版主 救命啊

captain-jack

大家谁有关于UDF定义材料的
我要把cross-WLF粘度模型加入FLUENT，但我编写的程序有问题，大家帮忙看一下 谢谢了
我觉得是我在调用C_W（cell，thread）时，有问题
运行时，FLUENT显示错误，见附近图
请教各位达人  我先谢谢了

#include"udf.h"
DEFINE_PROPERTY(cell_viscosity,cell,thread)
{
  real mu_lam;
  real nd,nc;
  real temp=C_T(cell,thread);
  real tempv=C_W(cell,thread);
  mu_lam=1.129e+014*exp(-34.03*(temp-373.15)/(51.6+(temp-373.15)));
  nc=mu_lam*tempv/258100;
  nd=mu_lam/(1+pow(nc,0.6));
  return nd;
}
DEFINE_PROPERTY(cell_density,cell,thread)
{
real Tm,tempP,tempB,tempV,tempT,tempD;
tempP=C_P(cell,thread);
Tm=428.15+4.127e-007*tempP;
tempT=C_T(cell,thread);
if(tempT>=Tm)
    {
      tempV=0.001045+6.759e-007*(tempT-428.15);
      tempB=1.309e+008*exp(tempT-428.15);
    }
else
  {
   tempV=0.001033+2.475e-007*(tempT-428.15);
   tempB=2.56e+008*exp(tempT-428.15);
  }
     tempD=pow(tempV*(1-log(1+tempP/tempB)),-1);
  
  
  return tempD;
}