[Madreed]

I am defining a new UDF (k-epsilon) model to ansys fluent , every time I hook the diffusivity as user defined in the material panel I can't initialize the case and it gives me MPI finalize with error.

Does any one can help me with this problem ?

The diffusion equation :
DEFINE_DIFFUSIVITY(ke_diffusivity, c, t, eqn)
{
real diff;
real mu = C_MU_L(c,t);
real mu_t = C_R(c,t)*C_MU*f_mu(c,t)*SQR(C_UDSI(c,t,TKE))/C_UDSI(c,t,TDR);
switch (eqn)
{
case TKE:
diff = mu_t/SIG_TKE + mu;
break;
case TDR:
diff = mu_t/SIG_TDR + mu;
break;
default:
diff = mu_t + mu;
}
return diff;