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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Supernova remnants and the physics of strong shock waves

This paper reports on a Workshop on Supernova Remnants and the Physics of Strong Shock Waves hosted by North Carolina State University at Raleigh, North Carolina, September 16-18, 1993. The workshop brought together observers, shock theorists, cosmic-ray specialists, and simulators to address the role supernova remnants can play in furthering our understanding of the complex plasma physics associated with collisionless shocks and particle acceleration. Over fifty scientists presented papers on various aspects of supernova remnants. In lieu of a proceedings volume, we present here a synopsis of the workshop, in the form of brief summaries of each workshop session.

Ellison, Donald C.↗
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↗
NASA NTRS1995

Electron distributions in quasi-perpendicular shocks: Simulations

One of the outstanding problems of collisionless shock physics is the microphysical description of electron thermalization and associated wave processes. Observations at the Earth's bow shock have shown various features in the electron velocity distribution such as a loss-cone, temperature anisotropies, and beams before a thermalized, flat-topped distribution is reached downstream. It is generally believed that in addition to the macroscopic fields, which determine the overall mapping of the distribution function, electromagnetic whistler waves and lower hybrid waves are responsible for the pitch-angle scattering and thermalization of the electrons. However, a genuine physical explanation requires analysis of the evolution of the electron distribution in the self-consistent wave fields. Particle simulations are ideally suited to do this, and include processes in which both electrons and ions are involved, such as cross-field streaming instabilities. We show results from l-D and 2-D simulations using an implicit code, which allows time steps large compared to the inverse electron plasma frequency. It is shown that in order to capture the correct propagation characteristics of the waves, a large mass ratio between the protons and electrons is generally necessary. Detailed diagnostics of the waves as well as the evolution of the electron velocity distribution will be presented and will be discussed in the frame work of satellite observations.

Pantellini, F.↗
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↗
NASA NTRS2022

Energy Partition at Collisionless Supercritical Quasiperpendicular Shocks

Collisionless shocks in astrophysical plasmas are important thermalizers, converting some of the incident flow energy into thermal energy, and non-thermalizers, partitioning that energy in unequal ways to different particle species, sub-populations thereof, and field components. This partition problem, or equivalently the shock equation of state, lies at the heart of shock physics. Here we employ systematically a framework to capture all the incident and downstream energy fluxes at two example traversals of the Earth’s bow shock by the Magnetospheric Multiscale Mission. Here and traditionally such data has to be augmented by information from other spacecraft, e.g., to provide more accurate measurements of the cold solar wind beam. With some care and fortuitous choices, the energy fluxes are constant, including instantaneous measurements through the shock layer. The dominant incident proton ram energy is converted primarily into downstream proton enthalpy flux, the majority of which is actually carried by a small fraction of suprathermal protons. Fluctuations include both real and instrumental effects. Separating these, resolving the solar wind beam, and other considerations point the way to a dedicated mission to solve this energy partition problem across a full range of plasma and shock conditions.

Steven J Schwartz↗
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↗
NASA NTRS2013

Launch Vehicle Debris Models and Crew Vehicle Ascent Abort Risk

For manned space launch systems, a reliable abort system is required to reduce the risks associated with a launch vehicle failure during ascent. Understanding the risks associated with failure environments can be achieved through the use of physics-based models of these environments. Debris fields due to destruction of the launch vehicle is one such environment. To better analyze the risk posed by debris, a physics-based model for generating launch vehicle debris catalogs has been developed. The model predicts the mass distribution of the debris field based on formulae developed from analysis of explosions. Imparted velocity distributions are computed using a shock-physics code to model the explosions within the launch vehicle. A comparison of the debris catalog with an existing catalog for the Shuttle external tank show good comparison in the debris characteristics and the predicted debris strike probability. The model is used to analyze the effects of number of debris pieces and velocity distributions on the strike probability and risk.

Gee, Ken↗
DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗
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↗
NASA NTRS2005

Fractionation and Accretion of Meteorite Parent Bodies

