Next Article in Journal
Advancing Green Hydrogen Purity with Iron-Based Self-Cleaning Oxygen Carriers in Chemical Looping Hydrogen
Previous Article in Journal
Prediction of Oxygen Evolution Activity for FeCoMn Oxide Catalysts via Machine Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyzing HDPE Thermal and Catalytic Degradation in Hydrogen Atmosphere: A Model-Free Approach to the Activation Energy

1
Centro de Química Estrutural, Institute of Molecular Sciences, Departamento de Engenharia Química, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
2
Departamento de Matemática Aplicada, Ciencia e Ingeniería de los Materiales y Tecnología Electrónica, Universidad Rey Juan Carlos, 28933 Madrid, Spain
3
Departamento de Engenharia Química, Instituto Superior de Engenharia de Lisboa, Instituto Politécnico de Lisboa, 1959-007 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Catalysts 2024, 14(8), 514; https://doi.org/10.3390/catal14080514
Submission received: 21 June 2024 / Revised: 2 August 2024 / Accepted: 3 August 2024 / Published: 9 August 2024
(This article belongs to the Section Industrial Catalysis)

Abstract

:
Despite the great interest in thermochemical processes for converting plastic waste into chemical feedstocks or fuels, their kinetics are still a less studied topic, especially under reductive conditions. In the present work, non-isothermal thermogravimetric analysis is used to study the thermal and catalytic conversion of HDPE promoted by parent and metal-based H-USY (15) and H-ZSM-5 (11.5) zeolites under a reducing hydrogen atmosphere. Additionally, the respective kinetic parameters (apparent activation energy, Ea, and frequency factor, A) were determined by applying two distinct model-free methods: Flynn–Wall–Ozawa (FWO) and Kissinger–Akahira–Sunose (KAS). The results showed that Ea of the thermal degradation of HDPE has an average value of 227 kJ/mol for both methods. In the presence of H-USY (15) and H-ZSM-5 (11.5) zeolites, Ea is strongly reduced and is highly dependent on conversion. In the case of H-USY (15), Ea varies from 78 to 157 kJ/mol for the KAS method and from 83 to 172 kJ/mol for the FWO method. Slightly lower values are reported for H-ZSM-5, with Ea values in the range of 53–122 kJ/mol for KAS and 61–107 kJ/mol for FWO. The presence and type of the metal source (Ni, Pt, or Pd) also affect the kinetic parameters of the reaction. The mean Ea values follow the order: Ni > Pt ≈ Pd for H-USY (15) or H-ZSM-5 zeolites. Accordingly, both parent and metal-based H-USY (15) and H-ZSM-5 zeolites can significantly reduce energy consumption in HDPE hydrocracking, thus promoting a more sustainable conversion of plastic waste.

1. Introduction

Over the past 50 years, plastics have provided and improved our lifestyle. However, changes in consumption and production patterns are leading to a rapid increase in plastic waste (PW) generation, which is becoming one of the most pressing environmental concerns [1]. The conversion of PW into resources is crucial for increasing efficiency and fostering circularity [2]. In recent years, chemical recycling and in particular, thermochemical processes, such as thermal and catalytic cracking, have emerged as promising options for converting PW into chemical feedstocks or fuels [3]. Thermal cracking (TC) is a widely used process but it suffers from two main problems: limited conversion of the feedstock at low temperatures and a wide carbon molecular weight distribution, leading to low selectivity for high-value products [4,5]. Catalytic cracking (CC), conversely, requires a catalyst [6] and allows the process to be carried out at lower temperatures, improving the yield and selectivity of the desired products [7]. Several authors have studied the TC and CC of polymers under an inert atmosphere, in particular polyethylene (PE) and polypropylene (PP), over different acidic catalysts [8], obtaining n-paraffins and olefins as the main products. However, to develop more efficient processes for the thermal and catalytic conversion of plastics, knowledge of the reaction kinetics is extremely important [9,10]. The kinetic parameters, apparent activation energy (Ea) and frequency factor (A), are closely related to the reaction and are highly dependent on the feedstock used, the reaction mechanism, the degree of conversion, and the kinetic study method chosen [11]. Thermogravimetric analysis (TGA) is one of the most widely used techniques to study the mechanism and determine the kinetic parameters (A and Ea) of thermal and catalytic cracking of HDPE under different conditions.
The activation energy values found in the literature for the degradation of HDPE in the absence of catalysts, and in an inert atmosphere, typically range from 204 to 473 kJ/mol, with most values clustering around 230–270 kJ/mol [10,12,13,14,15]. The use of catalysts, such as H-Al-MCM-41 [16], H-USY [17,18], H-Y [17,19], and H-ZSM-5 [17,18,19], can significantly lower the activation energy, with values ranging from 60 to 185 kJ/mol, indicating a more facilitated degradation process.
In turn, the chemical conversion of plastics by CC under hydrogen (H2) pressure—hydrocracking—offers many advantages when compared to alternative TC and CC routes. The ability to process a wide range of feedstocks, produce more stable and valuable products, and reduce the coke precursors responsible for deactivating the catalyst are the most significant [20,21]. However, to the best of our knowledge, there are no TGA kinetic studies for the conversion of HDPE conducted in a reducing atmosphere.
Therefore, the main objective of this work is to determine the apparent activation energy (Ea) for the thermal and the catalytic conversion of HDPE under H2 over parent and metal-based H-USY (15) and H-ZSM-5 (11.5) zeolites. Two integral iso-conversional methods, based on the analysis of multiple TGA curves measured at different heating rates, were applied, namely, Flynn-Wall-Ozawa and Kissinger–Akahira–Sunose. The impact of incorporating a metal function (Ni, Pt, or Pd) to the parent zeolites on the apparent Ea, associated with the multistep processes occurring in HDPE hydrocracking is examined and discussed.
Therefore, the primary objective of this study is to ascertain the apparent activation energy (Ea) associated with the thermal and catalytic conversion of HDPE under a H2 atmosphere. The catalytic systems under investigation include the H-USY (15) and H-ZSM-5 (11.5) zeolites. The impact of incorporating a metal function (Ni, Pt, or Pd) into the parent zeolites on the apparent Ea associated with the multistep process occurring in HDPE hydrocracking is examined and discussed. Two integral iso-conversional methods, based on the analysis of multiple TGA curves measured at different heating rates, were applied, namely, Flynn-Wall-Ozawa (FWO) and Kissinger–Akahira–Sunose (KAS)

2. Theoretical Approach

Several methods are used in the literature for non-isothermal solid state-kinetics calculations based on TGA data. These methods fall into two major categories: model-fitting and model-free or iso-conversional methods [22,23,24].
Model-fitting methods are applied to determine Ea and A from a single simulated curve [25]. In this case, TGA data are fitted to different reaction models (Table S1 of the Supporting Information) in order to obtain the best statistical fit and consequently the most suitable kinetic parameters (Ea and A). The main advantage of model-fitting methods, when applied to solid-state reactions, is their experimental simplicity, since only a single TGA experiment is required to determine the kinetic data. However, model-fitting methods also present several problems, the most important of which is the inability to unequivocally identify the unknown reaction model [26].
Conversely, model-free or iso-conversional methods allow the determination of the Arrhenius parameters without pre-selecting a reaction model and a reaction order. However, from an experimental point of view, this approach is much more time-consuming as it requires TGA curves at different heating rates to obtain a dependence of the kinetic parameters on conversion and/or temperature [27]. According to the International Confederation for Thermal Analysis and Calorimetry (ICTAC) Kinetics Committee [28], the use of multi-rate programs is recommended for the calculation of reliable kinetic parameters. This methodology not only avoids the risk of obtaining erroneous kinetics, due to an inappropriate reaction model but also allows the prediction of some mechanistic steps, since Ea can be estimated for a range of conversion values [29,30]. Several authors have contributed to the extension and consolidation of the knowledge in IC methods, such as Friedman [31], Flynn and Wall [32,33], Ozawa [34], Kissinger [35], and Vyazovkin [28,36], among others.
According to the authors listed above, the rate of a chemical reaction, in this specific case, the rate of HDPE degradation, can be described as follows:
d x d t = K T f x
where K(T) is the kinetic constant rate and f x is a function depending on the reaction mechanism involved. The term x is defined as the advancement of the reaction and is expressed in terms of conversion according to the following equation:
x = x 0 x t x 0 x f
where x 0 and x f are the initial and final solid mass, respectively; and x t the mass at a specific time t. In its turn, the temperature-dependent term, K(T) is generally denoted as:
K T = A   e x p ( E a R T )
where A is the frequency factor (min−1), Ea is the apparent Activation Energy (J/mol), R is the gas constant (8.314 J/K mol), and T is the absolute temperature. The generalized kinetic rate equation is expressed in terms of x and obtained by the combination of the three previous equations:
d x d t = K T f x = A   e x p ( E a R T ) f ( x )
For a dynamic TGA process, the heating rate can be defined as β =   d T /dt. Including β in the previous equation results in the following relationship:
d x d t = 1 β A   e x p ( E a R T ) f ( x )
Therefore, the integral form of the inverse of f x , which describes the reaction kinetic model, can be expressed as follows:
G x = 0 x d x f x = A β T 0 T exp E a R T d T = A E a β R   p ( u )
the term p(u) represents a temperature-dependent integral, without an analytical solution
p u = 0 u exp u u 2 d u
with
u = E a R T
However, Equation (6) can be solved by using different approximations for the p(u) function. Accordingly, several integral model-free iso-conversional methods have been developed for the estimation of the kinetic parameters (Ea and A).
For these methods, p(u) can be expressed by the following general approximate equation:
p u = e ( C   u D ) u B
where the values of B, C, and D depend on the method applied.
By replacing Equation (9) in Equation (6), linearizing and simplifying, it is possible to obtain a general expression that relates ln(βi/TBx,I) with 1/Tx,i:
ln β i T x ,   i B = ln A G ( x ) D + 1 B ln ( E a R ) C ( E a R T x , i )
This means that, at any constant conversion, a linear relationship is obtained by plotting ln(βi/TBx,I) against 1/Tx,i. Consequently, from the slope of the line [−CEa/R], it is possible to calculate the Ea value for the corresponding extent of conversion, while from the line interception [ln[A/G(x)] − D + (1 − B) ln(Ea/R)], a combination of the pre-exponential factor with G(x), [ln A/G(x)] [37] is obtained.
In this study, the Flynn–Wall–Ozawa (FWO) and Kissenger–Akahira–Sunose (KAS) methods were employed to ascertain Ea. In accordance with the FWO method [32], which employs Doyle’s approximation [38], the values of B, C, and D are 0, 1.0516, and 5.3305, respectively. In the case of the KAS method, which employs the more accurate approximation of Murray and White [39], the values of B, C, and D are 2, 1, and 0, respectively. Substituting B, C, and D by the respective values into Equation (10), yields the final expressions for each method, as outlined in Table 1.
It is important to note that the accuracy of integral iso-conversional methods is constrained by the underlying assumptions upon which they are based. Alternative differential iso-conversional methods, such as the Friedman (FR) method, have been developed to address these limitations. However, the application of the FR method necessitates numerical differentiation of typically noisy TGA data, which can give rise to other issues, including numerical instabilities and discontinuous Ea profiles.

3. Results

3.1. Catalyst Characterization

The PXRD diffractograms of the parent and metal-based H-USY (15) and H-ZSM-5 (11.5) zeolites are displayed in Figures S1 and S2 of the Supporting Information. For the H-USY zeolite, peaks at 2 θ = 6.2, 10.3, 12.1, 15.9, 18.9, 20.7, 24.0, and 27.5° are observed. Such peaks are characteristic of the FAU framework in USY zeolites [40]. In contrast, for H-ZSM-5 (11.5), diffraction peaks are observed at 2 θ = 7.9, 8.8, 14–17, and 23–25°, which are characteristic of the MFI framework in ZSM-5 zeolites [41] After incorporation of Ni, Pt, and Pd no signs of structural changes are observed, as the metal-based zeolites show the characteristic peaks of the parent structure. However, in the case of Ni-based H-USY (15) and H-ZSM-5 (11.5) new diffraction peaks are observed at 2 θ = 7.9 and 8.8, typical of Ni0 species. Pt0 and Pd0 reflection peaks are not observed in the PXRD diffractograms as the metal content (0.5 wt.%) is below the detection limit of the instrument. A slight decrease (≤10%) [Table 2] in the crystallinity of the parent systems is detected for all the metal-based systems, as a result of metal incorporation.
In terms of textural properties, the H-USY (15) zeolite exhibits a higher external surface area (Sext) and a higher microporous volume (Vmicro), due to its small crystal sizes (0.7 μm) and its microporous structure with a secondary mesoporosity, respectively. According to the literature [42], H-USY-type zeolites have a three-dimensional pore structure in which the basic sodalite structural units are assembled in spherical supercages with a diameter of 1.3 nm and an aperture window of 0.7 nm. Conversely, the HZSM-5 zeolite is composed of a three-dimensional structure interconnected by two channels with 10-membered rings with the following dimensions: 0.51 × 0.55 nm (sinusoidal channels) and 0.56 × 0.53 nm (straight channels) [43]. Upon metal incorporation a slight reduction in Sext and Vmicro is observed, especially for H-ZSM-5 (11.5) zeolites, indicating the deposition of Ni, Pt, and Pd particles on the external surface of the zeolite and/or inside its porous structure. Nevertheless, H-ZSM-5 (11.5) reveals a higher number of total acid sites, 1270 μmol/g, than H-USY (15) zeolite, 651 μmol/g. However, in both systems, the addition of a metallic source leads to a slight decrease in total acidity, probably due to the partial blockage of the pores opening by the metallic particles, which limits the access of the NH3 to the active sites. This fact is corroborated by N2 sorption results.

3.2. Thermogravimetric Analysis

