Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemistry - Chemical 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 109 records · Page 6

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005–1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2–3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of Indirect Quantification of 133mXe to Calibration of HPGe Detector

The INL Noble Gas Laboratory provides intercomparison samples for the noble gas analysis laboratories as part of the CTBTO PrepCom IMS. Xe-133m is one of the four relevant radionuclides in nuclear explosion monitoring. Without commercially available Xe-133m calibration standards laboratories must create and improve calibration methods. Improvements in calibration methods at the INL NGL benefit the CTBTO PrepCom through better certified values for Xe-133m intercomparison samples. Calibration of High Purity Germanium detectors for Xe-133m quantification is complicated by the coexistence of Xe-133 in samples under analysis. Xe-133 is typically produced in larger quantities, has higher gamma emission probabilities, and its gammas are detected more efficiently than Xe-133m. Xe-133m activity of samples can be indirectly inferred through the 133:133m activity ratio of a batch of material, and the Xe-133 counts in the assay of a small aliquot of the same material. This indirect quantification method can be leveraged to perform detector calibrations for quantification of Xe-133m. Calibrations can be performed by inferring the Xe-133m to certify the sample, and direct counting to determine detector efficiency. A comparison of method results will be shown.

133mXe↗

Radiolytic Gas Production from Aluminum Coupons (Alloy 1100 and 6061) in Helium Environments—Assessing the Extended Storage of Aluminum Clad Spent Nuclear Fuel

Corrosion of aluminium alloy clad nuclear fuel, during reactor operation and under subsequent wet storage conditions, promotes the formation of aluminium hydroxide and oxyhydroxide layers. These hydrated mineral phases and the chemisorbed and physisorbed waters on their surfaces are susceptible to radiation-induced processes that yield molecular hydrogen gas (H 2 ), which has the potential to complicate the long-term storage and disposal of aluminium clad nuclear fuel through flammable and explosive gas mixture formation, alloy embrittlement, and pressurization. Here, we present a systematic study of the radiolytic formation of H 2 from aluminium alloy 1100 (AA1100) and 6061 (AA6061) coupons in “dry” (~0% relative humidity) and “wet” (50% relative humidity) helium environments. Cobalt-60 gamma irradiation of both aluminium alloy types promoted the formation of H 2 , which increased linearly up to ~2 MGy, and afforded G-values of 1.1 ± 0.1 and 2.9 ± 0.1 for “dry” and “wet” AA1100, and 2.7 ± 0.1 and 1.7 ± 0.1 for “dry” and “wet” AA6061. The negative correlation of H 2 production with relative humidity for AA6061 is in stark contrast to AA1100 and is attributed to differences in the extent of corrosion and varying amounts of adsorbed water in the two alloys, as characterized using optical profilometry, scanning electron microscopy, Raman spectroscopy, and X-ray diffraction techniques.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reactive burn model calibration using high-throughput initiation experiments at sub-millimeter length scales

We report a first-of-its-kind model calibration was performed using Sandia National Laboratories’ high-throughput initiation (HTI) experiment for two types of vapor-deposited explosive films consisting of hexanitrostilbene (HNS) or pentaerythritol tetranitrate (PETN). These films exhibit prompt initiation, and they reach steady detonation at sub-millimeter length scales. Following prior work on HNS, we test the hypothesis of approximating these explosive films as fine-grained homogeneous solids with simple Arrhenius kinetics burn models. The model calibration process is described herein using a single-step as well as a two-step Arrhenius rate law, and it consists of systematic parameter sampling leading to a reduction in the model degrees of freedom. Multiple local minima are observed; results are given for seven different optimized parameter sets. Each model set is further evaluated in a two-dimensional simulation of the critical failure thickness for a sustained detonation. Overall, the two-step Arrhenius kinetics model captures the observed behavior for HNS; however, neither model produces a good fit to the PETN data. We hypothesize that the HTI results for PETN correspond to a heterogeneous response, owing to the smaller reaction zone of PETN compared to HNS (i.e., it does not homogenize the fine-grained hot spots as well). Future work should consider using the ignition and growth model for PETN, as well as other reactive burn models such as xHVRB, AWSD, PiSURF, and CREST.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disruption of electrostatic contacts in the HNH nuclease from a thermophilic Cas9 rewires allosteric motions and enhances high-temperature DNA cleavage

