Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High 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 469 records · Page 26

Visualization of Hemispherical Post‐Detonation Fireball Internal Structures

Hemispherical charges are initiated above a transparent plate in a half‐plane configuration, allowing optical access to the internal regions of the luminous fireball. High‐speed visualization enables characterization of the internal luminous structure, showing clearly the bright shell and dark core regions. Spectroscopic and pyrometric diagnostics are applied to these flows, generating quantitative information on the spatial distribution of temperature inside the fireball. Both cased and uncased charges can be examined using this methodology, and several different high explosives are tested. This approach can be useful in validation of detailed explosive fireball models. Initial comparisons with such models are presented.

detonation↗

New Developments in Proton Radiography at the Los Alamos Neutron Science Center (LANSCE)

An application of nuclear physics, a facility for using protons for flash radiography, was developed at the Los Alamos Neutron Science Center (LANSCE). Protons have proven far superior to high energy x-rays for flash radiography because of their long mean free path, good position resolution, and low scatter background. Although this facility is primarily used for studying very fast phenomena such as high explosive driven experiments, it is finding increasing application to other fields, such as tomography of static objects, phase changes in materials and the dynamics of chemical reactions. The advantages of protons are discussed, data from some recent experiments will be reviewed and concepts for new techniques are introduced.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Μeasurements of Sn Ejecta Particle-Size Distributions Using Ultraviolet In-line Fraunhofer Holography

Sn ejecta particle-size distributions from the break-up of a microsheet in a vacuum will be presented. The micro-sheet was created from a high explosive driven shock wave passing through a precision groove machined into a Sn sample. The particle-size measurements were carried out using an ultraviolet in-line Fraunhofer holography diagnostic. In conclusion, the diagnostic will be presented along with particle-size distributions over most of the ejecta velocities throughout the microsheet.

36 MATERIALS SCIENCE↗

Fate and transport of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and its degradation products in sedimentary and volcanic rocks, Los Alamos, New Mexico

High-explosive compounds including hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) were used extensively in weapons research and testing at Los Alamos National Laboratory (LANL). Liquid effluents containing RDX were released to an outfall pond that flowed to Cañon de Valle at LANL's Technical Area 16 (TA-16), resulting in the contamination of the alluvial, intermediate and regional groundwater bodies. Monitoring of groundwater within Cañon de Valle has shown persistent RDX in the intermediate perched zone located between 225 and 311 m below ground surface. Monitoring data also show detectable levels of RDX putative anaerobic degradation products. Batch and column experiments were conducted to determine the extent of adsorption-desorption and transport of RDX and its degradation products (MNX, DNX, and TNX) in major rock types that are within the RDX plume. All experiments in this paper were performed in the dark using water obtained from a well located at the center of the plume, which is fairly oxic and has a neutral pH of 7.5. Retardation factors and partitioning coefficient (K d ) values for RDX were calculated from batch experiments. Additionally, retardation factors and K d values for RDX and its degradation products were calibrated from column experiments using a one-dimensional transport model with equilibrium sorption (linear isotherm). Results from the column and batch experiments showed little to no sorption of RDX to the aquifer materials tested, with retardation factors ranging from 1.0 to 1.8 and K d values varying from 0 to 0.70 L/kg. Finally, results also showed no measurable differences between the transport properties of RDX and its degradation products.

54 ENVIRONMENTAL SCIENCES↗

Embedded smoothed particle hydrodynamics

In this paper, a hybrid method for embedding a meshfree smoothed particle hydrodynamics into mesh-based ALE hydro schemes is presented. The method relies on the embedded mesh framework presented in Puso et al. (2015, 2012) and utilizes SPH (Owen, 2010) in the solid foreground; when an appropriate damage or failure model is used with SPH, it becomes a natural way to model material fracture and cracking in high-explosive driven experiments. Overall, we provide various examples which validate the coupling between SPH and the embedded mesh method.

42 ENGINEERING↗

High pressure-high temperature phase diagram of an energetic crystal: Dihydroxylammonium 5,5'-bistetrazole-1,1'-diolate (TKX-50)

Raman spectroscopy was used to examine the high pressure-high temperature structural and chemical stability of an insensitive, high-performance energetic crystal – dihydroxylammonium 5,5'-bistetrazole-1,1'-diolate (TKX-50). The phase diagram was determined over 8 GPa and (293-760) K. Under isobaric heating, the melting/decomposition of TKX-50 is preceded by a transformation to two consecutive high-temperature intermediates; a lower-temperature intermediate – diammonium 5,5’-bistetrazole-1,1'-diolate, and a higher-temperature intermediate – dihydroxylammonium 5,5'-bistetrazolate and/or diammonium 5,5'-bistetrazolate. Pressure strongly increases the transition temperatures for these transformations and subsequent decomposition. As a result, significant increase in the chemical stability of TKX-50 and intermediates with pressure was attributed to a suppression of hydrogen-transfer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation and analysis of the 3D double shock-induced recompression of spall using a hybrid PAGOSA with FLIP+MPM

