Publications

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

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.