求大神看看能量方程源项udf问题出在哪里（源项由空化过程中物质汽化产生）

fyz1998

#include "udf.h"
DEFINE_SOURCE(energy1, cell, t, dS, eqn)
{
  Thread **pt;
  real re=199176;
  real MDv=1.138;
  real MDl=806.08;
  real Ps=C_P(cell, t);
real av=C_VOF(cell, pt[1]);
  real MDm=av*MDv+(1-av)*MDl;
  real RB=pow(3*av/((1-av)*4*3.14159*pow(10, 8)), 1/3);
  real Pv=0.0583*pow(C_T(cell,t), 4)-8.9018*pow(C_T(cell,t), 3)+517*pow(C_T(cell,t), 2)-16031*C_T(cell,t)+278654;
  real source;
  if(Ps<Pv)
    {
      source=-(3*av*(1-av)*MDv*MDl*pow(2*(Pv-Ps)/3/MDl, 1/2)/RB/MDm)*re;

      dS[eqn]=-(3*av*(1-av)*MDv*MDl*re/RB/MDm)*(1/(3*MDl*pow(2*(Pv-Ps)/3/MDl, 1/2)));
     }
  else
     {
       source=(3*av*(1-av)*MDv*MDl*pow(2*(Ps-Pv)/3/MDl, 1/2)/RB/MDm)*re;

       dS[eqn]=-(3*av*(1-av)*MDv*MDl*re/RB/MDm)*(1/(3*MDl*pow(2*(Ps-Pv)/3/MDl, 1/2)));
      }

  return source;


}

目前已经发现问题出在第九行（如何获得网格中气相体积/质量分数）。以上为错误的udf，通过添加能量方程源项后一步都算不了。但是把第九行气相分数改成常数可以计算，有没有大神教教小弟第九步到底怎么获得网格内的气相体积分数。能指明个大致方向就万分感谢了。