水蒸气蒸发UDF
《水蒸气蒸发UDF》由会员分享,可在线阅读,更多相关《水蒸气蒸发UDF(3页珍藏版)》请在装配图网上搜索。
1、对飞昂公司水一水蒸汽蒸发过程UDF程序段的解释#iiiclude udf.h包括常规宏#iiiclude sg_mphase.li包括体积分数宏 CVOF(C,T)define T_SAT 373 ,定义蒸发温度 100Cdefine LAT_HT l.e3定义蒸发潜热 J,KgDEFINE_SOURCE(liq_sic, cell, pn_th, dS, eqn)液相质量源项 UDFTluead *nux_th. *sec_th;,定义计算区线指针real m_dot_l;定义液相质量转移kg/(in2.s)mix_th = THREAD_SUPER_THREAD(pri_th);指向混合区
2、的主相即液相的指针sec_th = THREAD_SUB_THREAD(mix_th, 1);指向单相控制区的气相的指针,气相为第二 相if(C_T(cell, inix_th)=T_SAT)如果液相单元的温度高于蒸发温度,液相向气相的质量 质量转移= -0*C_VOF(cell, pri_th)*C_R(cell, pn_th)*I混合区液相温度-饱和温度I 0.1X1/被相XQ液相3fhbs(C_T(cell. mix_th) - T_SAT)/T_SAT;饱和温度mdot =dSeqn = -O.l*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T
3、_SAT)门T_SAT;定义源项对质量转移偏导s4)=a + b(/else (= 0*C_VOF(celL sec_th)*C_R(cell. sec_th)* fhbs(T_SAT-C_T(celLmix_th)/T_SAT;如果指向混合区液相的单元温度小于蒸发温度,气相向液相的质量转移,液相得mdot = 0.1 x V,相 x q 气相I混合区液相温度-饱和温度IX 饱和温度dSeqn = 0.;由于是气相向液相转移,所以液相的质量源项对质量转移的偏导为零return m dot 1;DEFINE_SOURCE(vap_src, cell. sec_th, dS, eqn)气相质量源项
4、 UDFTluead * mix_th, *pri_th;real m dot v;mix_th = THREAD_SUPER_THREAD(sec_th);指向混合区的第二相即气相的指针pn_th = THREAD_SUB_THREAD(mix_th, 0);指向单相控制区的液相的指针,液相为主相if(C_T(cell, niix_th)=T_SAT)如果混合区单元的温度高于蒸发温度,液相向气相的质量 质量转移m_do顷=0.1*C_VOF(cell, pri_th)*C_R(cell. pii_th)*fhbs(C_T(cell. mix_th) - T_SAT)/T_SAT;dSeqn
5、= 0.;由于是液相向气相转移,所以气相的质量源项对来自液相的质量转移的 偏导为零1else (m_dot_v = -0*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fhbs(T_SAT-C_T(celLmix_th)/T_SAT;如果指向混合区的单元温度小于蒸发温度,气相向液相的质量转移,气相失mdot = -O.lx 吃和 x p倾I混合区液相温度-饱和温度IX饱和温度dSeqn = -0.l*C_R(cell, sec_th)* fhbs(C_T(cell, mix_th) - T_SAT)/T_SAT;由于是气相向液相转移,所以气相的质量源项对自身的质量
6、转移的偏导不为零 1return m dot v;DEFINE_SOURCE(enrg*rc, cell, mix_th, dS, eqn)混合模型能量源项 UDFTluead *pri_th, *sec_th;real m_dot;pn_th = THREAD_SUB_THREAD(nux_tli, 0);指向混合区的液相的指针sec_th = THREAD_SUB_THREAD(nux_th. 1);指向混合区的气相的指针if(C_T(celL nux_th)=T_SAT)如果混合区的单元温度高于蒸发温度。质量转移由液相向 气相转移,吸热,质量转移量前有负号m_dot = -0*C_VOF
7、(cell, pri_th)*C_R(cell, pri_th)*fhbs(C_T(cell. mix_th) - T_SAT)T_SAT;dSeqn = -0.1 *C_VOF(cHl, pri_th)*C_R(cHl, pri_th)/T_SAT;else (相反,气相向液相转移则放热m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fhbs(T_SAT-C_T(celLmix_th)/T_SAT;dSeqn = 0.l*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;return LAT_HT *m_dot;(化潜热与质量转移率相乘得能量源项W/ni3
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。