In order to evaluate the decomposition of HDPE in the presence and absence of H-ZSM-5 (11.5) and H-USY (15) zeolites, the TGA profiles at a constant heating rate (10 °C/min) are displayed in Figure 1a, while the temperatures at which HDPE mass loss is 5, 50, and 95% (T5%, T50%, and T95%) are shown in Figure 1b.
According to the data, the HDPE thermal degradation occurs in one single step, starting at 433 °C and ending at 488 °C. Similar HDPE degradation profiles are reported in the literature under inert conditions, as well as identical degradation ranges (between 381 and 501 °C) [44,45,46]. Instead, the catalyst seems to play a crucial role in the degradation. In this case, a substantial decrease in the energy requirements is observed upon the addition of H-USY (15) and H-ZSM-5 (11.5). Although both catalysts reduce HDPE degradation temperatures, H-USY (15) is the one that allows a higher diminishment, around 162 °C on T5%. For the same mass loss, the H-ZSM-5 (11.5) allows for a reduction of 115 °C. According to the literature [47,48], the catalytic activity of zeolites in this type of reaction is strongly dependent on the acidity and the accessibility of the acid sites of the catalyst. Although the H-USY (15) zeolite exhibits a lower acid content compared to H-ZSM-5 (11.5), it has a three-dimensional channel system with large cavities, also known as supercages. These supercages, together with a higher external surface area, reduce the diffusional limitation constraints, and favor the access of polymer macromolecules to the active sites. These facts lead to an increase in cracking reactions and consequently to an easy degradation of HDPE [8].
Together with the physical and chemical features, the heating rate (β) used in the TGA experiments also has an important influence on HDPE degradation. The effect of β (5,10, 15, and 20 °C/min) on the HDPE TGA profiles and first-order derivative (DTG) curves for the thermal and catalytic systems is displayed in Figure 2.
The T5%, T50%, and T95% and the maximum weight loss temperature (Tm), obtained from the peak of the DTG curve, are summarized in Table 3.
Data reveal that in the absence of a catalyst, an increase in β delays the degradation rate of HDPE, leading to a shift of T5%, T50%, and T95% to higher temperatures. For the extreme cases of β (5 and 20 °C/min), the HDPE degradation range varies from 433 to 473 °C and from 463 to 501 °C, respectively. The DTG results also show that the maximum loss rate occurs at 468 °C (β = 5 °C/min) and 493 °C (β = 20 °C/min), indicating heat and mass transfer limitations during the thermal run for higher heating rates [49,50].
A shift in the mass loss curve to a higher temperature as the heating rate increased was also observed when both catalysts were used. In the case of H-USY (15), the degradation temperature range suffers a significant increase with β, varying from 255–398 °C at 5 °C/min to 292–431 °C at 20 °C/min. For H-ZSM-5 (11.5), a change of the same order is observed. For both catalytic systems, although the TGA profiles present an identical shape for different β values, the DTG data reveal the appearance of an additional peak at lower heating rates. In previous studies, the authors concluded that at sufficiently low decomposition rates, overlapping adjacent peaks can be separated due to lower heat transfer limitations. In this case, the different decomposition reactions that occur simultaneously at higher β can be identified separately at lower β [51].

3.3. Thermal and Catalytic Cracking Kinetics

Kinetic analysis plays a main role in reactor design and scale-up processes, providing insight into the reaction mechanism [52]. The objective of kinetic analysis is the interpretation of the experimentally determined kinetic parameters (Ea and A), each of which is associated with a fundamental theoretical concept. Ea represents the minimum amount of energy required for a reaction to occur and is strongly dependent on the nature and concentration of the reactants, the reaction conditions, and the presence or absence of a catalyst and its type [53]. In contrast, A is related to the frequency of vibrations of the activated complex.

3.3.1. Apparent Activation Energy (Ea)

The Ea values for both the thermal and catalytic degradations of HDPE under a reducing atmosphere were calculated by applying the FWO, and KAS methods to the previous TGA data at different β values.
Figure 3 depicts the fitting results obtained from the model-free methods plots for the thermal degradation of HDPE. The resulting Ea and A values and correlation coefficient data are listed in Table 4 for thermal and catalytic conditions.
Regarding the thermal data, correlation coefficients (R2) higher than 0.999 are observed, indicating an excellent fit of the method to the experimental data and an accurate estimation of Ea. Moreover, it is observed that the fitted lines are approximately parallel, suggesting a unification of multiple reactions mechanism [29,54]. In this case, the Ea varies between approximately 211 and 234 kJ/mol for both the FWO and KAS methods as the HDPE conversion increases. The mean Ea value for all methods is 227 kJ/mol.
To the best of our knowledge, this is the first study reported in the literature to employ non-isothermal methods to determine the Ea for the thermal and catalytic degradation of HDPE in a reductive H2 atmosphere. It should be noted that the present results fall within the range reported by other authors for the Ea of the thermal degradation of HDPE in an inert atmosphere (203–250 kJ/mol) [10,12,13,15,55].
When a catalyst is added to the reaction, the linear correlation coefficients are lower than those obtained under thermal conditions (0.895 ≤ R2 ≤ 0.973), indicating a less accurate estimation. This could be due to: (i) slight differences in the amount of the catalyst used to perform the TGA experiments (~2–5%), which could affect the TGA profiles; or, (ii) poorer model fit for more complex mechanisms. As expected, the presence of H-USY (15) and H-ZSM-5 (11.5) zeolites leads to a significant reduction in Ea, since the reaction rate increases. With regard to the H-USY (15) data, it can be observed that the Ea values vary with the extent of the reaction. The values range from 78 to 157 kJ/mol for the KAS methods and from 83 to 172 kJ/mol for the FWO method. In the case of the H-ZSM-5 (11.5) zeolite, the Ea values vary between 53 and 122 kJ/mol for the KAS method and between 60 and 107 kJ/mol for the FWO method. The differences between the ranges of Ea for these two methods are more obvious for this catalyst, probably due to the worst linear fitting obtained for the iso-conversion plots (Table 4). The average Ea values are, however, quite similar, at 77 kJ/mol and 81 kJ/mol, respectively.
It is also interesting to note that the Ea observed for the catalytic conversion of HDPE over H-ZSM-5 (11.5) zeolite is significantly lower than over H-USY (15). Although H-USY (15) exhibits more suitable structural and textural properties for the reaction, H-ZSM5 zeolite has a higher number of acid sites, which are the main responsible for the HDC reactions [18]. The Ea values obtained are in agreement with those reported in the literature for both systems [H-USY (15): 95–118 kJ/mol and H-ZSM-5 (11.5): 60–124 kJ/mol] under an inert atmosphere [17,18,19].
In addition, data show that an Ea dependency with the extent of reaction is observed either for thermal or catalytic processes, suggesting that in all cases, multi-step complex mechanisms are involved [11,23], as previously observed. Similar behavior was observed by other authors [9,56,57], who reported an increase in the effective Ea with the course of the reaction for PE, as a result of the change in the limiting step during the reaction. For the particular case of thermal run, it is well known that the degradation occurs via a random scission mechanism followed by a radical transfer process [46,58]. Nevertheless, according to Peterson et al. [56] and Vyazovkin et al. [57], the lower Ea value observed in the early stages of the process is most likely related to the ignition of the degradation process, occurring at thermally labile bonds such as chain branches, unsaturated compounds, peroxides, and carbonyl groups, formed during the manufacture, storage, and processing of HDPE. Once the weak links are consumed, the limiting step of the degradation shifts towards the degradation initiated by random scission, leading to an increase in Ea. Instead, and although the Ea also increases with the extent of conversion, a distinct trend is observed for the catalytic systems, suggesting that a different mechanism is involved. This is corroborated by literature data, indicating that in the presence of a catalyst, HDPE degradation proceeds by the carbocation mechanism [59].

3.3.2. Effect of Adding a Metallic Source (Ni, Pt, Pd) to H-USY (15) and H-ZSM-5 (11.5) on the Kinetic Parameters

