Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HMX”

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 37 records · Page 2

Modeling The Effects of Shock Pressure and Pore Morphology on Hot Spot Mechanisms in HMX

Here, we investigate the effects of shock pressure and pore morphology on the formation and growth of hot spots in HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Both non-reactive and reactive ALE3D simulations are used in these studies. Our non-reactive simulations show a viscous-dominated pore collapse mode at lower shock pressures (2–10 GPa) with shear band formation and a hydrodynamic-dominated mode at higher shock pressures (20-40 GPa) due to bulk melting. When normalized by bulk shock heating, viscous-dominated pore collapse modes are more efficient at generating hot spots. Pore morphology influences the post-collapse temperature distributions and reaction rate for a fixed pore area and shock pressure. We find that multiple surface pores at the binder-grain interface tend to react the fastest. Due to their upstream location in the HMX grain, the surface pores collapse sooner than interior pores; thus, the extent of reaction will generally favor these morphologies because they have more time to grow. In general, multiple smaller hot spots tend to react faster than a single larger hot spot because they accelerate one another's burning. The rank order of morphology effects, however, is not the same for non-reactive and reactive simulations. For example, while multiple surface pores produce the highest reaction rates they do not produce the highest (non-reactive) hot spot temperatures. In conclusion, our numerical studies provide insights on hot spot mechanisms in lieu of direct measurements and can be used to develop advanced shock initiation models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray Phase Contrast Imaging of the Impact of Multiple HMX Particles in a Polymeric Matrix

The initiation of high explosives (HEs) under shock loading lacks a comprehensive understanding: particularly at the particle scale. One common explanation is the hot spot theory, which suggests that energy in the material resulting from the impact event is localized in a small area causing an increase in temperature that can lead to ignition. This study focuses on the response of HMX particles (a common HE) within a polymer matrix (Sylgard-184 (R)), a simplified example of a polymer-bound explosive (PBX). These PBXs consist of multiple HMX particles in a single polymer-bound sample. A light gas gun was used to load the samples at impact velocities above 400 m/s. The impact events were visualized using X-ray phase-contrast imaging (PCI) allowing real-time observation of the impact event. The experiment used two different types of samples (multi-particle and two crystals) and found evidence of cracking and debonding in both sample types. Additionally, it was found that the multiple particle samples showed similar evidence of damage at lower velocities than that of single particle samples. This is an expected result as the multiple particles add additional interfaces for stress concentration and frictional heating.

42 ENGINEERING↗

Phase Contrast X‐Ray Imaging of the Collapse of an Engineered Void in Single‐Crystal HMX

Abstract Imaging the collapse of a single void that creates a hot spot initiation site in an otherwise defect‐free explosive is challenging given the spatial and temporal scales involved in explosive systems. This work presents our attempt to examine a single hot spot mode (void collapse) in single‐crystal octahydro‐l,3,5,7‐tetranitro‐l,3,5,7‐tetrazocine (HMX) embedded in Sylgard. The hot spot heating mechanisms involved with pore collapse include adiabatic heating, jetting, and viscoplastic dissipation. Quantifying the dynamics of a pore collapse is a crucial step to understanding which mechanisms dominate during ignition events. Our experiments were conducted with a single‐stage, light‐gas gun at Argonne National Laboratory's Advanced Photon Source, applying the phase contrast imaging technique while collecting high‐speed video. The details of HMX single crystal production, defect (pore) engineering, and sample construction, along with experimental results are presented here. These results demonstrate that detailed collapse dynamics can be obtained from homogeneous, single‐crystal explosives with this approach. Qualitative comparisons are made with simulation data which show good agreement in the transition between a quasi‐symmetric pore collapse and an asymmetric collapse with jetting across the pore as measured with normalized pore area and pore circularity.

42 ENGINEERING↗

Cohesive finite element modeling of the delamination of HTPB binder and HMX crystals under tensile loading

Accurately modeling the mechanical behavior of the polymer binders and the degradation of interfaces between binder and crystal is important to science-based understanding of the macro-scale response of polymer bonded explosives. The paper presents a description of relatively a simple bi-crystal HMX-HTPB specimen and associated tensile loading experiment including computed tomography imaging, the pertinent constitutive theory, and details of numerical simulations used to infer the behavior of the material during the delamination process. Within this work, mechanical testing and direct numerical simulation of this relatively simple bi-crystal system enabled reasonable isolation of binder-crystal interface delamination, in which the effects of the complicated thermomechanical response of explosive crystals were minimized. Cohesive finite element modeling of the degradation and delamination of the interface between a modified HTPB binder and HMX crystals was used to reproduce observed results from tensile loading experiments on bi-crystal specimens. Several comparisons are made with experimental measurements in order to identify appropriate constitutive behavior of the binder and appropriate parameters for the cohesive traction-separation behavior of the crystal-binder interface. This research demonstrates the utility of directly modeling the delamination between binder and crystal within crystal-binder-crystal tensile specimen towards characterizing the behavior of these interfaces in a manner amenable to larger scale simulation of polycrystalline PBX materials. One critical aspect of this approach is micro computed tomography imaging conducted during the experiments, which enabled comparison of delamination patterns between the direct numerical simulation and actual specimen. In addition to optimizing the cohesive interface parameters, one important finding from this investigation is that understanding and representing the strain-hardening plasticity of HTPB binder is important within the context of using a cohesive traction-separation model for the delamination of a crystal-binder system.

