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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

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 molecular 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↗

DOE OSTI2025

A High-Order 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 due to 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 reactive Euler equations encountered in high-speed combustion. 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 DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗

DOE OSTI2024

Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver

A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

A magnetic analog of pressure–strain interaction

We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.

Barbhuiya, M. Hasan (ORCID:0000000163301650)↗

DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

High-Fidelity, Low-Dissipation/Symmetry-Preserving Numerical Scheme for Solving the Euler Equations with Unstructured, Metric-Based Mesh Adaptation

This work presents an overview of a high-fidelity compressible Euler solver that utilizes the continuous Galerkin (CG) method with added artificial numerical diffusion for stabilization to solve a variety of unsteady and steady benchmark inviscid flow problems. This work shows that discretizing the Euler equations with this CG approach and first order basis functions produces a cost-effective stencil as well as simple well-posed boundary conditions. We show through convergence testing with manufactured solutions that the reduced stencil of CG, combined with the low amount of artificial diffusion required when using the stabilization method outlined in this work, leads to stable and highly accurate results for a variety of unsteady and steady applications. When combined with the adaptive mesh refinement approach used for many of the cases in this work, our results show that the flow solver achieves even more accurate results. A variety of inviscid flow cases are presented in this work, including transient 2D cases with complex shock structures and several steady 3D airfoils sections with a constant span.

Doetsch, Kevin [ORNL] (ORCID:0000000267051705)↗