Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock-compression”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Direct observation of diamond formation in a shock-compressed high explosive

Understanding the formation timescale and structure of carbonaceous reaction products is critical for modeling the high-pressure equation-of-state of organic materials. We use the National Ignition Facility to shock-compress polycrystalline TATB (C 6 H 6 N 6 O 6 ) samples to ~70–130 GPa and ~4000–5500 K, employing in situ nanosecond X-ray diffraction to probe reaction products and velocimetry to measure transmitted compression wave profiles. Our diffraction data is consistent with the formation of diamond over timescales less than ~60 ns. This represents carbon condensation from a molecular explosive on timescales three times faster than previously reported and the earliest observation of diamond produced from reacting TATB. Reactive flow simulations with explicit chemistry reproduce the observed temporal structure within wave profiles to inform the distribution of P-T states. These findings provide direct evidence of ultrafast diamond formation in a reactive system at extreme conditions and provide new constraints for models of shock and detonation chemistry.

Clarke, Samantha M. [Lawrence Livermore National L↗
DOE OSTI2025

Chemical decomposition in shock-compressed 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals: Time-resolved Raman spectroscopy

Understanding the evolution of chemical decomposition—at the molecular level—in shock-compressed insensitive high explosive (IHE) single crystals is an important need. Toward this goal, time-resolved Raman spectra were measured in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals shock compressed to stresses above 20 GPa. At 22 GPa and higher stresses, the Raman peak intensities were significantly reduced, showing the onset and extent of chemical changes. At 33 GPa, no Raman peaks were observed, suggesting complete decomposition of FOX-7. The 22 GPa onset from Raman data is in marked contrast to the much higher (32 GPa) onset determined previously from wave profile measurements—showing the significantly greater sensitivity of Raman measurements for gaining insight into shock-induced molecular-level changes in IHE single crystals. The close match between the complete loss of Raman peaks at 33 GPa and the 32 GPa reaction threshold determined from continuum data makes a good case that the two results are related and they reflect significant energy release due to decomposition.

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

Mechanics of pore array collapse and interaction in shock-compressed polymethyl methacrylate (PMMA)

Recent studies on dynamic pore collapse have revealed significant development of shear localization, which can lead to material failure in porous structures and hot spot generation in energetic materials. These findings have dramatically improved the understanding of failure mechanisms during pore collapse but also prompt further investigation of realistic porous materials. In particular, porous media consist of many pores and porous networks. Even in low-porosity materials, pores can form in close proximity during the manufacturing process, leading to the critical question of pore–pore interaction during collapse under dynamic loading conditions. This study investigates, via plate impact experiments coupled with high-speed internal digital image correlation and shadowgraphy techniques, the collapse of two pores in shock-compressed PMMA at stresses between 0.4 and 1 GPa. The results of these experiments provide new insights into shear localization in pore collapse, in addition to distinct interactions between pores. Shadowgraphy measurements reveal novel, direct visualization of shear band development and crack evolution from pore surfaces. Spacing between adiabatic shear bands is measured over a range of impact stresses and is predicted accurately by the Grady–Kipp model. Pore interactions are found to effect a transition in the impact stress threshold at which different failure mechanisms initiate and are also found to possibly influence preferential sites for shear cracking. Throughout the study, numerical and theoretical models are leveraged to understand shear localization behavior. The role of baroclinicity and wave interactions between the pores is used to elucidate interaction mechanisms between pores.

Lawlor, Barry P. [California Institute of Technolo↗
DOE OSTI2026

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2024

Sound speed measurements in shock-compressed cemented Tungsten carbides at pressures up to 100 GPa

To gain insights into thermodynamic states attained during shock compression of cemented tungsten carbide with 3.7 wt.% cobalt binder, we present results of longitudinal sound (release wave) speed measurements and their analysis at peak stresses up to 100 GPa (volumetric compression ratio ~ 15%). The sound speeds are determined using front-surface impact and release-wave overtake plate impact experimental configurations using laser interferometry. The measured sound speed data along with estimates for bulk sound speeds obtained using the fourth-order Birch-Murnaghan EoS and thermodynamics are used to determine the longitudinal moduli and shear moduli of shocked tungsten carbide at the various peak compression states attained in the experiments. Here, the longitudinal sound speeds were found to increase linearly with volume compression ratio from 6.97 ± 0.010 km/s at ambient conditions to 8.26 ± 0.156 km/s at a volume compression ratio of ~ 15%. The corresponding longitudinal elastic moduli also increase nearly linearly with the volume compression ratio but remain consistently lower than their theoretical predictions based on continuum models with no damage. Also, the sensitivity of shear moduli to pressure, as predicted by the Steinberg-Guinan model, is reduced substantially and the shear moduli of cemented WC with 3.7 wt.% Co remains nearly constant at ~ 310 GPa at the various peak compression stress states investigated in the present study.

36 MATERIALS SCIENCE↗
DOE OSTI2026

High-pressure melt dynamics in shock-compressed titanium

In this work we study the high-pressure melting behavior of titanium using laser-driven shock compression with in situ femtosecond x-ray diffraction and molecular-dynamics simulations based on a machine-learned interatomic potential. The MD simulations predict the solid-liquid coexistence on the Hugoniot in the ∼111−124GPa range. Experimentally, we observe the first evidence of liquid at 86 GPa. We also observe pronounced microstructural changes with pressure, with strong grain refinement associated with the emergence of liquid, within the solid-liquid coexistence (∼110−126GPa). Above 126 GPa, we observe the persistence of residual levels of highly textured crystalline Ti to ∼180GPa, well above the expected melt completion pressure. We discuss the accuracy that current laser-shock experimental platforms have at determining the melt onset and completion pressures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2025

Shock compression of additively manufactured high solids loaded polymer composites with aligned porosity

The effects of aligned porosity on the shock-compression response of additively manufactured (AM) composite structures with collinear and log-cabin geometries are investigated. The composites contain two inorganic and two organic particle populations and a UV-curable polymer binder. X-ray phase contrast imaging (X-PCI) is used as an in-material diagnostic in plate-impact experiments performed at the Dynamic Compression Sector of the Advanced Photon Source at the Argonne National Laboratory. Macroscale shock and particle velocities and microscale pore collapse velocities are measured from X-PCI images obtained from samples impacted in various orientations relative to the AM build directions. The linear shock and particle velocity relation for the structures with collinear build geometry shows directional differences, with the build direction exhibiting the highest slope and the filament print direction having the shallowest slope. Composites with log-cabin geometry, due to their slightly lower porosity and smaller sized pores, show higher shock velocities than the collinear structures at the same particle velocity. Microscale directionality effects are also observed, with the pore collapse velocity dependent on the pore configuration. Specifically, the pore collapse velocity is highest for pores elongated along the impact direction which is the case for the colinear structures impacted in the filament print direction, which is consistent with the directionality trend observed in the data for the shock and particle velocity equation of state. The pore collapse velocity approaches ∼80 % of the shock velocity and is several times the average bulk particle velocity, suggesting a highly hydrodynamic mechanism of pore collapse. Overall, both the macroscale response as measured by the shock and particle velocity relationship, and the microscale response as measured by the pore collapse velocity, are found to be anisotropic and dependent on the AM composite structure.

Wagner, K. B. [Georgia Institute of Technology, At↗
DOE OSTI2024

Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO 3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of theO-O and Mg-Si networks beyond Earth’s interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets withmasses nearly three times that of Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stability of the fcc phase in shocked nickel up to 332 GPa

Despite making up 5-20 wt.% of Earth’s predominantly iron core, the melting properties of elemental nickel at core conditions remain poorly understood, due largely to a dearth of experimental data. We present here an in situ X-ray diffraction study performed on laser shock-compressed samples of bulk nickel, reaching pressures up to ~ 500 GPa. Hugoniot states of nickel were targeted using a flat-top laser drive, with in situ X-ray diffraction data collected using the Linac Coherent Light Source. Rietveld methods were used to determine the densities of the shocked states from the measured diffraction data, while peak pressures were determined using a combination of measured particle velocities, shock transit times, hydrodynamic simulations, and laser intensity calibrations. We observed solid compressed face-centered cubic (fcc) Ni up to at least 332 ± 30 GPa along the Hugoniot—significantly higher than expected from the majority of melt lines that have been proposed for nickel. We also bracket the partial melting onset to between 377 ± 38 GPa and 486 ± 35 GPa.

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

Shock compression of FeOOH and implications for iron-water interactions in super-earth magma oceans

Iron(Fe)-water reactions in a magma ocean can influence water storage and density of planets. These reactions can form Fe-O-H phases, whose density, melting, and electronic properties at planetary interior conditions are important for informing planetary models. Here, we study natural goethite (α-FeOOH) that is shock-compressed along its principal Hugoniot. Analysis of our velocity interferometer system for any reflector (VISAR) results extends the equation of state to over 800 GPa. X-ray diffraction and VISAR reflectivity results indicate the onset of melting occurs at ~95 GPa with complete melting by 166 GPa, which may be relevant to low seismic velocity anomalies observed above the core-mantle boundary. Analysis of X-ray emission spectroscopy results up to 285 GPa shows the spin crossover of Fe, with dominantly low spin Fe above ~265 GPa in the melt, supporting formation of dense basal magma oceans in terrestrial planets. Using our measured FeOOH densities, we model planetary interiors up to 10 Earth masses. Assuming FeOOH forms via iron-water reactions, the radius decreases by up to 28%, while the density increases by up to 165% compared to the unreacted cases, providing an avenue to investigate water storage and evolution in super-Earths and sub-Neptunes.

Zhang, Yanyao [Stanford Univ., CA (United States)]↗
DOE OSTI2026

Compression rate effects on the Bi-III stability field in dynamic diamond anvil cell XFEL studies of bismuth

This study aims to resolve conflicting observations of bismuth’s high-pressure phase transitions across static, intermediate, and shock compression regimes. We probed its high-pressure structural sequence using the dynamic diamond anvil cell with and time-resolved X-ray diffraction with microsecond resolution at an X-ray free-electron laser. At room temperature and pressures up to 20 gigapascals, bismuth evolves through the same structural sequence previously identified under static compression. However, the transformation behavior differs in important ways under dynamic loading. In particular, the transition from the incommensurate intermediate-pressure phase to the high symmetry high pressure phase begins at pressures about 2 to 4 gigapascals lower than in static experiments, and the onset of this transformation depends on the compression rate. As a result, the stability field of the intermediate-pressure phase is reduced under rapid compression. Our results suggest that at sufficiently fast loading rates the intermediate phase may be bypassed entirely, consistent with previous shock-compression observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Shock compression of single-crystal austenitic FeCr18Ni12.5 stainless steel to 60 GPa

We measured the austenitic FeCr18Ni12.5 stainless steel Hugoniot as a function of crystallographic direction to approximately 60 GPa. We shock-compressed FeCr18Ni12.5 samples oriented along $\langle$100$\rangle$, $\langle$110$\rangle$, and $\langle$111$\rangle$ to mean stresses ranging 30.5–58.1 GPa via Ta plate impact in a large-bore powder gun and measured the free-surface velocities with laser interferometry. We unambiguously observed the largest post-shock free-surface velocity along $\langle$100$\rangle$ in each experiment, which consequently produced the lowest shock velocity along that orientation. However, the propagation of experimental uncertainties through the impedance matching scheme used to compute the shock velocity produced sufficient uncertainty overlap to preclude definitive conclusion of Hugoniot anisotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of warm dense matter by combining static and dynamic compression

Shock-compressing material samples that are precompressed with diamond-anvil cells (DACs) enables benchmarking equation of state measurements over a wide range of high-pressure conditions, in particular, warm dense matter relevant to planetary science, astrophysics, and inertial confinement fusion. These targets provide optical access to the shocked sample, so optical velocity interferometry and pyrometry are typical diagnostics. We present a scaling for the coupling of the laser drive to the shock pressure in DACs based on over 150 prior experiments that can be used to design future experiments. For initial static pressures ≤ 5 GPa, using a sapphire window mitigates the effect of photoionization compared to diamond windows and allows the measurement of shock pressures up to 1 TPa with optical diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Calibration and characterization of the line-VISAR diagnostic at the HED-HIBEF instrument at the European XFEL

In dynamic-compression experiments, the line-imaging Velocity Interferometer System for Any Reflector (VISAR) is a well-established diagnostic used to probe the velocity history, including wave profiles derived from dynamically compressed interfaces and wavefronts, depending on material optical properties. Knowledge of the velocity history allows for the determination of the pressure achieved during compression. Such a VISAR analysis is often based on Fourier transform techniques and assumes that the recorded interferograms are free from image distortions. In this paper, we describe the VISAR diagnostic installed at the HED-HIBEF instrument located at the European XFEL along with its calibration and characterization. It comprises a two-color (532, 1064 nm), three-arm (with three velocity sensitivities) line imaging system. We provide a procedure to correct VISAR images for geometric distortions and evaluate the performance of the system using Fourier analysis. We finally discuss the spatial and temporal calibrations of the diagnostic. As an example, we compare the pressure extracted from the VISAR analysis of shock-compressed polyimide and silicon.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗
DOE OSTI2026

Temperature and conductivity in shock compressed bridgmanite MgSiO 3 up to 2 TPa

The melting behavior and transport properties of MgSiO 3 at multi-megabar pressures remain poorly constrained despite their importance for high-pressure silicate physics. Here we report the first direct measurements of temperature and optical reflectivity in shock-compressed bridgmanite (MgSiO 3 ) using laser-driven decaying shock compression combined with velocimetry and optical pyrometry. Temperature and reflectivity data spanning approximately 4000–60 000 K were used to constrain the MgSiO 3 melting curve and to infer its electrical conductivity. We find that the MgSiO 3 melting curve becomes shallower than that of iron above 400 GPa, yielding lower melting temperatures in planetary mantles than predicted by several previous theoretical estimates. Across the solid-liquid transition, the inferred electrical conductivity increases significantly, reaching ∼2000 Ω cm −1 . These results provide experimental benchmarks for theoretical models of silicate melting and transport under extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