压力控制的相变UDF

#includeudf.h#includesg_mphase.h#defineLAT_HT1.e3DEFINE_SOURCE(air_src,cell,pri_th,dS,eqn)/*主相湿空气质量输运源项*/{Thread*mix_th,*sec_th;/*计算区域指针*/realm_v;/*主相质量输运量*/realmfv;/*湿空气中水蒸气的质量分数*/realt_air;/*湿空气的温度*/realhr;/*湿空气的含湿量*/realpa;/*湿空气的总压力*/realpa_v;/*湿空气的水蒸气分压*/realmfs;/*licl溶液的质量分数*/realt_sol;/*licl溶液的温度*/realps_w;/*licl溶液的水蒸气分压*/mix_th=THREAD_SUPER_THREAD(pri_th);/*混合区液相指针*/sec_th=THREAD_SUB_THREAD(mix_th,1);/*单相区气相指针*/mfv=C_YI(cell,pri_th,0);/*获得湿空气中水蒸气的质量分数,注意species列表中的顺序，水是第一个时序号才是0,dS要对该项求导*/hr=mfv/(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,0);/*获得溶液中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.0011080673*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);/*masstransferrateindehumidifier*/dS=0.00165*((pa*(0.378*mfv+0.622)-pa*mfv*0.378)/pow(0.378*mfv+0.622,2));/*湿空气质量源项对湿空气的质量转移求导，对湿空气中水蒸气质量分数mfv求导*/}elseif(pa_vps_w)/*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/{m_v=-0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w);/*masstransferrateindehumidifier*/dS=0;/*水蒸气由溶液向湿空气转移，湿空气质量源项对质量转移的偏导为0*/}returnm_v;}DEFINE_SOURCE(sol_src,cell,sec_th,dS,eqn)/*次相氯化锂溶液质量输运源项*/{Thread*mix_th,*pri_th;realm_s;/*次相溶液相质量输运量*/realmfv;realt_air;/*湿空气的温度*/realhr;/*湿空气的含湿量*/realpa;/*湿空气的总压力*/realpa_v;/*湿空气的水蒸气分压*/realmfs;/*licl溶液的质量分数*/realt_sol;/*licl溶液的温度*/realps_w;/*licl溶液的水蒸气分压*/mix_th=THREAD_SUPER_THREAD(sec_th);/*混合区液相指针*/pri_th=THREAD_SUB_THREAD(mix_th,0);/*单相区气相指针*/mfv=C_YI(cell,pri_th,0);/*获得湿空气中水蒸气的质量分数,注意species列表中的顺序，水是第一个时序号才是0*/hr=mfv/(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,0);/*获得溶液中licl的质量分数,注意species列表中的顺序，licl是第一个时序号才是0,dS要对该项求导*/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_vps_w)/*如果湿空气的水蒸气分压大于溶液的水蒸气分压，则水蒸气向溶液移动，冷凝放热，属于除湿过程*/{m_s=0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w);/*masstransferrateindehumidifier*/dS=0;/*水蒸气由湿空气转移到溶液，溶液的质量源项对质量转移的偏导为0*/}elseif(pa_vps_w)/*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/{m_s=-0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w);/*masstransferrateindehumidifier*/dS=-0.00165*(0.06183172-0.0022161346*mfs+0.0035782429*t_sol);/*dS：溶液的质量源项对质量转移求偏导，对溶液中水蒸气质量分数mfs求导*/}returnm_s;}DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn)/*混合相能量方程输运源项*/{Thread*pri_th,*sec_th;realm_e;/*相变质量*/realmfv;realt_air;/*湿空气的温度*/realhr;/*湿空气的含湿量*/realpa;/*湿空气的总压力*/realpa_v;/*湿空气的水蒸气分压*/realmfs;/*licl溶液的质量分数*/realt_sol;/*licl溶液的温度*/realps_w;/*licl溶液的水蒸气分压*/pri_th=THREAD_SUB_THREAD(mix_th,0);/*混合区液相指针*/sec_th=THREAD_SUB_THREAD(mix_th,1);/*单相区气相指针*/mfv=C_YI(cell,pri_th,0);/*获得湿空气中水蒸气的质量分数,注意species列表中的顺序，水是第一个时序号才是0*/hr=mfv/(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,0);/*获得溶液中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.0011080673*pow(mfs,2)-0.0035782429*t_sol*mfs;/*溶液水蒸气分压*/if(pa_vps_w)/*如果湿空气的水蒸气分压大于溶液的水蒸气分压，则水蒸气向溶液移动，冷凝放热，属于除湿过程*/{m_e=0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w);/*masstransferrateindehumidifier*/dS=0;}elseif(pa_vps_w)/*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/{m_e=-0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w);/*masstransferrateindehumidifier*/dS=0;}returnLAT_HT*m_e;}