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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Large N algebras and generalized entropy

We construct a Type II ∞ von Neumann algebra that describes the large N physics of single-trace operators in AdS/CFT in the microcanonical ensemble, where there is no need to include perturbative 1/N corrections. Using only the extrapolate dictionary, we show that the entropy of semiclassical states on this algebra is holographically dual to the generalized entropy of the black hole bifurcation surface. From a boundary perspective, this constitutes a derivation of a special case of the QES prescription without any use of Euclidean gravity or replicas; from a purely bulk perspective, it is a derivation of the quantum-corrected Bekenstein-Hawking formula as the entropy of an explicit algebra in the G → 0 limit of Lorentzian effective field theory quantum gravity. In a limit where a black hole is first allowed to equilibrate and then is later potentially re-excited, we show that the generalized second law is a direct consequence of the monotonicity of the entropy of algebras under trace-preserving inclusions. Finally, by considering excitations that are separated by more than a scrambling time we construct a “free product” von Neumann algebra that describes the semiclassical physics of long wormholes supported by shocks. We compute Rényi entropies for this algebra and show that they are equal to a sum over saddles associated to quantum extremal surfaces in the wormhole. Surprisingly, however, the saddles associated to “bulge” quantum extremal surfaces contribute with a negative sign.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Ray-tracing image simulations of transparent objects with complex shape and inhomogeneous refractive index

Optical images of transparent three-dimensional objects can be different from a replica of the object’s cross section in the image plane, due to refraction at the surface or in the body of the object. Simulations of the object’s image are thus needed for the visualization and validation of physical models. We report ray tracing image simulations that achieved high physical fidelity, reproducing optical behaviors and image features not rendered in previous studies. We replicated brightfield microscopy images of drops with complex shapes, and images of pressure and shock waves traveling inside them. For high physical fidelity, the simulations must replicate the spatial and angular distribution of illumination rays, and both the experiment and the simulation must be designed for accurate optical modeling. The simulations are highly sensitive to the properties of the drops and can be used to diagnose and refine fluid dynamics models. The simulated images can also be optimized to extract multiple 3D properties from experimental images. Compared to specialized single-shot 3D imaging methods, this approach has the advantage that it preserves the experimental simplicity, the high resolution, and the visual interpretability characteristic to basic optical imaging. The techniques introduced here are directly applicable to optical microscopy, so they can be used in other fields, such as microfluidics and biology, to expand the type and the accuracy of three-dimensional information that can be extracted from basic optical images.

Cavitation↗
DOE OSTI2024

Diamond under extremes

Diamond is, by virtue of the covalent bonding between atoms and the very strong carbon to carbon bonds, the hardest natural material. It has been a fascinating material since its discovery, first as a decorative gem and more recently, for its numerous industrial uses because of its extreme hardness, elastic modulus, and optical transparency. In recent years, it has become a preferred ablator for laser shock experiments, and this has led to its choice as the capsule material for fusion experiments at the National Ignition Facility. Further, this review covers both experimental and computational (including machine learning) advancements in research on diamond subjected extreme conditions of temperature and pressure. The synergy between shock and ramp loading experiments and atomic level simulations is proving to be powerful in advancing our understanding of diamond under extremes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Probing shocks interacting with radiation waves with the Radishock experiment

Both radiation flows and shocks have been extensively studied in the laboratory in the past few decades due to their critical roles in many astrophysical and high-energy density physics processes. In the Radishock experiment, a halfraum-powered radiation wave is driven into a low-density foam and interacts with an ablatively driven, counter-propagating shock. The interacting waves produce a spike in energy density with a temperature greater than the local temperature of the individual waves. As in the successful predecessor experiment, COAX, the primary diagnostic uses absorption spectroscopy at many locations down the cylindrical target, enabling a spatial temperature inference of the radiation wave and its interactions with the shock. Combined with a radiography diagnostic that is capable of imaging the shock and interaction features, we are able to study and inform model predictions of the interaction spike phenomenon. We describe the underlying physics behind the shock interactions with the radiation front and the implications of this experimental study for a broad range of astrophysical phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Collisionless ion-electron energy exchange in magnetized shocks

