奇怪：同一个UDF k-e湍流模型能用，k-w湍流模型不能用？？？？

shangdong

#include "udf.h"
real Re_k    = 6.0;
real C_0     = 13.0;
real C_s(cell_t c, Thread *t)
{return fabs(C_U(c,t)*(C_U(c,t)*C_DVDX(c,t)+C_V(c,t)*C_DVDY(c,t))-C_V(c,t)*(C_U(c,t)*C_DUDX(c,t)+C_V(c,t)*C_DUDY(c,t)))
/pow((C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t)),1.5);}
real R_l(cell_t c, Thread *t)
{return C_s(c,t)*sqrt(C_MU_L(c,t)/C_R(c,t)/C_D(c,t));}
real alph_1(cell_t c, Thread *t)
{return C_0*R_l(c,t)/(0.1+R_l(c,t));}
real Re_t(cell_t c, Thread *t)
{return C_R(c,t)*C_K(c,t)/C_MU_L(c,t)/C_D(c,t);}
real alph(cell_t c, Thread *t)
{return (0.025+alph_1(c,t)+Re_t(c,t)/Re_k)/(1+Re_t(c,t)/Re_k);}
/* Set the turbulent viscosity */
DEFINE_TURBULENT_VISCOSITY(user_mu_t,c,t)
{
real mu_t;
real rho = C_R(c,t);
real k = C_K(c,t);
real d = C_D(c,t);
real a = alph(c,t);
mu_t = a*rho*k/d;
return mu_t;
}
计算湍流粘性的
k-e湍流模型能用也能算，但k-w湍流模型只能编译，一算就出错。
请问斑竹看一下
多谢

Joywe

采用K_W模型的时候，你不能访问C_D(c,t),肯定又是fatal error 和四句话，呵呵....

youy624

程序中的real d = C_D(c,t);是指扩散系数，即k-epsilon模型中的epsilon，而k-omega模型中没有这个物理量的。故报错。