热非平衡效应和切应力耦合作用下波状液膜流的稳定性.pdfVIP

热非平衡效应与切应力耦合作用下 波状液膜流的稳定性1 王松岭，张润盘，叶学民，李春曦 华北电力大学动力工程系，河北省 保定市 071003 E-mail：zhengren_wu@163.com 摘摘 要要：：当液体薄膜流界面同时存在切应力和相变时，切应力和相变引起的非热平衡效应将 摘摘 要要：： 改变液膜的流动状态和界面力平衡状态，二者的相互耦合将使流体动力学特性和稳定性更为 复杂。基于完整边界条件，采用参数摄动法，建立了沿倾斜壁面下降的热非平衡效应与切应 力耦合作用下的二维液膜流动的Orr-Sommerfeld方程，得到波速和扰动增长率的表达式，着 重分析了相变和切应力的影响。结果表明：冷凝状态及其相变强度的增加使得液膜的稳定性 得到加强；而蒸发状态及其相变强度的增加导致液膜稳定性减弱；反向切应力有利于流动稳 定，正向切应力不利于流动稳定；热非平衡效应不仅影响扰动增长率，而且对临界波数和临 界雷诺数也有影响。 关键词关键词：：波动液膜 稳定性 相变 切应力 热非平衡效应 关键词关键词：： 中图分类号中图分类号：：O352 中图分类号中图分类号：： 1．引言．引言 ．．引言引言 液体薄膜流以高传热传质特性在传统工业和高技术领域备受关注，而流动稳定性是影 响其性能的重要因素之一。液膜的稳定性研究始于Benjamin[1]和Yih[2]对等温层流降膜的研 究；Kapitza[3]采用边界层积分法将N-S方程简化为三阶常微分方程(简称Kapitza方程) ，获得 了波动液膜的平均厚度和平均速度。Krantz等[4]对等温降膜采用数值计算方法求解了O-S方 程。之后，线性理论被广泛用于分析液膜流动的稳定性。Nguyen和Balakotaiah[5]提出了对简 化模型的改进模型，但分析结果仅适用于小雷诺数范围。基于边界层理论，Alekseenko等[6-7] [8] [9] 采用双波方程对切应力下的液膜稳定性进行了弱非线性研究。Yu等 和Brauner等 用不同 方法在中、高雷诺数下对等温降膜进行了线性稳定性分析。当液膜的汽液界面存在蒸发或 冷凝等相变过程，或存在界面切应力时，将直接影响液膜的表面流态，从而改变其流动稳 定性及动力学特征。Bankoff等[10]、Spindler等[11]、叶学民等[12-15]对蒸发或冷凝状态下的液膜 表面波的稳定性进行了研究，分析了不同因素的影响规律。Miesen和Bosersma[16]、Aktershev 和Alekseenko[17]、王松岭等[18]分析了切应力作用下波状液膜的动力学特征和流动稳定性，指 明正向和反向切应力的作用相反。研究表明，汽液界面处的相变和切应力对流体动力学特 征和稳定性具有强烈的非线性影响。 当汽液界面同时存在切应力和相变时，二者的相互耦合将使其非线性作用更为复杂化， 相变过程所引起的热非平衡效应不仅改变液膜的流动状态和界面平衡状态，而且界面处蒸发 量或冷凝量的大小对稳定性也有显著的影响[19] 。目前的研究仅限于考虑其中一种因素的影 响。本文通过建立热非平衡效应和切应力相互耦合作用下波状液膜流的Orr-Sommerfeld方程 （简称O-S方程），利用摄动法求得扰动增长率，进而分析在保证液膜不被撕裂的情况下的 热非平衡效应及其相变强度和切应力对流动稳定性的影响规律。 2 ．波状液膜流动模型．波状液膜流动模型 ．．波状液膜流动模型波状液膜流动模型 粘性液体受重力驱动，沿倾斜壁面形成波状下降液膜，液膜表面外是流动蒸汽，汽态分 子冷凝变成液态或液体分子蒸发，通过质量传递在汽液界面处形成热非平衡效应，如图1 所示。对于粘性、不可压缩牛顿流体，波动液膜的控制方程组和连续性方程为 ∂u 1 + u ⋅ ∇u= − ∇p + g+ν∆u (1) ∂t ρ 本课题得到高等学校博士学科点专项科研基金 （项目编号：