1-s2.0-S0021999118302699-main.pdf (2.55 MB)
0/0

An iterative interface reconstruction method for PLIC in general convex grids as part of a Coupled Level Set Volume of Fluid solver

Download (2.55 MB)
journal contribution
posted on 14.05.2018 by Maciej Skarysz, Andrew Garmory, Mehriar Dianat
Reconstructing the interface within a cell, based on volume fraction and normal direction, is a key part of multiphase flow solvers which make use of piecewise linear interface calculation (PLIC) such as the Coupled Level Set Volume of Fluid (CLSVOF) method. In this paper, we present an iterative method for interface reconstruction (IR) in general convex cells based on tetrahedral decomposition. By splitting the cell into tetrahedra prior to IR the volume of the truncated polyhedron can be calculated much more rapidly than using existing clipping and capping methods. In addition the root finding algorithm is designed to take advantage of the nature of the relationship between volume fraction and interface position by using a combination of Newton's and Muller's methods. In stand-alone tests of the IR algorithm on single cells with up to 20 vertices the proposed method was found to be 2 times faster than an implementation of an existing analytical method, while being easy to implement. It was also found to be 3.4–11.8 times faster than existing iterative methods using clipping and capping and combined with Brent's root finding method. Tests were then carried out of the IR method as part of a CLSVOF solver. For a sphere deformed by a prescribed velocity field the proposed method was found to be up to 33% faster than existing iterative methods. For simulations including the solution of the velocity field the maximum speed up was found to be approximately 52% for a case where 12% of cells lie on the interface. Analysis of the full simulation CPU time budget also indicates that while the proposed method has produced a considerable speed-up, further gains due to increasing the efficiency of the IR method are likely to be small as the IR step now represents only a small proportion of the run time.

Funding

M. Skarysz has been financially supported by a Jaguar Land Rover/UK Engineering and Physical Sciences Research Council (EPSRC) industrial CASE studentship ref 1689977 which is gratefully acknowledged by the authors.

History

School

  • Aeronautical, Automotive, Chemical and Materials Engineering

Department

  • Aeronautical and Automotive Engineering

Published in

Journal of Computational Physics

Citation

SKARYSZ, M., GARMORY, A. and DIANAT, M., 2018. An iterative interface reconstruction method for PLIC in general convex grids as part of a Coupled Level Set Volume of Fluid solver. Journal of Computational Physics, 368, pp.254–276.

Publisher

© Elsevier

Version

VoR (Version of Record)

Publisher statement

This work is made available according to the conditions of the Creative Commons Attribution 4.0 International (CC BY 4.0) licence. Full details of this licence are available at: http://creativecommons.org/licenses/ by/4.0/

Acceptance date

23/04/2018

Publication date

2018

Notes

This is an Open Access Article. It is published by Elsevier under the Creative Commons Attribution 4.0 International Licence (CC BY). Full details of this licence are available at: http://creativecommons.org/licenses/by/4.0/

ISSN

0021-9991

Language

en

Licence

Exports

Logo branding

Licence

Exports