动量源项 多相流

rouyuan朱

源项通过编译
但是一加上就算不了了
求大伙儿帮我瞅瞅啊！
万分谢

DEFINE_SOURCE(zmom_source,cell,mix_th,ds,eqn)
{
  Thread *pri_th, *sec_th;
  real zmom;
  real r;

  pri_th = THREAD_SUB_THREAD(mix_th,0); //指向混合区的液相
sec_th = THREAD_SUB_THREAD(mix_th,1);//指向混合区的气相


  r=g/(C_P(cell,sec_th)-C_P(cell,pri_th));

  if (C_VOF(cell,pri_th)>0)// 区域内全是液体，用液体的源项计算
  {
  zmom=r*n*frel*C_U(cell,mix_th)/(2*D);
  ds[eqn]=r*n*frel/(2*D);
  }
  else //液体体积分数为0，全是蒸汽，用蒸汽的源项计算
  {
  zmom=r*n*frev*C_U(cell,mix_th)/(2*D);
  ds[eqn]=r*n*frev/(2*D);
  }

  return zmom;

}