Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives

Transcription

1 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 7 Central European Journal of Energetic Materials,, 9(), 7-38 ISSN Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives Waldemar A. TRZCIŃSKI Military University of Technology, Faculty of Advanced Technology and Chemistry Kaliskiego, -98 Warsaw, Poland wtrzcinski@wat.edu.pl Abstract: Theoretical models proposed in the literature for the deflagration-todetonation transition (DDT) in cast explosives are evaluated for primary explosives (complex compounds) in this work. The one-dimensional model of burning (deflagration), consistent with the classical Chapman-Jouguet theory and a model of burning under the conditions of zero mass velocity behind the flame front are presented, and the physical phenomena accompanying the accelerating wave of flame in solid explosives are described. The results of calculations taken from the literature are presented for the cast high explosive (pentolite). The model of acceleration of the deflagration wave was used to estimate the time and distance at which the process of burning leads to the emergence of a shock wave in primary explosives. The influence of burning rate and the physical properties of an explosive on the distance of deflagration to detonation transition is analysed. Keywords: explosives, deflagration to detonation transition Introduction From the huge number of experimental studies of the process of deflagration to detonation transition (DDT) in the detonation of condensed explosives, the following conclusions can be drawn []: the most important difference between the DDT process in primary and secondary explosives is the speed of transition from deflagration to detonation in the case of primaries the process is much faster; because the difference is rather quantitative than qualitative, all the

2 8 W.A. Trzciński relationships derived for secondaries can also be applied to primaries; it is far easier to build a quantitative model for the DDT process in cast explosives than in porous materials. The first conclusion comes from the prime function of primary and secondary explosives. Since primaries need to have undergone a violent detonation, if only by being exposed to the heat pulse, whilst the detonation of secondary explosives is mainly initiated by the shock wave. It is commonly believed that the faster the DDT, the better the primary material is. The ease of the transition to detonation is also a weakness of primary explosives, because accidental ignition leads to detonation. The second conclusion results from the observation of many similarities between primary and secondary explosives in the DDT process. Both types of materials are subject to this process, but in the case of secondary explosives the transition takes longer and the charge needs to be enclosed in heavy confinement. The final conclusion follows from the results of the DDT process in cast explosives and good agreement between the experimentally determined time and distance of detonation development with data received from the one-dimensional model of Macek []. There is no similar model for the transition to detonation in the case of porous explosives. The experimental data confirm the existence of a convective flame front with high surface area in these materials and, consequently, the formation of a complex flow field, affecting the DDT process. Over the years, a consensus has been reached on the mechanism controlling the DDT process [3-6]. Convective flow and hot-gas ignition contribute only initial, low velocity phenomena. The formation of a compaction wave in undisturbed material begins a series of non-convective events that lead to the final shockto-detonation (type I DDT). Another mechanism was observed for low density charges. In this case, the outbreak of detonation is preceded by a rapid build-up of pressure at the point where the detonation occurs (type II DDT) [7-9]. There are complex numerical models in which the porous explosive is treated as a two-phase mixture of interacting continua, consisting of grains of reacting solid and the gaseous products of combustion [-]. At the moment, there is no simple physical model taking into account the complexity of the DDT in porous explosives. In this paper, the one-dimensional model proposed in [] and developed in [] is used to predict the path of transition of burning to detonation in primary explosives. In the first part, Macek s model of the development of deflagration to detonation in cast secondary explosives is described. A critical analysis of this model given in work [] is also presented. This is followed by a description of the one-dimensional model of burning (deflagration), consistent with the classical

3 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 9 CJ theory, and the physical phenomena accompanying the accelerating wave of flame in solid explosives is discussed. Details of models of the DDT analyzed in [] and the results of the calculations presented are considered. Finally, the acceleration model is used to estimate the time and distance at which the process of burning leads to the emergence of a detonation wave in primary explosives. The influence of the burning rate and the physical properties of the explosive on the distance of the DDT is analyzed. Macek s model for DDT If the deflagration wave is initiated at the bottom of a closed tube, it propagates into the explosive and generates a stream of gas in the opposite direction to the wave motion. Since these gases stop at the bottom of the tube, their kinetic energy changes in internal energy. Increasing pressure in the gas is the source of compression waves which interact with the wave and accelerate the movement of the deflagration flame (Figure ). Moving faster and faster, the combustion front generates in turn compression waves in the unreacted explosive before the flame front, which creates a complex, unsteady flow. The interaction of compression waves creates a shock wave inside the material. This wave initiates reactions leading to generation of a detonation wave. The amplitude of the shock wave which initiates detonation depends on the type of explosive. SW CP FF Time CW Explosive Figure. Distance DDT in a distance-time diagram: CP combustion products, FF flame front, CW compression waves, SW shock wave.

