板壁水膜波动流动数值研究 numerical study of water film wave flowing on wall.pdfVIP

第27卷第5期 核动力T程 V01．27．NO．5 2O06 年lO月 NuclearPOwer 0ct．20O6 Engineering 文章编号：0258-0926(2006)05—0029-04 板壁水膜波动流动数值研究 田瑞峰1，李兆俊1，张庆武2 (I．哈尔滨工程大学动力与核能工程学院．15000I；2．哈尔滨工程大学外语系，15000I) 摘要：板壁水膜流动特性是影响水膜二次携带的主要因素之一。通过数值方法对某一规格板壁表面水膜 波动流动特性进行研究。在对水膜表面进行受力分析的基础上，得出了水膜自由表面的动力学边界条件，建 立了板壁水膜波动流动的数学模型。连续性方程采用中心差分格式，动量方程采用迎风和中心差分混合的偏 心差分格式，采用simple算法对板壁水膜波动流动方程进行了求解，采用流体界面追踪法(VOF：v01． ume．of_Fluid)方法对水膜自由表面进行模拟追踪。计算结果表明，水膜在板壁中部的波动比其它部位剧烈， 另外也发现入口扰动频率对水膜波动的影响程度要大于扰动振幅对水膜波动的影响。 关键词：水膜；波动特性；vOF方法；数值模拟 中图分类号：TL7 文献标识码：A 1 引 言 的粘性流体。 流体在波形板汽．水分离器内的流动是复杂 (2)假设在汽．液界面处无相变发生。 的三维两相非定常流动。水滴在板壁上形成水膜， (3)板壁表面水膜在汽流携带力的作用下，沿 汽流横掠竖直平板水膜表面，水膜在重力和汽流 壁面成波状流动状态。 切应力作用下流动。在没有扰动的情况下，水膜 实际中，水膜厚度与流动方向尺度相比非常 保持稳定，以平滑状态流动。实际运行中，由于 小，可以认为水膜的流动为二维流动【l~3】。如图1 大量水滴的入射、汽流的压力、速度的脉动、波 所示。 形板的震动等，都可能引起水膜的波动。在汽流 作用下，水膜波动会发生几个转变。汽流速度足 够小时，汽．水界面处于相对稳定的状态；随着汽 流速度的增加，界面由于Kelvin．Helmholtz不稳 定性而变得不稳定。不稳定的水膜是形成水膜破 裂，造成二次携带的一个主要原因，所以研究板 壁水膜波动流动特性对二次携带有重要意义。 图l水膜流动示意图 Water 本文采用vOF方法对板壁水膜在气流切应 CoordinatesFilm onWall Fig．1 Flowing 力作用下的自由流动进行研究。计算结果表明， x一为水膜流动方向；一为水膜厚度方向 水膜在板壁中部的波动比其它部位剧烈；发现入 2．2基本方程 口扰动频率对水膜波动的影响程度要大于扰动振 (1)水膜流动连续性方程 幅对水膜波动的影响。 坐+坐：0 f1) 舐 瑟 2板壁水膜流动数学模型 (2)水膜流动动量方程 2．1基本假设 d甜 d“ d“ l∞ “．，d一“ d“、 …