24小时热门版块排行榜    

CyRhmU.jpeg
查看: 799  |  回复: 1

horstchen

银虫 (初入文坛)

[求助] fluent UDF给定进口速度,大家来找错了已有1人参与

如图所示为给定的进口速度,cpp文件如下:
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity, thread, position)
{
  real x[ND_ND]; /* this will hold the position vector */
  real xc,y;
  face_t f;
  begin_f_loop(f,thread)
{
  F_CENTROID(x, f, thread);
  xc=x[0];
  y=x[1];
  if(pow((xc*xc+(y-1)*(y-1))/(1.44*1.44),0.5)<=0.95&&pow((xc*xc+(y-1)*(y-1))/(1.44*1.44),0.5)>=0)
          F_PROFILE(f, thread, position)=10*pow(xc*xc+(y-1)*(y-1),0.5)/1.44;
  else
          F_PROFILE(f, thread, position)=0;
}
  end_f_loop(f, thread)
}


问题:为什么原点(0,0)的速度与周围的不太一致
根据
       if(pow((xc*xc+(y-1)*(y-1))/(1.44*1.44),0.5)<=0.95&&pow((xc*xc+(y-1)*(y-1))/(1.44*1.44),0.5)>=0)
          F_PROFILE(f, thread, position)=10*pow(xc*xc+(y-1)*(y-1),0.5)/1.44;

原点的速度应该与以(0,1)为圆心,半径为1的所有点速度是一样的
求各位虫朋友帮忙,看看哪儿出错了  
fluent UDF给定进口速度,大家来找错了
tangential velocity.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mwm621

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
horstchen: 金币+2, 有帮助, 单位应该没错,其他点的速度大小分布都符合 2014-07-04 09:00:56
我看小圆边缘与0.1是一致的,是不是楼主fluent里面的单位设置和udf中的不一致啊?楼主可以检查下哦~
2楼2014-07-04 08:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 horstchen 的主题更新
信息提示
请填处理意见