A study of planar Richtmyer-Meshkov instability in fluids with Mie-Grüneisen equations of state

We present a numerical comparison study of planar Richtmyer-Meshkov instability with the intention of exposing the role of the equation of state. Results for Richtmyer-Meshkov instability in fluids with Mie-Gruneisen equations of state derived from a linear shock-particle speed Hugoniot relationship (Jeanloz, J. Geophys. Res. 94, 5873, 1989; McQueen et al., High Velocity Impact Phenomena (1970), pp. 294–417; Menikoff and Plohr, Rev. Mod. Phys. 61(1), 75 1989) are compared to those from perfect gases under nondimensionally matched initial conditions at room temperature and pressure. The study was performed using Caltech’s Adaptive Mesh Refinement, Object-oriented C++ (AMROC) (Deiterding, Adaptive Mesh Refinement: Theory and Applications (2005), Vol. 41, pp. 361–372; Deiterding, “Parallel adaptive simulation of multi-dimensional detonation structures,” Ph.D. thesis (Brandenburgische Technische Universitat Cottbus, September 2003)) framework with a low-dissipation, hybrid, center-difference, limiter patch solver (Ward and Pullin, J. Comput. Phys. 229, 2999 (2010)). Results for single and triple mode planar Richtmyer-Meshkov instability when a reflected shock wave occurs are first examined for mid-ocean ridge basalt (MORB) and molybdenum modeled by Mie-Gruneisen equations of state. The single mode case is examined for incident shock Mach numbers of 1.5 and 2.5. The planar triple mode case is studied using a single incident Mach number of 2.5 with initial corrugation wavenumbers related by k_1 = k_2+k_3. Comparison is then drawn to Richtmyer-Meshkov instability in perfect gases with matched nondimensional pressure jump across the incident shock, post-shock Atwood ratio, post-shock amplitude-to-wavelength ratio, and time nondimensionalized by Richtmyer’s linear growth time constant prediction. Differences in start-up time and growth rate oscillations are observed across equations of state. Growth rate oscillation frequency is seen to correlate directly to the oscillation frequency for the transmitted and reflected shocks. For the single mode cases, further comparison is given for vorticity distribution and corrugation centerline shortly after shock interaction. Additionally, we examine single mode Richtmyer-Meshkov instability when a reflected expansion wave is present for incident Mach numbers of 1.5 and 2.5. Comparison to perfect gas solutions in such cases yields a higher degree of similarity in start-up time and growth rate oscillations. The formation of incipient weak waves in the heavy fluid driven by waves emanating from the perturbed transmitted shock is observed when an expansion wave is reflected.

[1]  D. Layzer,et al.  On the Instability of Superposed Fluids in a Gravitational Field. , 1955 .

[2]  I. Amdur,et al.  Kinetic Theory of Gases , 1959 .

[3]  R. D. Richtmyer.  Taylor instability in shock acceleration of compressible fluids , 1960 .

[4]  E. Meshkov.  Instability of the interface of two gases accelerated by a shock wave , 1969 .

[5]  J. N. Fritz,et al.  CHAPTER VII – THE EQUATION OF STATE OF SOLIDS FROM SHOCK WAVE STUDIES , 1970 .

[6]  A. J. Cable,et al.  High-velocity impact phenomena , 1970 .

[7]  K. Meyer,et al.  Numerical Investigation of the Stability of a Shock‐Accelerated Interface between Two Fluids , 1972 .

[8]  P. A. Thompson,et al.  Compressible‐Fluid Dynamics , 1973 .

[9]  David H. Rudy,et al.  A nonreflecting outflow boundary condition for subsonic navier-stokes calculations , 1980 .

[10]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

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

[12]  R. Jeanloz.  Shock wave equation of state and finite strain theory , 1989 .

[13]  K. Thompson.  Time-dependent boundary conditions for hyperbolic systems, II , 1990 .

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

[15]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .

[16]  Yang,et al.  Quantitative theory of Richtmyer-Meshkov instability. , 1993, Physical review letters.

[17]  Ravi Samtaney,et al.  Circulation deposition on shock-accelerated planar and curved density-stratified interfaces: models and scaling laws , 1994, Journal of Fluid Mechanics.

[18]  Qiang Zhang,et al.  Small amplitude theory of Richtmyer–Meshkov instability , 1994 .

[19]  Bradford Sturtevant,et al.  Experiments on the Richtmyer-Meshkov instability of an air/SF6 interface , 1995 .

