Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Empirical correlations between the function of entropy ( Z S ) and net artificial viscous work in a shock physics hydrocode

Entropy is a state variable that may be obtained from any thermodynamically complete equation of state (EOS). However, hydrocode calculations that output the entropy often contain numerical errors; this is not because of the EOS, but rather the solution techniques that are used in hydrocodes (especially Eulerian) such as convection, remapping, and artificial viscosity. Here, in this work, empirical correlations are investigated to reduce the errors in entropy without altering the solution techniques for the conservation of mass, momentum, and energy. Specifically, these correlations are developed for the function of entropy Z S , and they depend upon the net artificial viscous work, as determined via Sandia National Laboratories’ shock physics hydrocode CTH. These results are a continuation of a prior effort to implement the entropy-based CREST reactive burn model in CTH, and they are presented here to stimulate further interest from the shock physics community. Future work is planned to study higher-dimensional shock waves, shock wave interactions, and possible ties between the empirical correlations and a physical law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

SIMD Programming for the SMASH Shock Physics Code

Many modern CPUs that are available to the NNSA as mission computing resources support vector instruction sets. Making good use of vector instructions, referred to as “vectorization”, is often critical to getting the best performance from these CPUs. While other codes choose to rely on compiler auto-vectorization, the SMASH shock physics code chooses to leverage APIs for explicit vectorization. These APIs are similar to directly calling the CPU vendor’s vector intrinsics, with the additional benefit of being vendor-agnostic. This document explains what the SIMD APIs are and how to use them in developing SMASH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

CRADA Abstract 2025-01U-JWS01

CRADA abstract for collaboration with Washington State University on shock physics research. Public abstract at time of execution is required per DOE O 483.1B.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Learning robust parameter inference and density reconstruction in flyer plate impact experiments

Estimating physical parameters or material properties from experimental observations is a common objective in many areas of physics and material science. In many experiments, especially in shock physics, radiography is the primary means of observing the system of interest. However, radiography does not provide direct access to key state variables, such as density, which prevents the application of traditional parameter estimation approaches. Here we focus on flyer plate impact experiments on porous materials, and resolving the underlying parameterized equation of state (EoS) and crush porosity model parameters given radiographic observation(s). We use machine learning as a tool to demonstrate with high confidence that using only high impact velocity data does not provide sufficient information to accurately infer both EoS and crush model parameters, even with fully resolved density fields or a dynamic sequence of images. We thus propose an observable data set consisting of low and high impact velocity experiments/simulations that capture different regimes of compaction and shock propagation, and proceed to introduce a generative machine learning approach which produces a posterior distribution of physical parameters directly from radiographs. We demonstrate the effectiveness of the approach in estimating parameters from simulated flyer plate impact experiments, and show that the obtained estimates of EoS and crush model parameters can then be used in hydrodynamic simulations to obtain accurate and physically admissible density reconstructions. Finally, we examine the robustness of the approach to model mismatches, and find that the learned approach can provide useful parameter estimates in the presence of out-of-distribution radiographic noise and previously unseen physics, thereby promoting a potential breakthrough in estimating material properties from experimental radiographic images.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

X-ray microscopy and talbot imaging with the matter in extreme conditions X-ray imager at LCLS

The last decade has shown the great potential that X-ray Free Electron Lasers (FEL) have to study High Energy Density (HED) physics. Experiments at FELs have made significant breakthroughs in Shock Physics and Dynamic Diffraction, Dense Plasma Physics and Warm Dense Matter Science, using techniques such as isochoric heating, inelastic scattering, small angle scattering and X-ray diffraction. In addition, and complementary to these techniques, the coherent properties of the FEL beam can be used to image HED samples with high fidelity. We present new imaging diagnostics and techniques developed at the Matter in Extreme Conditions (MEC) instrument at Linac Coherent Light Source (LCLS) over the last few years. We show results in Phase Contrast Imaging geometry, where the X-ray beam propagates from the target to a camera revealing its phase, as well as in Direct Imaging geometry, where a real image of the sample plane is produced in the camera with a spatial resolution down to 200 nm. Last, we show an implementation of the Talbot Imaging method allowing both X-ray phase and intensity measurements change introduced by a target with sub-micron resolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Rotation‐Controlled Diurnal Evolution of Uranus' Asymmetric Bow Shock at Equinox

