Search NASA⌕ Search

SEARCH · Search NASA

Results for “Meshkov”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
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 OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2021

Hydrodynamic instability at impact interfaces and planetary implications

Abstract Impact-induced mixing between bolide and target is fundamental to the geochemical evolution of a growing planet, yet aside from local mixing due to jetting – associated with large angles of incidence between impacting surfaces – mixing during planetary impacts is poorly understood. Here we describe a dynamic instability of the surface between impacting materials, showing that a region of mixing grows between two media having even minimal initial topography. This additional cause of impact-induced mixing is related to Richtmyer-Meshkov instability (RMI), and results from pressure perturbations amplified by shock-wave refraction through the corrugated interface between impactor and target. However, unlike RMI, this new impact-induced instability appears even if the bodies are made of the same material. Hydrocode simulations illustrate the growth of this mixing zone for planetary impacts, and predict results suitable for experimental validation in the laboratory. This form of impact mixing may be relevant to the formation of stony-iron and other meteorites.

58 GEOSCIENCES↗
DOE OSTI2015

Stabilization of high-compression, indirect-drive inertial confinement fusion implosions using a 4-shock adiabat-shaped drive

Hydrodynamic instabilities and poor fuel compression are major factors for capsule performance degradation in ignition experiments on the National Ignition Facility. Using a recently developed laser drive profile with a decaying first shock to tune the ablative Richtmyer-Meshkov (ARM) instability and subsequent in-flight Rayleigh-Taylor growth, we have demonstrated reduced growth compared to the standard ignition pulse whilst maintaining conditions for a low fuel adiabat needed for increased compression. Using in-flight x-ray radiography of pre-machined modulations, the first growth measurements using this new ARM-tuned drive have demonstrated instability growth reduction of ∼4× compared to the original design at a convergence ratio of ∼2. Corresponding simulations give a fuel adiabat of ∼1.6, similar to the original goal and consistent with ignition requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Ab initio determination of the instability growth rate of warm dense beryllium-deuterium interface

Accurate knowledge about the interfacial unstable growth is of great importance in inertial confinement fusion. During implosions, the deuterium-tritium capsule is driven by laser beams or X-rays to access the strongly coupled and partially degenerated warm dense matter regime. At this stage, the effects of dissipative processes, such as diffusion and viscosity, have significant impact on the instability growth rates. Here, we present ab initio molecular dynamics simulations to determine the equations of state and the transport coefficients. Several models are used to estimate the reduction in the growth rate dispersion curves of Rayleigh-Taylor and Richtmyer-Meshkov instabilities with considering the presence of these dissipative effects. We show that these instability growth rates are effectively reduced when considering diffusion. The findings provide significant insights into the microscopic mechanism of the instability growth at the ablator-fuel interface and will refine the models used in the laser-driven hydrodynamic instability experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On shock driven jetting of liquid from non-sinusoidal surfaces into a vacuum

Other work employed Richtmyer-Meshkov theory to describe the development of spikes and bubblesfrom shocked sinusoidal surfaces. Here, we discuss the effects of machining different two-dimensional shaped grooves in copper and examine the resulting flow of the material after being shocked into liquid on release. For these simulations, a high performance molecular dynamics code, SPaSM, was used with machined grooves of kh 0 = 1 and kh 0 = 1/8, where 2h 0 is the peak-to-valley height of the perturbation with wavelength λ, and k = 2π/λ. The surface morphologies studied include a Chevron, a Fly-Cut, a Square-Wave, and a Gaussian. Furthermore, we describe extensions to an existing ejecta source model that better captures the mass ejected from these surfaces. We also investigate the same profiles at length scales of order 1 cm for an idealized fluid equation of state using the FLASH continuum hydrodynamics code. Our findings indicate that the resulting mass can be scaled by the missing area of a sinusoidal curve with an effective wavelength, λeff , that has the same missing area. Finally, our extended ejecta mass formula works well for all the shapes considered and captures the corresponding time evolution and total mass.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Multimode instability evolution driven by strong, high-energy-density shocks in a rarefaction-reflected geometry

