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

Detonation-induced transformation of graphite to hexagonal diamond

In this work, we explore the structural evolution of highly oriented pyrolytic graphite (HOPG) under detonation-induced shock conditions using in situ synchrotron x-ray diffraction in the ns timescale. We observe the formation of hexagonal diamond (lonsdaleite) at pressures above 50 GPa, in qualitative agreement with recent gas gun experiments. First-principles density functional calculations reveal that under uniaxial compression, the energy barrier for the transition toward hexagonal diamond is lower than that for cubic diamond. Finally, no indication of cubic diamond formation was observed up to >70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

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

High-Fidelity Analysis of ULF Wave Mode Structure Following Interplanetary Shock Compression of the Dayside Magnetopause Using MMS Multi-Point Observations

During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross-spectral analysis based on high-fidelity multi-probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode (|m| ≤ 3) magnetic fluctuations in the Pc4-5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump-on-tail distribution at m≈20 to a dominant peak at m≈50 in about 10 min. In addition, an interesting change of sign in m from negative to positive is observed as MMS crosses ~11 MLT pre-noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti-sunward propagation. Finally, for both shocks, the contribution of higher frequency waves (Pc-4 range compared with Pc-5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗
DOE OSTI2021

Time-resolved particle image velocimetry measurements of the turbulent Richtmyer–Meshkov instability

Experiments are presented on the Richtmyer–Meshkov instability (RMI) with a three-dimensional, multi-mode initial perturbation. The experiments use a vertical shock tube, where a stably stratified interface is formed between air and sulphur hexafluoride (SF 6 ) via counterflow. A perturbation is imposed at the interface by vertical oscillation of the gas column, forming Faraday waves. The interface is accelerated by a Mach 1.17 (in air) shock wave, and the development of the mixing region between the gases is investigated using particle image velocimetry. Following shock acceleration, a reflected shock wave from the bottom of the shock tube interacts with the mixing layer a second time (reshock). The experiment is initialized with both high and low amplitude perturbations to examine the effect of the perturbation amplitude on measured quantities. The instability growth exponent (θ) is determined from the kinetic energy field using the width of the mixing layer and the decay of kinetic energy, which are found to be in agreement when the flow is most strongly excited. A growth exponent of θ ≈ 0.5 is found for all cases except the high-amplitude reshocked regime (where θ ≈ 0.33). High-amplitude experiments exhibit the transitional outer Reynolds number (Re ≡ hh˙/ν >10 4 ) required for mixing transition following the incident shock, and both experiments are elevated well above this threshold following reshock. However, neither set of experiments meet the more stringent requirements proposed by Zhou et al. (Phys. Rev.E, vol. 67, issue 5, 2003) which include the time dependent aspect of the RMI, an observation which is also made when examining the spectra.

42 ENGINEERING↗
DOE OSTI2022

Shock compression behavior of stainless steel 316L octet-truss lattice structures

Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.

42 ENGINEERING↗
DOE OSTI2020

Thermal interfaces in dynamic compression experiments

Temperature is a difficult thermodynamic variable to measure in dynamic compression experiments. Optical pyrometry is a general-purpose technique for measuring temperature from a radiant surface, but that surface is often the interface between distinct materials with temperatures that vary spatially along the loading direction. This leads to a fundamental problem: how does the measured interface temperature relate to this temperature profile along the compression axis? Numerical analysis of loading history and thermal diffusion at these interfaces shows that seemingly subtle changes in experiment geometry can lead to very different temperature profiles. We compare these results to laboratory temperature measurements of shock-compressed tin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Plant heat stress: Concepts directing future research

Predicted increases in future global temperatures require us to better understand the dimensions of heat stress experienced by plants. In this review we highlight four key areas for improving our approach towards understanding plant heat stress responses. First, although the term "heat stress" is broadly used, that term encompasses heat shock, heat wave, and warming experiments, which vary in the duration and magnitude of temperature increase imposed. A greater integration of results and tools across these approaches is needed to better understand how heat stress associated with global warming will affect plants. Secondly, there is a growing need to associate plant responses to tissue temperatures. We review how plant energy budgets determine tissue temperature and discuss the implications of using leaf versus air temperature for heat stress studies. Third, we need to better understand how heat stress affects reproduction, particularly understudied stages such as floral meristem initiation and development. Fourth, we emphasize the need to integrate heat-stress recovery into breeding programs to complement recent progress in improving plant heat-stress tolerance. Taken together, we provide insights into key research gaps in plant heat stress and provide suggestions on addressing these gaps to enhance heat stress resilience in plants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Observation of breakdown wave mechanism in avalanche ionization produced atmospheric plasma generated by a picosecond CO 2 laser

Understanding the formation and long-timescale evolution of atmospheric plasmas produced by ultrashort, long-wavelength infrared (LWIR) pulses is an important but partially understood problem. Of particular interest are plasmas produced in air with a peak laser intensity ~10 12 W/cm 2 , the so-called clamping intensity observed in LWIR atmospheric guiding experiments where tunneling and multi-photon ionization operative at near-IR or shorter wavelengths are inoperative. We find that avalanche breakdown on the surface of aerosol (dust) particles can act to seed the breakdown of air observed above the 200 GW/cm 2 threshold when a train of 3 ps 10.6 μm laser pulses separated by 18 ps is used. The breakdown first appears at the best focus but propagates backward toward the focusing optic as the plasma density approaches critical density and makes forward propagation impossible. The velocity of the backward propagating breakdown can be as high as 10 9 cm/s, an order of magnitude greater than measured with ns pulse-produced breakdown, and can be explained rather well by the so-called breakdown wave mechanism. Transverse plasma expansion with a similar velocity is assisted by UV photoionization and is observed as a secondary longitudinal breakdown mechanism in roughly 10% of the shots. When a cm-size, TW power beam is propagated, interception of aerosol particles is guaranteed and several (40 cm -3 ) breakdown sites appear, each initially producing a near-critical density plasma. On a 10 ns–1 μs timescale, shockwaves from each site expand radially and coalesce to produce a large hot gas channel. The radial velocity of the expansion agrees well with the prediction of the blast wave theory developed for ultrafast atmospheric detonations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Characterizing a supernova’s standing accretion shock instability with neutrinos and gravitational waves

Here, we perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection P D as well as the false identification probability P FI . We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in P D (at P FI = 0.1) of SASI activities that is larger by up to ≈ 40% for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Magnetic field effects on laser energy deposition and filamentation in magneto-inertial fusion relevant plasmas

We report on experimental measurements of how an externally imposed magnetic field affects plasma heating by kJ-class, nanosecond laser pulses. The experiments reported here took place in gas cells analogous to magnetized liner inertial fusion targets. We observed significant changes in laser propagation and energy deposition scale lengths when a 12T external magnetic field was imposed in the gas cell. We find evidence that the axial magnetic field reduces radial electron thermal transport, narrows the width of the heated plasma, and increases the axial plasma length. Reduced thermal conductivity increases radial thermal gradients. This enhances radial hydrodynamic expansion and subsequent thermal self-focusing. Our experiments and supporting 3D simulations in helium demonstrate that magnetization leads to higher thermal gradients, higher peak temperatures, more rapid blast wave development, and beam focusing with an applied field of 12T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-dimensional Thomson scattering measurements of misaligned electron density and temperature gradients and associated Biermann battery produced fields

We present optical Thomson scattering measurements of electron density and temperature in high Mach number laser-driven blast waves in homogeneous gases. Taylor–Sedov blast waves are launched in nitrogen (N2) or helium (He) at pressures between 0.4 mTorr and 10 Torr by ablating a solid plastic target with a high energy laser pulse (10 J, 1012 W/cm2). Experiments are performed at high repetition rate (1 Hz), which allows one-dimensional and two-dimensional Thomson scattering measurements over an area of several cm2 by automatically translating the scattering volume between shots. Electron temperature and density in the blast wave fronts were seen to increase with increasing background gas pressure. Measured electron density and temperature gradients were used to calculate ∂B/∂t∝∇Te×∇ne. The experimentally measured ∂B/∂t showed agreement with the magnetic field probe (B-dot) measurements, revealing that magnetic fields are generated in the observed blast waves via the Biermann battery effect. The results are compared to numerical three-dimensional collisional magnetohydrodynamic simulations performed with FLASH, and are discussed in the context of spontaneous magnetic field generation via the Biermann battery effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