#include "udf.h"
#include "sg_mphase.h"
DEFINE_SOURCE(air_src, cell, pri_th, dS, eqn)
{
Thread *mix_th,*sec_th,* thi_th;
real m_dot_1;
mix_th = THREAD_SUPER_THREAD(pri_th);
sec_th = THREAD_SUB_THREAD(mix_th, 1);
thi_th = THREAD_SUB_THREAD(mix_th, 2);
m_dot_1 = 0;
dS = 0;
return (m_dot_1);
}
DEFINE_SOURCE(water_vap, cell, sec_th, dS, eqn)
{
Thread * mix_th, *thi_th; /*
混合相，第三相
*/
real m_dot_vap; /*
定义水汽传递率
*/
real T_w; /*
定义露点温度
*/
mix_th = THREAD_SUPER_THREAD(sec_th); /*
指向第二相
*/
thi_th = THREAD_SUB_THREAD(mix_th, 2); /*
指向第三相
*/
T_w=8.22+12.4*log(C_P(cell,sec_th))+1.9*log(C_P(cell,sec_th))*log(C_P(c
ell,sec_th));
if(C_T(cell, mix_th)<=T_w) {
m_dot_vap = -0.01*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_w-
C_T(cell,mix_th))/T_w; /*
低于露点水蒸气传递至液态水
*/
dS[eqn] = -0.01*C_R(cell, sec_th)*fabs(T_w-C_T(cell,mix_th))/T_w;
}
else {
m_dot_vap = 0.01*C_VOF(cell, thi_th)*C_R(cell,
thi_th)*fabs(C_T(cell,mix_th)-T_w)/T_w ; /*
高于露点液态水传递至水蒸气
*/
dS [eqn] = 0;
}
return (m_dot_vap);
}
DEFINE_SOURCE(water_liq, cell, thi_th, dS, eqn)
{
Thread *mix_th, *sec_th; /*
混合相，第二相
*/
real m_dot_liq; /*
定义水汽传递率
*/
real T_w; /*
定义露点温度
*/
mix_th = THREAD_SUPER_THREAD(thi_th); /*
指向第三相
*/
sec_th = THREAD_SUB_THREAD(mix_th, 1); /*
指向第二相
*/