We present an experiment using lasers to produce a shock pressure of >10 Mbar, which we then use to drive Richtmyer–Meshkov and Rayleigh–Taylor growth at a 2D multimode perturbed interface. Key features of this platform are that we can precisely reproduce the perturbation from iteration to iteration of the experiment, facilitating analysis, and that the lasers allow us to produce very strong shocks, creating a plasma state in the system. We also implement a Bayesian technique to analyze the multimode spectra. This technique enables us to draw quantitative conclusions about the spectrum, even in the presence of significant noise. For instance, we measure the signal contained in the seeded modes over time, as well as the transition of the initial growth rate of these modes into the overall saturation behavior of the spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Shock-driven discrete vortex evolution on a high-Atwood number oblique interface

Here, we derive a model describing vorticity deposition on a high-Atwood number interface with a sinusoidal perturbation by an oblique shock propagating from a heavy into a light material. Limiting cases of the model result in vorticity distributions that lead to Richtmyer-Meshkov and Kelvin-Helmholtz instability growth. For certain combinations of perturbation amplitude, wavelength, and tilt of the shock, a regime is found in which discrete, co-aligned, vortices are deposited on the interface. The subsequent interface evolution is described by a discrete vortex model, which is found to agree well with both RAGE simulations and experiments at early times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Proton radiography measurements and models of ejecta structure in shocked Sn

We present the experimental validation of ejecta source, mass and velocity models using proton radiography. Ejecta mass and velocity measurements were performed under 7 atm of He gas on Sn targets with 4 concentric sinusoidal features with radial wavelength of λ=2 mm and amplitude η 0 = 0.159 mm. 42 images at 275 and 400 ns intervals were obtained between 0 and 14 µs from the time of shock breakout. Volume density profiles for the spikes were determined through Abel inversion and spike, bubble and surface velocities were determined using proton radiography and also verified through Photon Doppler Velocimetry (PDV). Ejecta sheets evolve to a self-similar density distribution that depends on a scaling variable z/η s t where η˙s is the spike tip velocity, z is the distance from the free surface and t is the time after shock breakout. Both the density profiles and the time dependence of the mass per unit area in the evolving spikes are in good agreement with a Richtmyer-Meshkov instability based model for ejecta production and evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Turbulent mixing and transition criteria of flows induced by hydrodynamic instabilities

In diverse areas of science and technology, including inertial confinement fusion (ICF), astrophysics, geophysics, and engineering processes, turbulent mixing induced by hydrodynamic instabilities is of scientific interest as well as practical significance. Because of the fundamental roles they often play in ICF and other applications, three classes of hydrodynamic instability-induced turbulent flows—those arising from the Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities—have attracted much attention. ICF implosions, supernova explosions, and other applications illustrate that these phases of instability growth do not occur in isolation, but instead are connected so that growth in one phase feeds through to initiate growth in a later phase. Essentially, a description of these flows must encompass both the temporal and spatial evolution of the flows from their inception. Hydrodynamic instability will usually start from potentially infinitesimal spatial perturbations, will eventually transition to a turbulent flow, and then will reach a final state of a true multiscale problem. Indeed, this change in the spatial scales can be vast, with hydrodynamic instability evolving from just a few microns to thousands of kilometers in geophysical or astrophysical problems. These instabilities will evolve through different stages before transitioning to turbulence, experiencing linear, weakly, and highly nonlinear states. Here, the challenges confronted by researchers are enormous. The inherent difficulties include characterizing the initial conditions of such flows and accurately predicting the transitional flows. Of course, fully developed turbulence, a focus of many studies because of its major impact on the mixing process, is a notoriously difficult problem in its own right. In this pedagogical review, we will survey challenges and progress, and also discuss outstanding issues and future directions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Using cylindrical implosions to investigate hydrodynamic instabilities in convergent geometry

