请各位帮我看一下我的出口温度分布udf的问题

gougouwu

已知出口温度分布，所以按模板弄了个出口温度分布的udf，interpreted通过之后，但计算发现这个udf没起作用。另外出口边界形状为一抛物线形，请给位大侠帮我看看，这个udf有什么问题。谢谢了~~

#include "udf.h"

DEFINE_PROFILE(shape_temperature,thread,position)
{
         real t[ND_ND];
         real x;
                 face_t f;

         begin_f_loop(f,thread)
         {
                 F_CENTROID(t,f,thread);
                 x=t[1];
                 F_PROFILE(f,thread,position)=272.346+6150.5*x*x-1.65923e8*x*x*x*x+2.00406e12*x*x*x*x*x*x;
         }
         end_f_loop(f,thread)
}

znon

原帖由 gougouwu 于 2011-3-2 17:10 发表

                               
登录/注册后可看大图

已知出口温度分布，所以按模板弄了个出口温度分布的udf，interpreted通过之后，但计算发现这个udf没起作用。

#include "udf.h"
...

出口边界条件是什么？
可能选择的出口边界条件是通过计算得来的，这样你自定义的温度分布就没有什么意义。

gougouwu

出口温度分布是实验得到的，就是关于x的一个分布

gougouwu

本来打算设定出口为壁面边界条件，满足marangoni滑移，和给定热流密度，但是设为wall之后fluent默认法相速度为零，而我这里又必须有法向速度，所以只能选择设为velocity inlet，给定温度分布，x，y方向的速度。

znon

原帖由 gougouwu 于 2011-3-3 15:51 发表

                               
登录/注册后可看大图

本来打算设定出口为壁面边界条件，满足marangoni滑移，和给定热流密度，但是设为wall之后fluent默认法相速度为零，而我这里又必须有法向速度，所以只能选择设为velocity inlet，给定温度分布，x，y方向的速度。

不了解具体问题。
我个人理解，Fluent中wall也有可以滑移的。另外也可以设定为wall，然后通过UDF更改临近wall这一层网格的法向速度，但这个比较麻烦一点。

herbert7017

同问！解决了吗？