UFD动量源项求助帖

qwzx93

为了计算声场，定义了一个UDF动量源项，在某个zone内施加一个动量源。
x[1]的范围在0.066到0.036，x[0]的范围在0.066到0.049之间。但是每次计算出来的结果和没加源项的结果一样，请看一下UDF编译的有没有问题？



DEFINE_SOURCE(mhd_x,c,t,dS,eqn)
{
  real x[ND_ND];
  real source,E;
  C_CENTROID(x,c,t);
  E=22.627-654*(0.066-x[1])-1308*(x[0]-0.049);
  if(E<=3.)
  {source=dS[eqn]=0.;}
  else
  {source=288*E;
   dS[eqn]=0.;}
   return source;
}