Uranus possesses the most extreme magnetic and rotational geometry in the solar system, resulting in a uniquely dynamic and asymmetric interaction between its magnetosphere and the solar wind. Here we investigate the diurnal evolution of the Uranian bow shock (BS) at equinox using global multifluid magnetohydrodynamic simulations constrained by Voyager 2 observations. The BS structure and variability are quantified using the stand‐off distance, terminator distance, flaring parameter, and BS parameter, enabling a systematic assessment of its global asymmetry over one planetary rotation. Our results show that the BS exhibits pronounced diurnal asymmetry that is primarily controlled by planetary rotation. Even under steady upstream solar wind conditions, the BS undergoes periodic expansion and contraction, reflecting rotation‐driven reconfiguration of the planetary magnetospheric topology. This diurnal modulation represents a remarkably large asymmetry in planetary BS geometry within the solar system. These results establish planetary rotation as the intrinsic primary driver of Uranus' BS dynamics at equinox and provide essential constraints for the design and formulation of future space missions to the ice giants. These findings also offer a unique framework for understanding shock physics in the abundant population of ice‐giant exoplanets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
DOE OSTI2025

Modeling and Design of a Sub-Nyquist PDV Digitization System

Photon-Doppler Velocimetry (PDV) uses interferometry of a transmitted optical signal reflected from a fast-moving device under test (DUT) to generate a product signal at a doppler shifted frequency which is proportional to the velocity of the object being measured. The technique has applications in shock physics, where it is used to measure fast-moving objects over a short lifespan of travel, such as a bullet, shrapnel, or other shock accelerated object.

42 ENGINEERING↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2025

Evidence for an Early Formation of Serenitatis Basin at 4.25 Ga Shifts Lunar Chronology

Troctolite sample 76535, collected in Serenitatis basin during Apollo 17, formed at least 50 km deep, experienced maximum shock pressures of 6 GPa, and has a 40 Ar/ 39 Ar excavation age of 4.25 Ga. Previous work attributed 76535 to the South Pole-Aitken (SPA) basin, presumably dating the SPA-impact and constraining lunar bombardment history. Here we use the iSALE-2D shock-physics code and gravity inversion modeling to determine if instead the Serenitatis impact event excavated 76535. We find nearly 140,000 km 3 of material (∼2% of near-surface ejecta) matching the depth and pressure constraints of 76535 is displaced to the surface during crater collapse of a Serenitatis-like impact event. We conclude that the Serenitatis impact event possibly excavated 76535, redefining its age to 4.25 Ga, 300 My older than the consensus age based on Apollo 17 samples. This finding would provide an important anchor point where lunar chronology where bombardment flux is especially uncertain.

