1、#include udf.h 包括常規(guī)宏#include sg_mphase.h 包括體積分數(shù)宏CVOF(C,T)#define T_SAT 373 定義蒸發(fā)溫度100#define LAT_HT 1.e3 定義蒸發(fā)潛熱J/KgDEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) 液相質量源項UDF Thread *mix_th, *sec_th; 定義計算區(qū)域線指針 real m_dot_l; 定義液相質量轉移 kg/(m2.s) mix_th = THREAD_SUPER_THREAD(pri_th); 指向混合區(qū)的主相即液相的指針 sec_th =

2、THREAD_SUB_THREAD(mix_th, 1);指向單相控制區(qū)的氣相的指針，氣相為第二相 if(C_T(cell, mix_th)=T_SAT) 如果液相單元的溫度高于蒸發(fā)溫度，液相向氣相的質量質量轉移 m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dSeqn = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; 定義源項對質量轉移偏導 ， else m_dot_l = 0

3、.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SAT;如果指向混合區(qū)液相的單元溫度小于蒸發(fā)溫度，氣相向液相的質量轉移，液相得 dSeqn = 0.; 由于是氣相向液相轉移，所以液相的質量源項對質量轉移的偏導為零return m_dot_l; DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) 氣相質量源項UDF Thread * mix_th, *pri_th; real m_dot_v; mix_th = THREAD_SUPER_THREAD(sec_t

4、h); 指向混合區(qū)的第二相即氣相的指針 pri_th = THREAD_SUB_THREAD(mix_th, 0); 指向單相控制區(qū)的液相的指針，液相為主相 if(C_T(cell, mix_th)=T_SAT) 如果混合區(qū)單元的溫度高于蒸發(fā)溫度，液相向氣相的質量質量轉移 m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dSeqn = 0.; 由于是液相向氣相轉移，所以氣相的質量源項對來自液相的質量轉移的偏導為零 else m_dot_v = -0.1*C_V

5、OF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SAT; 如果指向混合區(qū)的單元溫度小于蒸發(fā)溫度，氣相向液相的質量轉移 dSeqn = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; 由于是氣相向液相轉移，所以氣相的質量源項對自身的質量轉移的偏導不為零 return m_dot_v; DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) 混合模型能量源項UDF Thread *pri_th, *s

6、ec_th; real m_dot; pri_th = THREAD_SUB_THREAD(mix_th, 0);指向混合區(qū)的液相的指針 sec_th = THREAD_SUB_THREAD(mix_th, 1); 指向混合區(qū)的氣相的指針 if(C_T(cell, mix_th)=T_SAT) 如果混合區(qū)的單元溫度高于蒸發(fā)溫度。質量轉移由液相向氣相轉移，吸 熱，質量轉移量前有負號 m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dSeqn = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; else 相反，氣相向液相轉移則放熱 m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SA