水沸騰算例UDF的詳細解析

1、【強烈推薦】水沸騰算例UDF的詳細解析【強烈推薦】水沸騰算例UDF的詳細解析#include udf.h /*包括常規(guī)宏*/#include sg_mphase.h /*包括體積分數(shù)宏*/ #define T_SAT 373 /*定義飽和溫度*/#define LAT_HT 1.e3 /*定義水蒸汽潛熱*/DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)/*定義液相源項*/Thread *mix_th, *sec_th; /*混合相、第二相 定義計算區(qū)域指針*/real m_dot_l; /*蒸發(fā)冷凝速率定義液相質(zhì)量轉(zhuǎn)移*/mix_th = THREA

2、D_SUPER_THREAD(pri_th); /*指向混合區(qū)的主相即液相的指針*/sec_th = THREAD_SUB_THREAD(mix_th, 1); /* 指向單相控制區(qū)的氣相的指針*/if(C_T(cell, mix_th)=T_SAT) m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;/*如果液相單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量轉(zhuǎn)移*/ dSeqn = -0.1*C_R(cell, pri_th)* fabs(C_T(cell, pri

3、_th) - T_SAT)/T_SAT; /*定義源項對質(zhì)量轉(zhuǎn)移偏導*/ else m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SAT; dSeqn = 0.; /*于是氣相向液相轉(zhuǎn)移，所以液相的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導為零*/ return m_dot_l;DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*定義氣相源項*/Thread * mix_th, *pri_th; /*混合相、第一相*/real m_dot_v;mix

4、_th = THREAD_SUPER_THREAD(sec_th); /*指向混合區(qū)的第二相即氣相的指針*/pri_th = THREAD_SUB_THREAD(mix_th, 0);/*指向單相控制區(qū)的液相的指針，液相為主相*/ if(C_T(cell, mix_th)=T_SAT) 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; /*如果混合區(qū)單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移*/ dSeqn = 0.; /*由于是液相向氣相轉(zhuǎn)移，所以氣相的

5、質(zhì)量源項對來自液相的質(zhì)量轉(zhuǎn)移的偏導為零*/ else m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SAT; dSeqn = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;/*由于是氣相向液相轉(zhuǎn)移，所以氣相的質(zhì)量源項對自身的質(zhì)量轉(zhuǎn)移的偏導不為零*/ return m_dot_v;DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) /*混合模型能量源項*/T

6、hread *pri_th, *sec_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) m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;/*如果混合區(qū)的單元溫度高于蒸發(fā)溫度。質(zhì)量轉(zhuǎn)移由液相向氣相轉(zhuǎn)移，吸熱，質(zhì)量轉(zhuǎn)移量前有負號*/ 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_SAT; dSeqn = -0.1*C_VO