Research Article  Open Access
Trung Toan Nguyen, Duc Nhan Phan, Van Thom Do, Hoang Nam Nguyen, "Kinetic Analysis of the Thermal Decomposition of PolymerBonded Explosive Based on PETN: ModelFitting Method and Isoconversional Method", Advances in Materials Science and Engineering, vol. 2020, Article ID 9260818, 12 pages, 2020. https://doi.org/10.1155/2020/9260818
Kinetic Analysis of the Thermal Decomposition of PolymerBonded Explosive Based on PETN: ModelFitting Method and Isoconversional Method
Abstract
This work investigates kinetics and thermal decomposition behaviors of pentaerythritol tetranitrate (PETN) and two polymerbonded explosive (PBX) samples created from PETN (named as PBXPN85 and PBXPP85) using the vacuum stability test (VST) and thermogravimetry (TG/DTG) techniques. Both modelfree (isoconversional) and modelfitting methods were applied to determine the kinetic parameters of the thermal decomposition. It was found that kinetic parameters obtained by the modified Kissinger–Akahira–Sunose method (using nonisothermal TG/DTG data) were close to those obtained by the isoconversional and modelfitting methods that use isothermal VST data. The activation energy values of thermal decomposition reactions were 125.6–137.1, 137.3–144.9, and 143.9–152.4 kJ·mol^{−1} for PBXPN85, PETN, and PBXPP85, respectively. The results demonstrate the negative effect of the nitrocellulosebased binder in reducing the thermal stability of single PETN, while the polystyrenebased binder seemingly shows no adverse influence on the thermal decomposition of PETN in our presented PBX compositions.
1. Introduction
Pentaerythritol tetranitrate (also called PETN), which is a popular nitrate explosive, is extensively used in civilian and military applications due to its high energy characteristics [1, 2]. The main drawback of PETN is that it is highly sensitive to impact and friction pulses, and thus there may be explosion risks in the processes of production, storage, and especially in the transportation of PETN. Therefore, PETN is often used in the form of a polymerbonded explosive (also called PBX) [3–5] to overcome this shortcoming. The PBX sample based on PETN is a composite explosive in which the PETN particles are covered and bonded by the polymeric binder. The polymeric binder in PBX composition can reduce the sensitivity of explosive to mechanical pulses and make the charge compression easier and safer [6–9].
Like all other highenergy materials, the PBX composition is also an unstable thermodynamic system [10]. Compared to the single explosive such as PETN, the thermal decomposition of PBX is generally more complex because of the incompatible reactions between explosives and the polymeric binders [11–14], and the safety and performance of PBX sample may change during aging. Therefore, the behaviors and kinetics of thermal decomposition of PBX samples should be evaluated and clarified to assess the safety and the thermal hazard potential during storage. The decomposition kinetics and thermal behaviors of several composite explosives based on PETN have been carefully examined, and it was shown that the presence of the additive materials exerts powerful impacts on the thermal decomposition of these explosives [15–19].
In this study, we investigate the thermal decomposition kinetics of two PETNbased PBX samples named as PBXPN85 and PBXPP85. These results were compared with those of a single PETN to evaluate the influence of binders on the thermal decomposition behaviors of explosives. Thermal analyses were conducted by thermogravimetry (TG/DTG) and vacuum stability test (VST), where the VST method (developed by STABIL, Czech Republic) is commonly used to examine the chemical stability, compatibility, and shelflife of energetic materials in recent years [10, 14, 20]. Based on TG/DTG curves and VST results, we determined kinetic parameters by using isoconversional and modelfitting methods according to the ICTAC kinetics committee recommendations [21]. It is noteworthy that while modelfitting and isoconversional methods are both applicable to isothermal data, the modelfitting method is generally not applicable to nonisothermal techniques since the kinetic parameters obtained for nonisothermal data utilizing modelfitting technique are highly variable [22].
2. Materials and Methods
2.1. Materials
We imported Class1 PETN with a melting temperature of over 139.0°C from South Korea. Good quality nitrocellulose (NC), that has a nitrogen content of 12.20%, was supplied from the Vietnam factory. Polystyrene (PS), which has an average molecular weight of 80,000 u, was produced in our laboratory. The widely used plasticizer, dioctyl phthalate (DOP), was utilized to make the binder from the PS and NC (the ratios of DOP/NC and DOP/PS were 3/1 and 2/1, respectively). Besides, some suitable solvents such as toluene and ethyl acetate are used in this study.
2.2. Experimental Techniques and Methods
2.2.1. Preparation of PBX Samples
PS and NC were plasticized by DOP (the ratios of DOP/NC and DOP/PS are 3/1 and 2/1, respectively), then dissolved in the suitable solvents (toluene for PS and ethyl acetate for NC), maintaining the polymer/solvent ratio of 1/15 (w/v), and kept for approximately 5 hours to obtain a homogeneous solution. The sample preparation method was carried out by mixing PETN crystals with the binders based on PS and NC. The mixing process took place on the heating device at 70°C for 30 minutes, and the PETN crystals were covered by the binder layers [6, 23]. After that, the mixtures were vacuumdried at 60 to 70°C for 5 hours to remove the solvents. Table 1 shows the compositions of two proposed PBX samples.

2.2.2. Experimental Techniques
The thermal decomposition behaviors of PBX samples were determined by employing thermogravimetry (TG/DTG) and vacuum stability test (VST). The experimental conditions are as follows.
TG/DTG analysis was conducted using a NETZSCH STA 409 PC/PG (NETZSCHGerätebau GmbH, Selb, Germany). PBX samples (around 5.0 mg) were put in an aluminum oxide crucible and heated at various heating rates of 3, 5, 7, and 10 K·min^{−1}. The sample vial was heated from 30 to 350°C. Experimental heating processes were carried on under a dynamic nitrogen atmosphere with a flow rate of 30 mL·min^{−1}.
The VST test was carried out in a STABIL apparatus (OZM Research, Pardubice, Czech Republic). During this test, PBX samples were heated isothermally at different temperatures of 135, 140, 145, and 150°C under vacuum pressure (at most 0.672 kPa). A pressure sensor and a computer are used to record the relationship of the released gas volume versus heating time. The sample mass was 20.0 mg, and we performed the test from 5 to 10 days (depending on the temperature test), until the end of the thermal decomposition reaction.
2.2.3. Kinetics Analysis
The kinetics of thermal decomposition reactions in solids can be expressed as the following equation:where t is time (s); φ is the extent of conversion; T is the temperature (K); and is a function of conversion representing the reaction model. For the TG tests, the extent of conversion φ is defined as follows:where m_{i}, m_{t}, and m_{f} are the initial sample mass, sample mass at any time t, and sample mass at the end of the decomposition, respectively. For the VST test, the value of φ is expressed as [13, 24]where p_{max} is the pressure at the end of decomposition, p_{t} is the pressure at any time t, and p_{0} is initial pressure. Following the Arrhenius equation, is computed bywhere R, T, E, and A are the universal gas constant, absolute temperature, activation energy, and preexponential factor, respectively.
(1) ModelFitting Method. Modelfitting methods employed different forms of reaction model for finding the referred decomposition kinetic parameters; some of these reaction models are illustrated in Table 2.

For isothermal conditions, equation (1) can be used in the integral form as follows [22, 24]:where is the integrated form of the reaction model . Replacing each reaction model from Table 1 into equation (5), the rate constant can be computed from the slope of the straight line when plotting against t. Then, we can base on and formulations in the reaction model to determine . For each selected reaction model, is determined at different temperatures, and the kinetic parameters can be calculated by the Arrhenius equation:
(2) Isoconversional (ModelFree) Method. The reaction model (as in equation (1)) in the isoconversional method is considered as being independent of temperature. Under the isothermal condition of the VST test, combining equations (4) and (5), we obtain [13, 22, 24]where the activation energy (E_{φ}) is estimated from the slope of the straight line when plotting (−lnt_{φ,j}) against (1/T_{j}). Thus, the dependencies of activation energy (E_{φ}) on the conversion of the thermal decomposition (φ) could be obtained.
At a linear heating rate and under nonisothermal condition, equation (1) can be described as the following:and Kissinger, Akahira, and Sunose introduced the modelfree method using the following linear equation [22, 26]:
For each extent of conversion φ, the values of E and A can be found from the slope and the intercept when plotting against 1/T. The function can be determined according to the methodology reported in the literature [27, 28].
3. Results and Discussions
3.1. The TG/DTG Results
TG/DTG curves of PETN, PBXPN85, and PBXPP85 at different heating rates (i.e., 3, 5, 7, and 10 K·min^{−1}) were recorded and are shown in Figure 1.
(a)
(b)
(c)
For comparison, the characteristic parameters of TG/DTG curves of PETN, PBXPN85, and PBXPP85 are summarized in Table 3.
 
