UDF更改湍流粘度

pc123

图为需要将湍流粘度ut修改成的公式，ρm为混合相密度，ρv为气相密度，ρl为液相密度。n取10，Cu取0.0845。
参照UDF说明书，利用DEFINE_TURBULENT_VISCOSITY编写了如下程序，在mixture模型，k-e RNG中导入UDF。但是导入后程序初始化直接出现闪退情况。
求大神帮忙看看问题出在哪，是语句有问题，还是这个宏就不能出现多相流之类的东西，或者不能这么使用在MIXTURE中啊？求大神解答。

#include "udf.h"
DEFINE_TURBULENT_VISCOSITY(user_mu_t,c,t)
{
real mu_t;
real rho_mix ;
real k ;
real d ;
Thread *liq ;
Thread *vap ;
real rho_l ;
real rho_v ;
real f;
real cmu = 0.0845;

liq=THREAD_SUB_THREAD(t,0);
vap=THREAD_SUB_THREAD(t,1);
rho_mix = C_R(c,t);
rho_v = C_R(c,vap);
rho_l = C_R(c,liq);
k = C_K(c,t);
d = C_D(c,t);
f = rho_v + (pow((rho_mix - rho_v),10.0))/(pow((rho_l - rho_v),9.0));
mu_t = (cmu*f*k*k)/d;
return mu_t;
}

pc123

公式如图