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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Equation-of-state, sound speed, and reshock of shock-compressed fluid carbon dioxide

Mechanical equation-of-state data of initially liquid and solid CO2 shock-compressed to terapascal conditions are reported. Diamond-sapphire anvil cells were used to vary the initial density and state of CO2 samples that were then further compressed with laser-driven shock waves, resulting in a data set from which precise derivative quantities, including Grüneisen parameter and sound speed, are determined. Reshock states are measured to 800 GPa and map the same pressure-density conditions as the single shock using different thermodynamic paths. The compressibility data reported here do not support current density-functional-theory calculations, but are better represented by tabular equation-of-state models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗
DOE OSTI2021

Sound speed measurements in lithium fluoride single crystals shock compressed to 168 GPa along [100]

The shock wave response of [100] lithium fluoride (LiF) single crystals at high stresses is of long-standing interest due to their extensive use as optical windows in dynamic compression experiments. The report of melting in shock compressed LiF single crystals between 134 GPa and 152 GPa – based on a single sound speed datum [J. Appl. Phys. 117, 045901 (2015)] – was surprising because good optical transmission was previously demonstrated in LiF shock compressed to ~200 GPa [J. Appl. Phys. 116, 033515 (2014)]. To address these apparent differences, we report on plate impact experiments on [100] LiF single crystals shock compressed to 168 GPa. Wave profiles were measured using laser interferometry to determine Hugoniot states and longitudinal sound speeds in shock compressed LiF. The measured Hugoniot states are in good agreement with those measured in previous studies. However, the measured sound speeds presented here show no evidence of melting up to 168 GPa. In particular, the abrupt drop reported previously in sound speed at 152 GPa was not observed in the present work. Furthermore, our results establish a lower bound of 168 GPa for the onset of melting in shock compressed LiF single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

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

Laser-driven collisionless shock acceleration of ions from near-critical plasmas

Here, this paper overviews experimental and numerical results on acceleration of narrow energy spread ion beams by an electrostatic collisionless shockwave driven by 1 μm (Omega EP) and 10 μm (UCLA Neptune Laboratory) lasers in near critical density CH and He plasmas, respectively. Shock waves in CH targets produced high-energy ~50 MeV protons (ΔE/E of ≤30%) and 314 MeV C 6+ ions (ΔE/E of ≤10%). Observation of acceleration of both protons and carbon ions to similar velocities is consistent with reflection of particles off the moving potential of a shock front. For shocks driven by CO 2 laser in a gas jet, ~30 MeV peak in He ion spectrum was detected. Particle-in-cell simulations indicate that regardless of the target further control over its density profile is needed for optimization of accelerated ion beams in part of energy spread, yield and maximum kinetic energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

In-Situ X-Ray Imaging High Strain Rate Compression of Laminate Al-Graphene Composite and Mechanical Property Characterization

Owing to the opaque nature of the laminated structures, traditional high-speed optical camera cannot be used to detect the dynamic process of sub-surface deformation. In this article, we report a study of using high speed X-ray imaging to study the high strain rate deformation in laminated Al structures. We used a Kolsky bar apparatus to apply dynamic compression and a high-speed synchrotron X-ray phase contrast imaging (PCI) setup to conduct the in situ X-ray imaging study. The in situ X-ray imaging captures the shock wave propagation in the laminated structures. After shock compression, we characterized the microstructures by using transmission electron microscopy (TEM), which demonstrates an increase of dislocation density. In conclusion, the micro-pillar compression tests show that the yield strength at 0.2% offset of laminated Al-graphene composite has a significant increase of 67%, from 30 to 50 MPa, compared to laminate Al after shock loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stress Waves in Solids [Slides]

The outline for SSDDS Lecture 3 includes: stress and strain, stress waves, waves at interfaces and boundaries, formation of shock waves, and problems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Novel method to control explosive shock sensitivity: Formulations, experiments and modeling

