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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

AWSD Reactive Burn Model for High Explosive LX‐14

ABSTRACT The results of an Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow calibration for the HMX‐based high explosive LX‐14 are presented. The parameters in the AWSD model are calibrated to experimental thermodynamic and gas gun data and to computational results from thermochemical calculations. There is no experimental rate stick data available for LX‐14; therefore, scaled experimental results from other PBX‐based high explosives are used in the calibration to fill this gap in data. Strong agreement is observed between the calibrated AWSD model and experimental data for LX‐14, including validation data that were not used in the calibration procedure. The developed model more accurately describes experimental shock‐to‐detonation results compared to several other reactive flow models for LX‐14 from the literature. The presented results illustrate that the AWSD model is capable of quantitatively describing the reactive burn of LX‐14.

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

Data‐Efficient Generation of Synthetic Microstructures of Polymer‐Bonded Energetic Material With Fine‐Tuned Stable Diffusion

Among current deep learning approaches for synthetic image generation, diffusion-based models stand out in terms of algorithmic stability and ability to retain high-fidelity image features with detailed resolution. Here, in this work, we employ Dreambooth, a method for fine-tuning Stable Diffusion, on X-ray CT images of microstructure of the polymer-bonded form (PBX) of a commonly used high explosive, Pentaerythritol tetranitrate (PETN), which yields generative models for creating synthetic PBX images. The models developed here represent five classes (or ‘lots’) of microstructures and demonstrate successful generation of images of each class with high fidelity, as verified by computed classification accuracy of ∼ 94% or higher. Data augmentation afforded by such image synthesis can be used to more reliably decipher underlying statistics, build processing-structure correlations, recognize off-normal structural anomalies, and identify age-related changes. Ideas related to converting image data into appropriate density mapping and performing mesoscale simulation or surrogate modeling of detonation are also discussed.

Dreambooth↗
DOE OSTI2026

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

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

Shock Response of Polymethyl Methacrylate (PMMA) Under Explosive Loading

Abstract The dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.

Torres, S. M. (ORCID:0009000287177206)↗
DOE OSTI2026

High sensitivity measurement of femtogram-level 238 Pu by thermal ionization mass spectrometry with correction of background 238 U

Low-concentration plutonium isotopic measurements provide a valuable fingerprint to identify possible material origins in the context of environmental contamination monitoring, environmental tracing, nuclear safeguards, nuclear forensics, and treaty monitoring. Thermal ionization mass spectrometry (TIMS) has long been recognized as one of the most sensitive techniques for plutonium isotopic measurement, but has not generally been utilized for 238 Pu determination, especially at low concentrations, due to an isobaric interference from 238 U, prevalent as background. Here, this work demonstrates a new analytical technique for correction of background 238 U from the 238 Pu counts collected during a TIMS analysis. The technique takes advantage of the higher ionization temperature of U relative to Pu and uses a 235 U tracer added after plutonium purification chemistry for 238 U semi-quantitative determination. This technique has been successfully applied to the determination of 238 Pu concentration and the 238 Pu/ 239 Pu ratio in sample types including single-isotope standards, isotopic certified reference materials, matrix-containing environmental reference materials, and simulated nuclear detonation debris. The results have been directly compared on a sample-by-sample basis to alpha spectrometry results. The detection limit by this new technique is 0.23 fg 238 Pu, which was previously unachievable by mass spectrometry and is possible in around 1 h of analysis time post-chemistry compared to weeks of required counting time by alpha spectrometry (at the same atom amount). We also present original, high precision data for Pu isotope concentrations in IAEA-384, Fangataufa Sediment, including 242 Pu results for the first time. The new technique allows measurement of a key isotope, 238 Pu, that was previously not measured in small environmental or nuclear forensics collections.

238Pu↗
DOE OSTI2024

Separation of gold from irradiated actinide material utilizing a 2- thenoyltrifluoroacetone extraction resin

Post-detonation nuclear forensics capabilities depend on the ability to rapidly isolate radionuclides to improve measurement quality. In this work an extraction chromatography resin was developed utilizing the extractants thenoyltrifluoroacetone and 1-octanol supported on Eichrom prefilter resin. Here the resin was tested in nitric and hydrochloric acid matrices. In nitric acid matrices the resin is able to extract zirconium, while in hydrochloric acid matrices it is possible to extract iron and gallium. In all acid conditions tested gold is retained but can be removed from the column with a 10% thiourea rinse.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗
DOE OSTI2026

