水蒸气蒸发UDF.docx

文档介绍：对飞昂公司水一水蒸汽蒸发过程UDF程序段的解释
#iiiclude ""'包括常规宏
#iiiclude ""'包括体积分数宏 CVOF(C,T)
^define T_SAT 373 ，定义蒸发温度 对飞昂公司水一水蒸汽蒸发过程UDF程序段的解释
#iiiclude ""'包括常规宏
#iiiclude ""'包括体积分数宏 CVOF(C,T)
^define T_SAT 373 ，定义蒸发温度 100°C
define LAT_HT '定义蒸发潜热 J，Kg
DEFINE_SOURCE(liq_sic, cell, pn_th, dS, eqn)'液相质量源项 UDF
Tluead *nux_th. *sec_th;，定义计算区线指针
real m_dot_l;'定义液相质量转移kg/()
mix_th = THREAD_SUPER_THREAD(pri_th);'指向混合区的主相即液相的指针
sec_th = THREAD_SUB_THREAD(mix_th, 1);指向单相控制区的气相的指针，气相为第二 相
if(C_T(cell, inix_th)>=T_SAT)'如果液相单元的温度高于蒸发温度，液相向气相的质量 质量转移
{
= -0』*C_VOF(cell, pri_th)*C_R(cell, pn_th)*
I混合区液相温度-饱和温度I —
fhbs(C_T(cell. mix_th) - T_SAT)/T_SAT;
饱和温度
mdot =
dS[eqn] = -*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)门T_SAT;
'定义源项对质量转移偏导s4)=a + b(/>
else (
= 0』*C_VOF(celL sec_th)*C_R(cell. sec_th)* fhbs(T_SAT-C_T(celLmix_th))/T_SAT;
'如果指向混合区液相的单元温度小于蒸发温度，气相向液相的质量转移，液相得
mdot = x V，[相 x q 气相
I混合区液相温度-饱和温度I
X 饱和温度
dS[eqn] = 0.;'由于是气相向液相转移，所以液相的质量源项对质量转移的偏导为零
return m dot 1;
DEFINE_SOURCE(vap_src, cell. sec_th, dS, eqn)'气相质量源项 UDF
{
Tluead * mix_th, *pri_th;
real m dot v;
mix_th = THREAD_SUPER_THREAD(sec_th);'指向混合区的第二相即气相的指针
pn_th = THREAD_SUB_THREAD(mix_th, 0);指向单相控制区的液相的指针，液相为主相
if(C_T(cell, niix_th)>=T_SAT)'如果混合区单元的温度高于蒸发温度，液相向气相的质量 质量转移
m_do顷=*C_VOF(cell, pri_th)*C_R(cell. pii_th)*
fhbs(C_T(cell.