Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Magnetohydrodynamic effects in a shock-accelerated gas cylinder

This work presents two-dimensional (2D) simulations on a cylindrical Richtmyer-Meshkov instability (RMI) in magnetohydrodynamics (MHD). Three studies are presented in an effort to quantify and qualify the evolution of the MHD RMI by varying the magnetic field orientation, strength of the magnetic field, and strength of the shock wave driving the instability. The orientations considered herein are either parallel or perpendicular to the shock wave motion. The second study varies the magnetic fields between 100, 250, and 500 gauss (G), while the third study considers incident shock wave Mach numbers M = 1.2, M = 1.66, and M = 2.2. These parameter ranges were selected to be easily achievable in experiments while the interface perturbation was selected such that its evolution is independent of either the shock wave or magnetic field orientations independently. It was found that the MHD RMI evolution is dependent upon the magnetic field orientation relative to the shock transit direction as well as their individual magnitudes. This is because the mechanism of suppression, attributed to Alfven waves, is a function of the magnetic field strength and the orientation of the magnetic field, while the mechanism of RMI evolution, baroclinic vorticity deposition, is a function of the Mach number. Stronger magnetic fields were found to provide greater mixing suppression and have significant effects on RMI-like interface morphology. Lastly, increasing the shock wave strength generated competing effects between higher RMI vorticity deposition and greater vorticity removal from the interface by faster Alfven waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Counterpropagating Radiative Shock Experiments on the Orion Laser

We present new experiments to study the formation of radiative shocks and the interaction between two counterpropagating radiative shocks. The experiments are performed at the Orion laser facility, which is used to drive shocks in xenon inside large aspect ratio gas cells. The collision between the two shocks and their respective radiative precursors, combined with the formation of inherently three-dimensional shocks, provides a novel platform particularly suited for the benchmarking of numerical codes. The dynamics of the shocks before and after the collision are investigated using point-projection x-ray backlighting while, simultaneously, the electron density in the radiative precursor was measured via optical laser interferometry. Modeling of the experiments using the 2D radiation hydrodynamic codes nym and petra shows very good agreement with the experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Laser pulse-length dependent ablation and shock generation in silicon at 5 × 10 14 W/ cm 2 intensities

The effect of laser pulse duration on energy coupling into a planar silicon target is investigated in experiments at the OMEGA-EP facility by varying the laser pulse length τ —spanning 3 orders of magnitude from 100 ps to 10 ns—while maintaining a constant peak laser intensity, I 0 = 5 × 10 14 W/ cm 2 . In theoretical models, the ablation pressure primarily scales for a given material with laser intensity and wavelength, which are all fixed variables here, allowing us to explore the specific role of laser pulse duration. Two-dimensional radiation-hydrodynamics simulations benchmarked with optical probing of the expanding plasma show that the pulse duration is critical for the ablation pressure to reach a steady state. Moreover, the pulse duration impacts shock decay and multiple wave effects, which strongly dictate the evolving shock profile that propagates within the laser-shocked target as ultimately measured by rear-surface diagnostics. The shock velocities inferred from the theoretical model, after considering shock decay, impedance matching, and shock Hugoniot, are found to be in good agreement with velocimetry measurements. However, discrepancies are observed with simulations for the shorter (0.1 ns) and longer (10 ns) pulse durations, which are respectively attributed to unaccounted contributions of kinetic absorption mechanisms and instabilities in simulations. Published by the American Physical Society 2024

Bailly-Grandvaux, M. (ORCID:0000000175294013)↗
DOE OSTI2019

Dalton’s and Amagat’s laws fail in gas mixtures with shock propagation

A shock propagating through a gas mixture leads to pressure, temperature, and density increases across the shock front. Rankine-Hugoniot relations correlating pre- and post-shock quantities describe a calorically perfect gas but deliver a good approximation for real gases, provided the pre-shock conditions are well characterized with a thermodynamic mixing model. Two classic thermodynamic models of gas mixtures are Dalton’s law of partial pressures and Amagat’s law of partial volumes. We measure post-shock temperature and pressure in experiments with nonreacting binary mixtures of sulfur hexafluoride and helium (two dramatically disparate gases) and show that neither model can accurately predict the observed values, on time scales much longer than that of the shock front passage, due to the models’ implicit assumptions about mixture behavior on the molecular level. However, kinetic molecular theory can help account for the discrepancy. Our results provide starting points for future theoretical work, experiments, and code validation.

42 ENGINEERING↗
DOE OSTI2023

Real-time (nanoseconds) determination of liquid phase growth during shock-induced melting

