Search NASA⌕ Search

SEARCH · Search NASA

Results for “shocks”

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 109 records · Page 6

Evolution of the mechanical properties of a cobalt-based alloy under thermal shocks

A cobalt-base alloy manufactured by hot isostatic pressing sintering (HIP) was investigated in the present work. A series of quenching thermal-shock experiments were carried out to study the effects of thermal-cycle shocks on this new cobalt-based alloy. Scanning Electron Microscope (SEM) and Energy Dispersion Spectrum (EDS) were employed to observe the evolution of microstructures under thermal-shock cycles. Results show that some floccus Co 3 W 3 C precipitated around the WC particles after several thermal-shock cycles. These precipitates and the carbide particles correspond to the locations with the high-stress concentration. The mechanical properties evolution under thermal-shock cycles was analyzed by Nanoindentation tests. The nanohardness presents no noticeable change with thermal-shock cycles. However, the reduced modulus demonstrates a decreasing trend with the thermal-shock cycles. The variation of the mechanical properties has an unestimated relationship with the residual stress and densities of dislocations caused by thermal shocks. Furthermore, the present work pointed out that the precipitates around the carbides are the places where the thermal fatigue crack initiated. How to deal with the precipitates will be the optimized way for this material.

36 MATERIALS SCIENCE↗

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗

Shock-ramp analysis test problem

Quasi-isentropic (ramp) compression is now a well-established experimental method and so are the analysis techniques to give Lagrangian sound speed, pressure, and density along the sample material's isentrope. A shock followed by ramp compression is a natural extension to investigate, for example, shock melt and refreeze on compression, or isentropes of states off the Hugoniot or principal isentrope. In practice, graded-density impactors produce initial shocks, compression by shaped laser pulses may be unable to produce a smooth pressure increase from zero, and incidental perturbations on the drive pulse may also give rise to shocks, so robust shock-ramp analysis methods will be needed. Appropriate analysis methods are needed for shock-ramp experiments, based on those for quasi-isentropic compression, and these require validation. This paper describes three different analyses of a shock-ramp test problem, including an assessment of their estimated errors. The methods tested were based on hydrodynamic characteristics or integration backward in space. All methods gave the known Lagrangian sound speed to within ~1%, and pressure and volume to within less than 2% and 1%, demonstrating that the analysis methods of isentropic compression experiments can be confidently extended to the analysis of shock and ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal (rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multiphase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(\sigma _{\rho /\langle \rho \rangle }^{2\Gamma }-1)^{1/2}(1-\sigma _{\rho /\langle \rho \rangle }^{-2})^{-1/2}\mathcal {M}^{-1}\Gamma ^{-1/2}$, with the density dispersion σρ/<ρ>, the turbulent Mach number $\mathcal {M}$, and the polytropic exponent Γ. Purely solenoidal and purely compressive driving correspond to b ~ 1/3 and b ~ 1, respectively. As a result, using simulations in which a shock is driven into a multiphase medium with structures of different sizes and Γ < 1, we find b ~ 1 for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbulence generation by shock interaction with a highly nonuniform medium

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, $\sigma_\rho/\bar{\rho}\sim1$) and strong shocks (shock Mach numbers $\bar{M}_s\gtrsim10$). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Magnitude and Waveform of Shock Waves Induced by X-ray Lasers in Water

The high energy densities deposited in materials by focused X-ray laser pulses generate shock waves which travel away from the irradiated region, and can generate complex wave patterns or induce phase changes. We determined the time-pressure histories of shocks induced by X-ray laser pulses in liquid water microdrops, by measuring the surface velocity of the microdrops from images recorded during the reflection of the shock at the surface. Measurements were made with ~30 µm diameter droplets using 10 keV X-rays, for X-ray pulse energies that deposited linear energy densities from 3.5 to 120 mJ/m; measurements were also made with ~60 µm diameter drops for a narrower energy range. At a distance of 15 µm from the X-ray beam, the peak shock pressures ranged from 44 to 472 MPa, and the corresponding time-pressure histories of the shocks had a fast quasi-exponential decay with positive pressure durations estimated to range from 2 to 5 ns. Knowledge of the amplitude and waveform of the shock waves enables accurate modeling of shock propagation and experiment designs that either maximize or minimize the effect of shocks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron–Ion Temperature Ratio in Astrophysical Shocks

Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T e /T p inferred from the Hα profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T e /T p for supernova remnant shocks using Hα profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the Hα determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks, and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence, and varying preshock conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

Insight into the Chemistry of PETN Under Shock Compression Through Ultrafast Broadband Mid-Infrared Absorption Spectroscopy

