[求助]FLUENT作水蒸气高速流动

mumu8888

不知道怎么把水蒸汽凝结过程用UDF加进去，请大家帮帮忙[br][br][以下内容由 mumu8888 在 2008年05月17日 09:44am 时添加] [br]
[br][br][以下内容由 mumu8888 在 2008年05月17日 09:44am 时添加] [br]

mumu8888

#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373
#define LAT_HT -1.e3
#define c 0.1

DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn)
{
  Thread *mix_th, *sec_th;
  real m_dot_l;
  mix_th = THREAD_SUPER_THREAD(pri_th);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);


if(C_T(cell, mix_th)<=T_SAT)            /*小于饱和温度且大于饱和压力，主相蒸汽冷凝*/
  {
    m_dot_l = -c*C_VOF(cell, pri_th)*C_R(cell, pri_th)*         /*主相质量减少*/
                fabs(C_T(cell, pri_th) - T_SAT)/ T_SAT;
    dS[eqn] = -c*C_R(cell, pri_th)*
      fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
                               }
  else {
    m_dot_l =c*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = 0.;
       }
   
  return m_dot_l;
}
DEFINE_SOURCE(liq_src, cell, sec_th, dS, eqn)
{
  Thread * mix_th, *pri_th;
  real m_dot_v;
  mix_th = THREAD_SUPER_THREAD(sec_th);
  pri_th = THREAD_SUB_THREAD(mix_th, 0);

  if(C_T(cell, mix_th)<=T_SAT)      /*小于饱和温度且大于饱和压力，辅相水滴冷凝，*/
  {  
    m_dot_v =c*C_VOF(cell, pri_th)*C_R(cell, pri_th)*    /*辅相质量增加*/
                fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
    dS[eqn] = 0.;
                               }
  else {
    m_dot_v = -c*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = -c*C_R(cell, sec_th)*
      fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;
       }
   
  return m_dot_v;
}
DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)
{
  Thread *pri_th, *sec_th;
  real m_dot;
  pri_th = THREAD_SUB_THREAD(mix_th, 0);
  sec_th = THREAD_SUB_THREAD(mix_th, 1);
  if(C_T(cell, mix_th)<=T_SAT) {
    m_dot = -c*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;
    dS[eqn] =  -c*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
                               }
   
  else {
    m_dot = c*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = c*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}
  return LAT_HT*m_dot;
}

mumu8888

上面是我自己改的
定饱和温度 不考虑饱和压力变化 和潜热变化的程序