Next Article in Journal
Inter-Month Nutrients Dynamic and Plant Growth in Calamagrostis angustifolia Community and Soil after Different Burning Seasons
Next Article in Special Issue
Investigation of Mineral Oil and CuO Mixed Synthetic Oil in Compression Ignition Engines: A Comparison of Physicochemical Attributes
Previous Article in Journal
Wildfire-Residential Risk Analysis Using Building Characteristics and Simulations to Enhance Structural Fire Resistance in Greece
Previous Article in Special Issue
Study on Pyrolysis Behaviors and Characteristics, Thermodynamics, Kinetics, and Volatiles of Single-Base Propellant, a Typical Energy-Containing Material
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetic Analysis of Thermal Decomposition of Polyvinyl Chloride at Various Oxygen Concentrations

1
Shaanxi Emergency Rescue and Aviation Forest Protection Center, Xi’an 710054, China
2
Forensic Science Institute, China People’s Police University, Langfang 065000, China
*
Author to whom correspondence should be addressed.
Fire 2023, 6(10), 404; https://doi.org/10.3390/fire6100404
Submission received: 23 September 2023 / Revised: 10 October 2023 / Accepted: 16 October 2023 / Published: 20 October 2023

Abstract

:
PVC plastic products are common combustible substances seen in fires, but their thermal degradation behavior under different oxygen concentrations has not been adequately studied. The thermal degradation behavior of PVC materials in atmospheres with different oxygen concentrations was analyzed via thermogravimetric–Fourier transform infrared spectroscopy (TG-FTIR). The TG results show that the thermal degradation process of PVC under a non-oxygenated atmosphere occurred in two stages, and the activation energies of the two stages were 130–175 KJ mol−1 and 230–320 KJ mol−1, respectively; under the oxygenated atmosphere, the thermal degradation process occurred in three stages. The activation energies of the three stages were 130–175 KJ mol−1, 145–510 KJ mol−1 and 75–190 KJ mol−1, respectively. And the reaction mechanism of the second stage of thermal degradation was changed from D-ZLT3 to En by the higher oxygen concentration. Infrared spectroscopy (FTIR) was used to analyze the pyrolysis process of PVC in the non-oxygenated atmosphere, and the eight major components were as follows, in descending order according to amount released: C-H stretching > HCl > C-Cl stretching > H2O > CO2 > C-H bending > C-H aliphatic bending > CH2. For the reaction of PVC at an oxygen concentration of 7%, the nine major components, in descending order according to amount released, were as follows: CO2 > HCl > H2O > CO > C-H stretching > C-Cl stretching > C-H aliphatic bending > C-H bending > CH2. For PVC reactions at oxygen concentrations of 14% and 21%, the five major components, in descending order according to amount released, were CO2 > HCl > CO > C-Cl stretching > H2O.

1. Introduction

Plastics play an important role in our daily lives and are used to manufacture various items, such as packaging, films, lids, bags and containers. PVC-based plastics are widely used in the manufacture of plastic products such as pipes, window frame materials, rigid films and cables and have excellent resistance to aging and high temperatures [1]. However, in a fire, PVC materials are exposed to high temperatures, and thermal radiation causes a rapid increase in the surface temperature of the material, contributing to pyrolysis and the release of pyrolysis gases. For example, in electrical fires, PVC is often used as an insulation material for wires, and the PVC material undergoes pyrolysis and combustion reactions, releasing large amounts of toxic gases [2]. In fact, PVC products (e.g., cables, pipes, floors, etc.) are often placed in confined spaces or small areas, resulting in a lack of sufficient oxygen on the surface of PVC products. This is likely to affect the pyrolysis and combustion reactions of PVC materials. Under a non-oxygenated atmosphere or low oxygen concentration, the combustion of plastics, wood, etc., will produce more toxic gases, such as CO, due to insufficient oxidation. For PVC materials, a typical polymer plastic, pyrolysis under different oxygen concentrations can produce different and complex pyrolysis products, some of which are toxic and harmful gases that can lead to serious injuries and economic losses and increase the difficulty of fire suppression. Therefore, it is crucial to study the variation in oxygen concentration in fires to understand the pyrolysis process of PVC materials and the effect of the released gases [3].
In recent years, scientists have studied the pyrolysis process and the pyrolysis products of PVC materials during fires. Thermogravimetric–infrared coupled analysis (TG-FTIR) can be used as the main method to study the pyrolysis process and pyrolysis products of substances under different oxygen concentrations. Benes et al. [4] used thermogravimetric analysis (TGA) combined with mass spectrometry and Fourier transform infrared (FTIR) spectroscopy to study this process. Thermal degradation of VC insulation using specimens was found to occur in three temperature ranges: 200–340 °C, 360–530 °C and 530–770 °C. And the PVC skeleton started to degrade in the range of 200~340 °C with the release of HCl, H2O, CO2 and benzene. Bhargava et al. [5] used the thermogravimetric analysis technique to carry out detailed kinetic analyses of PVC materials using a model-free approach and distributional reactivity models, and various distributional reactivity models were used to explain the thermal decomposition process of PVC polymers. Wang et al. [6] investigated the pyrolysis characteristics of PVC using TG-FTIR and found that the zero-order model and the sixth-order model were responsible for the two PVC sheath pyrolysis stages, respectively. Fourier transform infrared spectroscopy (FTIR) analysis showed that for PVC insulation, the release of seven major components were, in descending order: CO2 > C-H stretching > C-H bending > H2O > C-H aliphatic bending > CH2 deformation > C-Cl stretching. Wang et al. [7] used multiple thermogravimetric experiments in combination with FTIR analysis to determine the pyrolysis behavior of PVC. Possible reaction mechanisms were predicted using the generalized master plot method, and the relevant pre-exponential factors were calculated using the compensation effect. The analysis yielded changes in thermal degradation behavior attributable to changes in chemical composition, molecular structure, compositional ratios and various additives. In their study on pyrolysis and the combustion of substances under different oxygen concentrations in a fire, Fang et al. [8] investigated the pyrolysis and combustion process of wood under different oxygen concentrations using TG-FTIR analysis and analyzed in-depth the effect of oxygen concentration on the pyrolysis process and fugitive gases. Moreover, the apparent activation energy of pyrolysis and combustion varied linearly with oxygen concentration, and the effect of oxygen concentration on the combustion mechanism of wood was investigated. Currently, thermal analysis of PVC cable insulation materials is usually carried out for an air or nitrogen atmosphere. No extensive research has been carried out on how changes in oxygen concentration affect the thermal degradation properties of PVC materials.
Present-day studies on the pyrolysis of PVC using coupled TG-FTIR techniques at different oxygen concentrations are lacking, and the analyses are not yet complete. Based on previous studies, it is known that the thermogravimetric analysis method can effectively analyze the thermal weight loss process under continuous heating conditions [9,10,11]. Therefore, the heating rates of 10, 20, 30 and 40 min−1 were adopted to simulate the heating process of PVC materials in a real fire. The thermal degradation activation energy and reaction mechanism were analyzed using model-free and model-fitting methods, respectively, and the composition of the released gases was determined via characteristic spectroscopy. The aim of this study was to gain insight into the pyrolysis characteristics of polyvinyl chloride at different oxygen concentrations and determine how the pyrolysis and combustion processes are affected by changes in the oxygen concentration gradient. Semi-quantitative methods of infrared spectroscopy were also used to characterize the temperature dependence of the escaping gases. This provides more accurate data to support improved fire protection measures.

