24小时热门版块排行榜

返回列表

风_秋_黄

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 304
帖子: 13
在线: 9.3小时
虫号: 2466127
注册: 2013-05-15
专业: 多相流热物理学

[求助] fluent相变模拟（求大神指点）急用，谢谢已有1人参与

我做的平板流动换热，低温流体流过高温（超过沸点）平板，模拟平板温度变化以及流体相变，但是模拟的相变有气相但是气泡不长大，不知是哪里出问题了？
平板（有一定厚度）有个初温（大于液体沸点），会随着流体流过而下降，流体被加热，超过沸点会汽化。程序如下：
# include "udf.h"
# include "sg_mphase.h"
# define T_SAT 373.15
# define lam 0.1
# define LAT 2250
  DEFINE_PROPERTY(surface_tem,c,t)
{
real surface;
real temp=C_T(c,t);
surface=0.09537-2.24*(pow(10,-6))*temp-2.56*(pow(10,-7))*pow(temp,2);
return surface;
}

DEFINE_SOURCE(liq_src,cell,first,dS,eqn)
{
Thread *mixer,*second;
real m_dot_first;
mixer=THREAD_SUPER_THREAD(first);
second=THREAD_SUB_THREAD(mixer,1);
if(C_T(cell,mixer)>=T_SAT)
{
m_dot_first=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,mixer)-T_SAT)/T_SAT;
      dS[eqn]=-lam*C_R(cell,first)*fabs(C_T(cell,mixer)-T_SAT)/T_SAT;
}
else
{
m_dot_first=lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,mixer)-T_SAT)/T_SAT;
      dS[eqn]=0;
}
return m_dot_first;
}

DEFINE_SOURCE(vap_src,cell,second,dS,eqn)
{
Thread *mixer,*first;
real m_dot_second;
mixer=THREAD_SUPER_THREAD(second);
first=THREAD_SUB_THREAD(mixer,0);
if(C_T(cell,mixer)>=T_SAT)
{
m_dot_second=lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,mixer)-T_SAT)/T_SAT;
      dS[eqn]=0;
}
else
{
m_dot_second=-lam*C_VOF(cell,second)*C_R(cell,second)*fabs(C_T(cell,mixer)-T_SAT)/T_SAT;
      dS[eqn]=-lam*C_R(cell,second)*fabs(C_T(cell,second)-T_SAT)/T_SAT;
}
return m_dot_second;
}

DEFINE_SOURCE(enrg_src,cell,mixer,dS,eqn)
{
Thread *first,*second;
real m_dot;
first=THREAD_SUB_THREAD(mixer,0);
second=THREAD_SUB_THREAD(mixer,1);
if(C_T(cell,mixer)>=T_SAT)
{
m_dot=-lam*C_VOF(cell,first)*C_R(cell,first)*fabs(C_T(cell,first)-T_SAT)/T_SAT;
      dS[eqn]=-lam*C_VOF(cell,first)*C_R(cell,first)/T_SAT;
}
else
{
m_dot=lam*C_VOF(cell,second)*C_R(cell,second)*fabs(C_T(cell,mixer)-T_SAT)/T_SAT;
      dS[eqn]=-lam*C_VOF(cell,second)*C_R(cell,second)/T_SAT;
}
return LAT*m_dot;
}

无标题.png

回复此楼