Correlations for Predicting Detonation Temperature of Pure and Mixed CNO and CHNO Explosives
Total Page:16
File Type:pdf, Size:1020Kb
Indian Journal of Engineering & Materials Sciences Vol. 12, April 2005, pp. 158-164 Correlations for predicting detonation temperature of pure and mixed CNO and CHNO explosives M H Keshavarz Department of Chemistry, Malek-ashtar University of Technology, Shahin-shahr P.O. Box 83145/115, Islamic Republic of Iran Received 25 May 2004; accepted 10 January 2005 In this paper, detonation temperature of CNO and CHNO explosives is evaluated by the base of atomic composition and heat of formation. It is shown here how only atomic composition and condensed phase heat of formation of CNO and CHNO explosives are sufficient for reliable prediction of detonation temperature as compared to complicated computer codes. New correlations are introduced so that they can easily be applied for determining detonation temperature via specified different pathways of decomposition of explosives. Calculated detonation temperatures by this procedure for both pure and explosive formulations show good agreement with respect to measured and the BKWS-EOS predictions of detonation temperatures. IPC Code: C06B 25/00 The detonation parameters such as heat of detonation, Kelvin and hundreds of kbar, as well as at much lower detonation temperature, detonation velocity and temperatures and pressures obtained during expansion detonation pressure are a variety of performance of the reaction products. Different computer codes parameters for measuring the effectiveness of such as BKW1 and RUBY2 and latter's offspring different explosives. When an explosive is initiated to TIGER3, CHEQ4, and CHEETAH5 assume all of the detonation, a detonation wave has traversed through chemical equations for all possible species in the its medium so that the hot gases and solid residues left reaction product gases and solving these with thermo- behind continue to exert pressure. Since detonation chemical analogues. They can also estimate the reaction of an explosive is extremely fast, the heat isentropic expansion having the equilibrium energy liberated by detonation will raise the temperature of and gas quantities along with the Rankine-Hugoniot gases, which will in turn cause them to expand and jump equations. For example CHEETAH, a C version work on surroundings. Most of the work of an of the FORTRAN equilibrium code TIGER, is an explosive in detonation is performed by the equilibrium thermochemical code used to estimate detonation reaction products. detonation parameters for explosives, as well as parameters other conditions such as a constant volume Typical detonation products of high explosives explosion. The detonation parameters can usually be with the elements carbon, hydrogen, nitrogen and calculated by a computer code when the heat of oxygen at high pressures (p~1-100 GPa) and formation and the density of the explosive substance temperatures (T~1000-5000 K) simultaneously are are known and the equation of state (EOS) is CO, N2, CO2, H2O, and solid carbon. Releasing a assumed. There are several forms of EOSs such as large amount of energy and complex chemical Becker-Kistiakosky-Wilson (BKW)6, the Jacobs- reactions are initiated to sustain the detonation Cowperthwaite-Zwisler (JCZ)7,8, Kihara-Hikita- process. The application of the hydrodynamic theory Tanaka (KHT)9, Exp-610 and Lennard-Jones- for calculating the detonation properties requires Denvoshire (LJD)11. Of different EOSs, the BKW- knowledge of the equation of state of the system that EOS is used extensively to calculate detonation can accurately reflect the thermodynamic properties properties1. The BKWC-EOS5, BKWR-EOS12 and of multicomponent mixtures at several thousand BKWS-EOS13 are three different parameterizations of the BKW-EOS so that BKWS-EOS is one of the best _________ E-mail: [email protected]; [email protected] EOSs for predicting detonation temperatures. KESHAVARZ: PREDICTION OF DETONATION TEMPERATURE OF CNO AND CHNO 159 Predicting fairly accurate data, by simple empirical interacting with a detector can be used for measuring methods, are highly desired for calculating the various detonation temperatures. Since measurements in detonation parameters of explosives. There is a porous systems include the effects of shocked air or continuing need in the field of explosives for reliable perhaps low density explosive material jetting into predictions of detonation parameters. Especially voids rather than the brightness of the pure detonation knowledge of the detonation temperature is important products, void free systems such as liquid explosives or because it is difficult to measure accurately. This single crystal systems may be believed to be more study had as its original purpose to develop a simple accurate. Estimated absolute accuracies are ±100 K for reliable method for calculating detonation liquid explosives and ±200 K for solid explosives32. temperatures of high explosives with the elements Detonation temperature is maximum temperature carbon, hydrogen, oxygen and nitrogen as compared that can be obtained by assuming that heat of to complex computer codes that is usually used to decomposition of explosive is used entirely to heat the compute detonation parameters. The purpose of this products. Since decomposition reaction is very fast, it work was to correlate detonation temperature with the can be assumed that energy released is completely explosive's elemental composition and condensed transferred to detonation products during a very short phase heat of formation of the explosive. This work time. The approximate detonation temperature can be shows useful applications of derived correlations to obtained under adiabatic conditions from molar heat both pure and mixed explosives. The calculated capacities of detonation products and heat of detonation temperature will compare with measured detonation. Some pathways for ideal and less ideal and computed results of the BKWS-EOS. It should explosives are recently suggested20 which can be be noted that predicted results are comparable with given for CNO and CHNO explosives as follows: outputs of complex computer codes and the accuracy is not necessarily enhanced by greater complexity. da≤ Furthermore, an important interesting observation of cb⎛⎞ this work is that much simpler equations could be NCOC(s)H22++−+dad() ⎜⎟ written for desk calculation of detonation temperature, 22⎝⎠ with about the same reliability as one might attach to … (1a) the more complex computer output. d>a and Correlations of Detonation Temperatures b/2 ≥d-a Though detonation parameters may be measured experimentally or calculated from theory, theoretical cb⎛⎞ NCOC(s)H22++−+dad() ⎜⎟ calculations are more convenient and reduce the costs 22⎝⎠ associated with synthesis, test and evaluation of the materials. Numerous studies have been developed for … (1b) predicting the Chapman-Jouguet (C-J) detonation d>a+b/2 pressures and velocities of pure or mixture of 14-30 CaHbNcOd and different classes of explosives . To compare with d ≤ 2a+b/2 equilibrium calculations, detonation velocities are measured at various charge diameters and cb⎛⎞ ⎛ b ⎞ Nad22++−+⎜⎟HO ⎜ 2 ⎟ CO extrapolated to an "infinite diameter". Detonation 22⎝⎠ ⎝ 2 ⎠ velocities can typically be measured to within a few percent, meanwhile the detonation pressures ⎛⎞b +−−⎜⎟da CO2 … (1c) determined by indirect method span a range of 10- ⎝⎠2 20% that non-equilibrium effects in reaction zones may contribute to large uncertainty31. If the d>2a+b/2 measurement is taken behind the von Neumann spike cb2 db− and in front of C-J plane, the measured pressures in ⎛⎞ ⎛ ⎞ NHOCO22+++−⎜⎟aa 2 ⎜ ⎟ O 2 this situation may be higher than equilibrium 22⎝⎠ ⎝ 4 ⎠ calculations. The brightness of the detonation front … (1d) 160 INDIAN J. ENG. MATER. SCI., APRIL 2005 To obtain reliable correlations for predicting In the definition of the parameter of ΔT, Td is detonation temperature, experimental detonation detonation temperature in K and ΔHf is condensed temperatures were used to optimize the correlations. heat of formation of explosive in kJ/mol. The glossary Optimized correlations corresponding to Eqs (1a)- of compounds is given in Appendix A. Detonation (1d) can also be given as: temperature estimated by this method for many different underoxidized and overoxidized CNO and Δ−Hdf 529.4 CHNO explosives are given in Table 1 and compared Δ=T 10.95×− 10−−−333ab 0.1132 +×−× 13.35 10 c 99.1 10 dwith corresponding BKWS-EOS as well as the … (2a) measured values. As indicated in Table 1, the new hand calculated detonation temperature shows surprisingly very good agreement with experimental Δ−Hadf 943.4 + 1229.5 Δ=T and one of the best complicated EOS, namely BKWS- −+×+×+0.1914abcd 59.67 10−−33 16.87 10 0.2224 EOS. Comparison of calculated results with … (2b) experimental data and BKWS-EOS listed in Table 1 may be taken as appropriate validation tests of the Δ−Habd1158.3 − 252.3 + 847.9 Δ=T f introduced simple correlations for use with CNO and −−×+×+0.2964abcd 55.09 10−−33 18.66 10 0.1911 CHNO explosives. It is worthwhile to note that the … (2c) present method is exceedingly simple and at the same time gives results that are comparable predicted to the Δ+Hab625.2 − 142.8 Δ=T f complexities of the other methods involving the 59.05×−×+×+× 10−−−−3333abcd 43.81 10 18.66 10 20.36 10 equation of state of the products. Considering few … (2d) percent deviations generally attributed to experimental measurements and BKW-EOS of where detonation temperatures, the agreement between calculated and measured temperatures is also Δ=TTd −298 … (3) satisfactory. Table 1⎯Comparison of detonation temperature of the new correlations and BKWS-EOS with measured values a 13 13 Name Condensed phase ρ0 (g/cc) Texp (K) Td (K) TBKWS-EOS %Dev BKWS- %Dev 13 ΔHf (kJ/mol) (K) EOS new at 298 K TNT -62.81 1.64 3401 3720 1.45 3401 3780 1.36 3401 3790 1.00 3400 3401 3750 -10.29 -0.03 0.80 3401 3610 ABH 485.69 1.64 4984 4710 TATB -154.29 1.88 3251 3250 1.85 3251 3260 DATB -98.81 1.80 3506 3550 1.78 3506 3690 DIPAM -28.47 1.76 4212 4040 HNAB 284.3 1.60 4621 4620 HNS 78.3 1.60 4136 4150 1.70 4136 4120 NONA 114.72 1.70 4930 4510 TACOT 464.76 1.85 4041 4040 Contd.