Hypersonic Shock Wave-Boundary Layer Interaction Experiments on a Cone-Slice-Ramp over a Range of Enthalpies
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
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.
Shock experiments are widely used to understand the mechanical and electronic properties of matter under extreme conditions. However, after shock loading to a Hugoniot state, a clear description of the post-shock thermal state and its impacts on materials is still lacking. We used diffraction patterns from 100-fs x-ray pulses to investigate the temperature evolution of laser-shocked Al–Zr metal film composites at time delays ranging from 5 to 75 ns driven by a 120-ps short-pulse laser. We found significant heating of both Al and Zr after shock release, which can be attributed to heat generated by inelastic deformation. A conventional hydrodynamic model that employs (i) typical descriptions of Al and Zr mechanical strength and (ii) elevated strength responses (which might be attributed to an unknown strain rate dependence) did not fully account for the measured temperature increase, which suggests that other strength-related mechanisms (such as fine-scale void growth) could play an important role in thermal responses under shock wave loading/unloading cycles. Our results suggest that a significant portion of the total shock energy delivered by lasers becomes heat due to defect-facilitated plastic work, leaving less converted to kinetic energy. This heating effect may be common in laser-shocked experiments but has not been well acknowledged. High post-shock temperatures may induce phase transformation of materials during shock release. Another implication for the study is the preservability of magnetic records from planetary surfaces that have a shock history from frequent impact events.
Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.
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.
In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.
This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of $192$ $\mu$m. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier–Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant $C$ 2 = $2$ compared to the accepted value $C$ 2 = $2/3$. In conclusion, results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.
A material at high temperature can react or decompose. For an energetic material, the reaction is exothermic and releases chemical energy that would further increase the temperature. Under some circumstances, when a reaction is triggered, such a reaction can propagate and the material rapidly releases a large amount of energy giving rise to an explosion. Examples of such materials are aerosols, suspensions of solid particles or liquid droplets in a gas; such as coal dust, grain dust and fuel-air explosions. Frequently, explosions are due to accidents. A spectacularly destructive example is the recent explosion of a large quantity of ammonium nitrate (thousands of tons) in Beirut, Lebanon (August 2020); see for example Beirut explosion. Ammonium nitrate is used as a fertilizer. It and the aerosols are not considered to be explosives due to the limited conditions for which an explosion can occur. An aerosol gets the oxidizer from the surrounding air. Burning requires diffusion of the oxidizer to the particle surface where the reaction occurs. A large density of small particles is required for a fast enough reaction to support an explosion. In contrast, an explosive is an energetic material with both fuel and oxidizer mixed on a molecular scale (either premixed gases or within molecules of a solid). This allows fast enough reactions over a wide range of conditions to support a self-propagating reactive wave known as a detonation wave. A detonation wave can be controlled and an explosive used for useful purposes such as in mining, construction, demolition, explosive welding, argon flash lamp, pulsed power using a magnetic flux generator [see also Goforth et al., 2015], jet cutter with shaped charge, explosive art, and generating conditions to study the response of materials at high strain rates and high pressures [see for example, Marsh, 1980]. Explosives are also used in conventional munitions and nuclear weapons. The focus of this book is on the theory and phenomenology of solid high explosives (HEs); in particular, plastic-bonded explosives (PBXs). Some aspects of detonation wave theory are needed to interpret explosive data. Hence, the theory is presented before the detonation wave phenomenology. A familiarity with fluid flow, specifically the notion of shock waves and the shock loci are assumed. In the remainder of this chapter we give a brief overview on the basic properties of detonation waves and PBXs.
TATB (1,3,5-triamino-2,4,6-trinitrobenzene) is a powerful explosive whose dynamical behavior is difficult to study because TATB is so insensitive to initiation by shock waves. We used a tabletop microscope equipped with 0–4.5 km/s laser-launched flyer plates to study shock initiation of TATB, which was fabricated in the form of an array of hundreds of plastic-bonded explosive minicharges (X-TATB = 80% TATB + 20% Sylgard 182 polymer). The 4 ns shocks from the flyer plates were not effective in initiating TATB, but we also developed a two-layer array where flyers first initiated a plastic-bonded PETN (pentaerythritol tetranitrate) charge (X-PETN = 80% PETN + 20% Sylgard), which drove an initiating 25 ns shock into the X-TATB. Thermal emission from shocked X-TATB was used to measure time-dependent temperature profiles with a resolution of 2 ns and to produce high-speed (5 ns) videos. In X-TATB, flyer plates produced 2500–3500 K hot spots and combustion at 2500 K. With X-PETN initiators, X-TATB had 3500–4000 K hot spots and a powerful volume explosion lasting a few nanoseconds. Prospects for producing TATB detonations on a tabletop are discussed.
The engineered macro-pore foam provides a new way to study thermonuclear burn physics by utilizing capsules containing deuterated (D) foam and filling tritium (T) gas in the engineered macro-pores. The implosion of a thermonuclear capsule filled with an engineered macro-pore foam will be complex due to the interaction of a shock wave with the engineered macro-pores. It is our goal to quantify how substantially complex foam structures affect the shape of shock and bulk shock speed. A cylinder-shape shock tube experiment has been designed and performed at the Omega Laser Facility. In order to examine how a foam structure will affect shock propagation, we performed several tests varying (1) engineered macro-pore size, (2) average foam density, and (3) with/without neopentane (C 5 H 12 ) gas. X-ray radiographic data indicate that shock speed through engineered macro-pore foams depends strongly on average foam density and less on pore size. In this work, experimental shock propagation data helped guide two numerical simulation approaches: (1) a 2D simulation with homogenizing foams rather than explicitly simulating engineered macro-pores and (2) a 2D toroidal-pore approximation adopting a toroidal-tube geometry to model engineered macro-pores.
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.
Indirectly driven shock-tube experiments were performed on the Omega Laser Facility to evaluate the relative importance of hohlraum x ray and radiative shock preheat sources on a low-density foam. X rays emitted from the hohlraum and a subsequent shock wave are channeled into a low-density foam sample, which contains a plastic preheat-witness disk placed downstream of the foam. Simultaneous radiographic measurements of the shock speed in the foam and the expansion rate of the witness disk due to preheat allow for the observation of effects from the hohlraum's low-energy and high-energy x-ray spectrum. Here we showed, from simulations, that low-energy x rays from the hohlraum are preferentially absorbed near the ablator surface (where the hohlraum and the shock tube meet), while higher-energy x rays largely pass through the ablator and foam and are volumetrically absorbed by the witness disk. Reproducing the experimentally measured shock speed and expansion of the witness disk simultaneously, we extracted the temperature evolution of preheated foam from the simulation and evaluated the relative importance of preheat sources on a low-density foam from hohlraum x-ray radiation and radiative shock. We found that radiation from the shock front was more effective at preheating the low-density foam than the high-energy x rays from the hohlraum. This shock-tube preheat experiment is important for understanding the results of the MARBLE experiments at the National Ignition Facility because initial conditions of foam-filled MARBLE capsules are sensitive to preheat.
When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.
Temperature determination in shock compressed solids constitutes an important and long-standing scientific need. Since the reduction of Bragg diffracted peaks due to temperature increase, using the Debye-Waller factor, is well established, we examined the use of this approach to determine temperatures in shock compressed gold and platinum by representing the shocked state as a superposition of density and temperature changes. Comparison of the calculated and measured diffraction peaks did not show good agreement, because x-ray diffraction (XRD) profiles in the shock compressed state are not governed solely by density and temperature changes. XRD results are also influenced measurably by shock wave induced microstructural changes. Our results demonstrate that contributions from microstructural changes need to be incorporated and modeled in the theoretical analysis to use XRD measurements for reliable temperature determination in shock compressed solids.
We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 , and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.
Polymers are widely implemented materials often used for safety and security in our daily lives. However, we do not have a thorough understanding of the phase transformations of these materials when exposed to extreme condition, high strain rate, events. To study polymers in these severe environments, we performed plate impact experiments that impart a one-dimensional shock wave, at high pressures and temperatures, to observe the materials response. Generally, velocimetry is the experimental diagnostic used to capture the movement of high velocity interfaces to then calculate the mechanical variables associated with shock compression. However, when a material undergoes a phase transition, such as between two crystalline phases or melting to decomposition during shock compression, velocimetry does not explicitly indicate what structural change occurred. At the Dynamic Compression Sector at the Advanced Photon Source gas gun driven compression has been paired with x-ray diffraction, providing the ability to study both the shock parameters and structural changes concurrently. We explored high density polyethylene (HDPE) both on its reactants Hugoniot (<10.12 GPa) and within its transition region (24.3 GPa) where HDPE is hypothesized to decompose from reactants to products. With x-ray diffraction we observed both polymer chain compression and crystal structure reorientation from orthorhombic to monoclinic above 6 GPa, as well as either a possible phase change or chemical decomposition at 24.3 GPa. With these experiments we were able to determine the crystalline and amorphous Hugoniot loci at pressures below the polymer transition region.
Due to its numerous practical applications and intriguing phase transformation behavior, shape memory alloys (SMAs) have garnered significant research and development interests. In the past, most studies on the mechanical behavior of SMAs have been conducted under uniaxial stress loadings. Limited research on SMAs under shock loading has not provided conclusive results regarding their transformation behavior and transformation stress under such loading. Additionally, there is a lack of comprehensive understanding regarding the effects of different stress states on transformation behavior. The main objectives of this study are to address these issues. To achieve these objectives, a series of shock wave experiments were designed and conducted. Additionally, quasi-static and dynamic uniaxial stress experiments were carried out to establish a baseline for comparison. The results revealed that the transformation stress under dynamic uniaxial strain shock loading was approximately 1.92 GPa in contrast to 0.5 GPa (quasi-static) to 0.8 GPa (dynamic) observed in uniaxial stress loading. The transformation behavior exhibited noticeable rate sensitivity for both types of loading. There appeared to be a critical strain rate above which the austenite phase was driven to a metastable state. This estimated critical axial strain rate along the loading direction was approximately 2 × 10 3 /s–4 × 10 3 /s for uniaxial stress loading and approximately 2 × 10 6 /s for uniaxial strain loading. The apparent high transformation stress for uniaxial strain loading can likely be attributed to a combination of high-pressure confinement and high strain rate. Furthermore, determining their relative contributions remains an open issue.