Publications

Results 1–25 of 28
Skip to search filters

Prediction of spatial distributions of equilibrium product species from high explosive blasts in air

50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference 2014

Brundage, Aaron B.; Attaway, Stephen W.; Hobbs, Michael L.; Kaneshige, Michael J.; Boye, Lydia A.

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.

More Details

Modeling compressive reaction and estimating model uncertainty in shock loaded porous samples of hexanitrostilbene (HNS)

AIP Conference Proceedings

Brundage, Aaron B.; Gump, J.C.

Neat pressings of HNS powders have been used in many explosive applications for over 50 years. However, characterization of its crystalline properties has lagged that of other explosives, and the solid stress has been inferred from impact experiments or estimated from mercury porosimetry. This lack of knowledge of the precise crystalline isotherm can contribute to large model uncertainty in the reacted response of pellets to shock impact. At high impact stresses, deflagration-to-detonation transition (DDT) processes initiated by compressive reaction have been interpreted from velocity interferometry at the surface of distended HNS-FP pellets. In particular, the Baer-Nunziato multiphase model in CTH, Sandia's Eulerian, finite volume shock propagation code, was used to predict compressive waves in pellets having approximately a 60% theoretical maximum density (TMD). These calculations were repeated with newly acquired isothermal compression measurements of fineparticle HNS using diamond anvil cells to compress the sample and powder x-ray diffraction to obtain the sample volume at each pressure point. Hence, estimating the model uncertainty provides a simple method for conveying the impact of future model improvements based upon new experimental data. © 2012 American Institute of Physics.

More Details

Unreacted equation of state development and multiphase modeling of dynamic compaction of low density hexanitrostilbene (HNS) pressings

Proceedings - 14th International Detonation Symposium, IDS 2010

Brundage, Aaron B.

Compaction waves in porous energetic materials have been shown to induce reaction under impact loading. In the past, simple two-state burn models such as the Arrhenius Burn model have been developed to predict slapper initiation in Hexanitrostilbene (HNS) pellets; however, a more sophisticated, fundamental approach is needed to predict the shock response during impact loading, especially in pellets that have been shown to have strong density gradients. The intergranular stress measures the resistance to bed compaction or the removal of void space due to particle packing and rearrangement. A constitutive model for the intergranular stress is needed for closure in the Baer-Nunziato (BN) multiphase mixture theory for reactive energetic materials. The intergranular stress was obtained from both quasi-static compaction experiments and from dynamic compaction experiments. Additionally, historical data and more recently acquired data for porous pellets compacted to high densities under shock loading were used for model assessment. Predicted particle velocity profiles under dynamic compaction were generally in good agreement with the experimental data. Hence, a multiphase model of HNS has been developed to extend current predictive capability.

More Details

Microenergetic shock initiation studies on deposited films of PETN

AIP Conference Proceedings

Tappan, Alexander S.; Wixom, Ryan R.; Trott, Wayne T.; Long, Gregory L.; Knepper, Robert; Brundage, Aaron B.; Jones, David A.

Films of the high explosive PETN (pentaerythritol tetranitrate) up to 500-μm thick have been deposited through physical vapor deposition, with the intent of creating well-defined samples for shock-initiation studies. PETN films were characterized with microscopy, x-ray diffraction, and focused ion beam nanotomography. These high-density films were subjected to strong shocks in both the out-of-plane and in-plane orientations. Initiation behavior was monitored with high-speed framing and streak camera photography. Direct initiation with a donor explosive (either RDX with binder, or CL-20 with binder) was possible in both orientations, but with the addition of a thin aluminum buffer plate (in-plane configuration only), initiation proved to be difficult. Initiation was possible with an explosively-driven 0.13-mm thick Kapton flyer and direct observation of initiation behavior was examined using streak camera photography at different flyer velocities. Models of this configuration were created using the shock physics code CTH. © 2009 American Institute of Physics.

More Details

Mesoscale simulations of shock initiation in energetic materials characterized by three-dimensional nanotomography

AIP Conference Proceedings

Brundage, Aaron B.; Wixom, Ryan R.; Tappan, Alexander S.; Long, G.T.

Three-dimensional shock simulations of energetic materials have been conducted to improve our understanding of initiation at the mesoscale. Vapor-deposited films of PETN and pressed powders of HNS were characterized with a novel three-dimensional nanotomographic technique. Detailed microstructures were constructed experimentally from a stack of serial electron micrographs obtained by successive milling and imaging in a dual-beam FIB/SEM. These microstructures were digitized and imported into a multidimensional, multimaterial Eulerian shock physics code. The simulations provided insight into the mechanisms of pore collapse in PETN and HNS samples with distinctly different three-dimensional pore morphology and distribution. This modeling effort supports investigations of microscale explosive phenomenology and elucidates mechanisms governing initiation of secondary explosives. © 2009 American Institute of Physics.

More Details
Results 1–25 of 28
Results 1–25 of 28