圆柱体表面热流UDF求助

zhouchuanwen

我要给定圆柱体侧面热流的分布，是XY为底面，往Z方向延伸的圆管。这个热流分布是与周向的角度关的。依葫芦画瓢写了个UDF，不知道这能不能给出正确的侧面热流。主要疑问是在 F_CENTROID(x,f,thread);这个面中心给在哪了~~~~大家帮帮忙~~急~~

如果描述不清，下有附图...


#include "udf.h"
#define Q 933.7
DEFINE_PROFILE(flux, thread, position)
{
  real x[ND_ND];               
  real y;
  face_t f;

  begin_f_loop(f, thread)
    {
      F_CENTROID(x,f,thread);
      y = x[1];
      if(y>=-0.035 && y<-0.034235)
      F_PROFILE(f, thread, position) = 3921.57*Q*(y+0.035)+27.4*Q;
      if(y>=-0.034235 && y<-0.014124)
      F_PROFILE(f, thread, position) = 870*Q*(y+0.034235)+30.4*Q;
      if(y>=-0.014124 && y<0.0175)
      F_PROFILE(f, thread, position) = -1514.7*Q*(y-0.0175);
      else
      F_PROFILE(f, thread, position) = Q;
    }
  end_f_loop(f, thread)
}

[ 本帖最后由 zhouchuanwen 于 2011-11-22 14:44 编辑 ]

yulei97

个人感觉函数没有考虑z坐标是否合理，疑问中

zhouchuanwen

我也觉得有问题，只是不知道怎么把Z坐标考虑进去啊？能不能用极坐标什么的？