Journal of Physics: Conference Series PAPER OPEN ACCESS Model of non-ideal detonation of condensed high explosives To cite this article: E B Smirnov et al 2016 J. Phys.: Conf. Ser. 774 012076 View the article online for updates and enhancements. Related content - Joining of tubes by gas detonation forming Vahid Jenkouk, Sandeep Patil and Bernd Markert - Experimental measurement of the scaling of the diameter- and thickness-effect curves for ideal, insensitive, and non-ideal explosives Scott I Jackson and Mark Short - An Experimental Study of Spherically Imploding Detonations Kunio Terao and Yoshiaki Furuya This content was downloaded from IP address 148.251.232.83 on 11/01/2019 at 13:30
Model of non-ideal detonation of condensed high explosives E B Smirnov, O V Kostitsin, A V Koval and I A Akhlyustin Federal State Unitary Enterprise Russian Federal Nuclear Center Academician Zababakhin All-Russian Research Institute of Technical Physics, Vasilieva 13, Snezhinsk, Chelyabinsk Region 456770, Russia E-mail: ewgeny smirnov@mail.ru Abstract. The Zeldovich Neumann Doering theory of ideal detonation allows one to describe adequately the detonation of charges with near-critical diameter. For smaller diameters, detonation velocity can differ significantly from an ideal value expected based on equilibrium chemical thermodynamics. This difference is quite evident when using non-ideal explosives; in certain cases, this value can be up to one third of ideal detonation velocity. Numerical simulation of these systems is a very labor-consuming process because one needs to compute the states inside the chemical reaction zone, as well as to obtain data on the equation of state of high-explosive detonation products mixture and on the velocity of chemical reaction; however, these characteristics are poorly studied today. For practical purposes, one can use the detonation shock dynamics model based on interrelation between local velocity of the front and its local curvature. This interrelation depends on both the equation of state of explosion products, and the reaction velocity; but the explicit definition of these characteristics is not needed. In this paper, experimental results are analyzed. They demonstrate interrelation between the local curvature of detonation front and the detonation velocity. Equation of detonation front shape is found. This equation allows us to predict detonation velocity and shape of detonation wave front in arbitrary geometry by integrating ordinary differential equation for the front shape with a boundary condition at the charge edge. The results confirm that the model of detonation shock dynamics can be used to describe detonation processes in non-ideal explosives. The theory of ideal detonation adequately describes detonation of charges with diameter, which is close to ultimate. With smaller diameters, the detonation velocity may sufficiently differ from ideal value, which is expected on the basis of equilibrium chemical thermodynamics. This difference is widely developed for insensitive high explosives and in separate cases, may achieve up to 30% of ideal detonation velocity [1]. This is related to the fact that under detonation of the charges with finite diameters, detonation front cannot be flat. Actually due to limited character of chemical transformation velocity directly behind the flat front of detonation wave (DW), the flow is subsonic. That is why, lateral unloading waves that appear under compressed matter expansion to the sides, penetrate primary the front, decreasing pressure on it and thereby decreasing its velocity, first of all near the charge. Thus shock front of detonation wave gets a shape that is protuberant towards the detonation propagation, and after completion of some section of development, it gets stationary shape. Numerical modeling of such systems is rather labor intensive, as far as it is necessary to calculate the states inside chemical reaction zone and it is also necessary to have information on equations of state for a mixture of the high explosive (HE) and explosion products (EP) and Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Published under licence by Ltd 1
Figure 1. Scheme of experimental studies of detonation development and propagation. chemical reaction velocity that currently are studied poorly. For practical purposes, we can use a model of detonation shock dynamics (DSD), which is based on interconnection of local front velocity with local curvature of detonation front. The interconnection depends on both equation of state for EP mixture and reaction velocity, but it is not necessary to explicitly specify these characteristics. The work analyzes the results of experiments that demonstrate the interconnection of local curvature of detonation wave front with the value of detonation velocity for this section of front. Experimental investigation into the process of detonation wave development and propagation in cyclic samples out of low-sensitivity HE was performedin therange of diameter change (1 10)d cr (d cr critical diameter of detonation). Initiation was made by a flat and divergent wave. The scheme of experiments is shown in figure 1. Streak camera recorded time profiles of detonation wave arrival to the edge surface of investigated cylindrical samples, which were recalculated into the profiles of detonation waves with regard for the velocity. Detonation velocity was registered by electro-contact technique. During experiments, we determined the profiles of detonation wave front at the section of development [2, 3] and after detonation establishment [2, 3] in cylindrical samples with diameter 15, 20, 40, 60 and 120 mm, as well as the values of detonation velocity corresponding to them. Determination of the shape and velocity in the process of detonation development is demonstrated by the results (figure 2) of measuring the velocity by electro-contact gauges that are installed at various removal r (figure 1) from the axis of cylindrical charge. Diameter of HE samples was 120 mm. Initiation was made by a flat wave. Electro-contact gauges were installed on the axis and at a distance 20 and 40 mm from the axis. The results in figure 2 show that the velocity of each section of the front in the process of detonation development is different. The character of the velocity change for different front sections also differs. The gauges installed on the axis of cylindrical part (r = 0 mm) show the establishment of a flat detonation wave [2,3], whenthevelocity passes via maximum valuecorrespondingto 1.01D st. Thegauges shifted from the axis to the periphery of the charge demonstrate the character of detonation establishment, which is observed with initiation of divergent wave [2, 3]. Such character of the velocity change for various sections of the front is explained by the curvature of the latter under the impact of lateral rarefaction waves. In the course of investigations, it was established that the character of detonation development does not depend on HE charge diameter. This is confirmed by the data obtained under flat- 2
Figure 2. Detonation development in high-explosive charge of 120-mm diameter under flat-wave initiation. wave initiation of cylindrical parts with diameter 60 and 120 mm that are shown in figure 3. The results on the plot are given in relative detonation velocity of cylindrical HE charges, as a function of the ratio of charge length to its diameter. The results given in figure 3 show that, in the beginning, detonation velocity rapidly increases withthecharge length x, achieving at x d/2 the valueof stationary detonation velocity for this charge diameter D st (d). At x 1.5d, detonation velocity achieves its maximum value exceeding D st (d) by 1%. After that, the lateral loading gradually decreases the detonation velocity, and at the charge length higher than 3d, the detonation velocity is practically equal to D st. After passing the section of development, which is as a rule of 3d (d charge diameter), detonation velocity and shape are established, and in further, for this diameter, remain constant [2, 3]. In order to determine the front curvature, the obtained profiles were approximated by analytical dependence of the following form: z = f(z), (1) where r radial coordinate, which is measured from the charge axis, z axial coordinate, which is calculated from the point of the detonation wave front arrival to the charge edge (figure 2). The curvature of any three-dimensional surface is determined from corresponding combination of its first and second spatial derivatives by orthogonal directions [4]. The first main curvature by radial coordinate is expressed as follows: k 1 (r) = z (r) [1+z (r) 2 ] 3/2. (2) In this charge geometry, it is assumed that the profile of detonation wave front is symmetrical relatively to the axis, and this permits to express the second main curvature: k 2 (r) = sinθ r = sinarctan(z (r)) r = z (r) r[1+z (r) 2 ] 1/2. (3) 3
Figure 3. Detonation development under flat-wave initiation of parts with diameter of 60 and 120 mm. Figure 4. Scheme of established detonation in cylindrical HE charge. For average surface curvature, we can write k(r) = 1 ( 2 z (r) [1+z (r) 2 ] 3/2 + z ) (r) r[1+z (r) 2 ] 1/2. (4) 4
Figure 5. Processing of DW front shape by functional dependence (6). Obtained equation is differential equation of the second order. Specifying the concrete type of dependence k(d), it is possible to get analytical solution for this equation. The authors of [1,5,6] got approximated solution for this equation in the form of logarithm of Bessel function: z(r) = ln(j 0 (r)), (5) where J 0 (r) Bessel function of the first kind, which for differencing convenience with accuracy up to a constant, may be expressed via trigonometric functions. Finally, for the function z(r), the following expression was proposed: z(r) = n 1 [ ( a i (ln cos η π r )]) i, (6) 2 R where r radial coordinate; R radius of cylindrical charge; a i, η adjustable coefficients. Experimentally obtained image of detonation wave front was digitized with resolution by radial coordinate r = 0.07 mm. After that, obtained profile (figure 5) was processed by analytical dependence (5) adding free term into it, which was to include possible front asymmetry. Unlike investigators [1,6], who used three first elements of series (6), in this work, the number of series terms was increased up to five, that permitted to expand the range of this function 5
Figure 6. Dependence of detonation velocity on the front curvature. applicability. The results of processing are shown in figure 5. The same figure shows deviations of experimental data from approximating function. Deviation between experimental and calculated values, except the edge zone, does not exceed 0.05 mm, with recalculation for times, it gives the value 0.007 µs. After substituting the coefficients, which were determined as the result of approximation, in expression (9), and in analytical expressions for the first and second derivatives of function z(r) with regard for formula (6), the function k(r) of DW front curvature was determined. From figure 4, it is seen that normal component of detonation velocity along radial coordinate r is related to stationary detonation velocity for this diameter D st (in cylindrical geometry) in the following way: D n (r) = D st cosθ(r) = D st cosarctan(z (r)) = D st [1+z (r) 2 ] 1/2, (7) where θ(r) angle between normal component of detonation velocity vector and axis of the charge by radial coordinate. Obtained are two functions D n (r) and k(r) in parametric form give connection of normal component of detonation velocity with curvature along DW profile (figure 6). Dependence D n (k), which is the connection of normal component of detonation velocity with curvature along radial coordinate, is in good agreement with the data that were obtained earlier for various diameters [2, 3]. Theory of DSD assumes that normal detonation velocity D is determined by a full front curvature k, and that edge angle θ (angle between front normal and HE edge) is unique for each combination of HE and lining material. Using this model, we can obtain conventional equations that describe stationary two-dimensional profiles of detonation profile in HE charges, having the shape of a plate, cylinder and a ring. 6
Let us consider detonation of HE cylinder lined by inert material, which is detonated in stationary regime and detonation velocity in the center is D st. For normal component of detonation velocity, we can write D n (k) = D st 1+(z (r) 2 ), (8) where function z(r) describes the front profile, and r radius. According to (4) expression for the full front curvature will be written as follows: ( 1 z (r) 2 [1+z (r) 2 ] 3/2 + z ) (r) D st r[1+z (r) 2 ] 1/2 = k 1+(z (r) 2 ). (9) Substituting y = z(r) and transformation give ( ) 1(1+y dr dy = D st 2k y 1+y 2 r 2 ) 3/2, (10) 1+y 2 ( ) 1(1+y dz dy = 2 y k D st y 1+y 2 r 2 ) 3/2. (11) 1+y 2 Integrating of equations (10), (11) gives the profile of stationary detonation front z(r) for specified D st. Initial point is y = 0, r = 0, z = 0, and finite point is y = tan(θ). Thus D st definitely determines the profile of stationary detonation front z(r), which corresponds to several combinations of the cylinder radius R 0 and lining material (edge angle θ). Therefore, we got conventional differential equations those describe the profiles of stationary two-dimensional detonation fronts in HE cylinder. It was revealed that one and the same stationary profile of detonation front corresponds to several combinations of lining material and determining charge size. Equation for the shape of detonation wave front permits to predict detonation velocity and the shape of wave front in arbitrary geometry by integrating conventional differential equation with boundary condition in the form of angle θ on the charge boundary. References [1] Kennedy D L and Drive G B 1998 Detonation, Proc., Symp. (Int.), 11th [2] Smirnov E B, Loboiko B G, Kostitsin O V, Belenovsky Y A, Prosvirnin K M and Kiselev A N 2010 Lavrentyev Readings on Mathematics, Mechanics and Physics (Novosibirsk: Institute of Hydrodynamics SB RAS) [3] Smirnov E B, Averin A N, Loboiko B G, Kostitsin O V, Belenovsky Y A, Prosvirnin K M and Kiselev A N 2011 Scientific and Technical Workshop Physics of Burning and Combustion devoted to K I Stchelkin 100-year Anniversary (Snezhinsk: RFNC-VNIITF) [4] Aminov Y A (ed) 1987 Differential Geometry and Topology of Curves (Moscow: Nauka) [5] Bdzil J B 1981 J. Fluid Mech. 108 195 226 [6] Lambourn B D and Swift D C 1989 Detonation, Proc., Symp. (Int.), 9th 7