Hypersonic Shock Wave-Boundary Layer Interaction in a Shock Tunnel using Nitric Oxide Laser Induced Fluorescence
Explore the source record for details and available documents.
Jans, Elijah Rhea [Sandia National Laboratories (S↗
SEARCH · Search NASA
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.
Heat Loss Correction Factor for Fireball Yield Measurements
Fireball yield calculations are performed using Taylor’s radius equation. However, at early times, the x-ray diffusion process has a significant effect on the size of the fireball, and at late times the shock wave degenerates into an acoustic wave. Thus, earlier research determined that the range of validity for Taylor’s equation was in the scaled time regime of 0.004–0.008 s. So, when performing a yield calculation, the only values that are analyzed are those that lie within that range of validity. While Taylor’s equations offer an approximation of the scaled yield, they also rely on several assumptions. For example, the equations assume that the shock wave is adiabatic, when, in fact, the first light pulse on the nuclear test films, which occurs as the shock wave is forming, suggests that the shock wave is not adiabatic. Furthermore, the amount of heat loss from the shock wave that occurs prior to the scaled time regime of 0.004–0.008 s can change significantly as a function of air density. This may cause the shock wave radius to be smaller relative to the assumed adiabatic shock wave. Therefore, the objective of this project was to determine the correction factor for heat transfer, which should be applied in order to produce a more accurate approximation of the weapon’s yield.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
X-ray phase-contrast imaging of strong shocks on OMEGA EP
The ongoing improvement in laser technology and target fabrication is opening new possibilities for diagnostic development. An example is x-ray phase-contrast imaging (XPCI), which serves as an advanced x-ray imaging diagnostic in laser-driven experiments. In this work, we present the results of the XPCI platform that was developed at the OMEGA EP Laser-Facility to study multi-Mbar single and double shocks produced using a kilojoule laser driver. Two-dimensional radiation-hydrodynamic simulations agree well with the shock progression and the spherical curvature of the shock fronts. It is demonstrated that XPCI is an excellent method to determine with high accuracy the front position of a trailing shock wave propagating through an expanding CH plasma that was heated by a precursor Mbar shock wave. Furthermore, the interaction between the rarefaction wave and the shock wave results in a clear signature in the radiograph that is well reproduced by radiation-hydrodynamic simulations.
47 OTHER INSTRUMENTATION↗
Shock Response of Polymethyl Methacrylate (PMMA) Under Explosive Loading
Abstract The dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.
Torres, S. M. (ORCID:0009000287177206)↗
Temporally and spatially resolved X-ray densitometry in a shock tube
Gas densities were measured by x-ray absorption spectroscopy to examine spatial inhomogeneity in density behind incident and reflected shock waves in a miniature (12.7 mm bore) high repetition rate shock tube. Shock waves were generated in argon. The pressure and temperature behind the incident shock were P 2 ~ 1.72 bar and T 2 ~ 800 K, while those behind the reflected shock were P 5 ~ 6.5 bar and T 5 ~ 1480 K. Time-resolved line-of-sight transmission measurements using x-rays at 9 keV photon energy were made along various axial and transverse locations covering the entire shock tube cross section. The x-ray transmission at each location was converted to pathlength-integrated densities, which in some conditions could be further converted to pathlength-averaged densities. The measured gas densities were compared with those calculated by the normal shock wave relations, and good agreement was found in the preshock and immediately behind the incident shock regions. Similar agreement was found for reflected shock conditions very near the endwall. However, increasingly large differences between the measured and calculated gas densities were found as distance from the endwall increased. Furthermore, reconstruction by Abel inversion allowed time and radially resolved densities to be obtained from the transverse measurements. While the profiles measured in this experiment are of insufficient signal-to noise level and resolution to definitively assign boundary layer thicknesses, future improvements may be able to do this. The radially resolved densities reveal growing sidewall thermal boundary layers in gases behind the reflected shock. These reconstructions were compared with a model of the thermal boundary layer and are consistent with values of the thermal boundary layer thickness being between 0.25 mm and 1 mm at 0.7 ms after shock reflection.
42 ENGINEERING↗
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↗
The effects of electrically exploding gold bridgewires into inert and explosive powder beds
Abstract The particle velocity created in beds of both low-density inert sugar and explosive PETN as a function of distance from an exploding bridgewire was measured using optical velocimetry and a silvered PMMA window. As expected, more violent bridge-bursts (from a greater-stored-energy capacitive discharge unit) resulted in greater particle velocities and a better supported compaction wave in sugar. In all cases, ramp waves, not shocks, were observed in the inert sugar. Large window velocities were observed for very powerful bursts (up to 270 m/s), but bursts required for stochastic detonator operation conditions resulted in sugar/PMMA window velocities of only 8–10 m/s 0.85 mm from the bridge location. In contrast, after a distance of only 0.65 mm, a building shock wave was observed in PETN under both threshold and reliable firing conditions. Subsequently a hot-spot-driven shock-to-detonation (SDT) process was observed prior to full detonation. The measured buildup process accounts for $$\approx $$ ≈ 66% of the so-called excess transit time (ETT) between the observed and theoretical total function time for the particular exploding-bridge-wire (EBW) detonator studied. The remainder must occur in the powerful output pellet region. In contrast to a common understanding, the ETT is found to be a weak function of the discharge energy. Thus, the operation of the detonator after a bridge-burst energy-to-powder reaction transition process is found to be hot-spot-driven SDT in both the low- and high-density pellets.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Shock Compression of Materials: SSDDS Lecture 4
Lecture roadmap: Definitions, explanation of shock wave, properties of a shock wave, relations to equation of state and quasistatic compression, calculations, measurements, and problems.
36 MATERIALS SCIENCE↗
Imaging perturbed shock propagation in powders
A novel experimental methodology is presented to study the deviatoric response of powders in shock regimes. The powders are confined to a cylindrical wedge volume, and a projectile-driven shock wave with a sinusoidally varying front propagates through the powder. The perturbed shock wave exhibits a damping behavior due to irreversible processes of viscosity and strength (deviatoric) of the powder with propagation through increasing powder thicknesses. The inclined surface of the wedge is polished and coated to establish a diffuse surface suitable for reflecting incident laser light into a high-speed camera imaging at 5 MHz. Images of the contrast loss upon shock wave arrival at the observation surface are post-processed for qualitative and quantitative information. New data of shock damping behavior with parameters of perturbation wavelength and initial shock strength are presented for powders of copper, tantalum, and tungsten carbide as well as their mixtures. We present the first full-field images showing additional spatial disturbances on the perturbed shock front that appear dependent on particle material and morphology.
47 OTHER INSTRUMENTATION↗
A numerical study of particle jetting in a dense particle bed driven by an air-blast
In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.
42 ENGINEERING↗
Symmetry and scaling in one-dimensional compressible two-phase flow
Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
From Interstellar Cloud to Star to Laboratory: Frontier HEDP Studies of Magnetized Colliding Plasma Flows with Strong Radiative Cooling (Final Report)
High-speed flows of gas are the norm in astrophysics where exploding stars, winds from the regions around black holes and streams of galactic material create “supersonic flows”. As the flows collide, strong shock waves are formed where gas is compressed and heated. The regions behind these shock waves are of great interest to scientists. New stars can form in these shock regions in the context of galaxies. In the context of planets like our Earth, the shocks waves come from the interaction of the solar wind with our planets’ magnetic field causing Aurora and placing astronauts and satellites in jeopardy. Thus, understanding supersonic flows of magnetized plasmas is essential for progress in many fields of astronomy and space physics (relevant to planets and the Earth). Remarkably, we can also create these kinds of flows in the laboratory using Z-pinches and laser driven experiments that include magnetic fields. The exploration of these forms of “High Energy Density Physics” (HEDP) flows is also critical in order to achieve Inertial Confinement Fusion, which is a step towards developing a sustainable source of Fusion energy.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai
The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.
54 ENVIRONMENTAL SCIENCES↗
System and method for rapid, high throughput, high pressure synthesis of materials from a liquid precursor
The present disclosure relates to a system and method for synthesis of condensed, nano-carbon materials to create nanoparticles. In one embodiment the system may have a source of liquid precursor, a flow control element and a shock wave generating subsystem. The flow control element is in communication with the source of the liquid precursor and creates a jet of liquid precursor. The shock wave generating subsystem drives a shock wave through at least a substantial portion of a thickness of the jet of liquid precursor to sufficiently compress the jet of liquid precursor, and to increase a pressure and a temperature of the jet of liquid precursor, to create solid state nanoparticles.
Armstrong, Mike↗
Hypersonic Shock Wave-Boundary Layer Interaction on a Cone-Slice-Wedge Model in a Shock Tunnel
Explore the source record for details and available documents.
Wagner, Justin Lawrence [Sandia National Laborator↗
Generation of shock trains in free liquid jets with a nanosecond green laser
Shock wave trains in liquid jets were previously generated only by ablation with femtosecond x-ray lasers. Here we show that shock trains in water microjets can be also generated using nanosecond green laser pulses with 1- to 10-mJ energy. Furthermore, the ablation of 15-, 20-, 30-, and 70-μm water microjets opened a gap in the jets and launched an initial shock wave. Fully developed shock trains were observed in the 30- and 70-μm jets up to 250-ns delays, and these trains were also transmitted inside the nozzles. A few tens of nanoseconds after the pulse, the shock dynamics and its pressure became similar to the ones generated by x-ray lasers, with a more rapid pressure decay in thinner jets. At time delays exceeding 100 ns in the 30-μm jets, the leading shock pressure stabilized to an approximately constant pressure of 40 MPa. The energy density deposited in the jets was estimated at 30 MJ/cm 3 by comparing the jet gaps in the green and x-ray laser experiments, and matched previous estimates for optical ablation in water. Here, the pressure decay in the 30-μm jets was modeled based on the pressure decay observed in x-ray laser experiments.
74 ATOMIC AND MOLECULAR PHYSICS↗
Modeling shock-induced void collapse in single-crystal Ta systems at the mesoscales
Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.
36 MATERIALS SCIENCE↗
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↗