Melting of solids is a fundamental natural phenomenon whose pressure dependence has been of interest for nearly a century. However, the temporal evolution of the molten phase under pressure has eluded measurements because of experimental challenges. By using the shock front as a fiducial, we investigated the time-dependent growth of the molten phase in shock-compressed germanium. In situ x-ray diffraction measurements at different times (1 to 6 nanoseconds) behind the shock front quantified the real-time growth of the liquid phase at several peak stresses. These results show that the characteristic time for melting in shock-compressed germanium decreases from ~7.2 nanoseconds at 35 gigapascals to less than 1 nanosecond at 42 gigapascals. Our melting kinetics results suggest the need to consider heterogeneous nucleation as a mechanism for shock-induced melting and provide an approach to measuring melting kinetics in shock-compressed solids.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Technique for Forming Solid D2 and D-T Layers for Shock Timing Experiments at the National Ignition Facility

Capsule implosion experiments on the National Ignition Facility (NIF) are driven with a carefully tailored laser pulse that delivers a sequence of shocks to the ablator and fuel. In order to ensure the shocks converge at the desired position, the shock strength and velocity are measured in experimental platforms referred to as keyhole targets. We made shock measurements on capsules completely filled with liquid deuterium for the solid deuterium tritide (D-T) layer campaigns. Modeling has been used to extend these results to form an estimate of the shock properties in solid D-T layers. Furthermore, to verify and improve the surrogacy of the liquid-filled keyhole measurements, we have developed a technique to form a solid layer inside the keyhole capsule. The layer is typically uniform over a 400-μm-diameter area. This is sufficient to allow direct measurement of the shock velocity. This layering technique has been successfully applied to 13 experiments on the NIF. The technique may also be applicable to fast-igniter experiments since some proposed designs resemble keyhole targets. We discuss our method in detail and give representative results.

42 ENGINEERING↗
DOE OSTI2022

Quantum-cascade-laser-absorption-spectroscopy diagnostic for temperature, pressure, and NO X 2 Π 1/2 at 500 kHz in shock-heated air at elevated pressures

The design, validation, and application of a quantum-cascade-laser-absorption-spectroscopy diagnostic for measuring gas temperature, pressure, and nitric oxide (NO) in high-temperature air are presented. A distributed-feedback quantum-cascade laser (QCL) centered near 1976 c m − <#comment/> 1 was used to scan across two transitions of NO in its ground electronic state ( X 2 Π <#comment/> 1 / 2 ). A measurement rate of 500 kHz was achieved using a single QCL by: (1) performing current modulation through a bias-tee, and (2) targeting closely spaced transitions with a large difference in lower-state energy. The diagnostic was validated in a mixture of 95% argon and 5% NO, which was shock-heated to ≈ <#comment/> 2000 to 3700 K. The average mean percent differences between laser-absorption-spectroscopy (LAS) measurements and predictions from shock-jump relations for temperature, pressure, and NO mole fraction were 3.1%, 4.1%, and 6.5%, respectively. The diagnostic was then applied to characterize shock-heated air at high temperatures (up to ≈ <#comment/> 5500 K ) and high pressures (up to 12 atm) behind either incident or reflected shocks. The LAS measurements were compared to theoretical predictions from shock-jump relations, pressure sensors mounted in the wall of the shock tube, and equilibrium values of the NO mole fraction. The average mean percent differences between LAS measurements and their aforementioned reference values were 3.2%, 10.8%, and 10.4% for temperature, pressure, and NO mole fraction, respectively. Last, a comparison between a measured NO mole fraction time history and a time-stepped homogeneous reactor simulation performed using two different chemical kinetics mechanisms is presented.

Gilvey, Jonathan J. (ORCID:0000000317664383)↗
DOE OSTI2017

Transition of dislocation glide to shear transformation in shocked tantalum

A TEM study of pure tantalum and tantalum-tungsten alloys explosively shocked at a peak pressure of 30 GPa (strain rate: ~1 x 10 4 sec -1 ) is presented. While no ω (hexagonal) phase was found in shock-recovered pure Ta and Ta-5W that contain mainly a low-energy cellular dislocation structure, shock-induced ω phase was found to form in Ta-10W that contains evenly distributed dislocations with a stored dislocation density higher than 1 x 10 12 cm -2 . The TEM results clearly reveal that shock-induced α (bcc) → ω (hexagonal) shear transformation occurs when dynamic recovery reactions which lead the formation low-energy cellular dislocation structure become largely suppressed in Ta-10W shocked under dynamic (i.e., high strain-rate and high-pressure) conditions. A novel dislocation-based mechanism is proposed to rationalize the transition of dislocation glide to twinning and/or shear transformation in shock-deformed tantalum. Lastly, twinning and/or shear transformation take place as an alternative deformation mechanism to accommodate high-strain-rate straining when the shear stress required for dislocation multiplication exceeds the threshold shear stresses for twinning and/or shear transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of Pressure Gradients on the Initiation of PBX-9502 via Irregular (Mach) Reflection of Low Pressure Curved Shock Waves