T_{ot}—onset temperature of decomposition exotherm peak, T_{i}—the initial temperature of the mass loss, and T_{p}—the peak temperature of the DTG curve. 
It has been shown that only a single decomposition step is related to the thermal decomposition of a single PETN (i.e., onestep reaction model) in all samples. The peak temperatures in DTG curves for the thermal decomposition of PETN, PBXPN85, and PBXPP85 were observed in a range of 186–201, 183–199, and 189–204°C, respectively. Generally, the decomposition residue of all samples increases with the increasing heating rate, and the residue mass of PBXPN85 and PBXPP85 may be less than that of PETN because of the presence of inert ingredient (i.e., DOP and PS) in their formulation.
As seen, the thermal decomposition of PETN occurred at a higher temperature than PBXPN85 but lower than that of PBXPP85. DTG curves in Figure 1 also demonstrated the role of the binder on the thermal decomposition behaviors of PETN. Thermal stability and decomposition temperature of single PETN are clearly affected by the binders. Specifically, the energetic binder (i.e., NCbased binder) reduces the decomposition temperature of single PETN, while the nonenergetic binder (i.e., PSbased binder) hardly affects the thermal stability of PETN. Interestingly, the PBXPN85 sample also shows an early decomposition. This may be caused by the exothermic decomposition effect of NC and the release of the gaseous products in the early decomposition of PETN in the binder matrix [13]. Contrarily, the decomposition temperature of PBXPP85 has a slight increase because this PBX sample has less amount of PETN (compared to single PETN) in its composition [8].
According to equation (2), the corresponding curve of the conversion φ versus temperature T was obtained and is presented in Figure 2. It can be observed that all of the curves have the form of the sigmoid function (i.e., the “S” shape curve), which is often obtained in the thermal analysis of high energetic material [29].
(a)
(b)
(c)
As recommended by the ICTAC kinetics committee [21, 30], the modified Kissinger–Akahira–Sunose (KAS) technique was applied for investigating the relationship between the activation energy and the conversion fraction. The fitting lines at each conversion fraction are presented in Figure 3 by plotting against 1/T.
(a)
(b)
(c)
The activation energy E_{φ} and the preexponential factor A of the thermal decomposition of all samples were evaluated based on the slopes and intercepts of the plots shown in Figure 3 and are expressed in Table 4. The dependence of activation energy on the conversion by the modified KAS method is shown in Figure 4.

The mean activation energy (E_{φ}) values (calculated according to the modified KAS method with the conversion range of 0.3–0.7 [13, 16]) of single PETN, PBXPN85, and PBXPP85 are 144.9, 137.1, and 152.4 kJ·mol^{−1}, respectively.
To verify these results obtained by the KAS method, the kinetic parameters were recalculated using NETZSCH Thermokinetic software according to the ASTM E698 [31]. The E_{φ} values of single PETN, PBXPN85, and PBXPP85 are 138.2, 129.7, and 150.9 kJ·mol^{−1}, respectively. The results were close to those obtained from the KAS method with the maximum difference of approximately 5%, which indicates the accuracy and reliability of measurements.
Compared to single PETN, the introduction of NC in PBXPN85 formulation leads to the decrease of its activation energy because the thermal stability of the nitrate group in NC is less than that in PETN [32, 33]. For PBXPP85 formulation, besides the reducing mechanical sensitivity of the explosive, an inert binder system (i.e., DOP/PS) acts as a stabilizer by trapping the active thermal decomposition product [29, 34], resulting in the higher activation energy and thus higher stability versus single PETN.
3.2. The VST Results
Results of TG/DTG were verified by VST studies, which were conducted at several isothermal temperatures (e.g., 125, 130, 135, and 140°C). From the relationship of the volume of released gas versus heating time in VST tests, the conversion fraction values (φ) for each time (t) were computed according to equation (3), and the curve of the conversion fraction (φ) versus heating time (t) is plotted for all samples in Figure 5.
(a)
(b)
(c)
It is obvious from Figure 5 that the obtained curves of PBXPN85, PBXPP85, and PETN obey the logarithmic trend, and the decomposition rate of all samples increased when the temperature was increased. In the case of PBXPN85 and PETN, it is apparent that the conversion grew rapidly during the first 2000 minutes and the conversion rate is low in the next 8000 minutes until the samples were completely decomposed. Similarly, the conversion rate of PBXPP85 can also be derived from Figure 5. In the case of PBXPP85, the conversion also increased rapidly in the first 2000 minutes at 135 and 140°C, while the conversion took longer time when tested at 125 and 130°C. By comparing the results of all samples, it was found that PBXPP85 sample took a longer time for completing the conversion than PBXPN85 and PETN samples. For example, at an isothermal temperature 125°C, the time for completing the conversion of PBXPN85, PETN, and PBXPP85 was 9990, 10485, and 10970 min, respectively. The differences in the thermal stability of three samples were assessed by the activation energy employing the modelfree (isoconversional) and modelfitting methods.
3.2.1. Isothermal Isoconversional Method from VST Data
Next, we utilize VST data to find kinetics parameters of PBX samples. By drawing versus 1/T curve at each temperature according to equation (7) for each conversion φ, the kinetic thermal decomposition parameters of three samples were calculated and are expressed in Table 5. The dependence of activation energy on the conversion is presented in Figure 6.

These E_{φ} values are in the same range of the E_{φ} values calculated from TG/DTG data using the modified KAS method, thus demonstrating the reliability of the VST test in evaluating the E_{φ} values of studied samples by the isothermal isoconversional method. The E_{φ} values of PETN, PBXPN85, and PBXPP85 were 137.3, 129.7, and 143.9 kJ·mol^{−1}, respectively. Besides that, the values of PETN, PBXPN85, and PBXPP85 were 31.9, 28.0, and 35.0 min^{−1}, respectively.
3.2.2. Isothermal ModelFitting Method from VST Data
The modelfitting method includes two fits: the first fit for determination of the model rate constant that best fits the VST data and the second fit for determination of several kinetic parameters such as preexponential factor A from the Arrhenius equation and the activation energy E_{φ} [13, 22]. Fourteen reaction models from Table 1 were selected for fitting the conversion of the thermal decomposition. Table 6 shows the results for the preexponential factor, the activation energy, and squares of the correlation coefficients obtained by plotting the integrated form versus heating time t at each isothermal temperature.

From Table 6, it was observed that the kinetic parameters (specifically, preexponential factor and activation energy), which are derived from the VST data utilizing the modelfitting method, are essentially independent of the reaction model used. It means that there is an insignificant difference between the kinetic parameters (i.e., E_{φ} and A) calculated using different model reactions. For selecting the best model fitting, the correlation coefficient R^{2} was used as a parameter. It was found that the best model that describes the thermal decomposition of PETN is contracting volume (R3) and the best model for PBXPN85 is 3D diffusionJander (D3) while the best model for PBXPP85 is secondorder (F2).
On the other hand, we found that the activation energy values calculated by the isothermal isoconversional method are approximate to those derived from the isothermal modelfitting method using the VST technique. To make a comparison, kinetic parameters obtained by different methods are presented in Table 7.

