Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nuclear Weapon Blast Effects

Presentation Overview: Nuclear Weapon Blast Effects; Conventional Explosions; Energy Distribution; Thermal Effects; Fire Storm Effects; Blast and Shock Effects; Types of Bursts

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

Microcomputed X-Ray Tomographic Imaging and Image Processing for Microstructural Characterization of Explosives

Microstructural characterization of composite high explosives (HEs) has become increasingly important over the last several decades in association with the development of high fidelity mesoscale modeling and an improved understanding of ignition and detonation processes. HE microstructure influences not only typical material properties (e.g., thermal, mechanical) but also reactive behavior (e.g., shock sensitivity, detonation wave shape). A detailed nondestructive 3D examination of the microstructure has generally been limited to custom-engineered samples or surrogates due to poor contrast between the composite constituents. Highly loaded (>90 wt%) HE composites such as plastic-bonded explosives (PBX) are especially difficult. Here, we present efforts to improve measurement quality by using single and dual-energy microcomputed X-ray tomography and state-of-the-art image processing techniques to study a broad set of HE materials. Some materials, such as PBX 9502, exhibit suitable contrast and resolution for an automatic segmentation of the HE from the polymer binder and the voids. Other composite HEs had varying levels of success in segmentation. Post-processing techniques that used commercially available algorithms to improve the segmentation quality of PBX 9501 as well as zero-density defects such as cracks and voids could be easily segmented for all samples. Aspects of the materials that lend themselves well to this type of measurement are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Radiation Transport Two-temperature GRMHD Simulations of Warped Accretion Disks

In many black hole (BH) systems, the accretion disk is expected to be misaligned with respect to the BH spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the BH can tear the disk into multiple, independently precessing "sub-disks." This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGNs) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations in BHXRBs and changing-look phenomena in AGNs. Here, we present the first radiative two-temperature general relativistic magnetohydrodynamic (GRMHD) simulation of a strongly tilted (65°) accretion disk around an M BH = 10 M ⊙ BH, which tears and precesses. This leads to luminosity swings between a few percent and 50% of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation-pressure-dominated, the accretion disk is thermally stable over t ≳ 14,000 r g /c. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases tenfold and the electron temperature reaches T e ~ 10 8 –10 9 K. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to T e ~ 108 K. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The temperatures of ejecta transporting in vacuum and gases

In this work, we measure continuous thermal radiance from evolving clouds of liquid metal fragments ejected into vacuum, nonreactive, and reactive gas. We implement a model for the thermalization of the ejecta and gas and use this to constrain the absolute temperature of the ejecta cloud. This model enables further analyses of ejecta thermal behavior under a variety of conditions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Global Kinetic Modeling of the Intrabinary Shock in Spider Pulsars

Abstract Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and light curves, which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, F ν ∝ const; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by ∼0.8 rad; (3) the peak flux exceeds the one at the inferior conjunction by a factor of 10.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Isopropanol dehydration reaction rate kinetics measurement using H 2 O time histories

H 2 O formation during thermal decomposition of isopropanol was measured behind reflected shock waves at temperatures ranging from 1127 to 1621 K at an average pressure of 1.42 atm using a laser absorption technique. Of the five modern chemical kinetics models used to compare the H 2 O time histories, the model from Li et al (Combust. Flame 2019;207:171-185) showed the best overall agreement. Sensitivity and rate of production analyses using the Li et al model (as well as those from AramcoMech 3.0, CRECK, and Togbe et al (Energy Fuel. 2011;25:676-683)) showed unimolecular dehydration of isopropanol, iC 3 H 7 OH ⇌ H 2 O + C 3 H 6 (R1), is nearly the sole reaction controlling H2O production at early times, allowing for an a priori measurement of the forward rate constant k 1 . The Arrhenius expression k1 (s –1 ) = 2.60 × 10 13 exp(-31 120 K/T) was determined to represent best the data from this study. According to the models and previous experimental investigations, the pressure investigated is well within the high-pressure limit (HPL) for this reaction. Additional, higher-pressure experiments also confirmed the HPL assumption. An uncertainty analysis was performed by varying secondary reactions within their uncertainties and examining their effect on the overall prediction, establishing an uncertainty within ±20% for all but the highest temperature cases, which have a maximum uncertainty of ±40%. Experiments conducted with a radical trapper, toluene, showed little influence from radical chemistry, suggesting this estimated uncertainty is fairly conservative. Experimental data from Heyne et al (Z Phys Chem. 2015;229:881-907) were found to be in good agreement with the rate measurements from this study and, therefore, a second Arrhenius expression, k 1 (s –1 ) = 2.11 × 10 13 exp(-30 820 K/T), was found to represent both datasets well. This second expression has a larger temperature range of 976-1621 K. Here the present study provides the first high-temperature data collected for this reaction, adding to the limited data available for isopropanol in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A multi-wavelength streaked optical pyrometer for dynamic shock compression measurements above 2500 K

