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

Henning Scheufler*, Johan Roenby

*Corresponding author

Research output: Contribution to journalJournal articleResearchpeer-review

3 Citations (Scopus)

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 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.

Original languageEnglish
JournalJournal of Computational Physics
Volume383
Pages (from-to)1-23
Number of pages23
ISSN0021-9991
DOIs
Publication statusPublished - 15 Apr 2019

Fingerprint

Surface Reconstruction
Surface reconstruction
Advection
Volume Fraction
mesh
Volume fraction
Mesh
Prisms
Software packages
Computational fluid dynamics
Fluids
Distance Function
advection
Cartesian
Costs
VOF Method
Isosurface
Open Source Software
Unstructured Mesh
Prism

Keywords

  • Interface reconstruction
  • isoAdvector
  • Multiphase flow
  • Reconstructed Distance Function
  • Unstructured meshes
  • Volume of fluid

Cite this

@article{03f1a189c2a044cb9b8448f10f86144f,
title = "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∞ 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{\circledR}. The extension module and all test cases presented in this paper are released as open source.",
keywords = "Interface reconstruction, isoAdvector, Multiphase flow, Reconstructed Distance Function, Unstructured meshes, Volume of fluid",
author = "Henning Scheufler and Johan Roenby",
year = "2019",
month = "4",
day = "15",
doi = "10.1016/j.jcp.2019.01.009",
language = "English",
volume = "383",
pages = "1--23",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press",

}

Accurate and efficient surface reconstruction from volume fraction data on general meshes. / Scheufler, Henning; Roenby, Johan.

In: Journal of Computational Physics, Vol. 383, 15.04.2019, p. 1-23.

Research output: Contribution to journalJournal articleResearchpeer-review

TY - JOUR

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

AU - Scheufler, Henning

AU - Roenby, Johan

PY - 2019/4/15

Y1 - 2019/4/15

N2 - 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∞ 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.

AB - 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∞ 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.

KW - Interface reconstruction

KW - isoAdvector

KW - Multiphase flow

KW - Reconstructed Distance Function

KW - Unstructured meshes

KW - Volume of fluid

UR - http://www.scopus.com/inward/record.url?scp=85061001917&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2019.01.009

DO - 10.1016/j.jcp.2019.01.009

M3 - Journal article

AN - SCOPUS:85061001917

VL - 383

SP - 1

EP - 23

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -