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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE 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↗

NASA NTRS1995

The termination shock: Physical processes

A brief overview of the physics thought to determine the properties of the heliospheric termination shock is given and the central role of hot interstellar pick-up ions and cosmic rays is emphasized.

Zank, G. P.↗

NASA NTRS2012

SN 1987 A: A Unique Laboratory for Shock Physics

Supernova 1987 A has given us an unprecedented view of the evolution of the explosion debris and its interaction with circumstellar matter. The outer supernova debris, now expanding with velocities approx.8000 km/s, encountered the relatively dense circumstellar ring formed by presupernova mass loss in the early 1990s. The shock interaction is manifested by UV-optical "hotspots", an expanding X-ray ring, an expanding ring of knotty non-thermal radio emission, and a ring of thermal IR emission from silicate dust Recent ultraviolet observations of the emissions from the reverse shock and the ring with the HST/COS reveal new details about the shock interaction. Lyman alpha emission from the reverse shock is much stronger than H alpha and they have different emission morphologies, pointing to different emission mechanisms. The reverse shock was detected for the first time in C IV 1550. The N V to C IV brightness ratio indicates the N/C abundance ratio in the expanding debris is about 100X solar, about 3X N/C in the inner ring.

Sonneborn, George↗

NASA NTRS2012

SN 1987A: A Unique Laboratory for Shock Physics

Supernova 1987 A is the brightest and nearest supernova observed since Kepler's SN1604, and is the only one close enough to resolve and directly observe the temporal growth of the ejecta. Over the past 25 years, intensive observations across the electromagnetic spectrum with observatories on the ground (Australia Telescope Compact Array, Gemini-S, Magellan, VLT) and in space (IUE, KAO, CGRO, Hubble, Chandra, Spitzer, Herschel) have given us an unprecedented view of the evolution of the debris of the supernova and of its shock interaction with circumstellar matter. The outer supernova debris, now expanding with velocities -8000 km/s, encountered the relatively dense circumstellar ring formed by presupernova mass loss starting in 1994. The resulting shock interaction has been manifested by: rapidly brightening UV-optical "hotspots", an expanding X-ray ring. an expanding ring of knotty non-thermal radio emission, and a ring of thermal IR emission from silicate dust. The recent evolution of these emissions reveal new details about the shock interaction, circumstellar material, and the star that exploded. Certain critical problems about SN 1987 A, such as the still undiscovered compact object formed in the explosion and the structure of the central debris, require the capabilities of JWST.

Sonneborn, George↗

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↗

NASA NTRS1979

Molecule formation and infrared emission in fast interstellar shocks. I Physical processes

The paper analyzes the structure of fast shocks incident upon interstellar gas of ambient density from 10 to the 7th per cu cm, while focusing on the problems of formation and destruction of molecules and infrared emission in the cooling, neutral post shock gas. It is noted that such fast shocks initially dissociate almost all preexisting molecules. Discussion covers the physical processes which determine the post shock structure between 10 to the 4 and 10 to the 2 K. It is shown that the chemistry of important molecular coolants H2, CO, OH, and H2O, as well as HD and CH, is reduced to a relatively small set of gas phase and grain surface reactions. Also, the chemistry follows the slow conversion of atomic hydrogen into H2, which primarily occurs on grain surfaces. The dependence of this H2 formation rate on grain and gas temperatures is examined and the survival of grains behind fast shocks is discussed. Post shock heating and cooling rates are calculated and an appropriate, analytic, universal cooling function is developed for molecules other than hydrogen which includes opacities from both the dust and the lines.

Hollenbach, D.↗

NASA NTRS1991

The plasma physics of shock acceleration

The history and theory of shock acceleration is reviewed, paying particular attention to theories of parallel shocks which include the backreaction of accelerated particles on the shock structure. The work that computer simulations, both plasma and Monte Carlo, are playing in revealing how thermal ions interact with shocks and how particle acceleration appears to be an inevitable and necessary part of the basic plasma physics that governs collisionless shocks is discussed. Some of the outstanding problems that still confront theorists and observers in this field are described.

Jones, Frank C.↗

NASA NTRS2012

Atomic Physics of Shocked Plasma in Winds of Massive Stars

High resolution diffraction grating spectra of X-ray emission from massive stars obtained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key wind parameters. Modeling of resolved emission line velocity profiles allows us to derive independent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of Xray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray emission begins moderately close to the stellar surface and extends throughout the wind. Some outstanding questions remain, including the possibility of large optical depths in resonance lines, which is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place constraints on shock size, density, and velocity structure

Physics↗

NASA NTRS2013

Understanding the Flow Physics of Shock Boundary-Layer Interactions Using CFD and Numerical Analyses

Computational fluid dynamic (CFD) analyses of the University of Michigan (UM) Shock/Boundary-Layer Interaction (SBLI) experiments were performed as an extension of the CFD SBLI Workshop held at the 48th AIAA Aerospace Sciences Meeting in 2010. In particular, the UM Mach 2.75 Glass Tunnel with a semi-spanning 7.75deg wedge was analyzed in attempts to explore key physics pertinent to SBLI's, including thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses were 3D CFD simulations using the OVERFLOW flow solver, with additional quasi-1D simulations performed with an in house MATLAB code interfacing with the NIST REFPROP code to explore perfect verses non-ideal air. A fundamental exploration pertaining to the effects of particle image velocimetry (PIV) on post-processing data is also shown. Results from the CFD simulations showed an improvement in agreement with experimental data with key contributions including adding a laminar zone upstream of the wedge and the necessity of mimicking PIV particle lag for comparisons. Results from the quasi-1D simulation showed that there was little difference between perfect and non-ideal air for the configuration presented.

Friedlander, David J.↗

NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

NASA NTRS2013

Understanding the Flow Physics of Shock Boundary Layer Interactions Using CFD and Numerical Analyses

No abstract available

Friedlander, David↗

NASA NTRS1987

An analysis of artificial viscosity effects on reacting flows using a spectral multi-domain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, M. G.↗

NASA NTRS1989

Analysis of artificial viscosity effects on reacting flows using a spectral multidomain technique

Macaraeg, Michele G.↗

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↗

NASA NTRS2013

New Evidence for Efficient Collisionless Heating of Electrons at the Reverse Shock of a Young Supernova Remnant

Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect K beta (3p yields 1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe K alpha (2p yields 1s) emission dominated by a relatively highly ionized component. Comparisons with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe K alpha morphology from the Chandra observations. Since strong Fe K beta fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.

Fe K diagnostics↗

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↗