Research Papers: Fundamental Issues and Canonical Flows

Observation of the Development of Secondary Features in a Richtmyer–Meshkov Instability Driven Flow

[+] Author and Article Information
Tennille Bernard

Department of Mechanical Engineering,
University of New Mexico,
Albuquerque, NM 87131
e-mail: tenncb10@unm.edu

C. Randall Truman

Department of Mechanical Engineering,
University of New Mexico,
Albuquerque, NM 87131
e-mail: truman@unm.edu

Peter Vorobieff

Department of Mechanical Engineering,
University of New Mexico,
Albuquerque, NM 87131
e-mail: kalmoth@unm.edu

Clint Corbin

Department of Mechanical Engineering,
University of New Mexico,
Albuquerque, NM 87131
e-mail: clcorbin@unm.edu

Patrick J. Wayne

Department of Mechanical Engineering,
University of New Mexico,
Albuquerque, NM 87131
e-mail: pwayne@unm.edu

Garrett Kuehner

Department of Mechanical Engineering,
University of New Mexico,
Albuquerque, NM 87131
e-mail: garrett.kuehner@gmail.com

Michael Anderson

Illinois Rocstar, LLC,
Champaign, IL 61826
e-mail: mjanderson@illinoisrocstar.com

Sanjay Kumar

Associate Professor
Department of Mechanical Engineering,
University of Texas, Brownsville,
TX 78520
e-mail: sanjay.kumar@utb.edu

From here on, the word “spike” in the text will refer to this feature, which is driven by shock focusing. It is comprised of heavy gas entering the surrounding light gas, but distinct in its physical origin from the RMI “spike” according to traditional terminology.

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received December 31, 2013; final manuscript received June 5, 2014; published online September 10, 2014. Assoc. Editor: Samuel Paolucci.

J. Fluids Eng 137(1), 011206 (Sep 10, 2014) (6 pages) Paper No: FE-13-1759; doi: 10.1115/1.4027829 History: Received December 31, 2013; Revised June 05, 2014

Richtmyer–Meshkov instability (RMI) has long been the subject of interest for analytical, numerical, and experimental studies. In comparing results of experiment with numerics, it is important to understand the limitations of experimental techniques inherent in the chosen method(s) of data acquisition. We discuss results of an experiment where a laminar, gravity-driven column of heavy gas is injected into surrounding light gas and accelerated by a planar shock. A popular and well-studied method of flow visualization (using glycol droplet tracers) does not produce a flow pattern that matches the numerical model of the same conditions, while revealing the primary feature of the flow developing after shock acceleration: the pair of counter-rotating vortex columns. However, visualization using fluorescent gaseous tracer confirms the presence of features suggested by the numerics; in particular, a central spike formed due to shock focusing in the heavy-gas column. Moreover, the streamwise growth rate of the spike appears to exhibit the same scaling with Mach number as that of the counter-rotating vortex pair (CRVP).

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Richtmyer, R. D., 1960, “Taylor Instability in Shock Acceleration of Compressible Fluids,” Commun. Pure Appl. Math., 13(2), pp. 297–319. [CrossRef]
Meshkov, E. E., 1969, “Instability of the Interface of Two Gases Accelerated by a Shock Wave,” Izv. Akad. Nauk SSSR, Mekh. Zhidk. Gaza, 4(5), pp. 151–157. [CrossRef]
Strutt, J. W., 1883, “Investigation of the Character of the Equilibrium of an Incompressible Heavy Fluid of Variable Density,” Proc. London Math. Soc., 14, pp. 170–177. [CrossRef]
Taylor, G., 1950, “The Instability of Liquid Surfaces When Accelerated in a Direction Perpendicular to their Planes,” Proc. R. Soc. London, Ser. A, 201(1065), pp. 192–196. [CrossRef]
Lewis, D., 1950, “The Instability of Liquid Surfaces When Accelerated in a Direction Perpendicular to their Planes. Part II,” Proc. R. Soc. London, Ser. A, 202(1068), pp. 81–96. [CrossRef]
Zhang, Q., and Sohn, S.-I., 1997, “Nonlinear Theory of Unstable Fluid Mixing Driven by Shock Wave,” Phys. Fluids, 9, pp. 1106–1124. [CrossRef]
Saffman, P., and Meiron, D., 1989, “Kinetic Energy Generated by the Incompressible Richtmyer–Meshkov Instability in a Continuously Stratified Fluid,” Phys. Fluids A, 1, pp. 1767–1771. [CrossRef]
Jacobs, J., and Sheeley, J., 1996, “Experimental Study of Incompressible Richtmyer–Meshkov Instability,” Phys. Fluids, 8, pp. 405–415. [CrossRef]
Brouillette, M., 2002, “The Richtmyer–Meshkov Instability,” Annu. Rev. Fluid Mech., 34, pp. 445–468. [CrossRef]
Kane, J., Drake, R., and Remington, B., 1999, “An Evaluation of the Richtmyer–Meshkov Instability in Supernova Remnant Formation,” Astrophys. J., 511(1), pp. 335–340. [CrossRef]
Wu, C., and Roberts, P., 1999, “Richtmyer–Meshkov Instability and the Dynamics of the Magnetosphere,” Geophys. Res. Lett., 26(6), pp. 655–658. [CrossRef]
Goncharov, V., 1999, “Theory of the Ablative Richtmyer–Meshkov Instability,” Phys. Rev. Lett., 82, pp. 2091–2094. [CrossRef]
Fishbine, B., 2002, “Code Validation Experiments,” Los Alamos Research Quarterly, pp. 6–14.
Jacobs, J., Jones, M., and Niederhaus, C., 1995, “Experimental Studies of Richtmyer–Meshkov Instability,” 5th International Workshop on Compressible Turbulent Mixing, Stony Brook, NY, July 18–21, pp. 18–21.
Prestridge, K., Vorobieff, P., Rightley, P., and Benjamin, R., 2000, “Validation of an Instability Growth Model Using Particle Image Velocimetry Measurements,” Phys. Rev. Lett., 84(19), pp. 4353–4356. [CrossRef] [PubMed]
Jacobs, J., Klein, D., Jenkins, D., and Benjamin, R., 1993, “Instability Growth Patterns of a Shock-Accelerated Thin Fluid Layer,” Phys. Rev. Lett., 70(5), pp. 583–586. [CrossRef] [PubMed]
Budzinski, J. M., Benjamin, R. F., and Jacobs, J. W., 1994, “Influence of Initial Conditions on the Flow Patterns of a Shock-Accelerated Thin Fluid Layer,” Phys. Fluids, 6, p. 3510. [CrossRef]
Vetter, M., and Sturtevant, B., 1995, “Experiments on the Richtmyer–Meshkov Instability of an Air/SF6 Interface,” Shock Waves, 4(5), pp. 247–252. [CrossRef]
Anderson, M., 2011, “Oblique Shock Interactions With Perturbed Density Interfaces,” Ph.D. dissertation, The University of New Mexico, Albuquerque, NM.
Anderson, M., Vorobieff, P., Conroy, J., Truman, C., White, R., and Kumar, S., 2014, “An Experimental and Numerical Study of Shock Interaction With a Gas Column Seeded With Droplets,” Shock Waves (submitted).
Crepeau, J., Needham, C. E., and Hikida, S., 2001, “Second Order Hydrodynamic Automatic Mesh Refinement Code (SHAMRC): Volume I, Methodology,” Applied Research Associates, Inc.
Crepeau, J., Happ, H., Hikida, S., and Needham, C. E., 2001, “Second Order Hydrodynamic Automatic Mesh Refinement Code (SHAMRC): Volume II, User's Manual,” Applied Research Associates, Inc.
Johnson, E., 2009, “Planar and Oblique Shock Wave Interaction With a Droplet Seeded Gas Cylinder,” Master's thesis, The University of New Mexico, Albuquerque, NM.
Vorobieff, P., Anderson, M., Conroy, J., White, R., Truman, C. R., and Kumar, S., 2011, “Vortex Formation in a Shock-Accelerated Gas Induced by Particle Seeding,” Phys. Rev. Lett., 106(18), p. 184503. [CrossRef] [PubMed]
Vorobieff, P., Tomkins, C., Kumar, S., Goodenough, C., Mohamed, N., and Benjamin, R., 2004, “Secondary Instabilities in Shock-Induced Transition to Turbulence,” Advances in Fluid Mechanics, V, A.Mendes, M.Rahman, and C.Brebbia, eds., WIT Press, Boston, MA, pp. 139–148.
Orlicz, G., 2012, “Incident Shock Mach Number Effects on Richtmyer–Meshkov Mixing With Simultaneous Density and Velocity Measurements,” Ph.D. dissertation, The University of New Mexico, Albuquerque, NM.


Grahic Jump Location
Fig. 1

Schematic of experimental setup. Top: side view. The pressurized driver section (left) is separated from the driven section by a diaphragm that is ruptured by a puncturer (dashed lines, labeled P). As the shock (schematically shown) propagates through the driven section, its progress is monitored by two pressure transducers (PT1 and PT2). Glycol droplets suspended in gas (air, air–SF6 mix, or air–SF6-acetone mix) are injected into the test section vertically, forming a seeded cylindrical column. The flow is illuminated with a horizontal laser sheet through an optical window in the end of the runoff section. The flow is illuminated with several pulsed lasers forming a single laser sheet, illuminating the horizontal centerline cross section of the shock tube. A 45° mirror above the test section reflects the view of the illuminated plane toward the camera. Bottom: close-up of the laser-illuminated flow visualization plane in the test section as seen by the camera, showing the cross section of the initial conditions and a dynamic image with a CRVP.

Grahic Jump Location
Fig. 2

Comparison of flow visualization with glycol tracers (top) and numerical modeling (bottom) of RMI developing in a heavy-gas column with a cylindrical cross section and an initially diffuse interface between the heavy gas and the surrounding light gas (after [19]). Mach number is 1.7, Atwood number is 0.5. Shock direction is from left to right. The leftmost image in each row represents the initial conditions immediately before shock arrival from the left. Time intervals between subsequent dynamic images are 50 μs. Experimental images are inverted, so darker areas correspond to seeded (injected) flow.

Grahic Jump Location
Fig. 6

Growth rates of the CRVP (top row) and the spike (bottom row) as the function of time after shock acceleration and of downstream distance. Schematic above the top row of plots illustrates how the streamwise feature size was determined from the image (left), using spanwise-averaged intensity (right).

Grahic Jump Location
Fig. 3

PLIF image acquired at M = 1.7 about 600 μs after shock acceleration, with observable features labeled. Vertical image extent is 2 cm. The image is inverted (darker parts correspond to more fluorescence). A reflection off the bottom wall of the test section is also labeled.

Grahic Jump Location
Fig. 4

Visualization of the same PLIF image using two different palettes. The tone curve on the left emphasizes low-brightness features, and displays intensity levels (based on original 16-bit image) 1200–1600. On the right, intensity levels 1000–9000 are similarly displayed, resulting in a much shallower tone curve.

Grahic Jump Location
Fig. 5

Experimental flow visualization using PLIF and Mie scattering off droplets for three Mach numbers. Flow direction is from left to right. Experimental images are inverted, so darker areas correspond to flow seeded with droplets (label “Mie”) or marked with acetone tracer (label “PLIF”). Extent of imaged area is 10.09 cm. Individual image timings (with time t = 0 corresponding to shock reaching the center of the gas column) are labeled.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In