In the instance of multiple fragment impact on cased explosive, isolated curved shocks are generated in the explosive. These curved shocks propagate and may interact and form irregular or Mach reflections along the interaction loci, thereby producing a single shock that may be sufficient to initiate PBX-9501. However, the incident shocks are divergent and their intensity generally decreases as they expand, and the regions behind the Mach stem interaction loci are generally unsupported and allow release waves to rapidly affect the flow. The effects of release waves and divergent shocks may be considered theoretically through a “Shock Change Equation”.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

Experimental and Theoretical Investigation of Shock-Induced Reactions in Energetic Materials

In this work, shock-induced reactions in high explosives and their chemical mechanisms were investigated using state-of-the-art experimental and theoretical techniques. Experimentally, ultrafast shock interrogation (USI, an ultrafast interferometry technique) and ultrafast absorption spectroscopy were used to interrogate shock compression and initiation of reaction on the picosecond timescale. The experiments yielded important new data that appear to indicate reaction of high explosives on the timescale of tens of picoseconds in response to shock compression, potentially setting new upper limits on the timescale of reaction. Theoretically, chemical mechanisms of shock-induced reactions were investigated using density functional theory. The calculations generated important insights regarding the ability of several hypothesized mechanisms to account for shock-induced reactions in explosive materials. The results of this work constitute significant advances in our understanding of the fundamental chemical reaction mechanisms that control explosive sensitivity and initiation of detonation.

36 MATERIALS SCIENCE↗
DOE OSTI2017

On the Link between the Release of Solar Energetic Particles Measured at Widespread Heliolongitudes and the Properties of the Associated Coronal Shocks

Under the paradigm that the main agents in the acceleration of solar energetic particles (SEPs) are shocks initially driven by coronal mass ejections, we analyze whether the properties of the shocks in the corona inferred from combining extreme-ultraviolet (EUV) and white-light (WL) observations from multiple vantage points together with magnetohydrodynamic (MHD) simulations of the corona can be used to determine the release of SEPs into different regions of the heliosphere and hence determine the longitudinal extent of the SEP events. We analyze the SEP events observed on 2011 November 3, 2013 April 11, and 2014 February 25 over a wide range of heliolongitudes. MHD simulations provide the characteristics of the background medium where shocks propagate, in particular the Alfvén and sound speed profiles that allow us to determine both the extent of the EUV waves in the low corona and the fast magnetosonic Mach number ( M {sub FM}) of the shocks. The extent of the EUV waves in the low corona is controlled by this background medium and does not coincide with the extent of the SEP events in the heliosphere. Within the uncertainties of (i) the extent and speed of the shock inferred from EUV and WL images and (ii) the assumptions made in the MHD models, we follow the evolution of M {sub FM} at the region of the shock magnetically connected to each spacecraft. The estimated release times of the first SEPs measured by each spacecraft does not coincide with the time when the M {sub FM} at this region exceeds a given threshold.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Evidence for a Merger-induced Shock Wave in ZwCl 0008.8+5215 with Chandra and Suzaku

We present the results from new deep Chandra (~410 ks) and Suzaku (~180 ks) observations of the merging galaxy cluster ZwCl 0008.8+5215 (z = 0.104). Previous radio observations revealed the presence of a double radio relic located diametrically west and east of the cluster center. Using our new Chandra data, we find evidence for the presence of a shock at the location of the western relic, RW, with a Mach number ${{ \mathcal M }}_{{S}_{X}}={1.48}_{-0.32}^{+0.50}$ from the density jump. We also measure ${{ \mathcal M }}_{{T}_{X}}={2.35}_{-0.55}^{+0.74}$ and ${{ \mathcal M }}_{{T}_{X}}={2.02}_{-0.47}^{+0.74}$ from the temperature jump, with Chandra and Suzaku, respectively. These values are consistent with the Mach number estimate from a previous study of the radio spectral index, under the assumption of diffusive shock acceleration (${{ \mathcal M }}_{\mathrm{RW}}={2.4}_{-0.2}^{+0.4}$). Interestingly, the western radio relic does not entirely trace the X-ray shock. A possible explanation is that the relic traces fossil plasma from nearby radio galaxies that is reaccelerated at the shock. For the eastern relic we do not detect an X-ray surface brightness discontinuity, despite the fact that radio observations suggest a shock with ${{ \mathcal M }}_{\mathrm{RE}}={2.2}_{-0.1}^{+0.2}$. The low surface brightness and reduced integration time for this region might have prevented the detection. The Chandra surface brightness profile suggests ${ \mathcal M }\lesssim 1.5$, while the Suzaku temperature measurements found ${{ \mathcal M }}_{{T}_{X}}={1.54}_{-0.47}^{+0.65}$. Lastly, we also detect a merger-induced cold front on the western side of the cluster, behind the shock that traces the western relic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Electron Heating in the Transrelativistic Perpendicular Shocks of Tilted Accretion Flows

Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount T e 2 / T e 2 , ad − 1 ≃ 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1 / T i1 , which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1 / T i1 = 1, whereas heating by B -parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1 / T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.

