A 3-D Volume-of-Fluid advection method based on cell-vertex velocities for unstructured meshes

• Published in: Computers & fluids Vol. 94; pp. 14 - 29
• Main Authors: Jofre, Lluís, Lehmkuhl, Oriol, Castro, Jesús, Oliva, Assensi
• Format: Journal Article Publication
• Language: English
• Published: Elsevier Ltd 01-05-2014; Elsevier
• Subjects: 3-D unstructured mesh; Enginyeria mecànica; Fluid mechanics; Interface-capturing; Mecànica de fluids; Volume-of-Fluid; Àrees temàtiques de la UPC; 3-D unstructured mesh; Volume-of-Fluid; Interface-capturing
• ISSN: 0045-7930; 1879-0747
• DOI: 10.1016/j.compfluid.2014.02.001

•Volume-of-Fluid method for 3D unstructured meshes.•Volume-of-Fluid advection method based on cell-vertex velocities.•The new implementation minimizes over/underlapping.•One of the first attempts to report results for unstructured grids.•Results for cartesian meshes are also reported.
A new geometrical Volume-of-Fluid (VOF) method for capturing interfaces on three-dimensional (3-D) Cartesian and unstructured meshes is introduced. The method reconstructs interfaces as first- and second-order piecewise planar approximations (PLIC), and advects volumes in a single unsplit Lagrangian–Eulerian (LE) geometrical algorithm based on constructing flux polyhedrons by tracing back the Lagrangian trajectories of the cell-vertex velocities. In this way, the situations of overlapping between flux polyhedrons are minimized, consequently, the accuracy in the solution of the advection equation is improved by minimizing the creation of overshoots (volume fractions over one), undershoots (volume fractions below zero) and wisps (fluid in void regions or vice versa). However, if not treated carefully, the use of cell-vertex velocities may result in the construction of flux polyhedrons that contain nonplanar faces and that do not conserve volume. Therefore, this work explains in detail a set of geometric algorithms necessary to overcome these two drawbacks. In addition, the new VOF method is analyzed numerically on 3-D Cartesian and unstructured meshes, first, by reconstructing the interface of spherical geometries and, second, by evaluating the final advection result of a sphere placed in a rotation, shear and deformation field.