DPM模型 添加外力

keith08

本人在用fluent DPM模型模拟空气和水绕过圆柱体的流动，水滴为离散相，现根据文献要在水滴的方程中加入外力，
               Fx=-2.4*(ux-4)/ √t
                Fy=-2.4*(uy)/ √t
我用DPM的body-force来实现，程序如下
#include "udf.h"
#define TSTART 0.0   /* field applied at t = tstart   */
DEFINE_DPM_BODY_FORCE(particle_body_force, p, i)
{
         real bforce;
         if(P_TIME(p)>TSTART)
           {
            if(i==0) bforce=-2.4*(P_VEL(p)[0]-4)/sqrt (P_TIME(p));   
            else if(i==1) bforce=-2.4*(P_VEL(p)[1]-0)/sqrt (P_TIME(p));
           }
         else
            bforce=0.0;
         /* an acceleration should be returned */
         return (bforce/P_MASS(p));
}
运行显示程序没问题，但是计算水滴轨迹时，即使把max number of steps 调到最大，最后也显示水滴轨迹计算未完成(incomplete)，求各位大神指教！