On improving mass conservation properties of the hybrid particle-level-set method

Abstract For the computation of multi-phase flows level-set methods are an attractive alternative to volume-of-fluid or front-tracking approaches. For improving their accuracy and efficiency the hybrid particle-level-set modification was proposed by Enright et al. [Enright D, Fedkiw R, Ferziger J, Mitchell I. A hybrid particle-level-set method for improved interface capturing. J Comput Phys 2002;183:83–116]. In actual applications the overall properties of a level-set method, such as mass conservation, are strongly affected by discretization schemes and algorithmic details. In this paper we address these issues with the objective of determining the optimum alternatives for the purpose of direct numerical simulation of dispersed-droplet flows. We evaluate different discretization schemes for curvature and unit normal vector at the interface. Another issue is the particular formulation of the reinitialization of the level-set function which significantly affects the quality of computational results. Different approaches employing higher-order schemes for discretization, supplemented either by a correction step using marker particles (Enright et al., 2002) or by additional constraints [Sussman M, Almgren AS, Bell JB, Colella P, Howell LH, Welcome ML. An adaptive level set approach for incompressible two-phase flows. J Comput Phys 1999;148:81–124] are analyzed. Different parameter choices for the hybrid particle-level-set method are evaluated with the purpose of increasing the efficiency of the method. Aiming at large-scale computations we find that in comparison with pure level-set methods the hybrid particle-level-set method exhibits better mass-conservation properties, especially in the case of marginally resolved interfaces.

[1]  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..

[2]  Martin E. Weber,et al.  Bubbles in viscous liquids: shapes, wakes and velocities , 1981, Journal of Fluid Mechanics.

[3]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[4]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[5]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

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

[7]  David J. Torres,et al.  The point-set method: front-tracking without connectivity , 2000 .

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

[9]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

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

[11]  P. Colella,et al.  An Adaptive Level Set Approach for Incompressible Two-Phase Flows , 1997 .

[12]  Hendric C. Kuhlmann,et al.  Free surface flows , 1998 .

[13]  S. Osher,et al.  Level set methods: an overview and some recent results , 2001 .

[14]  Mark Sussman,et al.  A parallelized, adaptive algorithm for multiphase flows in general geometries , 2005 .

[15]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[16]  Oleg V. Vasilyev.  High Order Finite Difference Schemes on Non-uniform Meshes with Good Conservation Properties , 2000 .

[17]  W. Rider,et al.  Stretching and tearing interface tracking methods , 1995 .

[18]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[19]  D. Kothe.  Perspective on Eulerian Finite Volume Methods for Incompressible Interfacial Flows , 1998 .

[20]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[21]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[22]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

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

[24]  R. Fedkiw,et al.  USING THE PARTICLE LEVEL SET METHOD AND A SECOND ORDER ACCURATE PRESSURE BOUNDARY CONDITION FOR FREE SURFACE FLOWS , 2003 .

[25]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[26]  B. Engquist,et al.  Numerical approximations of singular source terms in differential equations , 2004 .

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

[28]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[29]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

[30]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .