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

Henning Scheufler*, Johan Roenby

*Kontaktforfatter

Publikation: Bidrag til tidsskriftTidsskriftartikelForskningpeer review

2 Citationer (Scopus)

Resumé

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.

OriginalsprogEngelsk
TidsskriftJournal of Computational Physics
Vol/bind383
Sider (fra-til)1-23
Antal sider23
ISSN0021-9991
DOI
StatusUdgivet - 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

Citer dette

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

I: Journal of Computational Physics, Bind 383, 15.04.2019, s. 1-23.

Publikation: Bidrag til tidsskriftTidsskriftartikelForskningpeer 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

VL - 383

SP - 1

EP - 23

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -