Exploding bridgewire detonators (EBWs) containing pentaerythritol tetranitrate (PETN) exposed to high temperatures may not function following discharge of the design electrical firing signal from a charged capacitor. Knowing functionality of these arbitrarily facing EBWs is crucial when making safety assessments of detonators in accidental fires. Orientation effects are only significant when the PETN is partially melted. The melting temperature can be measured with a differential scanning calorimeter. Nonmelting EBWs will be fully functional provided the detonator never exceeds 406 K (133 °C) for at least 1 h. Conversely, EBWs will not be functional once the average input pellet temperature exceeds 414 K (141 °C) for a least 1 min which is long enough to cause the PETN input pellet to completely melt. Functionality of the EBWs at temperatures between 406 and 414 K will depend on orientation and can be predicted using a stratification model for downward facing detonators but is more complex for arbitrary orientations. A conservative rule of thumb would be to assume that the EBWs are fully functional unless the PETN input pellet has completely melted.
Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. We present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.
Determining the thermal response of energetic materials at high densities can be difficult when pressure dependent reactions occur within the interior of the material. At high temperatures, reactive components such as hexahydro-l,3,5-tri-nitro-l,3,5-triazine (RDX), ammonium perchlorate (AP), and hydroxyl-terminated polybutadiene (HTPB) decompose and interact. The decomposition products accumulate near defects where internal pressure ultimately causes mechanical damage with closed pores transitioning into open pores. Gases are no longer confined locally; instead, they freely migrate between open pores and ultimately escape into the surrounding headspace or vent. Recently we have developed a universal cookoff model (UCM) coupled to a micromechanics pressurization (MMP) model to address pressure-dependent reactions that occur within the interior of explosives. Parameters for the UCM/MMP model are presented for an explosive and two propellants that contain similar portions of both aluminum (Al) and a binder. The explosive contains RDX and the propellants contain AP with no RDX. One of the propellants contains small amounts of curing catalysts and a burn modifier whereas the other propellant does not. We found that the cookoff behavior of the two propellants behave similarly leading and conclude that small amounts of catalysts or burn modifiers do not influence cookoff behavior appreciably. Kinetic parameters for the UCM/MMP models were obtained from the Sandia Instrumented Thermal Ignition (SITI) experiment. Validation is done with data from other laboratories.
Exploding bridgewire detonators (EBWs) containing pentaerythritol tetranitrate (PETN) exposed to high temperatures may not function following discharge of the design electrical firing signal from a charged capacitor. Knowing functionality of these arbitrarily facing EBWs is crucial when making safety assessments of detonators in accidental fires. Orientation effects are only significant when the PETN is partially melted. Here, the melting temperature can be measured with a differential scanning calorimeter. Nonmelting EBWs will be fully functional provided the detonator never exceeds 406 K (133 °C) for at least 1 h. Conversely, EBWs will not be functional once the average input pellet temperature exceeds 414 K (141 °C) for a least 1 min which is long enough to cause the PETN input pellet to completely melt. Functionality of the EBWs at temperatures between 406 and 414 K will depend on orientation and can be predicted using a stratification model for downward facing detonators but is more complex for arbitrary orientations. A conservative rule of thumb would be to assume that the EBWs are fully functional unless the PETN input pellet has completely melted.
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.
Understanding titanium particle combustion processes is critical not only for characterizing existing pyrotechnic systems but also for creating new igniter designs. In order to characterize titanium particle combustion processes, morphologies, and temperatures, simultaneous spatially-resolved electric field holography and imaging pyrometry techniques were used to capture post-ignition data at up to 7 kHz. Due to the phase and thermal distortions present in the combustion cloud, traditional digital in-line holography techniques fail to capture accurate data. In this work, electric field holography techniques are used in order to cancel distortions and capture the three-dimensional spatial locations and diameters of the particles. In order to estimate the projected surface temperatures of the titanium particles, an imaging pyrometry method that ratios emission at 750 and 850 nm is utilized. Using these diagnostics, joint statistics are collected for particle size, morphology, velocity, and temperature. Results show that, early in the combustion process, the titanium particles are primarily oxidized by potassium perchlorate inside the igniter cup, resulting in projected surface temperatures near 3000 K. Later in the process, the particles interact with ambient air, resulting in lower surface temperatures around 2400 K and the formation of flame zones. These results are consistent with adiabatic flame temperature predictions as well as particle morphology observations of a titanium core with a TiO2 surface. Late stage particle expansion, star fragmentation, and molten droplet breakup events are also observed using the time-resolved morphology and temperature diagnostics. These results illustrate the different stages of titanium particle combustion in pyrotechnic environments, which can be used to inform improvements in next-generation igniters.
Hobbs, Michael L.; Britt, Phillip F.; Hobbs, David T.; Kaneshige, Michael; Minette, Michael; Mintz, Jessica; Pennebaker, Frank M.; Parker, Gary R.; Rosenberg, David; Schwantes, Jon; Williams, Audrey; Pierce, Robert
Hobbs, Michael L.; Britt, Phillip F.; Hobbs, David T.; Kaneshige, Michael; Minette, Michael; Mintz, Jessica; Pennebaker, Frank M.; Parker, Gary R.; Pierce, Robert; Rosenberg, David; Schwantes, Jon; Williams, Audrey
Precise wording is important in every field of study, including operational procedures. Confusion in the wording “organic” and “inorganic” may have contributed to substitution of an organic kitty litter for an inorganic adsorbent used to prepare nuclear waste for disposal at an underground salt repository. Adsorbents prevent liquids like nitric acid from causing corrosion within the waste drums. However, combination of organic material with nitric acid can cause heat- and gas-generating reactions resulting in thermal runaway, rapid pressurization, and drum rupture. In 2014, waste Drum 68660 containing nitric acid-soaked organic kitty litter exploded and released transuranic waste into the repository. The cause of the accident was never identified. Here we show that the root cause of Drum 68660 igniting was restriction of the drum vent resulting in accelerated nitric acid chemistry, thermal runaway, and radiation dispersal.
In previous work, commercially available downward facing exploding bridgewire detonators (EBWs) were exposed to elevated temperatures. These detonators were then initiated using a firing set which discharged a high amplitude short duration electrical pulse into a thin gold bridgewire. Responses of the detonators were measured using photonic doppler velocimetry (PDV) and high-speed photography. A time delay of 2 μs between EBW initiation and first movement of an output flyer separated operable detonators from inoperable detonators or duds. In the current work, we propose a simple method to determine detonator operability from the calculated state of the detonator at the time the firing set is initiated. The failure criterion is based on the gap distance between the exploding bridgewire (EBW) and the adjacent initiating explosive within the detonator which is low-density pentaerythritol tetranitrate (PETN) that melts between 413–415 K (140–142 °C). The gap forms as PETN melts and flows to the bottom of the input pellet. Melting of PETN is modeled thermodynamically as an energy sink using a normal distribution spread over a temperature range between the onset temperature of 413 K and the ending temperature of 415 K. The extent of the melt is determined from the average temperature of the PETN. The PETN liquid is assumed to occupy the interstitial gas volume in the lower part of the input pellet. The vacated volume from the relocated liquid forms the gap between the EBW and the PETN. The remaining sandwiched layer consists of solid PETN particles and gas filling interstitial volume. We predict that a threshold gap between 17–27 μm separates properly functioning detonators from duds.
Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (∼406 K). The temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.
This is the Sandia report from a joint NSRD project between Sandia National Labs and Savannah River National Labs. The project involved development of simulation tools and data intended to be useful for tritium operations safety assessment. Tritium is a synthetic isotope of hydrogen that has a limited lifetime, and it is found at many tritium facilities in the form of elemental gas (T2). The most serious risk of reasonable probability in an accident scenario is when the tritium is released and reacts with oxygen to form a water molecule, which is subsequently absorbed into the human body. This tritium oxide is more readily absorbed by the body and therefore represents a limiting factor for safety analysis. The abnormal condition of a fire may result in conversion of the safer T2 inventory to the more hazardous oxidized form. It is this risk that tends to govern the safety protocols. Tritium fire datasets do not exist, so prescriptive safety guidance is largely conservative and reliant on means other than testing to formulate guidelines. This can have a consequence in terms of expensive and/or unnecessary mitigation design, handling protocols, and operational activities. This issue can be addressed through added studies on the behavior of tritium under representative conditions. Due to the hazards associated with the tests, this is being approached mainly from a modeling and simulation standpoint and surrogate testing. This study largely establishes the capability to generate simulation predictions with sufficiently credible characteristics to be accepted for safety guidelines as a surrogate for actual data through a variety of testing and modeling activities.
Transforming polymorphs, melting, and boiling are physical processes that can accelerate decomposition rates during cookoff of PETN and make measurements difficult. For example, splashing liquids from large bubbles filled with decomposition products clog pressure tubing in sealed experiments. Boil over can also extinguish thermal excursions in vented experiments making ignition difficult. For better measurements, we have modified the Sandia Instrumented Thermal Ignition (SITI) experiment to obtain better sealed and vented cookoff data for PETN by reducing the sample size and including additional gas space to prevent clogged tubing and boil over. Ignition times were not affected by 1) increasing the gas space by a factor of 3 in sealed SITI experiments or by 2) venting the decomposition gasses. That is, thermal ignition of PETN is not pressure dependent and the rate-limiting step during PETN decomposition likely occurs in the condensed phase. A simple decomposition model was calibrated using these observations and includes rate acceleration caused by melting and boiling. The model is used to predict internal temperatures, pressurization, and thermal ignition in a wide variety of experiments. The model is also used with SITI data to estimate the previously unreported latent enthalpy (5 J/g) associated with the α (PETN-I) to β (PETN-II) polymorphic phase transformation of PETN.
Cookoff experiments of powdered and pressed TATB-based plastic bonded explosives (PBXs) have been modeled using a pressure-dependent universal cookoff model (UCM) in combination with a micromechanics pressurization (MMP) model described in a companion paper. The MMP model is based on the accumulation of decomposition gases at nucleation sites that load the surrounding TATB crystals and binder. This is the first cookoff model to use an analytical mechanics solution for compressibility and thermal expansion to describe internal pressurization caused by both temperature and decomposition occurring within closed-pore explosives. This approach produces more accurate predictions of ignition time and pressurization within high-density explosives than simple equation-of-state models. The current paper gives details of the reaction chemistry, model parameters, predicted uncertainty, and validation using experiments from multiple laboratories with errors less than 6 %. The UCM/MMP model framework gives more accurate thermal ignition predictions for high density explosives that are initially impermeable to decomposition gases.
We have completed a series of both vented and sealed cookoff experiments of black powder and smokeless powder in our Sandia Instrumented Thermal Ignition (SITI) apparatus at bulk densities of 1078 and 729 kg/m3, respectively. The confining aluminum cylinder was ramped from room temperature to a set point temperature and then held at the setpoint temperature until ignition. The setpoint temperatures varied between 495 to 523 K for the black powder and 401 to 412 K for the more sensitive smokeless powder. The vented experiments show a significant delay in thermal ignition, indicating that the ignition is dependent on pressure. Post experimental debris shows greater violence for our smokeless powder experiments than our black powder experiments. A simplified universal cookoff model (UCM) was calibrated using the black powder and smokeless powder SITI data and used to predict pressurization and thermal ignition. The current work presents the first calibration of the UCM with a double base propellant. This work also presents the first pressure-dependent cookoff model for black powder and smokeless powder.
We have completed a series of vented and sealed cookoff experiments of the ε-polymorph of CL-20 in our Sandia Instrumented Thermal Ignition (SITI) apparatus using both powder and pressed pellets at nominal densities of 313 ± 8 kg/m3 and 1030 ± 4 kg/m3, respectively. The boundary temperature of our aluminum confinement cylinder was ramped in 10 minutes from room temperature to a prescribed set-point temperature ranging between 448 nd 468 K and held at the set-point temperature until ignition. A universal cookoff model (UCM) has been calibrated using the ε-CL-20 SITI data to predict pressurization and thermal ignition of ε-CL-20. The ignition model was validated by using one-dimensional time-to-explosion (ODTX) ignition data from a different laboratory. We found that a thirtyfold increase in the reaction rates due to liquefaction at 520 K could explain the high temperature ODTX cookoff data. The model gives a plausible explanation of why melting is important in fast cookoff events involving CL-20. Our model also gives support to 520 K as the liquefaction point of CL-20, which has different values in the literature.
The solubility of RDX (hexahydro-1,3,5-tri-nitro-1,3,5-triazine) in TNT (2,4,6-trinitrotoluene) at elevated temperatures is required to accurately predict the response of Comp-B3 (60:40 RDX:TNT) during accidents involving fire. As the temperature increases, the TNT component melts, the RDX partially dissolves in the liquid TNT, and the remaining RDX melts (203 ∘C) as the Comp-B thermally ignites. In the current work, we used a differential scanning calorimeter (DSC) to estimate the solubility of RDX in TNT at the melting point of RDX. Most DSC measurements of Comp-B3 do not show an RDX melt endotherm. The absence of an endotherm associated with the RDX melt has been interpreted as RDX being completely dissolved in TNT before reaching the melting point. We have observed that the endotherm is not absent, but is masked by exothermic reactions occurring at these elevated temperatures. We have inhibited the exothermic reactions by venting our DSC samples and measuring the RDX melt endotherm in our Comp-B3 samples at about 203 ∘C. Using the measured heat flow associated with the RDX melt and the latent melting enthalpy of RDX, we have approximated the solubility of RDX in TNT to be roughly 50–100 g RDX per 100 g TNT. The broad range is based on corrections for exothermic reactions occurring as the RDX melts.
Thermal conductivity has been determined for a variety of energetic materials (EMs) using finite element analysis (FEA) and cookoff data from the Sandia Instrumented Thermal Ignition (SITI) experiment. Materials studied include melt-cast, pressed, and low-density explosives. The low-density explosives were either prills or powders with some experiments run at pour density (not pressed). We have compared several of our thermal conductivities with those in the literature as well as investigated contact resistance between the confining aluminum and explosive, multidimensional heat transfer effects, and uncertainty in the thermocouple bead positions. We have determined that contact resistance is minimal in the SITI experiment, the heat transfer along the midplane is one-dimensional, and that uncertainty in the thermocouple location is greatest near the heated boundary. Our values of thermal conductivity can be used with kinetic mechanisms to accurately predict thermal profiles and energy dissipation during the cookoff of explosives.
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.
We have completed a series of small-scale cook-off experiments of ammonium nitrate (AN) prills in our Sandia Instrumented Thermal Ignition test at nominal packing densities of about 0.8 g/cm3. We increased the boundary temperature of our aluminum confinement cylinder from room temperature to a prescribed set-point temperature in 10 min. Our set-point temperature ranged from 508 to 538 K. The external temperature of the confining cylinder was held at the set-point temperature until ignition. We used type K thermocouples to measure temperatures associated with several polymorphic phase changes as well as melting and boiling. As the AN boiled, our thermocouples were destroyed by corrosion, which may have been caused by reaction of hot nitric acid (HNO3) with nickel to form nickel nitrate, Ni(NO3)2. Videos of the corroding thermocouples showed a green solution that was similar to the color of Ni(NO3)2. We found that ignition was imminent as the AN boiling point was exceeded. Ignition of the AN prills was modeled by solving the energy equation with an energy source due to desorption of moisture and decomposition of AN to form equilibrium products. A Boussinesq approximation was used in conjunction with the momentum equation to model flow of the liquid AN. We found that the prediction of ignition was not sensitive to small perturbations in the latent enthalpies.
We have used several configurations of the Sandia Instrumented Thermal Ignition (SITI) experiment to develop a pressure-dependent, four-step ignition model for a plastic bonded explosive (PBX 9407) consisting of 94 wt.% RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine), and a 6 wt.% VCTFE binder (vinyl chloride/chlorotrifluoroethylene copolymer). The four steps include desorption of water, decomposition of RDX to form equilibrium products, pressure-dependent decomposition of RDX forming equilibrium products, and decomposition of the binder to form hydrogen chloride and a nonvolatile residue (NVR). We address drying, binder decomposition, and decomposition of the RDX component from the pristine state through the melt and into ignition. We used Latin Hypercube Sampling (LHS) of the parameters to determine the sensitivity of the model to variation in the parameters. We also successfully validated the model using one-dimensional time-to-explosion (ODTX and P-ODTX) data from a different laboratory. Our SITI test matrix included 1) different densities ranging from 0.7 to 1.63 g/cm3, 2) free gas volumes ranging from 1.2 to 38 cm3, and 3) boundary temperatures ranging from 170 to 190 °C. We measured internal temperatures using embedded thermocouples at various radial locations as well as pressure using tubing that was connected from the free gas volume (ullage) to a pressure gauge. We also measured gas flow from our vented experiments. A borescope was included to obtain in situ video during some SITI experiments. We observed significant changes in the explosive volume prior to ignition. Our model, in conjunction with data observations, imply that internal accumulation of decomposition gases in high density PBX 9407 (90% of the theoretical maximum density) can contribute to significant strain whether or not the experiment is vented or sealed.
We combined the miniature cylinder expansion test (Mini-Cylex), with the Disk Acceleration Test (DAX) using Type K copper, Picatinny Liquid Explosive, and photonic Doppler velocimetry. We estimated the CJ state using plate reverberation methods at the test cap. We extracted velocities at 2, 7, and 10 volume expansions to fit Jones-Wilkins-Lee Equation of State and estimated Gurney velocity at the tube wall. The test also provides an additional method to estimate reaction products Hugoniot through knowledge of the copper test cap. Our experiments and simulations are within expected uncertainty. The test and the analysis effectively reduce costs while keeping or increasing fidelity.
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.
A diffusion-limited reaction model was calibrated for Al/Pt multilayers ignited on oxidized silicon, sapphire, and tungsten substrates, as well as for some Al/Pt multilayers ignited as free-standing foils. The model was implemented in a finite element analysis code and used to match experimental burn front velocity data collected from several years of testing at Sandia National Laboratories. Moreover, both the simulations and experiments reveal well-defined quench limits in the total Al + Pt layer (i.e., bilayer) thickness. At these limits, the heat generated from atomic diffusion is insufficient to support a self-propagating wave front on top of the substrates. Quench limits for reactive multilayers are seldom reported and are found to depend on the thermal properties of the individual layers. Here, the diffusion-limited reaction model is generalized to allow for temperature- and composition-dependent material properties, phase change, and anisotropic thermal conductivity. Utilizing this increase in model fidelity, excellent overall agreement is shown between the simulations and experimental results with a single calibrated parameter set. However, the burn front velocities of Al/Pt multilayers ignited on tungsten substrates are over-predicted. Possible sources of error are discussed and a higher activation energy (from 41.9 kJ/mol.at. to 47.5 kJ/mol.at.) is shown to bring the simulations into agreement with the velocity data observed on tungsten substrates. This higher activation energy suggests an inhibited diffusion mechanism present at lower heating rates.
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].
Al/Pt nanolaminates with a bilayer thickness (tb, width of an Al/Pt pair-layer) of 164 nm were irradiated with single laser pulses with durations of 10 ms and 0.5 ms at 189 W/cm2 and 1189 W/cm2, respectively. The time to ignition was measured for each pulse, and shorter ignition times were observed for the higher power/shorter pulse width. Videographic images of the irradiated area shortly after ignition show a non-uniform radial brightness for the longer pulse, while the shorter pulse shows uniform brightness. A diffusion-limited single step reaction mechanism was implemented in a finite element package to model the progress from reactants to products at both pulse widths. The model captures well both the observed ignition delay and qualitative observations regarding the non-uniform radial temperature.
We have previously developed a PBX 9501 cookoff model for the plastic bonded explosive PBX 9501 consisting of 95 wt% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazoncine (HMX), 2.5 wt% Estane® 5703 (a polyurethane thermoplastic), and 2.5 wt% of a nitroplasticizer (NP): BDNPA/F, a 50/50 wt% eutectic mixture bis(2,2-dinitropropyl)-acetal (BDNPA) and bis(2,2-dinitropropyl)-formal (BDNPF). This fivestep model includes desorption of water, decomposition of the NP to form NO2, reaction of the NO2 with Estane and HMX, and decomposition of HMX [1]. This model has been successfully validated with data from six laboratories with scales ranging from 2 g to more than 2.5 kg of explosive. We have determined, that the PBX 9501 model can be used to predict cookoff of other plastic bonded explosives containing HMX and an inert binder, such as LX-04 consisting of 85 wt% HMX and 15 wt% Viton A (vinylidine fluoride/hexafluoropropylene copolymer), LX-07 (90 wt% HMX and 10 wt% Viton A), LX- 10-0 (95 wt% HMX and 5 wt% Viton A), and LX-14 consisting of 95.5 wt % HMX and 4.5 wt% Estane® 5702-F1 (a polyurethane thermoplastic). Normally our cookoff models are verified using Sandia’s Instrumented Thermal Initiation (SITI) experiment. However, SITI data for LX-04, LX-07, LX-10-0, and LX-14 are not available at pressed density; although, some molding powder SITI data on LX-10-0 and LX-14 exists. Tarver and Tran [2] provide some one-dimensional time-to-explosion (ODTX) data for these explosives. The applicability of the PBX 9501 model to LX-04, LX-07, LX-10-0, AND LX-14 was made using this ODTX data [2]. The PBX 9501 model is applied to these other explosives by accounting for the correct amount of HMX in the explosive and limiting the NP reaction. We have found the PBX 9501 model to be useful for predicting the response of these PBXs to abnormal thermal environments such as fire.
We previously developed a PETN thermal decomposition model that accurately predicts thermal ignition and detonator failure [1]. This model was originally developed for CALORE [2] and required several complex user subroutines. Recently, a simplified version of the PETN decomposition model was implemented into ARIA [3] using a general chemistry framework without need for user subroutines. Detonator failure was also predicted with this new model using ENCORE. The model was simplified by 1) basing the model on moles rather than mass, 2) simplifying the thermal conductivity model, and 3) implementing ARIA’s new phase change model. This memo briefly describes the model, implementation, and validation.
We have used a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment to develop a pressure-dependent, five-step ignition model for a plastic bonded explosive (PBX 9501) consisting of 95 wt% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazoncine (HMX), 2.5 wt% Estane® 5703 (a polyurethane thermoplastic), and 2.5 wt% of a nitroplasticizer (NP): BDNPA/F, a 50/50 wt% eutectic mixture bis(2,2-dinitropropyl)-acetal (BDNPA) and bis(2,2-dinitropropyl)-formal (BDNPF). The five steps include desorption of water, decomposition of the NP to form NO2, reaction of the NO2 with Estane® and HMX, and decomposition of HMX. The model was fit using our experiments and successfully validated with experiments from five other laboratories with scales ranging from about 2 g to more than 2.5 kg of PBX. Our experimental variables included density, confinement, free gas volume, and temperature. We measured internal temperatures, confinement pressure, and ignition time. In some of our experiments, we used a borescope to visually observe the decomposing PBX. Our observations included the endothermic β–δ phase change of the HMX, a small exothermic temperature excursion in low-density unconfined experiments, and runaway ignition. We hypothesize that the temperature excursion in these low density experiments was associated with the NP decomposing exothermically within the PBX sample. This reactant-limited temperature excursion was not observed with our thermocouples in the high-density experiments. For these experiments, we believe the binder diffused to the edges of our high density samples and decomposed next to the highly conductive wall as confirmed by our borescope images.
Occasionally, our well-controlled cookoff experiments with Comp-B give anomalous results when venting conditions are changed. For example, a vented experiment may take longer to ignite than a sealed experiment. In the current work, we show the effect of venting on thermal ignition of Comp-B. We use Sandia's Instrumented Thermal Ignition (SITI) experiment with various headspace volumes in both vented and sealed geometries to study ignition of Comp-B. In some of these experiments, we have used a boroscope to observe Comp-B as it melts and reacts. We propose that the mechanism for ignition involves TNT melting, dissolution of RDX, and complex bubbly liquid flow. High pressure inhibits bubble formation and flow is significantly reduced. At low pressure, a vigorous dispersed bubble flow was observed.
Occasionally, our well-controlled cookoff experiments with Comp-B give anomalous results when venting conditions are changed. For example, a vented experiment may take longer to ignite than a sealed experiment. In the current work, we show the effect of venting on thermal ignition of Comp-B. We use Sandia's Instrumented Thermal Ignition (SITI) experiment with various headspace volumes in both vented and sealed geometries to study ignition of Comp-B. In some of these experiments, we have used a boroscope to observe Comp-B as it melts and reacts. We propose that the mechanism for ignition involves TNT melting, dissolution of RDX, and complex bubbly liquid flow. High pressure inhibits bubble formation and flow is significantly reduced. At low pressure, a vigorous dispersed bubble flow was observed.
A waste drum located 2150 feet underground may have been the root cause of a radiation leak on February 14, 2014. Information provided to the WIPP Technical Assessment Team (TAT) was used to describe the approximate content of the drum, which included an organic cat litter (Swheat Scoop®, or Swheat) composed of 100% wheat products. The drum also contained various nitrate salts, oxalic acid, and a nitric acid solution that was neutralized with triethanolamine (TEA). CTH-TIGER was used with the approximate drum contents to specify the products for an exothermic reaction for the drum. If an inorganic adsorbent such as zeolite had been used in lieu of the kitty litter, the overall reaction would have been endothermic. Dilution with a zeolite adsorbent might be a useful method to remediate drums containing organic kitty litter. SIERRA THERMAL was used to calculate the pressurization and ignition of the drum. A baseline simulation of drum 68660 was performed by assuming a background heat source of 0.5-10 W of unknown origin. The 0.5 W source could be representative of heat generated by radioactive decay. The drum ignited after about 70 days. Gas generation at ignition was predicted to be 300-500 psig with a sealed drum (no vent). At ignition, the wall temperature increases modestly by about 1°C, demonstrating that heating would not be apparent prior to ignition. The ignition location was predicted to be about 0.43 meters above the bottom center portion of the drum. At ignition only 3-5 kg (out of 71.6 kg total) has been converted into gas, indicating that most of the material remained available for post-ignition reaction.
Blast waves from an explosion in air can cause significant structural damage. As an example, cylindrically-shaped charges have been used for over a century as dynamite sticks for mining, excavation, and demolition. Near the charge, the effects of geometry, standoff from the ground, the proximity to other objects, confinement (tamping), and location of the detonator can significantly affect blast wave characteristics. Furthermore, nonuniformity in the surface characteristics and the density of the charge can affect fireball and shockwave structure. Currently, the best method for predicting the shock structure near a charge and the dynamic loading on nearby structures is to use a multidimensional, multimaterial shock physics code. However, no single numerical technique currently exists for predicting secondary combustion, especially when particulates from the charge are propelled through the fireball and ahead of the leading shock lens. Furthermore, the air within the thin shocked layer can dissociate and ionize. Hence, an appropriate equation of state for air is needed in these extreme environments. As a step towards predicting this complex phenomenon, a technique was developed to provide the equilibrium species composition at every computational cell in an air blast simulation as an initial condition for hand-off to other analysis codes for combustion fluid dynamics or radiation transport. Here, a bare cylindrical charge of TNT detonated in air is simulated using CTH, an Eulerian, finite volume, shock propagation code developed and maintained at Sandia National Laboratories. The shock front propagation is computed at early times, including the detonation wave structure in the explosive and the subsequent air shock up to 100 microseconds, where ambient air entrainment is not significant. At each computational cell, which could have TNT detonation products, air, or both TNT and air, the equilibrium species concentration at the density-energy state is computed using the JCZS2i database in the thermochemical code TIGER. This extensive database of 1267 gas (including 189 ionized species) and 490 condensed species can predict thermodynamic states up to 20,000 K. The results of these calculations provide the detailed three-dimensional structure of a thin shock front, and spatial species concentrations including free radicals and ions. Furthermore, air shock predictions are compared with experimental pressure gage data from a right circular cylinder of pressed TNT, detonated at one end. These complimentary predictions show excellent agreement with the data for the primary wave structure.
This report presents plots of specific heat, enthalpy, entropy, and Gibbs free energy for 1439 species in the JCZS2i database. Included in this set of species are 496 condensed-phase species and 943 gas-phase species. The gas phase species contain 80 anions and 112 cations for a total of 192 ions. The JCZS2i database is used in conjunction with the TIGER thermochemical code to predict thermodynamic states from ambient conditions to high temperatures and pressures. Predictions from the TIGER code using the JCZS2i database can be used in shock physics codes where temperatures may be as high as 20,000 K and ions may be present. Such high temperatures were not considered in the original JCZS database, and extrapolations made for these temperatures were unrealistic. For example, specific heat would sometimes go negative at high temperatures which fails the definition of specific heat. The JCZS2i database is a new version of the JCZS database that is being created to address these inaccuracies. The purpose of the current report is to visualize the high temperature extrapolations to insure that the specific heat, enthalpy, entropy, and Gibbs free energy predictions are reasonable up to 20,000 K.
A decomposition chemistry and heat transfer model to predict the response of removable epoxy foam (REF) exposed to fire-like heat fluxes is described. The epoxy foam was created using a perfluorohexane blowing agent with a surfactant. The model includes desorption of the blowing agent and surfactant, thermal degradation of the epoxy polymer, polymer fragment transport, and vapor-liquid equilibrium. An effective thermal conductivity model describes changes in thermal conductivity with reaction extent. Pressurization is modeled assuming: (1) no strain in the condensed-phase, (2) no resistance to gas-phase transport, (3) spatially uniform stress fields, and (4) no mass loss from the system due to venting. The model has been used to predict mass loss, pressure rise, and decomposition front locations for various small-scale and large-scale experiments performed by others. The framework of the model is suitable for polymeric foams with absorbed gases. Published by Elsevier Ltd.
A decomposition chemistry and heat transfer model to predict the response of removable epoxy foam (REF) exposed to fire-like heat fluxes is described. The epoxy foam was created using a perfluorohexane blowing agent with a surfactant. The model includes desorption of the blowing agent and surfactant, thermal degradation of the epoxy polymer, polymer fragment transport, and vapor-liquid equilibrium. An effective thermal conductivity model describes changes in thermal conductivity with reaction extent. Pressurization is modeled assuming: (1) no strain in the condensed-phase, (2) no resistance to gas-phase transport, (3) spatially uniform stress fields, and (4) no mass loss from the system due to venting. The model has been used to predict mass loss, pressure rise, and decomposition front locations for various small-scale and large-scale experiments performed by others. The framework of the model is suitable for polymeric foams with absorbed gases.
A case study is reported to document the details of a validation process to assess the accuracy of a mathematical model to represent experiments involving thermal decomposition of polyurethane foam. The focus of the report is to work through a validation process. The process addresses the following activities. The intended application of mathematical model is discussed to better understand the pertinent parameter space. The parameter space of the validation experiments is mapped to the application parameter space. The mathematical models, computer code to solve the models and its (code) verification are presented. Experimental data from two activities are used to validate mathematical models. The first experiment assesses the chemistry model alone and the second experiment assesses the model of coupled chemistry, conduction, and enclosure radiation. The model results of both experimental activities are summarized and uncertainty of the model to represent each experimental activity is estimated. The comparison between the experiment data and model results is quantified with various metrics. After addressing these activities, an assessment of the process for the case study is given. Weaknesses in the process are discussed and lessons learned are summarized.
A decomposition model has been developed to predict the response of removable syntactic foam (RSF) exposed to fire-like heat fluxes. RSF consists of glass micro-balloons (GMB) in a cured epoxy polymer matrix. A chemistry model is presented based on the chemical structure of the epoxy polymer, mass transport of polymer fragments to the bulk gas, and vapor-liquid equilibrium. Thermophysical properties were estimated from measurements. A bubble nucleation, growth, and coalescence model was used to describe changes in properties with the extent of reaction. Decomposition of a strand of syntactic foam exposed to high temperatures was simulated.
Response of removable epoxy foam (REF) to high heat fluxes is described using a decomposition chemistry model [1] in conjunction with a finite element heat conduction code [2] that supports chemical kinetics and dynamic radiation enclosures. The chemistry model [1] describes the temporal transformation of virgin foam into carbonaceous residue by considering breakdown of the foam polymer structure, desorption of gases not associated with the foam polymer, mass transport of decomposition products from the reaction site to the bulk gas, and phase equilibrium. The finite element foam response model considers the spatial behavior of the foam by using measured and predicted thermophysical properties in combination with the decomposition chemistry model. Foam elements are removed from the computational domain when the condensed mass fractions of the foam elements are close to zero. Element removal, referred to as element death, creates a space within the metal confinement causing radiation to be the dominant mode of heat transfer between the surface of the remaining foam elements and the interior walls of the confining metal skin. Predictions were compared to front locations extrapolated from radiographs of foam cylinders enclosed in metal containers that were heated with quartz lamps [3,4]. The effects of the maximum temperature of the metal container, density of the foam, the foam orientation, venting of the decomposition products, pressurization of the metal container, and the presence or absence of embedded components are discussed.
A Simple Removable Epoxy Foam (SREF) decomposition chemistry model has been developed to predict the decomposition behavior of an epoxy foam encapsulant exposed to high temperatures. The foam is composed of an epoxy polymer, blowing agent, and surfactant. The model is based on a simple four-step mass loss model using distributed Arrhenius reaction rates. A single reaction was used to describe desorption of the blowing agent and surfactant (BAS). Three of the reactions were used to describe degradation of the polymer. The coordination number of the polymeric lattice was determined from the chemical structure of the polymer; and a lattice statistics model was used to describe the evolution of polymer fragments. The model lattice was composed of sites connected by octamethylcylotetrasiloxane (OS) bridges, mixed product (MP) bridges, and bisphenol-A (BPA) bridges. The mixed products were treated as a single species, but are likely composed of phenols, cresols, and furan-type products. Eleven species are considered in the SREF model - (1) BAS, (2) OS, (3) MP, (4) BPA, (5) 2-mers, (6) 3-mers, (7) 4-mers, (8) nonvolatile carbon residue, (9) nonvolatile OS residue, (10) L-mers, and (11) XL-mers. The first seven of these species (VLE species) can either be in the condensed-phase or gas-phase as determined by a vapor-liquid equilibrium model based on the Rachford-Rice equation. The last four species always remain in the condensed-phase. The 2-mers, 3-mers, and 4-mers are polymer fragments that contain two, three, or four sites, respectively. The residue can contain C, H, N, O, and/or Si. The L-mer fraction consists of polymer fragments that contain at least five sites (5-mer) up to a user defined maximum mer size. The XL-mer fraction consists of polymer fragments greater than the user specified maximum mer size and can contain the infinite lattice if the bridge population is less than the critical bridge population. Model predictions are compared to 133-thermogravimetric analysis (TGA) experiments performed at 24 different conditions. The average RMS error between the model and the 133 experiments was 4.25%. The model was also used to predict the response of two other removable epoxy foams with different compositions as well as the pressure rise in a constant volume hot cell.
An efficient polymer mass loss and foam response model has been developed to predict the behavior of unconfined polyurethane foam exposed to fire-like heat fluxes. The mass loss model is based on a simple two-step mechanism using distributed reaction rates. The mass loss model was implemented into a multidimensional finite element heat conduction code that supports chemical kinetics and dynamic enclosure radiation. A discretization bias correction model was parameterized using elements with characteristic lengths ranging from 0.1 cm to 1 cm. Bias corrected solutions with these large elements gave essentially the same results as grid-independent solutions using 0.01-cm elements. Predictions were compared to measured decomposition front locations determined from real-time X-rays of 9-cm diameter, 15-cm tall cylinders of foam that were heated with lamps. The calculated and measured locations of the decomposition fronts were well within 1 cm of each other and in some cases the fronts coincided.
A Simple PolyUrethane Foam (SPUF) mass loss and response model has been developed to predict the behavior of unconfined, rigid, closed-cell, polyurethane foam-filled systems exposed to fire-like heat fluxes. The model, developed for the B61 and W80-0/1 fireset foam, is based on a simple two-step mass loss mechanism using distributed reaction rates. The initial reaction step assumes that the foam degrades into a primary gas and a reactive solid. The reactive solid subsequently degrades into a secondary gas. The SPUF decomposition model was implemented into the finite element (FE) heat conduction codes COYOTE [1] and CALORE [2], which support chemical kinetics and dynamic enclosure radiation using 'element death.' A discretization bias correction model was parameterized using elements with characteristic lengths ranging from 1-mm to 1-cm. Bias corrected solutions using the SPUF response model with large elements gave essentially the same results as grid independent solutions using 100-{micro}m elements. The SPUF discretization bias correction model can be used with 2D regular quadrilateral elements, 2D paved quadrilateral elements, 2D triangular elements, 3D regular hexahedral elements, 3D paved hexahedral elements, and 3D tetrahedron elements. Various effects to efficiently recalculate view factors were studied -- the element aspect ratio, the element death criterion, and a 'zombie' criterion. Most of the solutions using irregular, large elements were in agreement with the 100-{micro}m grid-independent solutions. The discretization bias correction model did not perform as well when the element aspect ratio exceeded 5:1 and the heated surface was on the shorter side of the element. For validation, SPUF predictions using various sizes and types of elements were compared to component-scale experiments of foam cylinders that were heated with lamps. The SPUF predictions of the decomposition front locations were compared to the front locations determined from real-time X-rays. SPUF predictions of the 19 radiant heat experiments were also compared to a more complex chemistry model (CPUF) predictions made with 1-mm elements. The SPUF predictions of the front locations were closer to the measured front locations than the CPUF predictions, reflecting the more accurate SPUF prediction of mass loss. Furthermore, the computational time for the SPUF predictions was an order of magnitude less than for the CPUF predictions.
A Chemical-structure-based PolyUrethane Foam (CPUF) decomposition model has been developed to predict the fire-induced response of rigid, closed-cell polyurethane foam-filled systems. The model, developed for the B-61 and W-80 fireset foam, is based on a cascade of bondbreaking reactions that produce CO2. Percolation theory is used to dynamically quantify polymer fragment populations of the thermally degrading foam. The partition between condensed-phase polymer fragments and gas-phase polymer fragments (i.e. vapor-liquid split) was determined using a vapor-liquid equilibrium model. The CPUF decomposition model was implemented into the finite element (FE) heat conduction codes COYOTE and CALORE, which support chemical kinetics and enclosure radiation. Elements were removed from the computational domain when the calculated solid mass fractions within the individual finite element decrease below a set criterion. Element removal, referred to as ?element death,? creates a radiation enclosure (assumed to be non-participating) as well as a decomposition front, which separates the condensed-phase encapsulant from the gas-filled enclosure. All of the chemistry parameters as well as thermophysical properties for the CPUF model were obtained from small-scale laboratory experiments. The CPUF model was evaluated by comparing predictions to measurements. The validation experiments included several thermogravimetric experiments at pressures ranging from ambient pressure to 30 bars. Larger, component-scale experiments were also used to validate the foam response model. The effects of heat flux, bulk density, orientation, embedded components, confinement and pressure were measured and compared to model predictions. Uncertainties in the model results were evaluated using a mean value approach. The measured mass loss in the TGA experiments and the measured location of the decomposition front were within the 95% prediction limit determined using the CPUF model for all of the experiments where the decomposition gases were vented sufficiently. The CPUF model results were not as good for the partially confined radiant heat experiments where the vent area was regulated to maintain pressure. Liquefaction and flow effects, which are not considered in the CPUF model, become important when the decomposition gases are confined.
Sensitivity/uncertainty analyses are not commonly performed on complex, finite-element engineering models because the analyses are time consuming, CPU intensive, nontrivial exercises that can lead to deceptive results. To illustrate these ideas, an analytical sensitivity/uncertainty analysis is used to determine the standard deviation and the primary factors affecting the burn velocity of polyurethane foam exposed to firelike radiative boundary conditions. The complex, finite element model has 25 input parameters that include chemistry, polymer structure, and thermophysical properties. The response variable was selected as the steady-state burn velocity calculated as the derivative of the burn front location versus time. The standard deviation of the burn velocity was determined by taking numerical derivatives of the response variable with respect to each of the 25 input parameters. Since the response variable is also a derivative, the standard deviation is essentially determined from a second derivative that is extremely sensitive to numerical noise. To minimize the numerical noise, 50-micron elements and approximately 1-msec time steps were required to obtain stable uncertainty results. The primary effect variable was shown to be the emissivity of the foam.
The decomposition of unconfined rigid polyurethane foam has been modeled by a kinetic bond-breaking scheme describing degradation of a primary polymer and formation of a thermally stable secondary polymer. The bond-breaking scheme is resolved using percolation theory to describe evolving polymer fragments. The polymer fragments vaporize according to individual vapor pressures. Kinetic parameters for the model were obtained from thermal gravimetric analysis (TGA). The chemical structure of the foam was determined from the preparation techniques and ingredients used to synthesize the foam. Scale-up effects were investigated by simulating the response of an incident heat flux of 25 W/cm2 on a partially confined 8.8-cm diameter by 15-cm long right circular cylinder of foam that contained an encapsulated component. Predictions of internal foam and component temperatures, as well as regression of the foam surface, were in agreement with measurements using thermocouples and X-ray imaging.
The decomposition of unconfined rigid polyurethane foam has been modeled by a kinetic bond-breaking scheme describing degradation of a primary polymer and formation of a thermally stable secondary polymer. The bond-breaking scheme is resolved using percolation theory to describe evolving polymer fragments. The polymer fragments vaporize according to individual vapor pressures. Kinetic parameters for the model were obtained from Thermal Gravimetric Analysis (TGA). The chemical structure of the foam was determined from the preparation techniques and ingredients used to synthesize the foam. Scale-up effects were investigated by simulating the response of an incident heat flux of 25 W/cm{sup 2} on a partially confined 8.8-cm diameter by 15-cm long right circular cylinder of foam which contained an encapsulated component. Predictions of center, midradial, and component temperatures, as well as regression of the foam surface, were in agreement with measurements using thermocouples and X-ray imaging.
Exponential-13,6 (EXP-13,6) potential pammeters for 750 gases composed of 48 elements were determined and assembled in a database, referred to as the JCZS database, for use with the Jacobs Cowperthwaite Zwisler equation of state (JCZ3-EOS)~l) The EXP- 13,6 force constants were obtained by using literature values of Lennard-Jones (LJ) potential functions, by using corresponding states (CS) theory, by matching pure liquid shock Hugoniot data, and by using molecular volume to determine the approach radii with the well depth estimated from high-pressure isen- tropes. The JCZS database was used to accurately predict detonation velocity, pressure, and temperature for 50 dif- 3 Accurate predictions were also ferent explosives with initial densities ranging from 0.25 glcm3 to 1.97 g/cm . obtained for pure liquid shock Hugoniots, static properties of nitrogen, and gas detonations at high initial pressures.
A database has been created for use with the Jacobs-Cowperthwaite-Zwisler-3 equation-of-state (JCZ3-EOS) to determine thermochemical equilibrium for detonation and expansion states of energetic materials. The JCZ3-EOS uses the exponential 6 intermolecular potential function to describe interactions between molecules. All product species are characterized by r*, the radius of the minimum pair potential energy, and {var_epsilon}/k, the well depth energy normalized by Boltzmann`s constant. These parameters constitute the JCZS (S for Sandia) EOS database describing 750 gases (including all the gases in the JANNAF tables), and have been obtained by using Lennard-Jones potential parameters, a corresponding states theory, pure liquid shock Hugoniot data, and fit values using an empirical EOS. This database can be used with the CHEETAH 1.40 or CHEETAH 2.0 interface to the TIGER computer program that predicts the equilibrium state of gas- and condensed-phase product species. The large JCZS-EOS database permits intermolecular potential based equilibrium calculations of energetic materials with complex elemental composition.
Rigid polyurethane foams are used as encapsulants to isolate and support thermally sensitive components within weapon systems. When exposed to abnormal thermal environments, such as fire, the polyurethane foam decomposes to form products having a wide distribution of molecular weights and can dominate the overall thermal response of the system. Decomposing foams have either been ignored by assuming the foam is not present, or have been empirically modeled by changing physical properties, such as thermal conductivity or emissivity, based on a prescribed decomposition temperature. The hypothesis addressed in the current work is that improved predictions of polyurethane foam degradation can be realized by using a more fundamental decomposition model based on chemical structure and vapor-liquid equilibrium, rather than merely fitting the data by changing physical properties at a prescribed decomposition temperature. The polyurethane decomposition model is founded on bond breaking of the primary polymer and formation of a secondary polymer which subsequently decomposes at high temperature. The bond breaking scheme is resolved using percolation theory to describe evolving polymer fragments. The polymer fragments vaporize according to individual vapor pressures. Kinetic parameters for the model were obtained from Thermal Gravimetric Analysis (TGA) from a single nonisothermal experiment with a heating rate of 20 C/min. Model predictions compare reasonably well with a separate nonisothermal TGA weight loss experiment with a heating rate of 200 C/min.
Rigid polyurethane foams are frequently used as encapsulants to isolate and support thermally sensitive components within weapon systems. When exposed to abnormal thermal environments, such as fire, the polyurethane foam decomposes to form products having a wide distribution of molecular weights and can dominate the overall thermal response of the system. Mechanical response of the decomposing foam, such as thermal expansion under various loading conditions created by gas generation, remains a major unsolved problem. A constitutive model of the reactive foam is needed to describe the coupling between mechanical response and chemical decomposition of foam exposed to environments such as fire. Towards this end, a reactive elastic-plastic constitutive model based on bubble mechanics describing nucleation, decomposition chemistry, and elastic/plastic mechanical behavior of rigid polyurethane foam has been developed. A local force balance, with mass continuity constraints, forms the basis of the constitutive model requiring input of temperature and the fraction of the material converted to gas. This constitutive model provides a stress-strain relationship which is applicable for a broad class of reacting materials such as explosives, propellants, pyrotechnics, and decomposing foams. The model is applied to a block of foam exposed to various thermal fluxes. The model is also applied to a sphere of foam confined in brass. The predicted mechanical deformation of the foam block and sphere are shown to qualitatively agree with experimental observations.
Determination of product species, equations-of-state (EOS) and thermochemical properties of high explosives and pyrotechnics remains a major unsolved problem. Although, empirical EOS models may be calibrated to replicate detonation conditions within experimental variability (5--10%), different states, e.g. expansion, may produce significant discrepancy with data if the basic form of the EOS model is incorrect. A more physically realistic EOS model based on intermolecular potentials, such as the Jacobs Cowperthwaite Zwisler (JCZ3) EOS, is needed to predict detonation states as well as expanded states. Predictive capability for any EOS requires a large species data base composed of a wide variety of elements. Unfortunately, only 20 species have known JCZ3 molecular force constants. Of these 20 species, only 10 have been adequately compared to experimental data such as molecular scattering or shock Hugoniot data. Since data in the strongly repulsive region of the molecular potential is limited, alternative methods must be found to deduce force constants for a larger number of species. The objective of the present study is to determine JCZ3 product species force constants by using a corresponding states theory. Intermolecular potential parameters were obtained for a variety of gas species using a simple corresponding states technique with critical volume and critical temperature. A more complex, four parameter corresponding state method with shape and polarity corrections was also used to obtain intermolecular potential parameters. Both corresponding state methods were used to predict shock Hugoniot data obtained from pure liquids. The simple corresponding state method is shown to give adequate agreement with shock Hugoniot data.
A summary of multidimensional modeling is presented which describes coupled thermals chemical and mechanical response of reactive and nonreactive materials. This modeling addresses cookoff of energetic material (EM) prior to the onset of ignition. Cookoff, lasting from seconds to days, sensitizes the EM whereupon combustion of confined, degraded material determines the level of violence. Such processes are dynamic, occurring over time scales of millisecond to microsecond, and thus more amenable for shock physics analysis. This work provides preignition state estimates such as the amount of decomposition, morphological changes, and quasistatic stress states for subsequent dynamic analysis. To demonstrate a fully-coupled thermal/chemical/quasistatic mechanical capability, several example simulations have been performed: (1) the one-dimensional time-to-explosion experiments, (2) the Naval Air Weapon Center`s (NAWC) small scale cookoff bomb, (3) a small hot cell experiment and (4) a rigid, highly porous, closed-cell polyurethane foam. Predictions compared adequately to available data. Deficiencies in the model and future directions are discussed.
Cookoff modeling of confined energetic materials involves the coupling of thermal, chemical and mechanical effects. In the past, modeling has focussed on the prediction of thermal runaway with little regard to the effects of mechanical behavior of the energetic material. To address the mechanical response of the energetic material, a constitutive submodel has been developed which can be incorporated into thermal-chemical-mechanical analysis. This work presents development of this submodel and its incorporation into a fully coupled one-dimensional, thermal-chemical-mechanical computer code to simulate thermal initiation of energetic materials. Model predictions include temperature, chemical species, stress, strain, solid/gas pressure, solid/gas density, yield function, and gas volume fraction. Sample results from a scaled aluminum tube filled with RDX exposed to a constant temperature bath at 500 K will be displayed. The micromechanical submodel is based on bubble mechanics which describes nucleation, decomposition, and elastic/plastic mechanical behavior. This constitutive material description requires input of temperatures and reacted fraction of the energetic material as provided by the reactive heat flow code, XCHEM, and the mechanical response is predicted using a quasistatic mechanics code, SANTOS. A parametric sensitivity analysis indicates that a small degree of decomposition causes significant pressurization of the energetic material, which implies that cookoff modeling must consider the strong interaction between thermal-chemistry and mechanics. This document consists of view graphs from the poster session.
An eXplosive CHEMical kinetics code, XCHEM was developed to solve the reactive diffusion equations associated with thermal ignition of energetic material. This method-of-lines code uses stiff numerical methods and adaptive meshing. Solution accuracy is maintained between multilayered materials consisting of blends of reactive components and/or inert materials. Phase change and variable properties are included in one-dimensional slab, cylindrical and spherical geometries. Temperature-dependent thermal properties was incorporated and modification of thermal conductivities to include decomposition effects are estimated using solid/gas volume fractions determined by species fractions. Gas transport properties are also included. Time varying temperature, heat flux, convective and thermal radiation boundary conditions, and layer to layer contact resistances are also implemented. The global kinetic mechanism developed at Lawrence Livermore National Laboratory (LLNL) by McGuire and Tarver used to fit One-Dimensional Time to eXplosion (ODTX) data for the conventional energetic materials (HMX, RDX, TNT, and TATB) are presented as sample calculations representative of multistep chemistry. Calculated and measured ignition times for explosive mixtures of Comp B (RDX/TNT), Octol, (HMX/TNT), PBX 9404 (HMX/NC), and RX-26-AF (HMX/TATB) are compared. Geometry and size effects are accurately modeled, and calculations are compared to experiments with time varying boundary conditions. Finally, XCHEM calculations of initiation of an AN/oil/water emulsion, resistively heated, are compared to measurements.
The Becker-Kistiakowsky-Wilson equation of state (BKW-EOS) has been calibrated over a wide initial density range near C-J states using measured detonation properties from 62 explosives at III total initial densities. Values for the empirical BKW constants {alpha}, {beta}, {kappa}, and {theta} were 0.5, 0.298, 10.5, and 6620, respectively. Covolumes were assumed to be invariant. Model evaluation includes comparison to measurements from 91 explosives composed of combinations of Al, B, Ba, C, Ca, Cl, F, H, N, 0, P, Pb, and Si at 147 total initial densities. Adequate agreement between predictions and measurements were obtained with a few exceptions for nonideal explosives. However, detonation properties for the nonideal explosives can be predicted adequately by assuming partial equilibrium. The partial equilibrium assumption was applied to aluminized composites of RDX, HMX, TNETB, and TNT to predict detonation velocity and temperature.
Variations of bed void fraction in a full-scale, reacting, fixed-bed coal gasifier have been deduced from measured axial pressure profiles obtained during gasification of seven coal types ranging from lignite to bituminous. Packed-bed pressure correlations were used to calculate the void fractions based on monotonic polynomial fits of measured pressure profiles. Insights into the fixed-bed combustion processes affected by the void distribution were obtained by a one-dimensional, steady-state, fixed-bed combustion model. Predicted temperature profiles from this model compare reasonably well to experimental data. The bed void distributions are not linear but are perturbed by vigorous reactions in the devolatilization and oxidation zones. Results indicate that a dramatic increase in temperature and associated gas release causes the bed to expand and the gas void space to increase. Increased void space localized in the combustion zone causes the steep temperature gradient to decrease and the location of the maximum temperature to shift. Also, large feed gas flow rates cause the void fraction in the ash zone to increase.