Energy partition between ions and electrons in collisionless shocks has been a long-standing unsolved fundamental physical question. Here, we show that kinetic simulations of moderate Alfv´enic Mach number, magnetized, collisionless shocks reveal rapid, faster-than-Coulomb, energy exchange between ions and electrons when the plasma is sufficiently magnetized. Using kinetic and multi-fluid models with counter-streaming ions, we identify resonances between electron whistler and ion magnetohydrodynamic waves that account for this rapid energy exchange.

High-energy-density plasmas↗
DOE OSTI2025

RaDIATE Collaboration Thermal shock studies

As next-generation accelerator target facilities (High Energy Physics, Spallation Sources, ... ) become increasingly more powerful and intense, high power target systems face key technical challenges, such as radiation damage and thermal shock. Those ultimately degrade the performance and lifetime of targets and have been identified as the leading cross-cutting challenges of high-power target facilities. In order to operate reliable beam-intercepting devices in the framework of energy and intensity increase for next generation accelerators, the RaDIATE Collaboration (Radiation Damage In Accelerator Target Environment), established in 2012 and managed by Fermilab, brings together existing expertise in nuclear material and accelerator targets from 20 international institutions, including CERN, to execute a coordinated strategy for high power targetry R&D. In this context, several thermal shock studies were performed at CERN's HiRadMat (High-Radiation to Materials) facility, that took a key step towards improving our knowledge on target damage tolerance. The first experiment HRMT-24, supported by EURCARD2 and completed in 2015, successfully validated the Johnson-Cook strength model developed at SwRI on beryllium S200FH (used for beam window material), providing a better confidence in simulating the thermal shock response of current and future S200FH beryllium components. HRMT-43, supported by ARIES and completed in 2018, tested various materials (Be, C, SiC, Si, Ti and ceramic nanofiber). It was a first and unique test which included pre-irradiated specimens from high energy proton beam irradiation to identify thermal shock response differences between non-irradiated and previously irradiated materials. Real-time measurement of dynamic thermomechanical response of graphite helped to benchmark numerical simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗
DOE OSTI2024

Formation of collisionless shocks driven by strongly magnetized relativistic electrons in the laboratory

In experiments performed with the OMEGA EP laser system, proton deflectometry captured magnetic field dynamics consistent with collisionless shock formation driven by strongly magnetized relativistic electrons. During laser-foil interactions, shocks can form as relativistic electrons and strong surface magnetic fields generated by a short-pulse laser impinge on a cooler plasma produced by a longer-pulse laser. Three-dimensional particle-in-cell simulations reproduce the magnetic draping and fast formation speeds measured in the experiment and reveal that this relativistic-electron-driven shock forms at an interface that is unstable to shear and streaming instabilities. The simulation results provide insight into the microphysics that may influence high-energy shocks observed in extreme astrophysical environments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Refractive Imaging of Air Shock Above Microscale Defects in Pentaerythritol Tetranitrate (PETN) Films

Physical vapor deposition (PVD) of high explosives can produce energetic samples with unique microstructure and morphology compared to traditional powder processing techniques, but challenges may exist in fabricating explosive films without defects. Deposition conditions and substrate material may promote microcracking and other defects in the explosive films. Here, we investigate effects of engineered microscale defects (gaps) on detonation propagation and failure for pentaerythritol tetranitrate (PETN) films using ultra-high-speed refractive imaging and hydrocode modelling. Observations of the air shock above the gap reveal significant instabilities during gap crossing and re-ignition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fast imaging of laser ablation of cerium: Dependence of plasma plume dynamics and structure on native oxide growth

We have imaged the evolution of plasma following nanosecond pulsed laser ablation of cerium in air with sub-10 ns resolution. We find the ablated plasma plume expansion and shape to be strongly dependent on the time after mechanical polishing of the sample surface. After an initial shock front common to both new and aged samples (velocities up to 16 km s ₋1 ), we observe a second front that is relatively localized to the surface vicinity of newly polished samples, but in the case of older samples it moves outward at velocities of up to 4 km s ₋1 . We attribute this behavior to the formation and growth of a native oxide layer on the order of hours after polishing. These results demonstrate that plasma imaging can be used to diagnose the presence of, and evaluate the extent of, thin surface oxides or other compounds that form soon after exposure to some atmosphere. These findings also highlight the need to explicitly consider the presence of such overlayers when modeling laser/metal interaction, when making measurements of vapor phase chemistry, and when analyzing the composition, phase, and morphology of solid particulates formed after the ablation of reactive metals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigating growth to detonation in vapor-deposited hexanitrostilbene and pentaerythritol tetranitrate films using high-throughput methods

