边界条件UDF求助

Cavalier2010

边界条件UDF求助

#include "udf.h"

DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
  real x[ND_ND];                /* this will hold the position vector */
  real y;
  face_t f;

  begin_f_loop(f, thread)
    {
      F_CENTROID(x,f,thread);
      y = x[1];
      F_PROFILE(f, thread, position) = 20 - y*y/(.0745*.0745)*20;
    }
  end_f_loop(f, thread)
}
各位大侠，上述是Fluent帮助中关于二维速度进口为抛物线的UDF，那当三维计算域中进口速度方程式为u = 20(1-r2/0.012)时，温度为T = 303+100 ((r/0.01)2-( r/0.01)4)，速度u与温度T都是径向坐标r的的函数，请问怎么改Fluent自带的这个例子，能使其进口速度u与温度T随径向坐标而变化，请明白的大侠指点。