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.

A diffuse interface approach for vector-valued PDEs on surfaces

Numerical Methods for Linear Diffusion Equations in the Presence of an Interface

An interface formulation of the Laplace-Beltrami problem on piecewise smooth surfaces

A Framework for Solving Parabolic Partial Differential Equations on Discrete Domains

Diffusion maps for embedded manifolds with boundary with applications to PDEs

A surface finite element method for the Navier-Stokes equations on evolving surfaces

An Efficient Numerical Method for Reaction-Diffusion Equation on the General Curved Surfaces.

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Global Well-Posedness

A high-order diffused-interface approach for two-phase compressible flow simulations using a discontinuous Galerkin framework

Analysis of the Crouzeix-Raviart Surface Finite Element Method for vector-valued Laplacians

A unified structure-preserving parametric finite element method for anisotropic surface diffusion

A High-Order Kernel Method for Diffusion and Reaction-Diffusion Equations on Surfaces

A diffuse-interface method for simulating two-phase flows of complex fluids

Diffuse Interface Model for Two-Phase Flows on Evolving Surfaces with Different Densities: Local Well-Posedness

A meshfree collocation scheme for surface differential operators on point clouds

Diffusion of tangential tensor fields: numerical issues and influence of geometric properties

A diffuse interface method for simulation-based screening of heat transfer processes with complex geometries

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