[求助]关于body_force的udf问题

rasca1

气液两相问题，需要给液体相加体力f＝Agsin（wt）
代码如下：
#include "udf.h"
#define ample 1.0
#define grav 9.81
#define omega 5.0
DEFINE_DPM_BODY_FORCE(particle_body_force,p,i)
{
int phase_domain_index;
real accel;
real current_time;
current_time=RP_Get_Real("flow-time");
Domain *subdomain;
Domain *mixture_domain;
subdomain = Get_Domain(3);
{
if(DOMAIN_ID(subdomain)==3)
{
if(i==0) accel=-ample*grav*sin(omega*current_time);
else if(i==1) accel=0.0;
}
return accel;
}
}
不知道是不是这样写的
加了以后，液面没有波动，不知道该怎么修改
多谢！

rasca1

为什么都没有人回呢？
烦请做过类似问题的xdjm多多指教了!