请帮忙看一个两相流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;
}
}
不知道是不是这样写的？
加了以后，液面没有波动
多谢！！！

topshen

气液体系模型不对吧，dmp不好用的吧。试试其他模型。试试源相修改。