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,513 records · Page 84

The Efficiency of Magnetic Field Amplification at Shocks by Turbulence

Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks (approx.100 micrograms) and in radio relics in galaxy clusters (approx. micrograms). We present high-resolution magnetohydrodynamic simulations of the interaction between pre-shock turbulence, clumping and shocks, to quantify the conditions under which turbulent dynamo amplification can be significant. We demonstrate numerically converged field amplification which scales with Alfven Mach number, B/B0 varies as MA, up to MA approx.150.This implies that the post-shock field strength is relatively independent of the seed field. Amplification is dominated by compression at low MA, and stretching (turbulent amplification) at high MA. For high MA, the B-field grows exponentially and saturates at equipartition with turbulence, while the vorticity jumps sharply at the shock and subsequently decays; the resulting field is orientated predominately along the shock normal (an effect only apparent in 3D and not 2D). This agrees with the radial field bias seen in supernova remnants. By contrast, for low MA, field amplification is mostly compressional, relatively modest, and results in a predominantly perpendicular field. The latter is consistent with the polarization seen in radio relics. Our results are relatively robust to the assumed level of gas clumping. Our results imply that the turbulent dynamo may be important for supernovae, but is only consistent with the field strength, and not geometry, for cluster radio relics. For the latter, this implies strong pre-existing B-fields in the ambient cluster outskirts.

lementary particles – magnetic fields – radiat↗

Comparison of the Coronal Mass Ejection Shock Acceleration of Three Widespread SEP Events During Solar Cycle 24

We studied three solar energetic particle (SEP) events observed on 14 August 2010, 3 November 2011, and 5 March 2013 by Solar Terrestrial Relations Observatory (STEREO) A, B, and near-Earth (L1) spacecraft with a longitudinal distribution of particles greater than 90 degrees. Using a forward modeling method combined with extreme ultraviolet and white-light images, we determined the angular extent of the shock, the time and location (cobpoint) of the shock intersection with the magnetic field line connecting to each spacecraft, and compute the shock speed at the cobpoint of each spacecraft. We then examine whether the observations of SEPs at each spacecraft were accelerated and injected by the spatially extended shocks or whether another mechanism such as cross-field transport is required for an alternative explanation. Our analyses results indicate that the SEPs observed at the three spacecraft on 3 November, STEREO B (STB) and L1 on 14 August, and the 5 March SEP event at STEREO A (STA) can be explained by the direct shock acceleration. This is consistent with the observed significant anisotropies, short time delays between particle release times and magnetic connection times, and sharp rises in the SEP time profiles. Cross-field diffusion is the likely cause for the 14 August SEP event observed by STA and the 5 March SEPs observed by STB and L1 spacecraft, as particle observations featured weak electron anisotropies and slow rising intensity profiles. Otherwise, the wide longitudinal spread of these SEP increases would require an existence of a circumsolar shock, which may not be a correct assumption in the corona and heliosphere.

Xie, H.↗

Delayed Shock-Induced Dust Formation in the Dense Circumstellar Shell Surrounding the Type IIn Supernova SN 2010jl

The light curves of Type IIn supernovae are dominated by the radiative energy released through the interaction of the supernova shock waves with their dense circumstellar medium (CSM). The ultraluminous Type IIn supernova SN2010jl exhibits an infrared emission component that is in excess of the extrapolated UV-optical spectrum as early as few weeks postexplosion. This emission has been considered by some as evidence for the rapid formation of dust in the cooling postshock CSM. We investigate the physical processes that may inhibit or facilitate the formation of dust in the CSM. When only radiative cooling is considered, the temperature of the dense shocked gas rapidly drops below the dust condensation temperature. However, by accounting for the heating of the postshock gas by the downstream radiation from the shock, we show that dust formation is inhibited until the radiation from the shock weakens as it propagates into the less dense outer regions of the CSM. In SN2010jl, dust formation can therefore only commence after day ∼380. Only the IR emission since that epoch can be attributed to the newly formed CSM dust. Observations on day 460 and later show that the IR luminosity exceeds the UV-optical luminosity. The postshock dust cannot extinct the radiation emitted by the expanding SN shock. Therefore, its IR emission must be powered by an interior source, which we identify as the reverse shock propagating through the SN ejecta. IR emission before day 380 must therefore be an IR echo from preexisting CSM dust.

Sarangi, Arkaprabha↗

MMS Observation of Inverse Energy Dispersion in Shock Drift Accelerated Ions

The four Magnetospheric Multiscale (MMS) spacecraft observed a ∼1 min burst of energetic ions (50–1000 keV) in the region upstream from the subsolar quasi-perpendicular bow shock on 6 December 2015. The composition, flux levels, and spectral indices of these energetic protons, helium, and oxygen ions greatly resemble those seen in the outer magnetosphere earlier while MMS crossed the magnetopause and differ significantly from those simultaneously observed far upstream by Advanced Composition Explorer (ACE). However, the event cannot be explained solely in terms of leakage from the magnetosphere. The strongly southward orientation of the interplanetary magnetic field (IMF) lines at the time of the event precludes any connection to the magnetosphere. This point is confirmed by the presence of energetic electrons, known to occur on magnetic field lines that graze the bow shock rather than connect to the magnetosphere. We suggest that the ions gradient drifted out of the nearby quasi-parallel foreshock and into the quasi-perpendicular bow shock. Each of the ion species exhibited an inverse energy dispersion. As predicted by models for shock drift acceleration, the energies of the ions increased as 𝜃(sub Bn), the angle between the IMF and the shock normal, increased. Finally, we note that a similar event was observed a few minutes later in the subsolar magnetosheath, indicating that such events can be swept downstream of the bow shock.

Lee, S. H.↗

Electron Scattering by Low-frequency Whistler Waves at Earth’s Bow Shock

Electrons are accelerated to nonthermal energies at shocks in space and astrophysical environments. While shock drift acceleration (SDA) has been considered a key process of electron acceleration at Earth’s bow shock, it has also been recognized that SDA needs to be combined with an additional stochastic process to explain the observed power-law energy spectra. Here, we show mildly energetic (∼0.5 keV) electrons are locally scattered (and accelerated while being confined) by magnetosonic-whistler waves within the shock transition layer, especially when the shock angle is large (θ(sub Bn) approximately equal or greater than 70°). When measured by the Magnetospheric Multiscale mission at a high cadence, ∼0.5 keV electron flux increased exponentially in the shock transition layer. However, the flux profile was not entirely smooth and the fluctuation showed temporal/spectral association with large-amplitude (δB/B ~ 0.3), low-frequency (approximately equal or less than 0.1 Ω(sub ce) where Ω(sub ce) is the cyclotron frequency), obliquely propagating (θ(sub kB) ~ 30°–60°, where θ(sub kB) is the angle between the wave vector and background magnetic field) whistler waves, indicating that the particles were interacting with the waves. Particle simulations demonstrate that, although linear cyclotron resonances with ∼0.5 keV electrons are unlikely due to the obliquity and low frequencies of the waves, the electrons are still scattered beyond 90° pitch angle by (1) resonant mirroring (transit-time damping), (2) non-resonant mirroring, and (3) subharmonic cyclotron resonances. Such coupled nonlinear scattering processes are likely to provide the stochasticity needed to explain the power-law formation.

Oka, M.↗

Energy Partition at Collisionless Supercritical Quasiperpendicular Shocks

Collisionless shocks in astrophysical plasmas are important thermalizers, converting some of the incident flow energy into thermal energy, and non-thermalizers, partitioning that energy in unequal ways to different particle species, sub-populations thereof, and field components. This partition problem, or equivalently the shock equation of state, lies at the heart of shock physics. Here we employ systematically a framework to capture all the incident and downstream energy fluxes at two example traversals of the Earth’s bow shock by the Magnetospheric Multiscale Mission. Here and traditionally such data has to be augmented by information from other spacecraft, e.g., to provide more accurate measurements of the cold solar wind beam. With some care and fortuitous choices, the energy fluxes are constant, including instantaneous measurements through the shock layer. The dominant incident proton ram energy is converted primarily into downstream proton enthalpy flux, the majority of which is actually carried by a small fraction of suprathermal protons. Fluctuations include both real and instrumental effects. Separating these, resolving the solar wind beam, and other considerations point the way to a dedicated mission to solve this energy partition problem across a full range of plasma and shock conditions.