Enabling probabilistic learning on manifolds through double diffusion maps

Here, we present a generative learning framework for probabilistic sampling that extends Probabilistic Learning on Manifolds (PLoM), which is designed to generate statistically consistent realizations of a random vector in a finite-dimensional Euclidean space, informed by a (representative) set of observations. In its original form, PLoM constructs a reduced-order probabilistic model by combining three main components: (a) kernel density estimation to approximate the underlying probability measure, (b) Diffusion Maps to characterize the manifold of the data, and (c) a reduced-order Itô Stochastic Differential Equation (ISDE) to sample from the learned distribution. However, its sampling dynamics are posed in the ambient space and the retained number of reduced coordinates is chosen by projection-reconstruction error. In practice, this often (i) requires more coordinates than the data’s intrinsic dimension to achieve stable sampling and (ii) lacks a smooth, basis-independent lifting back to the data domain; moreover, standard Diffusion Maps emphasize harmonic eigenfunctions and can miss non-harmonic latent structure. We address these limitations by decoupling geometry learning from sampling: a first Diffusion Maps pass identifies non-harmonic coordinates on which we formulate a full-order ISDE directly in the latent space, while Double Diffusion Maps captures multiscale geometric features and Geometric Harmonics (GH) learns a smooth lifting map to the ambient variables that is independent of the particular diffusion basis. This hybrid design preserves the system’s dynamical richness with a compact geometric representation and enables principled out-of-sample inference. The effectiveness and robustness of the proposed method are illustrated through two numerical studies: one based on data generated from two-dimensional Hermite polynomial functions and another based on high-fidelity simulations of a detonation wave in a reactive flow.

Double diffusion maps↗
DOE OSTI2025

Analysis of 127 Xe tracer measurements using a net counts method

A suite of measurement systems were deployed as part of the Physical Experiment 1 series of experiments, which involved detonating chemical explosives along with radionuclide tracers in an underground cavity, at the Nevada National Security Site (NNSS) in the United States. One of the radionuclide tracers, 127 Xe was released from the containment following the explosion and detected on a SAUNA Q B sampler situated approximately 3.5 km away. The system uses a beta-gamma coincidence detector system to measure fission product radioisotopes of xenon relevant to nuclear explosion monitoring. In this work we use the coincidence measurement data to analyse and interpret the results from the SAUNA Q B system, to calculate the measured 127 Xe activity concentration(s).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

U-Pu and Ba-Cs isotopic measurements on Trinitite by laser ablation sampling on the Neoma MC-ICP-MS

In this study we present the results of combined U-Pu and Ba-Cs isotope measurements obtained by laser ablation (LA) sampling of two glassy debris fragments (‘Trinitite’) from the world's first atomic bomb detonation conducted in New Mexico on July 16, 1945. Our primary goal in conducting these measurements was to understand whether examination of the U-Pu and Ba-Cs systematics by direct sampling (e.g. without any chemical separation or purification prior to isotope ratio measurement) could yield meaningful information that would differentiate the Trinitite fragments from glassy material lacking a nuclear fission signature. These measurements were conducted on a ThermoFisher Scientific Neoma multi collector – inductively coupled plasma – mass spectrometer (MC-ICP-MS), which is a relatively new MC-ICP-MS platform, so we also examine the behavior of these isotope systems in standards sampled in solution and via LA. Unsurprisingly, the measurements made on purified solutions of the U, Pu, and Ba isotopic standards produce high precision isotope ratios. Furthermore, this extends to the U-Pu measurements made by LA sampling, with the expected degradation in precision and accuracy related to matrix effects and signal intensity fluctuation. However, the Ba-Cs data acquired by LA is of low precision across all of the matrices examined and bears evidence of complex mass fractionation that will require further investigation to resolve. In total, our results indicate that the observed U-Pu isotope data are of sufficient quality to accurately constrain the U and Pu isotopic composition of glass containing sub-ppm levels of these elements which in turn could be used to differentiate glass containing anthropogenic fission products from natural glass whereas the Ba-Cs LA data cannot be used for this purpose until further methodological refinement is performed.

Ba-Cs↗
DOE OSTI2024

A new method for quantifying 64 Cu in nuclear debris samples