42 ENGINEERING↗

Importance of microstructural features in mechanical response of cast-cured HMX formulations

We report that cast-cure formulations of plastic-bonded explosives (PBX) present some advantages in manufacturing and application compared to traditional pressed slurry formulated materials. For example, these formulations can be cast to specific shapes and then cured in place, avoiding the need for machining precision parts. However, the microstructure of these materials can be greatly affected by the specific manufacturing process. Here, we evaluate the effect of minor changes to the formulation and manufacturing process on several cyclotetramethylene-tetranitramine (HMX) PBXs. The binders were based on hydroxyl-terminated polybutadiene (HTPB) and cast-cured using diphenylmethane diisocyanate (isonate) as a curing agent. We examined the materials using X-ray computed tomography (CT) imaging and uniaxial compression testing. The isonate content was found to significantly affect the modulus and strength of the binder. The presence of significant void content could be controlled by adding a centrifuging step during the curing process, but the resulting effect on mechanical properties was relatively minor. Finally, mesoscale simulations showed that differences in the mechanical strength of the binder were not sufficient to describe the differences observed in mechanical testing, indicating that the HMX-binder adhesion was also being changed by the manufacturing process.

36 MATERIALS SCIENCE↗

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↗

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↗

Numerical parameter optimizations of the Ignition and Growth model for a HMX plastic bonded explosive

We present here a novel method for optimizing parameters of the Ignition and Growth (I&G) reactive flow model for high explosives. The I&G model predicts the shock initiation response of explosives subjected to dynamic loading. However, calibrating the model is a time-consuming process because it has several parameters, typically requires multiple data sets, and there are no direct correlations between parameters and data. In this study, we couple a scalable optimization algorithm to simulations of shock initiation experiments in the multi-physics code ALE3D. We develop four I&G model parameter sets for the HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane)-based explosive LX-07 [90 wt. % HMX, 10 wt. % Viton A] based on minimizing the difference between calculations and measurements at the embedded pressure gauges in 1D gun shock experiments. Our study shows that a cost function based on both shock time of arrival and pressure pulse shape demonstrates the best agreement with experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ignition thresholds of aluminized HMX-based polymer-bonded explosives

Here, the ignition of aluminized HMX-based polymer-bonded explosives (PBXs) under shock loading is studied via mesoscale simulations. The conditions analyzed concern loading pulses of 20 nanoseconds to 0.8 microseconds in duration and impact piston velocities on the order of 400-1000 m/s or loading stresses on the order of 3-14 GPa. The sets of samples studied have stochastically similar microstructures consisting of a bimodal distribution of HMX grains, an Estane binder, and aluminum particles 50-100 µm in diameter. The computational model accounts for constituent elasto-viscoplasticity, viscoelasticity, bulk compressibility, fracture, interfacial debonding, internal contact, bulk and frictional heating, and heat conduction. The analysis focuses on the development of hotspots under different material settings and loading conditions. In particular, the ignition thresholds in the forms of the James relation and the Walker-Wasley relation and the corresponding ignition probability are calculated and expressed as functions of the aluminum volume fraction for the PBXs analyzed. Lastly, it is found that the addition of aluminum raises the ignition thresholds, causing the materials to be less sensitive. Dissipation and heating mechanism changes responsible for this trend are delineated.

36 MATERIALS SCIENCE↗

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CISME report on experiments with HMX powder, PBX 9501 prills, and heated PBX 9501

The ability to predict post-ignition explosive response is relevant to explosives safety, particularly weapons assembly/disassembly operations that must operate under a sufficiently conservative safety basis to meet risk mitigation criteria for inadvertent nuclear detonation (IND). In particular, hazard determinations hinge on whether the PBX 9501 main charge can ultimately attain detonation from an initial, non-shock, ignition stimulus. The CISME experiment uses a spherical explosive charge, which is thermally ignited in the center using a fiber-delivered infrared laser pulse. This ignition mechanism - i.e. a localized thermal hotspot - is representative of the ignition-causing insults typical in accident scenarios, including low-velocity impact, drop, skid, and punch events. In fiscal year 2018 CISME experiments were conducted with pristine PBX 9501 spheres, confined and unconfined. Sub-detonative, low-violence reactions were observed. Here we report experiments conducted since 2018. The explosive charge was varied: experiments were conducted with room temperature pour-density HMX powder and prills of PBX 9501, and heated, thermally damaged PBX 9501. Experiments were conducted in a confined configuration, with two 4 mm-thick stainless-steel hemispherical shells bolted around the sphere. Detonation, via DDT (deflagration-to-detonation transition) was observed in the 3-inch diameter powdered HMX and the 6-inch diameter heated and thermally-damaged PBX 9501 tests. No transition to DDT was observed in smaller diameters for those configurations, and no DDT was observed for pour-density prills of PBX 9501 of any diameter. The original test series demonstrated that the response of pristine, undamaged PBX 9501 to a central ignition is relatively benign. This new test series shows that increased porosity, lack of polymer binders and thermal damage can greatly increase the violence of the post-ignition response, up to and including detonation. These results bridge the explosive response of PBX 9501 from sub-detonative through detonation as a function of diameter and provide valuable test cases for modeling efforts that strive to capture the DDT threshold. By exploring explosive response in a more representative 3D geometry, we can also benchmark these new results to more plentiful DDT datasets generated from testing in 1D tube confinement (the historical test configuration)–the configuration upon which almost all DDT models were developed.

