24小时热门版块排行榜

返回列表

Tegrey

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 43.5
帖子: 29
在线: 14.3小时
虫号: 4489071
注册: 2016-03-10
专业: 内流流体力学

[求助] 关于冷凝UDF已有2人参与

找到的网上的关于水的冷凝UDF 如下

#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 373.15
#define LAT_HT 2260000
DEFINE_SOURCE(vap_src, cell, pri_th, dS, eqn)
{
  Thread *mix_th, *sec_th;
  real m_dot_v;
  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_v = 100*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, mix_th)-T_SAT)/T_SAT;
   dS[eqn] = 100*C_R(cell, sec_th)*fabs(C_T(cell, mix_th)-T_SAT)/T_SAT;
  }
  else
  {
   m_dot_v = -100*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
   dS[eqn] = 0.;
  }
return m_dot_v;
}
DEFINE_SOURCE(liq_src, cell, sec_th, dS, eqn)
{
  Thread *mix_th, *pri_th;
  real m_dot_l;
  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_l = 100*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_l = -100*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn] = -100*C_R(cell, sec_th)* fabs(C_T(cell, mix_th)-T_SAT)/T_SAT;
  }
  return m_dot_l;
}
DEFINE_MASS_TRANSFER(lic_gas_source,cell,thread,from_index,from_species_index,to_index,to_species_index)
{
Thread *gas, *liq;
real m_1g;
gas=THREAD_SUB_THREAD(thread,from_index);
liq=THREAD_SUB_THREAD(thread,to_index);
m_1g=0.;
if(C_T(cell,liq)>=T_SAT)
m_1g=100*C_VOF(cell,liq)*C_R(cell,liq)*fabs(C_T(cell,liq)-T_SAT)/T_SAT;
else
{
m_1g=100*C_VOF(cell,gas)*C_R(cell,gas)*fabs(T_SAT-C_T(cell,gas))/T_SAT;
}
return m_1g;
}
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 = -100*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(C_T(cell, mix_th)-T_SAT)/T_SAT;
dS[eqn] = -100*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;
  }
  else
  {
m_dot = 100*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn] = 100*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
  }
  return LAT_HT*m_dot;
}

有几个问题，第一在mass_tansfer中
if(C_T(cell,liq)>=T_SAT)
m_1g=100*C_VOF(cell,liq)*C_R(cell,liq)*fabs(C_T(cell,liq)-T_SAT)/T_SAT;
else
{
m_1g=100*C_VOF(cell,gas)*C_R(cell,gas)*fabs(T_SAT-C_T(cell,gas))/T_SAT;
}
m_1g这个的符号都是正号码？？

还有一个问题是应用这个UDF计算后，不应该低于373K的时候不应该就是液相了么，但是仿真算出来虽然有相变，
但是气相区域还存在许多低于373K的区域，，，  求UDF 大神指点

回复此楼