Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dynamic x-ray diffraction and nanosecond quantification of kinetics of formation of β -zirconium under shock compression

We report the atomic- and nanosecond-scale quantification of kinetics of a shock-driven phase transition in Zr metal. We uniquely make use of a multiple shock-and-release loading pathway to shock Zr into the β phase and to create a quasisteady pressure and temperature state shortly after. Coupling shock loading with in situ time-resolved synchrotron x-ray diffraction, we probe the structural transformation of Zr in the steady state. Our results provide a quantified expression of kinetics of formation of β-Zr phase under shock loading: transition incubation time, completion time, and crystallization rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fast imaging of laser ablation of cerium: Dependence of plasma plume dynamics and structure on native oxide growth

We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Gas density structure of supersonic flows impinged on by thin blades for laser-plasma accelerator targets

Density transition injection is an effective technique for controllably loading electrons into a trapped phase for laser plasma accelerators. One common technique to achieve the required fluid structure is to impinge a thin blade on the plume of a supersonic nozzle. Density transitions induced in this way are often assumed to be bow shocks and therefore sharp, but simulations and fluorescence measurements presented in this work show that in many cases of interest, the density transition accessible to a laser propagating transverse to the shock is an intercepting shock, and therefore, shock thickness and density vary with pressure, laser height, and blade position. The fluid dynamics of a supersonic nozzle impinged on by a thin, flat object are explored through simulations and relevant features are verified via planar laser-induced fluorescence measurements. Overall, the implications of the results for tuning electron beam injectors in laser plasma accelerators are discussed.

42 ENGINEERING↗
DOE OSTI2023

Modeling shock-induced void collapse in single-crystal Ta systems at the mesoscales

Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗
DOE OSTI2022

An Efficient B-Spline Lagrangian/Eulerian Method for Compressible Flow, Shock Waves, and Fracturing Solids

This study presents a new method for modeling the interaction between compressible flow, shock waves, and deformable structures, emphasizing destructive dynamics. Extending advances in time-splitting compressible flow and the Material Point Methods (MPM), we develop a hybrid Eulerian and Lagrangian/Eulerian scheme for monolithic flow-structure interactions. We adopt the second-order WENO scheme to advance the continuity equation. To stably resolve deforming boundaries with sub-cell particles, we propose a blending treatment of reflective and passable boundary conditions inspired by the theory of porous media. The strongly coupled velocity-pressure system is discretized with a new mixed-order finite element formulation employing B-spline shape functions. Shock wave propagation, temperature/density-induced buoyancy effects, and topology changes in solids are unitedly captured.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Ultrafast X-ray Diffraction Study of a Shock-Compressed Iron Meteorite above 100 GPa

Natural kamacite samples (Fe92.5Ni7.5) from a fragment of the Gibeon meteorite were studied as a proxy material for terrestrial cores to examine phase transition kinetics under shock compression for a range of different pressures up to 140 GPa. In situ time-resolved X-ray diffraction (XRD) data were collected of a body-centered cubic (bcc) kamacite section that transforms to the high-pressure hexagonal close-packed (hcp) phase with sub-nanosecond temporal resolution. The coarse-grained crystal of kamacite rapidly transformed to highly oriented crystallites of the hcp phase at maximum compression. The hcp phase persisted for as long as 9.5 ns following shock release. Comparing the c/a ratio with previous static and dynamic work on Fe and Fe-rich Fe-Ni alloys, it was found that some shots exhibit a larger than ideal c/a ratio, up to nearly 1.65. This work represents the first time-resolved laser shock compression structural study of a natural iron meteorite, relevant for understanding the dynamic material properties of metallic planetary bodies during impact events and Earth’s core elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Recreating Giants Impacts in the Laboratory: Shock Compression of MgSiO 3 Bridgmanite to 14 Mbar

Understanding giant impacts requires accurate description of how extreme pressures and temperatures affect the physical properties of the constituent materials. Here, we report shock experiments on two polymorphs of MgSiO 3 : enstatite and bridgmanite (perovskite) crystals.We obtain pressure-density shock equation of state to 14 Mbar and more than 9 g/cm 3 , a 40% increase in density from previous data on MgSiO 3 . Density-functional-theory molecular dynamics (DFT-MD) simulations provide predictions for the shock Hugoniot curves for bridgmanite and enstatite and suggest that the Grüneisen parameter decreases with increasing density. The good agreement between the simulations and the experimental data, including for the shock temperature along the enstatite Hugoniot reveals that DFT-MD simulations reproduce well the behavior of dense fluid MgSiO 3 .We also reveal a high optical reflectance indicative of a metal-like electrical conductivity which supports the hypothesis that magma oceans may contribute to planetary magnetic field generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Shock wave formation in the thermosphere by an earthgrazing fireball: Empirical evidence for volatile-enhanced hydrodynamic shielding

