24小时热门版块排行榜

返回列表

li2hebe

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 4
帖子: 8
在线: 4小时
虫号: 2926222
注册: 2014-01-11
专业: 多相流热物理学

[交流] fluent问题急求大神帮助！！！已有4人参与

最近在模拟饱和蒸汽倾斜管内凝结换热的例子，可是液相体积分数一直是在10e-7左右徘徊，udf就是用的网上搜到的蒸发冷凝模型，另外还有个质量输运方程，具体见下方，然后我设置气相为主相液相为次相，用的mix模型，有大神能帮我解答一下吗，要是有热心人能帮帮忙模拟一下最好了，小弟感激不尽，因为是初学者，所以好多地方一头雾水，谢谢大家了！管子模型是0.2m×10m，2维模型，倾角60°，管内压力15kpa恒定，蒸汽温度327K，壁温293K，入口质量边界条件0.1kg/s，出口压力边界条件出口混合物温度325K

#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 327
#define LAT_HT 2.e3

DEFINE_SOURCE(liq_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 = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS[eqn] = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

}  else {
m_dot_l = 0.1*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(vap_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 = 0.1*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 = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;

dS[eqn] = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;

}
  return m_dot_v;
}
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=-0.1*C_VOF(cell,liq)*C_R(cell,liq)*fabs(C_T(cell,liq)-T_SAT)/T_SAT;
else{
if(m_1g==0.)
m_1g=0.1*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 = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;
dS[eqn] =  -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
}
else {
m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
dS[eqn] = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;
}
  return LAT_HT*m_dot;
}

[ Last edited by li2hebe on 2014-2-2 at 17:00 ]

回复此楼