In this work, an innovative hybrid algorithm PAGOSA with FLIP + MPM is first presented by coupling PAGOSA with the particle FLIP + MPM, and applied to systematically exploring the 3D recompression of spall induced by double-shock waves in ductile materials as the few existing studies are currently limited to the pure longitudinal recompression of spall, and ignore the 3D effects. This algorithm solves the difficulties the grid-based methods encounter when applied to capturing the fracture in material. We first validate the capabilities of PAGOSA with FLIP + MPM to predict different fracture/fragmentation cases by solving two benchmark problems and comparing with the analytical solutions or the experiment results. The convergences and the infinity-norm errors are also investigated. Subsequently, the 3D effects are demonstrated by simulating the spallation in material using PAGOSA with FLIP + MPM and comparing with those in the longitudinal cases. Additionally, some factors that affect the spallation with 3D effects are also discussed. Finally, the recompressions of spall in material with/without 3D effects are systematically simulated and analyzed. The recompression of spall driven by a high explosive is also simulated to show the ability of PAGOSA with FLIP + MPM to handle this kind of problem. The conditions that can lead to a recompression of spall are explored. Also, numerical results show that PAGOSA with FLIP + MPM can accurately predict different fracture cases, that the 3D effects play an important role in fracture/fragmentation, and the recompression of spall is very sensitive to the occurrence conditions, shows great differences when considering 3D effects in real applications. Moreover, the present hybrid method can be easily extended to other grid-based techniques employed for fracture in materials.

3D effects↗

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗

Detonation propagation for shock-driven, subsonic and supersonic confiner flow

In this work, we study the compressible flow dynamics of two-dimensional, steady detonation wave propagation in a high explosive (HE) confined by aluminium (Al) or stainless steel (SS), outside of which is an air layer. We examine how the thickness of the confinement affects the subsonic detonation driving zone structure (DDZ) and the detonation speed (D 0 ), demonstrating a strong dependence on whether the oblique shock-driven flow in the confiner is supersonic, as for SS, or subsonic, as for Al. A characteristic path analysis is used to examine the information flow from the material boundaries through the supersonic flow regions in both the HE and confiner that can impact the sonic surfaces bounding the subsonic flow regions. It is shown that the nature of gas-dynamic wave reflection off the SS–air or Al–air boundary can significantly influence the DDZ and D 0 .

54 ENVIRONMENTAL SCIENCES↗

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗

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↗

Predicting Energetics Materials’ Crystalline Density from Chemical Structure by Machine Learning

To expedite new molecular compound development, a long-sought goal within the chemistry community has been to predict molecules’ bulk properties of interest a priori to synthesis from a chemical structure alone. In this work, we demonstrate that machine learning methods can indeed be used to directly learn the relationship between chemical structures and bulk crystalline properties of molecules, even in the absence of any crystal structure information or quantum mechanical calculations. We focus specifically on a class of organic compounds categorized as energetic materials called high explosives (HE) and predicting their crystalline density. An ongoing challenge within the chemistry machine learning community is deciding how best to featurize molecules as inputs into machine learning models—whether expert handcrafted features or learned molecular representations via graph-based neural network models—yield better results and why. We evaluate both types of representations in combination with a number of machine learning models to predict the crystalline densities of HE-like molecules curated from the Cambridge Structural Database, and we report the performance and pros and cons of our methods. Our message passing neural network (MPNN) based models with learned molecular representations generally perform best, outperforming current state-of-the-art methods at predicting crystalline density and performing well even when testing on a data set not representative of the training data. However, these models are traditionally considered black boxes and less easily interpretable. Here, to address this common challenge, we also provide a comparison analysis between our MPNN-based model and models with fixed feature representations that provides insights as to what features are learned by the MPNN to accurately predict density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High pressure–temperature phase diagram of 1,1-diamino-2,2-dinitroethylene (FOX-7)

In this study, the pressure–temperature (P–T) phase diagram of 1,1-diamino-2,2-dinitroethylene (FOX-7) was determined by in situ synchrotron infrared radiation spectroscopy with the resistively heated diamond anvil cell (DAC) technique. The stability of high-P–T FOX-7 polymorphs is established from ambient pressure up to 10 GPa and temperatures until decomposition. The phase diagram indicates two near isobaric phase boundaries at ~2 GPa (α → I) and ~5 GPa (I → II) that persists from 25 °C until the onset of decomposition at ~300 °C. In addition, the ambient pressure, high-temperature α → β phase transition (~111 °C) lies along a steep boundary (~100 °C/GPa) with a α–β–δ triple point at ~1 GPa and 300 °C. A 0.9 GPa isobaric temperature ramping measurement indicates a limited stability range for the γ-phase between 0.5 and 0.9 GPa and 180 and 260 °C, terminating in a β–γ–δ triple point. With increasing pressure, the δ-phase exhibited a small negative dT/dP slope (up to ~0.2 GPa) before turning over to a positive 70 °C/GPa slope, at higher pressures. The decomposition boundary (~55 °C/GPa) was identified through the emergence of spectroscopic signatures of the characteristic decomposition products as well as trapped inclusions within the solid KBr pressure media.