Hydrodynamic instabilities such as the Rayleigh–Taylor (RT) and Richtmyer–Meshkov instabilities disrupt inertial confinement fusion (ICF) implosions through the growth of 3D perturbations. Growth of these 3D imperfections at the interfaces of an ICF capsule during implosion lead to mixing between materials that is detrimental to performance. These instabilities have been studied extensively in planar geometry, but such experiments lack the effects of convergence in spherical implosions. While several studies have been performed in spherical geometry, these often lack a direct means to measure perturbation growth. Experiments in cylindrical geometry include convergence effects while maintaining direct diagnostic access. Although cylinders have less compression than spheres, they do provide an excellent platform to validate modeling for convergent geometries. The problem with previous cylindrical implosion experiments was that the convergence ratios were limited to ~4. With the National Ignition Facility (NIF), larger cylindrical targets can be driven to convergences of 10–15 while maintaining a large enough final diameter to measure perturbation growth. This paper reviews the design process used to both benchmark radiation hydrodynamics codes and enable 1D post-processed simulations to explore design space to separate compression effects from acceleration/deceleration RT instability. Results from 1D simulations suggest that cylindrical implosions on the NIF can produce high-convergence experiments to validate RT instability growth for ICF implosions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Equation-of-state measured via x-ray phase contrast imaging for Epon 828/DEA epoxy

Epoxies are a broad class of polymer materials often used as adhesive, structural or binding materials. Epon 828 is an epoxy resin that can be polymerized with a variety of curing agents with the choice of curing agent potentially having an effect on the resulting epoxy polymer’s material properties. In this study, the dynamic behavior of Epon 828 epoxy resin cured with diethanolamine (DEA) is investigated through a series of tamped Richtmyer-Meshkov instability (RMI) experiments measured with x-ray phase-contrast imaging. The measured shock and particle velocities are combined with data in the literature to calibrate Mie-Grüneisen equations-of-state (EOS) for portions and combinations of the collective dataset. The calibrated Mie-Grüneisen EOS are validated against particle velocity profiles extracted from published literature using the Eulerian hydrocode CTH. Here, the Mie-Grüneisen EOS fit to only the tamped RMI experimental data presented here most closely follows the particle velocity profile in the published literature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental evaluation of shear modulus scaling of dynamic strength at extreme pressures

Recent progress in the development of dynamic strength experimental platforms is allowing for unprecedented insight into the assumptions used to construct constitutive models operating in extreme conditions. In this work, we make a quantitative assessment of how tantalum strength scales with its shear modulus to pressures of hundreds of gigapascals through a cross-platform examination of three dynamic strength experiments. Specifically, we make use of Split–Hopkinson pressure bar and Richtmyer–Meshkov instability experiments to assess the low-pressure strain and strain rate dependence. Concurrent examination of magnetically driven ramp-release experiments up to pressures of 350GPa allows us to examine the pressure dependence. In this work, using a modern description of the shear modulus, validated against both ab initio theory and experimental measurements, we then assess how the experimentally measured pressure dependence scales with shear modulus. Furthermore, we find that the common assumption of scaling strength linearly with the shear modulus is too soft at high pressures and offer discussion as to how descriptions of slip mediated plasticity could result in an alternative scaling that is consistent with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Further developments on observations of the Taylor instability in a dusty plasma

We previously reported the observation of a hydrodynamic-type interchange instability occurring sporadically at the curved interface between two dust clouds of different densities, in a three-dimensional dusty plasma [Pacha et al., Phys. Plasmas 19, 014501 (2012)]. The mechanism for the instability was not conclusively identified. Further analysis of the video images has been performed, showing that the instability is associated with the appearance of large amplitude dust acoustic waves or shocks. Evidence of a dust acoustic shock at the boundary, which could drive a Richtmyer–Meshkov-type instability, will be presented. We also discuss the possibility, suggested by Avinash and Sen [Phys. Plasmas 22, 083707 (2015)], that the observed instability is due to a Rayleigh–Taylor instability driven by a component of the gravitational acceleration, which is perpendicular to the curved boundary.

Pacha, K. A.↗
DOE OSTI2021

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