用UDF中的DEFINE_PROFILE宏定义进流处速度剖面

Wind

[这个贴子最后由Wind在 2002/04/21 09:34pm 编辑]

用UDF中的DEFINE_PROFILE宏定义进流处速度剖面，3D。速度剖面是一曲面，沿高度呈指数变化，沿流向垂直的方向不变。其中用到F_CENTROID(x,f,t)宏，对于三维情况，参数x我选用二维数组，需要返回z坐标。刚开始看UDF,把manual中的例子改了一下，不知道这样写得对否，请不吝赐教。
/***********************************************************************/
/*UDF for specifying steady-state velocity profile boundary condition   */
/***********************************************************************/
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity,thread,index)
{
real x[ND_ND][ND_ND];          /* this will hold the position vector */
real z;
face_t f;
  
begin_f_loop(f,thread)
   {
    F_CENTROID(x,f,thread);
    z=x[1][1];
    F_PROFILE(f,thread,index)=40.*pow(z/200.,0.25);
   }
end_f_loop(f,thread)
}

Wind

编译出现问题，用过UDF的大侠请不吝赐教。

hongfw

x[ND_ND]这一数组不能改,它的意思为:x[0]代表x,x[1]代表y,x[2]代表z,在你的问题中一一对应就行!

Wind

3X。如果要定义一个3D问题，同时需要y,z的值，直接用y=x[1];z=x[2]就可以了吗？这一部分的内容在哪里有介绍？我看的是6.0的UDF Manual。

Wind

我参照你说的改了一下，编译通过了。但还有一个问题，在Fluent UDF Manual 5.4.2节中对F_CENTROID(x,f,t)有一个说明，x可以是1维，也可以是2维或三维数组。如果face f是一个面，那么我理解返回的x应是一个二维数组，那么这个二维数组的结构是怎样的？我还没有找到在哪里有说明。