|
发表于 2023-3-2 18:39:03
|
显示全部楼层
#include "udf.h"
#include "sg_mphase.h"
#define T_SAT 377 //饱和温度
#define VOF_SAT 0.5 //定义边界
#define LAT_HT 2.257e6//相变潜热1000j/kg
#define M_W 0.018 //水的摩尔质量18g/mol
#define g_c 8.314 //气体分子常数8.314j/(mol*k)
#define v_md 0.58979 //水蒸气的密度
//液体首相的传质
DEFINE_SOURCE(liq_src,cell,pri_th,dS,eqn)
{
Thread *mix_th, *sec_th;
real m_dot_l;
//mix_th
mix_th =THREAD_SUPER_THREAD(pri_th);
sec_th =THREAD_SUB_THREAD(mix_th,1);
//混合相温度导致相变
m_dot_l= -5*(C_T(cell,pri_th)-T_SAT)*v_md*LAT_HT*pow((1-C_VOF(cell,pri_th))*C_VOF(cell,pri_th),2)/(T_SAT*sqrt(2*3.14*g_c*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
mix_th =THREAD_SUPER_THREAD(sec_th);
pri_th =THREAD_SUB_THREAD(mix_th,0);
//混合相温度导致相变
m_dot_v= 5*(C_T(cell,pri_th)-T_SAT)*v_md*LAT_HT*pow((1-C_VOF(cell,pri_th))*C_VOF(cell,pri_th),2)/(T_SAT*sqrt(2*3.14*g_c*T_SAT));
//对体积分数求导
dS[eqn]= 0;
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);
//混合相温度导致相变
m_dot= -5*(C_T(cell,pri_th)-T_SAT)*v_md*LAT_HT*pow((1-C_VOF(cell,pri_th))*C_VOF(cell,pri_th),2)/(T_SAT*sqrt(2*3.14*g_c*T_SAT));
//对体积分数求导
dS[eqn]= 0;
//热量的传递
return LAT_HT*m_dot;
} |
|