Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. We found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.
Simultaneous data of the quasi-static compaction and electrical conductivity of porous, binary powder mixtures have been collected as a function of bulk density. The powder mixtures consist of a metal conductor, either titanium or iron, an insulator, and pores filled with ambient air. The data show a dependency of the conductivity in terms of relative bulk density and metal volume fraction on conductor type and conductor particle characteristics of size and shape. Finite element models using particle domains generated by discrete element method are used to simulate the bulk conductivity near its threshold while the general effective media equation is used to model the conductivity across the compression regime.
This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.
Validated models of melt cast explosives exposed to accidental fires are essential for safety analysis. In the current work, we provide several experiments that can be used to develop and validate cookoff models of melt cast explosives such as Comp-B3 composed of 60:40 wt% RDX:TNT. We present several vented and sealed experiments from 2.5 mg to 4.2 kg of Comp-B3 in several configurations. We measured pressure, spatial temperature, and ignition time. Some experiments included borescope images obtained during both vented and sealed decomposition. We observed the TNT melt, the suspension of RDX particles in the melt, bubble formation caused by RDX decomposition, and bubble-induced mixing of the suspension. The RDX suspension did not completely dissolve, even as temperatures approached ignition. Our results contrast with published measurements of RDX solubility in hot TNT that suggest RDX would be completely dissolved at these high temperatures. These different observations are attributed to sample purity. We did not observe significant movement of the two-phase mixture until decomposition gases formed bubbles. Bubble generation was inhibited in our sealed experiments and suppressed mixing.
In previous studies, we found that the nitroplasticizer in the HMX-based explosive PBX 9501 played a crucial role in cookoff, especially when predicting response in larger systems. We have recently completed experiments with a similar explosive, LX-14, that has a relatively nonreactive binder. We expected the ignition times for LX-14 to be longer than PBX 9501 since PBX 9501 has a more reactive binder. However, our experiments show the opposite trend. This paradox can be explained by retention of reactive gases within the interior of LX-14 by the higher strength binder resulting in faster ignition times. In contrast, the binder in PBX 9501 melts at low temperatures and does not retain decomposition gases as well as the LX- 14 binder. Retention of reactive gases in LX-14 may also explain the more violent response in oblique impact tests when compared to PBX 9501.
In this study, we have made reasonable cookoff predictions of large-scale explosive systems by using pressure-dependent kinetics determined from small-scale experiments. Scale-up is determined by properly accounting for pressure generated from gaseous decomposition products and the volume that these reactive gases occupy, e.g. trapped within the explosive, the system, or vented. The pressure effect on the decomposition rates has been determined for different explosives by using both vented and sealed experiments at low densities. Low-density explosives are usually permeable to decomposition gases and can be used in both vented and sealed configurations to determine pressure-dependent reaction rates. In contrast, explosives that are near the theoretical maximum density (TMD) are not as permeable to decomposition gases, and pressure-dependent kinetics are difficult to determine. Ignition in explosives at high densities can be predicted by using pressure-dependent rates determined from the low-density experiments as long as gas volume changes associated with bulk thermal expansion are also considered. In the current work, cookoff of the plastic-bonded explosives PBX 9501 and PBX 9502 is reviewed and new experimental work on LX-14 is presented. Reactive gases are formed inside these heated explosives causing large internal pressures. The pressure is released differently for each of these explosives. For PBX 9501, permeability is increased and internal pressure is relieved as the nitroplasticizer melts and decomposes. Internal pressure in PBX 9502 is relieved as the material is damaged by cracks and spalling. For LX-14, internal pressure is not relieved until the explosive thermally ignites. The current paper is an extension of work presented at the 26th ICDERS symposium [1].
It has long been observed that oxidation processes in metals tend to follow a parabolic rate law associated with the growth of a surface oxide layer. Here we observe that for certain titanium powders, the expected parabolic law (∝ t1/2) is recovered, yet for others, the exponent differs significantly. One explanation for this non-parabolic, anomalous diffusion arises from fractal geometry. Theo retical considerations indicate that the time response of diffusion-limited processes in an object closely follow a power-law in time (tn) with n=(E−D)/2, where E is the object's Euclidean dimension and D is its boundary's Hausdorff dimension. Non-integer, (fractal) values of D will result in n≠1/2. Finite element simulations of several canonical fractal objects were performed to verify the application of this theory; the results matched the theory well. Two different types of titanium powder were tested in isothermal thermogravimetric tests under dilute oxygen. Time-dependent mass uptake data were fit with power-law forms and the associated exponents were used to determine an equivalent fractal dimension. One Ti powder type has an implied surface dimension of ca. 2.3 to 2.5, suggesting fractal geometry may be operative. The other has a dimension near 2.0, indicating it behaves like traditional material.