Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock wave”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2020

The effects of mixture preburning on detonation wave propagation

Pressure gain combustion in the form of continuous detonations can provide a significant increase in the efficiency of a variety of propulsion and energy conversion devices. In this regard, rotating detonation engines (RDEs) that utilize an azimuthally-moving detonation wave in annular systems are increasingly seen as a viable approach to realizing pressure gain combustion. However, practical RDEs that employ non- premixed fuel and oxidizer injection need to minimize losses through a number of mechanisms, including turbulence-induced shock-front variations, incomplete fuel-air mixing, and premature deflagration. In this study, a canonical stratified detonation configuration is used to understand the impact of preburning on detonation efficiency. It was found that heat release ahead of the detonation wave leads to weaker shock fronts, delayed combustion of partially-oxidized fuel-air mixture, and non-compact heat release. Furthermore, large variations in wave speeds were observed, which is consistent with wave behavior in full-scale RDEs. Peak pressures in the compression region or near triple points were considerably lower than the theoretically-predicted values for ideal detonations. Analysis of the detonation structure indicates that this deflagration process is parasitic in nature, reducing the detonation efficiency but also leading to heat release far behind the wave that cannot directly strengthen the shock wave. As a result, this parasitic combustion leads to commensal combustion (heat release far downstream of the wave), indicating that it is the root cause of combustion efficiency losses.

42 ENGINEERING↗
DOE OSTI2022

Development of slurry targets for high repetition-rate x-ray free electron laser experiments

Combining an x-ray free electron laser with a high-power laser driver enables the study of equations-of-state, high strain-rate deformation processes, structural phase transitions, and transformation pathways as a function of pressure to hundreds of GPa along different thermodynamic compression paths. Future high repetition-rate laser operation will enable data to be accumulated at >1 Hz, which poses a number of experimental challenges, including the need to rapidly replenish the target. Here, we present a combined shock compression and an x-ray diffraction study on epoxy (50% vol.)-crystalline grains (50% vol.) slurry targets, which can be fashioned into extruded ribbons for high repetition-rate operation. For shock-loaded NaCl-slurry samples, we observe pressure, density, and temperature states within the embedded NaCl grains consistent with observations from shock-compressed single-crystal NaCl.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Sound speed measurements in shock compressed cemented tungsten carbide: Evolution of elastic moduli with damage at pressures to 100 GPa

The motivation of the present study is to gain insights into the evolution of elastic properties of cemented tungsten carbides (WC) shock compressed to 100 GPa. Seven plate impact experiments—two front surface impact and five release wave overtake—are conducted to make simultaneous measurements of Hugoniot states and longitudinal sound speeds in shocked WC with 3.7wt.% cobalt binder. The sound speeds along with estimates for bulk sound speeds, obtained using the Birch–Murnaghan EoS, are analyzed to determine the elastic moduli—longitudinal, bulk, and shear—as a function of Hugoniot stress. The longitudinal and bulk sound speeds at Hugoniot states of interest are found to increase linearly with longitudinal stress. Consistent with the increase in sound speeds, the longitudinal and bulk moduli also increase with Hugoniot stress; however, the increase in longitudinal modulus is modest when compared to predictions of theoretical models that account for pressure and temperature dependence of elastic moduli, but with no damage. The shear moduli remain nearly constant at ~318 GPa over the range of Hugoniot states investigated. These values are, however, much lower than those predicted by the Steinberg–Guinan model with no damage. Poisson’s ratio decreases initially from its ambient value of 0.208 to ~0.199 for Hugoniot stress ≤10 GPa indicating consolidation of the WC microstructure with low initial stress; however, with an increase in Hugoniot stress to ~100 GPa, Poisson’s ratio increases to ~0.317, indicating degradation of shear moduli with increasing stress. The product of density and Grüneisen parameter (ρΓ), after an initial spike, remains nearly constant for volumetric strains ≥0.07. The maximum average temperature rise is estimated to be ~286°C at the highest Hugoniot stress employed in the study.

36 MATERIALS SCIENCE↗
DOE OSTI2026

XFEL imaging techniques for high energy density and inertial fusion energy research at HED-HiBEF

