Advances in Materials Science and Engineering
Volume 2020 (2020), Article ID 9260818, 12 pages
https://doi.org/10.1155/2020/9260818

Research Article

Kinetic Analysis of the Thermal Decomposition of Polymer-Bonded Explosive Based on PETN: Model-Fitting Method and Isoconversional Method

Trung Toan Nguyen,¹ Duc Nhan Phan,¹ Van Thom Do,² and Hoang Nam Nguyen³

¹Faculty of Special Equipment, Le Quy Don Technical University, 236 Hoang Quoc Viet, Ha Noi 100000, Vietnam
²Faculty of Mechanical Engineering, Le Quy Don Technical University, 236 Hoang Quoc Viet, Hanoi 100000, Vietnam
³Modeling Evolutionary Algorithms Simulation and Artificial Intelligence, Faculty of Electrical & Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam

Correspondence should be addressed to Hoang Nam Nguyen; nguyenhoangnam@tdtu.edu.vn

Received 25 August 2019; Revised 30 May 2020; Accepted 17 June 2020; Published 11 July 2020

Academic Editor: Alain Portavoce

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.

Abstract.

This work investigates kinetics and thermal decomposition behaviors of pentaerythritol tetranitrate (PETN) and two polymer-bonded explosive (PBX) samples created from PETN (named as PBX-PN-85 and PBX-PP-85) using the vacuum stability test (VST) and thermogravimetry (TG/DTG) techniques. Both model-free (isoconversional) and model-fitting 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 non-isothermal TG/DTG data) were close to those obtained by the isoconversional and model-fitting 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⁻¹ for PBX-PN-85, PETN, and PBX-PP-85, respectively. The results demonstrate the negative effect of the nitrocellulose-based binder in reducing the thermal stability of single PETN, while the polystyrene-based 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 polymer-bonded 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 high-energy 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 PETN-based PBX samples named as PBX-PN-85 and PBX-PP-85. 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 shelf-life 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 model-fitting methods according to the ICTAC kinetics committee recommendations [21]. It is noteworthy that while model-fitting and isoconversional methods are both applicable to isothermal data, the model-fitting method is generally not applicable to non-isothermal techniques since the kinetic parameters obtained for non-isothermal data utilizing model-fitting technique are highly variable [22].

2. Materials and Methods

2.1. Materials

We imported Class-1 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 vacuum-dried at 60 to 70°C for 5 hours to remove the solvents. Table 1 shows the compositions of two proposed PBX samples.

Table 1: The composition of two PBX samples based on PETN.

Composition
Content of materials (wt%)

PETN
DOP
PS
NC

PBX-PP
85.00
10.00
5.00
—

PBX-PN
85.00
11.25
—
3.75

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 (NETZSCH-Gerä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⁻¹. 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⁻¹.

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 mi, mt, and mf 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 pmax is the pressure at the end of decomposition, pt is the pressure at any time t, and p₀ is initial pressure. Following the Arrhenius equation, is computed by

where R, T, E, and A are the universal gas constant, absolute temperature, activation energy, and pre-exponential factor, respectively.

(1) Model-Fitting Method. Model-fitting methods employed different forms of reaction model for finding the referred decomposition kinetic parameters; some of these reaction models are illustrated in Table 2.

Table 2: Different reaction models generally applied to describe the thermal decomposition in solids [21, 22, 24, 25].

Model no.
Reaction model

1
Avarami-Erofe’ev (A2)

[−ln]1/2

2
Avarami-Erofe’ev (A3)

[−ln]1/3

3
Avarami-Erofe’ev (A4)

[−ln]1/4

4
Power law (P2)
2φ(1/2)
φ(1/2)

5
Power law (P3)
3φ(2/3)
φ(1/3)

6
Power law (P4)
4φ(3/4)
φ(1/4)

7
1-D diffusion (D1)
1/(2φ)
φ²

8
2-D diffusion (D2)

[ln] + φ

9
3-D diffusion-Jander (D3)
[32/3]/2
[1/3]²

10
First-order (F1)

−ln

11
Second-order (F2)
²

12
Third-order (F3)
³
/2

13
Contracting area (R2)
21/2

14
Contracting volume (R3)
32/3

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 (Model-Free) 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/Tj). Thus, the dependencies of activation energy (Eφ) on the conversion of the thermal decomposition (φ) could be obtained.

At a linear heating rate and under non-isothermal condition, equation (1) can be described as the following:

and Kissinger, Akahira, and Sunose introduced the model-free 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, PBX-PN-85, and PBX-PP-85 at different heating rates (i.e., 3, 5, 7, and 10 K·min⁻¹) were recorded and are shown in Figure 1.

(a)

(b)

(c)

Figure 1: TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 at several heating rates.