2. Materials and Methods

2.1. Materials and Sample Preparation

This study used commercially available PVC composition materials. The sample material is used in large quantities in the manufacture of commercial PVC products (such as plumbing, electrical wiring, etc.) and was provided by Fengtai Polymer Material Co., Ltd. (Weifang, Shandong Province, China), as a reference. The particle diameter was 0.06–0.08 cm. Through oven drying at 80 °C for 12 h, the sample species comprising mostly water were removed, ensuring that the mechanism analysis was not affected by atmospheric water molecules.

2.2. TG-FTIR Tests

Thermogravimetric analysis is considered to be the most effective and primary method for evaluating the behavior of heat pyrolysis of substances in a fire scenario and for calculating their kinetic parameters. For this purpose, a thermogravimetric analyzer, the STA4499F5 thermogravimetric analyzer model from NETZSCH-Geratebau GmbH, instrument Co., Ltd. (Selb, Germany), was used. The prepared powder samples of PVC insulation powder were placed in an alumina crucible and weighed using an inbuilt TG balance. The sample mass was controlled at 5 ± 0.2 mg, the gas flow rate was 100 mL/min, and the atmosphere was filled with nitrogen and oxygen set at 21%, 14%, 7% and 0% oxygen. The heating rates were set to 10, 20, 30 and 40 Kmin−1 and the temperature range was 35–1000 °C to simulate the pyrolysis process of PVC under different oxygen concentrations in a fire scenario. Each test was repeated at least three times to verify the reproducibility of the results. The maximum deviation was less than ±3%.
TG-FTIR provides additional information on the emission characteristics of gases. The gaseous substances produced through the programmed heating in the TGA entered the Bruker Invenio S spectrometer via a heat transfer tube that can be heated up to 230 °C to prevent heat dissipation condensation of the decomposition products, thus allowing simultaneous STA-FTIR analysis as well as the combination of quantitative TG analysis and FTIR. FTIR has a spectral range of 650–4000 cm−1 with a resolution of 4 and 32 cm−1 scans per sample.

2.3. Kinetics Analysis

To obtain specific data on the thermal degradation induced by the pyrolysis reaction within the thermogravimetric analysis, with regular increases in temperature with time, the pyrolysis parameters obtained from the thermogravimetric analysis can be evaluated using non-isothermal kinetic methods [12]. The conversion rate of PVC insulation during pyrolysis can be expressed as
α = m 0 m t m 0 m f
where m0 = initial mass, mt = instantaneous mass and mf = final mass. The rate of transformation of a substance can be expressed as
d α d t = β d α d T = K ( T ) f ( α )
where t, T, β, f(x) and K(T) represent the time (min), temperature (K), ambient warming rate = dT/ (Kmin−1), the differential transformation function of the reaction mechanism and the rate constant of the sample, respectively, as described by the Arrhenius Law in Equation (3).
K ( T ) = A exp ( E R T )
where E, A and R represent the activation energy required for the reaction (kJ mol−1), the pre-exponential factor (s−1) and the gas constant (8.314 J mol−1 K−1), respectively. Combining Equations (2) and (3), with α = 0 and T being equal to the starting temperature T0, the integral form g(α) of the mechanism function can be expressed as
g ( α ) = 0 α d α f ( α ) = A β T 0 T exp ( E R T ) d T

2.3.1. Model-Free Method

The model-free method is a standard approach to estimating the activation energy of composites without knowing or assuming a reaction model [13]. The two commonly used model-free methods chosen for this study were the Kissinger–Akahira–Sunose (KAS) [14] and Friedman (FR) [15], models as shown in Table 1.

2.3.2. Model-Fitting Method

