为什么导入UDF后与没导入前算的结果一样
悬赏金币:5 zxf_19970610-zxf 时间:2020-11-20 17:06
#include "math.h"
#define n 17500
#define r 0.025
#define x0 1.3774e-5
#define y0 -5.9393e-6
#define v 0.0001
DEFINE_PROFILE(x_velcity,t,i)
{
face_t f;
real c[ND_ND];
real current_timestep;
current_timestep= RP_Get_Real("physical-time-step");
real x,y,angle,u_x,w;
begin_f_loop(f,t);
{
F_CENTROID(c,f,t);
x=c[0];
y=c[1];
w=(n*3.1415926)/30;
angle=atan(fabs(y-(y0+(v*current_timestep)))/fabs(x-(x0+(v*current_timestep))));
u_x=w*r*sin(angle);
if(y>(y0+(v*current_timestep)))
u_x=-u_x;
F_PROFILE(f,t,i)=u_x;
}
end_f_loop(f,t)
}
DEFINE_PROFILE(y_velcity,t,i)
{
face_t f;
real c[ND_ND];
real current_timestep;
current_timestep = RP_Get_Real("physical-time-step");
real x,y,angle,u_y,w;
begin_f_loop(f,t);
{
F_CENTROID(c,f,t);
x=c[0];
y=c[1];
w=(n*3.1415926)/30;
angle=atan(fabs(y-(y0+(v*current_timestep)))/fabs(x-(x0+(v*current_timestep))));
u_y=w*r*cos(angle);
if(x<(x0+(v*current_timestep)))
u_y=-u_y;
F_PROFILE(f,t,i)=u_y;
}
end_f_loop(f,t)
}
回答(0)