VOF法与LevelSet法的比较.pdfVIP

VOF法与Level Set法的比较 1 2 汪迎春 ，汪德爟 （1．2．南京河海大学环境学院，南京，210098） 摘 要：VOF法和Level Set法同属于界面捕获类算法，有很多相似之处。本文对两种方法进行了详细的介绍，并 用剪切流场进行了测试，比较分析了两种方法的优缺点。测试的结果表明，VOF法容易破碎，界面重构复杂，但有 些VOF实现方案能够模拟非常尖锐的界面，这一点上要比Level Set 法好。Level Set法不需要重构界面，模拟出 的界面形状要比VOF法光滑，整体效果要比VOF法好。 关键词：Level Set；VOF；界面 1 引言 在计算流体力学中，界面分层流问题是很普遍的。既有水和气强分层的自由面问题，也有冷热水弱 分层的内界面问题。这类问题的计算因为其复杂性成为当前的热点和难点。目前解决此类问题的数值方 法有很多。LevelSet 法与 VOF 法同属于界面捕获类方法，因为它们在欧拉网格下构造，都用一个标量函 数描述界面的几何特性，具有很强的拓扑处理能力，因而是两种很流行的界面流模拟方法。 VOF法在20世纪70 年代末由Hirt和Nichols 等最先提出。VOF法是80 年代国外最流行的方法，目 前国内对 VOF法的研究和应用也越来越多。 水平集方法（Level Set Method）由Osher和Sethian[1]于1988年提出后，引起了人们的广泛的注意。 水平集方法不需要显式地追踪运动界面，避免了对拓扑结构变化的处理，计算稳定，已经应用到诸多领 域，如材料学、流体力学、图象处理、计算机视觉、网格生成等。Level Set法是90年代国外很盛行的方 法，目前国内的研究和应用情况还不如VOF法充分。 2 两种方法介绍 2.1 VOF 方法 VOF方法的基本思想是，在欧拉网格系统上定义一个函数 f(流体体积份额量），根据每个网格内所含 某种物质的体积量来定义了在此网格上的值，然后用体积跟踪的方法求解方程： ∂ + ⋅=∇ 0 f u f ∂t VOF方法采用几何与代数相结合的方法，即在上一时刻重构界面，然后用代数方法算出下一时刻的 f 值。 VOF方法在代数处理阶段常采用贡献网格方法（donor-cell) ，通量修正方法（flux 一 corrected transport , FCT )，迎风方法。界面重构策略可分为分片常数体积跟踪方法和分片线性体积跟踪方法。 分片常数体积跟踪方法主要有两种：一种是 Hirt[2]等的VOF法，另一种是 Noh 和Woodward 的 SLIC 方 ——————————— 作者简介：汪迎春(1976-),男，南京河海大学环境学院，博士研究生 通讯地址：南京市西康路1号 河海大学七舍三单元506 E-MAIL：wangyingch@ 1 [3] 法 。SLIC方法中界面重构在同一时刻要求与某一坐标轴平行。而 Hirt方法是根据界面网格内 f 的大小， 结合相邻网格 f 的值，选择与之适合的坐标轴平行，相对接近于真实界面。分片常数体积跟踪方法重构 界面精度低，易产生非物理界面现象。分片线性体积跟踪法，最早是由Debar提出的，典型的有Youngs方 法[4]和Flair方法。在分片线性体积跟踪方法中 Youngs方法的界面重构是在一个网格内进行，界面不穿过 网格边（面）。在二维流体运动，网格内界面的法向由 9 个网格的 f 的值决定。界面的位置由法向和界 面所在网格的 f 的值及相邻网格的 f 的值确定，并且保证体