二维数值模拟蒸汽膜覆盖的液滴水力学特性.pdfVIP

314 中国科学 E 辑 工程科学 材料科学 2005, 35(3): 314~323 * 段日强 姜胜耀 张佑杰 贾海军 杨星团 ( , 102201) 清华大学核能与新能源技术研究院教育部先进反应堆工程与安全重点实验室 北京 摘要 采用移动粒子半隐式法(moving particle semi-implicit method, 简称MPS 法)直接数值模拟了蒸汽膜覆盖的液滴破裂, 求解了关于液滴、蒸汽和环境液体的 二维 方程 计算结果表明蒸汽膜的存在抑制了液滴破裂 液滴破裂 Navier-Stokes . , 的临界Weber 数随蒸汽膜厚度增加而增加. 液滴破裂过程可分为两个阶段, 液滴 变形与破裂主要发生于后一个阶段. 同时计算结果揭示, 除了液滴与蒸汽界面处 的Rayleigh-Taylor 不稳定性较强外, 蒸汽膜覆盖的液滴破裂机理与无蒸汽膜覆盖 的液滴破裂机理是一致的. 数值模拟结果的部分特征得到了相关实验的验证. 关键词 液滴破裂 蒸汽膜 数值模拟 粒子法 移动粒子半隐式法 蒸汽爆炸, 其特征表现为蒸汽快速大量的产生. 在轻水反应堆堆芯熔化严重 事故中, 由于熔融金属的温度远高于冷却水的饱和沸腾温度, 高温熔融金属与冷 却水的直接接触导致了蒸汽爆炸这一重要物理现象的发生. 当发生堆芯熔化事 故后 熔融金属以射流的形式射入冷却水中 首先破裂变成为相对较大的液滴 , , , 由于液滴与冷却水的接触面积有限此时尚不足以触发蒸汽爆炸, 触发蒸汽爆炸 的诱因是由于以上生成的相对大的熔融金属液滴的第二次细粒化而导致熔融金 属与冷却水的界面积大量快速增加. 关于熔融金属液滴第二次细粒化过程, 过去 许多研究者为了简化分析通常不考虑冷却水沸腾所产生的蒸汽影响而单纯地研 究无蒸汽膜覆盖的液滴破裂特性[1], 但合理的假设应该是熔融金属液滴被冷却水 沸腾所产生的蒸汽膜所覆盖. 关于无蒸汽膜覆盖的液滴破裂物理机理, 先前的研究者[2,3]对其已进行了非 常广泛的调查并已有了相当深入的了解, 但是关于蒸汽膜覆盖的液滴破裂机理, 由于以下几方面的原因, 到现在为止对其认识还非常有限: 从实验的角度看, 由 于蒸汽膜的存在导致了光学观察的模糊化以及甚至被阻挡, 因而难以得到清晰 2004-06-21 收稿, 2004-12-09 收修改稿 * ( : 国家自然科学基金资助项目 批准号 SCIENCE IN CHINA Ser. E Engineering Materials Science 第 3 期 段日强等: 二维数值模拟蒸汽膜覆盖的液滴水力学特性 315 的液滴破裂图像 而在数值模拟方面 由于该问题涉及到大界面变形的多流体多 ; , 相流动,