INFLUENCE OF INITIAL DENSITY ON THE REACTION ZONE FOR STEADY-STATE DETONATION OF HIGH EXPLOSIVES

Similar documents
Theory and Detonation Products Equations of State for a New Generation of Combined Effects Explosives. Dr. Ernest L. Baker

Los Alamos NATIONAL LABORATORY LA-UR Double Shock Initiation of the HMX Based Explosive EDC-37

Lecture 7 Detonation Waves

EQUATION OF STATE FOR MODELING THE DETONATION REACTION ZONE

n v molecules will pass per unit time through the area from left to

ISOCHORIC BURN, AN INTERNALLY CONSISTENT METHOD FOR THE REACTANT TO PRODUCT TRANSFORMATION IN REACTIVE FLOW

Application of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures

Application of Proton Radiography to High Energy Density Research

Detonation of Gas Particle Flow

STUDY OF PHYSICAL-CHEMICAL TRANSFORMATIONS IN DETONATION WAVE BY THE ELECTRIC CONDUCTIVITY METHOD

Presented at the 2000 RGD Conference, Sydney Australia. The Simulation of Detonations Using a Monte Carlo Method

Simulation of sympathetic detonation by a CIP Eulerian code

ISENTROPIC COMPRESSION OF HIGH EXPLOSIVES WITH THE Z ACCELERATOR

AA210A Fundamentals of Compressible Flow. Chapter 13 - Unsteady Waves in Compressible Flow

EXPERIMENTAL STUDY OF SHOCK WAVE INTERACTING PLANE GAS-PLASMA BOUNDARY

Blast wave attenuation by lightly destructable granular materials

Optimization of Wall Velocity Measurements Using Photonic Doppler Velocimetry (PDV) *)

Numerical investigation of shock wave dense particles cloud interaction

Shock compression of precompressed deuterium

APPLICATION OF DETONATION TO PROPULSION. S. M. Frolov, V. Ya. Basevich, and V. S. Aksenov