The reaction mechanism can be determined using the Málek method [16] and is expressed as follows
Y α = Z α Z 0.5 = f α g α f 0.5 g 0.5 = T α T 0.5 2 d α / d t α d α / d t 0.5
where ( d α / d t ) 0.5 and T 0.5 correspond to the conversion rate and temperature at a conversion rate of 0.5, respectively. The mathematical equations on the left-hand side of the equation are general theoretical models that are constructed. These models can be used to characterize the rates of reactions, substance transformations, etc. As shown on the right-hand side, experimental data are obtained by heating the sample and measuring the change in the rate of reaction. In order to determine the most appropriate model to describe the pyrolysis of the sample, we need to compare the experimental data with the theoretical model. If the model can accurately predict the experimental results, then we can assume that the model is most suitable for describing the pyrolysis process of the sample. Shown in Table 2 [17] is the general theoretical model corresponding to the formula. It is a general theoretical model corresponding to the formula used to describe the mathematical expression of the relationship between the change in the thermal properties of a substance in the process of thermal analysis, such as enthalpy change, heat capacity change and other factors, such as temperature and time.

3. Results and Discussion

3.1. Thermogravimetric Analysis

In order to investigate the pyrolysis process of PVC materials under different oxygen concentrations in fire conditions, the thermogravimetric loss behavior of the whole process of combustion and pyrolysis of PVC materials under different oxygen concentrations was analyzed by determining the TG-DTG curves of experimental samples with four sets of different oxygen concentrations (0–21%) with the heating rates 10–40 °C min−1, as shown in Figure 1 and Figure 2. According to a previous study, the pyrolysis process of PVC usually has two main weight loss peaks [18]. The first weight loss peak is the autocatalytic reaction of the dehydrochlorination of PVC material, which is the main weight loss stage of PVC pyrolysis. The second weight loss peak may be due to the breakage of conjugated polyene chains and interactions within the molecule [19]. And, as mentioned in Figure 1 and Figure 2, due to the participation of different concentrations of oxygen in the reaction, the TG and DTG curve with a high oxygen concentration shifts towards the low-temperature section in the later stage of the reaction. The appearance temperature and peak value of the first weight loss peak of PVC pyrolysis did not change significantly with the increase in oxygen concentration. However, the second weight loss peak changed significantly. For the second weight loss peak, the values were 462, 439, 438 and 438 at 10 °C min−1 for 0%, 7%, 14% and 21% oxygen concentrations, respectively, and 476, 448, 449 and 447 at 20 °C min−1. The values were 487, 453, 455 and 456 at 30 °C min−1. The values were 492, 462, 461 and 459 at 40 °C min−1, respectively. It can be observed that the characteristic temperature points associated with the second weight loss peak in PVC pyrolysis shifted to higher temperatures as the heating rate increased. However, the characteristic temperature points of the reaction remained essentially unchanged regardless of variations in oxygen concentrations. It is worth noting that there was a third distinct peak for PVC pyrolysis with oxygen, which was seen at the values 567, 542 and 532 at 10 °C min−1 for 7%, 14% and 21% oxygen concentrations, 600, 575 and 552 at 20 °C min−1, 632, 587 and 572 at 30 °C min−1 and 656, 612 and 585 at 40 °C min−1, respectively, and the third peak of PVC pyrolysis was sharper at higher oxygen concentrations.
The characteristic values of the main DTG peaks of PVC pyrolysis at different oxygen concentrations are shown in Table 3. It can be observed that the weight loss process of the dehydrochlorination reaction in the pyrolysis of PVC materials was not affected by the presence of oxygen. This is because oxygen reacts with the dehydrochlorination products to form water and other decomposition products. In this reaction, PVC is dechlorinated to form hydrochloric acid and other products. These products then react with oxygen to form water and other oxidation products. This is used to compensate for the elemental hydrogen and chlorine consumed in the oxidation reaction. As a result, the influence on the TG curve is not significant. Further analysis of this phenomenon will be discussed later. It is worth noting that there was more pyrolysis residue of PVC in the non-oxygenated atmosphere than at other oxygen concentrations. This indicates that substances react more completely when oxygen is involved in the reaction. However, the pyrolysis rate in the non-oxygenated atmosphere was always greater than that at different oxygen concentrations until the end of the reaction. This is because the pyrolysis process relies heavily on intermolecular thermal movement and the breaking of chemical bonds. PVC will crack more readily under anaerobic conditions due to the low intermolecular forces between PVC molecules. However, as the oxygen concentration increases, oxygen begins to participate in the reaction process. This introduces additional chemical reaction channels and energy dissipation pathways that reduce the rate of the dehydrochlorination reaction [20]. The reaction rate is greater in the presence of oxygen, where the pyrolysis of PVC molecules undergoes an oxidation reaction. The hydrocarbons in the molecular chain are oxidized by oxygen to form CO2 and H2O. Therefore, as can be seen in the graph, even a small amount of oxygen intervening in the reaction will change PVC from a fast pyrolysis reaction dominated by the dehydrochlorination reaction to a slower oxidation reaction dominated by oxygen.

3.2. Activation Energy Analysis

