如图所示为给定的进口速度,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 |