壓力控制的相變UDF

1、#include #include #define LAT_HTDEFINE_SOURCE(air_src,cell,pri_th,dS,eqn)Thread *mix_th,*sec_th;realreal質(zhì)量分?jǐn)?shù)*/realrealrealrealrealrealreal/*主相濕空氣質(zhì)量輸運(yùn)源項(xiàng)*/m_v;mfv;t_air;hr;pa;pa_v;mfs;t_sol;ps_w；/*濕空氣的溫度*/*濕空氣的含濕量*/*濕空氣的總壓力*/*濕空氣的水蒸氣分壓*/*licl溶液的質(zhì)量分?jǐn)?shù)*/*licl溶液的溫度*/*licl溶液的水蒸氣分壓*/*計(jì)算區(qū)域指針*/*主相質(zhì)量輸運(yùn)量*/*濕空氣中

2、水蒸氣的/*混合區(qū)液相指針*/*單相區(qū)氣相指針*/*獲得濕空氣中水蒸氣的質(zhì)量分mix_th=THREAD_SUPER_THREAD(pri_th);sec_th=THREAD_SUB_THREAD(mix_th,1);mfv=C_YI(cell,pri_th,0);數(shù)，注意species列表中的順序，水是第一個時序號才是0,dS要對該項(xiàng)求導(dǎo)*/hr=mfv/(1-mfv);/*由濕空氣中水蒸氣的質(zhì)量分?jǐn)?shù)算到濕空氣的含濕量*/pa=ABS_P(C_P(cell,pri_th),op_pres);/*獲取濕空氣的總壓力*/pa_v=(pa*hr)/(hr+;/*獲取濕空氣水蒸氣分壓*/* t_ai

3、r=C_T(cell,pri_th);/*濕空氣的溫度，單位為攝氏度*/mfs=C_YI(cell,sec_th,0);/*獲得溶液中l(wèi)icl的質(zhì)量分?jǐn)?shù)，注意species列表中的順序，licl是第一個時序號才是0,*/t_sol=C_T(cell,sec_th);/*溶液的溫度，單位為攝氏度*/ps_w=+0.*mfs+*pow(t_sol,2)+*pow(mfs,2)*t_sol*mfs;/*溶液水蒸氣分壓*/if(pa_vps_w)/*如果濕空氣的水蒸氣分壓大于溶液的水蒸氣分壓，則水蒸氣向溶液移動，冷凝放熱，屬于除濕過程*/m_v=*(pa*mfv/(1-mfv)/(mfv/(1-mfv

4、)+)-ps_w);/* masstransfer rate in dehumidifier */dS=*(pa*mfv+-pa*mfv*/pow*mfv+,2);/*濕空氣質(zhì)量源項(xiàng)對濕空氣的質(zhì)量轉(zhuǎn)移求導(dǎo)，對濕空氣中水蒸氣質(zhì)量分?jǐn)?shù)mfv求導(dǎo)*/else if(pa_vps_w)/*如果濕空氣的水蒸氣分壓大于溶液的水蒸氣分壓，則水蒸氣向溶液移動，冷凝放熱，屬于除濕過程*/m_s=*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+)-ps_w)；/* masstransfer rate in dehumidifier */dS=0；/*水蒸氣由濕空氣轉(zhuǎn)移到溶液，溶液的質(zhì)量源項(xiàng)對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為0*/else if(pa_vps_w)/*如果濕空氣的水蒸氣分壓大于溶液的水蒸氣分壓，則水蒸氣向溶液移動，冷凝放熱，屬于除濕過程*/m_e=*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+)-ps_w)；/* masstransfer rate in dehumidifier */dS=0；else if(pa_vps_w)/*如果濕空氣的水蒸氣分壓小于溶液的水蒸氣分壓，則水蒸氣向濕空氣移動，蒸發(fā)吸熱，屬于再生過程*/m_e=*(pa*mfv/(1-mfv)/(mfv/(