边界条件UDF求助

Cavalier2010

#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随径向坐标而变化，请明白的大侠指点。

JohnQiang

#include "udf.h"

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

  begin_f_loop(f, thread)
     {
       F_CENTROID(xx,f,thread);
       x = xx[0];
       y = xx[1];
        r =sqrt(x*x+y*y);
       F_PROFILE(f, thread, position) =20*(1-r*r/0.012);
     }
   end_f_loop(f, thread)
}
DEFINE_PROFILE(inlet_r_temperature, thread, position)
{
   real xx[ND_ND];                /* this will hold the position vector */
   real x,y,r;
   face_t f;

  begin_f_loop(f, thread)
     {
       F_CENTROID(xx,f,thread);
       x = xx[0];
       y = xx[1];
        r =sqrt(x*x+y*y);
       F_PROFILE(f, thread, position) =303+100* ((r/0.01)*(r/0.01)-( r/0.01));
     }
   end_f_loop(f, thread)
}
试试吧，
记得论坛里有过这样的帖子，你可以搜搜

Cavalier2010

大侠，牛人啊，请问F_CENTROID(xx,f,thread);
       x = xx[0];
       y = xx[1];
中的xx是什么意思?后面也出现了。

JohnQiang

原帖由 Cavalier2010 于 2012-7-23 18:52 发表

                               
登录/注册后可看大图

大侠，牛人啊，请问F_CENTROID(xx,f,thread);
       x = xx[0];
       y = xx[1];
中的xx是什么意思?后面也出现了。

惭愧，我也是菜鸟，只希望能帮到你。
xx[ND_ND]是数组，其中ND_ND是一个宏，对二维表示2，对三维表示3
xx是数组名字，自己定的。
详细的，请参考UDF帮助文件。