Search NASA⌕ Search

SEARCH · Search NASA

Results for “DETONATION”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗
DOE OSTI2026

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

The Effects of Shockwave Pressures on Ultrafast Vibrational Energy Transfer in BNFF, a Hydrogen-Free Energetic Material

Energy conversion in energetic materials from shock-wave-induced lattice compression to bond breaking critically depends on vibrational coupling and energy transfer between intra- and intermolecular vibrations, though the details of the mechanisms remain unknown. Herein, we indirectly tune the strength of intermolecular interactions in 3,4-bis(3-nitrofurazan-4-yl)furoxan (BNFF), a hydrogen-free energetic material characterized by van der Waals interactions, by applying high static pressure using a diamond anvil cell and monitoring vibrational energy transfer (VET) with ultrafast broadband infrared pump–probe spectroscopy. As BNFF is compressed from ambient pressure to 9 GPa, we find that VET accelerates by ∼ 0.9 ps/GPa. Density functional theory is applied in tandem with experiments to assign mode character and elucidate VET pathways. In conclusion, we find that furazan ring O–N–O vibrations, which are high-frequency detonation-relevant vibrational modes, experience increased sensitivity to lattice compression under shockwave pressures. These findings provide new mechanistic insight into how intermolecular interactions govern the rate and selectivity of VET.

Energy transfer↗
DOE OSTI2025

Quantitative Encapsulation and Homogeneity Assessment of Sol–Gel Based Nuclear Explosive Debris Simulants

Nuclear explosive debris simulants are an important material in training and validating aspects of post-detonation nuclear forensic processes. Realistic simulants should replicate several aspects of nuclear explosive debris such as the size, shape, color, density, and chemical and radiological properties. Silica particles produced via sol-gel synthesis have recently been found to successfully reproduce many of these parameters including the controllable incorporation of radionuclide content. However, to be useful as a benchmarking material for validation and verification of laboratory methodologies, radionuclide content from batch-to-batch must be reproducible. Here, in this work, we explore the variance in radionuclide distribution incorporated into sol-gel benchmarking materials with respect to sample subdivision. Results will help inform the sample sizes required to minimize variance between samples.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Direct observation of diamond formation in a shock-compressed high explosive

Understanding the formation timescale and structure of carbonaceous reaction products is critical for modeling the high-pressure equation-of-state of organic materials. We use the National Ignition Facility to shock-compress polycrystalline TATB (C 6 H 6 N 6 O 6 ) samples to ~70–130 GPa and ~4000–5500 K, employing in situ nanosecond X-ray diffraction to probe reaction products and velocimetry to measure transmitted compression wave profiles. Our diffraction data is consistent with the formation of diamond over timescales less than ~60 ns. This represents carbon condensation from a molecular explosive on timescales three times faster than previously reported and the earliest observation of diamond produced from reacting TATB. Reactive flow simulations with explicit chemistry reproduce the observed temporal structure within wave profiles to inform the distribution of P-T states. These findings provide direct evidence of ultrafast diamond formation in a reactive system at extreme conditions and provide new constraints for models of shock and detonation chemistry.