Senior Scientist Stuart J. Weidenschilling presents his final administrative report for the research program on which he was the Principal Investigator. The research program resulted in the following publications: 1) Particle-gas dynamics and primary accretion. J. N. Cuzzi and S. J . Weidenschilling. To appear in Meteorites and the Early Solar System 11 (D. Lauretta et a]., Eds.), Univ. Arizona Press. 2005; 2) Timescales of the solar protoplanetary disk. S. Russell, L. Hartmann, J . N. Cuzzi, A. Krot, M. Gounelle and S. J. Weidenschilling. To appear in Meteorites and the Early Solar System II (D. Lauretta et al., Eds.), Univ. Arizona Press, 2005; 3) Nebula evolution of thermally processed solids: Reconciling astrophysical models and chondritic meteorites. J. N. Cuzzi, F. J. Ciesla, M. I. Petaev, A. N. Krot, E. R. D. Scott and S . J. Weidenschilling. To appear in Chondrites and the Protoplanetary Disk (A. Krot et a]., Eds.), ASP Conference Series, 2005; 4) Possible chondrule formation in planetesimal bow shocks: Physical processes in the near vicinity of the planetesimal. L. L. Hood, F. J. Ciesla and S. J. Weidenschilling. To appear in Chondrites and the Protoplanetary Disk (A. Krot et al., Eds.), ASP Conference Series, 2005; 5) From icy grains to comets. In Comets II (M. Festou et al., Eds.), Univ. Arizona Press, pp. 97- 104, 2005; 6) Evaluating planetesimal bow shocks as sites for chondrule formation. F. J . Ciesla, L. L. Hood and S. J. Weidenschilling. Meteoritics & Planetary Science 39, 1809-1 821, 2004; and 7) Radial drift of particles in the solar nebula: Implications for planetesimal formation. Icarus 165, 438-442, 2003.

Weidenschilling, Stuart J.↗
DOE OSTI2022

Energetic Materials

Energetic materials comprise explosives, pyrotechnics, and propellants. The science of energetic materials is dedicated to developing a means to predict performance and safety characteristics with high fidelity. This is a particular challenge and is predicated on materials science and engineering, physics, chemistry, and dynamic response in extreme conditions. Fundamental elements of these complicated composite materials remain grand challenges—from the design of high-energy metastable molecules, to the engineering of composite formulations, to the processing parameters that link to safety and performance characteristics in as-yet undetermined ways. Key elements include crystalline mechanics, grain dynamics, multiphase interfaces, thermal and mechanical damage, and failure—all linked to multistep and high-rate chemistry and shock physics. A future revolution in our understanding and predictive capability for energetic materials behavior and responses is dependent upon sustained focus and advances in materials research and development.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
NASA NTRS2004

Chondrites and the Protoplanetary Disk, Part 2

Contents include the following: On the Dynamical Evolution of a Nebula and Its Effect on Dust Coagulation and the Formation of Centimeter-sized Particles. The Mineralogy and Grain Properties of the Disk Surfaces in Three Herbig Ae/Be Stars. Astrophysical Observations of Disk Evolution Around Solar Mass Stars. The Systematic Petrology of Chondrites: A Consistent Approach to Assist Classification and Interpretation. Understanding Our Origins: Formation of Sun-like Stars in H II Region Environments. Chondrule Crystallization Experiments. Formation of SiO2-rich Chondrules by Fractional Condensation. Refractory Forsterites from Murchison (CM2) and Yamato 81020 (CO3.0) Chondrites: Cathodoluminescence, Chemical Compositions and Oxygen Isotopes. Apparent I-Xe Cooling Rates of Chondrules Compared with Silicates from the Colomera Iron Meteorite. Chondrule Formation in Planetesimal Bow Shocks: Physical Processes in the Near Vicinity of the Planetesimal. Genetic Relationships Between Chondrules, Rims and Matrix. Chondrite Fractionation was Cosmochemical; Chondrule Fractionation was Geochemical. Chondrule Formation and Accretion of Chondrite Parent Bodies: Environmental Constraints. Amoeboid Olivine Aggregates from the Semarkona LL3.0 Chondrite. The Evolution of Solids in Proto-Planetary Disks. New Nickel Vapor Pressure Measurements: Possible Implications for Nebular Condensates. Chemical, Mineralogical and Isotopic Properties of Chondrules: Clues to Their Origin. Maximal Size of Chondrules in Shock-Wave Heating Model: Stripping of Liquid Surface in Hypersonic Rarefied Gas Flow. The Nature and Origin of Interplanetary Dust: High Temperature Components. Refractory Relic Components in Chondrules from Ordinary Chondrites. Constraints on the Origin of Chondrules and CAIs from Short-lived and Long-lived Radionuclides. The Genetic Relationship Between Refractory Inclusions and Chondrules. Contemporaneous Chondrule Formation Between Ordinary and Carbonaceous Chondrites. Chondrules and Isolated Grains in the Fountain Hills Bencubbinite. Implications of Chondrule Formation in a Gas of Solar Composition. Implications of Meteoritic Cl-36 Abundance for the Origin of Short-lived Radionuclides in the Early Solar System. Size Sorting and the Chondrule Size Spectrum. Comparative Study of Refractory Inclusions from Different Groups of Chondrites. In Situ Investigation of Mg Isotope Distributions in an Allende CAI by Combined LA-ICPMS and SIMS Analyses Photochemical Speciation of Oxygen Isotopes in the Solar Nebula.

