Acta Phys. -Chim. Sin. ›› 2013, Vol. 29 ›› Issue (08): 1623-1631.doi: 10.3866/PKU.WHXB201305282


Method to Determine the Kinetic Parameters of the Autocatalytic Decomposition Reaction and Critical Rate of Temperature Rise of Thermal Explosion of Energetic Materials from DSC Curves

HU Rong-Zu1, ZHAO Feng-Qi1, GAO Hong-Xu1, YAO Er-Gang1, ZHANG Hai2, WANG Yao3, Chang Xiang-Yu3, ZHAO Hong-An4   

  1. 1 Science and Technology on Combustion and Explosion Laboratory, Xi’an Modern Chemistry Research Institute, Xi’an 710065, P. R. China;
    2 Department of Mathematics/Institute of data analysis and computation chemistry, Northwest University, Xi’an 710069, P. R. China;
    3 Institute for Information Science and System Science, Xi’an Jiaotong University, Xi’an 710049, P. R. China;
    4 College of Communication Science and Engineering, Northwest University, Xi’an 710069, P. R. China
  • Received:2013-04-08 Revised:2013-05-27 Published:2013-07-09
  • Contact: HU Rong-Zu
  • Supported by:

    The project was supported by the National Natural Science Foundation of China (21173163).


To evaluate the thermal safety of energetic materials (EMs) using the critical rate of temperature rise of thermal explosion (dT/dt)Tb and obtain basic data used to calculate (dT/dt)Tb, nine expressions to calculate (dT/dt)Tb for EMs were derived from Semenov thermal explosion theory and nine autocatalytic reaction rate equations, dα/dt=Aexp(-E/RT)α(1-α) (I), dα/dt=Aexp(-E/RT)(1-α)n(1+Kcatα) (II), dα/dt=Aexp(-E/RT)[αa-(1-α)n] (III), dα/dt=A1exp(-Ea1/RT)(1-α)+A2exp(-Ea2/RT)α(1-α) (IV), dα/dt=A1exp(-Ea1/RT)(1-α)m+A2exp(-Ea2/RT)αn(1-α)p (V), dα/dt=Aexp(-E/RT)(1-α) (VI), dα/dt=Aexp(-E/RT)(1-α)n (VII), dα/dt=A1exp(-Ea1/RT)+A2exp(-Ea2/RT)(1-α) (VII) and dα/dt=A1exp(-Ea1/RT)+A2exp(-Ea2/RT)α(1-α) (IX), using reasonable hypotheses. A method to determine the kinetic parameters in the nine autocatalytic decomposition reaction rate equations and (dT/dt)Tb of EMs when autocatalytic decomposition converts into thermal explosion from difference scanning calorimetry (DSC) curves at different heating rates (β) was presented. The onset temperature (Te0) corresponding to heating rate β→0, the critical temperature of thermal explosion (Tb) and the conversion degree (αb) corresponding to Tb used in the calculation of (dT/dt)Tb were obtained from analysis of DSC curves. The kinetic parameters of the autocatalytic decomposition reactions in Eqs.(I) and (VI), and Eqs.(II)-(V), (VII)-(IX) were estimated by the linear least-squares method and trust region approach, respectively. The values of (dT/dt)Tb for EMs were obtained from basic DSC data. Results show that (1) under non-isothermal DSC conditions, the autocatalytic decomposition reaction of nitrocotton (NC) (13.54% N) can be described by the apparent empiric-order equation dα/dt=1015.82exp(-170020/RT)(1-α)1.11+1015.82exp(-157140/RT)α1.51(1-α)2.51; (2)The value of for NC (13.54% N) when autocatalytic decomposition converts into thermal explosion is 0.103 K·s-1.

Key words: Critical rate of temperature rise, Non-isothermal DSC, Thermal explosion, Nitrocotton, Trust region approach


  • O641