A boundary condition capturing method for incompressible flame discontinuities

Abstract In this paper, we propose a new numerical method for treating two-phase incompressible flow where one phase is being converted into the other, e.g., the vaporization of liquid water. We consider this numerical method in the context of treating discontinuously thin flame fronts for incompressible flow. This method was designed as an extension of the Ghost Fluid Method (1999, J. Comput. Phys.152, 457) and relies heavily on the boundary condition capturing technology developed in Liu et al. (2000, J. Comput. Phys.154, 15) for the variable coefficient Poisson equation and in Kang et al. (in press J. Comput. Phys.) for multiphase incompressible flow. Our new numerical method admits a sharp interface representation similar to the method proposed in Helenbrook et al. (1999, J. Comput. Phys.148, 366). Since the interface boundary conditions are handled in a simple and straightforward fashion, the code is very robust, e.g. no special treatment is required to treat the merging of flame fronts. The method is presented in three spatial dimensions, with numerical examples in one, two, and three spatial dimensions.

[1]  H. Markstein.  Nonsteady flame propagation , 1964 .

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

[3]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[4]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

[5]  Forman A. Williams,et al.  The Mathematics of Combustion , 1985, Frontiers in applied mathematics.

[6]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

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

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

[9]  D. Elsworth.  Computational Methods in Fluid Flow , 1993 .

[10]  Charles S. Peskin,et al.  Improved Volume Conservation in the Computation of Flows with Immersed Elastic Boundaries , 1993 .

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

[12]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

[13]  G. Son,et al.  Numerical Simulation of Film Boiling Near Critical Pressures With a Level Set Method , 1998 .

[14]  R. Fedkiw,et al.  The Ghost Fluid Method for de agration and detonation discontinuities , 1998 .

[15]  Gretar Tryggvason,et al.  A Front Tracking Method for the Motion of Premixed Flames , 1998 .

[16]  D. Juric,et al.  Computations of Boiling Flows , 1998 .

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

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

[19]  Brian T. Helenbrook,et al.  A Numerical Method for Solving Incompressible Flow Problems with a Surface of Discontinuity , 1999 .

[20]  Brian T. Helenbrook,et al.  The role of Landau-Darrieus instability in large scale flows , 1999 .

[21]  S. Welch,et al.  A Volume of Fluid Based Method for Fluid Flows with Phase Change , 2000 .

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

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

[24]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .