Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock instability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Velocity-space signatures of shock-drift acceleration at quasi-perpendicular collisionless shocks

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3 ≤ M A ≤ 15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6 ≤ M f ≤ 9.4) and shock-normal angle 45°≤ θ Bn ≤ 90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers M A ≳8 (M f ≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks.

Howes, Gregory G. [Univ. of Iowa, Iowa City, IA (U↗
DOE OSTI2022

Experiments on the single-mode Richtmyer–Meshkov instability with reshock at high energy densities

Here, the hydrodynamic instability growth of a reshocked single-mode interface between high energy density fluids is studied. A laser-driven shock wave is used to drive an initially solid, sinusoidal interface between a dense plastic (1.43 g/cc) and a light foam (≈ 0.110 g/cc). After the interface has grown to a nonlinear state where the amplitude is of order of the wavelength, it is reshocked. The reshock compresses the nonlinear perturbation, which then grows at about twice the rate. While the pre-reshock growth rate is sensitive to the initial amplitude and wavelength of the perturbation, the post-reshock growth rate is comparatively insensitive to the initial condition. Qualitatively, we observe that the perturbations are less coherent after reshock, consistent with the idea that having a reshock accelerates the transition to turbulence. We find that some memory of the initial condition remains, even after reshock at late time: it appears if the initial perturbations have large enough wavelengths, and the flow structure of size comparable to the initial wavelength persists through reshock. Our results agree with design simulations and are consistent with the phenomenology of reshock studies in conventional gaseous shock tubes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Turbulence generation by shock interaction with a highly nonuniform medium

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, $\sigma_\rho/\bar{\rho}\sim1$) and strong shocks (shock Mach numbers $\bar{M}_s\gtrsim10$). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing the performance of MagLIF with 3D MHD simulations

The Magnetized Liner Inertial Fusion (MagLIF) experimental platform at Sandia National Laboratories has realized a number of recent improvements in electrical current and laser preheat coupling but fusion yields are below expectations, based on modeling with clean (that is, simulations with no interfacial instabilities and no mix) two-dimensional (2D) magneto-hydrodynamics (MHD) simulations, by a factor of 5 or more. However, magnetized liner implosions are known to produce helical magneto-Rayleigh–Taylor (MRT) modes and complex stagnation structures that cannot be directly modeled in 2D. This paper presents the results of 3D HYDRA MHD simulations, including helical MRT, showing that degradation from these instabilities can readily reproduce experimental yields, but it is difficult to simultaneously match the yield along with inferred fuel temperature, pressure, and burn history (assessed with x-rays), which are larger in the experiments considered here. Additional analysis methods and future experiments are proposed to help address the discrepancies. The 3D simulations also show, without substantive improvements to stability, changing the applied Bz or preheat is unlikely to increase performance of MagLIF to the same degree as clean 2D simulations. Finally, the first 3D HYDRA simulations with the Hall term show the ab initio production of helical MRT and produce a clear change in stagnation morphology compared to pre-seeded simulations without Hall. In spite of the differences, the simulations still produce comparable fusion performance but with longer, lower power, neutron yield history.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Flow visualization in high-energy-density laboratory plasmas

Flow visualization is often essential for evaluating high-energy-density (HED) laboratory plasmas, examining their dynamics, and assessing design strategies. However, imaging at fine phenomenological flow scales is a major challenge. Here, a zone-plate–based radiography system tailored for applications in HED physics is demonstrated. The burst-mode radiography system combines a dual-flash laser-generated x-ray source, a high-magnification zone plate optic, and a multiple frame hybrid complementary metal-oxide-semiconductor (hCMOS) sensor for image detection. Operating at a magnification of 42 x⁠, a pair of radiographs is measured in rapid succession, each time gated by a 4.75 keV flash of x rays lasting 100 ps. The system achieves a spatial resolution below 3 μm. The potential of this approach for applications in HED physics is exemplified by observing small-scale spike morphologies at a pusher-foam interface impacted by a strong, laser-driven blast wave.

Flow instabilities↗
DOE OSTI2025

Optimizing the Smoothness and Thickness Uniformity of Thin-Film Parylene-N Vapor-Deposited Coatings for Inertial Confinement Fusion Experiments

Polymer coatings with submicrometer smoothness and constant thickness are a required component in a variety of inertial confinement fusion experiments. Smoothness is important for minimizing Rayleigh-Taylor-driven hydrodynamic instabilities, and uniform thickness is important for uniform shock propagation and shell convergence, both of which are critical phenomena that affect the experiment. The preferred polymer coating method is to vapor deposit the parylene-N polymer because it provides nominally smooth conformal coatings. As the coating thickness exceeds 5 µm, however, dome-shaped nodular growth defects develop and the thickness will vary by up to 17% over a distance of 3 cm. This study presents a deterministic method for achieving uniform film thicknesses with ±2% variability over 3 cm and a predictive method to control the thickness to within 5% of the desired value. A coating smoothness of ∼50 nm rms, measured over 40 000 µm 2 , was achieved by adding additional surfaces near the substrates. This additional area improved the thickness uniformity, an effect that is attributed to the low sticking coefficient of the parylene monomer.

chemical vapor deposition (CVD)↗
DOE OSTI2023

X-ray Pump-Probe Measurements using a Laser-Plasma Accelerator (Final Report)

This project funded one graduate student, Mario Balcazar, to perform experiments on high power laser facilities to study the application of X-rays generated by laser wakefield accelerated electrons for advanced radiography, in particular X-ray phase contrast imaging of hydrodynamic shocks/instabilities. The work made use of the DOE’s LaserNetUS facilities. We performed dynamic phase-contrast X-ray imaging of the shock generated by a 200 ps duration laser pulse (E = 1 J, I = 10 15 Wcm –2 ) with a 30 µm diameter liquid target, with unprecedented spatio-temporal resolution. This includes multi-shock generation within the liquid jet and plasma instability formation. Innovative electron-beam radiography was used to probe the laser-plasma interplay finding evidence of bilateral heating of the water followed by strong electric field generation. These measurements help explain some of the discrepancies between simulation and experiment and pave the way to better plasma diagnostic systems in HED and ICF physics experiments. This work was performed in collaboration with researchers from Michigan, LLNL, LBL and Sandia National Laboratories, Queen’s University Belfast and Imperial College.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2021

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗
DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks

The major goals of the project are to study three fundamental questions related to collisionless, non-relativistic, quasi-perpendicular shocks by combining continuum Vlasov-Maxwell simulations (the Gkeyll framework) and high resolution, in situ heliospheric observations from the MMS spacecraft. The specific questions are: 1. How and where are electrons energized and heated in a quasi-perpendicular shock, including the ultimate source of entropy production? 2. What role do waves, turbulence, and secondary instabilities play in the conversion of directed flow energy to electron energy? 3. How do the parameters of the shock determine the energization of electrons?

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
DOE OSTI2023

Soliton formation and topology manipulation of coupled spins via ultrafast re-magnetization

The major goal of the project was to explore the properties of magnetic order driven far out of equilibrium by optical excitations. These include fundamental questions related to the interplay of magnetic, structural and electronic degrees of freedom in materials that are optically excited. These problems break down in short term quenching of the magnetization as a results of the energy of the optical pulse studying how the energy of the optical pulse and the angular momentum of the magnetization flow between the different degrees of freedom. This is followed by the longer time re-emergence of the magnetism as the system cools. Within this goal we explored the formation of solitons, both with and without non-trivial topology, via rapid re-magnetization processes after optical-driven ultrafast demagnetization. Guided by theory, we predict that turbulence and modulational instabilities will drive the formation of solitons, including dispersive shock waves, magnon droplets, and skyrmions. In the case of topological defect formation, i.e. skyrmion generation, related processes have long been predicted based on general principles of symmetry-breaking phase transitions; the density of topological defects in a long-range-ordered phase can be controlled by varying the quench rate through the second-order phase transition, i.e. the Kibble-Zurek (KZ) mechanism. We, more broadly, had the goal to study a range of other emergent magnetic behaviors after optical excitation. While great attention has been given to ultrafast demagnetization, less is known of the subsequent spin dynamics and coupling to the lattice. What is increasingly appreciated is that ultrafast demagnetization leads to spin currents that can carry angular momentum from the rapidly demagnetized sample. These spin currents appear in many ways but fundamentally controls demagnetization, drives interactions between regions of the material, and control the spin resulting structure. They can be probed many ways including THz emission which, in turn, gives insight into the demagnetization processes. Simultaneously if there is magneto-elastic coupling ultrafast demagnetization can drive structural excitations that are not expected from thermal energy added to the lattice. We explored a range of complex phenomena that arise from photo-excitation of magnetic systems with the coupling of electronic, magnetic and structural degrees of freedom. We have made a number of fundamental discoveries that are reflected in our publications list with additional work still being prepared for publication. Highlights of this work include (i) ultra-efficient, nonlinear THz surface acoustics (ii) spin-current-mediated rapid magnon localization and coalescence, (iii) spin-wave soliton formation in ferromagnetic FePt nanoparticles, (iv) dynamic phonon coupling in elemental antiferromagnetic Cr, (v) THz emission from Co/Pt bilayers and FeRh/Pt bilayers, (vi) theoretical investigation of spin hydrodynamics, solitons and shock waves, and (vii) ultrafast perturbation of magnetic domains by optical pumping.

36 MATERIALS SCIENCE↗
DOE OSTI2026

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Understanding the evolution of liquid and solid microjets from grooved Sn and Cu samples using radiography

Experiments were performed on grooved Sn and Cu samples to study the temporal evolution of microjets. Jets were generated by the impact of gun-launched flyer plates against the back of grooved targets made from either Cu or Sn (groove depth of ~250 μm). The Hugoniot states in the various Sn targets encompassed conditions where solid phases are maintained throughout (7 and 16 GPa) and also conditions where melting occurs upon the release of compression (25 and 34 GPa); the transition occurs near a Hugoniot pressure of 23 GPa. Cu targets at 27 and 56 GPa provide comparisons in which the jets move at similar speeds but remain solid. In all cases, the spatial distribution of mass within the microjets was measured using high-speed synchrotron radiography. The result is a time history of the jet thickness profile from which quantities like total jet mass and jet velocity can be derived. In both the solid and liquid states, we generally observe that an increase in the shock strength leads to an increase in jet mass. However, this trend breaks down for Hugoniot states near the transition from continuously solid to melted-on-release. This is evidenced by the observation that there was no difference in the rate of mass flow in Sn jets at 16 and 25 GPa, while similar pressure jumps on either side of this range caused substantial changes in the jet mass. This contrasts with the behavior of smaller polishing defects that were present on the same samples (~1 μm deep). From these, no ejecta mass was detected below the melt boundary, but obvious microjets were generated once melting occurred. This indicates that crossing the bulk melt-on-release threshold can alternately promote or inhibit the flow of mass into microjets based on the amplitude of the initial perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