Accurate and efficient surface reconstruction from volume fraction data on general meshes

Abstract

Simulations involving free surfaces and fluid interfaces are important in many areas of engineering. There is, however, still a need for improved simulation methods. Recently, a new efficient geometric VOF method called isoAdvector for general polyhedral meshes was published. We investigate the interface reconstruction step of isoAdvector, and demonstrate that especially for unstructured meshes the applied isosurface based approach can lead to noisy interface orientations. We then introduce a novel computational interface reconstruction scheme based on calculation of a reconstructed distance function (RDF). By iterating over the RDF calculation and interface reconstruction, we obtain second order convergence of both the interface normal and position within cells even with a strict $L_{\infty }$ error norm. In 2D this is verified with reconstruction of a circle on Cartesian meshes and on unstructured triangular and polygonal prism meshes. In 3D the second order convergence is verified with reconstruction of a sphere on Cartesian meshes and on unstructured tetrahedral and polyhedral meshes. The new scheme is combined with the interface advection step of the isoAdvector algorithm. Significantly reduced absolute advection errors are obtained, and for CFL number 0.2 and below we demonstrate second order convergence on all the mentioned mesh types in 2D and 3D. The implementation of the proposed interface reconstruction schemes is straightforward and the computational cost is significantly reduced compared to contemporary methods. The schemes are implemented as an extension to the Computational Fluid Dynamics (CFD) Open Source software package, OpenFOAM. The extension module and all test cases presented in this paper are released as open source.

References in corpus (1)

• A Computational Method for Sharp Interface Advection

Cited by in corpus (20)

• Unstructured un-split geometrical Volume-of-Fluid methods -- A review
• varRhoTurbVOF: a new set of volume of fluid solvers for turbulent isothermal multiphase flows in OpenFOAM
• THINC scaling method that bridges VOF and level set schemes
• gVOF: An open-source package for unsplit geometric volume of fluid methods on arbitrary grids
• Deep learning of interfacial curvature: a symmetry-preserving approach for the volume of fluid method
• Influence of Cavity Geometry on the Bubble Dynamics of Nucleate Pool Boiling
• Computing volume fractions and signed distances from triangulated surfaces immersed in unstructured meshes
• Simulation and analytical modeling of high-speed droplet impact onto a surface
• An unstructured finite-volume level set / front tracking method for two-phase flows with large density-ratios
• Stabilizing the unstructured Volume-of-Fluid method for capillary flows in microstructures using artificial viscosity
• Inconsistencies in Unstructured Geometric Volume-of-Fluid Methods for Two-Phase Flows with High Density Ratios
• Numerical wetting benchmarks -- advancing the plicRDF-isoAdvector unstructured Volume-of-Fluid (VOF) method
• TwoPhaseFlow: An OpenFOAM based framework for development of two phase flow solvers
• Experimental and Numerical Study of Microcavity Filling Regimes for Lab-on-a-Chip Applications
• An unstructured geometrical un-split VOF method for viscoelastic two-phase flows
• Second-order accurate normal reconstruction from volume fractions on unstructured meshes with arbitrary polyhedral cells
• THINC-scaling scheme that unifies VOF and level set methods
• The influence of volumetric shrinkage on the metal solidification process under localized energy deposition
• A Deep Learning Algorithm for Piecewise Linear Interface Construction (PLIC)
• cfdmfFTFoam: A front-tracking solver for multiphase flows on general unstructured grids in OpenFOAM