Figure 3 shows the model-free activation energies, E, calculated using the KAS and Friedman methods with different conversion rates (Fitted correlation coefficient R2 values were calculated to be greater than 0.9). It was found that the E values calculated using the KAS and Friedman methods are similar, so they are not discussed separately. And it was clearly found that the activation energy value distribution of the splitting phenomenon, combined with the number of peaks in the previous section, were divided into two or three stages of pyrolysis [21]. In the figure, the splitting point is defined as the point at which there is an abrupt change in the pattern of activation energy change. The splitting point of PVC pyrolysis in the non-oxygenated atmosphere was α = 0.4, and the splitting points of PVC pyrolysis with oxygen participation were α = 0.5 and 0.8. The results show that the pyrolysis of PVC occurs in phases, regardless of whether there is oxygen participation. This is consistent with previous studies [22,23]. Generally speaking, the activation energies at different oxygen concentrations showed a relatively flat temperature trend in the first stage, with activation energy values in the range of 130–175 kJ mol−1, while in the second stage, the activation energies showed a wave growth trend due to the complexity and variability of the internal molecular structure reactions in the non-oxygenated atmosphere. In the second stage, due to the complexity of the internal molecular structure under the non-oxygenated atmosphere, the activation energies showed a wave growth trend, and the activation energies were in the range of 230–320 kJ mol−1 after a significant decrease at α = 0.7, while the activation energies required for oxidation reactions in pyrolysis reactions with different oxygen concentrations significantly increased and then rapidly decreased. The activation energies were in the range of 145–510 kJ mol−1, indicating that the oxidation reaction requires more energy to break the molecular bonds in the pyrolysis reaction with oxygen. Additionally, in the third stage, the breaking of this part of the bond energy led to the decomposition of the substances that could not be decomposed through pyrolysis. In this case, the hydrocarbon groups in the polymer were decomposed more fully, and the hard-to-decompose bonds, such as aromatic ring bonds, were ruptured to form hydrocarbons and other small groups. This shift in activation energy to the range of 75–190 kJ mol−1 signified a significant decrease in the activation energy. It is noteworthy that there was a sudden drop in E in the second stage with no oxygen. The reason for this is that the pyrolysis of PVC is specific due to its own material complexity. Different heating rates affect the reduction in activation energy. This observation is consistent with the TG and DTG curves shown in Figure 1 and Figure 2. It is noteworthy that the change in activation energy with the change in oxygen concentration was not significant. The reason for this may be that the activation energies of the reactants are usually related to the transition state structure and energy of the reaction. In the PVC pyrolysis reaction, with the change in oxygen concentration, the molecular structures of the reactants and products did not change significantly, and the structures of the transition states were relatively stable, resulting in an insignificant change in the activation energy.

3.3. Mechanism Function Analysis

