Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shocks”

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

Sublimation and oxidation measurements of graphite and carbon black at high temperatures in a shock tube using absorption imaging and thermal emission

Surface mass loss rates due to sublimation and oxidation at temperatures of 3000–7000 K have been measured in a shock tube for graphite and carbon black (CB) particles. Diagnostics are presented for measuring surface mass loss rates by diffuse backlit illumination extinction imaging and thermal emission. Here, the surface mass loss rate is found by regression fitting extinction and emission signals with an independent spherical primary particle assumption. Measured graphite sublimation and oxidation rates are reported to be an order of magnitude greater than CB sublimation and oxidation rates. It is speculated that the difference between CB and graphite surface mass loss rates is largely due to the primary particle assumption of the presented technique which misrepresents the effective surface area of an aggregate particle where primary particles overlap and shield inner particles. Measured sublimation rates are compared to sublimation models in the literature, and it is seen graphite shows fair agreement with the models while CB underestimates, likely a result of the particle shielding affect not being considered in the sublimation model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Styrene Thermal Decomposition and Its Reaction with Acetylene under Shock Tube Pyrolysis Conditions: an Experimental and Kinetic Modeling Study

Styrene is an important compound for polymer production and a key intermediate in gas-phase kinetics of polycyclic aromatic hydrocarbons (PAHs). For the first time, the pyrolysis of styrene with and without the presence of acetylene is investigated in a single-pulse shock tube coupled to gas chromatography and mass spectrometry. For each reaction system, quantitative speciation profiles are probed within the temperature range of 1100-1730 K, nominal pressure of 20 bar, and reaction duration of similar to 4 ms. A kinetic model is built to simulate the results. The model explains how styrene is consumed under high-pressure pyrolytic conditions, how the secondary chemistry of intermediate products affect subsequent PAH formation, and how acetylene addition alters the reaction pathways. Throughout the temperature range, styrene breakdown is dominated by the bimolecular interaction between styrene and hydrogen atom, which produces benzene+vinyl or phenyl and ethylene through the stabilization of 2-phenylethyl and its subsequent dissociation. As a result, large amounts of phenyl accumulate, which react with styrene to form C14H12 species while simultaneously releasing H atoms through addition/elimination reactions. The reactivity of fuel consumption is preserved by the regeneration of H atoms as chain carriers. Several C14H10 compounds are formed as a result of the following breakdown of the C14H12 isomers, particularly stilbene, 1,1-diphenyl ethylene, and 9-methyl-9H-fluorene. The presence of acetylene as a co-reactant with styrene allows the Hydrogen-Abstraction-Acetylene-Addition (HACA) pathway to proceed from phenyl radical to enhance the production of phenylacetylene at very low temperatures and acenaphthylene. This hinders the formation of C14H12 isomers, exclusive products from pure styrene dissociation by competing with the styrene+phenyl routes.

Kinetic Modeling↗
DOE OSTI2025

Laser-driven flash x-ray radiography of a shocked metallic foil

Characterizing hydrodynamic instability evolution in millimeter-scale, high-Z foils is crucial for understanding complex phenomena in high-energy-density physics. Here, we demonstrate a proof-of-concept, laser-driven flash x-ray radiography platform tailored for two-dimensional linear density mapping in shocked high-Z foils. Using chromium (Cr) foils with internal shockwaves (∼100 μm width), our platform achieves a spatial resolution of 59.8 ± 1.4 μm by employing a broadband x-ray source extending into the hundreds of keV range. The setup combines a compound parabolic concentrator cone with a tantalum wire target, a magnetic field to deflect residual transmitted electrons, and a copper casing to shield the sides and rear of the image plate pack. By varying the delay of the short-pulse beam driving the flash x-ray source, we resolve shockwave dynamics, specifically the velocity, position, width, and density profile, within the Cr foil. Reported experimental results are consistent with the corresponding hydrodynamics and radiation transport simulations, which accurately reproduce the measured electron and x-ray source terms. These developments enable the conversion of shockwave radiographs into two-dimensional density maps, enhancing interpretability for hydrodynamic instability evolution applications and validating the simulation approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2025

Laser Measurements and Modeling of Shock Tunnel Freestream Velocity and Multispecies Thermal Nonequilibrium

Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.

Aerodynamics↗
DOE OSTI2024

A Comprehensive Review on Finite Element Analysis of Laser Shock Peening

