找高手看一下关于FLUENT定义流场密度的UFD程序

qq347388344

想通过UDF来定义流场密度关于压力的函数，自己编写了一个
但运行起来老是错误
Error：Divergence detected in AMG solver：pressure correction

以下是我自己的程序，想请高手看一下程序编写是否符合规范
DEFINE_PROPERTY(density, c, t)
{
        real rho;
        real p;
        real x1, x2, x[ND_ND];
        real temp;
   
  C_CENTROID(x,c,t);
  x1 = x[0];
  x2 = x[1];
  p = C_P(c,t);
   
  temp = pow(x2/Hy,3.0)-3*x2/Hy;
  rho = 1/Kappa*(p-0.096*(pow(x2,2.0)-pow(x1,2.0))-Kappa/4.0*pow(x1,2.0)*pow(temp,2.0));
  return rho;
}