Steven J Schwartz↗

Examination of Mars2020 Shock-Layer Conditions via Infrared Laser Absorption Spectroscopy of CO2 and CO

Amid-infrared laser absorption diagnostic was deployed to study a simulated Mars2020 shock layer in the Electric Arc Shock Tube (EAST) facility at NASA’s Ames Research Center. Rapid RF-diplexing techniques enabled quantitative temperature and number density measurements of CO2 and CO with μs-resolution over an incident shock velocity range of 1.39 – 3.75 km/s. Two interband cascade lasers were utilized at 4.17 and 4.19 μm to resolve eight CO2 asymmetric stretch fundamental band (𝜈3) transitions from two different vibrational levels: 0000 (ground state) and 0110 (first excited bending mode). The probed rotational levels span across J” = 58 to J” = 140. Results are compared to DPLR simulations of the shock layer using kinetic mechanisms of Johnston et al. and Cruden et al. At shock velocities below 3.1 km/s, the agreement between the measurements and the Johnston mechanism is typically within 5% for temperature and within 10% for number density. At shock velocities above 3.1 km/s, the CO2 measurement becomes sensitive to a thin boundary layer and corrections of this effect are presented. On test cases with enough energy to dissociate CO2, a quantum cascade laser scanned the P(2, 20), P(0, 31), and P(3, 14) transitions of the CO fundamental band at 4.98 μm. CO formation rate is measured to be close to the Johnston kinetic mechanism at low velocities, and then trending towards the Cruden kinetic mechanism at high velocities. On a few low velocity test cases, rovibrational relaxation of the Martian atmosphere is probed with μs resolution.

Laser Absorption↗

Measurement and Analysis of the Extreme Physical Shock Environment Experienced by Crane-Mounted Radiation Detection Systems.

Potentially, radiation detectors at ports of entry could be mounted on container gantry crane spreaders to monitor cargo containers entering and leaving the country. These detectors would have to withstand the extreme physical environment experienced by these spreaders during normal operations. Physical shock data from the gable ends of a spreader were recorded during the loading and unloading of a cargo ship with two Lansmont SAVER 9X30 units (with padding) and two PCB Piezotronics model 340A50 accelerometers (hard mounted). Physical shocks in the form of rapid acceleration were observed in all accelerometer units with values ranging from 0.20 g’s to 199.99 g’s. The majority of the shocks for all the Lansmont and PCB accelerometers were below 50 g’s. The Lansmont recorded mean shocks of 21.83 ± 13.62 g’s and 24.78 ± 11.49 g’s while the PCB accelerometers experienced mean shocks of 34.39 ± 25.51 g’s and 41.77 ± 22.68 g’s for the landside and waterside units, respectively. Encased detector units with external padding should be designed to withstand at least 200 g’s of acceleration without padding and typical shocks of 30 g’s with padding for mounting on a spreader.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Systems and methods for shock-resistant memory devices

A shock-resistant memory device comprises a housing and a memory module. The memory module is disposed within the housing and surrounded by potting material to protect the memory module from damage during a shock event. The housing can include a port that accommodates a data connection between the memory module and a sensor from which data is desirably received by the memory module. During a shock event the connection between the memory module and the sensor may be severed, but data stored in the memory module can be retained in the memory module which is protected by the housing. To facilitate retrieval of the memory device subsequent to a shock event, a balloon can be affixed to the housing. The balloon can be configured to inflate subsequent to the shock event so that the shock-resistant memory device does not sink in water and to make the memory device more visible for recovery.

42 ENGINEERING↗

Drop Hammer Tests with Three Oleo Strut Models and Three Different Shock Strut Oils at Low Temperatures

Drop hammer tests with different shock strut models and shock strut oils were performed at temperatures ranging to -40 C. The various shock strut models do not differ essentially regarding their springing and damping properties at low temperatures; however, the influence of the different shock strut oils on the springing properties at low temperatures varies greatly.

LUBRICATING OILS, MINERAL↗

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↗