The catalytic conversion under a reducing atmosphere can proceed by distinct pathways depending on the catalyst used. In the presence of a monofunctional acidic catalyst, such as H-ZSM-5 (11.5) and H-USY (15), the reaction takes place by the carbocation mechanism, as mentioned before. However, when a metallic function is introduced on a monofunctional zeolite, the reaction mechanism changes. In this case, the metal favors dehydrogenation and hydrogenation reactions, while the remaining acidic centers are responsible for the cracking and isomerization reactions [60]. To evaluate the effect of the addition of a metal source on the catalytic degradation of HDPE in a H2 atmosphere, Ni, Pt, and Pd were added to H-USY (15) and H-ZSM-5 (11.5) zeolites. The TGA profiles are shown in Figure S3 of the Supporting Information.
Data show a positive effect on the degradation temperatures with the addition of Ni, Pt, and Pd either to H-USY (15) or to H-ZSM-5 (11.5), suggesting that the metallic function enhances the catalytic conversion of HDPE in the presence of a reducing H2 atmosphere. This is due to the formation of olefinic intermediates on the metallic sites, which are easily converted into carbocation intermediates, contributing to an increase in the cracking rate. During the catalytic conversion over a monofunctional catalyst, these olefinic intermediates are not formed, thus reducing the rate of carbocation formation. However, the reduction in HDPE degradation temperatures is more pronounced for the noble metal than for the transition metal catalysts, leading to lower degradation temperatures according to the following order: Pd < Pt < Ni. The promising results obtained with noble metal systems are probably due to two main facts: (i) higher hydrogenation ability, which promotes the dehydrogenation/hydrogenation reactions [61]; and (ii) smaller particle size and consequently, higher dispersion of the metal on the H-USY (15) and H-ZSM-5 (11.5) supports (Figure S4 and Table S1 of the Supporting Information). However, from an economic point of view, the use of a Ni metallic source could be more favorable.
In order to understand how the presence of a metal source and its nature affects the kinetic parameters of the process, the KAS method was applied to TGA data at different heating rates (β = 5, 10, 15, and 20 °C/min) using Ni-, Pt-, and Pd-based zeolites as catalysts. The HDPE degradation profiles at different β values are shown in Figure S3, while the resulting Ea and A values are summarized in Table 5.
As previously observed for parent zeolites, an increase in β leads to a shift in the degradation profiles of all the metal-based zeolites to higher temperatures, as a result of the reduced heat transfer and shorter exposure time of HDPE at a given temperature. Furthermore, higher heating rates also change kinetics, resulting in a slower degradation of HDPE macromolecules [24]. Data available in Table 5 also show that in the presence of bifunctional catalysts, lower Ea values are observed compared to the thermal run data, as previously observed for the monofunctional systems. In general, a slight increase in the average Ea value is observed for the metal-based catalysts compared to the respective parent zeolites. This is probably the result of two main aspects, namely, (i) the contribution of the dehydrogenation/hydrogenation reactions promoted by metal source to the mean Ea, and (ii) the reduction of the total acidity of H-USY (15) and H-ZSM-5 (11.5) zeolites upon the impregnation process, both of which are essential to the cracking and isomerization steps [60].
Another explanation was suggested by the studies of Ramdoss and Tarrer [62], who developed a kinetic model combining series and parallel reactions for the conversion of post-consumer plastics to light oil, heavy oil, gas, and coke, in a H2 atmosphere. Looking at the activation energies obtained for the various reaction pathways, it was found that the highest Ea values, 121.7 kJ/mol and 113.5 kJ/mol, correspond to the conversion of light and heavy liquids into gaseous fractions, respectively, while the lowest Ea value (11.8 kJ/mol) results from the conversion of heavier hydrocarbons into liquid fractions. Recent studies performed by our team show that there is a greater tendency to produce gaseous products over Ni-, Pt-, and Pd-based zeolites compared to their respective parent zeolites, which could thus contribute to the increase in Ea over the bifunctional catalysts. Miskolczi et al. [63], who calculated the Ea for the thermal and catalytic degradation of end-of-life vehicle plastic waste under a N2 atmosphere, over H-Y and H-ZSM-5 (11.5) loaded with distinct metals (Fe, Sn, Ce, Mg, Ni, Cu, and Zn), also found that the addition of a metal source to the parent zeolites results in an increase in the apparent Ea of the reaction for all cases, in the order Cu < Ce < Mg < Ni < Fe(III) < Fe(II) < Zn < Sn, thus, corroborating our data.

4. Materials and Methods

4.1. Materials

Commercial high-density polyethylene (HDPE, MW = 155,000 g/mol; MWD = 5.4; d = 0.95 g/cm3 and Tm = 140 °C), was supplied by Repsol (Sines, Portugal) without additives and used in powder form.
The zeolite catalysts, NH4-ZSM-5(11.5) (CBV-2314) and H-USY (15) (CBV-720), were supplied by Zeolyst (Conshohocken, PA, USA). The metal salts, nickel nitrate hexahydrate (Ni(NO3)2·6H2O), and tetraamine platinum (II) nitrate (Pt (NH3)4·(NO3)2), used for the preparation of Ni- and Pt-impregnated catalysts were supplied by Merck (Darmstadt, Germany), both with a purity of 99%. The palladium nitrate (Pd(NO3)2) solution used for the preparation of Pd-impregnated zeolites was supplied by Sigma-Aldrich (St. Louis, MO, USA).

4.2. Catalyst Preparation

NH4-ZSM-5 was calcined at 500 °C for 6 h under an airflow of 4 L·h−1.g−1 in order to transform the ammonium form to the protonic form (H-ZSM-5 (11.5)). A plateau at 200 °C, for 1 h, was carried out to remove the water and prevent the zeolite from steaming. A heating rate of 2 °C·min−1 was used.
The preparation of Ni-, Pt-, and Pd-based zeolites was performed by the incipient wetness impregnation method, using an aqueous solution of Ni(NO3)2·6H2O, Pt(NH3)4·(NO3)2 and Pd(NO3)2, respectively. The volume of the metal precursor solution used was calculated to match the volume of the zeolite pores and was added dropwise to the support. The impregnated samples were air-dried at 80 °C for 24 h and then calcined at 500 °C in an airflow of L·h−1·g−1. Prior to use, all impregnated catalysts were also pre-activated in a reactor with an H2 flow of 4 L·h−1·g−1 for 2 h, at 500 °C for Ni catalysts and at 450 °C for Pt and Pd catalysts.

4.3. Catalyst Characterization

The textural properties of the parent and metal-based catalysts were assessed by N2 sorption. Measurements were performed at −196 °C using an Autosorb IQ apparatus from Quantachrome (Boynton Beach, FL, USA). Prior to measurement, the materials were degassed under vacuum at 90 °C for 1 h and then heated at 350 °C for 5 h. The external surface area (Sext) and micropore volume (Vmicro) were calculated using the t-plot method, while the total pore volume was determined from the adsorbed volume of nitrogen at a relative pressure (P/P0) of 0.95. The difference between Vtotal and Vmicro gives the mesopore volume (Vmeso).
Powder X-ray diffraction (PXRD) analysis was used to obtain information on the phases present in the prepared catalysts, as well as to identify any damage to the zeolite materials upon metal incorporation. The PXRD experiments were carried out using a Bruker AXS Advance D8 diffractometer (Billerica, MA, USA) equipped with a 1D detector (SSD 160) and using a Ni filter, with a CuKα radiation source (λ = 1.5406 nm) and operating at 40 kV and 30 mA. The scanning range was from 5° to 80° (2 theta), with a step size of 0.03° and a step time of 2 s.
The density and strength distribution of acid sites of the zeolites were determined by ammonia temperature-programmed desorption (NH3-TPD). In these experiments, the sample was pre-treated under helium atmosphere at 350 °C for 1 h and then cooled to 125 °C and ammonia-saturated in a stream of 15% NH3/He at a flow of 30 mL/min, for 1 h. Prior to the desorption step, the sample was outgassed under helium at 125 °C for 30 min, to remove the physiosorbed ammonia. Chemisorbed NH3 was desorbed at a heating rate of 10 °C/min until 700 °C and the amount of ammonia desorbed in the effluent stream was detected by a thermal conductivity detector (TCD). The total acidity of the catalyst was obtained by integrating the area under the desorption curve.
TEM analysis was carried out to determine the size of Ni0, Pt0, and Pd0 particles on reduced samples. The experiments were carried out using an HRTEM 2010 JEOL LaB6 microscope (200 kV) (Tokyo, Japan). Data analysis was performed by using ImageJ software (Version 1.54). For this purpose, a minimum of 150 or 300 metallic particles were counted in each sample, depending on the metal source, transition or noble metal, respectively. The dispersion of metal-based catalysts, which is the ratio between the total number of metal atoms present on the surface of the catalyst and the total number of metal atoms, was determined assuming that Ni, Pt, and Pd are spherical particles, through the following equation:
D = 6/d.p,
where d is the mean particle size and p is the atomic density.

