Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock waves”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2025

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Unraveling the role of temperature on the onset of ejecta formation at atomic scales

Shock wave interactions with perturbations on a free surface can lead to the inversion and growth of the perturbation and eventual ejection of a jet of material, referred to as an ejecta microjet, from shocked surfaces. Here, this study carries out large-scale molecular dynamics (MD) simulations using single-crystal Cu and Sn systems with a pre-existing groove to characterize the localized gradients in temperatures and pressures generated that render the microjet formation. MD simulations are carried out for three loading orientations (along the [001], [011], and [111] directions) and shock pressures ranging from 16 GPa to 100 GPa to understand the role of temperatures generated at the groove vertex on the formation of ejecta. The simulations suggest that the interaction of the shock wave with the groove results in a localized increase in temperatures, leading to localized softening at the groove vertex and the generation of an ejecta microjet. For Cu systems, the simulations suggest that the orientation effects on shock wave structures, velocities, and localized softening affect the ejecta formation at low pressures. Jetting is only observed when temperatures at the groove vertex are high enough to induce localized softening (close to melting temperature), and the jet velocity increases with shock pressure. In contrast, the loading orientation is rendered inconsequential for Sn systems due to the melting of the material during shock compression at the pressures chosen. The jet velocities are similar, regardless of crystal orientation in Sn systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗
DOE OSTI2025

Absolute Hugoniot measurements in low-density plastic foams on the NIKE shock compression platform

The propagation of Mbar-range shock waves in low-density foam materials is of interest to target design in all approaches to inertial confinement fusion, high energy density physics, and laboratory astrophysics. We report absolute Hugoniot measurements for CH 1.72 N 0.086 O 0.37 plastic foams with 73 and 94 mg/cm 3 densities in the 32–107 km/s shock velocity range. The experiments were performed on the shock compression platform developed on the NIKE KrF laser facility at the U.S. Naval Research Laboratory. NIKE's 4 or 8-ns long flat-top laser pulses drive steady shock waves into foam targets at the ablative pressures of 1–7 Mbar. The propagation of the ablation and shock fronts is tracked continuously in time using streaked side-on monochromatic x-ray imaging radiography. The straight x−t trajectories of the shock and ablation fronts in the recorded streak images confirmed their steadiness. The SESAME tabulated equation of state predictions generally agrees with our Hugoniot data within a 95% confidence band. The experimental uncertainty of the evaluated shock density compression ratios remains large, indicating the need for more experiments and improved theoretical understanding of the strong shock propagation mechanisms in dry foams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of a collapsing gas bubble on the shock-to-detonation transition in liquid nitromethane

We studied the shock-induced collapse of butane gas bubbles in the homogeneous explosive nitromethane (NM) to investigate the effects of hot spot formation on the detonation process. A butane bubble was injected into a sample of NM, and a shock wave from a flat plate impactor compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and bubble collapse processes. A multiband optical fiber pyrometer measured the time-resolved thermal radiance, and we used the results and emissivity values extracted from spectral fits to estimate temperatures. We measured the characteristics of the shock-to-detonation transition in NM with and without a bubble. All experiments were performed at shock pressures near 8 GPa, where neat NM can detonate. A single bubble in this system was shown to sensitize NM, leading to a reduced run-to-detonation time. We used hydrodynamic modeling to predict shock wave propagation, the extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. We used a temperature-dependent Arrhenius burn model for simulations, and it yielded much better results than reactive burn models that depend only on pressure and density.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Optimizing infrasound observations for sample return capsule re-entry: Insights from OSIRIS-REx and Hayabusa2

The atmospheric entry of meteoroids presents a rare and unpredictable phenomenon, posing challenges for systematic observation and detailed characterization. Such events are nonetheless critical for advancing understanding of acoustic wave propagation, atmospheric structure, and entry dynamics. In contrast, sample return capsules (SRCs) from space missions follow well constrained re-entry trajectories, enabling planned observations of shock wave generation and propagation under controlled conditions. This study compares two SRC atmospheric entries, Hayabusa2 in 2020 and OSIRIS-REx in 2023, to assess how different infrasound array configurations influence shock wave detection and trajectory validation. Hayabusa2’s re-entry was monitored using a distributed network of 28 portable infrasound sensors across seven arrays in Woomera, Australia, permitting three-dimensional reconstruction of the trajectory and analysis of wave characteristics. For OSIRIS-REx, a compact four-sensor array deployed near Eureka Airport provided trajectory confirmation through arrival-time differences and back-azimuth estimates. Spectral and waveform analyses revealed differences in signal properties associated with variations in entry angle and velocity. The results illustrate both the strengths and the limitations of the deployed array configurations. The dense and distributed arrays during Hayabusa2’s re-entry enabled detailed trajectory reconstruction, whereas the compact array at Eureka primarily provided confirmation of signal coherence and back-azimuth consistency. These case studies highlight that even relatively small arrays, if located close to the predicted ground track, can still capture useful information on arrival direction and timing. Such insights provide practical guidance for planning future observational campaigns of SRC returns and other controlled atmospheric entries, and may also inform approaches to opportunistic observations of meteoroid events and other atmospheric acoustic phenomena.