INVESTIGATION INTO RESEARCH DATA ON HE SENSITIVITY TO IMPACT ( # 264)

Model of non-ideal detonation of condensed high explosives

Direct Simulation of Ultrafast Detonations in Mixtures

Laser Supported Detonation in Silica-based Optical Fibers

Constitutive Modelling of the Shock Behaviour of a Highly Porous Material

MEASUREMENT OF LOCAL POROSITY OF PARTICLES REINFORCED COMPOSITES WITH THE LASER OPTOACOUSTIC METHOD

KINETIC MODELING OF OXIDATION METHANE CONVERSION IN REGIME OF FILTRATION COMBUSTION WITH SUPERADIABATIC HEATING

About Some Features of a Magma Flow Structure at Explosive Volcano Eruptions

Numerical Analysis of the Deflagration to Detonation Transition in Primary Explosives

Numerical simulation of air blast waves

An Analytic Investigation of Hot-Spot Formation in Compressible Energetic Materials

Experimental and Computational Determination of Detonation Velocities

Chemical inhibiting of hydrogen-air detonations Sergey M. Frolov

Simulation of mixing of heterogeneous HE components

MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES

CRANFIELD UNIVERSITY DEFENCE COLLEGE OF MANAGEMENT AND TECHNOLOGY ENGINEERING SYSTEMS DEPARTMENT. PhD THESIS. Academic Year

EVIDENCE FOR THERMAL EQUILIBRIUM IN THE DETONATION OF HMX

On thermodynamic cycles for detonation engines

CHAPTER 5 TNT EQUIVALENCE OF FIREWORKS

A Vortex Model for Studying the Effect of Shock Proximity on Richtmyer-Meshkov Instability at High Mach Number

Thermodynamic evolution of phase explosion during high-power nanosecond laser ablation

Sound velocities in highly-oriented pyrolytic graphite shocked to 18. GPa: orientational order dependence and elastic instability

Effectiveness and Hydrodynamic Stability of Low-Entropy Laser-Driven Acceleration of Thin Foil.

A BURN MODEL BASED ON HEATING DUE TO SHEAR FLOW: PROOF OF PRINCIPLE CALCULATIONS. F. J. Zerilli, R. H. Guirguis, and C. S. Coffey

Chamber Development Plan and Chamber Simulation Experiments

THEORETICAL AND EXPERIMENTAL MESOSCALE STUDIES OF IMPACT- LOADED GRANULAR EXPLOSIVE AND SIMULANT MATERIALS

THE study of diffusion in gases and liquids always

INTERNATIONAL ATOMIC ENERGY AGENCY Division of Physical and Chemical Sciences Physics Section

Detonability of Ammonium Nitrate and Mixtures on Its Base

Raman Albedo and Deep-UV Resonance Raman Signatures of Explosives

Vibrational degrees of freedom in the Total Collision Energy DSMC chemistry model

Determining the Equivalent Explosive Effect for Different Explosives. Jerome B. Johnson U. S. Army Cold Regions Research and Engineering Laboratory

Shock Wave Boundary Layer Interaction from Reflecting Detonations

Phase-conjugation of high power molecular C02 and CO lasers radiation inside their active medium

Analitical criteria of homogeneous and heterogeneous detonation of liquid energetic materials

Fundamentals of explosive chemistry. Synopsis:

Modeling and Simulation of Fabrication, Transport, and Combustion of. Flow Environments

Determination of the Onset of Mechanical Destruction of Silver Azide Crystals Initiated by a Laser Pulse

Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture

Alternative speckle photography techniques for plastic deformation investigation

he science of high explosives is basically a coupling of chemistry and fluid mechanics. While each of these fields is in itself quite welldeveloped

Laser-Induced Explosion and Detonation in Gas-Particle and Gas-Droplet Mixtures

Density Field Measurement by Digital Laser Speckle Photography

THERMODYNAMIC THEORY OF NON-IDEAL DETONATION AND FAILURE

High-temperature and pressure aluminum reactions in carbon dioxide rich post-detonation environments

THERORY OF DETONATION WAVES 1 By J. VON NEUMANN

Surface Properties of EUVL. Laser Shock Cleaning (LSC)

Differential criterion of a bubble collapse in viscous liquids

Experimental Validation of Real Gas CO 2 Model near Critical Conditions

Detonation Structure

Target Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu

Molecular Dynamics Simulation of Diffusion of Gases and Liquids in Conditions of Phase Transition

Applying the Staggered Mesh Godunov (SMG) Method to Reactive Flows

Mass flow determination in flashing openings

Rarefied Gas Dynamics

The Reactions of Electrons and Holes Leading to the Metal Clusters Formation in Silver Halides: New Quantitative Characteristics 1

Theory and calibration of JWL and JWLB thermodynamic equations of state

UNIVERSITY OF CALGARY. Stability of Detonation Wave under Chain Branching Kinetics: Role of the Initiation Step. Michael Anthony Levy A THESIS

Preparation of the concerned sectors for educational and R&D activities related to the Hungarian ELI project

The role of diffusion at shear layers in irregular detonations

FRONTAL AND TILTED PDV PROBES FOR MEASURING VELOCITY HISTORY OF LASER- SHOCK INDUCED CALIBRATED PARTICLES

Size Segregation in the Brazil Nut Effect

FUNDAMENTALS OF CHEMISTRY Vol. II - Irreversible Processes: Phenomenological and Statistical Approach - Carlo Cercignani

An Introduction to Chemical Reactions, Gases, and Chemical Explosives

DEFLAGRATION AND DETONATION MODELING OF HETEROGENEOUS CONDENSED PI-IASE EXPLOSIVES. Joseph Ryan Peterson

Laser trigged proton acceleration from ultrathin foil

Ultrafast Shock-Induced Chemistry in Carbon Disulfide Probed with Dynamic. Kathryn E. Brown, Cynthia A. Bolme, Shawn D. McGrane, David S.

UNCLASSIFIED ADO ARMED SERVICES TECHNICAL INFORMAION AGENCY ARUNGIDN HALL STATIO ARLINGTFO 12, VIRGINIA UNCLASSIFIED

Ultrasonic Study of Binary Mixture of Acetone with Bromobenzene and Chlorobenzene at Different Frequencies


arxiv: v1 [cond-mat.other] 11 Sep 2008

SINGLE AND DOUBLE SHOCK INITIATION OF EDC37

Compressive Residual Stress Optimization in Laser Peening of a Curved Geometry

Strength effects in diamond under shock compression from 0.1 to 1 TPa

Physics of Explosions

Multicomponent diffusion in gases and plasma mixtures

Influence of the thickness absorbing film on the PETN ignition threshold by a laser pulse

Quantitative Analysis of Damage in PBX 9501 Subjected to a Linear Thermal Gradient

Cosmological 1st Order Phase Transitions: Bubble Growth & Energy Budget

On Analogy of 2D and 3D Combustible Mixture Flows

Transcription:

INFLUENCE OF INITIAL DENSITY ON THE REACTION ZONE FOR STEADY-STATE DETONATION OF HIGH EXPLOSIVES Alexander V. Utkin, Sergey A. Kolesnikov, Sergey V. Pershin, and Vladimir E. Fortov Institute of Problems of Chemical Physics RAS, Chernogolovka, Russia,4242, The laser interferometric system VISAR was used to investigate the detonation waves structure of pressed RDX, HMX, TNETB, and ZOX with different initial density. The experimental results are the surface velocity profiles of foils placed at the boundary between a HE sample and a water window. Critical initial densities ρ c at which the reaction zone structure changes crucially were found: Von Neumann spike was recorded if the density was less than the critical value, otherwise monotone pressure increase in the reaction zone was observed. The ρ c is equal to.72 g/cm,.84 g/cm,.56 g/cm, and.7 g/cm for RDX, HMX, TNETB and ZOX respectively. The results obtained for RDX demonstrate that the ρ c essentially depends on the sample structure and is determined not only by the HE particle size, but also by the pressing process. The results of this work can be explained by growth of the initial decomposition rate of an explosive with the increase of density if we assume that physicochemical transformations take place in a compression wave. INTRODUCTION According to the classical theory, the detonation wave consists of a shock jump and a chemical reaction zone, in which the pressure decreases and the matter expands, i.e. Von Neumann spike is shaped. The Chapman- Jouguet point, where particle velocity relative to the shock front is equal to the local sound velocity, ends a steady-state reaction zone. Inside of this zone a subsonic flow takes place and the conditions to supply chemical energy to the shock front are realized. Such flow structure allows us to establish the selection rule for detonation wave velocity. Numerous experimental data confirm the validity of this model for heterogeneous high explosives (HE). However, it was found 2 that in RDX and HMX at high initial density the pressure increases in the reaction zone and the spike does not form. Papers -5 contain indirect confirmation of existence of similar detonation wave structure in high density RDX and PETN: the authors stated that they did not observe the spike.

The detonation wave without Von Neumann spike does not correspond to the classical model. Moreover, it is not clear, whether the Chapman-Jouguet state will be reached, and what the selection rule of detonation velocity is in this case. To solve these key theoretical problems it is necessary () to understand if the observed phenomenon is the unique property of RDX and HMX, or it reflects the property common for powerful explosives, (2) to define the initial density at which the reaction zone structure changes crucially. For these purposes the experimental investigation of the reaction zone transformation under initial density increase in pressed RDX (C H 6 N 6 O 6 ), HMX (C 4 H 8 N 8 O 8 ), ZOX (C 6 H 8 N 0 O 6, Bis (2, 2, 2 Trinitroethyl - N - nitro) Ethylenediamine), and TNETB (C 6 H 6 N 6 O 4, Trinitroethyl trinitrobutyrate) was conducted. THE SCHEME OF EXPERIMENTS The scheme of experiments is shown in Figure. The detonation in the samples was initiated by a shock wave with amplitude ~4 GPa, created by the HE plane generator (). The diameter of the charges was 0 mm, the length changed from 40 up to 80 mm. Steadystate detonation wave was formed at these geometrical sizes. The wave profiles were registered by laser interferometer VISAR with 2 ns time resolution and an accuracy of the velocity determination was equal to ±0 m/s. The laser beam reflected from a 00-400 µm aluminum foil (2) placed between the charge and the water window (). As the result of the 2 experiment we have the velocity of the foil - water border, which represents all the details of the reaction zone structure in detonation wave. EXPERIMENTS AND RESULTS RDX and HMX RDX samples of different initial density were pressed from powder with a mean particle diameter of 80 µm. To obtain in the range of.50 g/cm to.74 g/cm, a small amount of acetone (less than wt%) was added into pressed RDX (RDX ). Higher density, up to.776 g/cm, was reached by increase of acetone quantity up to 0 wt% and by curing of the samples under pressure (RDX 2 ). TABLE. PARAMETERS OF EXPERIMENTAL SETUP FOR RDX. N, g/cm h Al, µm /l,mm/ Comments mm,5 400 0/45 spike, 70ns, RDX 2,5 200 0/45 spike, RDX,60 400 0/45 spike, 70ns, RDX 4,60 200 0/50 spike, RDX 5,69 200 0/40 spike, RDX 6,7 200 0/52 increase, 0 ns, RDX 7,72 200 0/40 spike, 40ns, RDX 2 8,776 200 0/40 spike, 0ns, RDX 2 9,68 200 0/45 spike, >40 ns, RDX 0,68 200 0/45 spike, >40 ns, RDX,74 200 0/40 increase, RDX 2,74 200 0/40 increase, RDX HE FIGURE. THE SCHEME OF EXPERIMENTS VISAR The experimental results for RDX are presented in Figures 2,. The initial density, foil thickness h Al, ratio of HE sample diameter to its length l, and comments concerning the reaction zone structure and characteristic time are given in Table. The

RDX,4 RDX : 2 =.776 g/cm =.72 g/cm 4 =.60 g/cm,2 2 =.5 g/cm 5 6 8 7 RDX : =.69 g/cm =.7 g/cm 0,20 FIGURE 2. THE VELOCITY OF AL FOR RDX AT LOW INITIAL DENSITY. determination accuracy of the characteristic time does not exceed ±5 ns. When the initial density is less than.72 g/cm (Figure 2) the velocity of the foil - water boundary decreases after the shock jump. Duration and amplitude of the spike are determined by Von Neumann spike in RDX. The subsequent velocity increase is caused by circulation of compression and rarefaction waves in Al foil. The attenuation of the velocity spike is recorded during its propagation through the foil. The spike amplitude decreases by 00-200 m/s as the foil thickness increases from 200 µm up to 400 µm. Therefore to find the exact value of the spike amplitude it is necessary to carry out experiments with thinner foils and to extrapolate their thickness to zero. But it was not the aim of this paper because we were interested in the qualitative change of the reaction zone structure when the initial density increases. Relatively thick foils can be used for that purpose. For initial densities of.5 g/cm and.60 g/cm (Figure 2) the velocity profiles show inflection at ~0 ns while the total spike duration is equal to 70 ns. That can reflect the two-step decomposition of RDX in detonation wave. Similar chemical reaction of TNT have been observed previously. 6 The fist rapid FIGURE. THE VELOCITY OF AL FOR RDX AT HIGH INITIAL DENSITY. reaction (80 ns) was interpreted as the formation of gas detonation products and the second stage over an additional 200 ns was attributed to diffusion controlled solid carbon particle formation. The velocity profiles inflection in Figure 2 may be the result of similar chemical reaction of RDX. The velocity spike duration changes insignificantly, but its amplitude notably drops as the density increases. When exceeds ρ c =.72 g/cm the situation changes crucially: instead of velocity decrease after the shock jump the monotonic increase is observed (curve 6 for =.7 g/cm in Figure ). It is necessary to note, that the particle velocity of explosion products increases with the increase of initial density except for critical density, where the abnormal velocity behavior is registered. For =.7 g/cm the particle velocity is the same as for.69 g/cm, whereas it has to increase by more than 50 m/s. It is possible to explain this by transition to underdriven detonation when Von Neumann spike disappears. Figure shows the experimental data for RDX 2. Here the velocity spike is registered for all initial densities. This result differs essentially from the result that one would

2,5 2,0,5 expect on the basis of the RDX data extrapolation to higher densities. The spike amplitude drops approximately twice when the increases from.72 up to.776 g/cm, and, by analogy with RDX, it is possible to expect that for RDX 2 ρ c.78 g/cm. This conclusion agrees with the data. 2 It can also be seen, that the particle velocity at.72 g/cm density (curve 7) is approximately 00 m/s higher than the velocity at.7 g/cm density for RDX, when the spike fades. That fact confirms that underdriven detonation is possible when the density exceeds ρ c. The obtained results demonstrate that the,4 fine RDX =,66 g/cm ρ0 =,48 g/cm ρ0 =,4 g/cm ρ0 =,2 g/cm ρ0 FIGURE 4. THE VELOCITY PROFILES FOR FINE RDX. h Al =200 µm. ρ c depends essentially on the sample structure and is determined not only by the RDX particle size, but also by the pressing process. It is known, 7 that at pressing many of the explosive particles are cracked and sheared. Pressing with a small quantity of acetone (RDX ) creates a lot of the potential centers of reaction, and the decomposition rate and explosive part reacting in a shock front increase. Therefore the detonation wave without Von Neumann spike is formed at >.72 g/cm. Pressing with a large quantity of acetone (RDX 2 ) gives the same density at smaller damage of RDX particles, which decreases the decomposition rate and the spike is recorded. It follows from these results that the reaction zone structure has to depend on both the initial density and the initial particle size of HE powder. To check this assumption the RDX samples were pressed from powder with a mean particle diameter of 5 µm. Figure 4 shows the experimental results for fine RDX. Maximum initial density did not exceed.66 g/cm, but from this figure we notice that the amplitude and duration of fine RDX are smaller than those of coarse RDX at the same density. It is conceivable that the critical density of fine RDX is below.72 g/cm. The most unusual result is the disappearance of,2 RDX RDX,2 9 =.69 g/cm 5 =.68 g/cm 2 0 6 FIGURE 5. EXPERIMENTAL DATA FOR RDX AT <ρ c. FIGURE 6. EXPERIMENTAL DATA FOR RDX AT ρ c.

,4 HMX,2 =.84 g/cm =.80 g/cm =.77 g/cm spike with the duration of ~0 ns is recorded at the =.77 g/cm. When the density is equal to ρ c =.84 g/cm the spike disappears, and the velocity remains constant after the shock jump. This result agrees well with the data given in paper, 2 where the pressure increase in the reaction zone was found at the initial density of.87 g/cm. TNETB and ZOX FIGURE 7. THE VELOCITY PROFILES FOR HMX. h Al =200 µm. Von Neumann spike in fine RDX at. g/cm, whereas it is recorded at densities both below and above this density. Reproducibility of experimental results is an important aspect for heterogeneous medium. In pressed HE the typical pore sizes change from one to few tens of microns, so one would be expect a marked difference between experimental data. But it follows from our experiments that it is not so. Figures 5, 6 show the results for initial density lower and higher than the critical density. A good reproducibility is the result of experimental setup. In the first place, the shock wave passes through the foil, the thickness of which significantly exceeds the pore sizes, and heterogeneity becomes negligible. In the second place, the foil velocity is recorded by laser interferometer which picks a signal off a spot with a diameter ~00 µm and averages the velocity too. Consequently, the replication of velocity profiles from shot to shot is good enough. The similar experiments were conducted for HMX. Samples were pressed from powder with a mean particle diameter of 50 µm. The quantity of acetone added in the HMX at pressing did not exceed 0.5 wt%. Figure 7 shows the experimental results. The velocity TNETB 8 is an explosive with low negative oxygen balance (-4.5 %) and density of a single crystal.89 g/cm. The detonation heat and sensitivity of TNETB are at the level of RDX. The experimental dependence of the detonation velocity D on the initial density can be approximated by linear relationship: D=.87+.92( -0.6), km/s, where the density is measured in g/cm. The TNETB was pressed without any solvents. TABLE 2. PARAMETERS OF EXPERIMENTAL SETUP FOR TNETB. h, N, /l, Comments g/cm Al µm mm/mm.48 200 0/42 spike, 50 ns 2.48 400 0/42 spike, 50 ns.5 200 0/42 spike, 40 ns 4.56 200 0/42 plateau 5.6 200 0/40 increase, 0 ns 6.6 200 0/80 increase, 0 ns 7.7 00 0/8 increase, >0 ns 8.7 200 0/8 increase, 0 ns The experimental results are presented in Figures 8, 9 and in the Table 2. At.48 g/cm initial density (curve ) the spike with ~50 ns duration is observed. The increasing of up to.5 g/cm (curve ) leads to decrease of spike amplitude, at =.56 g/cm (curve 4) the velocity spike disappears and behind the shock jump the parameters are constant. At further initial density increase, the velocity

2,7 2,7 TNETB 7 8 =.7 g/cm =.48 g/cm TNETB 5 6 4 =.6 g/cm =.5 g/cm =.56 g/cm 2, 2 2,,8 FIGURE 8. THE VELOCITY OF AL FOR TNETB. increases by ~00 m/s in ~0 ns time behind the jump. Good reproducibility of the experimental results is observed for the double length charge (curves 5, 6 in Figure 9). It means that steady-state detonation is realized in all the experiments. As well as for RDX, the abnormal change of particle velocity at critical density ρ c =.56 g/cm (Figure 9) was found: for ρ c the velocity is approximately 50 m/s lower than for.5 g/cm. It can be explained by transition to underdriven detonation at the moment of disappearance of Von Neumann ZOX ρ0=.7 g/cm,8 FIGURE 9. THE VELOCITY OF AL FOR TNETB NEARLY ρ c. spike. ZOX 9 is a powerful explosive completely balanced on oxygen, with the density of a single crystal.87 g/cm. The dependence of the detonation velocity D on the initial density can be approximated by the following relationship: D=4.0+.92( -0.6), km/s. Figure 0 shows the experimental results for ZOX. The disappearance of Von Neumann spike is watched at ρ c =.7 g/cm : after the jump the velocity remains practically constant. Whereas at smaller initial density the spike is registered. At =.5 g/cm its amplitude equals ~00 m/s, and duration is ~50 nc. The increase of initial density to.6 g/cm reduces amplitude and duration of the velocity spike approximately by 20 %. ρ0=.6 g/cm DISCUSSION 2,2 2,0 ρ0=.5 g/cm 0,20 FIGURE 0. THE VELOCITY OF AL FOR ZOX. h Al =200 µm. The results of this work confirm the possibility of detonation wave propagation without Von Neumann spike in a powerful HE. The reaction zone structure changes qualitatively at a critical initial density. It can be explained by growth of the initial decomposition rate of an explosive with the increase of density, if we assume that

Pressure a b 2 Volume FIGURE. PHASE TRAJECTORY IN THE PLANE PRESSURE VOLUME. physicochemical transformations take place in a compression wave. The influence of RDX charge structure on the critical density agrees well with this assumption. Besides when the pressure increases in the reaction zone, the final state of detonation products can be on the weak part of a detonation Huhoniot. Steady-state detonation wave without Von Neumann spike does not correspond to the classical model, which ignores the chemical reaction in a shock zone. Nevertheless, a strong coupling between the shock and reaction zones was found. 0 - In the book 2 the main attention is given to the fact, that after the shock jump the nonequilibrium HE Hugoniot is realized. Further there is a relaxation to an equilibrium HE Hugoniot and chemical reaction. Mathematical model of detonation phenomenon with the reaction in shock compression is analyzed in the book. It is proposed that the shock front thickness is produced by viscosity. Figure shows the 4 5 6 V 0 c results of these investigations in the plane pressure volume. According to the classical theory, shock jump compresses the material from the initial state to the point 2 on the Hugoniot of HE () and initiates the chemical reaction. Then the reaction products expand along the Rayleigh line (5) from the point 2 to the point 4 on the Hugoniot for the reaction products (). If we consider the HE decomposition inside of the shock front, the matter state changes steadily along the trajectory a, which does not coincide neither with the HE Hugoniot nor with the Rayleigh line. At low decomposition rate the steadystate detonation corresponds to the classical model. The Von Neumann spike amplitude decreases only and the part of explosive reacts in the shock front. However when the decomposition rate increases the structure of detonation wave changes significantly: the final state 4 on the detonation Hugoniot is reached as a result of monotonic pressure increase (trajectory b). At very fast decomposition rate the final state can belong to the weak part of the reaction product Hugoniot (point 6). Thus underdriven detonation can be realized. Curve c in figure shows the phase trajectory in this case. Attention must be paid to the next result: the selection rule of the detonation velocity gives a classical Chapman-Jouguet detonation when the final state belongs to the strong part of the reaction product Hugoniot (point 4), and it does not matter whether the Von Neumann spike exists or not. But if the underdriven detonation is realized the detonation velocity will be determined by the kinetics of chemical reaction and shock wave front, and not thermodynamics and gas dynamics only, as it takes place in the classical theory of detonation. ACKNOWLEDGEMENTS The work has been funded by Russian Fund for Basic Research, grant number 00-0-

208a and by the Program of Russian Academy of Science Physics and Chemistry of Extreme States of Matter. REFERENCES. Zeldovich, Ya.B. The Theory of Shock Waves and Introducing in Gas Dynamics, Academy Science of USSR, Moscow, 946, pp.45-98. 2. Ashaev, V.K., Doronin, G.S., Levin, A.D. Detonation Wave Structure of Condenced HE, The Physics of Combustion and Explosion, Vol. 24, No., 988, pp.95-0 (in Russian).. Voskoboinikov, I.M., Gogulya, M.F. Shock Front Illumination in Liquids nealy HE boudury, Chemical Physics, Vol., No. 7, 984, pp.06-04 (in Russian). 4. Tarver, C.M., Breithaupt, R.D., Kury, J.W. Detonation Waves in Pentaerythritol Tetranitrate, J.Appl.Phys., Vol. 8, No, 997, pp. 79-7202. 5. Lubyatinsky, S.N., Loboiko, B.G. Density Effect on Detonation Reaction Zone Length in Solid Explosives, In Shock Compression of Condensed Matter 997; AIP: New York, 998; Vol. 429, pp.74-746. 6. Kury, J.W., Breithaupt, R.D., Tarver, C.M. Detonation Waves in Trinitrotoluene, Shock Waves, Vol. 9, No 4, 999, pp. 227-27. 7. Burnside, N.J., Son, S.F., Asay, B.W., Skidmore, C.B., Particle Characterization of Pressed Granular HMX, in Shock Compression of Condensed Matter 997; AIP: New-York, 998; Vol. 429, pp.57-574. 8. Price, D. Dependence of Damage Effects upon Detonation Parameters of Organic High Explosives, Chemical Reviews, Vol. 59, No. 5, 959, pp. 80-825. 9. Haishan, D., Properties of Bis(2,2,2,- Trinitroethyl-N-Nitro) Ethylenediamine and Formulations Thereof, in X Symposium (International) on Detonation; AIP: New-York, 989, Vol. 429, pp.995-000. 0. Hirschfelder, J.O., Curtiss, C.F. Theory of Detonations. I. Irreversible Unimolecular Reaction, J. Chem. Phys., Vol. 28, No 6, 958, p. 0-47.. Linder, B., Curtiss, C.F., Hirschfelder, J.O. Theory of Detonations. II. Reversible Unimolecular Reaction, J. Chem. Phys., Vol. 28, No. 6, 958, pp. 47-5. 2. Dremin, A.N. Toward Detonation Theory, Springer, New York, 999, pp.2-20.. Williams, F.A. Combustion Theory, Addison-Wesley Publishing Company, London, 964, pp.94-2.