Laser shock peening (LSP) is a formidable cold working surface treatment that provides high-energy precision to enhance the mechanical properties of materials. This paper delves into the intricacies of the LSP process, offering insights into its methodology and the simulation thereof through the finite element method. This review critically examines various points, such as laser energy, overlapping of shots, effect of LSP on residual stress, effect of LSP on grain refinement, and algorithms for simulation extrapolated from finite element analyses conducted by researchers, shedding light on the nuanced considerations integral to this technique. As the significance of LSP continues to grow, the collective findings underscore its potential as a transformative technology for fortifying materials against mechanical stress and improving their overall performance and longevity. The discourse encapsulates the evolving landscape of the LSP, emphasizing the pivotal role played by finite element analysis in advancing our understanding and application of this innovative surface treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Improving Chirped Fiber Bragg Grating Resolution for Position-Sensitive Sensors in Shock- and Detonation-Driven Experiments

Chirped fiber Bragg gratings (CFBGs) are robust diagnostic sensors that are widely used to track detonation-driven and shock wave propagation. CFBGs are inscribed with a linearly chirped periodic index of refraction changes that alter the Bragg wavelength along the length of the probe. The light return of each individual Bragg element is captured by a detector at a unique time to map the full reflected spectrum. The CFBG spectrum is measured with a dispersive Fourier transform of the reflected light that temporally stretches the spectrum to increase spatial resolution and make a one-to-one map of the wavelength on a time axis. Here, we propose an improvement of CFBG temporal resolution by incorporating two co-linear laser pulses with orthogonal polarization states and a 5 ns time offset. The two separate signals were split and tracked by two separate detectors. An oscilloscope captured good separation in the signals, and two separate spectrograms were generated and interleaved in the post-processing of the data. This novel technique doubled the CFBG temporal resolution and led to a doubled location resolution. As a proof-of-concept of this technique, the resolution improvement was compared between standard CFBG measurements and the two polarization states method on a position-sensitive CFBG sensor. CFBG resolution doubling will advance sensor capabilities and will have a direct impact on improving capture and analysis in dynamic, high-explosive experiments.

42 ENGINEERING↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2025

Multi-messenger dynamic imaging of laser-driven shocks in water using a plasma wakefield accelerator

Understanding dense matter hydrodynamics is critical for predicting plasma behavior in environments relevant to laser-driven inertial confinement fusion. Traditional diagnostic sources face limitations in brightness, spatiotemporal resolution, and in their ability to detect relevant electromagnetic fields. In this work, we present a dual-probe, multi-messenger laser wakefield accelerator platform combining ultrafast X-rays and relativistic electron beams at 1 Hz, to interrogate a free-flowing water target in vacuum, heated by an intense 200 ps laser pulse. This scheme enables high-repetition-rate tracking the evolution of the interaction using both particle types. Betatron X-rays reveal a cylindrically symmetric shock compression morphology assisted by low-density vapor, resembling foam-layer-assisted fusion targets. The synchronized electron beam detects time-evolving electromagnetic fields, uncovering charge separation and ion species differentiation during plasma expansion – phenomena not captured by photons or hydrodynamic simulations. We show that combining both probes provides complementary insights spanning kinetic to hydrodynamic regimes, highlighting the need for hybrid physics models to accurately predict fusion-relevant plasma behavior.

imaging techniques↗
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 OSTI2025

Global heliospheric termination shock strength in the solar–interstellar interaction

A heliospheric termination shock (HTS) surrounds our Solar System at approximately 100 astronomical units from the Sun, where the expanding solar wind (SW) is compressed and heated before encountering the interstellar medium. HTS-accelerated particles govern the pressure balance with the interstellar medium, but little is known about the global properties of the HTS beyond in situ measurements from Voyager in only two directions of the sky. Here we fill this gap by extracting the HTS strength using particle-in-cell, test particle and magnetohydrodynamic simulations, constrained by Interstellar Boundary Explorer observations of energetic neutral atoms produced from HTS-accelerated particles. Our results reveal there is a higher compression near the poles during solar minimum compared with solar maximum due to the higher Mach number flow. North–south asymmetries arise from the disparate evolution of the polar coronal holes, while minimum compression near the flanks is probably due to SW slowing from mass loading over a greater distance to the HTS. The results imply a strong connection between the HTS strength and the SW and interstellar medium dynamics.

Astrophysical plasmas↗
DOE OSTI2024