4 W.A. Trzciński The mechanism of increase in pressure in the flame front depends on the type of fuel. Turbulence and possible diffusion of hot combustion products are the dominant phenomena in gases and porous solid explosives. In homogeneous explosives, these processes are absent, and the main mechanism of increase in pressure is accelerated by laminar burning [3]. The conditions necessary for the formation of a shock front in a burning explosive are strictly defined []. Firstly, it is necessary to limit the flow of the resulting combustion products to the rear. This provides increased pressure acting on the explosive and generates in it disturbances in the direction of propagation of the flame wave. The leading edge of these disturbances, according to the accelerating piston mechanism, transforms into a shock front. If expansion of combustion products is unrestricted, the pressure does not rise sufficiently to form a shock wave in the explosive. Secondly, the increase in pressure must be rapid. Since the compression waves are catching up and generate a shock wave inside the charge, the time needed to increase pressure from the baseline value to the value on the initiating shock wave must be lower than the time of passage of the sound wave by the charge. Macek [] proposed the first generation of experimental and theoretical methods for testing the DDT phenomenon in secondary explosives. He measured the time-dependence of the pressure in the combustion products of the explosive enclosed in a metal tube, and the time and place of the emergence of the detonation wave in the explosive. The dependence of the measured pressure p in the combustion products on the time t was approximated by the function: p = p e at () where p is the initial pressure and a is an experimentally determined constant. Macek assumed that the pressure in the combustion products is equal to the pressure in the explosive. Macek described continuous compression of the explosive by the modified Tait equation: n ρ c ρ p = () n ρ where ρ, ρ and c are mean density, initial density and velocity of sound in the explosive, and n is a constant. The movement of the compression waves in the explosive was calculated in [] by the method of characteristics. Macek assumed that n = 3, so that the

5 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives characteristics C+ in the compression wave are straight lines (Figure ). The point of intersection of the characteristics on which the difference in pressure is a few kbar, was assumed as the beginning of the formation of the shock wave of such an amplitude as to be able to boost the explosive to detonate. The initiated detonation wave travels through the material in the direction of propagation of the flame front (detonation) and in the opposite direction (retonation). Using such an approximation, Macek obtained good agreement between the calculated times and distances to detonation transition and measured values (time ~ μs, distance ~ cm) for the two high explosives diethylnitramine dinitrate (DINA) and pentolite (PETN/TNT 5/5) []. Macek also attempted to clarify the apparent exponential relationship between pressure and time behind the front, using a one-dimensional adiabatic model of burning. The dependence of the linear flame front speed U on the pressure p had the form: U = β p λ (3) where β and λ are constants. The author of work [] assumed that for the tested materials λ = and β = cm/sec kbar. The resulting time-dependence of the pressure matched equation () by a constant of proportionality, which takes into account the size of the flame front surface. Although Macek obtained quantitatively good agreement between the experimental data and the calculated results for two high explosives, according to the authors of work [] the theoretical model has weaknesses. Firstly, assuming that the pressures in the combustion products and the explosive are equal, Macek ignored the pressure drop at the deflagration wave front, which can be significant if the process reaches the deflagration CJ regime. According to the CJ theory, the pressure in unreacted material can be much higher than the pressure in combustion products. Secondly, the dependence of the linear burning velocity on pressure describes well the experimental dependence for only about kbar pressure. The actual relation between the burning rate and pressure in the process of an accelerating deflagration wave, which generates a detonation wave, is different from the burning rate law (3). The third weak point of Macek s analysis is the omission of momentum transferred to the reaction products at the deflagration wave front. In the case of a closed tube, the energy associated with the movement of particles in fact affects the acceleration of the flame front.

6 W.A. Trzciński Models of DDT taking into account the properties of the deflagration wave Figures and 3 show, on the pressure-volume and pressure-particle velocity planes, schematic Hugoniot curves for the reaction products of an explosive in the case of detonation and deflagration []. The upper curves and the CJ points correspond to detonation, while the lower curves and CJ points are related to the deflagration. The parameters of the CJ point correspond to a stationary combustion wave propagating in an explosive. detonation CJ Hugoniot adiabate of reaction products Pressure shock adiabate of explosive A deflagration initial state p, v CJ Specific volume Figure. Pressure-volume diagram for deflagration and detonation []. Pressure detonation CJ reaction produtcs unreacted explosive deflagration initial state CJ _ + Particle velocity Figure 3. Pressure-particle velocity diagram for deflagration and detonation [].

7 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 3 As previously stated, a deflagration wave initiated at the bottom of a closed tube moves along the tube and burns the explosive, generating the reaction products moving in the opposite direction. The reaction products are decelerated at the bottom of the tube and their kinetic energy transforms into internal energy. As a result, the pressure increases and compression waves propagate in the products, causing acceleration of the deflagration wave. The accelerating flame front, in turn, generates a compression wave in the explosive which leads to a complex unsteady flow ahead of the front. The initial state of the unreacted explosive is the compressed state. This complex movement associated with the accelerating deflagration wave is illustrated in Figure 4 using particle velocity-distance diagrams []. For simplicity, the deflagration wave front was presented as a discontinuity, although in reality it has a finite width. In Figures 4a and 4b, the mass velocity of the reaction products on the front of the deflagration wave is negative and the occurrence of a compression area is necessary to reach a value of zero at the bottom of the tube. Figure 4c shows a case in which the mass velocity drop at the deflagration wave front is exactly equal to the mass velocity of the compressed material before the front. In the next phase of acceleration of the deflagration wave (Figure 4d) mass velocity behind the wave front is positive and it is necessary to create an expansion area to satisfy the boundary condition at the bottom of the tube. However, when the compression zone behind the wave front no longer exists, then the mechanism of acceleration of the deflagration wave associated with the process of gas compression also vanishes. Although the exact solution of the complex unsteady flow generated by the deflagration wave in a closed tube is very difficult, there are two analytical models of the process of combustion wave acceleration which allow the upper and lower limits of combustion velocity to be estimated in a relatively simple way. For the detonation wave only the CJ point corresponds to a stable process, but in the case of the deflagration wave, the final state of the products can be identified by any point on the Hugoniot curve between points A and CJ (Figures and 3), namely the so-called weak deflagration branch. Point A corresponds to a hypothetical deflagration at constant pressure, during which the mass velocity does not change, but the CJ point is associated with deflagration, for which the combustion wave velocity, mass velocity and pressure reduction are the largest.

8 4 W.A. Trzciński compression (a) deflagration wave front compression compression (b) compression deflagration wave front Mass velocity compression (c) deflagration wave front u = compression (d) deflagration wave front decompression Figure 4. Distance from a closed tube bottom Particle velocity-distance diagrams for accelerating deflagration waves []. Taking into account the above-described mechanism of acceleration and deflagration wave properties, the authors of work [] have considered two analytical models describing the process of flame acceleration. In the first, proposed by Adams and Pack [4] to describe deflagration waves in gases, there is a compression area of the explosive terminated by the deflagration wave whose velocity is such that the momentum of compressed explosive is exactly equal to the momentum fall at the wave front and the mass velocity of reaction products

