A multiphase model for compressible flows with interfaces, shocks, detonation waves and cavitation

A compressible multiphase unconditionally hyperbolic model is proposed. It is able to deal with a wide range of applications: interfaces between compressible materials, shock waves in condensed multiphase mixtures, homogeneous two-phase flows (bubbly and droplet flows) and cavitation in liquids. Here we focus on the generalization of the formulation to an arbitrary number of fluids, and to mass and energy transfers, and extend the associated Godunov method. We first detail the specific problems involved in the computation of thermodynamic interface variables when dealing with compressible materials separated by well-defined interfaces. We then address one of the major problems in the modelling of detonation waves in condensed energetic materials and propose a way to suppress the mixture equation of state. We then consider another problem of practical importance related to high-pressure liquid injection and associated cavitating flow. This problem involves the dynamic creation of interfaces. We show that the multiphase model is able to solve these very different problems using a unique formulation. We then develop the Godunov method for this model. We show how the non-conservative terms must be discretized in order to fulfil the interface conditions. Numerical resolution of interface conditions and partial equilibrium multiphase mixtures also requires the introduction of infinite relaxation terms. We propose a way to solve them in the context of an arbitrary number of fluids. This is of particular importance for the development of multimaterial reactive hydrocodes. We finally extend the discretization method in the multidimensional case, and show some results and validations of the model and method.

[1]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[2]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[3]  Jiaquan Gao,et al.  How to prevent pressure oscillations in multicomponent flow calculations , 2000, Proceedings Fourth International Conference/Exhibition on High Performance Computing in the Asia-Pacific Region.

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

[5]  Qiang Zhang,et al.  Three-Dimensional Front Tracking , 1998, SIAM J. Sci. Comput..

[6]  S. SIAMJ.,et al.  A SIMPLE METHOD FOR COMPRESSIBLE MULTIFLUID FLOWS , 1999 .

[7]  C. Farhat,et al.  A method of finite element tearing and interconnecting and its parallel solution algorithm , 1991 .

[8]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[9]  De-kang Mao,et al.  A treatment of discontinuities for finite difference methods in the two-dimensional case , 1993 .

[10]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

[11]  N. Kolev,et al.  Fragmentation and coalescence dynamics in multiphase flows , 1993 .

[12]  M. Baer,et al.  A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials , 1986 .

[13]  R. Abgrall,et al.  A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows , 1999 .

[14]  Keh-Ming Shyue,et al.  An Efficient Shock-Capturing Algorithm for Compressible Multicomponent Problems , 1998 .

[15]  Joseph M. Powers,et al.  Theory of two-phase detonation—Part I: Modeling☆ , 1990 .

[16]  A. Dervieux,et al.  Multifluid incompressible flows by a finite element method , 1981 .

[17]  Richard Saurel,et al.  Numerical analysis of a ram accelerator employing two-phase combustion , 1996 .

[18]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[19]  E. F. Toro,et al.  Riemann-problem-based techniques for computing reactive two-phased flows , 1989 .

[20]  Lionel Sainsaulieu,et al.  Finite volume approximation of two phase-fluid flows based on an approximate Roe-type Riemann solver , 1995 .

[21]  E. Daniel,et al.  Droplet break-up through an oblique shock wave , 1996 .

[22]  F. Williams,et al.  A model for the deflagration of nitramines , 1988 .

[23]  Richard Saurel,et al.  A numerical study of weak shock wave propagation in a reactive bubbly liquid , 1996 .

[24]  D. Benson.  Computational methods in Lagrangian and Eulerian hydrocodes , 1992 .

[25]  Rémi Abgrall,et al.  Efficient numerical approximation of compressible multi-material flow for unstructured meshes , 2003 .

[26]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[27]  Ronald Fedkiw,et al.  An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows , 1999 .

[28]  L. van Wijngaarden,et al.  One-Dimensional Flow of Liquids Containing Small Gas Bubbles , 1972 .

[29]  M. Yuen,et al.  Heat-transfer measurements of evaporating liquid droplets , 1978 .

[30]  P. B. Butler,et al.  Shock Development and Transition to Detonation Initiated by Burning in Porous Propellant Beds , 1981 .

[31]  R. Menikoff,et al.  The Riemann problem for fluid flow of real materials , 1989 .

[32]  Steven F. Son,et al.  Two-phase modeling of DDT: Structure of the velocity-relaxation zone , 1997 .

[33]  Rémi Abgrall,et al.  Proposition de méthodes et modèles eulériens pour les problèmes à interfaces entre fluides compressibles en présence de transfert de chaleur , 2002 .

[34]  Richard Saurel,et al.  A mechanistic model for shock initiation of solid explosives , 1999 .

[35]  S. Bankoff,et al.  Propagation of pressure waves in bubbly mixtures , 1984 .

[36]  Smadar Karni,et al.  Hybrid Multifluid Algorithms , 1996, SIAM J. Sci. Comput..

[37]  Christophe Morel,et al.  The local volumetric interfacial area transport equation: derivation and physical significance , 1999 .

[38]  D. Drew,et al.  Theory of Multicomponent Fluids , 1998 .

[39]  M. Larini,et al.  Ignition and growth of a detonation by a high energy plasma , 1992 .

[40]  D. Bestion,et al.  The physical closure laws in the CATHARE code , 1990 .

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

[42]  J. Edwards,et al.  Ausm schemes and extensions for low mach and multiphase flows , 1999 .

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

[44]  Randall J. LeVeque,et al.  Two-Dimensional Front Tracking Based on High Resolution Wave Propagation Methods , 1996 .

[45]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .