《压力控制的相变UDF》由会员分享，可在线阅读，更多相关《压力控制的相变UDF（6页珍藏版）》请在金锄头文库上搜索。
1、#i nclude udf.h#in elude sg_mphase.h #defi ne LAT_HT 1.e3DEFINE_SOURCE(air_src,cell,pri_th,dS,eq n) /* 主相湿空气质量输运源项 */ Thread *mix th,*sec th;/* 计算区域指针 */real m v;/*主相质量输运量*/real mfv;/*湿空气中水蒸气的质量分数*/real t air;/*湿空气的温度*/real hr; /*湿空气的含湿量*/real pa;/*湿空气的总压力*/real pa v; /*湿空气的水蒸气分压*/real mfs; /*licl溶液
2、的质量分数*/real t sol; /*licl 溶液的温度 */real ps_w; /*licl溶液的水蒸气分压*/mix th=THREAD SUPER THREAD(pri th);/* 混合区液相指针*/sec_th=THREAD_SUB_THREAD(mix_th,1);/* 单相区气相指针*/mfv=C Yl(cell,pri th,0);/*获得湿空气中水蒸气的质量分数，注意species列表中的顺序，水是第一个时序号才是0,dS要对该项求导*/hr=mfv/(1-mfv);/*由湿空气中水蒸气的质量分数算到湿空气的含湿量*/pa=ABS P(C P(cell,pri th)
3、,op pres); /* 获取湿空气的总压力 */ pa v=(pa*hr)/(hr+0.622);/*获取湿空气水蒸气分压 */* t air=C T(cell,pri th); /*湿空气的温度，单位为摄氏度*mfs=C YI(cell,sec th,O);/*获得溶液中licl的质量分数，注意species列表中的顺序，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.001
4、1080673*pow(mfs,2)-0.0035782429*t_sol*mfs; /* 溶液水蒸气分压 */if(pa_vps_w) /*如果湿空气的水蒸气分压大于溶液的水蒸气分压，则水蒸气向溶液移动，冷凝放热，属于除湿过程*/m_v=0.00165*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+0.622)-ps_w);/* mass transfer ratein dehumidifier */dS=0.00165*(pa*(0.378*mfv+0.622)-pa*mfv*0.378)/pow(0.378*mfv+0.622,2); /* 湿空 气质量源项对湿空气的质量转移
5、求导，对湿空气中水蒸气质量分数mfv 求导*/else if(pa_vps w)/*如果湿空气的水蒸气分压大于溶液的水蒸气分压，则水蒸气向溶液移动，冷凝放热，属于除湿过程*/m_s=0.00165*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+0.622)-ps_w);/* mass tran sfer rate indehumidifier */else if(pa_vvps_w) /*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/ m_s=-0.00165*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+0.622)-ps
6、_w);/* mass transfer ratein dehumidifier */dS=-0.00165*(0.06183172-0.0022161346*mfs+0.0035782429*t_sol); /*dS :溶液的质量源项对质量转移求偏导，对溶液中水蒸气质量分数mfs求导*/return m_s;DEFINE_SOURCE(e nrg_src,cell,mix_th,dS,eq n) /* 混合相能量方程输运源项 */Thread *pri th,*sec th;real m_e;/*相变质量*/real mfv;real t air;/*湿空气的温度*/real hr; /*湿
7、空气的含湿量*/real pa;/*湿空气的总压力*/real pa v; /*湿空气的水蒸气分压*/real mfs; /*licl溶液的质量分数*/real t sol; /*licl 溶液的温度 */real ps_w; /*licl溶液的水蒸气分压*/pri th=THREAD SUB THREAD(mix th,0);/* 混合区液相指针*/sec_th=THREAD_SUB_THREAD(mix_th,1);/* 单相区气相指针*/mfv=C_YI(cell,pri_th,O);/*获得湿空气中水蒸气的质量分数，注意species列表中的顺序，水是第一个时序号才是0*/hr=mfv
8、/(1-mfv);/*由湿空气中水蒸气的质量分数算到湿空气的含湿量*/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,O);/*获得溶液中licl的质量分数，注意species列表中的顺序，licl是第一个时序号才是 0*/t sol=C T(cell,sec th); /*溶液的温度，单位为摄氏度*/ps_w=1.1491997+0.11
9、242725*t_sol-0.06183172*mfs+0.0014521543*pow(t_sol,2)+0.0011080673*pow(mfs,2)-0.0035782429*t_sol*mfs; /* 溶液水蒸气分压 */if(pa vps w)/*如果湿空气的水蒸气分压大于溶液的水蒸气分压，则水蒸气向溶液移动，冷|琵放热，属于除湿过程*/、匕M凝 m_e=0.00165*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+0.622)-ps_w);/* mass transfer ratein dehumidifier */ dS=O;else if(pa_vps_w) /*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/m_e=-0.00165*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+0.622)-ps_w);/* mass transfer ratein dehumidifier*/dS=0;return LAT_HT*m_e;