Allosteric signaling within multidomain proteins is a driver of communication between spatially distant functional sites. Understanding the mechanism of allosteric coupling in large multidomain proteins is the most promising route to achieving spatial and temporal control of the system. The recent explosion of CRISPR-Cas9 applications in molecular biology and medicine has created a need to understand how the atomic level protein dynamics of Cas9, which are the driving force of its allosteric crosstalk, influence its biophysical characteristics. In this study, we used a synergistic approach of nuclear magnetic resonance (NMR) and computation to pinpoint an allosteric hotspot in the HNH domain of the thermostable GeoCas9. We show that mutation of K597 to alanine disrupts a salt-bridge network, which in turn alters the structure, the timescale of allosteric motions, and the thermostability of the GeoHNH domain. This homologous lysine-to-alanine mutation in the extensively studied mesophilic S. pyogenes Cas9 similarly alters the dynamics of the SpHNH domain. We have previously demonstrated that the alteration of allostery via mutations is a source for the specificity enhancement of SpCas9 (eSpCas9). Hence, this may also be true in GeoCas9.

Chemistry↗

Effect of PBX Prill Size on the Compressive Properties of PBX 9502

Plastic-Bonded Explosive (PBX) prills are agglomerates a few millimeters in diameter composed of High Explosive crystals non-uniformly distributed in and surrounded by a polymer binder. When pressed at elevated temperatures at or above the binder melt temperature the polymer flows to uniformly coat the HE crystals to create a well-consolidated compaction. In reality it does so imperfectly, such that x-ray tomographic scans often look like a collection of prills mashed together. It is speculated that the larger the prill, the more non-uniform the binder distribution, but this has never been proven. Even if that were not so, the larger the prill the farther binder must flow toward the middle to homogenize during pressing. Thus, the degree of homogenization depends in part on prill size. The degree to which binder flows during pressing in turn affects the void distribution within pressed charges, which one suspects will affect shock sensitivity and material strength. In this paper we explore, for ~30 PBX 9502 formulation batches using the same TATB powder lot, how prill size effects the compressive mechanical properties.

36 MATERIALS SCIENCE↗

Energetic Copolymers From LLM-105 and Aliphatic Isocyanates

Energetic polymers typically feature fuel-rich backbones with pendant oxidizing explosophores, which are used to modify pressure and temperature characteristics in energetic formulations. However, these pendant explosophores generally increase sensitivity and decrease the thermal stability of the polymer in the condensed phase. Direct polymerization of insensitive high explosives (IHEs) bearing polymerizable functionalities, such as amines, provides a synthetic platform for deriving tunable energetic materials. Here, this work demonstrates the first direct polymerization of the IHE 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) with aliphatic isocyanate comonomers to yield an energetic copolyurea (PUa1) and an energetic copolyurea network (PUa2). 1 H and 13 C nuclear magnetic resonance (NMR) and Fourier transform infrared spectroscopy (FTIR) confirm copolyurea synthesis. PUa1 exhibited branching reactions commonly found in polyurea syntheses, and these side reactions were suppressed in PUa2 through use of a polyfunctional isocyanate. Differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) demonstrate the energetic IHE based polyurea preserves the energetic decomposition of the parent monomer LLM-105 in PUa1 (1203 J g −1 ) and PUa2 (687 J g −1 ), albeit with reduced thermal stability with peak decomposition temperatures of 235°C and 233°C, respectively. This work presents a new approach for generating energetic polymers from IHE cores with tunable energetic, thermal, and structural characteristics.

Chemistry↗

TIGER, A thermodynamics Equilibrium Tool for Explosives Update: Improved Solver

The thermodynamic equilibrium code known as TIGER has been in use since the 1970s and is designed to calculate the performance of energetic materials during detonation at high temperatures and pressures. The original TIGER code utilized a limited database consisting of 12 gaseous and 3 condensed constituents, which were composed of the elements carbon (C), hydrogen (H), nitrogen (N), oxygen (O), and aluminum (Al). In contrast, the more modern JCZS3 database features a significantly larger dataset, including 756 gaseous species, 189 positive and negative ions, and 496 condensed constituents derived from 62 different elements. While this expanded product species database allows for the exploration of a broader range of problems relevant to our laboratory, it also introduces challenges related to convergence, particularly when both liquid and solid species coexist under a vapor dome. In this work, we present an improved TIGER solution technique aimed at accurately determining the equilibrium state for systems that can produce a diverse array of products, including both solid and liquid condensed species. This report presents the status of the TIGER solver as of the end of FY25. Ongoing efforts are focused on enhancing the solver, and the report outlines several planned improvements. We anticipate that further developments will require additional documentation in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-Temperature-Stress Superposition of PBX 9502 Compressive Creep Data for Lifetime Predictions