As mentioned above, the pyrolysis process of PVC can be analyzed in two or three stages, and a single model-free approach cannot adequately represent the mechanism function of the PVC pyrolysis process. Therefore, for different stages of the pyrolysis process, the corresponding reaction mechanisms need to be determined [24]. Using the Málek method, experimental and theoretical calculations were plotted, as shown in Figure 4, which shows the experimental (point) and theoretical (line) main plots of the conversion α at different oxygen concentrations. In order to facilitate the comparison of the mechanism functions for different oxygen concentrations, the mechanisms of the second and third stages of the pyrolysis reaction involving oxygen were analyzed together with the second stage of the non-oxygenated atmosphere reaction, as Figure 4b. It is worth noting that the mechanism function of the first stage of PVC pyrolysis did not change significantly after oxygen participation in the reaction. However, the second stage of pyrolysis had a gradual change in the mechanism function due to the increasing concentration of oxygen.
From the figure, it can be inferred that for the pyrolysis reaction mechanism under the condition of no oxygen, the overall reaction mechanism was most similar to the exponent power models (E2). The reaction mechanisms under the reaction with oxygen participation were all close to the reaction order models F2 (0 ≤ α ≤ 0.5) and E2 (0.5 ≤ α ≤ 1). The changes in the reaction mechanism proved that oxygen participation in PVC pyrolysis changes the decomposition reaction from one dominated by dehydrochlorination in a non-oxygenated atmosphere to an oxidation reaction dominated by oxygen at the mechanistic level. In the second stage, the mechanism function was clearly affected by the oxygen concentration. The mechanism function was closest to that of the 3D Zhuravlev–Leskin–Tempelman (D-ZLT3) model in a non-oxygenated atmosphere and at a 7% oxygen concentration, to E2 at a 14% oxygen concentration and to F2 at a 21% oxygen concentration. These three theoretical models represent the diffusion model, reaction order model and exponent power model, implying that the oxygen concentration changes the dependence between the rate of warming and the reactants. The main reason for this difference is the increased demand for oxygen in the second stage, and a higher oxygen concentration will have more influence on the mechanism of PVC pyrolysis. The higher the concentration of oxygen, the greater the impact on the pyrolysis mechanism of PVC.
From the above discussion, we have gained insights into the thermal cracking behavior of PVC at varying concentrations of oxygen. A simulation of PVC pyrolysis was executed using NETZSCH Kinetics Neo, utilizing the heat loss data calculated earlier (as depicted in Figure 5). Notably, the R2 value of each model’s fit exceeded 0.9, underscoring the strong applicability of the fit. The reason why the R2 value at 7% and 14% oxygen was lower than that at 0% and 21% oxygen is because an insufficient oxidation reaction will be affected by the heating rate, which leads to an error in the calculated factor. This, in turn, influences the rate of temperature rise, causing inaccuracies in the computed value of A. In conclusion, NETZSCH Kinetics Neo simulation of the pyrolysis process of PVC under various oxygen concentrations demonstrated its capacity to accurately construct applicable models, thus affirming the reliability of the experiment.
The differences in the pyrolysis process of PVC at different oxygen concentrations were analyzed in detail above. The gases formed during the pyrolysis process can be analyzed via infrared spectroscopy, which provides a scientific analytical method with which to study the quality, nature and quantity of functional groups under different oxygen concentrations [25]. In Figure 6, the three-dimensional TG-FTIR spectrum of the pyrolysis products of PVC under different oxygen concentrations at a heating rate of 40 Kmin−1 are shown. Among them, the 3D TG-FTIR spectrum in the non-oxygenated atmosphere was significantly different from the spectrum seen for the pyrolysis reaction involving oxygen. The peak temperatures in the anoxic atmosphere were concentrated at 400–600 °C, while the peak temperatures of the pyrolysis reaction with oxygen were scattered at 350–800 °C. This indicates that the production of fugitive gases from the PVC pyrolysis reaction involving oxygen was widely distributed in different temperature bands, and the peak of fugitive gases for this reaction was significantly higher than that of fugitive gases in the non-oxygenated atmosphere [26]. However, the peak height did not represent the amount of fugitive gas produced, and the vibration frequencies of different peaks were not the same. The process of fugitive gas production under different oxygen concentrations will be analyzed in detail in the following.
Figure 7 shows the FTIR absorption spectra of the PVC material at a heating rate of 40 Kmin−1. Table 4 lists the main vibrational frequency bands of the gases from the PVC material. It can be seen that the main components of the gases were similar for different oxygen concentrations [27]. For the PVC material, there were obvious absorption bands at 4000~3500 cm−1 and 1800~1300 cm−1, which correspond to O-H stretching vibration and O-H bending vibration, respectively, showing the decomposition of oxygen-containing groups to form H2O. In addition, the characteristic peak at 2934 cm−1 indicates the presence of a stretching mode. The characteristic peak at 2934 cm−1 indicates the presence of C-H bonds in the stretching mode, which was due to the back-and-forth movement of carbon and hydrogen atoms in the group along the bonding direction. CO2 was located at 2400–2260 cm−1, and the peaks of CO2 were very sharp in the reaction involving oxygen. In addition, CO was located at 2173 and 2120 cm−1, representing the stretching vibration and bending vibration of the CO molecule. The peak signals for C-H aliphatic bending were weaker near 1458 cm−1. As a result of the volume, bond length and bond energy of C-C and C-H bonds in hydrocarbons, they interacted and caused bending of carbon atoms around neighboring carbon–hydrogen bonds. In addition, the characteristic absorption band of CH2 deformation was observed near 1338 cm−1. C-H bending near Cl atoms formed an absorption band at about 1240 cm−1, which refers to the carbon and hydrogen atoms in the group moving back and forth in a direction perpendicular to the bond. The signals found at 852 cm−1 and 660 cm−1 can be attributed to C-Cl stretching, which occurred due to a change in the relative plane angle of the C-Cl bond. Combined with the studies summarized above, the main pyrolytic components of the measured PVC insulation in a helium atmosphere were found to be HCl, H2O, CO2, CO, C-H stretching, C-H aliphatic bending, CH2 deformation, C-H bending and C-Cl stretching.
Figure 8 shows the regularization curves of the temperature evolution of the main pyrolysis components of PVC at different oxygen concentrations at a temperature increase rate of 40 Kmin−1. It was found that the gas of PVC with a temperature increase in the non-oxygenated atmosphere presents two peaks, while the reaction with oxygen participation presents more than two peaks. In PVC pyrolysis, this is the most intuitive phenomenon seen in the oxygen participation reaction regarding the composition of the gas.
For the PVC pyrolysis reaction in the non-oxygenated atmosphere, the first peak, which was mainly a sharp peak formed by the volatilization of HCl gas and the stretching vibration of the C-Cl bond, was in the temperature range of 380–420 °C. The second peak position was in the temperature range of 440–560 °C, and the changes in C-H stretching and vibration, volatilization of the remaining HCl gas, aliphatic bending of the C-H bond, CH2 deformation and C-H bending had a similar trend. For the PVC pyrolysis reaction involving oxygen, the changes brought about by the change in oxygen concentration were obvious. In particular, the change in the CO2 curves at a 7% oxygen concentration showed three peaks and a temperature range of 190–750 °C compared to those at an oxygen concentration above 14%. This is consistent with the weight loss temperature range of the DTG curve. The extra peaks were due to the increased demand for oxygen by the reactants, as evidenced by the increased CO production in the temperature range between the peaks. The most significant changes to occur due to the increasing oxygen concentration were seen in CO2, CO, HCl, C-H stretching and C-Cl stretching, with CO and HCl stretching being the main components of the toxic gases that contribute to the pyrolysis smoke produced by the burning of PVC in fires [28]. Overall, in conjunction with the TG-FTIR analysis process described above, for the pyrolysis process of PVC in the non-oxygenated atmosphere, the seven major components, in descending order according to amount released, were C-H stretching > HCl > C-Cl stretching > H2O > CO2 > C-H bending > C-H aliphatic bending. The reason for this is that the non-oxygenated atmosphere, as the dominant pyrolysis reaction, leads to weaker C-H bonds, which are easily broken due to high temperatures. Thus, C-H stretching increases the vibration frequency due to the high temperature, releasing a large number of hydrocarbon radicals, which triggers a chain reaction that promotes the decomposition of PVC and the release of other components. For the reaction of PVC at a 7% oxygen concentration, the nine main components, in descending order according to amount released, were CO2 > HCl > H2O > CO > C-H stretching > C-Cl stretching > C-H aliphatic bending > C-H bending > CH2 deformation. For PVC, the five main components in the reactions at 14% and 21% oxygen concentrations were, in descending order according to amount released, CO2 > HCl > CO > C-Cl stretching > H2O. The rest of the components did not have significant peaks in the infrared spectrograms due to the difference in oxygen concentration.

4. Conclusions

