壓力控制的相變udf

1、#include “udf.h“#include “sg_mphase.h“#define LAT_HT 1.e3DEFINE_SOURCE(air_src,cell,pri_th,dS,eqn) /*主相濕空氣質(zhì)量輸運源項*/{Thread *mix_th,*sec_th; /*計算區(qū)域指針*/real m_v; /*主相質(zhì)量輸運量*/real mf

2、v; /*濕空氣中水蒸氣的質(zhì)量分?jǐn)?shù)*/real t_air; /*濕空氣的溫度*/real hr; /*濕空氣的含濕量*/real pa; /*濕空氣的總壓力*/real pa_v; /*濕空氣的水蒸氣分壓*/real mfs; /*licl 溶液的質(zhì)量分?jǐn)?shù)*/real t_sol; /*licl 溶液的溫度*/real ps_w; /*licl 溶液的水蒸

3、氣分壓*/mix_th=THREAD_SUPER_THREAD(pri_th); /*混合區(qū)液相指針*/sec_th=THREAD_SUB_THREAD(mix_th,1); /*單相區(qū)氣相指針*/mfv=C_YI(cell,pri_th,0); /*獲得濕空氣中水蒸氣的質(zhì)量分?jǐn)?shù),注意species 列表中的順序，水是第一個時序號才是 0,dS 要對該項求導(dǎo)*/hr=mfv/(1-mfv);

4、 /*由濕空氣中水蒸氣的質(zhì)量分?jǐn)?shù)算到濕空氣的含濕量*/pa=ABS_P(C_P(cell,pri_th),op_pres); /*獲取濕空氣的總壓力*/pa_v=(pa*hr)/(hr+0.622); /*獲取濕空氣水蒸氣分壓*//* t_air=C_T(cell,pri_th); /*濕空氣的溫度，單位為攝氏度*/mfs=C_YI(cell,sec_th,0); /*獲得溶液中 licl 的質(zhì)量分?jǐn)?shù),注意 specie

5、s 列表中的順序，licl 是第一個時序號才是 0,*/t_sol=C_T(cell,sec_th); /*溶液的溫度，單位為攝氏度*/ps_w=1.1491997+0.11242725*t_sol-0.06183172*mfs+0.0014521543*pow(t_sol,2)+0.0011080673*pow(mfs,2)-0.0035782429*t_sol*mfs; /*溶液水蒸氣分壓*/if(pa_v>ps_w)

6、 /*如果濕空氣的水蒸氣分壓大于溶液的水蒸氣分壓，則水蒸氣向溶液移動，冷凝放熱，屬于除濕過程*/{m_v=0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w); /* mass transfer rate in dehumidifier */dS=0.00165*((pa*(0.378*mfv+0.622)-pa*mfv*0.378)/pow(0.378*mfv+0.622,2)

7、); /*濕空氣質(zhì)量源項對濕空氣的質(zhì)量轉(zhuǎn)移求導(dǎo)，對濕空氣中水蒸氣質(zhì)量分?jǐn)?shù) mfv 求導(dǎo)*/}else if(pa_vps_w) /*如果濕空氣的水蒸氣分壓大于溶液的水蒸氣分壓，則水蒸氣向溶液移動，冷凝放熱，屬于除濕過程*/{m_s=0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w); /* mass transfer rate in dehumidifier */dS=

8、0; /*水蒸氣由濕空氣轉(zhuǎn)移到溶液，溶液的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為 0*/}else if(pa_v<ps_w) /*如果濕空氣的水蒸氣分壓小于溶液的水蒸氣分壓，則水蒸氣向濕空氣移動，蒸發(fā)吸熱，屬于再生過程*/{m_s=-0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w); /* mass transfer rate in dehumidifier */dS=-0

9、.00165*(0.06183172-0.0022161346*mfs+0.0035782429*t_sol); /*dS：溶液的質(zhì)量源項對質(zhì)量轉(zhuǎn)移求偏導(dǎo)，對溶液中水蒸氣質(zhì)量分?jǐn)?shù) mfs 求導(dǎo)*/}return m_s;}DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) /*混合相能量方程輸運源項*/{Thread *pri_th,*sec_th;real m_e; /*相變質(zhì)量*/real

10、 mfv;real t_air; /*濕空氣的溫度*/real hr; /*濕空氣的含濕量*/real pa; /*濕空氣的總壓力*/real pa_v; /*濕空氣的水蒸氣分壓*/real mfs; /*licl 溶液的質(zhì)量分?jǐn)?shù)*/real t_sol; /*licl 溶液的溫度*/real ps_w; /*licl 溶液的水蒸氣分壓*/pri_th=THREAD_SUB_THREAD(mix_th,0);

11、 /*混合區(qū)液相指針*/sec_th=THREAD_SUB_THREAD(mix_th,1); /*單相區(qū)氣相指針*/mfv=C_YI(cell,pri_th,0); /*獲得濕空氣中水蒸氣的質(zhì)量分?jǐn)?shù),注意species 列表中的順序，水是第一個時序號才是 0*/hr=mfv/(1-mfv); /*由濕空氣中水蒸氣的質(zhì)量分?jǐn)?shù)算到濕空氣的含濕量*/pa=ABS_P(C_P(cell,pri_th