9 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 5 is zero (Figure 4c). In this model, there must be no rarefaction or compression waves in the reaction products to meet the boundary condition at the bottom of the tube. Since the change in particle velocity with distance behind the front of the combustion wave is zero, the pressure change in the time on the left border must also be zero. Although in this model the law of conservation of momentum is fulfilled, it is not physically correct because it is not possible in this model to accelerate the flame and there is no mechanism to increase pressure at the bottom of the tube. However, the model of Adams and Pack allows the lower flame velocity and the pressure drop at the front of the accelerating deflagration wave to be designated, since for each wave propagating with a faster speed there is a decrease in mass velocity at its front and it is possible to realize a mechanism of flame acceleration. In the second model it is assumed that deflagration is always the CJ wave with the greatest possible velocity of propagation and with the biggest drop in pressure at the front. This model allows the mechanism of flame acceleration in accordance with the scheme shown in Figure 4 to be described. It was used for the first time by Troshin to study the DDT process in gases [5]. An interesting state of deflagration, where the two models provide the same mass velocity and pressure in the reaction products, is the case shown in Figure 4c. For the CJ deflagration, this state corresponds to the case when the stage of flame acceleration and generation of new compression waves ends. Table shows, as calculated in [], the parameters of this state for TNT and pentolite. The values obtained for the compression wave, 95 kbar for TNT and 9.5 kbar for pentolite, are too big to be a criterion for DDT, because for pentolite its sensitivity to the shock wave is about kbar, and cast TNT can be initiated by a shock wave with an amplitude of 6 kbar. The transition from burning to detonation in pentolite is relatively easy, in TNT it is impossible. Thus, there is another criterion of flame acceleration which determines whether the transition from deflagration to detonation occurs in a cast explosive. Table. CJ deflagration wave parameters of zero-mass velocity behind the wave front [] Deflagration wave parameters Cast TNT Pentolite Shock pressure in explosive [kbar] CJ deflagration velocity [mm/μs] CJ pressure after deflagration [kbar] In the paper [] the discussed models of acceleration of flame were used to calculate the transition from burning to detonation in pentolite. To determine

10 6 W.A. Trzciński the time and distance of creation of a shock wave in the explosive ahead of the flame front, the method of characteristics was used. The time-dependence of pressure in the combustion products of pentolite for the pressure below kbar was described by Macek s relationship [] with the parameters p =.8 kbar, a =. /μs, and for the pressure above kbar by the relation given in [6]: p = (kbar)e.8856(t-.5) (4) where p is the pressure behind the front of the deflagration wave. Time t is expressed in microseconds. To describe the physical properties of the explosive, the modified Tait equation () with n = 3 was applied for isentropic compression. Using the relations valid along the set of negative characteristics C, the speed of sound c, mass velocity u, density ρ and the internal energy e of the compressed explosive in front of the deflagration wave were determined along the straight positive characteristics C + as a function of the pressure p : /3 3 p ρ = ρ + ρ c /3 3 p u = c + ρ c / 3 3 p u + c = c + ρ c p c ρ e e = (8) ρ ρ (5) (6) (7) where ρ, c are the parameters of equation (), e = c /. The laws of conservation of mass, momentum and energy flow at the deflagration wave front can be written in the following form: ρ (w u ) = ρ (w u ) (9) p + ρ (w u ) = p + ρ (w u ) ()

11 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 7 ( w u ) p ( w u ) p e q = e + + () ρ ρ where the index refers to the parameters behind the deflagration front, w denotes the velocity of the deflagration wave and q is the heat of reaction. The physical properties of the gaseous reaction products are described by a polytropic gas model, for which the relations are true: p = ρ k const () ρ D p CJ = (3) k + where k is the adiabatic exponent of the reaction products. In the paper [] its value was proposed to derive from the relationship linking the parameters of CJ detonation: ρ D p CJ = (4) k + where D is the detonation velocity. In the case of the Adams and Pack model, the mass velocity u behind the front of the combustion wave is zero. Then, from the conservation laws (9)- (), the deflagration wave velocity as a function of parameters marked with an index can be determined: w p k 3 ( )( ) ( ) p u w + k q + e + u k q + e + + = (5). u ρ ρ After determining the value of the deflagration wave velocity from equation (5) for a given pressure p, parameters p and ρ can be calculated from equations (9)-() assuming u =. In the paper [], the following data were used for pentolite: ρ =.67 g/cm 3, c =.43 mm/μs, q = 6 cal/g, k =.598. The results of calculations for the deflagration waves are presented in Table.