The imaging platform developed at the High Energy Density-Helmholtz International Beamline for Extreme Fields (HED-HiBEF) instrument at the European X-ray Free Electron Laser (XFEL) and its applications to HED and fusion related research are presented. The platform combines the XFEL beam with the high-intensity short-pulse laser ReLaX and the high-energy nanosecond-pulse laser DiPOLE-100X. The spatial resolution is better than 500 nm and the temporal resolution of the order of 50 fs. The influence of the XFEL source in the x-ray imaging method is discussed. Free-propagation x-ray phase contrast imaging and Talbot-Lau imaging setups are shown. We show examples of blast waves and converging cylindrical shocks in aluminum, resonant absorption measurements of specific charged states in copper with ReLaX and planar shocks in polystyrene material generated by DiPOLE-100X. For the first time, we show the application of Talbot-Lau interferometry to convergent cylindrical shocks as well as resonant absorption processes. We also discuss the possibilities introduced by combining this imaging platform with a kJ-class laser.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Few-femtosecond resolved imaging of laser-driven nanoplasma expansion

Abstract The free expansion of a planar plasma surface is a fundamental non-equilibrium process relevant for various fields but as-yet experimentally still difficult to capture. The significance of the associated spatiotemporal plasma motion ranges from astrophysics and controlled fusion to laser machining, surface high-harmonic generation, plasma mirrors, and laser-driven particle acceleration. Here, we show that x-ray coherent diffractive imaging can surpass existing approaches and enables the quantitative real-time analysis of the sudden free expansion of laser-heated nanoplasmas. For laser-ionized SiO 2 nanospheres, we resolve the formation of the emerging nearly self-similar plasma profile evolution and expose the so far inaccessible shell-wise expansion dynamics including the associated startup delay and rarefaction front velocity. Our results establish time-resolved diffractive imaging as an accurate quantitative diagnostic platform for tracing and characterizing plasma expansion and indicate the possibility to resolve various laser-driven processes including shock formation and wave-breaking phenomena with unprecedented resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗
DOE OSTI2023

Kinetic characterization of 3D magnetic reconnection: A transformative step (Final Report)

We investigated magnetic reconnection in the Earth’s magnetopause, magnetotail, and the bow shock, by means of fully kinetic simulations and analysis of space data obtained by NASA’s Magnetospheric Multiscale mission. We elucidated signatures of particle energization in particle distribution functions in the vicinity of reconnection X-line, properties of reconnection in shock-driven turbulence, and wave excitations and nonlinear structures. We advanced the understanding of magnetic reconnection from the viewpoint of particle kinetics in the diffusion region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Go with the Radflow: Thermal radiation experiments in the high energy-density regime [Slides]

Los Alamos National Laboratory is a nuclear weapons laboratory supporting our nation's defense. In support of the mission is a high energy-density physics program in which we design and execute experiments to study radiation-hydrodynamics phenomena and improve the predictive capability of our large-scale multi-physics software codes on our supercomputers. The Radflow Project is maturing a unique spectroscopic measurement that yields a spatially dependent supersonic radiation wavefront profile. The spectroscopic measurement is much more constraining than the typically used radiography, which measures the density changes due to the shock as the radiation wave goes subsonic. Recently, the spectroscopic measurement has been applied to advanced targets with a detailed heterogeneity representing a single realization of a stochastic medium. The Radflow Project also conducts opacity experiments to explain the current conundrum of the experimental iron opacity differing from theory, which is important in modeling the sun. More than that, though, we need to be able to answer these science questions and be able to validate our theory and codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Shock-induced melting of [100] lithium fluoride: Sound speed and Hugoniot measurements to 230 GPa

Although [100] lithium flouride (LiF) is the most widely used optical window material in dynamic compression experiments, its high stress (>100 GPa) shock compression response, including melting, is not well understood. To address this need, we measured wave profiles in plate impact experiments to determine the Hugoniot states and longitudinal sound speeds in [100] LiF crystals shock compressed to 231 GPa. The measured peak states are fitted well by a linear shock velocity-particle velocity relation, providing an accurate determination of the LiF Hugoniot curve to significantly higher stresses than previous experiments. The longitudinal sound speeds show a near linear increase with density compression to 182 GPa. Between 182 GPa and 195 GPa, the sound speed and the longitudinal modulus decrease abruptly, due to shock-induced melting. The increasing sound speeds and moduli at higher stresses suggest that shock compressed LiF is fully liquid at 195 GPa and above, allowing determination of the Gruneisen parameter for liquid LiF. Here, the melt stress determined here differs from that predicted by current multiphase equations of state for LiF. Our results provide important insight into the high stress solid and liquid states of shock compressed LiF and point to the need for an improved multiphase equation of state at high pressures and high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Global Hybrid Simulations of Interaction Between Interplanetary Rotational Discontinuity and Bow Shock/Magnetosphere: Can Ion-Scale Magnetic Reconnection be Driven by Rotational Discontinuity Downstream of Quasi-Parallel Shock?

