Search NASA⌕ Search

SEARCH · Search NASA

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 OSTI2023

Hybrid Particle Method for Computational Shock Physics

A long-standing area of research for Eulerian shock wave physics codes has been the treatment of strength and damage for materials. Here we present a method that will aid in the analysis of strength and failure in shock physics applications where excessive diffusion of critical variables can occur and control the solution outcome. Eulerian methods excel for large deformation simulations in general but are inaccurate in capturing structural behavior. Lagrangian methods provide better structural response, but finite element meshes can become tangled. Therefore, a technique for merging Lagrangian and Eulerian treatments of material response, within a single numerical framework, was implemented in the Multiple Component computational shock physics hydrocode. The capability is a Lagrangian/Eulerian Particle Method (LEPM) that uses particles to interface a Lagrangian treatment of material strength with a more traditional Eulerian treatment of the Equation of State (EOS). Lagrangian numerical methods avoid the advection diffusion found in Eulerian methods, which typically strongly affects strength constitutive law internal variables, such as equivalent plastic strain, porosity and/or damage. The Lagrangian capability enhances existing capabilities and permits accurate predictions of high rate, large deformation and/or shock of mechanical structures.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2023

Shock physics and shadowgraphic measurements of laser-produced cerium plasmas

Shadowgraphic measurements are combined with theory on gas-dynamics to investigate the shock physics associated with nanosecond laser ablation of cerium metal targets. Time-resolved shadowgraphic imaging is performed to measure the propagation and attenuation of the laser-induced shockwave through air and argon atmospheres at various background pressures, where stronger shockwaves characterized by higher propagation velocities are observed for higher ablation laser irradiances and lower pressures. The Rankine-Hugoniot relations are also employed to estimate the pressure, temperature, density, and flow velocity of the shock-heated gas located immediately behind the shock front, predicting larger pressure ratios and higher temperatures for stronger laser-induced shockwaves.

42 ENGINEERING↗
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 OSTI2023

Final Scientific Report on: Magnetized Shock Physics and Convergent Flows

The objective of this project was to study magnetized shocks in plasmas using the Omega Laser Facility to drive a solid foil into a magnetized gas, which would be preionized by the x-rays from the laser-solid interaction. The MIFEDS (magneto-inertial fusion electrical discharge system) provided the magnetic field. Magnetized shocks are of interest in magneto-inertial fusion, such as MagLIF and magnetized spherical implosions with both direct and indirect drive, space physics, heliophysics, and astrophysics. In magneto-inertial fusion the shocks of interest are collisional and perpendicular to the magnetic field lines, since the magnetic field needs to be compressed. In space, the shocks of interest are typically collisionless, due to the low plasma density, and both perpendicular and parallel to the magnetic field. The marginally collisionless regime between these is also of interest in both laboratory, fusion, and space plasmas, but has barely been investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Imaging velocity interferometer system for any reflector (VISAR) diagnostics for high energy density sciences

Two variants of optical imaging velocimetry, specifically the one-dimensional streaked line-imaging and the two-dimensional time-resolved area-imaging versions of the Velocity Interferometer System for Any Reflector (VISAR), have become important diagnostics in high energy density sciences, including inertial confinement fusion and dynamic compression of condensed matter. Here, we give a brief review of the historical development of these techniques, then describe the current implementations at major high energy density (HED) facilities worldwide, including the OMEGA Laser Facility and the National Ignition Facility. We illustrate the versatility and power of these techniques by reviewing diverse applications of imaging VISARs for gas-gun and laser-driven dynamic compression experiments for materials science, shock physics, condensed matter physics, chemical physics, plasma physics, planetary science and astronomy, as well as a broad range of HED experiments and laser-driven inertial confinement fusion research.

47 OTHER INSTRUMENTATION↗
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 OSTI2023

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2023

ExactPack: A python library of exact analytic solutions

Verification of multi-physics simulation software against problems with known analytic or semi-analytic solutions is an important aspect of research into a wide variety of fields involving the motion of fluids, shock physics and other dynamic material properties. Previous work comparing simulation results against analytic solutions has been ad-hoc, with developers frequently writing their own analytic solvers. This has resulted in a large amount of duplicated effort. The python library ExactPack has been developed as a collection of analytic and semi-analytic solvers to a variety of multi-physics problems, providing a consistent API to a set of well-tested solver implementations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