PVC materials at different oxygen concentrations were analyzed using the thermogravimetric–infrared coupling (TG-FTIR) method in the range of 10–40 Kmin−1. The thermogravimetric processes and gas components of PVC materials at different oxygen concentrations were obtained from the analysis. The main conclusions are as follows:
(1)
Different oxygen concentrations and heating rates in a fire affected the pyrolysis process of PVC materials. The increase in heating rate caused the TG curve to shift to the high-temperature section and led to a higher DTG peak. Oxygen concentration mainly affected the second weight loss peak of the PVC pyrolysis process. When oxygen was involved in the pyrolysis reaction, the second weight loss peak changed to two weight loss peaks. The prolongation of pyrolysis and the increase in oxygen concentration led to the gradual acceleration of the oxidation reaction of PVC.
(2)
The activation energies of PVC materials in oxygenated and non-oxygenated atmospheres increased with an increase in the conversion rate in the first stage, and the activation energies were in the range of 130–175 kJ mol−1. The activation energies in the oxygenated atmosphere showed a wave-type trend in the first stage, and the activation energies were in the range of 230–320 kJ mol−1; the activation energies in the second stage increased significantly, and then decreased rapidly, in the range of 130–510 kJ mol−1. In the third stage, the activation energies rose sharply before falling quickly, ranging from 75 to 510 kJ mol−1. The activation energies of PVC in the oxygen-containing atmosphere did not change significantly as a result of the increase in oxygen concentration, and the overall trends of the activation energies were similar.
(3)
The effect of oxygen concentration on the mechanism function of PVC pyrolysis was mainly a change in the reaction mechanism of the second pyrolysis stage from D-ZLT3 to E at higher concentrations of oxygen. But there was no significant effect on the reaction mechanism in the first pyrolysis stage.
(4)
The pyrolysis of PVC in a non-oxygenated atmosphere was investigated using infrared spectroscopy (FTIR), and the eight major released components, in descending order according to amount released, were C-H stretching > HCl > C-Cl stretching > H2O > CO2 > C-H bending > C-H aliphatic bending > CH2. For PVC at a 7% oxygen concentration, the nine major released components, in descending order according to amount released, were CO2 > HCl > H2O > CO > C-H stretching > C-Cl stretching > C-H aliphatic bending > C-H bending > CH2. For the PVC reaction at 14% and 21% oxygen concentrations, the five major released components, in descending order according to amount released, were: CO2 > HCl > CO > C-Cl stretching > H2O.

Author Contributions

Conceptualization, Y.W., writing—original draft preparation, S.Y., writing—review and editing, P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Data Availability Statement

