Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Isostructural Phase Transition of Fe 2 ⁢O 3 under Laser Shock Compression

We present in situ x-ray diffraction and velocity measurements of Fe 2⁢ O 3 under laser shock compression at pressures between 38–122 GPa. None of the high-pressure phases reported by static compression studies were observed. Instead, we observed an isostructural phase transition from 𝛼−Fe 2⁢ O 3 to a new 𝛼′−Fe 2 ⁢O 3 phase at a pressure of 50–62 GPa. The 𝛼′−Fe 2⁢ O 3 phase differs from 𝛼−Fe 2 ⁢O 3 by an 11% volume drop and a different unit cell compressibility. We further observed a two-wave structure in the velocity profile, which can be related to an intermediate regime where both 𝛼 and 𝛼′ phases coexist. Furthermore, density functional theory calculations with a Hubbard parameter indicate that the observed unit cell volume drop can be associated with a spin transition following a magnetic collapse.

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

Measurement of turbulent velocity and bounds for thermal diffusivity in laser shock compressed foams by x-ray photon correlation spectroscopy

Experimental benchmarking of transport coefficients under extreme conditions is required for validation of differing theoretical models. To date, measurement of transport properties of dynamically compressed samples remains a challenge with only a limited number of studies able to quantify transport in high pressure and temperature matter. x-ray photon correlation spectroscopy utilizes coherent x-ray sources to measure time correlations of density fluctuations, thus providing measurements of length and timescale-dependent transport properties. Here, we present a first-of-a-kind experiment to conduct x-ray photon correlation spectroscopy in laser shock compression experiments. We report measurement of the turbulent velocity in the wake of a laser driven supersonic shock and place an upper bound on thermal diffusivity in a solid density plasma on nanosecond timescales.

High-energy-density plasmas↗
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

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

Shock compression of diamond single crystals to 120 GPa: Refractive index and nonlinear photoelasticity

The optical response of transparent solids at extreme conditions is important for both fundamental science and many applications. Strong transparent solids are of particular interest for use as optical windows in dynamic compression experiments. Due to diamond’s exceptional strength and optical properties, laser-driven shock experiments and plate impact experiments were carried out to examine the diamond optical response for shock wave compression along two different crystal orientations. Using laser interferometry at 532 nm and 1550 nm wavelengths, optical transparency was observed and refractive indices were determined for [100] diamond at stresses up to 119 GPa and for [111] diamond at stresses up to 87 GPa. From these results, the nonlinear photoelastic response for [100] and [111] diamond was determined, revealing significant dependence on both crystal orientation and laser wavelength. To enable [100] diamond as an interferometry window in dynamic compression experiments, the requisite window corrections were determined for 532 nm and 1550 nm wavelengths. Because of diamond’s excellent x-ray transparency, the present findings will be particularly useful for incorporating [100] diamonds as windows in dynamic compression experiments involving x-ray diffraction or other x-ray diagnostics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1985

Shock Compression of Liquid Helium to 56 GPa (560) Kbar

Shock-wave data are presented for liquid helium which has been compressed to densities up to five times greater than the normal liquid. The helium was heated to temperatures up to 21,000 K, while the maximum pressure attained was 56 GPa. The properties of helium and hydrogen are important for modeling the giant planets Saturn and Jupiter where these elements are the major constituents. Conditions on Saturn are of particular interest because studies have suggested that this planet has an internal energy source which is associated with unmixing and gravitational separation the hydrogen-helium fluid at pressures below 1 TPa. The existence of this phase transition depends very sensitively on the hydrogen and helium equation of state. In the experiments, strong shock waves were generated by the impact of planar projectiles into cryogenic specimen holders.

Nellis, W. J.↗
NASA NTRS1984

Shock compression of liquid helium to 56 GPa (560 kbar)

Shock-wave data are presented for liquid helium which has been compressed to densities up to five times greater than the normal liquid. The helium was heated to temperatures up to 21,000 K, while the maximum pressure attained was 56 GPa. The properties of helium and hydrogen are important for modeling the giant planets Saturn and Jupiter where these elements are the major constituents. Conditions on Saturn are of particular interest because studies have suggested that this planet has an internal energy source which is associated with unmixing and gravitational separation of the hydrogen-helium fluid at pressures below 1 TPa. The existence of this phase transition depends very sensitively on the hydrogen and helium equation of state. In the experiments, strong shock waves were generated by the impact of planar projectiles into cryogenic specimen holders.

Nellis, W. J.↗
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2026

Phase Transitions of Eutectic High Entropy Alloy AlCoCrFeNi 2.1 Under Shock Compression

High entropy alloys (HEAs) are a new class of metals that exhibit unique mechanical performance. Among HEAs, additively manufactured eutectic high entropy alloys (AM EHEAs) have recently emerged as candidate materials for use in extreme conditions due to their simultaneous high strength and ductility. However, the deformation and structural evolution of AM EHEAs under conditions of high pressure have not been well characterized, limiting their use in extreme applications. Dynamic compression experiments and molecular dynamics simulations are presented to study the structural evolution of AM EHEA AlCoCrFeNi 2.1 when compressed to pressures up to 400 GPa. In situ X-ray diffraction measurements capture the appearance of face-centered cubic and body-centered cubic phases at different pressure conditions, with pure- and mixed-phase regions. Furthermore, understanding the phase stability and structural evolution of the AM EHEA offers new insights to guide the development of high-performance complex materials for extreme conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Shock compression and adiabatic release of lunar fines from Apollo 17