In this work, a high-throughput experimental setup was used to characterize initiation threshold and growth to detonation in the explosives hexanitrostilbene (HNS) and pentaerythritol tetranitrate (PETN). The experiment sequentially launched an array of laser-driven flyers to shock samples arranged in a 96-well microplate geometry, with photonic Doppler velocimetry diagnostics to characterize flyer velocity and particle velocity at the explosive–substrate interface. Vapor-deposited films of HNS and PETN were used to provide numerous samples with various thicknesses, enabling characterization of the evolution of growth to detonation. One-dimensional hydrocode simulations were performed with reactions disabled to illustrate where the experimental data deviate from the predicted inert response. Prompt initiation was observed in 144 μm thick HNS films at flyer velocities near 3000 m/s and in 125 μm thick PETN films at flyer velocities near 2400 m/s. This experimental setup enables rapid quantification of the growth of reactions in explosive materials that can reach detonation at sub-millimeter length scales. These data can subsequently be used for parameterizing reactive burn models in hydrocode simulations, as discussed in Paper II [D. E. Kittell, R. Knepper, and A. S. Tappan, J. Appl. Phys. 131, 154902 (2022)].

36 MATERIALS SCIENCE↗
DOE OSTI2024

High strain-rate strength response of single crystal tantalum through in-situ hole closure imaging experiments

The properties of crystalline materials often depend on directionality and operating conditions. Specifically, the strength of materials can depend anisotropically on crystal direction and the loading condition. To probe these effects, a preliminary series of high strain-rate (> 105/s) strength plate-impact hole closure experiments were performed on high purity single crystal Tantalum cubes. The orientation of the single crystals with respect to impact/loading were varied to provide data to inform crystal plasticity modeling efforts. The experiments consist of in-situ high-resolution X-ray radiographic imaging of the hole collapse under dynamic compression conditions to infer the material strength via its resistance to closure at increasing levels of plastic strain. The experiments are compared against hydrocode simulation predictions. Here, a comparison with simple elastic perfectly plastic strength model predictions is presented to elucidate the response of the different crystal orientations at high strain-rate and large plastic strains.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Effects of Climate Change on Interregional Electricity Market Dynamics on the U.S. West Coast

The United States (U.S.) West Coast power system is strongly influenced by variability and extremes in air temperatures (which drive electricity demand) and streamflows (which constrain hydropower production). As hydroclimate changes across the West Coast, a combination of forces may work in tandem to make its bulk power system more vulnerable to physical reliability issues and market price shocks. In particular, a warmer climate is expected to increase summer cooling (electricity) demands and shift the average timing of peak streamflow (hydropower production) away from summer to the spring and winter, depriving power systems of hydropower when it is needed the most. Here, we investigate how climate change could alter interregional electricity market dynamics on the West Coast, including the potential for hydroclimatic changes in one region (e.g. Pacific Northwest (PNW)) to “spill over” and cause price and reliability risks in another (e.g. California). We find that the most salient hydroclimatic risks for the PNW power system are changes in streamflow, while risks for the California system are driven primarily by changes in summer air temperatures, especially extreme heat events that increase peak system demand. Altered timing and amounts of hydropower production in the PNW do alter summer power deliveries into California but show relatively modest potential to impact prices and reliability there. Instead, it is future climate conditions (extreme heat) in California that may exert a stronger significant influence on prices and reliability in the PNW, especially if California continues to rely on its northern neighbor for imported power to meet higher summer demands.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Scaling arguments for a plasma experiment relevant to an interplanetary coronal mass ejection

We present the design of a plasma experiment scaled to be relevant to an interplanetary coronal mass ejection (ICME) traveling from the Sun to the Earth. The experiment would launch a compact torus of plasma (the scaled analog of the coronal mass ejected by the Sun) into a magnetized background plasma (the scaled analog of the interplanetary medium, the diffuse magnetized plasma that fills the Solar System). We present the complete process of scaling the ICME system to a laboratory experiment: beginning with the ICME features that we aim to replicate with a laboratory experiment, explaining how these desired features are translated into specific design constraints, and presenting a set of experimental parameters that meet those constraints. We find that the experiment would be possible with plasma conditions typically achieved in large-scale facilities such as the Big Red Ball at the University of Wisconsin or the LArge Plasma Device at the University of California, Los Angeles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