Data is available in a publicly accessible repository.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, J.; Sun, L.; Ma, C.; Qiao, Y.; Yao, H. Thermal Degradation of Pvc: A Review. Waste Manag. 2016, 48, 300–314. [Google Scholar] [CrossRef]
  2. Babrauskas, V. Mechanisms and Modes for Ignition of Low-Voltage, Pvc-Insulated Electrotechnical Products. Fire Mater. 2006, 30, 151–174. [Google Scholar] [CrossRef]
  3. Wang, Z.; Wei, R.; Wang, X.; He, J.; Wang, J. Pyrolysis and Combustion of Polyvinyl Chloride (Pvc) Sheath for New and Aged Cables Via Thermogravimetric Analysis-Fourier Transform Infrared (Tg-Ftir) and Calorimeter. Materials 2018, 11, 1997. [Google Scholar] [CrossRef] [PubMed]
  4. Beneš, M.; Milanov, N.; Matuschek, G.; Kettrup, A.; Plaček, V. Thermal Degradation of Pvc Cable Insulation Studied by Simultaneous TG-FTIR and TG-EGA methods. J. Therm. Anal. Calorim. 2004, 78, 621–630. [Google Scholar] [CrossRef]
  5. Bhargava, A.; van Hees, P.; Andersson, B. Pyrolysis Modeling of Pvc and Pmma Using a Distributed Reactivity Model. Polym. Degrad. Stab. 2016, 129, 199–211. [Google Scholar] [CrossRef]
  6. Wang, C.; Liu, H.; Zhang, J.; Yang, S.; Zhang, Z.; Zhao, W. Thermal Degradation of Flame-Retarded High-Voltage Cable Sheath and Insulation Via Tg-Ftir. J. Anal. Appl. Pyrolysis 2018, 134, 167–175. [Google Scholar] [CrossRef]
  7. Wang, Z.; Xie, T.; Ning, X.; Liu, Y.; Wang, J. Thermal Degradation Kinetics Study of Polyvinyl Chloride (Pvc) Sheath for New and Aged Cables. Waste Manag. 2019, 99, 146–153. [Google Scholar] [CrossRef] [PubMed]
  8. Fang, M.X.; Shen, D.K.; Li, Y.X.; Yu, C.J.; Luo, Z.Y.; Cen, K.F. Kinetic Study on Pyrolysis and Combustion of Wood under Different Oxygen Concentrations by Using Tg-Ftir Analysis. J. Anal. Appl. Pyrolysis 2006, 77, 22–27. [Google Scholar] [CrossRef]
  9. Ma, W.; Rajput, G.; Pan, M.; Lin, F.; Zhong, L.; Chen, G. Pyrolysis of Typical Msw Components by Py-Gc/Ms and Tg-Ftir. Fuel 2019, 251, 693–708. [Google Scholar] [CrossRef]
  10. Sun, C.; Li, C.; Tan, H.; Zhang, Y. Synergistic Effects of Wood Fiber and Polylactic Acid During Co-Pyrolysis Using Tg-Ftir-Ms and Py-Gc/Ms. Energy Convers. Manag. 2019, 202, 112212. [Google Scholar] [CrossRef]
  11. Volli, V.; Gollakota, A.R.K.; Shu, C.M. Comparative Studies on Thermochemical Behavior and Kinetics of Lignocellulosic Biomass Residues Using Tg-Ftir and Py-Gc/Ms. Sci. Total Environ. 2021, 792, 148392. [Google Scholar] [CrossRef] [PubMed]
  12. Vyazovkin, S.; Burnham, A.K.; Criado, J.M.; Pérez-Maqueda, L.A.; Popescu, C.; Sbirrazzuoli, N. Ictac Kinetics Committee Recommendations for Performing Kinetic Computations on Thermal Analysis Data. Thermochim. Acta 2011, 520, 1–19. [Google Scholar] [CrossRef]
  13. Vyazovkin, S.; Wight, C.A. Model-free and model-fitting approaches to kinetic analysis of isothermal and nonisothermal data. Thermochim. Acta 1999, 341, 53–68. [Google Scholar] [CrossRef]
  14. Kissinger, H.E. Reaction kinetics in differential thermal analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  15. Friedman, H.L. Kinetics of Thermal Degradation of Char-Forming Plastics from Thermogravimetry. Application to a Phenolic Plastic. J. Polym. Sci. Part C Polym. Symp. 1954, 6, 183–195. [Google Scholar] [CrossRef]
  16. Málek, J. The kinetic analysis of non-isothermal data. Thermochim. Acta 1992, 200, 257–269. [Google Scholar] [CrossRef]
  17. Bai, F.; Guo, W.; Lü, X.; Liu, Y.; Guo, M.; Li, Q.; Sun, Y. Kinetic Study on the Pyrolysis Behavior of Huadian Oil Shale Via Non-Isothermal Thermogravimetric Data. Fuel 2015, 146, 111–118. [Google Scholar] [CrossRef]
  18. Wu, J.; Chen, T.; Luo, X.; Han, D.; Wang, Z.; Wu, J. Tg/Ftir Analysis on Co-Pyrolysis Behavior of Pe, Pvc and Ps. Waste Manag. 2014, 34, 676–682. [Google Scholar] [CrossRef] [PubMed]
  19. Zhu, H.M.; Jiang, X.G.; Yan, J.H.; Chi, Y.; Cen, K.F. Tg-Ftir Analysis of Pvc Thermal Degradation and Hcl Removal. J. Anal. Appl. Pyrolysis 2008, 82, 1–9. [Google Scholar] [CrossRef]
  20. Aboulkas, A.; El harfi, K.; El Bouadili, A. Thermal degradation behaviors of polyethylene and polypropylene. Part I: Pyrolysis kinetics and mechanisms. Energy Convers. Manag. 2010, 51, 1363–1369. [Google Scholar] [CrossRef]
  21. Vyazovkin, S.; Burnham, A.K.; Favergeon, L.; Koga, N.; Moukhina, E.; Pérez-Maqueda, L.A.; Sbirrazzuoli, N. Ictac Kinetics Committee Recommendations for Analysis of Multi-Step Kinetics. Thermochim. Acta 2020, 689, 178597. [Google Scholar] [CrossRef]
  22. Çepelioğullar, Ö.; Pütün, A.E. Thermal and Kinetic Behaviors of Biomass and Plastic Wastes in Co-Pyrolysis. Energy Convers. Manag. 2013, 75, 263–270. [Google Scholar] [CrossRef]
  23. Sánchez-Jiménez, P.E.; Perejón, A.; Criado, J.M.; Diánez, M.J.; Pérez-Maqueda, L.A. Kinetic model for thermal dehydrochlorination of poly (vinyl chloride). Polymer 2010, 51, 3998–4007. [Google Scholar] [CrossRef]
  24. Xu, L.; Jiang, Y.; Wang, L. Thermal decomposition of rape straw: Pyrolysis modeling and kinetic study via particle swarm optimization. Energy Convers. Manag. 2017, 146, 124–133. [Google Scholar] [CrossRef]
  25. Gao, N.; Li, A.; Quan, C.; Du, L.; Duan, Y. Tg–Ftir and Py–Gc/Ms Analysis on Pyrolysis and Combustion of Pine Sawdust. J. Anal. Appl. Pyrolysis 2013, 100, 26–32. [Google Scholar] [CrossRef]
  26. Wang, Y.; Kang, N.; Lin, J.; Lu, S.; Liew, K.M. On the Pyrolysis Characteristic Parameters of Four Flame-Retardant Classes of Pvc Sheathless Cable Insulation Materials. J. Anal. Appl. Pyrolysis 2023, 170, 105901. [Google Scholar] [CrossRef]
  27. Xie, Q.; Zhang, H.; Tong, L. Experimental study on the fire protection properties of PVC sheath for old and new cables. J. Hazard. Mater. 2010, 179, 373–381. [Google Scholar] [CrossRef] [PubMed]
  28. Singh, R.; Pant, D. Polyvinyl Chloride Degradation by Hybrid (Chemical and Biological) Modification. Polym. Degrad. Stab. 2016, 123, 80–87. [Google Scholar] [CrossRef]
