非热平衡效应与应力耦合作用下冷凝液膜动力特性.pdfVIP

非热平衡效应与应力耦合作用下冷凝液膜非热平衡效应与应力耦合作用下冷凝液膜1 非热平衡效应与应力耦合作用下冷凝液膜非热平衡效应与应力耦合作用下冷凝液膜 的动力特性的动力特性 的动力特性的动力特性 王松岭，张润盘，叶学民，马辉 华北电力大学电站设备状态监测与控制教育部重点实验室，河北省保定市(071003) E-mail：zhengren_wu@163.com 摘摘 要要：：基于完整边界条件，采用参数摄动法，建立了沿倾斜壁面下降的非热平衡效应与切 摘摘 要要：： 应力耦合作用下的二维液膜流动的Orr-Sommerfeld 方程，得到波速和扰动增长率的表达式， 着重分析了相变和切应力的影响。结果表明：冷凝状态及其相变强度的增加使得液膜的稳定 性得到加强；而蒸发状态及其相变强度的增加导致液膜稳定性减弱；反向切应力有利于流动 稳定，正向切应力不利于流动稳定；非热平衡效应不仅影响扰动增长率，而且对临界波数和 临界雷诺数也有影响。 关键词关键词：：液膜 冷凝 质量跃变流率 表面波 关键词关键词：： 中图分类号中图分类号：：O352 中图分类号中图分类号：： 1．引言．引言 ．．引言引言 液体薄膜流的稳定性直接关系到其高传热传质特性在传统工业和高技术领域的应用。液 膜的稳定性分析始于Benjamin[1]和Yih[2]对等温层流降膜的研究。Kapitza[3]采用边界层积分方 法将Navier-Stokes方程简化为3阶常微分方程(后称此为Kapitza方程) ，获得了液膜在波动情况 下的平均厚度和平均速度。之后，线性理论曾一度被用于分析液膜流动的稳定性[4] 。基于边 界层理论，Alekseenko [5,6]采用双波方程对切应力下的液膜稳定性进行了弱非线性研究。Yu[7] 和Brauner[8]用不同的方法在不同雷诺数范围内对等温降膜进行了线性稳定性分析。 Bank0ff[9]、Spindler[10] 等对蒸发或冷凝状态下的液膜表面波进行了研究，并分析了不同因素 [11] 的影响。顾汉洋 通过在连续液相流动中注入气体产生单气泡，采用双平行探针技术和摄 像方法对汽液流中气泡的形态特征进行了实验研究，系统考察了液体流速、管道倾角以及气 泡长度对气泡形态特征的影响。叶学民[12,13]推导了在蒸发、等温和冷凝状态下普遍适用的线 性二维液膜表面波扰动空间演化方程，分析了有界面切应力作用下的永久破断特性。而考虑 冷凝和切应力同时存在情况下的液膜稳定性分析，几乎没有文献记载。 本文直接从包含丰富流动特性的Navier-Stokes方程着手进行理论建模，基于完整的边界 条件，推导出包含Orr-Sommerfeld方程，利用参数摄动法得到冷凝和切应力同时存在情况下 波速C 的表达式，并分析了雷诺数、波数、韦伯数、质量跃变流率、倾角等参数对液膜动力 特性的影响。 2 ．冷凝液膜流动模型．冷凝液膜流动模型 ．．冷凝液膜流动模型冷凝液膜流动模型 粘性液体受重力驱动，沿倾斜壁面形成波状的下降液膜，液体自由表面附近是流动蒸汽。 由于非热平衡效应的影响下，汽态分子会冷凝下来变成液态，通过质量传递进入液膜流动 中(如图1)。根据实际情况，液膜的流动可认为是二维流动。对于粘性、不可压缩牛顿流体， 波动液膜的控制方程组和连续性方程为： 本课题得到高等学校博士学科点专项科研基金 （项目编号：20040079004 ）的资助。 -1- ∂u 1 + u ⋅ ∇u= − ∇p + g+ν∆u (1)