Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock timing”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗
DOE OSTI2025

Fractal Scaling of Explosively Driven Product Gases

ABSTRACT Characterization of the interface between explosive product gases and ambient air in an explosion is a complicated task due to the turbulent mixing and inherently three‐dimensional expansion of the interface. This study aims to quantify the evolution of the interface as a temporally varying Hausdorff dimension. Two test series were conducted with Composition C‐4 charges with masses of 105 and 880 g. Imaging data were collected from the time of detonation until shock wave detachment using ultra‐high‐speed cameras. Gas cloud profiles were extracted using automated image processing algorithms, and the Hausdorff dimension of these two‐dimensional slices of the gas cloud was then estimated using boxcounting algorithms. When scaled with standard gas dynamic nondimensional scalings, the Hausdorff dimension of all explosive events appears to collapse towards a single curve. The fireball was initially nonfractal and began to develop fractal properties as the shock wave separated from the detonation products. Artificial perturbation of the charged surface had no detectable impact on the evolution of the Hausdorff dimension in the early development of the fireball outside of error, despite visible phenomenological differences in the early development of mixing on the fireball surface.

42 ENGINEERING↗
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Atmospheric clearing of solid-particle debris using femtosecond filaments

Through nonlinear self-focusing, femtosecond pulses can propagate several kilometers beyond diffraction limits, forming an ionization channel in air known as a laser filaments. It has been demonstrated that in the wake of the filament, aerosols can be effectively cleared to improve the transmission of subsequent laser pulses or secondary light sources, pertinent to applications in atmospheric sensing. However, the current understanding of aerosol clearing is founded on interactions with droplets to simulate fogs and clouds and thus does not extend to solid particles or atmospheric debris. Using optical trapping, we isolate both graphite and silica microparticles and directly measure the subsequent displacement caused by the filament using time-resolved shadowgraphy. The shock wave from the filament is demonstrated to propel particles away from the filament, directly contributing to atmospheric debris clearing. Particles exposed to the laser light in either the intense filament core or the surrounding energy reservoir are axially displaced along the beam path. Furthermore, it is found that the optomechanical properties of the particle largely influence the axial displacement induced by laser exposure through mechanisms such as radiation pressure, mass ejection from ablation or optical damage, and particle deagglomeration.

42 ENGINEERING↗
DOE OSTI2025

​​Hydrogen Transportation and Distributed Energy Systems Seismic Risk Assessment for Cascadia Subduction Zone Airport Facilities​

Portland International Airport in Oregon is exploring operating a fleet of 28 fuel cell electric buses to support airport operations and possibly provide backup power during outages. In this report, we evaluate the risk of hydrogen deployment at an airport considering the potential for seismic activity. We present a seismic risk assessment for pressurized piping for a generic hydrogen refueling station at the airport. Only pressurized piping was considered in the risk analysis because of its vulnerability and availability of fragility information unlike the other components. Seismic capacity of other components should be included as more information becomes available. To characterize the seismic hazard, we considered the initiating event frequency of a Cascadia Subduction Zone earthquake event based on data from the literature. Seismic stressors that translate to the site from an earthquake event are expressed as peak ground acceleration. Site disturbance is characterized as a function of soil conditions, represented by different shear wave velocity values, which affect how seismic waves propagate and impact structures. Pipe fragility curves as a function of seismic capacity were used to correlate ground acceleration and failure probability. These fragility characteristics were combined with site-specific soil conditions to calculate the probability of pipe rupture during a seismic event. It should be noted that multiple simultaneous failures due to an earthquake event as a common cause are not considered. Some epistemic uncertainties such as self-exciting shocks, aftershocks, and the time variant nature of the ground motion are not considered either.

08 HYDROGEN↗
DOE OSTI2025

Autoignition Delay Time Measurements of Dimethyl Ether at 110 Bar Inside a Shock Tube

Fuel mixtures of propane and dimethyl ether have been proposed as potential diesel substitutes for use in compression-ignition engines. There is an interest in phasing out diesel fuel because its combustion byproducts are harmful to the environment. Propane and dimethyl ether burn cleaner and produce significantly less carbon dioxide than diesel. The combustion of these fuels at engine-relevant conditions must be characterized to create an engine capable of running on propane and dimethyl ether. In this work, IDT measurements of stoichiometric neat DME diluted in 75% argon were completed for temperatures spanning 800-1200K at 110 bar. A shock tube will be used for all experiments in this work. Experimental IDTs were then compared to simulated values using several well-established mechanisms in the literature.

Guenther, Aaron C. [University of Central Florida]↗
DOE OSTI2026

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗
DOE OSTI2020

Influence of the Interplanetary Magnetic Field Cone Angle on the Geometry of Bow Shocks

The frozen-in interplanetary magnetic field (IMF) in the solar wind is one of the most important parameters affecting the Earth’s space weather. In the early studies of the IMF’s influence on space weather, significant effects of the north–south component of the IMF B {sub Z} were emphasized, while the radial component of the IMF B {sub X} was largely ignored. However, the IMF near the Earth is not always dominated by the north–south component of the IMF B {sub Z}, and the radial component of the IMF B {sub X} also plays an important role. However, while the effects of the IMF B {sub X} (cone angle) on the magnetopause have been studied in recent years, there has been much less effort to quantify the B {sub X} (cone angle) effects on the bow shock. In this paper, using the bow shock crossing data from multiple satellites, we investigate the IMF cone angle effect on the dayside and nightside of the bow shock. Our results show that under the radial IMF condition, the dayside of the bow shock is located closer to the Earth than the average. At the same time, on the nightside, the bow shock is farther away from the Earth than the average. The mechanism explaining the bow shock location under the radial IMF is not completely understood. We believe that the magnetosonic Mach number and unusual conditions of the magnetosheath, especially for low dynamic pressure, play an important role. In the future, more work is needed to describe the reactions of the Earth’s magnetosphere to different IMF orientations, especially to the radial IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Shock compression of additively manufactured high solids loaded polymer composites with aligned porosity

