一种基于非结构化网格的VOF方法.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：096040 一种基于非结构化网格的VOF方法 黄萌，陈斌* （西安交通大学动力工程多相流国家重点实验室，西安，710049 ） （Tel：029 E-mail: chenbin@mail.xjtu.edu.cn） 摘要：本文基于三角形非结构化网格，发展出一种类似于结构化网格中 SLIC （Simple Line Interface Calculation ）方法的界面重构手段。该方法为纯欧拉方法，根据位于两相界面网格的相邻网格流体信息 在该网格内部构造一条平行于某边的界面线，通过界面线与网格边界围成的几何形状计算流体输运， 采用修正系数保证了全局和局部的质量守恒。使用均匀平直速度场中的同心双运动界面、旋转速度场 中的马蹄形界面和剪切流场中的圆形运动界面模型等三个经典算例检测了方法的效果，其中前两个算 例结果较令人满意，但剪切运动的处理结果还不太理想，必须进一步发展更高精度的界面重构方式。 关键词：两相流，数值模拟，非结构化网格，VOF 0 前言 研究气液/液液两相流体时，必须考虑相界面在流体运动中的形变和位移。现有的界 [1] 面捕获方法可以分为 Front Tracking 法和 Front Capturing 法两类 。Front Tracking 法在 Euler 网格上进行离散，在网格中布置若干标记点（Marker ）或流体质点（Particle ），通 过对这些标记点或质点的 Lagrangian 追踪来模拟界面的变化。利用这种方法可以显式追 踪界面，从而更精确地计算界面的曲率，但是界面的重构非常复杂，难以处理界面拓扑 结构发生变化的情形。Front capturing 方法包括 VOF （Volume of Fluid ）方法[2]和 Level set 方法[3]，这类方法是完全的 Eulerian 方法，可以比较容易地处理复杂的物质界面以及界 面拓扑结构发生变化的情形（如合并、交叉、破碎等）。同Level set 方法相比，VOF 方 法通过求解关于 VOF 函数的标量输运方程将复杂的界面处理为内部的运动边界，无需 特殊处理即可达到处理拓扑结构变化的鲁棒性和计算的简单高效。 VOF 方法的概念最早由 Hirt 和 Nichols[4]提出，采用流体体积比函数 F 来表示每个 网格中指定相流体所占的体积。VOF 函数随流体运动发生改变，根据每个网格对应的F 函数值可以判断两相界面的位置并通过某种方式重构两相界面。他们提出的施主受主原 理很大程度上影响了后来其他研究者的改进方法，其基本思想直到现在仍然很受重视， 后人主要从方程差分格式和自由面流体传输等方面着手，结合自由面重构格式对 VOF 方法进行改进，发展出了 SLIC[5]、FLAIR[6]、PLIC[7]等各种界面重构方法。到现在为止， 结构化网格上的 VOF 方法已经相当成熟，如由 Youngs 发展的 PLIC 方法可以构造出清 *基金项目：国家自然科学基金项目）；教育部新世纪优秀人才计划项目（NCET-07-0661 ）； 教育部科学技术研究重点项目（107101） 晰准确的两相界面。 然而结构化网格的适用范围有限。当计算区域不规则的时候，非结构化网格可以适 应各种形状的计算区域，发展基于非结构化网格的 VOF 方法，可以有效模拟复杂计算 区域中的多相流动，比结构化网格具有更加普遍的适用性。当前国际国内针对非结构化 网格上的VOF 方法刚刚开展初步研究，一般采用欧拉-拉格朗日相结合的界面重构方式， 如 JI Chun-ning et. al[8]用这一方法成功进行了简单涡和复杂变形场的计算，Xiaofeng Yang et. al[9]详细描述了