The values in Table 7 demonstrate the role of the binder on the initiation activation energy of PETN. The energetic binder based on NC reduces the activation energy (i.e., PBXPN85 has lower thermal stability than that of PETN), while the introduction of the nonenergetic material based on PS increases the activation energy of the explosive.
The ranges of E_{φ} values obtained by different methods were 125.6–137.1 kJ·mol^{−1} for PBXPN85, 137.3–144.9 kJ·mol^{−1} for PETN, and 143.9–152.4 kJ·mol^{−1} for PBXPP85. These ranges are narrow, proving the equivalence of these methods with each other. It is also noteworthy that E_{φ} of PETN obtained in this work is in good agreement with that reported by Pouretedal et al. [15] (135.1–136.9 kJ·mol^{−1}) and by Künzel et al. [16] (137.4 ± 3.4 kJ·mol^{−1}).
According to Brill et al. [35], the E_{φ} and values—determined from isothermal techniques in the same conditions—usually compensate each other. They found that all of pairs that lie in the regression line represent legitimate measurement and reliable results. We have plotted all E_{φ} values of PETN, PBXPN85, and PBXPP85 versus corresponding values (presented in Table 7) in Figure 7.
(a)
(b)
It can be observed the linear relationship of E_{φ} with values, which confirm the accuracy of the kinetic parameters obtained by the isoconversional method and modelfitting method using the VST test.
4. Conclusions
The thermal decomposition behavior and kinetics of PETN, PBXPN85 (85% PETN, 15% NCbased binder), and PBXPP85 (85% PETN, 15% PSbased binder) were examined by TG/DTG and VST techniques using modelfitting and isoconversional methods. The activation energy values computed by the nonisothermal technique (modified KAS method) were close to those calculated by modelfitting and isothermal isoconversional methods, ranging from 125.6–137.1 kJ·mol^{−1} for PBXPN85, 137.3–144.9 kJ·mol^{−1} for PETN, and 143.9–152.4 kJ·mol^{−1} for PBXPP85. The results indicate that the NCbased binder accelerates the thermal decomposition of PBX, thus reducing the activation energy of single PETN. One more valuable conclusion of this study is the influence of the PSbased binder. Specifically, the presence of the PSbased binder causes no adverse effect on the thermal stability of the main explosive in PBX composition.
Data Availability
The data used to support the findings of this study are included in the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
DVT gratefully acknowledges the support of Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant no. 107.022018.30.
References
 J. P. Agrawal, High Energy Materials: Propellants, Explosives and Pyrotechnics, John Wiley & Sons, Hoboken, NJ, USA, 2010.
 T. Urbanski, Chemistry and Technology of Explosives, vol. 2, Pergamon Press, New York, NY, USA, 1965.
 A. Elbeih, J. Pachman, S. Zeman, and Z. Akstein, “Replacement of PETN by bicycloHMX in semtex 10. Problem of mechatronic,” Armament, Aviation, Safety Engineering, vol. 1, pp. 7–16, 2010. View at: Google Scholar
 H.S. Kim and B.S. Park, “Characteristics of the insensitive pressed plastic bonded explosive, DXD59,” Propellants, Explosives, Pyrotechnics, vol. 24, no. 4, pp. 217–220, 1999. View at: Publisher Site  Google Scholar
 S. Moore, M. Schantz, and W. MacCrehan, “Characterization of three types of semtex (H, 1A, and 10),” Propellants, Explosives, Pyrotechnics, vol. 35, no. 6, pp. 540–549, 2010. View at: Publisher Site  Google Scholar
 A. Elbeih, S. Zeman, M. Jungova, and Z. Akstein, “Effect of different polymeric matrices on the sensitivity and performance of interesting cyclic nitramines,” Central European Journal of Energetic Materials, vol. 9, no. 2, pp. 131–138, 2012. View at: Google Scholar
 A. Elbeih, S. Zeman, and J. Pachman, “Effect of polar plasticizers on the characteristics of selected cyclic nitramines,” Central European Journal of Energetic Materials, vol. 10, no. 3, pp. 339–350, 2013. View at: Google Scholar
 N. T. Toan, P. D. Nhan, and V. H. Phuong, “Thermal decomposition and shelflife of PETN and PBX based on PETN using thermal methods,” Vietnam Journal of Science and Technology, vol. 56, no. 3, pp. 303–311, 2018. View at: Publisher Site  Google Scholar
 A. Elbeih, M. M. Mohamed, and T. Wafy, “Sensitivity and detonation characteristics of selected nitramines bonded by sylgard binder,” Propellants, Explosives, Pyrotechnics, vol. 41, no. 6, pp. 1044–1049, 2016. View at: Publisher Site  Google Scholar
 B. Vogelsanger, “Chemical stability, compatibility and shelf life of explosives,” CHIMIA International Journal for Chemistry, vol. 58, no. 6, pp. 401–408, 2004. View at: Publisher Site  Google Scholar
 Q.L. Yan, S. Zeman, and A. Elbeih, “Recent advances in thermal analysis and stability evaluation of insensitive plastic bonded explosives (PBXs),” Thermochimica Acta, vol. 537, pp. 1–12, 2012. View at: Publisher Site  Google Scholar
 Q.L. Yan, S. Zeman, P. E. Sánchez Jiménez, F.Q. Zhao, L. A. PérezMaqueda, and J. Málek, “The effect of polymer matrices on the thermal hazard properties of RDXbased PBXs by using modelfree and combined kinetic analysis,” Journal of Hazardous Materials, vol. 271, pp. 185–195, 2014. View at: Publisher Site  Google Scholar
 M. AbdElghany, A. Elbeih, and S. Hassanein, “Thermal behavior and decomposition kinetics of RDX and RDX/HTPB composition using various techniques and methods,” Central European Journal of Energetic Materials, vol. 13, no. 3, pp. 714–735, 2016. View at: Publisher Site  Google Scholar
 T. Nguyen, D. Phan, D. Nguyen, V. Do, and L. Bach, “The chemical compatibility and adhesion of energetic materials with several polymers and binders: an experimental study,” Polymers, vol. 10, no. 12, p. 1396, 2018. View at: Publisher Site  Google Scholar
 H. R. Pouretedal, S. Damiri, M. Ravanbod, M. Haghdost, and S. Masoudi, “The kinetic of thermal decomposition of PETN, Pentastite and Pentolite by TG/DTA nonisothermal methods,” Journal of Thermal Analysis and Calorimetry, vol. 129, no. 1, pp. 521–529, 2017. View at: Publisher Site  Google Scholar
 M. Künzel, Q.L. Yan, J. Šelešovský, S. Zeman, and R. Matyáš, “Thermal behavior and decomposition kinetics of ETN and its mixtures with PETN and RDX,” Journal of Thermal Analysis and Calorimetry, vol. 115, no. 1, pp. 289–299, 2013. View at: Publisher Site  Google Scholar
 K.S. Jaw and J.S. Lee, “Thermal behaviors of PETN base polymer bonded explosives,” Journal of Thermal Analysis and Calorimetry, vol. 93, no. 3, pp. 953–957, 2008. View at: Publisher Site  Google Scholar
 C.C. Huang, M.D. Ger, and S.I. Chen, “Study on thermal decomposition of Pentolites by modified vacuum stability apparatus and differential scanning calorimetry,” Propellants, Explosives, Pyrotechnics, vol. 17, no. 5, pp. 254–259, 1992. View at: Publisher Site  Google Scholar
 S. K. Bhattacharia, J. Nunley, and B. L. Weeks, “New insights into kinetics of PETN decomposition from the product and reactant point of view: an investigation with mass spectrometry and differential scanning calorimetry,” Thermochimica Acta, vol. 617, pp. 38–43, 2015. View at: Publisher Site  Google Scholar
 W. P. C. De Klerk, M. A. Schrader, and A. C. Van der Steen, “Compatibility testing of energetic materials, which technique?” Journal of Thermal Analysis and Calorimetry, vol. 56, no. 3, pp. 1123–1131, 1999. View at: Publisher Site  Google Scholar
 S. Vyazovkin, A. K. Burnham, J. M. Criado, L. A. PérezMaqueda, C. Popescu, and N. Sbirrazzuoli, “ICTAC kinetics committee recommendations for performing kinetic computations on thermal analysis data,” Thermochimica Acta, vol. 520, no. 12, pp. 1–19, 2011. View at: Publisher Site  Google Scholar
 S. Vyazovkin and C. A. Wight, “Isothermal and nonisothermal kinetics of thermally stimulated reactions of solids,” International Reviews in Physical Chemistry, vol. 17, no. 3, pp. 407–433, 1998. View at: Publisher Site  Google Scholar
 Q.L. Yan, S. Zeman, F.Q. Zhao, and A. Elbeih, “Nonisothermal analysis of C4 bonded explosives containing different cyclic nitramines,” Thermochimica Acta, vol. 556, pp. 6–12, 2013. View at: Publisher Site  Google Scholar
 A. Elbeih, M. AbdElghany, and T. Elshenawy, “Application of vacuum stability test to determine thermal decomposition kinetics of nitramines bonded by polyurethane matrix,” Acta Astronautica, vol. 132, pp. 124–130, 2017. View at: Publisher Site  Google Scholar
 A. Singh, G. Kaur, C. Sarkar, and N. Mukherjee, “Investigations on chemical, thermal decomposition behavior, kinetics, reaction mechanism and thermodynamic properties of aged TATB,” Central European Journal of Energetic Materials, vol. 15, no. 2, pp. 258–282, 2018. View at: Publisher Site  Google Scholar
 H. E. Kissinger, “Variation of peak temperature with heating rate in differential thermal analysis,” Journal of Research of the National Bureau of Standards, vol. 57, no. 4, pp. 217–221, 1956. View at: Publisher Site  Google Scholar
 G. I. Senum and R. T. Yang, “Rational approximations of the integral of the Arrhenius function,” Journal of Thermal Analysis, vol. 11, no. 3, pp. 445–447, 1977. View at: Publisher Site  Google Scholar
 V. Georgieva, D. Zvezdova, and L. Vlaev, “Nonisothermal kinetics of thermal degradation of chitosan,” Chemistry Central Journal, vol. 6, no. 1, p. 81, 2012. View at: Publisher Site  Google Scholar
 Q.L. Yan, S. Zeman, and A. Elbeih, “Thermal behavior and decomposition kinetics of Viton A bonded explosives containing attractive cyclic nitramines,” Thermochimica Acta, vol. 562, pp. 56–64, 2013. View at: Publisher Site  Google Scholar
 Q.L. Yan, S. Zeman, J. Š mand/Cro, R. Svoboda, and A. Elbeih, “Thermal behavior and decomposition kinetics of Formexbonded explosives containing different cyclic nitramines,” Journal of Thermal Analysis and Calorimetry, vol. 111, no. 2, pp. 1419–1430, 2012. View at: Publisher Site  Google Scholar
 ASTME69811, Standard Test Method for Arrhenius Kinetic Constants for Thermally Unstable Materials, American Society for Testing and Materials Philadelphia, Philadelphia, PA, USA, 1979.
 J. Kimura, “Chemiluminescence study on thermal decomposition of nitrate esters (PETN and NC),” Propellants, Explosives, Pyrotechnics, vol. 14, no. 3, pp. 89–92, 1989. View at: Publisher Site  Google Scholar
 C.C. Huang, M.D. Ger, Y.C. Lin, and SunIChen, “Thermal decomposition of mixtures containing nitrocellulose and pentaerythritol tetranitrate,” Thermochimica Acta, vol. 208, pp. 147–160, 1992. View at: Publisher Site  Google Scholar
 N. P. Loginov and S. N. Surkova, “Effectiveness of the action of stabilizers in explosive compositions under mechanical loading,” Combustion, Explosion, and Shock Waves, vol. 42, no. 1, p. 100, 2006. View at: Publisher Site  Google Scholar
 T. B. Brill, P. E. Gongwer, and G. K. Williams, “Thermal decomposition of energetic materials. 66. Kinetic compensation effects in HMX, RDX, and NTO,” The Journal of Physical Chemistry, vol. 98, no. 47, pp. 12242–12247, 1994. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Trung Toan Nguyen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.