1、Fluent 蒸發(fā)相變 模擬 UDF經(jīng)過幾天的不懈折騰，終于找到一個較為完成的用于fluent蒸發(fā)相變模擬的udf的一個程序。而且注釋相對完整。#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/Kg DEFINE_SOURCE

2、(liq_src, cell, pri_th, dS, eqn)     /液相質(zhì)量源項UDF Thread *mix_th, *sec_th;    /定義計算區(qū)線指針 real m_dot_l;     /定義液相質(zhì)量轉(zhuǎn)移 kg/(m2.s) mix_th = THREAD_SUPER_THREAD(pri_th);        /指向混合區(qū)的主相即液相的指針 sec_th = THREAD_SUB_THREAD(mix_th, 1);  

3、     /指向單相控制區(qū)的氣相的指針，氣相為第二相 if(C_T(cell, mix_th)>=T_SAT)           /如果液相單元的溫度高于蒸發(fā)溫度，液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移                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

4、;   dSeqn = -0.1*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT;  /定義源項對質(zhì)量轉(zhuǎn)移偏導(dǎo)        else                           m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T

5、_SAT-C_T(cell,mix_th)/T_SAT;   /如果指向混合區(qū)液相的單元溫度小于蒸發(fā)溫度，氣相向液相的質(zhì)量轉(zhuǎn)移，液相得                   dSeqn = 0.;/由于是氣相向液相轉(zhuǎn)移，所以液相的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零       return m_dot_l; DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /氣相質(zhì)量源項UDF Thread

6、* mix_th, *pri_th; real m_dot_v; mix_th = THREAD_SUPER_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)移      &#

7、160;             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.;由于是液相向氣相轉(zhuǎn)移，所以氣相的質(zhì)量源項對來自液相的質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零            else  

8、;                       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;                /如果指向混合區(qū)的單元溫度小于蒸發(fā)溫度，氣相向液相的質(zhì)量轉(zhuǎn)移，氣相失         

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

10、ead *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)移量前有負號               m_dot = -0.1*

11、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         /相反，氣相向液相轉(zhuǎn)移則放熱                 m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th)/T_SAT;