Thermal Decomposition Behavior and Thermal Safety of Nitrocellulose with Different Shape CuO and Al/CuONanothermites.

Bamboo leaf-like CuO(b) and flaky-shaped CuO(f) were prepared by the hydrothermal method, and then combined with Al nanoparticles to form Al/CuO(b) and Al/CuO(f) by the ultrasonic dispersion method. The phase, composition, morphology, and structure of the composites were characterized by X-ray powder diffraction (XRD), transmission electron microscopy (TEM), scanning electron microscopy (SEM), and energy scattering spectrometer (EDS). The compatibility of CuO, Al/CuO and nitrocellulose (NC) was evaluated by differential scanning calorimetry (DSC). The effects of CuO and Al/CuO on the thermal decomposition of NC were also studied. The results show that the thermal decomposition reactions of CuO-NC composite, Al/CuO-NC composite, and NC follow the same kinetic mechanism of Avrami-Erofeev equation. In the cases of CuO and Al/CuO, they could promote the O-NO2 bond cleavage and secondary autocatalytic reaction in condensed phase. The effects of these catalysts have some difference in modifying the thermolysis process of NC due to the microstructures of CuO and the addition of Al nanopowders. Furthermore, the presence of Al/CuO(f) can make the Al/CuO(f)-NC composite easier to ignite, whereas the composites have strong resistance to high temperature. Compatibility and thermal safety analysis showed that the Al/CuO had good compatibility with NC and it could be used safely. This contribution suggests that CuO and Al/CuO played key roles in accelerating the thermal decomposition of NC.


Introduction
Nanosize metal oxides exhibit excellent electrical, optical, magnetic, and catalytic properties, because they have high specific surface area and surface energy, and active sites. Therefore, the preparation and application of nanosize metal oxides have attracted widespread attention. Recently, nanosize CuO has been extended studies for the potential applications in ion batteries [1,2], gas sensors [3], catalysts [4,5], magnetoelectric effects [6], etc. It has become the typical representative of nanosize transition metal oxides. In the field of explosives and solid propellants [7], nanosize CuO, as an important catalyst, has been applied for promoting combustion. There are a lot of reports regarding the catalytic effects of nanometal oxides to the main components of solid propellants [8]. CuO can decrease the decomposition peak temperature and activity energy of ammonium perchlorate (AP), and increase the decomposition reaction rate and releasing heat of AP [9,10]. When adding 2.3. Preparation of Flaky-shaped CuO(f) and Al/CuO(f) 0.47 g of CuCl 2 ·2H 2 O was dissolved in 20 mL deionized water, and then 10 mL NaOH solution (1 mol·L −1 ) was added dropwise. After stirring for 30 min, the mixture was poured into a hydrothermal reaction vessel and reacted at 180 • C for 16 h. After the reaction is completed and natural cooling afterwards, the precipitate was separated and then washed with anhydrous ethanol and deionized water for several times. Subsequently, the precipitate was dried in an oven at 60 • C. Finally, the flaky-shaped CuO(f) was obtained.
The Al/CuO(f) nanothermite was prepared by the ultrasonic dispersion method after the CuO(f) were mixed up with Al nanopowders according to the mole rate of 1.33:1 for Al:CuO [36,37].

Samples Characterization
An FEI Quanta 400 (FEI Co., Hillsboro, OR, USA) field-emission environment scanning electron microscope (SEM) characterized the morphology and structure of the sample. The acceleration voltage was 30 kV and current was 4 A. OXFORD INCAIE350 energy-dispersive X-ray spectroscopy (EDS) from OXFORD Instruments INC (Oxford, UK) was used to roughly examine the composition of the sample. The discharge voltage was 4-10 kV and the distance between the electrodes was exactly 1 mm. The morphology and size of the sample were investigated via transmission electron microscopy (TEM) and high-resolution TEM using HITACHI H-7650B (Tokyo, Japan). The operating voltage was 80 kV and resolution was 0.2 nm (lattice image). The purity and phase structure of the sample were confirmed by powder X-ray diffraction (XRD) analysis on a Rigaku D/MAX-3C (Tokyo, Japan) X-ray powder Diffractometer. The radiation source was Cu K α (λ = 1.5418 Å) at 40 kV and 40 mA, and the range of 2θ was 10 • to 80 • .