Astronomy & Astrophysics↗
DOE OSTI2026

The Role of Three-dimensional Effects on Ion Injection and Acceleration in Perpendicular Shocks

Understanding the conditions that enable particle acceleration at nonrelativistic collisionless shocks is essential to unveil the origin of cosmic rays. We employ 2D and 3D hybrid simulations (with kinetic ions and fluid electrons) to explore particle acceleration and magnetic field amplification in nonrelativistic perpendicular shocks, focusing on the role of shock drift acceleration and its dependence on the shock Mach number. We perform an analysis of the ion injection process and demonstrate why efficient acceleration is only observed in 3D. In particular, we show that ion injection critically depends on the “porosity” of the magnetic turbulence in the downstream region near the shock, a property describing how easily the postshock region allows particles to traverse it and return upstream without being trapped. This effect can only be properly captured in 3D. Additionally, we explore the impact of numerical resolution on ion energization, highlighting how resolving small-scale turbulence—on scales below the thermal ion gyroradius—is essential for accurately modeling particle injection. Overall, our results emphasize the necessity of high-resolution 3D simulations to capture the fundamental microphysics driving particle acceleration at perpendicular shocks.

Orusa, Luca [Princeton Univ., NJ (United States); ↗
DOE OSTI2017

Effects of Preionization in Radiative Shocks. I. Self-consistent Models

In this paper we treat the preionization problem in shocks over the velocity range 10 < v {sub s} < 1500 km s{sup −1} in a self-consistent manner. We identify four distinct classes of solutions controlled by the value of the shock-precursor parameter, Ψ=Q/v{sub s}, where Q is the ionization parameter of the UV photons escaping upstream. This parameter determines both the temperature and the degree of ionization of the gas entering the shock. In increasing velocity, the shock solution regimes are cold neutral precursors (v {sub s} ≲ 40 km s{sup −1}), warm neutral precursors (40 ≲ v {sub s} ≲ 75 km s{sup −1}), warm partly ionized precursors (75 ≲ v {sub s} ≲ 120 km s{sup −1}), and fast shocks in which the preshock gas is in photoionization equilibrium and is fully ionized. The main effect of a magnetic field is to push these velocity ranges to higher values and to limit the postshock compression. In order to facilitate comparison with observations of shocks, we provide a number of convenient scaling relationships for parameters, such as postshock temperature, compression factors, cooling lengths, and Hβ and X-ray luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

Abstract In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ALMA-SZ detection of a galaxy cluster merger shock at half the age of the universe

We present ALMA measurements of a merger shock using the thermal Sunyaev–Zel’dovich (SZ) effect signal, at the location of a radio relic in the famous El Gordo galaxy cluster at z≈0.9. Multi-wavelength analysis in combination with the archival Chandra data and a high-resolution radio image provides a consistent picture of the thermal and non-thermal signal variation across the shock front and helps to put robust constraints on the shock Mach number as well as the relic magnetic field. We employ a Bayesian analysis technique for modeling the SZ and X-ray data self-consistently, illustrating respective parameter degeneracies. Combined results indicate a shock with Mach number M=2.4{sub −0.6}{sup +1.3}, which in turn suggests a high value of the magnetic field (of the order of 4–10 μG) to account for the observed relic width at 2 GHz. At roughly half the current age of the universe, this is the highest-redshift direct detection of a cluster shock to date, and one of the first instances of an ALMA-SZ observation in a galaxy cluster. It shows the tremendous potential for future ALMA-SZ observations to detect merger shocks and other cluster substructures out to the highest redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counter- propagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scat- tering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. We compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