[求助]两相流中的udf问题

rasca1

我用vof模型求解气液两相流
需要用DEFINE_DPM_BODY_FORCE给液体加体积力
代码如下：
#include "udf.h"
DEFINE_DPM_BODY_FORCE(particle_body_force,p,i)
{
real flow_time=RP_Get_Real("flow-time");
real accel;
if(i==0) accel=0.001*sin(5*flow_time);
else if(i==1) accel=0.0;
return(accel);
}
请问一下，需要添加什么才能使得该体积力作用到液体所在相？
非常感谢