12 8 W.A. Trzciński Table. The results of calculations for deflagration waves for pentolite [] Adams and Pack model p [kbar] p [kbar] u [mm/μs] w [mm/μs] u +c [mm/μs] The calculated values of the deflagration wave velocity are approximately 3 times greater than the velocity obtained from equation (3) describing the dependence of the linear burning rate on pressure. It was observed in work [7] that the deflagration wave velocities in explosives are of the order of - mm/μs, which correlates well with the results in the fourth column of Table. This supports the critical remarks made by the authors of [] on Macek s extrapolation of experimental data obtained from measurements of the linear combustion velocity above the pressure of kbar. Data from Table also show that for pressures up to kbar, the pressures p and p differ slightly, which means that the model of Adams and Pack in this range is similar to the model of constant pressure deflagration. The motion trajectory of the combustion wave front can be calculated from the formula: x w () t w() τ = t dτ (6) where w(t) is determined from the dependence of w on p from Table, and the time-dependence of p, according to the relationship (4). Development of the combustion wave front and compression waves in pentolite on a plane x-t are shown in Figure 5.

13 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 9 4 Figure kbar c u + c 4 kbar shock formation kbar kbar kbar Characteristics of compression waves generated by the combustion of pentolite obtained from the Adams and Pack model. The characteristics of C + (u + c) intersect in an area approximately 3 mm from the bottom of the tube after a period of approximately 9 μs after the initiation of combustion of pentolite. It is assumed that at the point of intersection of the characteristics, the shock wave is generated which is able to initiate the explosive to detonation. The parameters of the transition of combustion into detonation estimated from the Adams and Pack model were compared with the results of calculations made on the basis of the CJ deflagration model applied to the explosive compressed by isentropic compression waves. At the flame wave front the CJ condition is fulfilled w = u + c, i.e.: ( w u ) p = k ρ. 7) Substituting relation (7) in equations (9)-(), the following expressions for the pressure and density behind the combustion wave front are obtained: ( ) ( ) ( q + e ) ( )( ) + p ρ p k ρ q + e ρ k q + e ( )( ) k q + e ( ) ( )( ) + p ρ k q + e + ρ p p + k p ρ k q + e + p ρ p = (8) ρ = (9)

14 3 W.A. Trzciński The deflagration velocity is then calculated from the relationship: p ρ w = u + k () ρ ρ The calculated results for pentolite are shown in Table 3 and in Figure 6. For a low initial pressure p, the pressure p behind the CJ deflagration wave front is much lower than the pressure estimated by the Adams and Pack model, and the flame wave velocity is nearly doubled. That means that the CJ deflagration wave simultaneously moves in the explosive at a greater distance. The time of emergence of the shock wave is about 8 μs and it is slightly shorter than the time calculated for the Adams and Pack model (around 9 μs), and the distance to shock formation increases (about 4 mm in Figure 6 to about 3 mm in Figure 5). Since the applied theoretical models allow the calculation of upper and lower limits for the deflagration to detonation transition, a critical time and distance of the transition will be bracketed by these results, if the isentropic compression of pentolite can be described by the equation () with n = 3. Table 3. The results of calculations for the deflagration waves for pentolite [] the CJ deflagration model p p u w [kbar] [kbar] [mm/μs] [mm/μs]

15 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 3 shock formation Figure kbar u + c c 4 kbar kbar kbar kbar Characteristics of compression waves generated by the combustion of pentolite for the CJ deflagration model. Modeling of the deflagration to detonation transition in primary explosives In the paper [] the models described above were used to calculate the distance and time of the passage from deflagration to detonation in the primary explosive lead azide. On the basis of experimental literature data, the timedependence of the pressure p for the deflagration wave in lead azide was obtained. To describe the physical properties of crystals of lead azide, the Tait equation and experimental data obtained from the relationship between pressure and mass velocity on the front of a shock wave were applied. The latter was used to calculate the mass velocity of compressed lead azide: p u = ρ c () where ρ = 4.6 g/cm 3, c =.5 mm/μs. In this case, the authors of [] did not give details of the DDT problem. Meanwhile, the use of data from the Hugoniot curve to describe the isentropic compression of the explosive is a coarse approximation. In addition, there is some ambiguity in the model associated with the calculation of the mass velocity of the compressed explosive by the use of the formula () and the velocity of sound, internal energy and density using the Tait model ().

