模型如图所示，一个简单管道，坐标系原点位于入口面的圆心。设置入口处速度随圆半径的增大而程抛物线分布，也即靠近管壁出速度为0，圆心出速度最大，设为0.1m/s。

udf代码修改如下：

#include "udf.h"

DEFINE_PROFILE(inlet_x_velocity, thread, position)

{

   real x[ND_ND];   二维中该数组有两个元素，存储x，y坐标；三维模型，该数组存储x,y,z三个坐标。

   real r,R;

   face_t f;

   R=0.0075;   管道半径为7.5mm

   begin_f_loop(f,thread)

   {

    F_CENTROID(x, f, thread);

    r=sqrt(x[1]*x[1]+x[2]*x[2])/R;  x[1],x[2]分别为入口面上face单元质心的y坐标值和z坐标值

    F_PROFILE(f, thread, position)=0.1*(1.0-r*r);

   }  

   end_f_loop(f,thread)

}

设置求解为瞬态求解，液体为液态水，k-e模型，求解时间10s，结果如下所示：

速度分布图

From

http://blog.sina.com.cn/s/blog_15ae0ebcb0102wof8.html