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 37 records · Page 2

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

CTH v.13.0

SAND2023-06647O The latest version of the CTH suite of computer codes is a flexible software system designed to treat a wide range of shock wave propagation and material motion phenomena in 1-D, 2-D, or 3-D. CTH has models suitable for most conditions encountered in shock physics including material strength, fracture, distended materials, high explosives, and a variety of boundary conditions. Finite-volume analogs of the Lagrangian equations of momentum and energy conservation are employed with continuous rezoning to construct Eulerian differencing. Adaptive mesh refinement is available for maximizing resolution in regions of interest. Full Lagrangian numerical methods are available for mechanical or structural response.

Schumacher, Shane↗

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↗

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↗

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↗

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↗

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↗

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↗