16 3 W.A. Trzciński The calculated distance to the transition from deflagration to detonation in the lead azide crystal was about mm, which is in good agreement with experimental data obtained for lead azide crystals of different thicknesses. In crystals with a thickness of less than mm, there was no DDT phenomenon, whilst in thicker crystals this process was observed. The calculated combustion velocities for both models (Adams and Pack and CJ) are an order of magnitude smaller than the velocity of.7 mm/μs measured by Chaudhri and Field [8]. Despite the objections presented earlier in [] this model was also used to calculate the process of accelerating deflagration in pressed charges of lead azide. The following data were adopted: ρ = 3.4 g/cm 3, c =.3 mm/μs. In this case, the distance to the creation of the shock wave was.9 mm. Porosity causes reduced velocity of sound in the pressed charge and hence the intersection of the relevant characteristics takes place earlier. Due to the limited amount of experimental data for primary explosives, an attempt was undertaken to study the influence of the burning rate of these materials on the process of DDT in a typical primary explosive. In spite of previous limits, as in [], models developed for cast explosives were used. The following parameters for the primary explosive were assumed: ρ =.4 g/cm 3, c =. mm/μs. The proposed value of ρ corresponds to the density of pressed complex compounds tested in the Department of Explosives at the Military University of Technology [9], and the velocity of sound is close to the velocity in pressed lead azide []. The heat of combustion q = cal/g was derived from work []. It was assumed that the adiabatic exponent of the combustion products k is equal to 3.. The value of p from equation () was assumed to be. kbar. p [kbar] a = a = a = 3 a = Figure 7. t [µs] Time-dependence of pressure in the combustion products for different values of a (Eq. ()).

17 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 33 Initially, the influence of the burning rate of the explosive on the DDT process was investigated. The burning rate of a confined explosive is represented by the dependence of pressure on time (Eq. ). Figure 7 shows the change in pressure for different values of a. The pressure builds up to a value of about 4 kbar in the period from to 6 μs for changes in a from 5 to. To perform the calculations for a typical primary explosive the models of flame acceleration discussed in the preceding paragraph were used. Table 3 presents the parameters for the deflagration wave as a function of the pressure in the compressed explosive. The Adams and Pack and CJ models were used in the calculations. The range of changes of the pressure p corresponds to the amplitude of shock waves leading to detonation of primary explosives (e.g. 6 kbar in the case of lead azide []). As previously stated, the flame front velocity w and pressure p in the combustion products behind the wave front calculated for the models from Adams and Pack and CJ are the limiting values of these parameters for the deflagration wave in the tested explosive. Table 3. p [kbar] Results of calculations for deflagration waves in the primary explosive u u +c Adams and Pack model (u = ) CJ model [mm/μs] [mm/μs] p w p w [kbar] [mm/μs] [kbar] [mm/μs] Figures 8- show the x-t plane images of the transition from deflagration to detonation for various values of the exponent a in Eq. (). A slight difference is observed in the distances from the tube bottom to the points of creation of a shock wave for both models used. The distance of the deflagration to detonation transition in the tested porous explosive decreases from about to 6 mm when the parameter a is changed from to 5. This change in the parameter a means a significant increase in rate of pressure rise in the reaction products (Figure 7). The range of changes in distance of intersection of the characteristics obtained corresponds to the run distances to detonation in primary explosives. For example, the DDT in porous lead azide takes place over a distance of -3 mm depending on the duration of the loading pressure (. and 3.5 μs) with an amplitude of

18 34 W.A. Trzciński 6-9 kbar []. In turn, a complex primary explosive [Cu (C H 4 N 4 ) 3 ] (ClO 4 ) triggers the detonation of PETN, if the length of the charge of initiating explosive is at least 4.5 mm, which means that at this distance burning must develop into detonation [9] kbar 8 kbar 6 kbar shock creation 3 kbar c Adams and Pack model Figure 8. a = kbar 8 kbar c shock formation 3 kbar a = 6 kbar CJ model Characteristics of compression waves generated by the combustion of the primary explosive enclosed in a tube for the exponent a = kbar 3 kbar 8 kbar shock formation c 6 kbar a = Adams and Pack model,,5,,5,,5 3, 3,5 4, 4,5 5, 5,5 6, Figure kbar 8 kbar shock formation c a = 6 kbar 3 kbar CJ model,,5,,5,,5 3, 3,5 4, 4,5 5, 5,5 6, Characteristics of compression waves generated by the combustion of the primary explosive enclosed in a tube for the exponent a =. 4, 3,5 3,,5, kbar 8 kbar shock formation,5 3 kbar c, 6 kbar a = 3,5 Adams and Pack model,,,5,,5,,5 3, 3,5 4, 4, 3,5 3,,5, kbar 8 kbar shock formation,5 c, a = 3,5 3 kbar 6 kbar CJ model,,,5,,5,,5 3, 3,5 4, Figure. Characteristics of compression waves generated by the combustion of the primary explosive enclosed in a tube for the exponent a = 3.

