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

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Effects of Turbulence on Heliosheath Ions and Implications for Energetic Neutral Atoms

The distribution of ions in the heliosheath—the region between the heliospheric termination shock and the heliopause—is important for understanding remote observations of energetic neutral atoms (ENAs). The ion distributions were estimated previously based on hybrid simulations of the heating and evolution of solar wind and interstellar pickup ions across the solar wind termination shock, but these estimates only provide the distributions near the shock. In this work, we use self-consistent hybrid kinetic simulations to investigate the effects of turbulence on ion distributions in the heliosheath. The simulations are compared against Voyager observations, constraining the feasible amplitude and compressibility of turbulence. We find that the heating due to turbulent dissipation can lead to a significant increase in the temperature of thermal solar wind ions. Both turbulent velocity fluctuations and the heating of solar wind ions increase the charge-exchange source for ENAs at low energies (around 100 eV), where current ENA models underpredict observations by more than an order of magnitude. However, the effects of turbulence are likely not strong enough to fully explain these discrepancies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamic response of 17-4 stainless steel as a function of manufacturing method and heat treatment

We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.

Compressive stress↗
DOE OSTI2025

Impact Testing Theories and Techniques Employed at the TA-III Mechanical Shock Complex

We present theories and techniques to support safe, predictable, and efficient impact testing at the Sandia National Laboratories TA-III Mechanical Shock Complex (MSC), located in Building 6570. The two testing instruments used at the MSC are the 20-inch actuator, which is a launcher that propels sleds down an indoor track, and an indoor 6-in bore compressed-gas powered cannon. The theories and techniques presented derive from basic mechanical engineering principles but are significantly informed and validated by testing experience at the MSC.

42 ENGINEERING↗
DOE OSTI2026

Improving indirect-drive ablator performance through pulse shaping and material selection in double shell implosions

The double shell path to inertially confined volumetric thermonuclear burn relies on efficient kinetic energy transfer from an outer ablator shell to an inner shell containing DT fuel at implosion velocities near 250 km/s. In order to experimentally realize adequate compression, target imperfections and assembly features must be tightly controlled. Pulse shape and ablator material are critical design choices that impact all aspects of the implosion. In this article, we evaluate three designs: an Al ablator with single-shock radiation drive, a CH ablator with single-shock drive, and a CH ablator with reduced-adiabat 2-shock drive. All designs use the same W inner (pusher) shell overcoated with a glow discharge polymer tamper and filled with liquid-density DT. One-dimensional HYDRA simulations revealed that the plastic ablator design required the use of the 2-shock drive to match the performance of the Al design. Two-dimensional simulations resolving broadband surface roughness up to modes of several hundred placed only on the ablator showed plastic to greatly outperform Al. The plastic ablator designs did exhibit higher levels of inner shell instability growth when surface roughness was only placed at the pusher/tamper interface, which was mitigated with higher foam cushion densities. A unique and additional benefit of a plastic ablator is also discussed, where heating of the ablator by Au M-band radiation is predicted to rapidly expand the hemi-shell surfaces, closing the equatorial joint gap ahead of the first shock. Taken together, the properties of the plastic ablator for indirect-drive double shell implosions appear promising in removing several long-standing fabrication challenges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effect on compression of lowering the design adiabat in the SQ-n campaign

A set of experiments in the SQ-n campaign attempts to increase the compression in layered implosions by reducing the design adiabat from 3.0 down to ∼ 2.6. The experiments utilize nanocrystalline high density carbon ablators driven by a ramped laser pulse drive. Contrary to expectations, the lower design adiabat layered implosion reaches the same level of compression, within the measurement uncertainty, as the previous set of SQ-n layered implosions at the higher design adiabat. Here, we examine potential explanations for this result, including incomplete melting or refreezing of the ablator after the first shock passes into the fuel layer, reproducibility, and dopant areal density. We also discuss additional potential paths to reaching lower adiabat including increased dopant areal density, ramped dopant profiles, and increased initial capsule size which are also based on experimental campaign results from the National Ignition Facility.

Physics - Plasma physics↗
DOE OSTI2025