The microstructure of a heterogeneous high explosive (HE) affects its shock initiation sensitivity and detonation performance. Alteration of the void content and/or void structure (i.e., bulk heating or mechanical damage) therefore changes the shock initiation behavior. Controlling and predicting the change in shock sensitivity after an HE has undergone microstructural changes addresses an important and challenging goal for the design and understanding of novel energetic material formulations. In this study, we aimed to develop HE systems to precisely tune shock sensitivity by a thermal treatment prior to use. Specifically, we incorporated a small fraction (1 wt% or less) of thermally expandable microspheres (TEMs) during the formulation process of various plastic-bonded explosives (PBX). TEMs typically consist of a thermoplastic acrylonitrile shell (10-50 µm diameter), which encapsulates an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol% (irreversibly). Here, we present our progress on shock sensitivity comparisons of HE formulations doped with TEMs after heating and expansion. Experiments, using Composition C-4 thermally cycled to 120°C, showed an increase in shock sensitivity. Additionally, mesoscale modeling revealed that the TEM itself does not act as a hotspot, but instead has a secondary effect on the response of nearby voids. We conclude that TEMs indeed do provide a method to tune shock sensitivity. Further, our results emphasize the non-uniformity of shock waves at the mesoscale, and that upstream defects influence the reactivity of downstream defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shock behavior of materials

The physics governing the response of materials to shock loading has many applications in various fields, ranging from Earth and planetary science to engineering and ballistics. Elucidating the behavior of materials experiencing the passage of shock waves motivates active research worldwide, including experimental developments with in situ and time-resolved measurements, modern computational approaches, and multi-scale simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

A Computational Study of the Thermodynamic Conditions Leading to Autoignition in Nanosecond Pulsed Discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results toward resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model (BRM). This can afford the use of very detailed chemical kinetic information. The full nonequilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

42 ENGINEERING↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker

Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes hot enough that the x-ray opacity falls significantly. Additionally, if the system includes a Lagrangian marker such that the mass within the marker is known, this additional information can be used to constrain the opacity as well as the Hugoniot state. In the limit that the opacity changes only on shock heating, and not significantly on subsequent isentropic compression, the opacity of the shocked material can be determined uniquely. More generally, it is necessary to assume the form of the variation of opacity with isentropic compression or to introduce multiple marker layers. Alternatively, assuming either the equation of state or the opacity, the presence of a marker layer in such experiments enables the non-assumed property to be deduced more accurately than from the radiographic density reconstruction alone. An example analysis is shown for measurements of a converging shock wave in polystyrene at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Listening for the Landing: Seismic Detections of Perseverance's Arrival at Mars With InSight

The entry, descent, and landing (EDL) sequence of NASA's Mars 2020 Perseverance Rover will act as a seismic source of known temporal and spatial localization. We evaluate whether the signals produced by this event will be detectable by the InSight lander (3,452 km away), comparing expected signal amplitudes to noise levels at the instrument. Modeling is undertaken to predict the propagation of the acoustic signal (purely in the atmosphere), the seismoacoustic signal (atmosphere-to-ground coupled), and the elastodynamic seismic signal (in the ground only). Our results suggest that the acoustic and seismoacoustic signals, produced by the atmospheric shock wave from the EDL, are unlikely to be detectable due to the pattern of winds in the martian atmosphere and the weak air-to-ground coupling, respectively. However, the elastodynamic seismic signal produced by the impact of the spacecraft's cruise balance masses on the surface may be detected by InSight. The upper and lower bounds on predicted ground velocity at InSight are 2.0 × 10 -14 and 1.3 × 10 -10 m s -1 . The upper value is above the noise floor at the time of landing 40% of the time on average. The large range of possible values reflects uncertainties in the current understanding of impact-generated seismic waves and their subsequent propagation and attenuation through Mars. Uncertainty in the detectability also stems from the indeterminate instrument noise level at the time of this future event. A positive detection would be of enormous value in constraining the seismic properties of Mars, and in improving our understanding of impact-generated seismic waves.

