压力控制的相变UDF.docxVIP

下载本文档

113
0
约4.83千字
约 4页
2017-08-23 发布于重庆
举报
版权申诉

压力控制的相变UDF.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

压力控制的相变UDF

#include udf.h#include sg_mphase.h#define LAT_HT 1.e3DEFINE_SOURCE(air_src,cell,pri_th,dS,eqn)? ? /*主相湿空气质量输运源项*/{? ? Thread *mix_th,*sec_th;? ?? ?? ?? ?? ?? ? /*计算区域指针*/? ? real m_v;? ?? ?? ?? ?? ?? ?? ???/*主相质量输运量*/? ? real mfv;? ?? ?? ?? ?? ?? ?? ?? ?? ? /*湿空气中水蒸气的质量分数*/? ? real t_air;? ???/*湿空气的温度*/? ? real hr;? ? /*湿空气的含湿量*/? ? real pa;? ? /*湿空气的总压力*/? ? real pa_v;??/*湿空气的水蒸气分压*/? ? real mfs;??/*licl溶液的质量分数*/? ? 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_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+0t_sol-0mfs+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);? ?/* 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));??/*湿空气质量源项对湿空气的质量转移求导，对湿空气中水蒸气质量分数mfv求导*/? ?? ? }? ? else if(pa_vps_w)? ?/*如果湿空气的水蒸气分压小于溶液的水蒸气分压，则水蒸气向湿空气移动，蒸发吸热，属于再生过程*/? ?? ? {? ?? ?? ? m_v=-0.00165*(((pa*mfv/(1-mfv))/(mfv/(1-mfv)+0.622))-ps_w);? ?/* mass transfer rate in dehumidifier */? ?? ?? ? dS=0;??/* 水蒸气由溶液向湿空气转移，湿空气质量源项对质量转移的偏导为0*/? ?? ?? ???}? ?return m_v;}DEFINE_SOURCE(sol_src,cell,sec_th,dS,eqn)? ? /*次相氯化锂溶液质量输运源项*/{? ? Thread *mix_th,*pri_th;? ? real m_s;? ? /*次相溶液相质量输运量*/? ? real mfv;? ? real t_air;? ???/*湿空气的温度*/? ? real hr;? ? /*湿空气的含湿量*/? ? real pa;? ? /*湿空气的总压力*/? ? real pa_v;??/*湿空气的水蒸