elbow管道的fluent模拟练习

入口速度要小，过大会出现涡流和回流，与实际例子不相符

#include "udf.h"

DEFINE_PROFILE(x_velocity,thread,nv)

{

  float x[3]; /* an array for the    coordinates */

  float y;

  face_t f;   /* f is a face    thread index */

  begin_f_loop(f, thread)

  {

     F_CENTROID(x,f,thread);

     y = x[1];

     F_PROFILE(f, thread, nv) = 20.*(1.-  y*y/(.0745*.0745));

  }

  end_f_loop(f, thread)

}

 

 

DEFINE_PROFILE 宏允许定义x_velocity函数

–所有的UDFs 以 DEFINE_ 宏开始

–x_velocity 将在 GUI中出现

–thread 和 nv DEFINE_PROFILE 宏的参数，
分别用来识别域和变量

–begin_f_loop宏通过thread指针，对所有的面f循环

F_CENTROID宏赋单元位置向量给 x[]

F_PROFILE 宏在面 f上施加速度分量

代码以文本文件保存
inlet_bc.c