This work presents a concept and implementation of a calibrated streaked spectral pyrometer (SSP) temperature diagnostic used in dynamically driven shock experiments on a two-stage gas gun. This system relies upon measuring the total system response using a tunable monochromator, a NIST-traceable calibrated power meter, and a SSP. The diagnostic performance is validated against previously measured temperatures of shock driven z-cut quartz at 99 GPa and 93 GPa. The results are found to agree with the literature to within 5% and are discussed in this manuscript. The experimental setup utilizes measurements from multiple SSP systems per sample, providing several independent measurements and substantially increasing confidence in the extrapolated shock driven sample temperature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

VIPIR: A High-Throughput Drop-Weight Impact Instrument for Imaging the Initiation and Propagation of Reactions in Energetic Materials

Characterizing the handling safety and sensitivity of explosives has been a challenging area of study for over 60 years. Historically one of the most accessible and widely utilized experiments has been the drop-weight impact test, which involves dropping a weight on a small sample sandwiched between two anvils. Because this experiment generally only utilizes sound thresholds to determine whether or not a sample reacted, the physical and chemical properties governing sensitivity remain convolved. Better understanding of chemical and material characteristics is needed to give the chemistry and engineering communities a predictive tool to determine the handling sensitivity of explosives prior to pursuing expensive and potentially hazardous synthesis and formulation operations. Here, we are developing a high throughput drop tower instrument capable of imaging the deformation and flow of energetic materials during impact and the resulting thermal ignition and propagation events. This instrument is based on previous designs but has been modified for higher throughput and tailorable modifications in the future. Herein, we present key design features that enable high-speed visible and thermal imaging of explosive initiation by sub-shock impacts, as well as preliminary results in which ignition sites were observed in an erythritol tetranitrate sample.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

The role of sample–window coupling in shock temperature experiments

Radiometric studies of opaque materials under shock compression probe a sample–window interface that may be very different than the bulk state. Subtle experimental details, particularly an initial sample–window gap vs an adhesively bonded stack, have tremendous influence on measured temperature. This work experimentally tests gapped vs bonded targets and theoretically investigates observed differences in the measured temperature. The results indicate that perturbations created by a thin adhesive are smaller than those suggested by a previous work, whereas sample–window gaps alter the measured state more significantly than that previously realized.

Adhesives↗
DOE OSTI2020

Step-patterned survivorship curves: Mortality and loss of equilibrium responses to high temperature and food restriction in juvenile rainbow trout ( Oncorhynchus mykiss )

While survivorship curves typically exhibit smooth declines over time, step-patterned curves can occur with multiple stressors within a life stage. To explore this process, we examined the effects of heat (24°C) and food restriction on juvenile rainbow trout (Oncorhynchus mykiss Walbaum) in challenge experiments. We observed step-patterned survivorship curves determined by mortality and loss of equilibrium (LOE) endpoints. To examine the cause of heterogeneity in the stress responses from early to late mortality and LOE, we measured indices of energetic reserves. The step transition in the survivorship curves, the peak mortality rates, and start of when individuals reached a critical energetic threshold (14% dry mass; 4.0 kJ∙g -1 energy) all occurred at around days 10–15 of the challenge. The coherence in these temporal patterns suggest heterogeneity in the cohort stress responses, in which an early subgroup died from heat stress and a late subgroup died from starvation. Thus, their endpoint sensitivities resulted in step-patterned survivorship curves. We discuss the implications of the study for understanding effects of multiple stressors on population heterogeneity and note the possible significance of stress response selection under climate change in which heat stress and food limitations occur in concert.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

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↗
DOE OSTI2023

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Melting and density of MgSiO 3 determined by shock compression of bridgmanite to 1254GPa

The essential data for interior and thermal evolution models of the Earth and super-Earths are the density and melting of mantle silicate under extreme conditions. Here, we report an unprecedently high melting temperature of MgSiO 3 at 500 GPa by direct shockwave loading of pre-synthesized dense MgSiO3 (bridgmanite) using the Z Pulsed Power Facility. We also present the first high-precision density data of crystalline MgSiO 3 to 422 GPa and 7200 K and of silicate melt to 1254 GPa. The experimental density measurements support our density functional theory based molecular dynamics calculations, providing benchmarks for theoretical calculations under extreme conditions. The excellent agreement between experiment and theory provides a reliable reference density profile for super-Earth mantles. Furthermore, the observed upper bound of melting temperature, 9430 K at 500 GPa, provides a critical constraint on the accretion energy required to melt the mantle and the prospect of driving a dynamo in massive rocky planets.

58 GEOSCIENCES↗
DOE OSTI2026

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗
DOE OSTI2023

Structures of strong shocks in low-density helium and neon gases

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture (90% helium, 10% neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than ~350 µm with the helium gas mixture but less than 30 µm in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Additionally, both effects explain the longer precursor layer in the helium shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