Hydrodynamic shielding is a theoretically well-established but observationally elusive and experimentally difficult-to-replicate phenomenon with implications that extend far beyond meteor physics. Rare earthgrazing meteoroids with infrasound signatures that penetrate to the ground can be used to probe hydrodynamic shielding that leads to strong shock formation at high altitude. Here, we report the first coordinated optical and multi-station infrasound observations of a centimeter-scale earthgrazing fireball that generated sustained cylindrical line shock at thermospheric altitudes near 92 km. The event was recorded by numerous optical stations and three infrasound arrays, allowing trajectory reconstruction, ablation behavior, acoustic source localization, and shock characteristics. Optical observations indicate early mechanical erosion and ablation/evaporation at exceptionally low dynamic pressure, consistent with a cometary or a porous, volatile-bearing CM chondritic object. Independent infrasound detections localize shock generation to multiple points along a 164 km trajectory segment near perigee. Weak-shock modeling yields a consistent blast radius of ∼30 m, implying an acoustic-equivalent source size far exceeding the physical dimensions of the ∼45 g nucleus. We demonstrate that classical gas dynamics and ablation-driven hydrodynamic shielding alone cannot account for these observations under ambient thermospheric conditions. We show that volatile release provides the additional flow-field density enhancement required to amplify hydrodynamic shielding, reduce the effective local Knudsen number, and sustain a shock envelope capable of radiating detectable infrasound. Furthermore, these results demonstrate that small, volatile-rich meteoroids can transiently establish continuum-like flow in rarefied environments.

Astrometry↗
DOE OSTI2025

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗
DOE OSTI2025

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. In conclusion, we present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Microphysically modified magnetosonic modes in collisionless, high-β plasmas

With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-β collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma β is increased. For mode amplitudes |δB/B 0 |≳0.3, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies |δB/B 0 |≳2β -1 . The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-β plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Verification of a Reduced-Order Extension Spring Model

Adopting reduced order models (ROMs) of springs lowers the computational cost of stronglink simulations. However, ROMs introduce currently unquantified error to such analyses. This study addresses that lack of data by comparing a hexahedral mesh to a commonly used ROM beam mesh. Two types of analyses were performed, a quasi-static displacement-controlled pull and a haversine shock, examining basic spring properties as well as dynamics and stress/strain data. Both tests showed good similarities between the hexahedral and beam meshes, especially when comparing reaction force and stress trends and maximums. Equivalent plastic strain results were not quite as favorable, indicating that the beam model may be less likely to correctly predict spring failure. Despite reducing computation times by over 48 hours in all shock cases, appropriate use of the ROM should carefully balance this advantage with its reduction in accuracy, especially when examining spring failure and outputting variables such as equivalent plastic strain.

42 ENGINEERING↗
DOE OSTI2022

High-precision inversion of dynamic radiography using hydrodynamic features

While radiography is routinely used to probe complex, evolving density fields in research areas ranging from materials science to shock physics to inertial confinement fusion and other national security applications, complications resulting from noise, scatter, complex beam dynamics, etc. prevent current methods of reconstructing density from being accurate enough to identify the underlying physics with sufficient confidence. In this work, we show that using only features that are robustly identifiable in radiographs and combining them with the underlying hydrodynamic equations of motion using a machine learning approach of a conditional generative adversarial network (cGAN) provides a new and effective approach to determine density fields from a dynamic sequence of radiographs. In particular, we demonstrate the ability of this method to outperform a traditional, direct radiograph to density reconstruction in the presence of scatter, even when relatively small amounts of scatter are present. Our experiments on synthetic data show that the approach can produce high quality, robust reconstructions. We also show that the distance (in feature space) between a testing radiograph and the training set can serve as a diagnostic of the accuracy of the reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Behavior of soda-lime silicate glass under laser-driven shock compression up to 315 GPa

Shock experiments give a unique insight into the behavior of matter subjected to extremely high pressures and temperatures. Understanding the behavior of materials under such extreme conditions is key to modeling material failure and deformation dynamics under impact. While studies on pure silica are extensive, the shock behavior of other commercial silicates that contain additional oxides has not been systematically investigated. To better understand the role of composition in the dynamic behavior of silicates, we performed laser-driven dynamic compression experiments on soda-lime glass (SLG) up to 315 GPa. Using the accurate pulse shaping offered by the long pulse laser system at the Matter in Extreme Conditions end-station at the Linac Coherent Light Source, SLG was shock compressed along the Hugoniot to multiple pressure-temperature points. Velocity Interferometer System for Any Reflector was used to measure the velocity and determine the pressure inside the SLG. The U s –u p relationship obtained agrees well with the previous parallel plate impact studies. Within the error bars, no transformation to the crystalline phase was observed up to 70 GPa, which is in contrast to the behavior of pure silica under shock compression. Our studies show that the glass composition strongly influences the shock compression behavior of the silicate glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding the plasticity contributions during laser-shock loading and spall failure of Cu microstructures at the atomic scales

A hybrid atomic-scale and continuum modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline (nc) Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. Here, the simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.

36 MATERIALS SCIENCE↗