For comparison, the characteristic parameters of TG/DTG curves of PETN, PBX-PN-85, and PBX-PP-85 are summarized in Table 3.

Table 3: Kinetic parameters from TG/DTG data of PBX-PN-85, PBX-PP-85, and PETN.

Material
β (K·min⁻¹)
TG/DTG curve

Tot (°C)
Ti (°C)
Tp (°C)
Mass loss (%)
Residue (%)

PETN
3.0
165.5
130.2
186.0
98.3
1.7

5.0
169.6
132.0
191.1
97.2
2.8

7.0
172.7
133.2
196.0
96.3
3.7

10.0
178.9
135.6
200.8
96.0
4.0

PBX-PN-85
3.0
162.3
127.1
183.9
95.3
4.7

5.0
166.8
129.8
188.2
94.4
5.6

7.0
169.0
130.9
194.5
94.0
6.0

10.0
175.3
133.5
198.8
93.1
6.9

PBX-PP-85
3.0
168.6
136.2
189.4
93.7
6.3

5.0
172.5
138.0
194.2
93.1
6.9

7.0
177.0
139.9
197.4
91.5
8.5

10.0
183.2
142.5
203.6
91.9
8.1

Tot—onset temperature of decomposition exotherm peak, Ti—the initial temperature of the mass loss, and Tp—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., one-step reaction model) in all samples. The peak temperatures in DTG curves for the thermal decomposition of PETN, PBX-PN-85, and PBX-PP-85 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 PBX-PN-85 and PBX-PP-85 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 PBX-PN-85 but lower than that of PBX-PP-85. 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., NC-based binder) reduces the decomposition temperature of single PETN, while the nonenergetic binder (i.e., PS-based binder) hardly affects the thermal stability of PETN. Interestingly, the PBX-PN-85 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 PBX-PP-85 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)

Figure 2: φ versus T curves of PETN, PBX-PN-85, and PBX-PP-85 at several heating rates.

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)

Figure 3: Isoconversional plots at selected conversion degree of PETN, PBX-PN-85, and PBX-PP-85 according to the modified KAS method.

The activation energy Eφ and the pre-exponential 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.

Table 4: Kinetic parameters of PETN and PBXs by the modified KAS method.

φ
PETN
PBX-PN-85
PBX-PP-85

Eφ (kJ⋅mol⁻¹)
ln A, (min⁻¹)
R²
Eφ (kJ⋅mol⁻¹)
ln A, (min⁻¹)
R²
Eφ (kJ⋅mol⁻¹)
ln A, (min⁻¹)
R²

0.1
133.8
24.7
0.9955
124.5
21.3
0.9835
139.3
27.8
0.9880

0.2
138.3
26.3
0.9915
136.3
25.4
0.9905
150.4
32.5
0.9950

0.3
141.3
27.2
0.9850
136.5
24.9
0.9945
152.6
33.2
0.9915

0.4
145.2
28.4
0.9884
138.6
26.8
0.9964
152.1
32.4
0.9835

0.5
145.3
28.6
0.9935
137.9
24.8
0.9924
155.5
34.1
0.9980

0.6
147.1
29.9
0.9828
136.6
24.6
0.9920
153.2
32.2
0.9909

0.7
145.8
28.9
0.9845
135.8
24.7
0.9960
148.5
31.3
0.9945

0.8
141.4
27.5
0.9899
137.9
26.4
0.9949
147.6
29.8
0.9952

0.9
141.8
27.7
0.9854
127.2
22.6
0.9935
146.7
27.8
0.9960

Mean
144.9
28.6

137.1
25.2

152.4
32.6

Figure 4: The dependence of activation energy on the conversion by the modified KAS method.

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, PBX-PN-85, and PBX-PP-85 are 144.9, 137.1, and 152.4 kJ·mol⁻¹, 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, PBX-PN-85, and PBX-PP-85 are 138.2, 129.7, and 150.9 kJ·mol⁻¹, 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 PBX-PN-85 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 PBX-PP-85 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)

Figure 5: Conversion fractions (φ) versus heating time t for PBX-PN-85, PBX-PP-85, and PETN.

It is obvious from Figure 5 that the obtained curves of PBX-PN-85, PBX-PP-85, and PETN obey the logarithmic trend, and the decomposition rate of all samples increased when the temperature was increased. In the case of PBX-PN-85 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 PBX-PP-85 can also be derived from Figure 5. In the case of PBX-PP-85, 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 PBX-PP-85 sample took a longer time for completing the conversion than PBX-PN-85 and PETN samples. For example, at an isothermal temperature 125°C, the time for completing the conversion of PBX-PN-85, PETN, and PBX-PP-85 was 9990, 10485, and 10970 min, respectively. The differences in the thermal stability of three samples were assessed by the activation energy employing the model-free (isoconversional) and model-fitting 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.