Measurement of Thermal Decomposition Properties and Compatibility
The thermal behavior and kinetic analysis of the samples were performed while using differential scanning calorimetry (DSC). The effects of CuO and Al/CuO on the thermal decomposition properties of nitrocellulose (NC) were investigated using DSC (Q2000, TA Co., New Castle, DE, USA) under an N 2 atmosphere at a flow rate of 50 mL·min −1 . The test temperature range was room temperature to 300 • C, and the heating rate was 5.0, 10.0, 15.0, 20.0, 25.0, and 30.0 • C·min −1 . The mass of each sample was about 0.22 mg and the sample cell is an aluminum crucible. The CuO-NC and Al/CuO-NC composites were obtained by physical mixing in an agate mortar at room temperature. The mass ratio of CuO: NC for the CuO-NC composites, and Al/CuO: NC for the Al/CuO-NC composites was 1: 19.
The DSC thermal analysis method was used to assess the compatibility of different CuO and Al/CuO nanothermites with NC. The results come from the above-mentioned DSC test at the heating rate of 10.0 • C·min −1 . The evaluation standard for compatibility of ingredients can refer to [29,30].  Figure 1 shows the SEM images of CuO prepared at different reaction temperatures. When the reaction temperature is 80 • C, the obtained CuO is spindle shaped in the whole, as shown in Figure 1a. That is thick in the middle and thin in both ends. There are many stripe protrusions in the middle spindle shaped structure, and the ends of spindle shaped structure are curly. When the reaction temperature is 120 • C, the obtained CuO is bamboo leaf-shaped (Figure 1b). There are a few protrusions on the surface, and the leaves are thin. When the reaction temperature is 160 • C, the obtained CuO is in the boat form ( Figure 1c). Its size becomes smaller than the structure that is shown in Figure 1a,b. When the reaction temperature is 160 • C, the size of obtained CuO is the smallest, and there are raised fragments on the surface (Figure 1d). These experimental results show that the size of the obtained CuO gradually decreases with the temperature increases and its morphology also obviously changed. The DSC thermal analysis method was used to assess the compatibility of different CuO and Al/CuO nanothermites with NC. The results come from the above-mentioned DSC test at the heating rate of 10.0 °C·min −1 . The evaluation standard for compatibility of ingredients can refer to [29,30]. Figure 1 shows the SEM images of CuO prepared at different reaction temperatures. When the reaction temperature is 80 °C, the obtained CuO is spindle shaped in the whole, as shown in Figure  1a. That is thick in the middle and thin in both ends. There are many stripe protrusions in the middle spindle shaped structure, and the ends of spindle shaped structure are curly. When the reaction temperature is 120 °C, the obtained CuO is bamboo leaf-shaped (Figure 1b). There are a few protrusions on the surface, and the leaves are thin. When the reaction temperature is 160 °C, the obtained CuO is in the boat form (Figure 1c). Its size becomes smaller than the structure that is shown in Figure 1a and 1b. When the reaction temperature is 160 °C, the size of obtained CuO is the smallest, and there are raised fragments on the surface (Figure 1d). These experimental results show that the size of the obtained CuO gradually decreases with the temperature increases and its morphology also obviously changed.  Figure 2 shows the SEM images of CuO prepared at 120 °C with different reaction times. When the reaction time is 2 h, the structure of obtained CuO is similar to the "shrimp head" with multiple curled "shrimp beard" at the end, as shown in Figure 2a. When the reaction time is 4 h, most of the obtained CuO structures are bamboo leaf-shaped ( Figure 2b). There are a lot of protrusions on the surface, and the leaves are thin. When the reaction time is 8 h, the obtained CuO is also bamboo leafshaped ( Figure 2c). There are a few protrusions at the middle region of surface, and the two ends are flaky-shaped. As the reaction time extend to 12 h (Figure 2d) or 16 h (Figure 2e), the morphology of the obtained CuO is found not significantly different from the CuO that is shown in Figure 2c. However, the agglomeration degree of CuO is obviously increased. This shows that the morphology of the obtained CuO is changed obviously with the extension of the reaction time. As the reaction time exceeds 8 h, the morphology is changed very little and it has significant aggregate.  Figure 2 shows the SEM images of CuO prepared at 120 • C with different reaction times. When the reaction time is 2 h, the structure of obtained CuO is similar to the "shrimp head" with multiple curled "shrimp beard" at the end, as shown in Figure 2a. When the reaction time is 4 h, most of the obtained CuO structures are bamboo leaf-shaped ( Figure 2b). There are a lot of protrusions on the surface, and the leaves are thin. When the reaction time is 8 h, the obtained CuO is also bamboo leaf-shaped ( Figure 2c). There are a few protrusions at the middle region of surface, and the two ends are flaky-shaped. As the reaction time extend to 12 h (Figure 2d) or 16 h (Figure 2e), the morphology of the obtained CuO is found not significantly different from the CuO that is shown in Figure 2c. However, the agglomeration degree of CuO is obviously increased. This shows that the morphology of the obtained CuO is changed obviously with the extension of the reaction time. As the reaction time exceeds 8 h, the morphology is changed very little and it has significant aggregate. TEM and SEM analysis were performed to probe the elemental composition and morphology of the bamboo leaf-like CuO(b) obtained by reacting at 120 °C for 8 h and corresponding Al/CuO nanothermites, the results are exhibited in Figure 3. It can be seen from the TEM image of bamboo leaf-like CuO(b) (Figure 3a) that the maximum leaf of CuO(b) is about 1.6 μm in length, 300 nm in width, and 57 nm in thickness. The leaf surface is not smooth. There are many folds and a few leaves have protrusions in the middle. The protrusion is about 25 nm wide and about 90 nm long. Figure 3b is a high-resolution transmission electron microscope (HRTEM) image of the CuO and its corresponding fast Fourier transform (FFT) pattern. It is found that the lattice fringes with a lattice distance of 0.27 nm appeared at this position, corresponding to the (110) plane. Figure 3c is the selected area electron diffraction (SAED) image of the CuO sample. From Figure 3c, one can see that there are some clear scattered points, which indicate that the CuO leaves are single crystal structure and well-crystallized. Figure 3d is the SEM image of Al/CuO nanothermites. The spherical Al nanopowders and bamboo leaf-like CuO stick to each other, as shown in Figure 3d. Additionally, there are some agglomerate for the Al nanopowders because of their small size. TEM and SEM analysis were performed to probe the elemental composition and morphology of the bamboo leaf-like CuO(b) obtained by reacting at 120 • C for 8 h and corresponding Al/CuO nanothermites, the results are exhibited in Figure 3. It can be seen from the TEM image of bamboo leaf-like CuO(b) (Figure 3a) that the maximum leaf of CuO(b) is about 1.6 µm in length, 300 nm in width, and 57 nm in thickness. The leaf surface is not smooth. There are many folds and a few leaves have protrusions in the middle. The protrusion is about 25 nm wide and about 90 nm long. Figure 3b is a high-resolution transmission electron microscope (HRTEM) image of the CuO and its corresponding fast Fourier transform (FFT) pattern. It is found that the lattice fringes with a lattice distance of 0.27 nm appeared at this position, corresponding to the (110) plane. Figure 3c is the selected area electron diffraction (SAED) image of the CuO sample. From Figure 3c, one can see that there are some clear scattered points, which indicate that the CuO leaves are single crystal structure and well-crystallized. Figure 3d is the SEM image of Al/CuO nanothermites. The spherical Al nanopowders and bamboo leaf-like CuO stick to each other, as shown in Figure 3d. Additionally, there are some agglomerate for the Al nanopowders because of their small size.

