Charge Characterization for Blast Effects Testing using Explosively Driven Blast Tubes
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Sensitivity to initiation is an important characteristic of high explosives. Many chemical and physical properties are fine-tuned to deliver an energetic material that behaves according to the user’s specifications exactly. In some cases, it is advantageous to study the effects that non energetic additives have on the performance of an explosive material. Following on the promising results obtained in other experiments, a study is conducted on the effects of hollow glass micro-balloons (bubbles, spheres), or HGMB, on a material’s sensitivity. Threshold test data is compared between several formulations of explosives. The control formulation is an RDX-based explosive with a non-energetic binder, and the test formulations are identical in composition, but with the addition of HGMB, up to 6% by volume. Safety testing of the new HGMB-doped explosive formulations is conducted to ensure they can be produced and handled safely prior to intentional detonation. This experiment utilizes an attenuation card of varying thicknesses between the initiator and the test article to attenuate the shock imparted, and a witness plate is analyzed for each test article to determine a GO or NO GO result. The data collected and the conclusions of this study affect how additives like HGMB are utilized in future explosive sensitivity threshold experiments.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The prevailing method to produce plastic‐bonded explosive (PBX) molding powder, or “prills”, is a complex, multiphase, and bespoke process that was developed by the high explosives (HEs) manufacturing industry several decades ago. This work demonstrates the utility of a simpler, widely‐used mechanical process—extrusion‐spheronization—to produce PBX prills. We begin by detailing precautions taken to enable safe remote operation of extrusion‐spheronization equipment with HE. We then perform a study investigating the effect of lacquer solvent composition on the particle packing, pressed density, and compressive strength properties of a 95 wt.% TATB/5 wt.% polymer binder formulation akin to PBX 9502. It was found that increased composition of low vapor pressure solvents caused prolonged retention of the solvent, resulting in tackier materials that would agglomerate and form larger prills. The larger prills also led to lower poured density and tapped density of HE prills and compressive strength of pressed PBX articles. The samples prepared with a 75% propyl acetate/25% butyl acetate lacquer solvent composition exhibited the highest compressive strength. However, it is believed that the prill packing and compressive strength properties are primarily driven by the prill size rather than the chemical composition of the lacquer itself. Extrusion‐spheronization remains a promising method to reliably and repeatably produce HE prills that is less sensitive to feedstock or process variation than traditional methods.
Hot spots are spatial regions of intense energy localization that govern initiation of secondary high explosives. Studies that characterize or compare simulated hot spots are frequently either qualitatively descriptive or resort to quantitative distribution functions that neglect stochastic variations and spatial correlations—effects that are also neglected in common comparison tests like the Kolmogorov–Smirnov test. To this end, we develop an image distinguishability analysis (IDA) test based on principal component (PC) analysis that makes pixel-by-pixel comparisons between small, for example, O(<10), image data sets. The IDA test makes comparisons through a generalized distance metric in the PC space and a test statistic that is derived to calculate mathematical equation-values. Here, we derive a statistical distribution and criticality criterion to determine whether images are distinguishable from established baselines. We apply the IDA test on images generated from molecular dynamics simulations of hot spots from pore collapse in TATB to assess scale invariance in the complex patterns of hot spots that form in a representative high explosive crystal. The IDA test shows that TATB hot spot spatial temperature fields and their derived temperature histograms exhibit scale-invariant features over specific intervals of shock orientation, strength, and initial pore diameter. However, the IDA test also shows that qualitatively different conclusions regarding invariance can be reached depending on whether the hot spot is treated as a spatially correlated field as opposed to a distribution function that lacks spatial information.
A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.
In this numerical study, we investigate the effects of layered high explosive (HE) charge design on Richtmyer–Meshkov instability (RMI) in metal plates with sinusoidal surface features. The detonation wave from the HE induces a shock in the metal target that subsequently interacts with the surface features; this results in vortex formation and ultimately RMI. Here, we seek to modify RMI by altering the detonation wave characteristics. The modification is investigated in a twofold manner: first, by varying the initial design of the unconfined surface of the target and second, by varying the charge design and composition. Within a limited scope of this design space, a wide variety of behaviors related to RMI growth are observed. Mechanistic actions, including exaggerated front curvature, behind these modifications are proposed. Charge designs, which modify RMI the most for a select target design, are then presented.
The shock-induced reaction of high explosives to gaseous and solid reaction products is a rapid, complex process. Understanding solid reaction product structure and formation kinetics is essential in determining the high-pressure equation of state of these multicomponent systems. We use the National Ignition Facility (NIF) to shock compress ~500-µm thick pressed powder high-explosive TATB samples to 70-135 GPa and collect in situ structural data on detonation byproducts using the TARDIS X-ray diffraction diagnostic. Velocimetry is used to record the transmitted compression wave profile, which is well described by Cheetah hydrocode simulations coupled with a reactive flow model, providing strong evidence of reaction in the TATB sample. While an unambiguous determination of the product phases was not possible owing to the low signal-to-noise quality of the diffraction signal, X-ray diffraction data of the product phases formed within the first 50 ns of the reaction is most consistent with a mixture of amorphous products and crystalline hexagonal diamond.
A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.
Mechanical properties are of interest for many plastic‑bonded explosive (PBX) materials with tensile properties being of particular interest. Direct tensile measurements using dogbone-shaped samples are considered the gold standard, but they are fairly large, making testing more costly and less desirable from a safety perspective. We investigated whether the measured tensile strength depends on the dogbone specimen size, which to our knowledge, has not been reported in the literature for PBX materials. Understanding this should inform the feasibility of employing smaller samples and how sample size should be considered when comparing PBX dogbone values in the literature. The TATB-based PBX dogbone sample size was varied by (a) scaling all dimensions proportionally and (b) varying only the length of the samples. It was observed that the measured tensile peak stress (strength) was a function of the sample size, and was more dependent on the diameter (cross-sectional area) than the length of the samples. Since peak stress is calculated as peak force normalized to the diameter of the sample, one might not expect an explicit diameter dependence for the peak stress. Therefore, these results suggest there may be an additional strengthening effect as the sample diameter is increased.
Carbon particulates generated during detonation depend upon high explosive type, composition, and detonation conditions. Although explosive composition greatly affects particulates, the focus of this work is on how detonation geometries that induce much higher temperatures and pressures in the high explosive lead to differing particulate morphologies. In this study, two geometries were used: Detonations were initiated in Composition B cylinders at one end in conventional detonations and initiated at both ends to produce colliding detonations. Each of these detonations was observed on the sub-μs timescale using fast radiography capturing images of the front moving through the cylinder, and colliding detonation fronts in real-time. These imaging experiments were complemented with time-resolved small-angle x-ray scattering (SAXS) experiments that were able to observe and determine the varying condensed carbon morphologies at different locations and times in each detonation. The detonations could be timed in such a way that the spatial and temporal dependence of the carbon morphology could be superimposed onto radiography images collected at the same point in time. The complementary approach is able to show that the carbon condensates are much larger when formed in the elevated temperature and pressure conditions near the location of colliding detonation fronts. Thermochemical modeling suggests that these larger particulates form either in the diamond phase or on the liquidus line of the carbon phase diagram. The increase in size observed by SAXS may correlate well with the increased residence time deeply in the diamond phase. Finally, these particulates can be described as nano-sized phases with some surface texture or otherwise near-surface intra-particle heterogeneity.
ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.
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.
Various small-scale experiments were performed to provide data for developing a model to predict the thermal response of LLM-105 over a wide range of conditions. The thermal decomposition of LLM-105 was studied as a function of sample mass, confinement of volatile products, and preheating time in both isothermal and ramped heating experiments. The thermal decomposition of LLM-105 is a two-step process, as shown by the two exothermic peaks in the heat flow profiles, which were fitted to two nth-order autocatalytic reaction models with a similar activation energy of ∼289 kJ/mol. The magnitude and shape of these peaks varied with sample mass and confinement. Increasing sample mass enhanced the second exotherm with respect to the first one, while increasing the level of confinement promoted a transition from a sublimation-dominated regime towards thermal decomposition. The effect of LLM-105 particle size on the rate of weight loss was evident for open-pan experiments, where bigger particles sublimed at lower temperatures than smaller particles. Thermal response and solid residue composition of LLM-105 samples were analyzed following preheating for different durations. Longer preheating times caused a shift of the second exotherm to lower temperatures and a decrease in the reaction enthalpy, confirming that LLM-105 decay is a consecutive reaction mechanism, probably autocatalytic. In conclusion, the kinetic model derived from ramped experiments was validated against the measured LLM-105 fraction remaining and the enthalpy remaining of the solid residue as a function of preheating times and showed good agreement.
Explore the source record for details and available documents.
A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.
We demonstrate a strongly thresholded response in cyclotrimethylene trinitramine (RDX) when it is cylindrically shocked using a novel waveguide geometry. Using ultrafast single-shot multi-frame imaging, we demonstrate that <100 μm diameter single crystals of RDX embedded in a polymer host deform along preferential planes for >100 ns after the shock first arrives in the crystal. We use in situ imaging and time-resolved photoemission to demonstrate that short-lived chemistry occurs with complex deformation pathways. Using scanning electron microscopy and ultra-small-angle X-ray scattering, we demonstrate that the shock-induced dynamics leave behind porous crystals, with pore shapes and sizes that change significantly with shock pressure. A threshold pressure of ~12 GPa at the center of convergence separated the single-mode planar crystal deformations from the chemistry-coupled multi-plane dynamics at higher pressures. Our observations indicate preferential directions for deformation in our cylindrically shocked system, despite the applied stress along many different crystallographic planes.