The effects of aligned porosity on the shock-compression response of additively manufactured (AM) composite structures with collinear and log-cabin geometries are investigated. The composites contain two inorganic and two organic particle populations and a UV-curable polymer binder. X-ray phase contrast imaging (X-PCI) is used as an in-material diagnostic in plate-impact experiments performed at the Dynamic Compression Sector of the Advanced Photon Source at the Argonne National Laboratory. Macroscale shock and particle velocities and microscale pore collapse velocities are measured from X-PCI images obtained from samples impacted in various orientations relative to the AM build directions. The linear shock and particle velocity relation for the structures with collinear build geometry shows directional differences, with the build direction exhibiting the highest slope and the filament print direction having the shallowest slope. Composites with log-cabin geometry, due to their slightly lower porosity and smaller sized pores, show higher shock velocities than the collinear structures at the same particle velocity. Microscale directionality effects are also observed, with the pore collapse velocity dependent on the pore configuration. Specifically, the pore collapse velocity is highest for pores elongated along the impact direction which is the case for the colinear structures impacted in the filament print direction, which is consistent with the directionality trend observed in the data for the shock and particle velocity equation of state. The pore collapse velocity approaches ∼80 % of the shock velocity and is several times the average bulk particle velocity, suggesting a highly hydrodynamic mechanism of pore collapse. Overall, both the macroscale response as measured by the shock and particle velocity relationship, and the microscale response as measured by the pore collapse velocity, are found to be anisotropic and dependent on the AM composite structure.

Wagner, K. B. [Georgia Institute of Technology, At↗
DOE OSTI2022

Shock Tube and Flame Speed Measurements of 2,4,4-Trimethyl-1-Pentene: A Co-Optima Biofuel

Abstract The combustion of 2,4,4-trimethyl-1-pentene (diisobutylene, C8H16), which is a biofuel and a component of surrogate fuels, is examined in this work. Carbon monoxide time–histories and ignition delay times are collected behind reflected shock waves utilizing a shock tube and mid-infrared laser absorption spectroscopy. Measurements were obtained near 10 atm pressure during stoichiometric oxidation of 0.15%C8H16/O2/Ar. Simulated results from chemical kinetic models are provided, and sensitivity analyses are used to discuss differences between models for both ignition delay times and carbon monoxide formation. In addition, laminar burning speeds are obtained at 1 atm, 428 K, and equivalence ratios, phi, between 0.91 and 1.52 inside a spherical chamber facility. Measured burning speeds are found to be less than that of ethanol over the equivalence ratio span. Burning speed measurements are compared to predictions of chemical kinetic mechanisms and are in agreement for the richest conditions; however, at lean conditions, the model predicts a far slower-burning speed. The maximum burning speed occurs at an equivalence ratio of 1.08 with a magnitude of 0.70 m/s. The current work provides the crucial experimental data needed for assessing the feasibility of this biofuel and for the development of future combustion chemical kinetics models.

Energy & Fuels↗
DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2026

Development of an ultrafast wide-field microscope for single-shot laser-driven shock and ablation studies

We build an ultrafast wide-field microscope that captures single-shot events at near-optical diffraction limit for time-resolved studies of ablation and shock compression. Here, we provide insights on how the imaging and spectroscopic modalities of our single-shot microscope are impacted by the depth of field of a microscope objective. At high spatial resolution, the depth of field can critically impede these experiments and novel strategies must be implemented to acquire high-quality optical micrographs at high throughput. Here, we highlight how custom lithography-prepared samples enable us to rectify optical artifacts that arise during sample motion. With a complete knowledge of how the mechanical and optical components couple together, we can then extract critical spatial and temporal features of our ablation and shock compression experiments. For example, during the ablation processes, we show that shock waves created in metal thin films form Newton’s interference rings, whose evolution can be recorded with picosecond temporal resolution and at near optical diffraction-limit resolution.

Ultrafast microscopy↗
DOE OSTI2020

Single-pulse (100 ps) extended x-ray absorption fine structure capability at the Dynamic Compression Sector

Determining real-time changes in the local atomistic order is important for a mechanistic understanding of shock wave induced structural and chemical changes. However, the single event and short duration (nanosecond times) nature of shock experiments pose challenges in obtaining Extended X-ray Absorption Fine Structure (EXAFS) measurements - typically used for monitoring local order changes. Here, we report on a new single pulse (~100 ps duration) transmission geometry EXAFS capability for use in laser shock-compression experiments at the Dynamic Compression Sector (DCS), Advanced Photon Source. We used a flat plate of highly oriented pyrolytic graphite (HOPG) as the spectrometer element to energy disperse x-rays transmitted through the sample. It provided high efficiency with approximately 15% of the x-rays incident on the HOPG reaching an x-ray area detector with high quantum efficiency. This combination resulted in good signal-to-noise ratio (~10 3 ), energy resolution of ~10 eV at 10 keV, EXAFS spectra covering 100s of eV, and good pulse to pulse reproducibility of our single pulse measurements. Ambient EXAFS spectra for Cu and Au are compared to reference spectra, validating our measurement system. Comparison of single pulse EXAFS results for ambient and laser shocked Ge(100) shows large changes to the local structure of the short lived state of shocked Ge. The current DCS EXAFS capability can be used to perform single pulse measurements in laser shocked materials from approximately 9-13 keV. These EXAFS developments will be available to all users of the DCS.

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