Table 5: Kinetic parameter of PETN and PBXs by the isothermal isoconversional method.

φ
PETN
PBX-PN-85
PBX-PP-85

Eφ (kJ⋅mol⁻¹)
ln A, (min⁻¹)
R²
Eφ (kJ⋅mol⁻¹)
ln A, (min⁻¹)
R²
Eφ (kJ⋅mol⁻¹)
ln A, (min⁻¹)
R²

0.1
130.2
29.4
0.9835
123.9
24.1
0.9966
135.6
31.1
0.9767

0.2
137.4
31.8
0.9830
127.2
26.3
0.9921
138.4
32.5
0.9773

0.3
137.9
32.0
0.9962
128.6
27.3
0.9938
145.6
34.9
0.9897

0.4
139.2
32.5
0.9949
132.3
28.9
0.9939
147.3
35.7
0.9890

0.5
142.6
33.7
0.9889
132.3
29.2
0.9924
147.9
36.2
0.9883

0.6
143.6
34.0
0.9930
135.9
30.5
0.9967
147.5
36.4
0.9885

0.7
136.8
32.0
0.9854
137.5
31.1
0.9858
145.5
36.1
0.9900

0.8
136.6
31.9
0.9868
128.7
28.6
0.9885
145.3
36.5
0.9857

0.9
131.2
30.2
0.9831
121.4
26.6
0.9941
142.0
36.2
0.9879

Mean
137.3
31.9

129.7
28.0

143.9
35.0

Figure 6: The dependence of activation energy on the conversion using the VST technique.

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, PBX-PN-85, and PBX-PP-85 were 137.3, 129.7, and 143.9 kJ·mol⁻¹, respectively. Besides that, the values of PETN, PBX-PN-85, and PBX-PP-85 were 31.9, 28.0, and 35.0 min⁻¹, respectively.

3.2.2. Isothermal Model-Fitting Method from VST Data

The model-fitting 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 pre-exponential 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 pre-exponential factor, the activation energy, and squares of the correlation coefficients obtained by plotting the integrated form versus heating time t at each isothermal temperature.

Table 6: Kinetic parameters of PETN and PBXs calculated from the model-fitting method.

Model
PETN
PBX-PN-85
PBX-PP-85

Eφ (kJ mol⁻¹)
ln A, (min⁻¹)
R²
Eφ (kJ mol⁻¹)
ln A, (min⁻¹)
R²
Eφ (kJ mol⁻¹)
ln A, (min⁻¹)
R²

1
137.4
32.4
0.9818
122.0
27.6
0.9828
144.4
34.6
0.9896

2
138.5
32.5
0.9815
121.6
27.2
0.9830
144.5
34.4
0.9896

3
138.3
32.3
0.9827
121.9
27.1
0.9821
144.5
34.0
0.9892

4
139.1
33.6
0.9900
123.5
28.7
0.9846
146.4
35.8
0.9908

5
138.5
33.0
0.9890
122.8
28.2
0.9833
146.2
35.4
0.9907

6
134.1
31.4
0.9845
122.9
27.9
0.9838
146.1
35.1
0.9904

7
139.4
33.0
0.9838
123.9
28.3
0.9853
146.9
35.3
0.9906

8
139.1
32.7
0.9935
123.9
28.1
0.9844
146.9
35.1
0.9908

9
139.4
32.3
0.9862
123.7
28.5
0.9862
147.6
35.8
0.9905

10
139.8
33.3
0.9754
119.7
27.2
0.9855
146.7
35.2
0.9886

11
140.6
32.6
0.9860
125.6
28.0
0.9918
148.5
34.9
0.9913

12
139.7
34.4
0.9853
124.9
29.8
0.9853
148.2
36.6
0.9914

13
141.3
36.2
0.9846
124.8
31.1
0.9846
150.6
36.9
0.9940

14
142.6
38.2
0.9834
129.6
34.2
0.9834
152.4
40.0
0.9903

From Table 6, it was observed that the kinetic parameters (specifically, pre-exponential factor and activation energy), which are derived from the VST data utilizing the model-fitting 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² 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 PBX-PN-85 is 3-D diffusion-Jander (D3) while the best model for PBX-PP-85 is second-order (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 model-fitting method using the VST technique. To make a comparison, kinetic parameters obtained by different methods are presented in Table 7.

Table 7: Kinetic parameters of PETN and PBXs obtained by several methods.

Sample
TG/DTG using KAS method
Vacuum stability test (VST)

Isoconversional method
Model-fitting method

Eφ (kJ·mol⁻¹)
(min⁻¹)
Eφ (kJ·mol⁻¹)
(min⁻¹)
Eφ (kJ·mol⁻¹)
(min⁻¹)

PETN
144.9
28.6
137.3
31.9
139.1
32.7

PBX-PN-85
137.1
25.2
129.7
28.0
125.6
28.0

PBX-PP-85
152.4
32.6
143.9
35.0
150.6
36.9

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., PBX-PN-85 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⁻¹ for PBX-PN-85, 137.3–144.9 kJ·mol⁻¹ for PETN, and 143.9–152.4 kJ·mol⁻¹ for PBX-PP-85. 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⁻¹) and by Künzel et al. [16] (137.4 ± 3.4 kJ·mol⁻¹).

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, PBX-PN-85, and PBX-PP-85 versus corresponding values (presented in Table 7) in Figure 7.