celestial mechanics: orbit determination↗
DOE OSTI2025

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
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 OSTI2025

Array-Based Seismic Measurements of OSIRIS-REx’s Re-Entry

The return home of the OSIRIS-REx spacecraft in September 2023 marked only the fifth time that an artificial object entered the Earth’s atmosphere at interplanetary velocities. Although rare, such events serve as valuable analogs for natural meteoroid re-entries; enabling study of hypersonic dynamics, shock wave generation, and acoustic-to-seismic coupling. Here, in this study, we report on the signatures recorded by a dense (100 m scale) 11-station array located almost directly underneath the capsule’s point of peak atmospheric heating in northern Nevada. Seismic data are presented, which allow inferences to be made about the shape of the shock wave’s footprint on the surface, the capsule’s trajectory, and its flight parameters.

58 GEOSCIENCES↗
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 OSTI2026

Discovery of a nearby radio relic in the low-mass, merging cluster Abell 4067

Shock waves generated during cluster mergers offer a powerful probe of how large-scale structure grows and evolves in the Universe. As part of the MeerKAT-South Pole Telescope (SPT) survey, we report the discovery of a single arc-like radio relic in the galaxy cluster Abell 4067 ($z=0.099$), one of the lowest-mass clusters known to host such a structure. MeerKAT UHF-band (0.58–1.09 GHz) observations reveal a relic with a largest linear size of $\sim 1.48\pm 0.02$ Mpc, located at a projected distance of 0.95 Mpc from the cluster centre. XMM–Newton X-ray data show that the relic’s position and orientation relative to the intracluster medium (ICM) elongation are consistent with a merger-driven shock-wave scenario. The relic has an estimated radio power of $3.10\pm 0.03\times 10^{24}$ W Hz$^{-1}$ at 150 MHz. When placed in the $P_{150\, \mathrm{MHz}}$–$M_{500}$ scaling relation, the Abell 4067 relic appears less luminous compared to relics in more massive clusters, suggesting an association with weak merger shocks. This finding supports the idea that relics in low-mass clusters may form through less energetic merger events, leading to weak merger shocks. The latter is supported by the absence of a detectable central radio halo in Abell 4067, which reinforces the idea that luminous radio haloes are not a universal outcome of cluster mergers and highlights the role of cluster mass, merger energetics, and evolutionary stage in shaping diffuse radio emission in the ICM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multiwavelength Afterglow Analysis of GRB 221009A: Unveiling the Evolution of a Cooling Break in a Wind-like Medium

Gamma-ray bursts (GRBs) are the most energetic explosions in the Universe, and their afterglow emission provides an opportunity to probe the physics of relativistic shock waves in an extreme environment. Several key pieces for completing the picture of GRB afterglow physics are still missing, including jet properties, the emission mechanism, and particle acceleration. Here, we present a study of the afterglow emission of GRB 221009A, the most energetic GRB ever observed. Using optical, X-ray, and gamma-ray data up to approximately 2 days after the trigger, we trace the evolution of the multiwavelength spectrum and the physical parameters behind the emission process. The broadband spectrum is consistent with the synchrotron emission emitted by relativistic electrons with its index of p = 2.29 ± 0.02. We identify a break energy at keV and an exponential cutoff at GeV in the observed multiwavelength spectrum. The break energy increases in time from $16.0^{+7.1}_{-4.9}$ keV at 0.65 days to $46.8^{+25.0}_{-15.5}$ keV at 1.68 days, favoring a stellar-wind-like profile of the circumburst medium with k = 2.4 ± 0.1 as in ρ(r) ∝ r -k . The high-energy attenuation at around 0.4 to 4 GeV is attributed to the maximum of the particle acceleration in the relativistic shock wave. This study confirms that the synchrotron process can explain the multiwavelength afterglow emission and its evolution.

79 ASTRONOMY AND ASTROPHYSICS↗