19 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 35,5,,5, kbar shock formation 8 kbar c,5 3 kbar a = 5 6 kbar Adams and Pack model,,,5,,5,,5 3,,5,,5 kbar 8 kbar shock formation, c,5 a = 5 6 kbar 3 kbar CJ model,,,5,,5,,5 3, Figure. Characteristics of compression waves generated by the combustion of the primary explosive enclosed in a tube for the exponent a = 5. We can compare how the deflagration front velocity evolves over time in the primary explosive for the Adams and Pack model and the CJ model. The calculated results are shown in Figure. In the Adams and Pack model the front velocity is lower for the same time compared with the CJ model. However, it should be noted that in the first case, the mass velocity behind the wave front is zero (Figure 4c) and the whole released reaction energy is transmitted to the explosive. Therefore, generation of a compression wave of a given pressure ahead of the front requires a lower velocity than in the CJ deflagration model, in which a part of the energy is consumed to give a negative momentum of the gaseous combustion products (Figure 4b).,,,8,8 w [mm/µs],6,4, Adams and Pack model, t [µs] a = a = a = 3 a = 4 w [mm/µs],6,4,, t [µs] CJ model Figure. The change in time of the wave flame velocity in a primary explosive for various values of the exponent a. a = a = a = 3 a = 4 An analysis of the position of the intersection of the characteristics in Figures 8- shows that the distance from the charge bottom to the intersection point is changed relatively little for a high variation in the rate of pressure rise. In order to check how the properties of the primary explosive, density and

20 36 W.A. Trzciński velocity of sound, affect the location of this point, calculations for a material with a higher density were performed. For this calculation the complex material was assumed in the form of a crystal of density ρ =. g/cm 3 [9] and velocity of sound c =. mm/μs (as in crystals of lead azide). The burning parameters were identical to the parameters of the porous complex explosive tested earlier. Figure 3 shows wave images of the explosive compression by the accelerating deflagration wave for both the models used. The assumed dependence of pressure behind the front corresponds to the parameter a = (Figure 7). A comparison of the location of points of intersection of the characteristics in Figures 9 and 3 shows that the distance to detonation in a crystalline material with a high density is more than twice that in the case of a porous material kbar 6 kbar shock formation 8 kbar 3 kbar c c =, mm/µs ρ =, g/cm a = Adams and Pack model kbar 8 kbar shock formation c 3 kbar c =, mm/µs 6 kbar ρ =, g/cm a = CJ model Figure 3. Characteristics of compression waves generated by the combustion of the primary explosive of high density. Summary The application of the model of accelerating deflagration in cast explosives to a description of the transition from combustion to detonation in primary explosives allowed us to obtain results consistent with experimental data in quantitative terms. Thus, the estimation of the velocity of a flame wave front which leads to initiation of a primary explosive to detonation is possible. However, a more accurate assessment of the front velocity and the passage from burning to detonation will be achievable only on the basis of experimental measurements of the time-dependence of the pressure in reaction products behind the deflagration wave front. In the case of porous explosives more precise prediction of the DDT parameters requires the use of complex theoretical models, which better describe the phenomena of accelerating deflagration in these materials. This in turn leads to the need to solve the DDT problem by application of the numerical methods used in fluid dynamics, for example, the method of differential schemes.

