#include "udf.h" //包括常规宏
#include "sg_mphase.h" // 包括体积分数宏 CVOF(C,T)
#define T_SAT 373 //定义蒸发温度 100℃
#define LAT_HT 1.e3 //定义蒸发潜热 J/Kg
DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) //液相质量源项 UDF
{
Thread *mix_th, *sec_th; //定义计算区线指针 real m_dot_l; //定义液相质量转移 kg/(m2.s)
mix_th = THREAD_SUPER_THREAD(pri_th); //指向混合区的主相即液相的指针
sec_th = THREAD_SUB_THREAD(mix_th, 1); //指向单相控制区的气相的指针,气相为第二相
if(C_T(cell, mix_th)>=T_SAT) //如果液相单元的温度高于蒸发温度,液相向气 相的质量质量转移
{