4.4. Thermogravimetric Analysis Procedure

The samples used for the TGA experiments were prepared by compression molding technique. The HDPE and catalyst, both in powder form, were first mechanically mixed in a mass ratio of 9/1 and then heated in a hydraulic press at 140 °C and 3 tons of pressure for 5 min [64].
The thermogravimetric experiments were performed using a Setaram TGA 92-16.18 thermobalance (Austin, TX, USA) under a H2 atmosphere, at a continuous flow rate of 30 mL/min. The temperature varied between 20 and 700 °C, at different heating rates (2, 5, 10, 15 °C/min). To avoid the presence of oxygen, a nitrogen purge of 15 min was performed before each experiment.

5. Conclusions

In this work, the kinetic parameters of the thermal and catalytic conversions of HDPE under a H2 atmosphere were determined by non-isothermal TGA experiments, applying two different model-free methods, namely Flynn–Wall–Ozawa and Kissinger–Akahira–Sunose. The results show that the thermal degradation of HDPE at the lowest heating rate takes place in the range of 433–473 °C, which means that a high energy input is required. Nevertheless, the addition of H-USY (15) and H-ZSM-5 (11.5) zeolites has a positive effect on the degradation process, allowing the temperature at which mass loss is 5% to be reduced by 162 and 115 °C, respectively.
The kinetic data for HDPE thermal degradation show that Ea values increase with conversion (from 211 to 234 kJ/mol) and that an average value of 227 kJ/mol is obtained regardless of the model-free method used (FWO or KAS).
In the presence of H-USY (15) and H-ZSM-5 (11.5) catalysts, Ea suffers a significant decrease and becomes highly dependent on conversion. However, the average Ea value obtained for the FWO and KAS methods is very similar. For H-USY (15), it ranges from 116 to 122 kJ/mol, while for H-ZSM-5, it is much lower and varies between 77 and 81 kJ/mol. The variation in Ea with the extent of reaction is an indication of the presence of a multi-step complex reaction mechanism.
The addition of a metal source to the H-USY (15) and H-ZSM-5 (11.5) zeolites also leads to a significant reduction in the apparent Ea of the reaction compared to the thermal process. Conversely, a slight increase in Ea is observed compared to the respective parent systems. This is probably related to the contribution of the dehydrogenation/hydrogenation reactions to the overall mechanism and thus, to the resulting Ea.
In conclusion, both parent and metal-based H-USY (15) and H-ZSM-5 (11.5) zeolites can greatly reduce the energy barrier and corresponding energy consumption for HDPE hydrocracking compared to thermal conditions, thus promoting a more sustainable conversion of plastic waste.
As a final point, it is also important to recall that HDC is a complex reaction, involving a multi-step mechanism. The elucidation of the mechanisms underlying complex processes is very challenging, and extensive kinetic studies and model simulations are needed to gain a deeper understanding. In this preliminary investigation on the kinetics of HDPE hydrocracking, a meaningful estimate of the average value for the apparent Ea associated with the multistep process of each of the thermal and catalytic systems studied was obtained. However, to gain deeper mechanistic insights into the HDC reaction, it would be beneficial to apply advanced iso-conversional methods (AIC), which use numerical integration over small temperature intervals as a component of Ea evaluation, in association with a model-free methodology to determine lnA. Further studies are in progress to address this topic.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/catal14080514/s1, Figure S1: PXRD diffractograms of parent and metal-based (Ni, Pt, Pd) H-USY (15) (A) and identification of Ni0, Pt0 and Pd0 species on XRD diffractograms (B); Figure S2: PXRD diffractograms of parent and metal-based (Ni, Pt, Pd) H-ZSM-5 (11.5) (A) and identification of Ni0, Pt0 and Pd0 species on XRD diffractograms (B); Figure S3- TGA profiles at distinct β (5, 10, 15 and 20 °C/min) for Ni-, Pt-, and Pd-based HUSY (A,C,E) and H-ZSM-5 (B,D,F) zeolites. Figure S4: TEM images of Ni-, Pt-, and Pd-based H-USY (15) and H-ZSM-5 (11.5) zeolites. Table S1: Metallic properties of metal-based H-USY (15) and H-ZSM-5 zeolites (11.5).

Author Contributions