Bjonnes, Evan [Lawrence Livermore National Laborat↗
DOE OSTI2026

Validation of 266 nm UV Raman thermometry using the alpha-to-beta quartz transition

Ultraviolet Raman scattering is a promising diagnostic for temperature measurements relevant to dynamic compression experiments, where improved signal strength and reduced background are critical. In this work, we experimentally validate Raman thermometry at 266 nm using crystalline quartz heated through the α to β phase transition in a tube furnace. Quartz provides a useful benchmark because of its well-characterized transition temperature near 846 K at ambient pressure. Temperatures inferred from the Stokes to anti-Stokes intensity ratio are compared with thermocouple measurements and with the observed phase transition. Single-shot Raman temperatures track the thermocouple closely, with an empirical single-shot scatter of approximately ±15 K under the present conditions. Temporal averaging reduces shot-to-shot scatter and improves agreement between the Raman temperatures and thermocouple readings to within approximately 8 K, while the transition is identified within 5 K of the accepted value. Measurements are performed using the fourth harmonic of a nanosecond pulsed Nd:YAG laser at 266 nm. The shorter excitation wavelength provides enhanced Raman scattering efficiency and reduced visible background relative to conventional visible excitation. Comparison with 532 nm Raman measurements acquired in the same experimental geometry demonstrates substantially improved signal-to-noise ratio for the UV configuration. These results establish 266 nm Raman thermometry as a reliable and accurate temperature diagnostic under controlled heating conditions and provide a quantitative benchmark for future application in dynamic compression and shock physics experiments.

Chemistry - Chemical explosives↗
DOE OSTI2024

Reconstructing Richtmyer–Meshkov instabilities from noisy radiographs using low dimensional features and attention-based neural networks

We develop an ML-based approach for density reconstruction based on transformer neural networks. This approach is demonstrated in the setting of ICF-like double shell hydrodynamic simulations wherein the parameters related to material properties and initial conditions are varied. The new method can robustly recover the complex topologies given by the Richtmyer-Meshkoff instability (RMI) from a sequence of hydrodynamic features derived from radiographic images corrupted with blur, scatter, and noise. A noise model is developed to characterize errors in extracting features from synthetic radiographs of the simulated density field. The key component of the network is a transformer encoder that acts on a sequence of features extracted from noisy radiographs. This encoder includes numerous self-attention layers that act to learn temporal dependencies in the input sequences and increase the expressiveness of the model. This approach is shown to exhibit an excellent ability to accurately recover the RMI growth rates, despite the gas-metal interface being greatly obscured by radiographic noise. Our approach can be applied in a broad array of fields involving shock physics and material science.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Single-shot imaging with randomized structured illumination at a free electron laser

Stroboscopic nanoscale imaging with free electron laser light is revolutionizing our understanding of fast dynamics in heterogeneous systems. The short wavelength of X-ray and extreme ultraviolet radiation makes it possible to achieve nanoscale resolution, while resonance with atomic transitions gives access to electronic and magnetic degrees of freedom. Here, we report on our implementation of a recently developed imaging method, randomized probe imaging, at a free electron laser. The advantage of randomized probe imaging over existing methods is its compatibility with extended and strongly scattering samples. Our implementation delivers robust single-shot reconstructions at up to a full-pitch resolution of 400 nm over a field of view with a 40 µm diameter. We also demonstrate single-shot imaging of magnetic domain structures using circular dichroism at resonance, paving the way to future time-resolved studies of magnetic dynamics, shock physics, and the dynamics of collective electronic phases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Four Channel Time Multiplexed Photonic Doppler Velocimetry using an Optical Switch

Photonic Doppler Velocimetry (PDV) is a diagnostic commonly used in shock physics and dynamic compression experiments to reliably get velocity information from experiments. In PDV systems, a common method of reducing experimental cost is to use time and frequency multiplexing to increase the number of PDV probes. With time multiplexing, interference between probes is a frequent problem. In this report, we look at using a high-speed optical switch to reduce this interference, including measuring the amount of interference generated to determine if it has the potential to affect experiments and integrating a time multiplexing system into an experiment. We find that, when applied to PDV systems, there is approximately (-23.4 ± 0.9) dB of interference measured in the short time Fourier transform between switch inputs. When an optical switch based time multiplexing system was integrated into a dynamic compression experiment, the system was able to successfully combine the signals from four different PDV probes onto a single optical cable without unacceptable levels of interference in the spectrogram. An optical switch based time multiplexing system appears to be a promising method for reducing the cost of fielding larger numbers of PDV probes in an experiment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

The right conditions for high-precision dynamic temperature and heat capacity measurement via pyrometry and conductivity

The pursuit of accurate bulk temperature T under extreme conditions has been a long-standing goal of the high pressure science community, complicated by a lack of data to inform models. To reach these extremely high-pressure, high-temperature (high P − T) conditions, a combination of dynamic and heated static experiments (e.g., diamond or gem anvil cel experiments) are used. For example, in a diamond anvil cell (DAC) experiment, a sample placed in the DAC is first pressurized. Following pressurization, the sample T is increased either by heating the entire DAC (usually using resistive heating, and limited to ∼1000K) or by applying intense laser power to the sample surfaces. In a dynamic experiment, the process of pressurizing the sample also heats it. In the case of shock physics experiments, such heating is substantial, easily reaching thousands of Kelvin; in our work we have seen T ∼17000K. Most methods of measuring temperature at ambient are not compatible with experiments under these high-pressure, high-temperature conditions: thermocouples break, melt, or have conductivity properties that differ from ambient where they are calibrated; thermometers would melt; both are too slow. As a result most methods are based on non-contact techniques such as x-ray diffraction broadening, neutron scattering, or optical methods. Of these, optical methods using the visible and near-infrared region of the spectrum are the most commonly used as the sources and detectors are readily available. In the case of optical methods the optical depth, and therefore the measurement location, is limited to the surface. When a window or anvil material is used, heat flows from the sample into the window/anvil. Likewise, if the sample undergoes a change in thermodynamic state, such as expansion upon release, different T may be expected. As a result, the surface or apparent temperature T app measurement will differ from the bulk or interior temperature that is desired. This surface measurement must be related to the bulk measurement using thermal transport models and material models. While it is tempting to conclude that one should just use x-ray methods that directly probe the interior, even these methods have been shown to depend on thermal transport and material models. Regardless of the method used to create the high P − T condition, therefore, we must understand the role of thermal transport and material models upon our interpretation of the T measurement, as well as the errors and uncertainties associated with the choice of models used in the analysis. This is a substantial area of research and this paper is by no means a complete survey of the relevant sources of uncertainty. For example, we have yet to begin to address alternate transport models in a detailed manner (e.g., Tan-Ahrens), or the many models that use additional layers to approximate melting, turbulence, or epitaxial phenomena). Likewise, we have not explored the impact upon uncertainty of thermal models that use temperature-dependent thermal transport coefficients, or the wide range of material models that can be applied. Instead, this paper focuses on using one simple model, the Urtiew-Grover model, to understand the sources of error in T measurement so that we may identify how best to focus future research efforts to return the best improvements and avoid working on over-optimizing a single type of measurement. To this end, we work through some of the best and worst case scenarios for T measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Optical Test Bench (SYSM 5620 Final Report)