Improving high temperature resilience of fiber sensor embedded smart components through laser shock peening

This study explores the use of laser shock peening (LSP) to enhance material properties and high-temperature performance of fiber-sensor-fused smart parts fabricated by additive manufacturing (AM) methods. Using embedded fiber sensors as distributed strain gauges, the study demonstrates that LSP can induce compressive strains of up to 130 µε on fiber embedded 1-mm below metal surfaces. The electron backscatter diffraction (EBSD) analysis shows that, with optimized LSP parameters, the metallic matrix undergoes substantial microstructural refinement, resulting in denser structures. Thermal cycling tests showed that the LSP process can increase fiber slippage temperatures by more than 50 °C. This work shows that the LSP process is an effective room-temperature process for enhancing both surface quality and increasing fiber slippage threshold under both thermal and mechanical stress.

Zhong, Shuda [University of Pittsburgh, PA (United↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2025

Dynamic compression effects of H 2 ⁡O in a dynamic diamond anvil cell: Origin of metastable ice VII and its crystal growth kinetics

We report on the structural verification of metastable ice VII solidifying in the phase space of ice VI at 1.80 GPa at room temperature. Using time-resolved (TR) x-ray diffraction and TR ruby luminescence paired with high-speed microphotography utilizing a dynamic diamond anvil cell, an initial compression rate range from 0.12 to 95.84 GPa/s was explored. The solidification pressure of metastable ice VII has a potential sigmoidal dependence upon compression rate with a turnover compression rate of ∼80 GPa/s. The preferred crystallization of ice VII in the stability field of ice VI is due to the increased nucleation rate of ice VII over ice VI at 1.77 GPa that is driven by the surface energy difference between the liquid and solid phases along with the change in Gibbs free energy of solidification. The dynamic pressure-volume–compression behaviors of ice phases (VI and VII) show a lattice stiffening in both phases, especially during the compression loading. It is also found that the compression rate greatly affects the solid-solid phase transition between ice VI and VII but does not affect the liquid-solid transition between water and ice VI as much. Lastly, a third phase transition was found to occur after metastable ice VII transforms into high-density amorphous (HDA) ice, which could be a disordered hydrogen-bonded network configuration of ice VII forming out of HDA ice facilitated by the decoupling of the oxygen movement and reorientation of the H 2⁡ O molecule. These results demonstrate the complexity of a seemingly simple molecule H 2⁡ O, how it can readily change its static properties with the modification of (de)compression rate, and highlight the need to use multiple TR structural and spectroscopic probes at higher time resolutions to realize the most comprehensive understanding.

Chemical bonding↗
DOE OSTI2025

Sound Velocity and Grüneisen Parameter in Shock‐Melted Silica at Deep Earth Pressures

Silica is a primary component of rocky planet interiors and its melt properties are important for understanding planetary formation and differentiation, magma oceans, and the deep mantle. Although well understood in the solid state, the high-pressure behavior of liquid silica is poorly constrained at lower mantle pressures. Using laser interferometry to measure shock wave profiles, we report measured stress-density states and longitudinal sound speeds in shock-synthesized stishovite, from fused silica staring material, and across the solid-liquid phase boundary up to 154 GPa. Our results constrain completion of melt at 80 GPa and show that at pressures relevant to the deep mantles of Earth-sized, rocky planets, the Grüneisen parameter for liquid silica increases with compression. This finding is consistent with a continuous increase in Si-O coordination above six for liquid silica at core-mantle boundary relevant pressures and temperatures.

Ocampo, I. K. [Princeton University, NJ (United St↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonseeded linewise temperature measurements by resonantly ionized photoemission thermometry in a Mach 4 Ludwieg tube

A one-dimensional (1D) thermometry using oxygen-tagging resonantly ionized photoelectron thermometry (O 2 RIPT) was employed to investigate thermal gradients within a Mach 4 Ludwieg tube. The Ludwieg tube is pulsed with a test duration of approximately 100 ms, providing a cold supersonic flow at Mach 4 ideal for studying aerothermal effects. This study focused on measuring freestream temperatures, capturing shock-induced heating behind a detached bow shock from a blunt cylinder, and resolving sharp temperature variations across a bow shock generated by a cylinder. The O 2 RIPT technique produced strong emission signals extending approximately 4 cm long, demonstrating its capability for precise temperature measurements in high-speed wind tunnel environments. The results confirm that O 2 RIPT is well-suited for applications in large-scale aerodynamic testing facilities, particularly in regions with strong compression effects, enabling the resolution of sharp thermal gradients. This method presents a promising solution for thermometry in dynamic flow conditions relevant to various experimental ground-test facilities.

McCord, Walker (ORCID:0000000179747550)↗
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 OSTI2025

Integrated simulations of premagnetized and self-magnetizing dynamic screw pinch-driven MagLIF

Magnetically driven implosions such as in magnetized liner inertial fusion (MagLIF) on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that dynamically redistribute liner mass during implosion, limiting fusion fuel compression and confinement, which ultimately degrades performance. Driving the implosion with an initially helical drive field that dynamically shifts the direction of the magnetic field surrounding the liner (i.e., a dynamic screw pinch, DSP) is a method proposed to mitigate MRTI in-flight and improve target performance. In DSPs, the axial drive magnetic field component implodes the liner and diffuses through the shocked, melted liner material into the fuel throughout the implosion. Liners can be designed to enable enough axial magnetic flux to diffuse through the liner material to effectively magnetize the fuel region without the need of an initial axial magnetic field (i.e., from external field coils). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing drive magnetic fields composed of axial and azimuthal components (a helical drive field). These simulated DSP-driven MagLIF targets demonstrate improved fuel conditions and thermonuclear yield compared to a traditional MagLIF target implosion. Synthetic x-ray radiography of the imploding liner material and x-ray emission images of the fuel region at the time of peak neutron yield rate indicate superior implosion morphology for DSP-MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale mol­ecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗
DOE OSTI2025

Femtosecond temperature measurements of laser-shocked copper deduced from the intensity of the x-ray thermal diffuse scattering

We present 50-fs, single-shot measurements of the x-ray thermal diffuse scattering (TDS) from copper foils that have been shocked via nanosecond laser ablation up to pressures above ∼135 GPa. We hence deduce the x-ray Debye–Waller factor, providing a temperature measurement. The targets were laser-shocked with the DiPOLE 100-X laser at the High Energy Density endstation of the European X-ray Free-Electron Laser. Single x-ray pulses, with a photon energy of 18 keV, were scattered from the samples and recorded on Varex detectors. Despite the targets being highly textured (as evinced by large variations in the elastic scattering) and with such texture changing upon compression, the absolute intensity of the azimuthally averaged inelastic TDS between the Bragg peaks is largely insensitive to these changes, and allowing for both Compton scattering and the low-level scattering from a sacrificial ablator layer provides a reliable measurement of $T/Θ^2_D$, where Θ D is the Debye temperature. We compare our results with the predictions of the SESAME 3336 and LEOS 290 equations of state for copper and find good agreement within experimental errors. We, thus, demonstrate that single-shot temperature measurements of dynamically compressed materials can be made via thermal diffuse scattering of XFEL radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Formation of a spherical plasma liner for plasma-jet-driven magneto-inertial fusion

Plasma-jet-driven magneto-inertial fusion is an alternative approach to controlled nuclear fusion, which aims to utilize a line-replaceable dense plasma liner as a repetitive spherical compression driver. In this experiment, first measurements of the formation of a spherical argon plasma liner formed from 36 discrete pulsed plasma jets are obtained on the Plasma Liner Experiment. Properties including liner uniformity and morphology, plasma density, temperature, and ram pressure are assessed as a function of time throughout the implosion process and indicate an apparent transition from initial kinetic inter-jet interpenetration to collisional regime near stagnation times, in accordance with theoretical expectation. A lack of primary shock structures between adjacent jets during flight implies that arbitrarily smooth liners may be formed by way of corresponding improvements in jet parameters and control. The measurements facilitate the benchmarking of computational models and understanding the scaling of plasma liners toward fusion-relevant energy density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