A high-resolution Godunov method for compressible multi-material flow on overlapping grids

A numerical method is described for inviscid, compressible, multi-material flow in two space dimensions. The flow is governed by the multi-material Euler equations with a general mixture equation of state. Composite overlapping grids are used to handle complex flow geometry and block-structured adaptive mesh refinement (AMR) is used to locally increase grid resolution near shocks and material interfaces. The discretization of the governing equations is based on a high-resolution Godunov method, but includes an energy correction designed to suppress numerical errors that develop near a material interface for standard, conservative shock-capturing schemes. The energy correction is constructed based on a uniform-pressure-velocity flow and is significant only near the captured interface. A variety of two-material flows are presented to verify the accuracy of the numerical approach and to illustrate its use. These flows assume an equation of state for the mixture based on the Jones-Wilkins-Lee (JWL) forms for the components. This equation of state includes a mixture of ideal gases as a special case. Flow problems considered include unsteady one-dimensional shock-interface collision, steady interaction of a planar interface and an oblique shock, planar shock interaction with a collection of gas-filled cylindrical inhomogeneities, and the impulsive motion of the two-component mixture in a rigid cylindrical vessel.

[1]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[2]  R. Abgrall.  How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

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

[4]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[5]  J. O. Hallquist,et al.  Modeling two-dimensional shock initiation and detonation-wave phenomena in PBX 9404 and LX-17 , 1981 .

[6]  Rémi Abgrall,et al.  A Simple Method for Compressible Multifluid Flows , 1999, SIAM J. Sci. Comput..

[7]  Smadar Karni,et al.  Viscous shock profiles and primitive formulations , 1992 .

[8]  R B Pember,et al.  A Dynamically Adaptive Arbitrary Lagrangian-Eulerian Method for Solution of the Euler Equations , 2003 .

[9]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[10]  W. Henshaw,et al.  An adaptive numerical scheme for high-speed reactive flow on overlapping grids , 2003 .

[11]  J. O. Hallquist,et al.  Modeling short pulse duration shock initiation of solid explosives , 1985 .

[12]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[13]  William D. Henshaw,et al.  A Scheme for Conservative Interpolation on Overlapping Grids , 1994, SIAM J. Sci. Comput..

[14]  W. Henshaw,et al.  Composite overlapping meshes for the solution of partial differential equations , 1990 .

[15]  Paul Glaister,et al.  An approximate linearised Riemann solver for the Euler equations for real gases , 1988 .

[16]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[17]  Joseph M. Powers,et al.  Two-phase viscous modeling of compaction of granular materials , 2004 .

[18]  Donald W. Schwendeman,et al.  Numerical shock propagation in non-uniform media , 1988 .

[19]  H. S. Udaykumar,et al.  An Eulerian method for computation of multimaterial impact with ENO shock-capturing and sharp interfaces , 2003 .

[20]  Patrick Jenny,et al.  Correction of Conservative Euler Solvers for Gas Mixtures , 1997 .

[21]  Keh-Ming Shyue,et al.  A wave-propagation based volume tracking method for compressible multicomponent flow in two space dimensions , 2006, J. Comput. Phys..

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

[23]  Smadar Karni,et al.  Multicomponent Flow Calculations by a Consistent Primitive Algorithm , 1994 .

[24]  C. M. Tarver,et al.  Phenomenological model of shock initiation in heterogeneous explosives , 1980 .

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

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

[27]  M. Berger.  ON CONSERVATION AT GRID INTERFACES. , 1987 .

[28]  B. M. Dobratz,et al.  Properties of chemical explosives and explosive simulants , 1972 .

[29]  Tariq D. Aslam,et al.  NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS , 2002 .

[30]  William D. Henshaw,et al.  Moving overlapping grids with adaptive mesh refinement for high-speed reactive and non-reactive flow , 2006, J. Comput. Phys..

[31]  Tariq D. Aslam,et al.  High resolution numerical simulation of ideal and non-ideal compressible reacting flows with embedded internal boundaries , 1997 .