An accurate conservative level set/ghost fluid method for simulating turbulent atomization

This paper presents a novel methodology for simulating incompressible two-phase flows by combining an improved version of the conservative level set technique introduced in [E. Olsson, G. Kreiss, A conservative level set method for two phase flow, J. Comput. Phys. 210 (2005) 225-246] with a ghost fluid approach. By employing a hyperbolic tangent level set function that is transported and re-initialized using fully conservative numerical schemes, mass conservation issues that are known to affect level set methods are greatly reduced. In order to improve the accuracy of the conservative level set method, high order numerical schemes are used. The overall robustness of the numerical approach is increased by computing the interface normals from a signed distance function reconstructed from the hyperbolic tangent level set by a fast marching method. The convergence of the curvature calculation is ensured by using a least squares reconstruction. The ghost fluid technique provides a way of handling the interfacial forces and large density jumps associated with two-phase flows with good accuracy, while avoiding artificial spreading of the interface. Since the proposed approach relies on partial differential equations, its implementation is straightforward in all coordinate systems, and it benefits from high parallel efficiency. The robustness and efficiency of the approach is further improved by using implicit schemes for the interface transport and re-initialization equations, as well as for the momentum solver. The performance of the method is assessed through both classical level set transport tests and simple two-phase flow examples including topology changes. It is then applied to simulate turbulent atomization of a liquid Diesel jet at Re=3000. The conservation errors associated with the accurate conservative level set technique are shown to remain small even for this complex case.

[1]  M. Herrmann.  A Domain Decomposition Parallelization of the Fast Marching Method , 2003 .

[2]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[3]  D. Chopp.  Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[4]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[5]  Robert D. Falgout,et al.  hypre: A Library of High Performance Preconditioners , 2002, International Conference on Computational Science.

[6]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

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

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

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

[10]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

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

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

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

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

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

[16]  C. Pierce,et al.  Progress-variable approach for large-eddy simulation of turbulent combustion , 2001 .

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

[18]  C. Weber.  Zum Zerfall eines Flüssigkeitsstrahles , 1931 .

[19]  Marcus Herrmann,et al.  A balanced force refined level set grid method for two-phase flows on unstructured flow solver grids , 2008, J. Comput. Phys..

[20]  Heinz Pitsch,et al.  High order conservative finite difference scheme for variable density low Mach number turbulent flows , 2007, J. Comput. Phys..

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

[22]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[23]  P. Moin,et al.  Effects of the Computational Time Step on Numerical Solutions of Turbulent Flow , 1994 .

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

[25]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

[26]  T. Ménard,et al.  Coupling level set/VOF/ghost fluid methods: Validation and application to 3D simulation of the primary break-up of a liquid jet , 2007 .

[27]  Stéphane Popinet,et al.  A front-tracking algorithm for accurate representation of surface tension , 1999 .

[28]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[29]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics , 2007, J. Comput. Phys..

[30]  Gunilla Kreiss,et al.  A conservative level set method for two phase flow II , 2005, J. Comput. Phys..

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

[32]  Andrea Prosperetti,et al.  Motion of two superposed viscous fluids , 1981 .