The plastic-bonded explosive (PBX) 9502 is a viscoelastic, high-solids loaded polymer bound composite comprised of 95 weight% (wt%) TATB explosive crystals and 5 wt% FK-800 polymer binder. The uniaxial quasi-static mechanical properties have been studied and characterized extensively over 25+ years, including creep. Creep is the strain evolution over time as a result of an asymmetric load typically at temperatures above ambient. Creep measurements to failure can take excessively long periods of time depending on various factors including applied stress, temperature, and loading rate that render a single test to failure impossible. Using the principles of time-temperature-stress superposition (TTSSP), short-term creep tests at various applied stresses and temperatures can be used for long-term creep predictions. The TTSSP model shifts short-term test data using two factors, one for temperature and stress, creating a master curve at a reference test condition. The master curve can be shifted to predict the creep response at any desired test condition. PBX 9502 compression specimens were tested at three temperatures and various applied stresses. This work discusses the TTSSP analysis of PBX 9502 compressive creep data to create a single comprehensive master curve.

36 MATERIALS SCIENCE↗

Deacetylation and Mechanical Refining (DMR) and Deacetylation and Dilute Acid (DDA) Pretreatment of Corn Stover, Switchgrass, and a 50/50 Corn Stover/Switchgrass Blend

The most economical feedstock available to future lignocellulosic biorefineries may be a feedstock blend. Corn stover, switchgrass, and a 50:50 blend of corn stover and switchgrass were subjected to three different dilute alkali deacetylation conditions in a 90 L reactor prior to either mechanical refining (DMR) in a disc refiner followed by a Szego mill or dilute acid pretreatment (DDA) at five different conditions in a 4 L batch steam explosion reactor. The forty-five DDA and nine DMR slurries were subjected to high solids enzymatic hydrolysis at an enzyme loading of 15 mg protein per gram cellulose using Novozymes Cellic CTec3 (DDA) or a 4:1 mix of Novozymes Cellic CTec3:HTec 3 (DMR). The highest glucose yield was achieved using DMR processed corn stover from the high severity deacetylation condition (92 °C and 100 kg NaOH/ODMT) at 95.9% ± 0.8%, while the greatest monomeric glucose yields for switchgrass were between 75% to 80%. The glucose yields from saccharification of the 50:50 corn stover/switchgrass blend was approximately the average of the two individual feedstocks for both DDA and DMR processes. Statistical analysis showed that for a given feedstock type dilute alkali deacetylation severity had a significant impact (p < 0.05) on enzymatic hydrolysis yields, and for a given deacetylation severity the feedstock type had a significant effect on glucose conversion with the addition of switchgrass resulting in lower glucose yields. A blending strategy to maximize enzymatic hydrolysis sugar yields while reducing carbohydrate solubilization during the alkaline treatment step was proposed that allows deacetylation of feedstocks with different recalcitrance at optimal conditions for that feedstock, then blending prior to enzymatic hydrolysis.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗

Flammable Gas Generation and Control at the Idaho Completion Project Legacy Combustible Gas Generation - 20191

On April 11, 2018, four drums containing transuranic waste at the Idaho National Laboratory underwent over-pressurization, ejecting their lids and spreading radiological waste within a facility. An investigation has found that waste in the drums generated methane gas, which contributed to the event. Subsequent to the investigation, the potential for drums to have methane and other flammable analytes whose concentrations could approach or exceed the lower flammability limit (LFL) and the adequacy of the controls to prevent or mitigate a possible deflagration was evaluated. An extensive review of the historical records was performed to determine how many drums exceed the LFL. The historical record identified a small quantity of drums that exceeded the LFL for xylene, hydrogen, and methane. The primary codified applicable code or standard for handling drums is 29 CFR 1910.120 (j) 'Handling Drums and Containers.' It is used throughout the commercial and government sector. This code is integrated throughout the Fluor Idaho safety management programs and procedures to control the Standard Industrial Hazards (SIHs) associated with drum and container handling. Example requirements include: - Drums and containers are inspected, and their integrity assured prior to being moved. - Site operations are organized to minimize the amount of drum or container movement. - When there is a reasonable possibility of flammable atmospheres being present, material handling equipment and hand tools are of the type to prevent sources of ignition. Drums and containers under pressure, as evidenced by bulging or swelling, are not moved until the cause for excess pressure is determined and appropriate containment procedures have been implemented to protect employees from explosive relief of the drum. Drums which exceed the WIPP FGA limit have NCRs associated with them. Once a drum fails flammable gas sampling, an NCR is initiated. The NCR drives placement into NCR dense pack rows. A standing order drives placement of the FGA failures to a single planer segregation location. These drums are not allowed for further processing or shipment until the NCR has been cleared. Once in the single planer rows, follow-up FGA testing is performed by CCP. The NCR is dispositioned by follow-up sampling occurring after adequate time has been provided for the gas to diffuse through the filter assembly, or, in some cases, the drums are repackaged into additional drums to reduce the source term. To date, the ICP project has performed approximately 1.5 million drums moves without a deflagration event due to combustible gas generation. Flammable gas generation in a transuranic waste drum is not unique to the ICP but is common across the DOE complex. Based on the experience at Idaho, application of industry standards is sufficient to control the risk of drum deflagration due to drum movements. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Energetic Basis and Design of Enzyme Function Demonstrated Using GFP, an Excited-State Enzyme