Figure 1. Total mass loss (TG) versus temperature curves for PVC with different oxygen concentrations: (a) β = 10 Kmin−1, (b) β = 20 Kmin−1, (c) β = 30 Kmin−1, (d) β = 40 Kmin−1.
Figure 1. Total mass loss (TG) versus temperature curves for PVC with different oxygen concentrations: (a) β = 10 Kmin−1, (b) β = 20 Kmin−1, (c) β = 30 Kmin−1, (d) β = 40 Kmin−1.
Fire 06 00404 g001
Figure 2. Derivative total mass loss (DTG) versus temperature curves for PVC with different oxygen concentrations: (a) β = 10 Kmin−1, (b) β = 20 Kmin−1, (c) β = 30 Kmin−1, (d) β = 40 Kmin−1.
Figure 2. Derivative total mass loss (DTG) versus temperature curves for PVC with different oxygen concentrations: (a) β = 10 Kmin−1, (b) β = 20 Kmin−1, (c) β = 30 Kmin−1, (d) β = 40 Kmin−1.
Fire 06 00404 g002
Figure 3. Activation energy of thermal degradation of PVC versus conversion α as determined using the (a) KAS method and (b) Friedman method.
Figure 3. Activation energy of thermal degradation of PVC versus conversion α as determined using the (a) KAS method and (b) Friedman method.
Fire 06 00404 g003
Figure 4. Experimental and theoretical master plot results:(a) stage I, (b) stage II–III.
Figure 4. Experimental and theoretical master plot results:(a) stage I, (b) stage II–III.
Fire 06 00404 g004
Figure 5. Fitted and experimental curves of conversion rate versus temperature at different oxygen concentrations.(a) O2-0%, (b) O2-7%, (c) O2-14%, (d) O2-21%.
Figure 5. Fitted and experimental curves of conversion rate versus temperature at different oxygen concentrations.(a) O2-0%, (b) O2-7%, (c) O2-14%, (d) O2-21%.
Fire 06 00404 g005
Figure 6. Three-dimensional mapping of the FTIR spectrum of the gas at a heating rate of 40 Kmin−1: (a) O2-0%, (b) O2-7%, (c) O2-14%, (d) O2-21%.
Figure 6. Three-dimensional mapping of the FTIR spectrum of the gas at a heating rate of 40 Kmin−1: (a) O2-0%, (b) O2-7%, (c) O2-14%, (d) O2-21%.
Fire 06 00404 g006aFire 06 00404 g006b
Figure 7. Absorption spectra at the heating rate of 40 Kmin−1.
Figure 7. Absorption spectra at the heating rate of 40 Kmin−1.
Fire 06 00404 g007
Figure 8. Temperature evolution patterns of the main pyrolysis components of PVC: (a) O2-0%, (b) O2-7%, (c) O2-14%, (d) O2-21%.
Figure 8. Temperature evolution patterns of the main pyrolysis components of PVC: (a) O2-0%, (b) O2-7%, (c) O2-14%, (d) O2-21%.
Fire 06 00404 g008aFire 06 00404 g008b
Table 1. Common model-free functional equations.
Table 1. Common model-free functional equations.
Model-Free MethodExpressions
Kissinger–Akahira–Sunose (KAS) ln β T 2 = ln A R E g α E R T
Friedman (FR) ln d α d t = ln β d α d T = ln f α A E R T
Table 2. Common kinetic mechanism (mode) functions [17].
Table 2. Common kinetic mechanism (mode) functions [17].
NumberModelDifferential Form f ( α ) Integral Form G ( α )
Diffusion model
11D diffusion D1 1 / 2 α 1 α2
22D diffusion–Valensi D-V2 ln 1 α 1 α + 1 α ln 1 α
33D diffusion–Jande r D-J3 6 1 α 2 / 3 [ 1 1 α 1 / 3 ] 1 / 2 [ 1 1 α 1 / 3 ] 1 / 2
43D Zhuravlev–Leskin–Tempelman D-ZLT3 3 / 2 1 α 4 / 3 1 α 1 / 3 1 1 1 α 1 / 3 1 2
Sigmoidal rate equations
5Avarami–Erofeev A2 1 / 2 1 α ln 1 α 1 ln 1 α 2
Reaction order models
6Second-order chemical reaction F2 1 α 2 1 α 1 1
Exponent power models
7Second-order E2 1 / 2 α ln α 2
Geometrical contraction models
8Contracting area R2 2 1 α 1 / 2 1 1 α 1 / 2
Table 3. DTG parameters for PVC material under different oxygen concentrations.
Table 3. DTG parameters for PVC material under different oxygen concentrations.
AtmosphereHeat Rate
/°C min−1
Tset/°CTfinal/°CDTGpeak1
/% min−1
DTGpeak2
/% min−1
DTGpeak3
/% min−1
O2-0%1027550620.33.9/
2028652335.77.5/
3028953553.010.9/
4029254468.714.2/
O2-7%1027071519.23.72.3
2028279536.07.04.3
3028581452.18.34.4
4028582663.812.15.6
O2-14%1026864218.42.82.9
2027966137.55.84.8
3028368952.18.56.3
4028969266.110.38.8
O2-21%1027061020.02.93.27
2028063635.05.06.6
3028566351.47.57.9
4028574070.511.78.4
Table 4. Band assignments and corresponding wave numbers for PVC.
Table 4. Band assignments and corresponding wave numbers for PVC.
Description of VibrationsWave Numbers (cm−1)
HCl3150–2500
H2O4000–3500, 1800–1300
C-H stretching2934
CO22400–2260
CO2173, 2120
aliphatic bending of C-H1458
deformation of CH21338
C-H bending1240
stretching of C-Cl660, 852
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, S.; Wang, Y.; Man, P. Kinetic Analysis of Thermal Decomposition of Polyvinyl Chloride at Various Oxygen Concentrations. Fire 2023, 6, 404. https://doi.org/10.3390/fire6100404

AMA Style

Yang S, Wang Y, Man P. Kinetic Analysis of Thermal Decomposition of Polyvinyl Chloride at Various Oxygen Concentrations. Fire. 2023; 6(10):404. https://doi.org/10.3390/fire6100404

Chicago/Turabian Style

Yang, Shuo, Yong Wang, and Pengrui Man. 2023. "Kinetic Analysis of Thermal Decomposition of Polyvinyl Chloride at Various Oxygen Concentrations" Fire 6, no. 10: 404. https://doi.org/10.3390/fire6100404

APA Style

Yang, S., Wang, Y., & Man, P. (2023). Kinetic Analysis of Thermal Decomposition of Polyvinyl Chloride at Various Oxygen Concentrations. Fire, 6(10), 404. https://doi.org/10.3390/fire6100404

Article Metrics

Back to TopTop