Single-step reinitialization and extending algorithms for level-set based multi-phase flow simulations

Abstract We propose efficient single-step formulations for reinitialization and extending algorithms, which are critical components of level-set based interface-tracking methods. The level-set field is reinitialized with a single-step (non iterative) “forward tracing” algorithm. A minimum set of cells is defined that describes the interface, and reinitialization employs only data from these cells. Fluid states are extrapolated or extended across the interface by a single-step “backward tracing” algorithm. Both algorithms, which are motivated by analogy to ray-tracing, avoid multiple block-boundary data exchanges that are inevitable for iterative reinitialization and extending approaches within a parallel-computing environment. The single-step algorithms are combined with a multi-resolution conservative sharp-interface method and validated by a wide range of benchmark test cases. We demonstrate that the proposed reinitialization method achieves second-order accuracy in conserving the volume of each phase. The interface location is invariant to reapplication of the single-step reinitialization. Generally, we observe smaller absolute errors than for standard iterative reinitialization on the same grid. The computational efficiency is higher than for the standard and typical high-order iterative reinitialization methods. We observe a 2- to 6-times efficiency improvement over the standard method for serial execution. The proposed single-step extending algorithm, which is commonly employed for assigning data to ghost cells with ghost-fluid or conservative interface interaction methods, shows about 10-times efficiency improvement over the standard method while maintaining same accuracy. Despite their simplicity, the proposed algorithms offer an efficient and robust alternative to iterative reinitialization and extending methods for level-set based multi-phase simulations.

[1]  Nikolaus A. Adams,et al.  Adaptive multi-resolution method for compressible multi-phase flows with sharp interface model and pyramid data structure , 2014, J. Comput. Phys..

[2]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[3]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[4]  Nikolaus A. Adams,et al.  Wavelet-based adaptive multi-resolution solver on heterogeneous parallel architecture for computational fluid dynamics , 2011, Computer Science - Research and Development.

[5]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[6]  J. Haas,et al.  Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities , 1987, Journal of Fluid Mechanics.

[7]  Nikolaus A. Adams,et al.  Anti-diffusion method for interface steepening in two-phase incompressible flow , 2011, J. Comput. Phys..

[8]  Nikolaus A. Adams,et al.  A conservative sharp interface method for incompressible multiphase flows , 2015, J. Comput. Phys..

[9]  P. Koumoutsakos,et al.  A Lagrangian particle level set method. , 2005 .

[10]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[11]  W. Q. Tao,et al.  A coupled volume-of-fluid and level set (VOSET) method for computing incompressible two-phase flows , 2010 .

[12]  Nikolaus A. Adams,et al.  Scale separation for multi-scale modeling of free-surface and two-phase flows with the conservative sharp interface method , 2015, J. Comput. Phys..

[13]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[14]  Xiangyu Hu,et al.  An interface interaction method for compressible multifluids , 2004 .

[15]  H. Choi,et al.  A direct reinitialization approach of level‐set/splitting finite element method for simulating incompressible two‐phase flows , 2011 .

[16]  Theo G. Theofanous,et al.  High-fidelity interface tracking in compressible flows: Unlimited anchored adaptive level set , 2007, J. Comput. Phys..

[17]  Oriol Lehmkuhl,et al.  A coupled volume-of-fluid/level-set method for simulation of two-phase flows on unstructured meshes , 2016 .

[18]  Bo Yu,et al.  A coupled volume-of-fluid and level set (VOSET) method based on remapping algorithm for unstructured triangular grids , 2017 .

[19]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[20]  Nikolaus A. Adams,et al.  On the HLLC Riemann solver for interface interaction in compressible multi-fluid flow , 2009, J. Comput. Phys..

[21]  P. Smereka,et al.  A Remark on Computing Distance Functions , 2000 .

[22]  Theo G. Theofanous,et al.  Direct numerical simulation of interfacial instabilities: A consistent, conservative, all-speed, sharp-interface method , 2013, J. Comput. Phys..

[23]  Frank Losasso,et al.  A fast and accurate semi-Lagrangian particle level set method , 2005 .

[24]  Elias Balaras,et al.  Topology preserving advection of implicit interfaces on Cartesian grids , 2015, J. Comput. Phys..

[25]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

[26]  T. Aslam A partial differential equation approach to multidimensional extrapolation , 2004 .

[27]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[28]  S. Zalesak Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[29]  H. Martin Bücker,et al.  Parallel re-initialization of level set functions on distributed unstructured tetrahedral grids , 2011, J. Comput. Phys..

[30]  W. Dehnen,et al.  Improving convergence in smoothed particle hydrodynamics simulations without pairing instability , 2012, 1204.2471.

[31]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[32]  James J. Quirk,et al.  On the dynamics of a shock–bubble interaction , 1994, Journal of Fluid Mechanics.

[33]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[34]  Tai Wang,et al.  A coupled volume-of-fluid and level set (VOSET) method on dynamically adaptive quadtree grids , 2013 .

[35]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[36]  Theo G. Theofanous,et al.  Adaptive characteristics-based matching for compressible multifluid dynamics , 2006, J. Comput. Phys..

[37]  Nikolaus A. Adams,et al.  A conservative interface method for compressible flows , 2006, J. Comput. Phys..

[38]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[39]  Frederick Stern,et al.  A new volume-of-fluid method with a constructed distance function on general structured grids , 2012, J. Comput. Phys..

[40]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[41]  E. Rouy,et al.  A viscosity solutions approach to shape-from-shading , 1992 .

[42]  Kai Schneider,et al.  An adaptive multiresolution scheme with local time stepping for evolutionary PDEs , 2008, J. Comput. Phys..