42 ENGINEERING↗

Shock and Reaction in Granular Bed of HMX-Aluminum Powders

Ignition and material response properties of aluminumized HMX heterogeneous explosive mixtures were explored in a series of planar impact experiments performed over multiple years. This work expands on previous work studying material response to impact in single-component HMX granular materials. The addition of nanometric aluminum is shown to affect the ignition sensitivity and growth to reaction from impact. The gas gun test results are presented here varying parameters of particle size, shock strength, and aluminum mass fraction.

36 MATERIALS SCIENCE↗

The Combustion of HMX

The burn rate of HMX was measured at high pressures (p more than 1000 psi). The self deflagration rate of HMX was determined from 1 atmosphere to 50,000 psi. The burning rate shows no significant slope breaks.

Boggs, T. L.↗

Critical analysis of nitramine decomposition data: Activation energies and frequency factors for HMX and RDX decomposition

A summary of a literature review on thermal decomposition of HMX and RDX is presented. The decomposition apparently fits first order kinetics. Recommended values for Arrhenius parameters for HMX and RDX decomposition in the gaseous and liquid phases and for decomposition of RDX in solution in TNT are given. The apparent importance of autocatalysis is pointed out, as are some possible complications that may be encountered in interpreting extending or extrapolating kinetic data for these compounds from measurements carried out below their melting points to the higher temperatures and pressure characteristic of combustion.

Schroeder, M. A.↗

Investigating deformation and mesoscale void creation in HMX based composites using tomography based grain scale FEM

The microstructure of plastic bonded explosives (PBXs) significantly affects their macroscale mechanical characteristics. Imaging and modeling of the mesoscale constituents allows for a detailed examination of the deformation of mechanically loaded PBXs. In this study, a pair of explosive HMX crystals were bound with an HTPB based polymer binder and imaged using micro Computed Tomography ($µ$CT). Delamination of the polymer binder from the crystals was induced via macro-scale application of a tensile load. Cohesive parameters for simulation of the crystal/binder interface were determined by comparing numerical and experimental results of the delamination of this bi-crystal system. Simulations showed macro-scale plastic behavior is accommodated through both crystal-binder delamination and plastic deformation of the binder itself. Imaging and simulation of this bi-crystal system demonstrated the complex interplay of material constitutive behavior and grain-scale geometry. Successful parameterization of this relatively simple mechanical system is a step toward a better understanding void nucleation in similarly composed polycrystalline composites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dependence of the Elastic Stiffness Tensors of PETN, α‐RDX, γ‐RDX, ϵ‐RDX, ϵ‐CL‐20, DAAF, FOX‐7, and β‐HMX on Hydrostatic Compression

Abstract The dependence of the components of the elastic stiffness tensors (or elastic constants) of the organic explosives PETN, RDX, CL‐20, DAAF, FOX‐7, and HMX on hydrostatic pressure up to 10 GPa have been computed using dispersion‐corrected density functional theory. We report the evolution of lattice parameters and the non‐zero stiffnesses for the tetragonal, orthorhombic, and monoclinic crystal symmetries. Linear and quadratic dependencies of the components of the elastic stiffness tensors on volumetric compression and hydrostatic pressure are tabulated for use in single crystal plasticity models.

36 MATERIALS SCIENCE↗

Thermal Gradient Effects on Local Hotspot Ignition in 1,3,5,7‐Tetranitro‐1,3,5,7‐tetrazocane (HMX)

Understanding hotspot ignition, growth, and criticality, as well as the timescales of each, is crucial for parameterizing mesoscale and continuum‐level models that rely on a statistical understanding of hotspots. However, these models often consider hotspots to have uniform temperatures or for hotspots of a given temperature and size to always behave the same. Therefore, using molecular dynamics simulations, we assess the influence of thermal distribution effects on hotspot local ignition and time to ignition in 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane (HMX) nanoscale hotspots. Finally, by assessing hotspots with a gradient driven from an initial two‐temperature core–shell setup, we show that small increases in the core temperature under a constant average temperature can lead to order‐of‐magnitude effects on reaction and local ignition timescales.

36 MATERIALS SCIENCE↗