Quantifying 64 Cu in post-detonation nuclear debris samples can provide important diagnostic information regarding the structural materials used within a nuclear device. However, this task is challenging due to the weak gamma emissions associated with the decay of 64 Cu, its short half-life (12.701 h), and the presence of interfering fission product radioisotopes. Large quantities of debris sample are generally needed to accurately quantify 64 Cu, which can be problematic in sample-limited scenarios where other radiometric analyses are required. Herein, we present a new method for the separation of 64 Cu from solutions of mixed fission products and demonstrate the quantification of its activity through use of gas-flow proportional beta counting. Here, the new method was validated through a series of rigorous tests and was shown to improve the detection limit of 64 Cu by over two orders of magnitude, from 2.5 × 10 6 to 1.3 × 10 4 atoms/sample for 100 min measurements.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Rapid analysis of 237 Np and Pu isotopes in unseparated sample matrices using ICP-MS/MS

Inductively coupled plasma tandem mass spectrometry (ICP-MS/MS) is an emerging technique for measuring actinide isotopes when assessing pre- and post- detonation nuclear material. In this study, ICP-MS/MS was investigated for direct Np and Pu quantitation in unseparated, dissolved bulk soil matrices. To achieve this, purified nitric oxide (NO) was investigated for the reactivity of Th, Np, U, Pu, Am, and Cm. Purifying NO prior to the collision reaction cell (CRC) results in increased sensitivity and allows for higher gas flows to be utilized for the removal of interferences. Here, the interference from uranium hydrides was mitigated to less than 3.85 x 10 -11 . This method was demonstrated on standard reference materials which were measured for 237 Np and 238,239,240 Pu in dilute sample digestions. The 238 Pu measurement was validated by spiking into a standard reference material and was accurately measured with an excess of 85000 of 238 U.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

147 Nd Quantification Using HSCCC-Purified Samples

Quantifying the fission product 147 Nd in nuclear debris samples is an important component of post-detonation nuclear forensics. The most accurate quantifications are obtained when Nd is purified from all other fission products, actinides, activation products, and environmental matrix contained within the debris. In this study, a recently developed method for Nd purification was tested, purifying 147 Nd from solutions of mixed fission products using high-speed counter-current chromatography (HSCCC). Importantly, the new method allowed for faster elution of Nd from the column as compared with established high performance liquid chromatography (HPLC) methods, and resulted in accurate/precise 147 Nd quantification by gamma-ray spectrometry. While the up-front equipment costs associated with HSCCC may be higher, its operational costs are on par with those of HPLC (solvents, extractants, power). Gas-flow proportional beta decay counting revealed contamination from the nearest neighbor lanthanide 143 Pr (a gamma-silent radioisotope) in the HSCCC-purified samples, but the activity contribution from 147 Nd could still be quantified. Remarkably consistent elution profiles were observed for the HSCCC method, spanning rare earth element (REE) loadings of more than 10 orders of magnitude (tracer to mmol quantities). In conclusion, the reliability and speed of the new method suggest utility for the rapid separation and quantification of 147 Nd in unknown samples.

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

Assessing Gas-Phase Ion Reactivity of 50 Elements with NO and the Direct Application for 239 Pu in Complex Matrices Using ICP-MS/MS

Understanding the reactivity of metal cations with various reaction gases in ICP-MS/MS is important to determine the best gas to use for a given analyte/interference pair. In this study, nitric oxide (NO) was investigated as the reaction gas following previous experimental designs. The reactions with 50 elements were investigated to examine periodic trends in reactivity, validate theoretical modeling of reaction enthalpies as a method to screen reactant gases, and provide a baseline for potential in-line gas separation methods. ICP-MS/MS studies involving actinides are typically limited to Th, U and Pu, with analyses of Np and Am rarely reported in the literature. To date, only two previous methods have investigated the use of NO in ICP-MS/MS analyses. To showcase the utility of NO, a method was developed to measure 239 Pu in the presence of environmental matrix constituent and other actinides, like what could be expected from post-detonation debris, with no chemical separation prior to analysis. 239 Pu + was reacted to form 239 Pu 16 O + , eliminating interferences derived from the sample matrix by measuring 239 Pu + intensity at m/z = 255 ( 239 Pu 16 O + ). To validate NO for 238 U 1 H + interference removal in environmental matrices, standard reference materials were diluted to 1 mg per gram of solution and spiked to 0.05 pg/g of 239 Pu and 1 µg/g 238 U (Pu/U = 5 x 10 -8 ). Measured 239 Pu concentrations were within 6% of the spiked value. In conclusion, these results demonstrate that reliable 239 Pu measurements can be made at levels relevant to nuclear forensics without the need for extensive chemical matrix separation prior to analysis.

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

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
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↗