(a)

(b)

Figure 7: The linear dependence of Eφ with of all samples obtained from the VST test.

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 model-fitting method using the VST test.

4. Conclusions

The thermal decomposition behavior and kinetics of PETN, PBX-PN-85 (85% PETN, 15% NC-based binder), and PBX-PP-85 (85% PETN, 15% PS-based binder) were examined by TG/DTG and VST techniques using model-fitting and isoconversional methods. The activation energy values computed by the non-isothermal technique (modified KAS method) were close to those calculated by model-fitting and isothermal isoconversional methods, ranging from 125.6–137.1 kJ·mol⁻¹ for PBX-PN-85, 137.3–144.9 kJ·mol⁻¹ for PETN, and 143.9–152.4 kJ·mol⁻¹ for PBX-PP-85. The results indicate that the NC-based 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 PS-based binder. Specifically, the presence of the PS-based 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.02-2018.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 bicyclo-HMX in semtex 10. Problem of mechatronic,” Armament, Aviation, Safety Engineering, vol. 1, pp. 7–16, 2010.

H.-S. Kim and B.-S. Park, “Characteristics of the insensitive pressed plastic bonded explosive, DXD-59,” Propellants, Explosives, Pyrotechnics, vol. 24, no. 4, pp. 217–220, 1999.

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.

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.

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.

N. T. Toan, P. D. Nhan, and V. H. Phuong, “Thermal decomposition and shelf-life of PETN and PBX based on PETN using thermal methods,” Vietnam Journal of Science and Technology, vol. 56, no. 3, pp. 303–311, 2018.

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.

B. Vogelsanger, “Chemical stability, compatibility and shelf life of explosives,” CHIMIA International Journal for Chemistry, vol. 58, no. 6, pp. 401–408, 2004.

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.

Q.-L. Yan, S. Zeman, P. E. Sánchez Jiménez, F.-Q. Zhao, L. A. Pérez-Maqueda, and J. Málek, “The effect of polymer matrices on the thermal hazard properties of RDX-based PBXs by using model-free and combined kinetic analysis,” Journal of Hazardous Materials, vol. 271, pp. 185–195, 2014.

M. Abd-Elghany, 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.

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.

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 non-isothermal methods,” Journal of Thermal Analysis and Calorimetry, vol. 129, no. 1, pp. 521–529, 2017.

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.

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.

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.

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.

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.

S. Vyazovkin, A. K. Burnham, J. M. Criado, L. A. Pérez-Maqueda, C. Popescu, and N. Sbirrazzuoli, “ICTAC kinetics committee recommendations for performing kinetic computations on thermal analysis data,” Thermochimica Acta, vol. 520, no. 1-2, pp. 1–19, 2011.

S. Vyazovkin and C. A. Wight, “Isothermal and non-isothermal kinetics of thermally stimulated reactions of solids,” International Reviews in Physical Chemistry, vol. 17, no. 3, pp. 407–433, 1998.

Q.-L. Yan, S. Zeman, F.-Q. Zhao, and A. Elbeih, “Noniso-thermal analysis of C4 bonded explosives containing different cyclic nitramines,” Thermochimica Acta, vol. 556, pp. 6–12, 2013.

A. Elbeih, M. Abd-Elghany, 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.

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.

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.

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.

V. Georgieva, D. Zvezdova, and L. Vlaev, “Non-isothermal kinetics of thermal degradation of chitosan,” Chemistry Central Journal, vol. 6, no. 1, p. 81, 2012.

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.

Q.-L. Yan, S. Zeman, J. Š mand/Cro, R. Svoboda, and A. Elbeih, “Thermal behavior and decomposition kinetics of Formex-bonded explosives containing different cyclic nitramines,” Journal of Thermal Analysis and Calorimetry, vol. 111, no. 2, pp. 1419–1430, 2012.

ASTM-E-698-11, 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.

C.-C. Huang, M.-D. Ger, Y.-C. Lin, and Sun-IChen, “Thermal decomposition of mixtures containing nitrocellulose and pentaerythritol tetranitrate,” Thermochimica Acta, vol. 208, pp. 147–160, 1992.

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.

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.