Morphology and Structure
We present their SEM images and energy-dispersive X-ray spectroscopy (EDS) patterns to analyze the morphology and composition of the above flaky-shaped CuO(f) and its corresponding Al/CuO(f) nanothermites samples, as shown in Figure 4a-d. It can be seen from Figure 4a that the shape of as-prepared CuO is flaky-shaped, the surface is not smooth, and it is covered with depressions and protrusions. Figure 4c is the SEM image of Al/CuO(f) nanothermites sample obtained by ultrasonic dispersion method with Al nanopowders and flaky-shaped CuO(f). Some spherical aluminum nanoparticles with smaller particle sizes are coated on the surface of some CuO fragments, as shown in Figure 4c. The self-agglomeration phenomenon of the Al nanopowders is weak due to the ultrasonic effect. The EDS and XRD results of flaky-shaped CuO(f) (Figures 4b and 5d) show that the sample obtained by the hydrothermal method is the pure CuO phase. Figure 4d is the EDS pattern of the Al/CuO(f) nanothermites. In combination with its XRD characterization results (Figure 5c Figure 4g is a HRTEM image of flaky-shaped CuO(f) and its corresponding FFT pattern (insert). It is found that the lattice fringes with a lattice distance of 0.27 nm appeared at this position, corresponding to the (110) plane. Figure 4h is a SAED pattern of the flaky-shaped CuO(f) sample. The clear scattered points presented in Figure 4h indicate that the flaky-shaped CuO(f) is a single crystal structure with well-crystallized. We present their SEM images and energy-dispersive X-ray spectroscopy (EDS) patterns to analyze the morphology and composition of the above flaky-shaped CuO(f) and its corresponding Al/CuO(f) nanothermites samples, as shown in Figure 4a to 4d. It can be seen from Figure 4a that the shape of as-prepared CuO is flaky-shaped, the surface is not smooth, and it is covered with depressions and protrusions. Figure 4c is the SEM image of Al/CuO(f) nanothermites sample obtained by ultrasonic dispersion method with Al nanopowders and flaky-shaped CuO(f). Some spherical aluminum nanoparticles with smaller particle sizes are coated on the surface of some CuO fragments, as shown in Figure 4c. The self-agglomeration phenomenon of the Al nanopowders is weak due to the ultrasonic effect. The EDS and XRD results of flaky-shaped CuO(f) (Figure 4b and 5d) show that the sample obtained by the hydrothermal method is the pure CuO phase. Figure 4d is the EDS pattern of the Al/CuO(f) nanothermites. In combination with its XRD characterization results (Figure 5c), it can be obtained that the Al/CuO(f) sample is a mixture of Al and CuO. Because there is no other characteristic diffraction peaks, except the Al and CuO. Figure 4e is a TEM image for the flaky-shaped CuO(f). It can be seen that the fragments of flaky-shaped CuO(f) are irregular in shape and have wrinkles on the surface. The thickness is about 90 nm and the dispersion is good. From the enlarged view of the flaky-shaped CuO(f) (Figure 4f), we can see that the folds on the fragments are clear. The largest fragment in size is approximately 0.75 μm wide and about 2.5 μm long, and the protrusions on the surface are about 0.13 μm wide. Figure 4g is a HRTEM image of flaky-shaped CuO(f) and its corresponding FFT pattern (insert). It is found that the lattice fringes with a lattice distance of 0.27 nm appeared at this position, corresponding to the (110) plane. Figure 4h is a SAED pattern of the flaky-shaped CuO(f) sample. The clear scattered points presented in Figure 4h indicate that the flaky-shaped CuO(f) is a single crystal structure with well-crystallized. The XRD was used to analyze the phase microstructure, and the results are shown in Figure 5.  The XRD was used to analyze the phase microstructure, and the results are shown in Figure 5.

Effect of CuO and Al/CuO on Thermal Decomposition of NC
The composite materials of CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC were prepared by mixing NC with CuO or Al/CuO to analyze the effect of bamboo leaf-like CuO(b), flakyshaped CuO(f) and their corresponding Al/CuO(b), Al/CuO(f) nanothermites on the thermal decomposition properties of NC. Figure 6 shows their SEM images. The CuO(b)-NC composite has a long rod shape and a rough surface. Around the rod shape structures, there are some CuO(b) particles. As can be seen from the magnified SEM image of CuO(b)-NC (Figure 6b), the surface of the NC is rough, and the bamboo leaf-shaped CuO(b) is adherent to the surface. The Al/CuO(b)-NC composite showed to be rod-shaped, and the bulk Al/CuO(b) grains are dispersed on the region near the NC short fibers, as shown in Figure 6c. It can be seen from Figure 6d that the Al/CuO(b) is adherent to  long rod shape and a rough surface. Around the rod shape structures, there are some CuO(b) particles. As can be seen from the magnified SEM image of CuO(b)-NC (Figure 6b), the surface of the NC is rough, and the bamboo leaf-shaped CuO(b) is adherent to the surface. The Al/CuO(b)-NC composite showed to be rod-shaped, and the bulk Al/CuO(b) grains are dispersed on the region near the NC short fibers, as shown in Figure 6c. It can be seen from Figure 6d that the Al/CuO(b) is adherent to the surface of the NC fibers, and there is some agglomeration on the spherical Al nanopowders and CuO fragments. The shape of CuO(f)-NC composite is rod-shaped or block-shaped with different length and size, as shown in Figure 6e. The magnified image of CuO(f)-NC (Figure 6f) shows that the surface of the NC is rough and a large number of CuO fragments adherent on the surface. The SEM images of Al/CuO(f)-NC is similar to the Al/CuO(b)-NC, and the structure is a short rod shape and the surface is rough, as shown in Figure 6g. There are also some agglomerations on the spherical Al nanopowders and CuO fragments (see Figure 6h).

Effect of CuO and Al/CuO on Thermal Decomposition of NC
The composite materials of CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC were prepared by mixing NC with CuO or Al/CuO to analyze the effect of bamboo leaf-like CuO(b), flakyshaped CuO(f) and their corresponding Al/CuO(b), Al/CuO(f) nanothermites on the thermal decomposition properties of NC. Figure 6 shows their SEM images. The CuO(b)-NC composite has a long rod shape and a rough surface. Around the rod shape structures, there are some CuO(b) particles. As can be seen from the magnified SEM image of CuO(b)-NC (Figure 6b), the surface of the NC is rough, and the bamboo leaf-shaped CuO(b) is adherent to the surface. The Al/CuO(b)-NC composite showed to be rod-shaped, and the bulk Al/CuO(b) grains are dispersed on the region near the NC short fibers, as shown in Figure 6c. It can be seen from Figure 6d that the Al/CuO(b) is adherent to the surface of the NC fibers, and there is some agglomeration on the spherical Al nanopowders and CuO fragments. The shape of CuO(f)-NC composite is rod-shaped or block-shaped with different length and size, as shown in Figure 6e. The magnified image of CuO(f)-NC (Figure 6f) shows that the surface of the NC is rough and a large number of CuO fragments adherent on the surface. The SEM images of Al/CuO(f)-NC is similar to the Al/CuO(b)-NC, and the structure is a short rod shape and the surface is rough, as shown in Figure 6g. There are also some agglomerations on the spherical Al nanopowders and CuO fragments (see Figure 6h). °C, 209.2 °C, and 209.3 °C. When compared with the NC, the thermal decomposition peak temperatures of CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC are lower than that of NC at 0.5 °C, 0.1 °C, 0.5 °C, and 0.4 °C , respectively. These results indicate that there is no reaction at a low temperature between NC and other Al and CuO reactants, and the compatibility of the bamboo leaf-like CuO(b), flaky-shaped CuO(f), Al/CuO(b), Al/CuO(f) with NC is good. Therefore, the composites can be used as a component in the preparation of propellants and explosives. To obtain a better understanding of the effect of the thermite reaction between Al nanopowders and CuO to the thermal decomposition of NC, the thermal reaction characteristics of Al/CuO(b) and Al/CuO(f) were investigated by DSC. Figure 8 shows the DSC curves for Al/CuO(b) and Al/CuO(f) nanothermites at a heating rate of 10 °C·min −1 . One can see that the exothermic peak temperature of and Al/CuO(f) nanothermites are almost the same, but the releasing heat per unit mass of Al/CuO(b) (1153 J·g 1 ) is obviously higher than the Al/CuO(f) (681.5 J·g 1 ). The weak endothermic peaks were observed at about 654 °C from the DSC curves of Al/CuO(b) and Al/CuO(f) nanothermites. These are the melting peak of Al. The molten aluminum continues to react with CuO. The releasing heat per unit mass of the Al/CuO(b) after the melting of aluminum is approximately 101.0 J·g 1 , and the Al/CuO(f) is very little. It can also see that the temperature of main exothermic peak for Al/CuO(b) and Al/CuO(f) nanothermites is far lower than the melting temperature of aluminum. The main exothermic peak temperatures of the thermite reaction are also much higher than the thermal decomposition of CuO-NC and Al/CuO-NC, as shown in Figure 7. This can further explain there being no reaction at low temperature between NC and other Al or Al/CuO. To obtain a better understanding of the effect of the thermite reaction between Al nanopowders and CuO to the thermal decomposition of NC, the thermal reaction characteristics of Al/CuO(b) and Al/CuO(f) were investigated by DSC. Figure 8 shows the DSC curves for Al/CuO(b) and Al/CuO(f) nanothermites at a heating rate of 10 • C·min −1 . One can see that the exothermic peak temperature of and Al/CuO(f) nanothermites are almost the same, but the releasing heat per unit mass of Al/CuO(b) (1153 J·g 1 ) is obviously higher than the Al/CuO(f) (681.5 J·g 1 ). The weak endothermic peaks were observed at about 654 • C from the DSC curves of Al/CuO(b) and Al/CuO(f) nanothermites. These are the melting peak of Al. The molten aluminum continues to react with CuO. The releasing heat per unit mass of the Al/CuO(b) after the melting of aluminum is approximately 101.0 J·g 1 , and the Al/CuO(f) is very little. It can also see that the temperature of main exothermic peak for Al/CuO(b) and Al/CuO(f) nanothermites is far lower than the melting temperature of aluminum. The main exothermic peak temperatures of the thermite reaction are also much higher than the thermal decomposition of CuO-NC and Al/CuO-NC, as shown in Figure 7. This can further explain there being no reaction at low temperature between NC and other Al or Al/CuO.

Non-isothermal Kinetic Analysis
In order to explore the reaction mechanism of the intense exothermic decomposition process of NC, CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC, the thermal decomposition reaction kinetics was investigated by the non-isothermal DSC method. The kinetic parameters, apparent activation energy (Ea), pre-exponential factor (A), and kinetic model f( α) were obtained.

Non-isothermal Kinetic Analysis
In order to explore the reaction mechanism of the intense exothermic decomposition process of NC, CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC, the thermal decomposition reaction kinetics was investigated by the non-isothermal DSC method. The kinetic parameters, apparent activation energy (E a ), pre-exponential factor (A), and kinetic model f (α) were obtained. From the non-isothermal DSC curves at different heating rates 5.0, 10.0, 15.0, 20.0, 25.0, and 30.0 • C·min −1 , one can obtained the values of the extent of conversion (α) to corresponding temperature (T) by integrating the peak area of DSC curves at different heating rates. The E α can be obtained from the slope of the liner plot of lgβ i versus T i in the isoconversional Flynn-Wall-Ozawa's method (Equation (1)) [38]. Subsequently, the values of E α to α can be obtained by repeating the procedure for a set of different α.
where, β is heating rate (K·min −1 ); A is pre-exponential factor (s −1 ); α is the extent of conversion; E α is activation energy (J·mol −1 ) at α; R is universal gas constant (8.314 J·mol −1 ·K −1 ); G(α) is the integral form of the reaction model; and, T is temperature (K). Some of the reaction models used in the non-isothermal kinetic analysis are listed in the literature [39]. The isoconversional methods (i.e. Flynn-Wall-Ozawa's method) do not need to know the reaction model, and the more accurate activation energy can also be obtained by using a set of non-isothermal curves under different heating rates [40,41]. Figure 9 shows the dependences of E α to α for NC, different CuO-NC, and Al/CuO-NC. The isoconversional results of E α to α are important in detecting and treating the multistep kinetics. The E α -α curves of the decomposition process for the five samples have almost the same characteristic, as can be seen from Figure 9. The activation energies of NC, different CuO-NC, and Al/CuO-NC have little changes in the range of 0.100~0.750. Additionally, the ranges can be selected to calculate the non-isothermal reaction kinetics parameters.  7)) and one differential method (Kissinger (Equation 8)) were employed. By substituting the α-T data and forty-one types of kinetic model functions into the integral and differential equations, the apparent activation energy (Ea), pre-exponential factor (A), and the most probable kinetic model f( α) were obtained by the logical choice method with the best linear correlation coefficient (r) [27]. Generally, the normal range of Ea and lg(A/s −1 ) for energetic materials is approximately 80 to 250 kJ•mol −1 and 7 to 30, respectively. The values of Ea, lg(A/s −1 ), and corresponding linear correlation coefficient (r) that was obtained by different methods at different heating rates are listed in Table 1 and Table 2. The values of Ea and lg(A/s −1 ) obtained from each single In order to obtain the kinetic parameters, six integral methods (MacCallum-Tanner (Equation (2)), Šatava-Šesták (Equation (3)), Agrawal (Equation (4)), General integral (Equation (5)), Universal integral (Equation (6)), Ozawa (Equation (7)) and one differential method (Kissinger (Equation (8))) were employed. By substituting the α-T data and forty-one types of kinetic model functions into the integral and differential equations, the apparent activation energy (E a ), pre-exponential factor (A), and the most probable kinetic model f (α) were obtained by the logical choice method with the best linear correlation coefficient (r) [27]. Generally, the normal range of E a and lg(A/s −1 ) for energetic materials is approximately 80 to 250 kJ·mol −1 and 7 to 30, respectively. The values of E a , lg(A/s −1 ), and corresponding linear correlation coefficient (r) that was obtained by different methods at different heating rates are listed in Tables 1 and 2. The values of E a and lg(A/s −1 ) obtained from each single non-isothermal DSC curve is in good agreement with the calculated values that were obtained by Kissinger's method and Ozawa's method. Therefore, one can conclude that the reaction mechanism of the intense exothermic decomposition process of NC, CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC is classified as Avrami-Erofeev equation: f (α) = 3/2(1−α)[−ln(1−α)] 1/3 (differential form) and G(α) = [−ln(1−α)] 2/3 (integral form).
RT e (or p) ln where, β is heating rate (K·min .1 ); A is pre-exponential factor (s −1 ); E a is apparent activation energy (J·mol −1 ); R is universal gas constant (8.314 J·mol −1 ·K −1 ); G(α) is the integral form of the reaction model; α is the extent of conversion; T is temperature (K); T 0 is the initial temperature (K) at which DSC curve deviates from the baseline of the non-isothermal DSC curve; T e is the onset temperature (K); T p is the peak temperature (K); E eO (J·mol −1 ); and, E pO (J·mol −1 ) are the apparent activation energy obtained from T e and T p by Ozawa's method, respectively; E K is the apparent activation energy (J·mol −1 ) obtained from T p by Kissinger's method. Table 1. Calculated values of kinetic parameters of decomposition reaction for NC, bamboo leaf-like CuO(b)-NC, and Al/CuO(b)-NC.  Substituting the E a and A values (listed in Tables 1 and 2)

Thermal Safety Analysis
The simple model derived earlier in the frame of Semenov's thermal explosion theory [43] was established for estimating the thermal ignition temperature (T be0 ) and the critical temperatures of thermal explosion (T bp0 ) of as-prepared NC-based complexes. All of the calculations are based on the differential scanning calorimetry (DSC) experiments and the safety parameters, such as the thermal conductivity, particle size, mechanical properties, and pressure dependence of melting point et al., are not involved [44,45].
The values of T e0 andT p0 can be calculated in accordance with the Equation (15) in order to obtain the self-accelerating decomposition temperature (T SADT ) of as-prepared samples. Then the value of T SADT can be obtained by the Equation (16). The values of T e0 (T SADT ) and T p0 for NC, CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC are listed in Table 3.
T ei(or pi) = T e0(or p0) + aβ i +bβ i 2 +cβ i 3 (15) where T e is the onset temperature; T p is the peak temperature; β is the heating rate; T e0 andT p0 are the onset and peak temperatures corresponding to β→0, respectively; and, a, b, and c are the polynomial coefficients; i= 1, 2, ···, 6. T SADT =T e0 (16) Table 3. Calculated values of kinetic parameters of decomposition reaction for NC, CuO(b)-NC, Al/CuO(b)-NC, CuO(f)-NC, and Al/CuO(f)-NC. The T be0 can be obtained substituting the values of E eO (listed in Tables 1 and 2) andT e0 (see Table 3) into Equation (16). The T bp0 can also be obtained by substituting E pO (listed in Tables 1 and 2) andT p0 (see Table 3) into the Equation (17). The values of T be0 and T bp0 for NC, CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC are also listed in Table 3. Generally, the value of T b (T be0 or T bp0 ) is one of the most important evaluation parameters for thermal safety, which represent the degree of difficulty of the transition from thermal decomposition to thermal explosion. The higher the value of T b , the transition can take place easier.
T be0 or bp0 = E eO or pO − E 2 eO or pO − 4E eO or pO RT e0 or p0 2R (17) While using T = T p0 , E a = E K , and A = A K , the values of activation entropy (∆S ), activation enthalpy (∆H ), and activation of activation free energy (∆G ) of the main exothermic decomposition reaction of the NC, CuO(b)-NC, CuO(f)-NC, Al/CuO(b)-NC, and Al/CuO(f)-NC are obtained by Equation (18) to Equation (20), as listed in Table 3.
According to the previous studies [27,28,46], the thermal decomposition reaction of NC is a typical competition reaction between the O-NO 2 bond rupture and the decomposition of polymer skeleton products decomposition. The NO 2 gas is the initial decomposition product in the first step of high-temperature pyrolysis, similar to RDX. However, the NO 2 that is produced by RDX pyrolysis is derived from the rupture of the N-NO 2 bond [47]. However, the O-NO 2 bond rupture is deemed to be the first step, resulting in the release of NO 2 . The NO 2 stagnates in the polymer skeleton and then reacts with the RO• radical or its decomposition products to produce the NO, NO 2 , CO 2 , CO, H 2 O, N 2 O, HCHO, HCOOH, etc. Finally, the secondary autocatalytic reaction is significantly strengthened.
From the above calculation results, it could be derived that the addition of CuO and Al/CuO nanothermites can reduce the apparent activation energy (E a ), pre-exponential factor (A), onset temperature (T e0 ), thermal decomposition peak temperature (T p0 ), critical thermal ignition temperature (T be0 ), and the critical temperatures of thermal explosion (T bp0 ), as compared with single-component NC. This indicates that the CuO and Al/CuO nanothermites can accelerate the thermal decomposition of NC. When compared with the CuO-NC composite system, the E a , A, T e0 , T p0 , and T be0 values of the Al/CuO-NC composite system are all decreased. It indicates that the addition of Al nanopowders can increase the reactive sites and then accelerate the thermal decomposition of NC because of the high specific surface area of Al nanopowders or CuO. The physical adsorption might occur between the Al nanopowders or CuO and decomposition product of nitrocellulose. This may further effect the autocatalytic decomposition reaction of nitrocellulose.
From Tables 1 and 2, one can find that the E a values of the thermal decomposition processes of Al/CuO(b)-NC and Al/CuO(f)-NC are 181.1 kJ·mol −1 and 180.6 kJ·mol −1 , and the A values are 10 17.4 s −1 and 10 17.3 s −1 , respectively. The little difference of E a and A between the two nanothermites, indicating that the addition of Al/CuO(b) and Al/CuO(f) may reduce the thermal decomposition energy barrier of the composites and make it easy to decompose. However, the self-accelerated decomposition temperature (T SADT ) of Al/CuO(f)-NC is 152.9 • C, which is 16.4 • C lower than that of Al/CuO(b)-NC. Additionally, the thermal decomposition peak temperature (T p0 ) of Al/CuO(f)-NC is 7.4 • C higher than the Al/CuO(b)-NC. This shows that the decomposition reaction of Al/CuO(f)-NC is stronger than that of Al/CuO(b)-NC, but the time-consumption of the decomposition process of Al/CuO(f)-NC is long. From the results of the critical thermal ignition temperature (T be0 ) and the critical temperatures of thermal explosion (T bp0 ), one can get that the T be0 of Al/CuO(f)-NC is 17.3 • C lower than the Al/CuO(b)-NC, but T pe0 of Al/CuO(f)-NC is higher than the Al/CuO(b)-NC. This indicates that the Al/CuO(f)-NC is more easily ignited, but its decomposition rate is relatively slow. Accordingly, the Al/CuO(f)-NC has a much higher thermal stability than that of the Al/CuO(b)-NC.

Conclusions
Two kinds of nanosize CuO with different morphologies and sizes were prepared via the hydrothermal method by adjusting the reaction temperature and time. The ultrasonic composite method was used to combine the aluminum nanopowders with the as-prepared bamboo leaf-shaped CuO(b) and flaky-shaped CuO(f) to obtain the corresponding Al/CuO(b) and Al/CuO(f) nanothermites. Based on the non-isothermal decomposition kinetics, the catalytic effect of the different morphologies of CuO and corresponding nanothermites on the thermal decomposition properties of NC was investigated. The additions of CuO(b), Al/CuO(b), CuO(f), and Al/CuO(f) do not change the thermal decomposition mechanism of NC. The catalytic effect of Al/CuO nanothermites to NC is better than the CuO, and the Al/CuO(b)-NC, Al/CuO(f)-NC is easier to decompose than the CuO(b)-NC and CuO(f)-NC. The results of compatibility and thermal safety analysis show that the CuO and Al/CuO catalysts have good compatibility with NC, and the catalysts can be used safely. However, the agglomeration of Al/CuO(f) is relatively serious. The self-accelerated decomposition temperature (T SADT ) and the critical thermal ignition temperature (T be0 ) of Al/CuO(f)-NC are the lowest, while the thermal decomposition peak temperature (T p0 ) and the critical temperatures of thermal explosion (T bp0 ) are higher. It shows that the Al/CuO(f)-NC composite material is more easily ignited after being integrated with NC and Al/CuO(f), and the thermal stability during the thermal decomposition process is also better. The CuO(b), Al/CuO(b), CuO(f), and Al/CuO(f) as the catalysts has a wide application prospect in solid propellants.