Clarke, Samantha M. [Lawrence Livermore National L↗
DOE OSTI2026

Violent mergers revisited: The origin of the fastest stars in the Galaxy

Binary systems of two carbon-oxygen white dwarfs are one of the most promising candidates for the progenitor systems of Type Ia supernovae. Violent mergers, where the primary white dwarf ignites when the secondary white dwarf smashes into it while being disrupted on its last orbit, were the first double degenerate merger scenario proposed that ignites dynamically. However, violent mergers likely contribute only a few percent to the total Type Ia supernova rate and do not yield normal Type Ia supernova light curves. Here we revisit the scenario, simulating a violent merger with better methods and, in particular, a more accurate treatment of the detonation. We find good agreement with previous simulations but with one critical difference: The secondary white dwarf being disrupted and accelerated towards the primary white dwarf and impacted by its explosion does not fully burn, and its core survives as a bound object. The explosion leaves behind a 0.16 M ⊙ star travelling 2800 km/s, making it an excellent (and so far the only) candidate to explain the origin of the fastest observed hypervelocity stars. We also show that before the explosion, 5 × 10 −3 M ⊙ of material predominantly consisting of helium, carbon, and oxygen had already been ejected at velocities above 1000 km/s. Finally, we argue that if a violent merger made the hypervelocity stars D6-1 and D6-3 and violent mergers require the most massive primary white dwarfs in binaries of two carbon-oxygen white dwarfs, there has to be a much larger population of white dwarf mergers with slightly lower mass primary white dwarfs. Because this population likely represents ≫10% of the Type Ia supernovae rate, it can essentially only give rise to normal Type Ia supernovae.

Astronomy and AstroPhysics↗
DOE OSTI2024

An exploration of anomalous electrical noise in shocked cyclotrimethylenetrinitramine (RDX)-based explosives

Gas gun shock experiments on cyclotrimethylenetrinitramine (RDX)-based explosive compositions that employ embedded gauge particle velocity tracers have noted a significant amount of electrical noise when compared to other explosive formulations. This paper reexamines previously published embedded gauge data on Cyclotols (60–80 wt. % RDX) to quantify the electromagnetic behavior of these materials. The primary observation is a fourfold increase in the electrical noise when Cyclotols are shocked above 4.22 ± 0.08 GPa. Electromagnetic gauge noise is also observed within particle velocity traces in reactive growth and off-Hugoniot shocks, although at higher pressures than the direct shock case, suggesting a temperature- or kinetically dependent transition. In all cases, the electrical noise disappears upon detonation. By comparing with the static high-pressure phase diagram of RDX, we interpret this change in electromagnetic behavior to be a change in the RDX crystal structure to a piezoelectric phase, although it is uncertain whether the γ or ε phase is responsible for the observed behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A generalization of the shock invariant relationship

Shock invariant relationship, which was conceived for inert shock waves to derive the 4th power relationship between shock pressure and maximum strain rate, is generalized for reactive shock waves such as Chapman–Jouget detonation and shock-induced vaporization. The generalization, based on the first-order reaction models, is a power function relationship between overall dissipated energy (Δe dis ) and reaction time Δτ such that Δe dis Δτ 1/α = constant, where the power coefficient α is found to be in the range of 2/3–4. Experimental data, though scarce, are consistent with the generalization. Implication of the generalization for inert shocks is also considered and suggests a broad range of the 4th power coefficient including an inequality equation that constrains the shock and particle velocity relationship.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗
DOE OSTI2026

Machine Learning–Augmented Laser-Induced Breakdown Spectroscopy for Spectral Discrimination of Iron Oxalates

Enhanced characterization and phase identification of post-PUREX Pu Oxalates (PuOXA) are pivotal for nonproliferation and pre-detonation nuclear forensics. Despite significant advances in the characterization of PuO 2 samples, little is known about the impact of both the chemical structure and oxidation states of PuOXA (i.e., Pu(III) and Pu(IV)) have on optical emission signatures. Here, we demonstrate the analytical capabilities of laser-induced breakdown spectroscopy (LIBS) applied to Fe(II) and Fe(III) oxalate samples as surrogates for PuOXA, highlighting the discriminating features in the LIBS emission spectra arising from differences in the oxidation states within mixed FeOXA samples. We report the enhancement of spectral feature selection using Principal Component Analysis (PCA), which enables the analytical superiority of machine learning algorithms such as Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest Regression (RFR) over conventional univariate techniques for phase discrimination and chemometric analysis. Cluster analysis revealed how both matrix effects and laser ablation influence cluster separability by introducing spectral artifacts that misdirect the maximization of variance. PCA-selected emission lines were used in the regression models, demonstrating that both univariate and multivariate linear regression models (i.e., PLSR and SVR) can achieve acceptable performance, with machine learning models outperforming conventional calibration regressions. Furthermore, the application of non-linearly activated PCA-selected emission lines illustrates how simplifying the data while retaining captured variance enables the use of less complex and more computationally efficient models. Furthermore, this is particularly evident in the underperformance of RFR, which suffers from increased computational costs and overfitting owing to its high complexity.

Oxalates↗
DOE OSTI2026

Mechanistic Investigation of Co(II) Extraction by TODGA to Aid Nuclear Forensic Separations

Nuclear forensic (NF) analysis supports law enforcement inquiries by analyzing evidence tainted with radioactive substances. Separation techniques can be used to identify and quantify actinides and fission 15 products in post-detonation (PD) debris. Environmental transition metals, also present in PD residues, have been observed to impact critical isotope extractions. For example, radio stable cobalt (Co), ubiquitous in urban environments, particularly in corrosion-resistant alloys, paint-drying agents, dyes, and pigments, can impact the separation of important actinides and fission products. The presented work aims to elucidate the chemistry Q2 governing Co extraction in samples pertinent to PDNF. Chemistry between Co and N,N,N,‘N’-tetraoctyl diglyco- 20 lamide (TODGA), the ligand present in the commercial chromatographic resin diglycolamide (DGA), were studied via solvent and chromatographic extraction and spectroscopic analyses. These results indicate that a tetrahedral Co(II) species is extracted by TODGA from highly acidic (>5M HCl) solutions via a spontaneous entropy-driven reaction. Furthermore, extraction trends in varied acid concentrations are consistent between solvent extraction and chromatographic extraction methods.

activation product↗
DOE OSTI2024

Underwater unexploded ordnance discrimination based on intrinsic target polarizabilities – A case study

Seabed unexploded ordnance that resulted partly from the high failure rate among munitions from more than 80 years ago and from decades of military training and testing of weapons systems poses an increasing concern all around the world. Although existing magnetic systems can detect clusters of debris, they are not able to tell whether a munition is still intact requiring special removal (e.g. in situ detonation) or is harmless scrap metal. The marine environment poses unique challenges, and transferring knowledge and approaches from land to a marine environment has not been easy and straightforward. On land, the background soil conductivity is much lower than the conductivity of the unexploded ordnance and the electromagnetic response of a target is essentially the same as that in free space. For those frequencies required for target characterization in the marine environment, the seawater response must be accounted for and removed from the measurements. The system developed for this study uses fields from three orthogonal transmitters to illuminate the target and four three-component receivers to measure the signal arranged in a configuration that inherently cancels the system's response due to the enclosing seawater, the sea–bottom interface and the air–sea interface for shallow deployments. The system was tested as a cued system on land and underwater in San Francisco Bay – it was mounted on a simple platform on top of a support structure that extended 1 m below and allowed the diver to place metal objects to a specific location even in low-visibility conditions. The measurements were stable and repeatable. Furthermore, target responses estimated from marine measurements matched those from land acquisition, confirming that the seawater and air–sea interface responses were removed successfully. Thirty-six channels of normalized induction responses were used for the classification, which was done by estimating the target principal dipole polarizabilities. Our results demonstrated that the system can resolve the intrinsic polarizabilities of the target, with clear distinctions between those of symmetric intact unexploded ordnance and irregular scrap metal. The prototype system was able to classify an object based on its size, shape and metal content and correctly estimate its location and orientation.

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

YOLO11 to SAM2 pipeline for feature extraction from nuclear test films

The response to the effects of nuclear detonations is supported by models that describe the evolution of the nuclear fireball and cloud and the associated transport of active debris. Validation of those descriptions relies on data from the nuclear test operations. Video records of those events offer a rich source of information that was exploited to a limited extent in historic analyses. Computer vision and machine learning techniques are powerful tools that can be used to increase the number of measurements that can be obtained from those films. In this work, we apply computer vision techniques to automatically track the temporal evolution of the nuclear fireball. In particular, we apply You Only Look Once 11 (YOLO11) and Segment Anything Model 2 (SAM2) in combination with minimal human intervention to digitized versions of the original nuclear test films. As part of the proposed workflow, the YOLO11 model is applied to films to determine bounding boxes for the fireball within each frame. These are then used as inputs to SAM2, which uses image segmentation to determine the fireball boundaries and their temporal evolution. We assess the accuracy of our approach by using it to determine the energy released during the Trinity nuclear test and comparing the results with previous analyses based on manual measurements.

Van Exel, Kimberly [ORNL] (ORCID:0009000877463894)↗
DOE OSTI2024

Codes for "Shallow Soil Response to a Buried Chemical Explosion with Geophones and Distributed Acoustic Sensing" DAG - 01101646

The codes reproduce the figures of the manuscript entitled "Shallow Soil Response to a Buried Chemical Explosion with Geophones and Distributed Acoustic Sensing" submitted to Journal of Geophysical Research - Solid Earth. Geophone data and Distributed acoustic sensing (DAS) data recorded during the Phase II of the The Source Physics Experiment (SPE) along a fiber-optic cable offshore were processed to understand the response of the shallow subsurface to an explosion. This Ground-based Nuclear Detonation Detection (GNDD), Low Yield Nuclear Monitoring (LYNM), and Source Physics Experiment (SPE) research was funded by the National Nuclear Security Administration, Defense Nuclear Nonproliferation Research and Development (NNSA DNN R&D).

Viens, Loic↗
DOE OSTI2024

ezECM

Allows the user to fit a classical event categorization matrix (ECM) model, and a novel Bayesian event categorization matrix model, both of which are used for nuclear detonation detection.

Koermer, Scott↗
DOE OSTI2026

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