fluent

    昨天分享了fluent udf help文件中关于入口velocity的例子的实现，帮助文件中例子是二维模型，为了更好地理解的udf代码的编写，今天做了个最简单的三维模型的例子。

模型如图所示，一个简单管道，坐标系原点位于入口面的圆心。设置入口处速度随圆半径的增大而程抛物线分布，也即靠近管壁出速度为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，结果如下所示：

速度分布图

中心平面上的动态压力分布图

经过二维和三维模型的udf实例实现，现在对fluent udf编程有了初步的认识，赶快尝试你的分析求解吧。