Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗
DOE OSTI2021

Self-consistent, high-order spatial profiles in a model for two-fluid turbulent mixing

A Reynolds-averaged Navier-Stokes model is presented with the property that it admits self-consistent, high-order spatial profiles in simulations of two-fluid turbulent mixing layers. Whereas previous models have been limited by the assumption of a linear mixing profile, the present paper relaxes this assumption and, as a result, is shown to achieve much better agreement with experimental profiles. Similarity analysis is presented to derive constraints on model coefficients to enforce desired self-similar growth rates that are fully consistent with the high-order spatial profiles. Through this similarity analysis, it is shown that care must be taken in model construction, as it is possible to construct certain terms in such a way as to leave growth rates unconstrained. This model, termed the k-Φ-L-a-V model, is then applied in simulations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz mixing layers. Overall, these simulations confirm that the expected growth parameters are recovered and high-order spatial profiles are maintained.

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

Reduced mixing in inertial confinement fusion with early-time interface acceleration

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. In conclusion, this characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Magnetohydrodynamic effects in a shock-accelerated gas cylinder

This work presents two-dimensional (2D) simulations on a cylindrical Richtmyer-Meshkov instability (RMI) in magnetohydrodynamics (MHD). Three studies are presented in an effort to quantify and qualify the evolution of the MHD RMI by varying the magnetic field orientation, strength of the magnetic field, and strength of the shock wave driving the instability. The orientations considered herein are either parallel or perpendicular to the shock wave motion. The second study varies the magnetic fields between 100, 250, and 500 gauss (G), while the third study considers incident shock wave Mach numbers M = 1.2, M = 1.66, and M = 2.2. These parameter ranges were selected to be easily achievable in experiments while the interface perturbation was selected such that its evolution is independent of either the shock wave or magnetic field orientations independently. It was found that the MHD RMI evolution is dependent upon the magnetic field orientation relative to the shock transit direction as well as their individual magnitudes. This is because the mechanism of suppression, attributed to Alfven waves, is a function of the magnetic field strength and the orientation of the magnetic field, while the mechanism of RMI evolution, baroclinic vorticity deposition, is a function of the Mach number. Stronger magnetic fields were found to provide greater mixing suppression and have significant effects on RMI-like interface morphology. Lastly, increasing the shock wave strength generated competing effects between higher RMI vorticity deposition and greater vorticity removal from the interface by faster Alfven waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Rayleigh–Taylor Instability With Varying Periods of Zero Acceleration

We present our findings from a numerical investigation of the acceleration-driven Rayleigh–Taylor Instability, modulated by varying periods without an applied acceleration field. It is well known from studies on shock-driven Richtmyer–Meshkov instability that mixing without external forcing grows with a scaling exponent as ≈ t 0.20-0.28 When the Rayleigh–Taylor Instability is subjected to varying periods of “zero” acceleration, the structural changes to the mixing layer remain remarkably small. After the acceleration is re-applied, the mixing layer quickly resumes the profile of development it would have had if there had been no intermission. As a result, this behavior contrasts in particular with the strong sensitivity that is found to other variable acceleration profiles examined previously in the literature.

42 ENGINEERING↗
DOE OSTI2023

Transition in helium bubble strengthening of copper from quasi-static to dynamic deformation

Damage from low-temperature irradiation and the subsequent degradation of materials performance pose significant challenges for the storage of radioactive materials and for peripheral components in some nuclear reactor designs. Fully understanding the mechanical behavior of such materials requires test data for strain rates in both the quasi-static (< 10/s) and dynamic (>> 10/s) regimes. While dynamic testing has generally been avoided in the past for neutron irradiated (contamination concerns) and ion irradiated (insufficient volume) materials, surface-sensitive Richtmyer-Meshkov instability (RMI) tests were used in the present work to overcome these limitations. Here, nanopillar compression, nanoindentation, and RMI testing data from a helium implanted surface layer (~10 µm thick) were compiled to explore the effects of helium bubbles on the materials strength of high-purity copper at strain rates of 0.001/s – 10 8 /s. While nano-mechanical testing revealed increases in yield strength and hardness with increasing helium dose from 1000 to 4000 appm He, RMI indicated no significant changes in strength as compared to unimplanted copper. Here, this discrepancy in behavior was rationalized through a combination of recent literature and follow-on molecular dynamics (MD) simulations, leading to the conclusion that the nanoscale helium bubbles acting as dispersed barriers to dislocation motion at quasi-static strain rates collapse under shock loading and cease to be effective barriers at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
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 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↗