Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular explosives”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Phosgene Neutralization Chemistry for the Explosive Destruction System (EDS)

This report summarizes activities at Sandia National Laboratories as part of the Explosive Destruction System (EDS) Phase 3 (P3) System design. An exploration of chemical neutralization strategies of phosgene was conducted for safe disposal of recovered mortars and M79 1000 lb. bombs filled with carbonyl dichloride "phosgene" or "CG agent" (molecular formula = COCl 2 ). The incumbent strategy utilized aqueous sodium hydroxide was found to be the worst-case scenario, producing enough CO 2 gas that would cause an unacceptable pressure and temperature spike. Several chemical neutralization strategies were evaluated based on criteria set by the operating envelope of the P3 design. In the end, it was determined that a pure solution of N-methyl ethanolamine (MeEA) or 90% monoethanolamine (MEA)(aqueous) provided the best balance reaction profile, cost, and safety.

36 MATERIALS SCIENCE↗
NASA NTRS2006

The 3D Structure of Eta Carinae's Nebula: A Definitive Picture from High-Dispersion Near-IR Spectra

High resolution long-slit spectra obtained with the Phoenix spectrograph on Gemini South provide our most accurate probe of the 3D structure of the Homunculus Nebula around Eta Carinae. Emission from molecular hydrogen at 2.122 microns traces a very thin outer skin, which contains the vast majority of the more than 10 solar masses of material in the nebula. This emission, in turn, yields our first definitive picture of the exact shape of the nebula, plus the latitude dependence of the mass-loss rate, speed, kinetic energy, shell thickness, and other properties associated with Eta Car's 19th century explosion. This will be critical for testing any models for the outburst mechanism. A preliminary analysis suggests that explosion from a critically rotating star was the dominant mechanism that gave rise to both the bipolar shape of the nebula and the production of its equatorial disk. [Fe II] emission in the near IR traces a geometrically thicker but less massive shell found on the inner surface of the H2 skin --- this is either a reverse shock that decelerates Eta Car's wind or a warm PDR. [Fe Ill emission also clarifies the structure of an inner "Little Homunculus" seen previously in HST/STlS spectra. Comparing these two tracers of cool molecular gas and warm partially-ionized gas resolves some significant confusion about the complex structure noted in previous studies.

Smith, N.↗
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↗
NASA NTRS1998

Spatially Resolved Spectroscopy of the SNR IC443

IC 443 is a supernova remnant of intermediate age, i.e. a few thousand years. It is especially interesting because part of its periphery is expanding into a molecular cloud while other sections are expanding into a typical interstellar medium of much lower density. Since the evolution of a supernova remnant through its various phases is affected by the density of the medium it expands into with the reasonable assumption that the supernova explosion was approximately symmetric we have an opportunity to observe a single object in two phases simultaneously. It was observed by ASCA in April, 1993 for a short period during the PV phase and more thoroughly in a 42 ksec exposure in March, 1994. The latter measurement provides most of the results that have been reported. Most of the analysis took place after the grant ended but is included here for completeness. The data was sent simultaneously to US and Japanese Pls. We worked independently. The software set of FTOOLs was used to construct images and spectra. They were judged to be rather unintuitive and not at all user friendly. I found I was using one FTOOL to read the header to obtain information that would only be provided to another FTOOL. The Japanese investigators were more successful. They analyzed the data and published results more rapidly. The scientific results summarized below are based primarily on their publications. Since IC 443 is an interesting example of a middle aged SNR in which a variety of processes are occurring it is one of a class. IC 443 exhibits shell-like emission in hard X-rays and extended soft X-rays with thin thermal spectra. It resembles SN 1006 in these respects. IC 443 contains hard X-rays in a semi-circular shell surrounding the thermal component. The total hard X-ray flux in the ASCA FOV is only a half of the Ginga hard component; which suggests that the hard X-rays are not confined only in the shell but some are extended larger than the ASCA FOV of eq 1 degree diameter. Japanese investigators examined the spatial structure of the thermal component and analyzed the GIS spectra with a non-equilibrium plasma model, and found no systematic variation of the interstellar absorption across the remnant. Evidence for shock acceleration of cosmic rays to high energies (10 TeV) was found by Keohane. X-ray imaging spectroscopy with ASCA reveals two regions of particularly hard emission: an unresolved source embedded in an extended emission region, and a ridge of emission coincident with the southeastern rim. Both features are located on part of the radio shell where the shock wave is interacting with molecular gas, and together they account for a majority of the emission at 7 keV. Though we would not have noticed it a priori, the unresolved feature is coincident with one resolved by the ROSAT HRI. The ASCA measurements were combined with higher energy data from the XTE and GRO missions and with radio and TeV gamma-ray data to produce a nonthermal multiwavelength spectrum for IC 443 which was fit with a cosmic ray interaction model. This model calculates the cynchrotron, bremsstrahlung, invers Compton, and neutral pion decay emission produced by locally accelerated cosmic ray interacting with ambient matter, soft photon fields, and magnetic fields.

Gorenstein, P.↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

An Updated Picture of Presolar History from Short-lived Radioactive Isotopes and Inferences on the Birth of the Sun

We examine the origin of the short-lived radionuclides (SLRs; defined as having half-lives between 0.1 and 100 Ma) present in the early solar system (ESS) by investigating how predictions of their abundances in the interstellar medium (ISM) from steady-state equilibrium relate to their ESS values. For this, we take into account the non-negligible time t iso elapsed between the isolation of the presolar molecular cloud and the formation of the ESS, during which the SLRs decayed freely. We also consider the alternative scenario in which the presolar molecular cloud remained partially mixed with the ISM, with a mixing timescale t mix . We find that the ESS abundances of 107 Pd and 182 Hf produced by slow neutron captures (s-process) and of 53 Mn and 60 Fe produced by explosive nucleosynthesis can be consistently explained within these scenarios. Their required tiso is 9–12 Ma and their required tmix is 11–14 Ma (with one potential exception of t mix = 38 Ma), depending on galactic uncertainties, such as the galactic star formation history and efficiency and the star-to-gas mass ratio. Another s -process SLR, 205 Pb, has a more uncertain ESS value and falls within only some of these time values. The same applies to the SLRs produced by the p -process ( 92 Nb and 146 Sm), depending on the latter’s half-life. In agreement with previous studies, we find that the ESS abundances of the rapid neutron-capture isotopes ( 129 I, 244 Pu, and 247 Cm) and of the most short-lived radionuclides ( 26 Al, 36 Cl, and 41 Ca) cannot be explained by assuming steady-state equilibrium in the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

Latest developments in cryogenic safety

The Cryogenic Safety Manual, sponsored by the British Cryogenics Council, was published over 10 years ago. A new updated version is now available. Some general aspects of cryogenic safety are highlighted, and attention is drawn to some of the more unusual hazardous situations. An awareness of the physical properties of the cryogenic fluids being dealt with is important in directing attention to hazardous situations which may arise. Because of this, the more important properties of the cryogenic fluids are given, such as molecular weight, boiling point and freezing point. From these properties, hazardous situations can be deduced. There are hidden dangers that are not always easy to spot. Some of the unexpected hazards, most of which have led to deaths, are: asphyxiation (anoxia), frost bites and hypothermia, explosions, and combustion. The aim of this publication is to help bring about increased safety in the production and use of cryogenic products through a deeper appreciation of the scientific, technological and administrative steps which must be made if accidents, some fatal, are to be voided in the future.

Webster, T. J.↗
NASA NTRS1982

Detonation propulsion experiments and theory

Test data are presented for the use of a single detonation of explosives in long-cone, short-cone, straight, and firing-plug nozzles to provide propulsion in a simulated Jupiter atmosphere, as well as the ambient gases N, CO2 and He. The long-cone nozzle yielded a progressive increase with ambient pressure for the higher molecular weight gases CO2 and N, while the lower molecular weight He and simulated Jupiter atmosphere showed a specific pulse decrease with increasing ambient pressure. The short-plug nozzle yielded a small specific impulse reduction with increasing ambient pressure, and its results were found to be nearly independent of ambient gas molecular weight. All data gathered are analyzed by using first principles, approximate blast wave theory predictions, and two-dimensional numerical calculations. Rarefaction and oscillatory wave phenomena are found to significantly influence specific impulse.

Back, L. H.↗
NASA NTRS1979

The depletion of interstellar gaseous iron

The Copernicus UV telescope was used to measure equivalent widths of interstellar Fe II resonance lines toward 55 early-type stars; the measurements permit the determination of Fe II column densities. The depletion of interstellar gaseous iron was obtained by combining these measurements with the results from a previous atomic and molecular hydrogen survey program; the derived depletions refer mostly to matter in H I regions. As an example, the nearly normal gaseous iron abundance in the distant high-latitude intermediate-velocity cloud toward HD 93521 is consistent with the idea that these clouds are produced by galactic supernova explosions.

Savage, B. D.↗
DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computational study of TATB

We have performed a computational study of TATB. The study is composed of two parts, one where we perform static T = 0 K calculations with several different DFT functionals, to investigate structural properties and obtain a cold curve. Even though the functionals used in this study give poor results at ambient pressure, they perform much better at compression, when the dominance of the van der Waals’ forces in the binding is replaced by more normal interactions. In particular AM05 gives good results, mainly because it doesn’t include any “faulty” van der Waals’. Indeed the van der Waals’ corrected functional (PBE+D2) results shows the worst performance on structural properties. For the cold curve the PBE+D2 functional gives the best results compared to experiments (room temperature) but since the structural properties were bad with this functional we cannot for sure say that this is in any way a validation. It is shown that scaling the AM05 results with a common factor reproduces the experimental data. We leave the question on why to a follow up study but speculate that maybe the experiment and our calculations do not represent the same system set-up. From the structural properties we identify a site in the TATB lattice that contains a high concentration of oxygen and hydrogen. The consequences of this finding, if any, are left to a future project. In the second part of the study we perform extensive high-quality DFT-MD (molecular dynamics with DFT forces) calculations to obtain points on the Hugoniot. This part of the study did not get finished within this project and a more extensive summary of this part will be given in a follow up report at a future date. However, we show that obtaining calculated data is essential for equations of state development since the existing experimental Hugoniot data barely deviates from room temperature isotherm data and thus give no information of the high temperature dependency of materials properties. Our conclusion so far in the project is that being able to calculate materials properties is essential for high explosives since experiments are unable to give information in the high temperature part of phase space. While blind use of DFT cannot (yet) give unambiguous answers, we can still explore qualitatively and use such information to guide us in EOS development. However, the main conclusion is that we need better functionals in DFT to easily extract needed properties from calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase

Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shock wave bending around a dusty plasma void

We report on experimental observations of the bending of a dust acoustic shock wave around a dust void region. This phenomenon occurs as a planar shock wavefront encounters a compressible obstacle in the form of a void whose size is larger than the wavelength of the wave. As they collide, the central portion of the wavefront, that is the first to touch the void, is blocked while the rest of the front continues to propagate, resulting in an inward bending of the shock wave. The bent shock wave eventually collapses, leading to the transient trapping of dust particles in the void. Subsequently, a Coulomb explosion of the trapped particles generates a bow shock. The experiments have been carried out in a Direct Current (DC) glow discharge plasma, where the shock wave and the void are simultaneously created as self-excited modes of a three-dimensional dust cloud. The salient features of this phenomenon are reproduced in molecular dynamics simulations, which provide valuable insights into the underlying dynamics of this interaction.

Dust acoustic waves↗
DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Type Ia supernovae: Pulsating delayed detonation models, IR light curves, and the formation of molecules

We computed optical and infrared light curves of the pulsating class of delayed detonation models for Type Ia supernovae (SNe Ia). It is demonstrated that observations of the IR light curves can be used to identify subluminous SNe Ia by testing whether secondary maxima occur in the IR. Our pulsating delayed detonation models are in agreement with current observations both for subluminous and normal bright SN Ia, namely SN1991bg, SN1992bo, and SN1992bc. Observations of molecular bands provide a test to distinguish whether strongly subluminous supernovae are a consequence of the pulsating mechanism occurring in a high-mass white dwarf (WD) or, alternatively, are formed by the helium detonation in a low-mass WD as was suggested by Woosley. In the latter case, no carbon is left after the explosion of low-mass WDs whereas a log of C/O-rich material is present in pulsating delayed detonation models.

Hoflich, Peter↗
DOE OSTI2026

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗
DOE OSTI2024

Distinguishing the XUV-induced Coulomb explosion dynamics of iodobenzene using covariance analysis

The primary and secondary fragmentation dynamics of iodobenzene following its ionization at 120 eV were determined using three-dimensional velocity map imaging and covariance analysis. Site-selective iodine 4d ionization was used to populate a range of excited polycationic parent states, which primarily broke apart at the carbon-iodine bond to produce I+ with phenyl or phenyl-like cations (C n H x + or C n H, with n = 1 – 6 and x = 1 – 5). The molecular products were produced with varying degrees of internal excitation and dehydrogenation, leading to stable and unstable outcomes. This further allowed the secondary dynamics of C 6 H x 2+ intermediates to be distinguished using native-frame covariance analysis, which isolated these processes in their own centre-of-mass reference frames. Furthermore, the mass resolution of the imaging mass spectrometer used for these measurements enabled the primary and secondary reaction channels to be specified at the level of individual hydrogen atoms, demonstrating the ability of covariance analysis to comprehensively measure the competing fragmentation channels of aryl cations, including those involving intermediate steps.

74 ATOMIC AND MOLECULAR PHYSICS↗