Loughborough University
Leicestershire, UK
LE11 3TU
+44 (0)1509 263171
Loughborough University

Loughborough University Institutional Repository

Please use this identifier to cite or link to this item: https://dspace.lboro.ac.uk/2134/32948

Title: An iterative interface reconstruction method for PLIC in general convex grids as part of a Coupled Level Set Volume of Fluid solver
Authors: Skarysz, Maciej
Garmory, Andrew
Dianat, Mehriar
Keywords: Interface reconstruction
PLIC
VOF
CLSVOF
Root-finding
Issue Date: 2018
Publisher: © Elsevier
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.
Abstract: 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.
Description: 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/
Sponsor: 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.
Version: Published
DOI: 10.1016/j.jcp.2018.04.044
URI: https://dspace.lboro.ac.uk/2134/32948
Publisher Link: https://doi.org/10.1016/j.jcp.2018.04.044
ISSN: 0021-9991
Appears in Collections:Published Articles (Aeronautical and Automotive Engineering)

Files associated with this item:

File Description SizeFormat
1-s2.0-S0021999118302699-main.pdfPublished version2.61 MBAdobe PDFView/Open

 

SFX Query

Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.