基于LevelSet方法对油水和气水两相界面的数值模拟.pdfVIP

27 第2 期 应 用 力 学 学 报 Vol.27 No.2 2010 年6 月 CHINESE JOURNAL OF APPLIED MECHANICS Jun.2010 文章编号 1000- 4939(2010) 02-0298-05 基于Level Set 方法对油水和 气水两相界面的数值模拟 王琳琳 田辉 李国君 （西安交通大学 710049 西安） 摘要：为研究两相界面迁移特性，基于Level Set 方法，建立了求解非定常不可压缩两相界面流动 的数值方法。计算中使用结构化网格采用Level Set 函数捕捉两相界面。通过对经典算例的模拟， 验证了数值方法对界面捕捉的有效性和精确性。模拟了油滴在水中上升、变形，与油层融合的过 程，研究了气泡在产 、发展、脱离阶段的变形机理，和表面张力系数对气泡形状的影响。计算 结果发现，表面张力系数越大，气泡在发展阶段持续的时间就越长，膨胀的程度也越大，并在脱 离时刻，气泡的体积也越大，为进一步研究两相界面迁移特性提供了新的途径。 关键词：Level Set 方法；数值模拟；两相流 中图分类号：TK 121 文献标识码： A 的相函数，并以此相函数为基础重构界面，求解相 1 引 言 函数后经界面重构，得到新的界面位置；并且在保 证一定的精度下，消耗较少的计算资源。在处理三 两相界面流动现象广泛存在于各个领域，如能 维问题，以及对界面的融合、破碎现象的描述时， 源与动力、石油和化工领域等。其具有许多独特的 Level Set 方法显示出比VOF 方法更大的优势。1988 流体力学性质，研究这些性质对流体力学的发展及 年，Osher 和 Sethian 提出了相界面随时间运动的 [1] 工程应用都有重要的意义。由于相间质量、动量、 Level Set 方法 ，其重要思想是引入Level Set 函数 能量交换主要发生在两相界面上，因此对两相流界 f ，将移动的相界面对应为零水平集，只要确定了 面迁移特性的研究格外受到关注，提出了多种数值 零水平集就可以确定界面的位置。这种方法自提出 方法，这些方法主要分为两类，即界面追踪类方法 以来，在界面的演化、流体力学、燃烧、图像处理 和界面捕获类方法。以PIC、MAC 等算法为代表的 [2] 和材料力学等领域都得到了广泛的应用 。 界面追踪类方法是通过拉格朗日法，采用一系列离 本文基于Level Set 方法，建立了一套求解两相 散的标记点，具有较高的精度。但这类方法对于拓 不可压缩流体相界面非定常迁移特性的数值方法。 扑结构变化大的问题处理起来比较困难，而且消耗 在对经典算例验证的基础上，模拟并研究了油滴与 较多的计算资源。以VOF、Level Set 等算法为代表 油层融合和气泡形成、发展、脱离的过程，分析了 的界面捕获类方法是用欧拉法构造表征界面特性 表面张力系数对气泡形状的影响。 来稿日期：2009-03-06 修回日期：20 10-05-25 第一作者简介：王琳琳，女，1981 年生，西安交通大学能源与动力工程学院，博士生；研究方向——多相流流动。 E-mail : linlinwang.724@stu.xjt u.edu.cn 第2 期 王琳琳，等 基于Level Set 方法对油水