Conceptualization, C.S.C., M.R.R., J.M.S. and A.F.; methodology, C.S.C. and A.F.; validation, C.S.C., A.F. and M.M.; formal analysis, C.S.C.; investigation, C.S.C., A.F. and M.M.; writing—original draft preparation, C.S.C.; writing—review and editing, C.S.C., M.R.R. and J.M.S.; visualization, C.S.C.; supervision, M.R.R. and J.M.S.; project administration, M.R.R. and J.M.S.; funding acquisition, M.R.R. and J.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundação para a Ciência e Tecnologia (FCT) [Project UIDB/00100/2020 (DOI 10.54499/UIDB/00100/2020), Project UIDP/00100/2020 (DOI 10.54499/UIDP/00100/2020), Project LA/P/0056/2020 (DOI 10.54499/LA/P/0056/2020), and (ERANETMED/0003/201). Cátia S. Costa is grateful for the PhD fellowship of the FCT-CATSUS program (PD/BD/12862/2017). The work is also part of the TED2021-129688B project, funded by MCIN/AEI/10.13039/501100011033 and by the European Union “NextGenerationEU”/PRTR. The content of this publication does not reflect the official opinion of the European Union. Responsibility for the information and views expressed in this paper lies entirely with the authors.

Data Availability Statement

Data are contained within the article or Supplementary Material.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mrowiec, B. Plastics in the circular economy (CE). Environ. Prot. Nat. Resour. J. Inst. Environ. Prot. Res. Inst. 2018, 29, 16–19. [Google Scholar] [CrossRef]
  2. Al-Salem, S.M.; Lettieri, P.; Baeyens, J. Recycling and recovery routes of plastic solid waste (PSW): A review. Waste Manag. 2009, 29, 2625–2643. [Google Scholar] [CrossRef] [PubMed]
  3. Al-Salem, S.M.; Antelava, A.; Constantinou, A.; Manos, G.; Dutta, A. A review on thermal and catalytic pyrolysis of plastic solid waste (PSW). J. Environ. Manag. 2017, 197, 177–198. [Google Scholar] [CrossRef] [PubMed]
  4. Mohanraj, C.; Senthilkumar, T.; Chandrasekar, M. A review on conversion techniques of liquid fuel from waste plastic materials. Int. J. Energy Res. 2017, 41, 1534–1552. [Google Scholar] [CrossRef]
  5. Panda, A.K. Thermo-catalytic degradation of different plastics to drop in liquid fuel using calcium bentonite catalyst. Int. J. Ind. Chem. 2018, 9, 167–176. [Google Scholar] [CrossRef]
  6. Ratnasari, D.K.; Nahil, M.A.; Williams, P.T. Catalytic pyrolysis of waste plastics using staged catalysis for production of gasoline range hydrocarbon oils. J. Anal. Appl. Pyrolysis 2017, 124, 631–637. [Google Scholar] [CrossRef]
  7. Almeida, D.; de Fátima Marques, M. Thermal and Catalytic Pyrolysis of Polyethylene Plastic Waste. Polimeros 2015, 26, 1–8. [Google Scholar] [CrossRef]
  8. Socci, J.; Osatiashtiani, A.; Kyriakou, G.; Bridgwater, T. The catalytic cracking of sterically challenging plastic feedstocks over high acid density Al-SBA-15 catalysts. Appl. Catal. A Gen. 2019, 570, 218–227. [Google Scholar] [CrossRef]
  9. Khedri, S.; Elyasi, S. Kinetic analysis for thermal cracking of HDPE: A new isoconversional approach. Polym. Degrad. Stab. 2016, 129, 306–318. [Google Scholar] [CrossRef]
  10. Kumar, S.; Singh, R.K. Pyrolysis kinetics of waste high-density polyethylene using thermogravimetric analysis. Int. J. ChemTech Res. 2014, 6, 131–137. [Google Scholar]
  11. Kumar, A.; Mylapilli, S.V.P.; Reddy, S.N. Thermogravimetric and kinetic studies of metal (Ru/Fe) impregnated banana pseudo-stem (Musa acuminate). Bioresour. Technol. 2019, 285, 121318. [Google Scholar] [CrossRef] [PubMed]
  12. Kayacan, İ.; Doğan, Ö.M. Pyrolysis of Low and High Density Polyethylene. Part I: Non-isothermal Pyrolysis Kinetics, Energy Sources, Part A Recover. Util. Environ. Eff. 2008, 30, 385–391. [Google Scholar] [CrossRef]
  13. Encinar, J.M.; González, J.F. Pyrolysis of synthetic polymers and plastic wastes. Kinetic study. Fuel Process. Technol. 2008, 89, 678–686. [Google Scholar] [CrossRef]
  14. Fernandes, G.J.T.; Fernandes, V.J.; Araujo, A.S. Catalytic degradation of polyethylene over SAPO-37 molecular sieve. Catal. Today 2002, 75, 233–238. [Google Scholar] [CrossRef]
  15. Budrugeac, P. Theory and practice in the thermoanalytical kinetics of complex processes: Application for the isothermal and non-isothermal thermal degradation of HDPE. Thermochim. Acta 2010, 500, 30–37. [Google Scholar] [CrossRef]
  16. Araújo, A.S.; Fernandes, V.J.; Araujo, S.A.; Ionashiro, M. Kinetic evaluation of the pyrolysis of high density polyethylene over H-AlMCM-41 material. Stud. Surf. Sci. Catal. 2002, 141, 473–478. [Google Scholar] [CrossRef]
  17. Street, S.; Centre, E.T.; Street, S. Catalytic degradation of high density polyethylene: An evaluation of mesoporous and microporous catalysts using thermal analysis. Thermochim. Acta 1997, 294, 65–69. [Google Scholar]
  18. Garforth, A.; Yeuh-Hui, L.I.N.; Sharratt, P.; Dwyer, J. Catalytic Polymer Degradation for Producing Hydrocarbons over Zeolites. Stud. Surf. Sci. Catal. 1999, 121, 197–202. [Google Scholar] [CrossRef]
  19. Coelho, A.; Costa, L.; Marques, M.D.M.; Fonseca, I.; Lemos, M.A.; Lemos, F. Using simultaneous DSC/TG to analyze the kinetics of polyethylene degradation-catalytic cracking using HY and HZSM-5 zeolites. React. Kinet. Mech. Catal. 2010, 99, 5–15. [Google Scholar] [CrossRef]
  20. Munir, D.; Abdullah; Piepenbreier, F.; Usman, M.R. Hydrocracking of a plastic mixture over various micro-mesoporous composite zeolites. Powder Technol. 2017, 316, 542–550. [Google Scholar] [CrossRef]
  21. Munir, D.; Irfan, M.F.; Usman, M.R. Hydrocracking of virgin and waste plastics: A detailed review. Renew. Sustain. Energy Rev. 2018, 90, 490–515. [Google Scholar] [CrossRef]
  22. Khawam, A.; Flanagan, D.R. Role of isoconversional methods in varying activation energies of solid-state kinetics: II. Nonisothermal kinetic studies. Thermochim. Acta 2005, 436, 101–112. [Google Scholar] [CrossRef]
  23. Ashraf, A.; Sattar, H.; Munir, S. A comparative applicability study of model-fitting and model-free kinetic analysis approaches to non-isothermal pyrolysis of coal and agricultural residues. Fuel 2019, 240, 326–333. [Google Scholar] [CrossRef]
  24. Murichan, N.; Cherntongchai, P. Kinetic Analysis of Thermal Degradation of Polyolefin Mixtures. Int. J. Chem. Eng. Appl. 2014, 5, 169–175. [Google Scholar] [CrossRef]
  25. Sánchez-Jiménez, P.E.; Pérez-Maqueda, L.A.; Perejón, A.; Criado, J.M. Clarifications regarding the use of model-fitting methods of kinetic analysis for determining the activation energy from a single non-isothermal curve. Chem. Cent. J. 2013, 7, 25. [Google Scholar] [CrossRef] [PubMed]
  26. Janković, B.; Adnadević, B.; Jovanović, J. Application of model-fitting and model-free kinetics to the study of non-isothermal dehydration of equilibrium swollen poly (acrylic acid) hydrogel: Thermogravimetric analysis. Thermochim. Acta 2007, 452, 106–115. [Google Scholar] [CrossRef]
  27. Ravari, F.; Noori, M.; Ehsani, M. Thermal Stability and Degradation Kinetic Studies of PVA/RGO Using the Model-fitting and Isoconversional (model-free) Methods. Fibers Polym. 2019, 20, 472–480. [Google Scholar] [CrossRef]
  28. Vyazovkin, S.; Chrissa, K.; Laura, M.; Lorenzo, D.; Koga, N.; Pijolat, M.; Roduit, B.; Sbirrazzuoli, N.; Josep, J. Thermochimica Acta ICTAC Kinetics Committee recommendations for collecting experimental thermal analysis data for kinetic computations. Thermochim. Acta 2014, 590, 1–23. [Google Scholar] [CrossRef]
  29. Yao, F.; Wu, Q.; Lei, Y.; Guo, W.; Xu, Y. Thermal decomposition kinetics of natural fibers: Activation energy with dynamic thermogravimetric analysis. Polym. Degrad. Stab. 2008, 93, 90–98. [Google Scholar] [CrossRef]
  30. Siddiqui, M.N.; Antonakou, E.V.; Redhwi, H.H.; Achilias, D.S. Kinetic analysis of thermal and catalytic degradation of polymers found in waste electric and electronic equipment. Thermochim. Acta 2019, 675, 69–76. [Google Scholar] [CrossRef]
  31. 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. 1964, 6, 183–195. [Google Scholar] [CrossRef]
  32. Flynn, J.H.; Wall, L.A. General Treatment of the Thermogravimetry of Polymers. J. Res. Natl. Bur. Standars-Phys. Chem. 1996, 70, 344–346. [Google Scholar] [CrossRef] [PubMed]
  33. Flynn, J.H. The “temperature integral”—Its use and abuse. Thermochim. Acta 1997, 300, 83–92. [Google Scholar] [CrossRef]
  34. Ozawa, T. A New Method of Analyzing Thermogravimetric Data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef]
  35. Kissinger, H.E. Variation of peak temperature with heating rate in differential thermal analysis. J. Res. Natl. Bur. Stand. 1956, 57, 217. [Google Scholar] [CrossRef]
  36. 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]
  37. Kumar, S.A.P.; Nagarajan, R.; Prasad, K.M.; Anand, B.; Murugavelh, S. Thermogravimetric study and kinetics of banana peel pyrolysis: A comparison of “model-free” methods. Biofuels 2019, 13, 129–138. [Google Scholar] [CrossRef]
  38. Doyle, C.D. Estimating isothermal life from thermogravimetric data. J. Appl. Polym. Sci. 1962, 6, 639–642. [Google Scholar] [CrossRef]
  39. Murray, P.; White, J. Kinetics of the thermal dehydration of clay. Part IV. Interpretation of the differential thermal analysis of the clay minerals. Trans. Br. Ceram. Soc. 1955, 54, 204–238. [Google Scholar]
  40. Wong, S.; Ngadi, N.; Abdullah, T.A.T.; Inuwa, I.M. Catalytic Cracking of LDPE Dissolved in Benzene Using Nickel-Impregnated Zeolites. Ind. Eng. Chem. Res. 2016, 55, 2543–2555. [Google Scholar] [CrossRef]
  41. Hamidzadeh, M.; Ghassemzadeh, M.; Tarlani, A.; Sahebdelfar, S. The effect of hydrothermal impregnation of Ni, Co, and Cu on HZSM-5 in the nitrogen oxide removal. Int. J. Environ. Sci. Technol. 2018, 15, 93–104. [Google Scholar] [CrossRef]
  42. Awala, H.; Gilson, J.P.; Retoux, R.; Boullay, P.; Goupil, J.M.; Valtchev, V.; Mintova, S. Template-free nanosized faujasite-type zeolites. Nat. Mater. 2015, 14, 447–451. [Google Scholar] [CrossRef]
  43. Caldeira, V.P.S.; Santos, A.G.D.; Oliveira, D.S.; Lima, R.B.; Souza, L.D.; Pergher, S.B.C. Polyethylene catalytic cracking by thermogravimetric analysis: Effects of zeolitic properties and homogenization process. J. Therm. Anal. Calorim. 2017, 130, 1939–1951. [Google Scholar] [CrossRef]
  44. Abbasa, A.S.; Saber, M.G. Kinetics of Thermal Pyrolysis of High-Density Polyethylene. Iraqi J. Chem. Pet. Eng. 2018, 19, 13–19. [Google Scholar] [CrossRef]
  45. Zhao, D.; Wang, X.; Miller, J.B.; Huber, G.W. The Chemistry and Kinetics of Polyethylene Pyrolysis: A Process to Produce Fuels and Chemicals. ChemSusChem 2020, 13, 1764–1774. [Google Scholar] [CrossRef]
  46. Anene, A.F.; Fredriksen, S.B.; Sætre, K.A.; Tokheim, L.A. Experimental study of thermal and catalytic pyrolysis of plastic waste components. Sustainability 2018, 10, 3979. [Google Scholar] [CrossRef]
  47. Costa, C.S.; Ribeiro, M.R.; Silva, J.M. Hydrocracking of HDPE over H-USY and mesostructure acidic materials: Accessibility and acidity effects. In Proceedings of the 8th International Symposium on Energy from Biomass and Waste—Venice 2020, Venice, Italy, 16–19 November 2020. [Google Scholar]
  48. Costa, C.S.; Thi, H.D.; Van Geem, K.M.; Ribeiro, M.R.; Silva, J.M. Assessment of acidity and the zeolite porous structure on hydrocracking of HDPE. Sustain. Energy Fuels 2022, 6, 3611–3625. [Google Scholar] [CrossRef]
  49. Shahbeig, H.; Nosrati, M. Pyrolysis of municipal sewage sludge for bioenergy production: Thermo-kinetic studies, evolved gas analysis, and techno-socio-economic assessment. Renew. Sustain. Energy Rev. 2020, 119, 109567. [Google Scholar] [CrossRef]
  50. Das, P.; Tiwari, P. Thermal degradation kinetics of plastics and model selection. Thermochim. Acta 2017, 654, 191–202. [Google Scholar] [CrossRef]
  51. Al-Salem, S.M.; Sharma, B.K.; Khan, A.R.; Arnold, J.C.; Alston, S.M.; Chandrasekaran, S.R.; Al-Dhafeeri, A.T. Thermal degradation kinetics of virgin polypropylene (PP) and PP with starch blends exposed to natural weathering. Ind. Eng. Chem. Res. 2017, 56, 5210–5220. [Google Scholar] [CrossRef]
  52. Mortezaeikia, V.; Tavakoli, O.; Khodaparasti, M.S. A review on kinetic study approach for pyrolysis of plastic wastes using thermogravimetric analysis. J. Anal. Appl. Pyrolysis 2021, 160, 105340. [Google Scholar] [CrossRef]
  53. Dhaundiyal, A.; Singh, S.B.; Hanon, M.M.; Rawat, R. Determination of Kinetic Parameters for the Thermal Decomposition of Parthenium hysterophorus. Environ. Clim. Technol. 2018, 22, 5–21. [Google Scholar] [CrossRef]
  54. Volli, V.; Purkait, M.K. Physico-chemical properties and thermal degradation studies of commercial oils in nitrogen atmosphere. Fuel 2014, 117, 1010–1019. [Google Scholar] [CrossRef]
  55. Sinfrônio, F.S.M.; Santos, J.C.O.; Pereira, L.G.; Souza, A.G.; Conceição, M.M.; Fernandes, V.J.; Fonseca, V.M. Kinetic of thermal degradation of low-density and high-density polyethylene by non-isothermal thermogravimetry. J. Therm. Anal. Calorim. 2005, 79, 393–399. [Google Scholar] [CrossRef]
  56. Peterson, J.D.; Vyazovkin, S.; Wight, C.A. Kinetics of the thermal and thermo-oxidative degradation of polystyrene, polyethylene and poly(propylene). Macromol. Chem. Phys. 2001, 202, 775–784. [Google Scholar] [CrossRef]
  57. Vyazovkin, S.; Sbirrazzuoli, N. Isoconversional kinetic analysis of thermally stimulated processes in polymers. Macromol. Rapid Commun. 2006, 27, 1515–1532. [Google Scholar] [CrossRef]
  58. Singh, B.; Sharma, N. Mechanistic implications of plastic degradation. Polym. Degrad. Stab. 2008, 93, 561–584. [Google Scholar] [CrossRef]
  59. Lin, Y.H.; Hwu, W.H.; Ger, M.D.; Yeh, T.F.; Dwyer, J. A combined kinetic and mechanistic modelling of the catalytic degradation of polymers. J. Mol. Catal. A Chem. 2001, 171, 143–151. [Google Scholar] [CrossRef]
  60. Weitkamp, J. Catalytic Hydrocracking-Mechanisms and Versatility of the Process. ChemCatChem 2012, 4, 292–306. [Google Scholar] [CrossRef]
  61. Tao, M.; Meng, X.; Xin, Z.; Bian, Z.; Lv, Y.; Gu, J. General Synthesis and characterization of well dispersed nickel-incorporated SBA-15 and its high activity in syngas methanation reaction. Appl. Catal. A Gen. 2016, 516, 127–134. [Google Scholar] [CrossRef]
  62. Ramdoss, P.K.; Tarrer, A.R. High-temperature liquefaction of waste plastics. Fuel 1998, 77, 293–299. [Google Scholar] [CrossRef]
  63. Miskolczi, N.; Juzsakova, T. Preparation and application of metal loaded ZSM-5 and y-zeolite catalysts for thermo-catalytic pyrolysis of real end of life vehicle plastics waste. J. Energy Inst. 2019, 92, 118–127. [Google Scholar] [CrossRef]
  64. Costa, C.S.; Muñoz, M.; Ribeiro, M.R.; Silva, J.M. A thermogravimetric study of HDPE conversion under a reductive atmosphere. Catal. Today 2020, 379, 192–204. [Google Scholar] [CrossRef]