Shock initiation of the HMX-based explosive PBX 9012: Experiments, uncertainty analysis, and unreacted equation-of-state

In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.

36 MATERIALS SCIENCE↗

Shock initiation and hot spots in plastic-bonded 1,3,5-triamino-2,4,6-trinitrobenzene (TATB)

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock recompaction of spall damage

Spall fracture is a high strain-rate damage phenomenon associated with shock or impulsive loading events. When a material that has been subjected to shock compression is allowed to release, rarefaction waves propagate into the sample and reduce the internal stress to zero. If multiple rarefaction waves intersect, they generate tension which, if sufficient, can nucleate voids in the material. It has been observed in several works investigating spall fracture that although the shock-wave profiles suggested spall occurred, imaging of the recovered sample revealed no voids or cracks. In this study, we aim to determine whether a second shock event could recompact existing spall damage, and if so, what form does the microstructure at the recompaction interface have? Through a series of gas-gun flyer-plate impact experiments, we demonstrate that modest shock stresses of 2 GPa–3 GPa are enough to both fully compact a damaged copper target back to a state of zero porosity and, furthermore, drive recrystallization of the interface such that there is a new bond formed where the free surfaces were brought together.

42 ENGINEERING↗

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Micro-strains, local stresses, and coherently diffracting domain size in shock compressed Al(100) single crystals

In situ x-ray diffraction (XRD) measurements and their analysis in Al single crystals shock compressed along the [100]-direction were utilized to examine shock wave induced microstructural heterogeneities. High-resolution XRD line profiles for the 200, 400, and 600 Al peaks were measured in uniaxial strain compression states to either 5.6 or 11.7 GPa and partial stress release to 3.5 or 6.6 GPa, respectively. Broadening of the XRD line profiles was analyzed to determine the magnitude of the longitudinal micro-strain distribution (0.195% and 0.28% full width at half maximum for 3.5 and 6.6 GPa stresses, respectively) and the size of coherently diffracting domains (CDDs) (0.125 and 0.07 μm for 3.5 and 6.6 GPa stresses, respectively). From the longitudinal micro-strain distributions, the distribution of local stress differences (or stress deviators) was obtained in the shocked state. The full width at half maximum of this distribution, a measure of the local stress inhomogeneities, is greater than half of the macroscopic stress difference for both 3.5 and 6.6 GPa peak stresses, suggesting considerable variation in local stress deviators. The CDD sizes determined here are comparable to characteristic length scales for defect-free regions determined from defect density measurements in post-shock recovery experiments. As a result, the present work represents an important step in understanding material microstructure and inhomogeneities in the shock-compressed state.

36 MATERIALS SCIENCE↗

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

Shock-induced chemical decomposition and overdriven detonation in hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) single crystals

Understanding the differences in the shock compression and detonation response of insensitive high explosives (IHEs) and conventional HEs (CHEs) is a long-standing need in HE science and technology. Having previously examined 1,1-diamino-2,2-dinitroethene (FOX-7) IHE single crystals, the shock and detonation response of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX)—a widely used CHE—was determined using wave profile measurements in ~250 μm thick single crystals shock compressed to 63 GPa. In marked contrast to FOX-7, RDX single crystals shocked along the [100] and [111] orientations showed wave profile features consistent with chemical decomposition onset at 15 GPa. These features were more pronounced for [100] RDX, suggesting a higher decomposition rate compared to [111] RDX. At 51 GPa and above, flat-topped single waves were observed for both orientations, showing the classic Chapman–Jouget (C–J) detonation response in which the decomposition is completed within the detonation front. The Hugoniot states and sound speeds determined for the detonation products were similar for both [100] and [111] orientations, showing that the overdriven detonation response for RDX does not depend on crystal orientation. The C–J pressure for RDX single crystals (35 GPa)—determined experimentally—is comparable to that of FOX-7. However, compared to FOX-7, chemical decomposition onset for RDX occurs at much lower pressures and the overdriven C–J detonation response occurs at higher pressures. Lastly, the present findings constitute the first experimental comparison of the shock and detonation response of conventional and insensitive HE single crystals over a broad pressure range below and above the C–J pressure.

36 MATERIALS SCIENCE↗