[20]  Ravi Samtaney,et al.  On initial‐value and self‐similar solutions of the compressible Euler equations , 1996 .

[21]  J. Jacobs,et al.  Experimental study of incompressible Richtmyer–Meshkov instability , 1996 .

[22]  E. Puckett,et al.  A High-Order Godunov Method for Multiple Condensed Phases , 1996 .

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

[24]  Qiang Zhang,et al.  AN ANALYTICAL NONLINEAR THEORY OF RICHTMYER-MESHKOV INSTABILITY , 1996 .

[25]  Hazak.  Lagrangian formalism for the Rayleigh-Taylor instability. , 1996, Physical review letters.

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

[27]  M. A. Jones,et al.  A membraneless experiment for the study of Richtmyer–Meshkov instability of a shock-accelerated gas interface , 1997 .

[28]  K. Nishihara,et al.  Asymptotic growth in the linear Richtmyer–Meshkov instability , 1997 .

[29]  Qiang Zhang.  Analytical Solutions of Layzer-Type Approach to Unstable Interfacial Fluid Mixing , 1998 .

[30]  Björn Sjögreen,et al.  The Convergence Rate of Finite Difference Schemes in the Presence of Shocks , 1998 .

[31]  Marilyn Schneider,et al.  Richtmyer–Meshkov instability growth: experiment, simulation and theory , 1999, Journal of Fluid Mechanics.

[32]  L. Chambers.  Linear and Nonlinear Waves , 2000, The Mathematical Gazette.

[33]  Keh-Ming Shyue,et al.  A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Grüneisen equation of state , 2001 .

[34]  J. G. Wouchuk.  Growth rate of the Richtmyer–Meshkov instability when a rarefaction is reflected , 2001 .

[35]  J. G. Wouchuk.  Growth rate of the linear Richtmyer-Meshkov instability when a shock is reflected. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[36]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[37]  S. Abarzhi.  A new type of the evolution of the bubble front in the Richtmyer–Meshkov instability , 2002 .

[38]  V N Goncharov,et al.  Analytical model of nonlinear, single-mode, classical Rayleigh-Taylor instability at arbitrary Atwood numbers. , 2002, Physical review letters.

[39]  J. Jacobs,et al.  PLIF flow visualization and measurements of the Richtmyer–Meshkov instability of an air/SF6 interface , 2002, Journal of Fluid Mechanics.

[40]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[41]  Grégoire Allaire,et al.  A five-equation model for the simulation of interfaces between compressible fluids , 2002 .

[42]  Sung-Ik Sohn.  Simple potential-flow model of Rayleigh-Taylor and Richtmyer-Meshkov instabilities for all density ratios. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  Ralf Deiterding,et al.  Parallel adaptive simulation of multi-dimensional detonation structures , 2003 .

[44]  J. Jacobs,et al.  Experimental study of the Richtmyer–Meshkov instability of incompressible fluids , 2003, Journal of Fluid Mechanics.

[45]  Julian C. Cummings,et al.  A Virtual Test Facility for the Simulation of Dynamic Response in Materials , 2002, The Journal of Supercomputing.

[46]  Ralf Deiterding,et al.  Construction and Application of an AMR Algorithm for Distributed Memory Computers , 2005 .

[47]  J. Jacobs,et al.  Experiments on the late-time development of single-mode Richtmyer–Meshkov instability , 2005 .

[48]  Oleg Schilling,et al.  Effects of WENO flux reconstruction order and spatial resolution on reshocked two-dimensional Richtmyer-Meshkov instability , 2006, J. Comput. Phys..

[49]  Ralf Deiterding,et al.  A low numerical dissipation patch-based adaptive mesh refinement method for large-eddy simulation of compressible flows , 2007, J. Comput. Phys..

[50]  D. Pullin,et al.  Startup process in the Richtmyer-Meshkov instability , 2007 .

[51]  Britton J. Olson,et al.  Rayleigh-Taylor shock waves , 2007 .

[52]  V. Gregory Weirs,et al.  Adaptive Mesh Refinement - Theory and Applications , 2008 .

[53]  D. Pullin,et al.  Shock-resolved Navier–Stokes simulation of the Richtmyer–Meshkov instability start-up at a light–heavy interface , 2009, Journal of Fluid Mechanics.

[54]  G. M. Ward,et al.  A hybrid, center-difference, limiter method for simulations of compressible multicomponent flows with Mie-Grüneisen equation of state , 2010, J. Comput. Phys..