Quantifying motion blur by imaging shock front propagation with broadband and narrowband X-ray sources

Time-integrated radiography using MeV Bremsstrahlung X-ray sources is the norm for imaging during system-level testing of components and structures under dynamic condition. One source of error in the analysis of the time-integrated radiography data sets stems from motion blur which smears out sharp interfaces to a greater degree with longer exposure times, which become necessary to provide sufficient signal-to-noise with low X-ray penetration of objects of interest. To quantify motion blur, a 1D shock wave through PMMA was investigated experimentally at The Dynamic Compression Sector at The Advanced Photon Source (DCS@APS) with tapered broadband and 25.46 ± 1.06 keV narrowband X-rays. Four cameras with different exposure times were used for each experiment to compare the effect that exposure time has on motion blur. In addition, our methodology to accurately simulate motion blur in terms of transmission and shape is presented and compared to our experimental results and quantified. There is a high level of agreement between the experimental and simulation results across the range of data sets investigated in this study with a percent difference range of 0.29–1.31% for the four shots. The methodology of this work serves as a steppingstone towards a physically validated model that could be used in conjunction with experimental results to deconvolve physical parameters, densities, and interfaces of interest in a way that would not be possible with experimental results alone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Data-driven picosecond X-ray imaging for quantitative plasma-induced shock characterization

Imaging dynamic events, especially shockwave behavior, is key to advancing high-energy-density (HED) research. Recent advances in fourth- and fifth-generation X-ray light sources allow for high-resolution imaging of fast phenomena, but limited beam time necessitates maximizing data acquisition. We present a benchtop-scale pulsed plasma device submerged in liquid heptane, capable of generating dynamic events at rates exceeding 10 Hz, supporting the field’s data-driven goals by producing large, high-quality imaging datasets. Using X-ray phase contrast imaging (XPCI) at the Advanced Photon Source, we imaged weak shockwaves (Mach ~ 1.2) in heptane interacting with plasma-induced cavitation bubbles, causing deviation from Rankine-Hugoniot behavior; to our knowledge, this represents the first direct imaging of such interaction. Our quantitative analysis offers insight into weak shock phenomena and energy-focusing applications in pulsed plasmas. These results highlight the potential for large datasets to advance dynamic HED research at current light source facilities, and have implications for fields such as inertial confinement fusion, plasma-enhanced chemical processing, and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Abrupt reduction in shipping emission as an inadvertent geoengineering termination shock produces substantial radiative warming

Abstract Human activities affect the Earth’s climate through modifying the composition of the atmosphere, which then creates radiative forcing that drives climate change. The warming effect of anthropogenic greenhouse gases has been partially balanced by the cooling effect of anthropogenic aerosols. In 2020, fuel regulations abruptly reduced the emission of sulfur dioxide from international shipping by about 80% and created an inadvertent geoengineering termination shock with global impact. Here we estimate the regulation leads to a radiative forcing of $$+0.2\pm 0.11$$ + 0.2 ± 0.11 Wm −2 averaged over the global ocean. The amount of radiative forcing could lead to a doubling (or more) of the warming rate in the 2020 s compared with the rate since 1980 with strong spatiotemporal heterogeneity. The warming effect is consistent with the recent observed strong warming in 2023 and expected to make the 2020 s anomalously warm. The forcing is equivalent in magnitude to 80% of the measured increase in planetary heat uptake since 2020. The radiative forcing also has strong hemispheric contrast, which has important implications for precipitation pattern changes. Our result suggests marine cloud brightening may be a viable geoengineering method in temporarily cooling the climate that has its unique challenges due to inherent spatiotemporal heterogeneity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Strength and thermal shock resistance of porous reaction-bonded silicon nitride by direct nitriding of binder jet additively manufactured silicon

The direct nitridation of binder jet additively manufactured silicon to produce porous silicon nitride was explored. The Taguchi Design of Experiments method was used to systematically study the effect of time, temperature, and nitrogen flow rate on the direct nitriding of binder jet additively manufactured silicon to produce porous reaction-bonded silicon nitride (RBSN). Highly porous (61%) samples with high conversion (99%) to silicon nitride were achieved using a thermal profile that utilized two isothermal holds. Further, it was shown that these samples exhibited thermal shock resistance and a low coefficient of thermal expansion. These results indicate that RBSN is readily achievable through additive manufacturing.

Ceramics↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