1. Introduction
Solid propellant gas generators (SPGGs) are an emerging technology in the field of Halon alternative research, which has potential applications for fire suppression in dry bays and nacelles [
1,
2]. Fire suppression agents can be expelled under the reaction of solid propellant, including 1,1,1,2,3,3,3-heptafluoropropane, pentafluoroethane, or potassium carbonate. SPGGs have some exceptional properties, which include easy storage, rich in inert gases, and no liquid component, all of which significantly enhance its applicability in aviation firefighting.
According to the Next Generation Fire Suppression Technology Program [
3], green energetic materials (GEMs) are widely used as fuels in SPGGs. Examples of GEMs include 5-aminotetrazole (5AT), bitetrazole (BT), guanidinium bitetrazole (GBT), triaminoguanidinim nitrate (TAGN), and bis(aminotetrazolyl)tetrazine (BTATZ). Among the studied fuels, 5AT (as shown in
Figure 1) is the most readily available, and its commercial application prospect is the broadest. In contrast, many of the other fuels considered have limited commercial availability. As an environmentally friendly green energetic material, 5-amino-1h-tetrazole (5AT/CH
3N
5) has become the top candidate for use in solid propellants, rocket engines, airbags, and fire extinguishing areas. What makes 5AT extremely attractive is its surprisingly potential energy properties, e.g., high nitrogen content (82.3%), high formation enthalpy (208.7 kJ mol
−1), moderate mechanical, and thermal sensitivity. Once the mixtures of 5AT and oxidizers (strontium nitrate, potassium nitrate, etc.) are ignited, violent redox reactions will immediately take place with the generation of large amounts of gases and heat, creating desirable impetus. Undoubtedly, the pyrolysis behavior of thermal hazard substances usually plays a crucial role in burning rates and gas production performance for propellants, pyrotechnics, and explosives. Pyrolysis of 5AT is also regarded as the fundamental step to generate gaseous products which support fuel ignition, burning, and energy release. To elucidate the mechanism governing propulsive performance, it appears essential to analyze the thermal degradation behavior of 5AT.
Pyrolysis analysis is an effective method to study the properties of energetic materials, which is a crucial chemical conversion method that can help us to understand the thermal and combustion behaviors. Thermogravimetry (TG) technology can exhibit the generation of products in different phases in a given temperature range and the kinetic parameters can be obtained through the TG data, which is of great importance for the prediction and understanding of the pyrolysis and combustion characteristics of energetic materials and can further expand its application ranges. Thermogravimetry combined with Fourier transform infrared spectroscopy (TG-FTIR) is another method to investigate these properties, which can convert the inherent energy inside the molecule of energetic materials to the solid, gas and liquid products, which allow us determine the decomposition pathways of the energetic materials and aid their application, management, and computer simulations.
However, 5AT-based propellants have some limitations in their combustion and thermal decomposition behaviors that hinder their development [
4,
5]. These include the high exhaust temperature and the unstable burning rate. To improve its combustion performance, different chemical compounds have been incorporated into 5AT-based propellants [
6,
7]. Research has indicated that transition metal oxides (TMOs) can be employed in minute quantities in solid propellants to enhance their combustion characteristics. Examples of these TMOs include Fe
2O
3, ZnO, CuO, NiO, and TiO
2. The pyrolysis property of combustible materials is a fundamental aspect of material combustion [
8]. Therefore, the thermal stability of catalyzed 5AT is essential for its safe handling, storage, and usage. Thus, to study the pyrolysis process and thermal behavior of 5AT with the addition of additives is necessary for its practical application. The objective of this study is to examine how nano-sized NiO affects the thermal decomposition characteristics of 5AT.
In previous studies, the effects of TMO addition have been conducted on traditional propellant ingredients [
9,
10]. Prajakta et al. [
9] and Wang et al. [
10] conducted studies on the effect of two different catalytic agents containing nano-copper oxide and nano-sized NiO on the degradation behavior of AP-based propellants, respectively. Prajakta et al. [
9] found that the degradation temperature and the activation energy decreased with the addition of nano-copper oxide. Wang et al. [
10] proposed that the degradation temperature of AP decreased by 93 °C after adding 2% nano-sized NiO. Paletsky et al. [
11] studied the degradation kinetic characteristics of 5AT at high heating rates (~100 °C s
−1), and distinguished two routes of 5AT decomposition. Levchik et al. studied the energies of atomic bonding of 5AT tautomers during its decomposition process [
12]. A review of the literature reveals that the kinetic parameters, thermodynamic parameters, and thermal safety parameters of 5AT with catalysts have not been extensively studied.
This article presents a comprehensive analysis of the thermostability of pure 5AT and its mixture with nano-sized NiO. The study used thermogravimetry (TG), differential scanning calorimetry (DSC), and Fourier transform infrared (FTIR) spectroscopy to evaluate the thermal kinetics, thermodynamic parameters, and thermal safety. The dynamic thermal TG-DSC techniques were used to obtain the mass and heat flow variations with the increase in temperature. The weight-loss curves were used to evaluate the thermal kinetics and thermodynamic parameters, while the heat flow curves were used to evaluate the thermal safety. The 3D spectrogram was plotted using FTIR analysis to identify the evolved gas distribution during pyrolysis. In this study, NiO nanoparticles were investigated due to their catalytic effects on the degradation process of 5AT. This may have implications for the research field of solid propellant gas generators in the aviation industry.
2. Theoretical Considerations
An examination of the degradation kinetics and pyrolysis mechanism of 5AT was conducted using three iso-conversional methods (FWO, FR, and KAS). These approaches can estimate the activation energy as an outcome of the conversion process without assuming any decomposition model. This is provided that at least three different heating rates are used. Below are the detailed descriptions of each method. A general expression for degradation kinetics can be found in equation [
13].
where α,
T, and
β denote the conversion rate, temperature, and heating rate constant (
), respectively.
A,
E, and
f(α) are the kinetic triplets, which refer to the pre-exponential factor, the activation energy, and the kinetic function, respectively. R is the gas constant (8.314 J mol
−1 K
−1). In terms of weight loss, the conversion coefficient
α is calculated as follows [
14]:
Here, the instant mass of a sample at time t, Wt, represents the difference between its initial weight, W0, and its final weight, Wf, which is the weight of the sample after it has been completely decomposed.
2.1. FR Method [15]
As a result of Friedman’s method, activation energy is determined by TG data for different heating rates, and the mathematical equation may be expressed as follows:
By mathematically plotting ln(dα/dt) vs. 1/T, a linear relation indicating a slope equal to −E/R can be obtained.
2.2. FWO Method [16]
In terms of the FWO method, it is an integral approach that does not depend on the degradation mechanism. It can be expressed in logarithmic form as follows:
An estimate of E can be obtained by linear regression of lnβ vs. 1/T at a fixed conversion with a slope of 0.4567E/RT.
2.3. KAS Method [17]
The KAS method, similar to the FWO method, is an integral method. The following is the equation for KAS:
Plotting ln(β/T2) vs. 1/T allows the calculation of E for each degree of conversion values.
4. Results and Discussion
4.1. Thermal Degradation Behavior
TG and DTG profiles for 5AT–1 and 5AT–2 with a characteristic heating rate are shown in
Figure 2 and
Figure 3.
Table 2 lists the characteristic temperatures of different mass-loss stages at different conversion rates. Comparing
Figure 2 with
Figure 3, the trends of the mass loss curves were similar, while both the degradation of 5AT–1 and 5AT–2 exhibited four degradation stages, which can be seen from the four peaks in DTG. The pyrolysis process is attributed to the decomposition of crystalline water inside the 5AT molecular structure when the temperature is below 100 °C.
However, from
Table 2, it can be found that the characteristic temperatures at different conversions for 5AT–1 were obviously lower than those for 5AT–2, indicating that nano-sized NiO particles could reduce the decomposition temperature. This result implies that the thermostability of 5AT–2 has been slightly reduced with NiO addition, even though a similar decomposition process occurs.
4.2. Reaction Kinetics of Isoconversions
To calculate the activation energy, three heating rates and three iso-conversional kinetic methods (FR, FWO, and KAS approaches) were employed to analyze the thermal degradation characteristics. According to Equations (3)–(5), the linear regression parameters against 1/
T can be calculated, respectively using the three thermal kinetic methods (FR, FWO, and KAS). The activation energy for each conversion was calculated based on the linear relationship between variables.
Table 3 shows the activation energy results and the linear correlation coefficients for 5AT degradation.
From
Table 3, it can be seen that the square of the correlation coefficient at higher conversions were much smaller than 1 when the kinetic calculations were performed on the sample doped with nano-sized NiO, which may be responsible for the pyrolysis process being dominated by the reaction temperature at the beginning of the pyrolysis reaction. As the reaction continues, the catalytic effect of NiO plays an important role in the later pyrolysis process, which promotes the pyrolysis reaction and more gas products are rapidly generated; this causes the instability of the pyrolysis process in the later stage with an increased conversion rate. Therefore, the square of the correlation coefficient at higher conversions were much smaller due to the catalytic effect of NiO and the accelerated generation of gas products.
The mean value of E for 5AT–1 was 161.1–171.6 kJ mol−1, while that for 5AT–2 was 115.0–121.9 kJ mol−1 (varying with the kinetic methods). The derived E of 5AT–2 was much lower than that of 5AT–1, indicating that the nano-sized NiO significantly affected 5AT activity. As a result, the 5AT decomposition with NiO required a lower temperature and less energy.
The correlation relationships between the activation energy (
E) obtained from the three kinetic methods (FR, FWO, and KAS methods) and the conversion (
α) are shown in
Figure 4,
Figure 5 and
Figure 6. Generally speaking, the activation energy values of 5AT–2 were mostly lower than those of 5AT–1 in the conversion (
α) range from 0.1 to 0.9. Kinetic parameters are influenced by multiple factors, such as reaction time, temperature, reaction rate, conversion, etc. Moreover, the kinetic parameters also can be influenced by each other. As the thermal decomposition reaction progresses, the reaction time, temperature, and conversion are increasing, but the reaction rate is not a constant and it causes the activation energy and conversion rate to keep fluctuating, which further affects the instability of reaction.
The differences between the FR, KAS, and FWO methods, according to the results in
Figure 4,
Figure 5 and
Figure 6, may be attributed to the formula’s calculation principle. Specifically, the independent variables are all associated with temperature (1/T), but the dependent variables are different for these three methods. For the FR method, the dependent variable (ln(d
α/d
t)) is related to time (t) and conversion (
α), which is only directly influenced by the accuracy of the experimental recordings of time (t) and conversion (
α) accompanied by no approximations for the formulation’s reasoning and expression. For the KAS and FWO methods, both the formulations’ reasoning processes and expressions are adjusted by several approximations. Therefore, there is a difference between these three methods due to the formula’s calculation principle and reasoning process. However, no matter what method is used, the mean activation energy calculation results had similar variation trends and effectively reflected the catalytic effect of NiO.
4.3. Effects of Heating Rate on DSC Curves
Figure 7 and
Figure 8 show the DSC curves of 5AT–1 and 5AT–2 at 150–400 °C using four different heating rates. The decomposition temperature of both 5AT–1 and 5AT–2 was shifted to higher temperatures.
Three characteristic temperatures were obtained from the DSC profiles, i.e.,
T0,
Tonset, and
Tp; the determination method was same as that used in previous studies [
18].
Figure 9 shows the characteristics of the temperature changes (
To,
Tonset, and
Tp). With the growth of the heating rate, the characteristic temperatures increased. On the other hand, the
To,
Tonset, and
Tp of 5AT–2 were lower than those of 5AT–1.
4.4. Kinetic Parameters Derived from DSC Curves
The thermal decomposition of the two 5AT samples was investigated using the ASTM method E698 [
19] to determine the Arrhenius parameters. To obtain the parameter value of
A, the decomposition was deemed to follow the first-order kinetics. As a result, the first-order kinetics model fit the overall observed decomposition behavior of both 5AT samples over this temperature range, as evidenced using the linear regression method. The linear representation of ln(
β/
Tp2) vs. 1/
Tp made it possible to determine
E from the linear slope, as shown in
Figure 10. The logarithm of
A, lg
A, can be derived from the numerical relationship from ASTM E698:
As for the thermodynamic parameters of the activation reaction, the thermal kinetic parameters were also investigated. Shown in Equations (7)–(11), ∆
G≠, ∆
H≠, and ∆
S≠ are free energy, enthalpy, and entropy of the activation, respectively.
A summary of the characteristic kinetic parameters is presented in
Table 4. Compared to 5AT–1, 5AT–2 remained lower in activation energy. The positive values of ∆
G≠ imply that the propellants require heating to undergo exothermic decomposition. Moreover, the higher thermodynamic value of 5AT–1 indicates a higher heat demand for exothermic decomposition.
4.5. Thermal Safety Studies
When the heating rate approaches zero, i.e.,
, the values of
T00,
Te0, and
Tp0 can be calculated using Equation (10). Equation (11) [
18] can be used to estimate the temperature of self-accelerating decomposition.
where
b,
c,
d, and
e are coefficients.
In Equation (12),
TTIT is the thermal ignition temperature [
20]. Alternatively, the thermal explosion critical temperature (
Tb) can be determined by incorporating
E and
Tp0 into the equation below.
All the thermal safety values of the propellants are listed in
Table 5. We can see that both
TSADT and
TTIT of 5AT–1 were higher than those of 5AT–2, indicating that the NiO addition into 5AT resulted in a reduced resistance to heat and thermal safety of the 5AT propellant. Meanwhile, the critical value of
Tb for 5AT–2 was lower than that of 5AT–1, indicating that the process of thermal explosion for 5AT–2 occurred more readily.
4.6. Gas Product Analysis Using Fourier Transform Infrared Spectroscopy
Following the acquisition of all FTIR data, a typical 2D spectrogram for the two samples were plotted from the TG-FTIR. The temperature and wavenumber of 5AT–1 and 5AT–2 were used to identify the related evolved gas information. By comparing the two spectrograms, it was observed that the addition of nano-sized NiO generated new gases.
As can be seen in
Figure 11 and
Figure 12, two spectrum band peaks within 3400–3250 cm
−1 suggest that there were symmetrical stretching vibrations of primary amine (–NH
2) groups within the organic compounds. The absorption bands at 1650–1500 and 3336 cm
−1 could be attributed to the in-plane bending vibrations and stretching vibrations of secondary amine (–NH) groups. The presence of the absorption band at 2150–2000 cm
−1 was considered to indicate the existence of hydrazoic acid (HN
3). Additionally, another remarkable spectrum band at 2410–2200 cm
−1 was associated with the formation of a cyanogroup (–CN).
NiO exerts different catalytic effects on different groups. For the –NH and –CN groups, the initial evolving temperature is moved forward with the addition of NiO, and the evolving temperature range of HN3 is extended due to the presence of NiO. However, the terminal temperature of –NH2 is shortened when adding NiO.
Previous research had confirmed that the generation of CH2N2 can be detected at 250~500 °C during the pyrolysis process. Therefore, the –NH group is associated with the generation of CH2N2. The gas product of NH2CN can be detected when the temperature is over 500 °C. The –CN is most likely to belong to the characteristic group of NH2CN. The two spectrum band peaks within the 3400–3250 cm−1 range are the characteristic absorption peaks of HCN.
Based on the above analysis, there were no new gas products generated when adding nano-sized NiO, but different catalytic effects on different groups were obvious. Specifically, the initial evolving temperature of CH2N2 was detected in advance (about 50 °C earlier) with the presence of NiO. Moreover, the generation of HCN followed a similar variation trend. The evolving temperature range of HN3 was extended from 330 °C to 430 °C with the same initial temperature. In particular, the temperature range of NH2CN was shortened from 780 °C to 580 °C under the same terminal temperature.