The National Ignition Facility (NIF) is a research facility that focuses on High Energy Density Science and Inertial Confinement Fusion. In the NIF, 192 laser beams are directed into a miniscule Target Capsule to generate an implosion – this is called a NIF Shot. Many NIF Users conduct experiments using different Targets and laser inputs to generate reactions that are of interest to the scientific community. Those reactions are monitored using various Target Diagnostic instruments, including Shot Time Photography – which provides crucial images of the Target as the implosion occurs. A subset of these NIF Users is less interested in the reaction itself and more interested in subjecting materials or devices to the intense environment that these reactions can create. Using the Target as a source of x-rays, neutrons, or blast loading for particular materials has been significant to the shock physics community and to those who work in extreme design spaces. However, the Material and Radiation Effects team – that performs these experiments on behalf of the NIF Users – has not been able to make good use of the existing Shot Time Photography Diagnostics. The STP cameras are focused on the Targets rather than the MRE Customers’ materials, the current system is not oriented to have a good view of the MRE Customers’ materials, and the cameras are limited in their capabilities (black and white, poor resolution, etc.). As MRE experiments increase in quantity and importance to the scientific community, the desire for visual monitoring of the experiments to substantiate models and inform other data is also rising.

42 ENGINEERING↗
DOE OSTI2026

Variation in surface properties as error sources in dynamic temperature measurements

Most temperature measurement methods rely upon the assumption of a single uniform source as part of the analysis in inferring an observed temperature T from the surface of a material. This is true for the radiometric and spectroscopic methods we rely upon in shock physics. However, there are obvious cases where this is a poor assumption that can cause large errors in analysis, such as the measurement of emission from a surface composed of more than one material, such as a rusted surface of Fe and Fe 2 O 3 , or a multi-material surface that could be produced through alloying methods (e.g., mokume gane) or 3D printing. Here we perform an initial examination of the effects of varied surface properties upon the analysis using idealized test data to constrain when we may have better confidence in the resulting calculations. We also examine the structures of the fits to test whether, in principle, comparing the fit with ideal gray bodies would reveal such a structure and enable the researcher to identify when a two-temperature fitting model should be used instead of a single- T model or a wavelength-dependent emissivity model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multi-frame x-ray radiography and image tracking for quantification of expansion in laser-driven tin ejecta microjets

One regime of experimental particle-laden flow study involves ejecta microjets—often defined as a stream of micrometer-scale particles generated through shock interaction with a non-uniform surface and generally travel above 1 km/s. In order to capture the change in characteristics as a function of propagation time, we apply a multi-frame x-ray radiography platform to observe and track the jet transport dynamics. A synchrotron x-ray source allows us to perform quantitative analyses and comparisons between the eight images captured by the imaging system. Observation of a single jet through time allows the use of a cross correlation algorithm to independently track various regions within the jet and quantify the jet expansion over time using normalized area and normalized areal density values. Through a comparison with the calculated values of ballistic transport, these findings show less expansion than expected for ballistically transporting particles. In conclusion, this work combines multi-frame synchrotron radiography with image tracking to establish a foundation for future studies on jet transport and particle interaction dynamics.

Sun, Yuchen [Lawrence Livermore National Laborator↗