VOF方法中自由面重构算法的研究.pdf

VOF 方法中自由面重构算法的研究 班长英 河海大学理学院，（210098 ） E-mail:banying_2004@163.com 摘要：针对处理自由面问题的流体体积函数（VOF）法提出了一种新的自由面重构的算法， 该方法是在单个网格内用斜线段近似自由面，根据边界网格的流体体积分数及边界线与网格 线的交点坐标列出对应于每个边界网格的方程，通过求解方程组即可准确构造出自由面的具 体位置和形状。给出了具体算例并用不同的方法进行数值模拟，结果表明：该方法明显优于 其他方法，它能更好的保持流体的界面形状，证实了新方法的有效性和可行性。 关键词：VOF 方法；自由面重构；数值模拟 A New Free Surface Reconstruction Method in VOF Ban Chang-ying (College of Sciences, Hohai University, Nanjin g 210098, China) Abstract: A two-dimensional free surface reconstruction approach for VOF(the volume of fluid) method was presented in this paper．It was using sloped line segment to reconstruct free surface in a boundary mesh.The equations were listed out according to the fractional volume of fluid and the intersection of the borderline and mesh.Then the surface was reconstructed by solving the equations.This paper also provided results of some sample models.The computational results indicated that this method was more efficient than other methods for reconstructing free surface. Key word:VOF method;free surface reconstruction；numerical simulation 1．引言 当前，运动界面问题的数值计算和模拟在 自然科学和工程计算中成为越来越热门的课 题，并承担着越来越重要的任务，原因是它涉及到了各种理论和技术领域的数值研究和应用， 例如机械、能源、电子、动力、水利和建筑[1−3]等。由于运动界面的复杂性和整体位置的不 断变化，所以精确模拟运动界面的难度非常大。目前处理运动界面问题比较理想的方法是 VOF(volume of fluid)法，它以存储量小、应用简单、追踪的锐利性好、精细而得到广泛应用。 VOF 方法的基本原理是通过研究网格单元中流体和网格体积比函数 F 来确定自由面， 追踪流体的变化。在一个网格单元内，若 F=1,则说明该网格充满流体；若 F=0,则该网格为 空网格；若 0F1，则该网格为自由面网格。由于 VOF 法追踪的是网格中的流体体积，而 非追踪自由面流体质点的运动（如 MAC 法，边界元法），因此可以处理自由面重入等强非 线性现象。 VOF 自由面追踪方法包括两部分研究内容：根据给定流速及重构好的 自由液面确定下 一时刻 F 值的自由面输运算法和根据给定的 F 值确定自由液面位置的 自由面重构算法； 本文仅讨论 自由面重构算法。施主受主法 [4 ]是利用水平和垂直自由液面单元近似模拟自由 面，当水平相邻自由面网格的界面倾斜角都大于 45 度时，网格间将产生非物理的空隙，为 了使VOF 法能够精确地定量