Figure 1. TGA profiles (a) and temperature at which mass loss is 5, 50, and 95% (T5%, T50%, and T95%) (b), in presence and absence of H-USY (15) and H-ZSM-5 (11.5) zeolites.
Figure 1. TGA profiles (a) and temperature at which mass loss is 5, 50, and 95% (T5%, T50%, and T95%) (b), in presence and absence of H-USY (15) and H-ZSM-5 (11.5) zeolites.
Catalysts 14 00514 g001
Figure 2. TGA and DTG profiles for thermal and catalytic degradation of HDPE under a H2 atmosphere without catalyst (a) and over H-USY (15) (b) and H-ZSM-5 (11.5) (c).
Figure 2. TGA and DTG profiles for thermal and catalytic degradation of HDPE under a H2 atmosphere without catalyst (a) and over H-USY (15) (b) and H-ZSM-5 (11.5) (c).
Catalysts 14 00514 g002
Figure 3. Linear fitting curves under distinct conversion levels for model-free methods applied to the HDPE thermal degradation TGA data.
Figure 3. Linear fitting curves under distinct conversion levels for model-free methods applied to the HDPE thermal degradation TGA data.
Catalysts 14 00514 g003
Table 1. Model-free methods applied in this study.
Table 1. Model-free methods applied in this study.
MethodExpressionPlot
Flynn-Wall-Ozawa (FWO) ln β = ln A E a R G ( x ) 5.311 1.052 E a R T ln β vs. 1/T
Kissinger–Akahira–Sunose (KAS) ln β T 2 = ln A R E a G ( x ) E a R T ln β T 2 vs. 1/T
Table 2. Physical and chemical properties of parent and metal-based H-USY (15) and H-ZSM-5 (11.5) zeolites.
Table 2. Physical and chemical properties of parent and metal-based H-USY (15) and H-ZSM-5 (11.5) zeolites.
CatalystCrystallinity
(%)
Sext
(m2/g)
Vmicro
(cm3/g)
Vmeso
(cm3/g)
Vtotal
(cm3/g)
Weak Acidity
(μmol/g)
Strong Acidity
(μmol/g)
Total Acidity
(μmol/g)
H-USY (15)1001890.250.230.48332 319 651
Ni (5)/H-USY (15)97.41700.240.240.48288 245 533
Pt (0.5)/H-USY (15)93.71840.240.240.48364 256 620
Pd (0.5)/H-USY (15)93.71800.240.230.47324 301 625
HZSM-5 (11.5)1001140.130.100.23740 530 1270
Ni (5)/H-ZSM-5 (11.5)91.2780.120.070.19670 500 1170
Pt (0.5)/H-ZSM-5 (11.5)95.7610.100.050.15650 530 1180
Pd (0.5)/H-ZSM-5 (11.5)97.6770.100.060.17650 500 1150
Table 3. Temperature at which mass loss is 5, 50, and 95% (T5%, T50%, and T95%) calculated from the TGA curves and maximum degradation temperature (Tm) determined from DTG profiles.
Table 3. Temperature at which mass loss is 5, 50, and 95% (T5%, T50%, and T95%) calculated from the TGA curves and maximum degradation temperature (Tm) determined from DTG profiles.
Sampleβ (°C/min)Temperature (°C)
T5%T50%T95%Tm
HDPE5433463473468
10433478488480
15455484494488
20463489501493
HDPE + H-USY (15)5255354398392
10271370410399
15289382418408
20292394431418
HDPE + H-ZSM-5 (11.5)5283362397389
10318407428421
15325413430432
20329421443442
Table 4. Ea, ln(A/G(x)), and R2 in function of the extent of conversion (x) determined by distinct model-free methods.
Table 4. Ea, ln(A/G(x)), and R2 in function of the extent of conversion (x) determined by distinct model-free methods.
xFWOKAS
Ea (kJ/mol)ln(A/G(x))R2Ea (kJ/mol)ln(A/G(x))R2
HDPE0.2212.533.81.000211.233.51.000
0.3220.934.91.000220.034.81.000
0.4226.935.80.999226.235.60.999
0.5229.236.00.999228.635.90.999
0.6232.036.30.999231.536.20.999
0.7232.736.30.999232.336.20.999
0.8233.036.20.999233.036.10.999
0.9234.136.30.999233.736.10.999
Av.227.735.70.999227.135.60.999
HDPE+
H-USY (15)
0.283.816.10.91978.314.40.899
0.393.117.40.97687.715.80.971
0.4104.719.10.99899.517.60.997
0.5115.020.50.988110.219.30.985
0.6126.022.20.979121.621.10.975
0.7137.023.90.970133.022.90.965
0.8148.825.70.967145.324.90.961
0.9172.229.60.988157.523.80.810
Av.122.621.80.973116.620.00.946
HDPE+
H-ZSM-5 (11.5)
0.260.610.90.88953.37.80.849
0.363.511.10.90056.18.60.865
0.469.511.90.90462.29.60.873
0.576.613.00.90969.610.80.882
0.684.314.20.91577.512.20.892
0.791.215.30.91884.713.40.897
0.898.516.40.92692.214.70.908
0.9106.917.70.929122.217.10.994
Av.81.413.80.91177.211.80.895
Table 5. Ea, ln(A/G(x)), and R2 as a function of the conversion (x) determined by the KAS method for metal (Ni, Pt, and Pd)-based H-USY (15) and H-ZSM-5 (11.5) zeolites. * [Av-Average].
Table 5. Ea, ln(A/G(x)), and R2 as a function of the conversion (x) determined by the KAS method for metal (Ni, Pt, and Pd)-based H-USY (15) and H-ZSM-5 (11.5) zeolites. * [Av-Average].
xNi/H-USY (15)Pt/H-USY (15)Pd/H-USY (15)
Ea
(kJ/mol)
ln(A/G(x))Ea
(kJ/mol)
ln(A/G(x))Ea
(kJ/mol)
ln(A/G(x))
0.2153.430.4219.345.1202.243.8
0.3149.728.8182.536.1180.435.4
0.4141.326.6165.331.8161.232.3
0.5134.324.8153.428.9149.729.2
0.6131.423.9144.026.6144.528.5
0.7132.923.9136.624.7132.226.7
0.8137.724.5131.523.4123.525.3
0.9159.628.2130.022.7127.624.5
Av.*142.526.4157.829.9152.630.7
xNi/H-ZSM-5 (11.5)Pt/H-ZSM-5 (11.5)Pd/H-ZSM-5 (11.5)
Ea
(kJ/mol)
ln(A/G(x))Ea
(kJ/mol)
ln(A/G(x))Ea
(kJ/mol)
ln(A/G(x))
0.292.217.870.412.068.315.3
0.395.318.172.212.070.116.1
0.4101.919.077.812.873.416.1
0.5109.220.285.414.077.215.1
0.6117.821.693.115.283.215.4
0.7126.623.099.916.387.216.1
0.8133.424.1107.517.590.517.1
0.9147.326.4115.318.898.218.7
Av.*115.521.390.214.887.016.2
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