21 Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives 37 References [] Tarver C.M., Goodale T.C., Shaw R., Cowperthwaite M., Deflagration-to- Detonation Transition Studies for Two Potential Isomeric Cast Primary Explosives, Sixth Symposium (International) on Detonation, Office of Naval Research, White Oak, 976, p. 3. [] Macek A., Transition from Deflagration to Detonation in Cast Explosives, J. Chem. Phys., 959, (), [3] McAfee J.M., Asay B.W., Campell A.W., Ramsay J.B., Deflagration to Detonation in Granular HMX, Ninth Symposium (International) on Detonation, Office of Naval Research, Portland, 989, p. 65. [4] McAfee J.M., Asay B.W., Bdzil J.B., Deflagration-to-Detonation in Granular HMX: Ignition, Kinetics, and Shock Formation, Tenth International Detonation Symposium, Office of Naval Research, Boston, 993, p. 76. [5] Parker G.R., Dickson P., Asay B.W., McAfee J.M., DDT of Hot, Thermally Damaged PBX 95 in Heavy Confinement, Fourteenth International Detonation Symposium, Coeur d Alene Resort,, ID 484. [6] Gifford M.J., Proud W.G., Field J.E., Observation of type II deflagration-todetonation transitions, CP6, Shock of Condensed Matter, American Institute of Physics,. [7] Field J.E., Walley S.M., Pround W.G., Balzer J.E., Gifford M.J., Grantham S.G., Greenaway M.W., Siviour C.R., The Shock Initiation and High Strain Rate Mechanical Charcterization of Ultrafine Energetic Powders and Compositions, Mat. Res. Soc. Symp. Proc., vol. 8, AA5.4., Materials Research Society, 4. [8] Walley S.M., Field J.E., Greenaway M.W., Crystal Sensitivieties of Energetic Materials, Materials Science and Technology, 6, (4), 4. [9] Pround W.G., Walley S.M., Willianson D.M., Collins A.L., Addiss J.W., Recent Trends in Research on Energetic Materials in Cambridge, Cent. Eur. J. Energ. Mater., 9, 6(), 67-. [] Bdzil J.B., Menikoff R., Son S.F., Kapila A.K., Steward D.S., Two-phase Modeling of Deflagration-to-Detonation Transition in Granual Materials: A Critical Examination Of Modeling Issues, Physics of Fluids, 999, (), [] Kapila A.K., Menikoff R., Bdzil J.B. Son S.F., Steward D.S., Two-phase Modeling of Deflagration-to-Detonation Transition in Granual Materials: Reduced equations, Physics of Fluids,, 3(), [] Narin B., Ozyoruk Y., Ulas A., Application of Parallel Processing to Numerial Modeling of Two-phase Deflagration-to-Detonation (DDT) Transition Phenomenon, in: Parallel Computation Fluid Dynamics (I.H. Tuncer et al.) 7; DOI:.7/ _5, Springer-Verlag Berlin Heidelberg 9. [3] Fogelzang A.E., Sinditski V.P., Egorshev V.Y., Serushkin V.V., Effect of Structure of Energetic Materials on Burning Rate, Proceedings of Symposium on Decomposition, Combustion, and Detonation Chemistry of Energetic Materials, Boston, 995, Materials Research Society, vol. 48, p. 5.

22 38 W.A. Trzciński [4] Adams G.K., Pack D.C, Some Observations on The Problem of Transition Between Deflagration and Detonation, Seventh Symposium (International) on Combustion, Butterworths, London, 959, p. 8. [5] Troshin Ya.K., The Generalized Hugoniot Adiabatic Curve, Seventh Symposium (International) on Combustion, Butterworths, London, 959, p [6] Price D., Wehner J.F., ThE Transition From Burning to Detonation in Cast Explosives, Combustion and Flame, 97, 9, 49. [7] Gibson R.W., Macek A., Flame Fronts and Compression Waves During Transition from Deflagration to Detonation In Solids, Eighth Symposium (International) on Combustion, Williams and Wilkins, Baltimore, 96, p [8] Chaudhri M.M., Field J.E., Deflagration in Single Crystals of Lead Azide, Fifth Symposium (International) on Detonation, Office of Naval Research, Pasadena, 97, p. 3. [9] Cudziło S., Complex compounds with explosive properties, Private Communication, 9. [] Fridrich M., Gilvez-Ruiz J.C., Klapotke T.M., Mayer, P. Weber B., Weigand J.J., BTA Copper Complexes, Inorganic Chemistry, 5, 44, [] Davies F.W., Zimmerschied A.B., Borgardt F.G., Avrami L., The Hugoniot and Shock Initiation Threshold of Lead Azide, Sixth Symposium (International) on Detonation, Office of Naval Research, White Oak, 976, p. 389.