UDF-浮力

zisejingling

大家帮忙看一下下面这个UDF程序：是关于浮力的，即离散相气泡在连续相水中上升过程中受到的浮力：
（用离散相模型） Fb=pl(水的密度)*Vb（气泡体积）*g(重力加速度)
#include "udf.h"
#define g 9.81
DEFINE_DPM_BODY_FORCE(buoyancy,tp,i)
{
/* get the cell and thread that the particle is currently in */
cell_t c=RP_CELL(&(tp->cCell));
Thread *t=RP_THREAD(&(tp->cCell));
Particle *p;
real Vb=M_PI*pow(P_DIAM(p),3.0)/6.0,Fb;
P_RHO(p) = 998.2;
Fb=P_RHO(p)*Vb*g;
return Fb;
}
不知道对不对，大家看一下！
谢谢了！