我编写了fluent内置Schnerr and Sauer空化模型的udf文件，但是导入udf的结果和使用fluent内置Schnerr and Sauer空化模型的结果不一样，希望各位大佬帮我看看到底有什么问题。

模型公式：

#include "udf.h"

#include "sg_mphase.h" 

#define P_v 30000

DEFINE_MASS_TRANSFER(water_to_vapor,cell,thread,from_index,from_species_index,to_index,to_species_index)

{

   real  m_dot;

  Thread *gas, *liq;

  liq = THREAD_SUB_THREAD(thread, from_index);

  gas = THREAD_SUB_THREAD(thread, to_index);

  m_dot=0.0;

  if(C_P(cell, thread)<=P_v)

     m_dot=3108.725549*C_R(cell,liq)*C_R(cell,gas)/C_R(cell,thread)*sqrt(2.0/3.0*(P_v-C_P(cell,thread))/C_R(cell,liq));

  else

     m_dot=-621.74511*C_R(cell,liq)*C_R(cell,gas)/C_R(cell,thread)*sqrt(2.0/3.0*(C_P(cell,thread)-P_v)/C_R(cell,liq));

  return( m_dot);

}

还有

可以用c_v += 1.0*(1.0+0.1/20*fabs(0.5*(C_DUDX(cell,thread)+C_DUDY(cell,thread)+C_DVDX(cell,thread)+C_DVDY(cell,thread))))表示吗