Costa, C.S.; Fernandes, A.; Munoz, M.; Ribeiro, M.R.; Silva, J.M. Analyzing HDPE Thermal and Catalytic Degradation in Hydrogen Atmosphere: A Model-Free Approach to the Activation Energy. Catalysts 2024, 14, 514. https://doi.org/10.3390/catal14080514

AMA Style

Costa CS, Fernandes A, Munoz M, Ribeiro MR, Silva JM. Analyzing HDPE Thermal and Catalytic Degradation in Hydrogen Atmosphere: A Model-Free Approach to the Activation Energy. Catalysts. 2024; 14(8):514. https://doi.org/10.3390/catal14080514

Chicago/Turabian Style

Costa, Cátia S., A. Fernandes, Marta Munoz, M. Rosário Ribeiro, and João M. Silva. 2024. "Analyzing HDPE Thermal and Catalytic Degradation in Hydrogen Atmosphere: A Model-Free Approach to the Activation Energy" Catalysts 14, no. 8: 514. https://doi.org/10.3390/catal14080514

APA Style

Costa, C. S., Fernandes, A., Munoz, M., Ribeiro, M. R., & Silva, J. M. (2024). Analyzing HDPE Thermal and Catalytic Degradation in Hydrogen Atmosphere: A Model-Free Approach to the Activation Energy. Catalysts, 14(8), 514. https://doi.org/10.3390/catal14080514

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop