请教用udf写一个离散相和液相动量交换的源项

huguihua

各位老师、师兄、师姐：
    我想用udf写一个离散相和液相动量交换的源项，源项方程如附件所示：     
为液体的动力粘度系数， 为粒子的密度， 为粒子的直径， 为相对雷洛数， 和 分别为粒子和液体的速度， 为阻力系数（0.44）, 为粒子的质量流率, 为时间尺度(time step).我编的程序是：
include"udf.h"
#define VISC 1.003e-06
#define mp 0.068
#define Cd 0.44
DEFINE_DPM_SOURCE(dpm_source,c,t,S,strength,p)
{real Re,T;
Re=fabs(C_U(c,t)-P_VEL(p)[0])*P_DIAM(p)/VISC;
T=P_RHO(p)*pow(P_DIAM(p),2)*24/(18*VISC*C_R(c,t)*Cd*Re);
S->momentum_s[3]-=(C_U(c,t)-P_VEL(p)[0])*mp*P_DT(p)/T;}
我用compiled udf编译后没有出现错误，可是当我在边界条件fluid中输入dpm_source后进行初始化并跌代计算时，就出现了错误：get_udf_function dpm_source has wrong type:5!=19.我想请教各位我是不是编的程序有问题？请各位帮忙我编出正确的程序来！感激不尽！

hocolt

get_udf_function dpm_source has wrong type:5!=19!
我也遇到了同样的错误，个位大侠，能不能指教以下，还有搂住，我想你的问题肯定解决了吧，你能不能给我们说说呀！！
谢谢