58 GEOSCIENCES↗
DOE OSTI2020

(U) PRAD0697 & PRAD0698: Complex Loading of CeO₂ Powder

Cerium(IV) oxide (CeO₂) powder is shock compressed using the Precision High Energy-density Liner Implosion eXperiment (PHELIX) platform. Experimental results are compared against several modeling approaches. Compaction behavior is best captured with a P-∝ model, which calculates CeO₂ powder bulk densities within 80-99% of experimental values but overpredicts densi cation at the cylindrical target's outer radius and center by up to 20%. Preliminary calculations suggest that accuracy could be increased with the inclusion of a coupled strength model. Several common computational modeling approaches for the shock compression response of granular materials and the magnetohydrodynamic (MHD) force upon the impactor/liner in pulsed power compression experiments are investigated and analyzed for their validity. The Bi-linear Ramp, P-∝ PACXP, and P-∝ Menikoff-Kober continuum compaction models are calibrated to planar impact Hugoniot data for CeO₂ powder and used to predict the powder's shock compaction response under non-planar shock wave compression. MHD calculations of the PHELIX pulsed power driver are performed using an idealized resistor-inductor-capacitor (RLC) circuit calibrated to previous experiments. All simulations are performed using the LANL code FLAG. Two validation experiments are computationally designed using the calibrated compaction and circuit models, executed using the PHELIX platform on CeO₂ targets with initial porous densities of 3.95 and 4.03 g/cm³, measured with proton radiography, and analyzed against the model predictions. The two P-∝ models more accurately describe CeO₂ powder densi cation than the Bi-linear Ramp model. However, the two P-∝ models overpredict bulk density of the shock compressed CeO₂ powder by up to 20% when the appropriate impact velocities are applied. MHD calculations for both validation experiments underpredict liner impact velocities by 4-11% when using the idealized RLC circuit model calibrated to previous experiments. Compensating underpredictions of impact velocity and overpredictions of powder densication lead to a false accuracy in pre-shot calculations compared to experimental data. To improve correlation between simulations and experiments, the following improvements are suggested: 1. A coupled strength model for CeO₂ powder that updates strength as a function of porosity and applied stress. 2. An improved MHD circuit model that more accurately captures the PHELIX machine.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A study of explosive-induced fracture in polymethyl methacrylate (PMMA)

The fracture response of geologic materials is of interest for applications, including geothermal energy harnessing and containment of underground explosions. To better understand the explosively induced fracture response of geomaterials, polymethyl methacrylate (PMMA) was used as a transparent rock surrogate to allow imaging of internal shock propagation and fracture growth processes. Experiments were conducted using high-speed shadowgraphy and photon Doppler velocimetry (PDV), which were compared to numerical simulations. Experiments measured fractures produced in 304.8mm × 304.8mm × 304.8mm PMMA cubes with two simultaneously initiated detonators. The cubes were subjected to varying amounts and directions of externally applied uniaxial stresses, including no stress, 2 MPa stress, and 20 MPa stress. The fracture radius as a function of time was extracted from the high-speed videos. Post-test images of the PMMA cubes aided in the determination of three-dimensional effects not directly imaged by the cameras. The surface velocity history and the shock response captured in PDV and the high-speed videos were compared to the simulated explosive-induced shock response. The simulation results indicate that the shock drives the fracture for the first 20 μs corresponding to a fracture radius of approximately 15 mm in the experiments. The gas-driven fracture extent was estimated analytically using an equilibrium stress distribution calculated after the shock wave propagation through the sample. Reduction in the gas pressure due to the leakage of the explosive products through the crack as a function of time was accounted for. In conclusion, the estimated fracture lengths were in agreement with the experimentally observed fracture lengths.

15 GEOTHERMAL ENERGY↗