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 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 OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2022

The ultrafast pixel array camera system and its applications in high energy density physics

Diagnostics in high energy density physics, shock physics, and related fields are primarily driven by a need to record rapidly time-evolving signals in single-shot events. These measurements are often limited by channel count and signal degradation issues on cable links between the detector and digitizer. Here we present the Ultrafast Pixel Array Camera (UPAC), a compact and flexible detector readout system with 32 waveform-recording channels at up to 10 Gsample/s and 1.8 GHz analog bandwidth. The compact footprint allows the UPAC to be directly embedded in the detector environment. A key enabling technology is the PSEC4A chip, an eight-channel switch-capacitor array sampling device with up to 1056 samples/channel. The UPAC system includes a high-density input connector that can plug directly into an application-specific detector board, programmable control, and serial readout, with less than 5 W of power consumption in full operation. We present the UPAC design and characterization, including a measured timing resolution of ~20 ps or better on acquisitions of sub-nanosecond pulses with minimal system calibrations. Example applications of the UPAC are also shown to demonstrate operation of a solid-state streak camera, an ultrafast imaging array, and a neutron time-of-flight spectrometer.

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

On the determination of shocked surface velocities

Laser Doppler velocimetry is commonly fielded in shock physics experiments to diagnose the times of arrival, magnitudes, and durations of shock impulses as they arrive at shock loaded surfaces. However, in the presence of ejecta, the resulting frequency noise often obscures the surface and subsequent impulses at later times. In this context, we evaluate a dominant wave analysis technique applied to velocimetry data from explosively shocked Sn subjected to a second shockwave, on which ejecta are emitted subsequent to both the first and subsequent shockwave impulses. The analysis reveals surface velocities, the times of arrival of subsequent impulses, and their durations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