压力控制的相变UDF

1、#include #include #define LAT_HTDEFINE_SOURCE(air_src,cell,pri_th,dS,eqn)Thread *mix_th,*sec_th;realreal质量分数*/realrealrealrealrealrealreal/*主相湿空气质量输运源项*/m_v;mfv;t_air;hr;pa;pa_v;mfs;t_sol;ps_w；/*湿空气的温度*/*湿空气的含湿量*/*湿空气的总压力*/*湿空气的水蒸气分压*/*licl溶液的质量分数*/*licl溶液的温度*/*licl溶液的水蒸气分压*/*计算区域指针*/*主相质量输运量*/*湿空气中

2、水蒸气的/*混合区液相指针*/*单相区气相指针*/*获得湿空气中水蒸气的质量分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+;/*获取湿空气水蒸气分压*/* t_ai

3、r=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=+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);/*湿空气质量源项对湿空气的质量转移求导，对湿空气中水蒸气质量分数mfv求导*/else if(pa_vps_w)/*如果湿空气的水蒸气分压大于溶液的水蒸气分压，则水蒸气向溶液移动，冷凝放热，属于除湿过程*/m_s=*(pa*mfv/(1-mfv)/(mfv/(1-mfv)+)-ps_w)；/* masstransfer rate in dehumidifier */dS=0；/*水蒸气由湿空气转移到溶液，溶液的质量源项对质量转移的偏导为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)/*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/m_e=*(pa*mfv/(1-mfv)/(mfv/(