竖直管内降膜流动气液两相运动数值模拟 numerical simulation on gas-liquid two-phase motion of falling film in vertical tube.pdfVIP

■ 第7卷第1期 热 科 学 与 技 术 V01．7NO．1 2008年3月 ofThermalScienceand Mar．2008 Journal Technology 文章编号：1671-8097(2008)01—0005-06 竖直管内降膜流动气液两相运动数值模拟 刘 升， 郝英立 (东南大学能源与环境学院洁净煤发电及燃烧技术教育部重点实验室，江苏南京 210096) 摘要：对较高流速下强迫竖直降膜进行了模拟。建立了描述竖直管内强迫降膜的物理模型和数学模型，采用 RNGk-e模型描述管内气体和液体的复杂湍流流动过程．采用VOF方法对气液相界面进行追踪。通过观察流 动过程中的液膜前端速度矢量变化、整个流场中的压力分布变化过程。对降膜前端的“托举”现象进行了分 析．对竖直管内液膜形成过程中的两相复杂运动进行了分析研究，指出了避免。托举”现泉出现的条件． 关键词：降膜；湍流；速度矢量；VOF方法；托举 中图分类号：TKl24文献标识码：A 高条件下，对竖直管内液膜从形成到稳定这一非 0 引 言 定常强迫降膜流动的数值模拟方法进行了探索， 降膜流动是在化工、能源领域常见的流动形 对相界面的波动和两相中的速度分布和压力分 式，研究和掌握其运动规律可以有效地组织相关 布，尤其是其液膜前端区域两相速度分布进行了 的过程操作，有利于生产的安全高效。众多学者对 详细的分析研究。 降膜的流动特征和换热机理从实验研究和数值计 1 数学模型 算两个方面进行了大量卓有成效的研究工作。 本文以三维大直径圆管内无相互渗透气一液 Liu和Gollub[1咭]通过实验研究了相界面的 of 稳定性和相界面波动产生和变化规律；Sekine和两相顺流降膜为研究对象，采用VOF(volume KumadaE3]通过实验观测了液膜的流动特性；fluid)方法，追踪气液相界面。图l给出了模拟对 象及其边界条件示意图。 Seban和FaghriE4]则研究了表面波动对输运特性 的影响。 水流进口 由于气液间自由界面的存在及其气液问复杂 的相互作用导致的随机变形、波动等复杂运动现 象，对降膜的建模和求解变得非常困难。自Bach 和Villadsen[53成功模拟了相界面上的波动形成 和液膜中的流场分布以来，降膜的流动、传热和传 质特性得到了广泛的关注，但是已有的模拟研究 工作大部分为液相雷诺数较低的大空间、稳定降 膜流动(或流动的稳定段)，采用周期性边界条件 或给定初始速度分布的方法来简化模拟。部分学 者对液相雷诺数较高的情况进行了模拟，但是其 在模拟中对气相仍然采用大空间条件处理，不考 图l 物理模型和边界条件 虑气相速度，这与工程实际情况尤其是动力工程 model condition Fig．1Physicala粤dboundary 中的实际应用有较大差别。本文在