We report the past decades have witnessed an explosion of de novo protein designs with a remarkable range of scaffolds. It remains challenging, however, to design catalytic functions that are competitive with naturally occurring counterparts as well as biomimetic or nonbiological catalysts. Although directed evolution often offers efficient solutions, the fitness landscape remains opaque. Green fluorescent protein (GFP), which has revolutionized biological imaging and assays, is one of the most redesigned proteins. While not an enzyme in the conventional sense, GFPs feature competing excited-state decay pathways with the same steric and electrostatic origins as conventional ground-state catalysts, and they exert exquisite control over multiple reaction outcomes through the same principles. Thus, GFP is an “excited-state enzyme”. Herein we show that rationally designed mutants and hybrids that contain environmental mutations and substituted chromophores provide the basis for a quantitative model and prediction that describes the influence of sterics and electrostatics on excited-state catalysis of GFPs. As both perturbations can selectively bias photoisomerization pathways, GFPs with fluorescence quantum yields (FQYs) and photoswitching characteristics tailored for specific applications could be predicted and then demonstrated. The underlying energetic landscape, readily accessible via spectroscopy for GFPs, offers an important missing link in the design of protein function that is generalizable to catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrogen Transport in a Model 9979 Shipping Package with Inner Convenience Cans

Radiolytic hydrogen production and accumulation inside containment packages is a concern at any facility responsible for their packaging, storage, transportation, and/or disposal. When hydrogen gas accumulates to concentrations above the Lower Flammability Limit (LFL) which is 4% or 40,000 ppm in air, the possibility of a deflagration or explosion increases. This concern persists over the course of the package lifetime which is unlimited when disposed of by burial or in permanent repositories. Here, we report on a numerical model used to predict the concentration of hydrogen within each layer of a Model 9979 package containing a convenience can assembly. Simulations show the hydrogen concentration to always be highest in the inner convenience can containing the radioactive source. When the radioactive source is within the Los Alamos National Laboratory (LANL) Packaging Limits, the hydrogen concentration is shown to remain well below the LFL at all times including packaging, storage, transportation, and disposal. A hydrogen transport model is presented for a Model 9979 package system containing a nested arrangement of convenience cans, which are tin oxide coated steel cans of various sizes with a slip-lid assembly. The inner convenience can contains the radioactive source material along with an unknown quantity of incidental water acquired from humid air or processing. While visible organic materials such as paper and plastics were purposely excluded from the inner can, it is not possible to claim the wastes are entirely organic free. The inner convenience can is tape sealed and placed into a plastic bag which is horsetail closed (i.e., twisted and taped). The bagged can is placed into an outer convenience can that is also tape sealed. The can assembly is then placed into the 30 gallon drum and subsequently placed inside the 55 gallon drum in the 9979 package. Here we assume hydrogen gas is produced in the inner convenience can from alpha radiolysis of water at a rate dependent on the quantity of uranium isotopes and water present. The hydrogen transport model was used to calculate hydrogen accumulations within the package’s five layers at different times and conditions. These simulations serve two purposes; (i) to build confidence in the model by comparing predicted values to measured values, and (ii) to check the steady state hydrogen concentrations that are approached at long times in the package’s lifetime. Model simulations were compared to gas samples taken from the 30 gallon drum after storage at LANL’s Chemistry and Metallurgy Research (CMR) building for around 500 days. Hydrogen concentration calculations over much longer periods (i.e., more than 270 years) included extreme storage durations, transportation at extreme cold temperatures, and disposal of packages assuming different average temperatures.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