In this study, thin films of pentaerythritol tetranitrate (PETN) were shock compressed using the laser driven shock apparatus at Los Alamos National Laboratory (LANL). Two spectroscopic probes were available to this apparatus: visible white light transient absorption spectroscopy (VIS) from 400 to 700 nm and mid-infrared transient absorption spectroscopy (MIR) from 1150 to 3800 cm –1 . Important PETN vibrational modes are the symmetric and antisymmetric NO 2 stretches at 1280 and 1650 cm –1 , respectively, as well as CH stretches at ~2900 cm –1 . Shock strength was varied from approximately 3 to 55 GPa to span from the chemically unreactive regime to the regime in which fast chemical reaction took place on the 250 ps time scale of the measurements. VIS and MIR results suggest irreversible chemistry was induced in PETN at pressures above 30 GPa. At lower shock pressures, the spectroscopy showed minimal changes attributable to pressure induced effects. Under the higher-pressure reactive conditions, the frequency region at the antisymmetric NO 2 stretch mode had a significantly increased absorption while the region around the symmetric NO 2 stretch did not. No observable increased absorption occurred in the higher frequency regions where CH-, NH-, and OH- bond absorptions would be observed. A broad absorption appeared on the shoulder at the red-edge of the CO 2 vibrational band around 2200 cm –1 . In addition to the experiments, reactive molecular dynamics were carried out under equivalent shock conditions to correlate the evolution of the infrared spectrum to molecular processes. The simulations show results consistent to experiments up to 30 GPa but suggest that NO and NO 2 related features provided the strongest contributions to the shocked infrared changes. Proposed mechanisms for shocked PETN chemistry are analyzed as consistent or inconsistent with the data presented here. Our experimental data suggests C≡O or N 2 O bond formation, nitrite formation, and absence of significant hydroxyl or amine concentrations in the initial chemistry steps in PETN shocked above 30 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock generated vorticity in spark discharges

Spark discharges induce a complex flow field consisting of a shock wave at early times (~ 1 μs), a pair of vortex rings, and a hot gas kernel. The vortex rings entrain ambient gas into the hot gas kernel and control its cooling and expansion. In this work, we investigate the shock wave's contribution in producing the vortex ring vorticity. We analyze high-speed (700 kHz) schlieren images of the shock wave for a range of electrical energies to measure the shock properties and estimate shock velocity and curvature. These measurements are combined with a model to calculate the vorticity generated. The measurements show that the highest vorticity is generated near the peak shock curvature location, and the shock curvature and strength increase with electrical energy deposited. A comparison of the vorticity estimated from the model to vorticity measurements from stereoscopic particle image velocimetry shows the results to be statistically equivalent. Furthermore, this suggests that the shock curvature and velocity contribute to the vortex rings induced by spark discharges.

42 ENGINEERING↗

Probabilistic Assessment of Damage from Same Shock Response Spectra due to Variations in Damping

Interpretation of field data from shock tests and subsequent assessment of product safety margins via laboratory testing are based on the shock response spectra (SRS). The SRS capture how a single degree of freedom (SDOF) structure responds to the shock at differing frequencies and, therefore, no longer contain the duration or other temporal parameters pertaining to the shock. A single duration can often be included in the technical specification or in the recreation of acceleration vs. time history from the specified SRS; however, there is little basis for that beyond technical judgment. The loss of such temporal information can result in the recreated SRS being the same while its effect on a system or component can be different. This paper attempts to quantify this deficiency as well as propose a simple method of capturing damping from shock waves that can allow the original waveform to be more accurately reconstructed from the SRS. In this study the decay rate associated with various frequencies that comprise the overall shock was varied. This variation in the decay rate leads to a variation in the acceleration vs. time history, which can be correlated to a “Damage Index” that captures the fatigue damage imparted to the object under shock. Several waveforms that have the same SRS but varying rates of decay for either high- or low-frequency components of the shock were investigated. The resulting variation in stress cycles and Damage Index is discussed in the context of the lognormal distribution of fatigue failure data. It is proposed that, along with the SRS, the decay rate is also captured to minimize the discrepancy between field data and representative laboratory tests.

36 MATERIALS SCIENCE↗

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data Based Shock Test Specification

It is vital that avionic packages used for testing and certifying the reliability and safety of U.S. nuclear weapons with platform aircraft survive exposure to shock environments during transportation and delivery. The objective of this research was to characterize the response to these transportation shock environments delivering accurate shock test specifications in order to set laboratory programming material and device certification rigor. Responses to shock events were analyzed in the frequency domain via the Shock Response Spectrum (SRS). Shocks were then grouped based on respective behavior of maximum response accelerations which were pseudorandomly resampled and compared to test data to form test specifications based on the MinerPalmgren hypothesis. In addition to discovering significant over testing in current shock specifications, a new systematic, data-driven approach to designing shock specifications was formulated.

42 ENGINEERING↗

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗

Insight into the chemistry of TNT during shock compression through ultrafast absorption spectroscopies

Thin films of trinitrotoluene (TNT) were shock compressed using the ultrafast laser shock apparatus at Los Alamos National Laboratory. Visible (VIS) and mid-infrared (MIR) transient absorption spectroscopies were simultaneously performed to probe for electronic and vibrational changes during shock compression of TNT. Additionally, three shock pressures (16 GPa, 33 GPa, and 45 GPa) were selected to observe no reaction, incipient reaction, and strongly developed reactions for TNT within the experimental time scale of <250 ps. Negligible absorption changes in MIR or VIS absorptions were observed at 16 GPa. At 33 GPa, MIR absorptions in the 3000 cm -1 –4000 cm -1 range were observed to increase during the shock and continue to increase during the rarefaction, in contrast to the VIS absorption measurements, which increased during the shock and almost fully recovered during rarefaction. At 45 GPa, both VIS and MIR absorptions were strong and irreversible. The intense and spectrally broad MIR absorptions were attributed to short lived intermediates with strong, spectrally broad absorptions that dominate the spectral response. The MIR and VIS absorption changes observed at 33 GPa and 45 GPa were credited to shock induced chemistry, most likely including the formation of a very broad hydrogenic stretch feature. The results from these experiments are consistent with the chemical mechanisms that include O–H or N–H formation such as CH 3 oxidation or C–N homolysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