Ion-scale magnetic reconnection has been observed downstream of the terrestrial quasi-parallel (Q-∥) shock. Whether it is driven by interplanetary discontinuities or turbulent Q-∥ shock, however, is unclear. Using three-dimensional global hybrid simulation, we investigate the generation of magnetic reconnection downstream of the Q-∥ shock, while an interplanetary rotational discontinuity (RD) is launched to the bow shock. Cases with various solar wind Alfvén Mach numbers, M A = 3.0 to 8, and propagation directions of the RD are presented. The propagation direction n is assumed to be in the GSE xz plane and pointing earthward, with n = (-sin(θ 12 /2),0,-cos(θ 12 /2)), where θ 12 is the angle between the upstream (B1) and downstream (B2) magnetic fields across the transmitted RD. It is found that magnetic reconnection occurs inside the RD downstream of the Q-∥ shock, forming flux ropes extending along the dawn-dusk direction about tens of ion inertial lengths. Large-amplitude low-frequency waves originated from the Q-∥ shock lead to the bending and squeezing of the field lines around the RD, which play an important role in triggering reconnection inside the RD. As the RD impacts the dayside magnetopause, magnetopause reconnection takes place between the field lines behind the RD and geomagnetic field lines. Nevertheless, no reconnection is found downstream of the Q-∥ shock itself or outside the RD in the magnetosheath. The existent and structure of reconnection in the magnetosheath are found to strongly depend on the parameters M A and n. Our simulation shows that ion-scale magnetic reconnection is driven by an external driver in the form of the compression of an RD around the bow shock and in the magnetosheath, rather than caused by the turbulent Q-∥ shock alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Proto-neutron star convection and the neutrino-driven wind: implications for the r-process

The neutrino-driven wind from proto-neutron stars is a proposed site for r-process nucleosynthesis, although most previous work has found that a wind heated only by neutrinos cannot produce the third r-process peak. However, several groups have noted that introducing a secondary heating source within the wind can change the hydrodynamic conditions sufficiently for a strong r-process to proceed. One possible secondary heating source is gravito-acoustic waves, generated by convection inside the proto-neutron star. As these waves propagate into the wind, they can both accelerate the wind and shock and deposit energy into the wind. Additionally, the acceleration of the wind by these waves can reduce the total number of neutrino captures and thereby reduce the final electron fraction of the wind. In neutron rich conditions, all of these effects can make conditions more favourable for r-process nucleosynthesis. Here, we present a systematic investigation of the impact of these convection-generated gravito-acoustic waves within the wind on potential nucleosynthesis. We find that wave effects in the wind can generate conditions favourable for a strong r-process, even when the energy flux in the waves is a factor of 10 –4 smaller than the total neutrino energy flux and the wind is marginally neutron rich. Nevertheless, this depends strongly on the radius at which the waves become non-linear and form shocks. We also find that both entropy production after shock formation and the acceleration of the wind due to stresses produced by the waves prior to shock formation impact the structure and nucleosynthesis of these winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗
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 OSTI2023

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Peculiarities of planar shockwave interaction with air–water interface and solid target

We present the results of the research on underwater electrical explosion of planar copper wire arrays, accompanied by the generation of a planar shock, using the synchrotron-based phase-contrast radiography imaging capabilities of the ID19 beamline at the European Synchrotron Radiation Facility. It is shown that the interaction of a strong shock with an air–water interface leads to additional acceleration of the wire expansion due to a rarefaction wave along with the emergence of micro-jets. In the case of a target placed above the array, we observed the formation of a cavity between the array and the target due to the interaction of two rarefaction waves. Furthermore, the results of two-dimensional hydrodynamic simulations of the wire explosion and the interaction of the generated shock with the water–air and water-target interfaces showed good agreement with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Subnanosecond phase transition dynamics in laser-shocked iron

Iron is one of the most studied chemical elements due to its sociotechnological and planetary importance; hence, understanding its structural transition dynamics is of vital interest. By combining a short pulse optical laser and an ultrashort free electron laser pulse, we have observed the subnanosecond structural dynamics of iron from high-quality x-ray diffraction data measured at 50-ps intervals up to 2500 ps. We unequivocally identify a three-wave structure during the initial compression and a two-wave structure during the decaying shock, involving all of the known structural types of iron (α-, γ-, and ε-phase). In the final stage, negative lattice pressures are generated by the propagation of rarefaction waves, leading to the formation of expanded phases and the recovery of γ-phase. Our observations demonstrate the unique capability of measuring the atomistic evolution during the entire lattice compression and release processes at unprecedented time and strain rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