An experimental investigation was conducted with the objective to obtain quantitative bounds, in terms of shock pressure and hence meteoroid impact velocity, concerning the conditions required to compact and lithify lunar fines. The measured pressure-particle velocity release states provide a basis for the determination of the approximate values of shock pressure associated with various postshock volumes and temperatures concomitant with solid-state vitrification and thermal melting. The implications of the obtained results for the study of regolith evolution are discussed.

Ahrens, T. J.↗
DOE OSTI2024

Pressure-induced transformation of Nb 2 O 5 under shock compression from first principles

Ab initio molecular dynamics (AIMD) simulations were carried out to investigate the equation of state of Nb 2 O 5 and its pressure-density relationship under shock conditions. Here, the focus of this study is on the monoclinic B-Nb 2 O 5 (C2/c) polymorph. Enthalpy calculations from AIMD trajectories at 300 K show that the pressure-induced transformation between the thermody-namically most stable crystalline monoclinic parent phase H-Nb 2 O 5 (P2/m) and B-Nb 2 O 5 occurs at ~1.9 GPa. This H→B transition is energetically more favorable than the H→L(Pmm2) pressure-induced transition recently observed at ~5.9-9.0 GPa. The predicted shock properties of Nb 2 O 5 polymorphs are also compared to their Nb and NbO 2 counterparts to assess the impact of niobium oxidation on shock response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Optical transparency of LiF (100) shock compressed to ∼360 GPa

High-purity [100] lithium fluoride (LiF) is the most widely used optical window in dynamic compression experiments due to its wide bandgap and well-characterized mechanical response. Recent plate-impact experiments established the [100] LiF Hugoniot to ∼230 GPa and demonstrated shock-induced melting onset at 182 GPa with complete melting by 195 GPa; theoretical models predict LiF optical transparency to nearly 900 GPa. To experimentally examine the optical transparency of [100] LiF at higher pressures and in the liquid state, laser-driven shock experiments were performed at peak stresses ranging from 223 to 363 GPa. Optical response was examined by measuring the particle velocity histories at the Kapton/LiF interface using laser interferometry at 532 and 1550 nm wavelengths; in-material particle velocities were obtained using established refractive-index corrections. Continuous photonic Doppler velocimetry fringes were observed across the entire stress range, demonstrating that LiF remains transparent to 1550 nm light throughout the multi-megabar regime investigated. At 532 nm, fringe visibility depended on the reflector coating: aluminum mirrors provided signals to ∼235 GPa, while gold mirrors extended this limit to 270.5 GPa, indicating that the shorter-wavelength response is likely sensitive to experimental configuration rather than to the loss of LiF transparency. Continued optical transparency to at least 360 GPa indicates that shock-melted LiF does not display bandgap closure over the stress range explored. Furthermore, these results provide direct experimental constraints on the high-pressure optical response and establish LiF (100) as a robust optical window material for laser-driven dynamic compression experiments approaching 400 GPa.

Renganathan, P. [Argonne National Laboratory (ANL)↗
NASA NTRS1985

Hybrid approach to steady transonic normal shock-compressible laminar boundary layer interactions over airfoils with suction

Transonic airfoil flow is to a large degree affected by viscous-inviscid interactions. Among them a key role is played by the boundary layer interaction with the shock wave embedded in the flow field and the interaction of the boundary layer with the sustained adverse pressure gradients. The effects of these interactions can be controlled or suppressed by introducing surface mass transfer such as suction. This investigation deals strictly with the study of shock-laminar boundary layer interactions including distributed mass transfer. The present calculations agree with experiments on a swept LFC wing at low speeds which indicates that with a suitable choice of the extent of the airfoil surface over which suction is applied and adjustment of the suction velocity, full-chord laminar flow can be maintained and separation can be prevented completely.

Ram, R. B.↗
NASA NTRS1973

Shock compression of a recrystallized anorthositic rock from Apollo 15

Hugoniot measurements on 15,418, a recrystallized and brecciated gabbroic anorthosite, yield a value of the Hugoniot elastic limit (HEL) varying from 45 to 70 kbar as the final shock pressure is varied from 70 to 280 kbar. Above the HEL and to 150 kbar, the pressure-density Hugoniot is closely described by a hydrostatic equation of state constructed from ultrasonic data for single-crystal plagioclase and pyroxene. Above 150 kbar, the Hugoniot states indicate that a series of one or more shock-induced phase changes are occurring in the plagioclase and pyroxene. From Hugoniot data for both the single-crystal minerals and the Frederick diabase, we infer that the shock-induced high-pressure phases in 15,418 probably consists of a 3.71 g/cu cm density, high-pressure structure for plagioclase and a 4.70 g/cu cm perovskite-type structure for pyroxene.

Ahrens, T. J.↗
DOE OSTI2026

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