水蒸汽蒸發(fā)過程UDF程序段注釋

1、Fluent蒸發(fā)相變模擬UDF經(jīng)過幾天的不懈折騰，終于找到一個(gè)較為完成的用于fluent蒸發(fā)相變模擬的udf的一個(gè)程序。而且注釋相對完整。#include udf.h/包括常加宏#include sg_mphase h” / 包括體積分?jǐn)?shù)宏 CVOF(C, T)#define T.SAT 373/定義蒸發(fā)溫度 100Cdefine LAT.HT 1. e3/定義蒸發(fā)潛熱 J/KgDEFINE.SOURCE (liq_src, cell, pri.th, dS, eqn) /液相質(zhì)雖源項(xiàng) UDF Thread *mix_th, *sec_th;/定義計(jì)算區(qū)線指針real m_dot_l;/定義液

2、相質(zhì)量轉(zhuǎn)移kg/(m2, s)mix.th = THREAD_SUPER_THREAD(pri_th);/指向混合區(qū)的主相即液相的指針m.dot.l = -0. l*C_V0F(cell, pri_th)*C_R(cell, pri.th)* fabs(C_T(cell, mix.th) - T.SAT)/T.SAT;sec_th = THREAD_SUB_THREAD(mix_th, 1);/指向單相控制區(qū)的氣相的指針，氣相為第二相dSeqn = -0. l*C_R(cell, pri_th)*fabs(C_T(cell, mix.th) - T.SAT)/T.SAT;/定義源項(xiàng)對質(zhì)量轉(zhuǎn)移偏

3、導(dǎo)elsem_dot_l = 0. l*C_V0F(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell, mix_th)/T_SAT;/如柬指向混合區(qū)液而的單元溫度小于蒸麥溫度.氣相向液相的質(zhì)量轉(zhuǎn)移，液相得dSeqn = 0.于是氣相向液相轉(zhuǎn)移.所以液相的質(zhì)雖源項(xiàng)對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零return m dot 1:rDEFINE.SOURCE (vap.src, cell, sec_th, dS, eqn) /氣相質(zhì)雖源項(xiàng) UDFThread * mix_th, *pri_th;real m dot v;mix.th = THREAD_SUPE

4、R_THREAD(sec_th); /指向混合區(qū)的第二相即氣相的指針pri.th = THREAD_SUB_THREAD(mix_th, 0);指向單相控制區(qū)的液相的指針，液相為主相if(C_T(cell, mix.th)=T_SAT)如果混合區(qū)單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)雖質(zhì)量轉(zhuǎn)移m.dot.v = 0. l*C_V0F(cell, pri_th)*C_R(cell, pri.th)*fabs(C_T(cell, mix.th) - T.SAT)/T.SAT;dSeqn = 0.;由于是液相向氣相轉(zhuǎn)移，所以氣相的質(zhì)雖源項(xiàng)對來自液相的質(zhì)最轉(zhuǎn)移的偏導(dǎo)為零elsem.dot.v = -0

5、. l*C_V0F(cell, sec_th)*C_R(cell, sec.th)* fabs(T_SAT-C_T(cell, mix_th)/T_SAT;/如果指向混合區(qū)的機(jī)元溫度小于蒸發(fā)溫度，氣相向液相的質(zhì)雖轉(zhuǎn)移，氣相失dSeqn = -0. l*C_R(cell, sec.th)* fabs(C_T(cell, mix.th) - T.SAT)/T_SAT;/由于是氣相向液相轉(zhuǎn)移，所以氣相的質(zhì)覺源項(xiàng)對自身的質(zhì)雖轉(zhuǎn)移的偏導(dǎo)不為零 return m dot v;DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) 混合模型能雖源項(xiàng) CDFThread

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

7、. l*C_V0F(cell, pri_th)*C_R(cell, pri_th)/T_SAT;else/相反，氣相向液相轉(zhuǎn)移則放熱m.dot = 0. l*C_VOF(cell, sec_th)*C_R(cell, sec.th) *fabs(T_SAT-C_T (cell, mix.th)/T_SAT;dSeqn = 0. l*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;return LAT_HT*m_dot; /氣化潛熱與質(zhì)最轉(zhuǎn)移率相乘得能雖源項(xiàng)阡m3最終成果如下:1 06e-061 016-069.526-078 99e-O75 468-077.936-07 7 41e-07 6 88巳07 635875.B2e-O7 5 29e-07766-07 4.236-07 3.706-07 3 178-07 2 6