Source record↗
NASA NTRS2002

Shocks in Dense Clouds in the Vela Supernova Remnant: FUSE

We have obtained 8 LWRS FUSE spectra to study a recently identified interaction of the Vela supernova remnant with a dense cloud region along its western edge. The goal is to quantify the temperature, ionization, density, and abundance characteristics associated with this shock/dense cloud interface by means of UV absorption line studies. Our detection of high-velocity absorption line C I at +90 to +130 km/s with IUE toward a narrow region interior to the Vela SNR strongly suggests the Vela supernova remnant is interacting with a dense ISM or molecular cloud. The shock/dense cloud interface is suggested by (1) the rarity of detection of high-velocity C I seen in IUE spectra, (2) its very limited spatial distribution in the remnant, and (3) a marked decrease in X-ray emission in the region immediately west of the position of these stars where one also finds a 100 micron emission ridge in IRAS images. We have investigated the shock physics and general properties of this interaction region through a focussed UV absorption line study using FUSE spectra. We have FUSE data on OVI absorption lines observed toward 8 stars behind the Vela supernova remnant (SNR). We compare the OVI observations with IUE observations of CIV absorption toward the same stars. Most of the stars, which are all B stars, have complex continua making the extraction of absorption lines difficult. Three of the stars, HD 72088, HD 72089 and HD 72350, however, are rapid rotators (v sin i less than 100 km/s) making the derivation of absorption column densities much easier. We have measured OVI and CIV column densities for the "main component" (i.e. the low velocity component) for these stars. In addition, by removing the H2 line at 1032.35A (121.6 km/s relative to OVI), we find high velocity components of OVI at approximately 150 km/s that we attribute to the shock in the Vela SNR. The column density ratios and magnitudes are compared to both steady shock models and results of hydrodynamical SNR modeling. We find that the models require the shock to be relatively slow (approximately 100 - 170 km/s) to match the FUSE data. We discuss the implications of our results for models of the evolution of the Vela SNR.

Nichols, Joy↗
DOE OSTI2022

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Project DarkStar: Vision for LLNL in 2030

DarkStar was a Strategic Initiative (FY2021-FY2024) to investigate applications of Artificial Intelligence (AI) and Machine Learning (ML) to scientific problems of complex hydrodynamics, shockwave physics and energetic materials. The research focused on physics and engineering design as a process that can be tremendously accelerated through merging AI with advanced physics simulation on exascale-class platforms, and to experimentally validate this revolutionary new approach through dynamic materials campaigns. A central thread of scientific inquiry was in the application of AI to enable human understanding of how to control hydrodynamic instability (which has impacts to areas such as inertial confinement fusion) via engineering features and time-dependent sources. Motivated by an unfinished line of research started by Dr. Johnny von Neumann, AI-enabled simulation approaches were developed that allowed DarkStar researchers to uncover several ground-breaking discoveries regarding hydrodynamic instability, including how to completely suppress Richtmyer-Meshkov instability (RMI). These S&T discoveries, along with other advances, have shown the way for an entirely new approach to time-dependent problems known as inverse design – the idea that complex systems can be developed directly from a final state that is to be achieved and resolve the initial design via satisfying several constraints simultaneously via AI/ML. Through experimental campaigns conducted across a wide range of facilities in the NNSA complex (the High Explosive Application Facility at LLNL, the Dynamic Compression Sector/Advanced Photon Source at Argonne National Lab, and Special Technologies Laboratory at MSTS) the radical new AI/ML approach to engineering complex material dynamics was verified, establishing a new field of study within the realm of shock physics. As advanced manufacturing capabilities continue to develop, the great importance of inverse design as a means to apply that technology effectively for NNSA missions will feature prominently over this decade. DarkStar has positioned NNSA as a world-leader in this newly emerging cross-disciplinary area of AI methods for advanced physics simulation and pioneered multiple novel approaches that have enabled the broader scientific community. By allowing us to see past the horizon, to 2030 and beyond, DarkStar has illuminated the vast potential of AI/ML to impact a wide range of new national security missions and, consequently, multiple areas of further research have already emerged across the NNSA and DOD complex.

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