36 MATERIALS SCIENCE↗

Ultrafast Photodissociation Dynamics of Nitromethane

Nitromethane (NM), a high explosive (HE) with low sensitivity, is known to undergo photolysis upon ultraviolet (UV) irradiation. The optical transparency, homogeneity, and extensive study of NM make it an ideal system for studying photodissociation mechanisms in conventional HE materials. The photochemical processes involved in the decomposition of NM could be applied to the future design of controllable photoactive HE materials. In this work, the photodecomposition of NM from the nπ* state excited at 266 nm is being investigated on the femtosecond time scale. UV femtosecond transient absorption (TA) spectroscopy and excited state femtosecond stimulated Raman spectroscopy (FSRS) are combined with nonadiabatic excited state molecular dynamics (NA-ESMD) simulations to provide a unified picture of NM photodecomposition. The FSRS spectrum of the photoproduct exhibits peaks in the NO 2 region and slightly shifted C–N vibrational peaks pointing to methyl nitrite formation as the dominant photoproduct. A total photolysis quantum yield of 0.27 and an nπ* state lifetime of ~20 fs were predicted from NA-ESMD simulations. Predicted time scales revealed that NO 2 dissociation occurs in 81 ± 4 fs and methyl nitrite formation is much slower having a time scale of 452 ± 9 fs corresponding to the excited state absorption feature with a decay of 480 ± 17 fs observed in the TA spectrum. Lastly, although simulations predict C–N bond cleavage as the primary photochemical process, the relative time scales are consistent with isomerization occurring via NO 2 dissociation and subsequent rebinding of the methyl radical and nitrogen dioxide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase diagram and decomposition of 1,1-diamino-2,2-dinitroethene single crystals at high pressures and temperatures

The high pressure-high temperature (HP-HT) phase diagram and decomposition of FOX-7, central to understanding its stability and reactivity, were determined using optical spectroscopy and imaging measurements in hydrostatically compressed and heated single crystals. Boundaries between various FOX-7 phases (α, α’, β, γ, and ε) and melting/decomposition curves were established up to 10 GPa and 750 K. Main findings are: (i) a triple point is observed between α, β, and γ phases ~ 0.6 GPa and ~ 535 K, (ii) previously suggested δ phase is not a new phase but is partly decomposed γ phase, (iii) the α-α’ transition takes place along an isobar, whereas the α’-ε transition pressure decreases with increasing temperature, and (iv) melting/decomposition temperatures increase rapidly with pressure, with an increase in the slope at the onset of the α’-ε transition. Our results differ from the recently reported HP-HT phase diagram for nonhydrostatically compressed polycrystalline FOX-7. In addition, the observed interplay between melting and decomposition suggests the suppression of melting with pressure. Our FTIR measurements at different pressures to 3.5 GPa showed similar decomposition products, suggesting similar decomposition pathways irrespective of the pressure. Lastly, the present results provide new insights into the structural and chemical stability of an important insensitive high explosive (IHE) crystal under well-defined HP-HT conditions.

36 MATERIALS SCIENCE↗

Comprehensive End-To-End Design of Novel High Energy Density Materials: III. Fused Heterocyclic Energetic Compounds

In contrast to synthesis and experimental characterization that are usually expensive and time-consuming, a coupled combination of ab initio and group additive theoretical methods represents a powerful alternative to computationally design and characterize new materials with targeted properties. Heterocyclic energetic materials are considered attractive candidates to replace conventional high explosives, such as PETN, RDX, HMX, and TNT. Heterocycles enable an intriguing opportunity for combinatorial design of new classes of energetic materials with targeted performance and sensitivity properties and offer ways for an accurate analysis of structure–property correlations of the materials. By using a computational strategy that links methods of quantum chemistry and semiempirical statistical analysis, we designed, modeled, and characterized a series of novel linear and fused explosive compounds. Proposed materials with best combinations of performance and sensitivity parameters can be now considered as potential candidates for synthesis and experimental characterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation into the Synthesis of Several TATB-Inspired Energetic Materials

The synthesis and characterization of several TATB-inspired energetic materials was explored. Despite their structural similarity to TATB, their behavior toward traditional synthetic pathways deviates significantly from that of the parent material. This change in expected reactivity is discussed and an efficient alternative synthesis is reported. The development of a simple and effective route to these TATB-type compounds opens the door to a potential new class of insensitive high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