壁面运动函数

france

我的模型是由一个圆柱面和一个偏心圆柱面组成，其中内圆柱面是绕自身轴线旋转运动的同时沿轴向做直线运动。不知道该用动网格编写UDF还是用PROFILE边界条件来写。下面是我自己编写的，轴向速度是常数，所以我就将旋转运动转化到x，y方向，形成三个方向的运动。请大神指点。
#include "udf.h"
DEFINE_PROFILE(xv_profile,t,i)
{
face_t f;
real R=4.6e-3;
real W=1;
real pi=3.1415926;
real u;
real flow_time=RP_Get_Real("flow_time");
begin_f_loop(f,t)
{
u=W*R*cos((180*W*flow_time)/pi)
F_PROFILE(f,t,i)=u;
}
end_f_loop(f,t);
}

DEFINE_PROFILE(yv_profile,t,i)
{
face_t f;
real R=4.6e-3;
real W=1;
real pi=3.1415926;
real v;
real flow_time=RP_Get_Real("flow_time");
begin_f_loop(f,t)
{
v=W*R*sin((180*W*flow_time)/pi)
F_PROFILE(f,t,i)=v;
}
end_f_loop(f,t);
}