Next Article in Journal
Carbonized Aminal-Linked Porous Organic Polymers Containing Pyrene and Triazine Units for Gas Uptake and Energy Storage
Previous Article in Journal
Phormidium ambiguum and Leptolyngbya ohadii Exopolysaccharides under Low Water Availability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Stability and Non-Isothermal Kinetic Analysis of Ethylene–Propylene–Diene Rubber Composite

1
Department of Physics, College of Science, King Faisal University, P.O. Box 400, Al-Ahsa 31982, Saudi Arabia
2
Department of Physics, Faculty of Science, Cairo University, Giza 12613, Egypt
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(8), 1890; https://doi.org/10.3390/polym15081890
Submission received: 18 February 2023 / Revised: 10 April 2023 / Accepted: 11 April 2023 / Published: 14 April 2023

Abstract

:
The purpose of this study was to investigate the thermal stability and the decomposition kinetics of ethylene–propylene–diene monomer (EPDM) composite samples loaded with and without lead powder (50, 100, and 200 phr lead) using thermogravimetric analysis (TGA). TGA was carried out at different heating rates (5, 10, 20, and 30 °C/min) under inert conditions in the temperature range of 50–650 °C. Lead addition did not significantly change the onset temperature or peak position corresponding to the maximum decomposition rate of the first derivative of the TGA curve (DTGA) (onset at about 455 °C and Tm at about 475 °C). Peak separation for the DTGA curves indicated that the main decomposition region for EPDM, the host rubber, overlapped the main decomposition region for volatile components. The decomposition activation energy (Ea) and pre-exponent factor (A) were estimated using the Friedman (FM), Kissinger–Akahira–Sunose (KAS), and Flynn–Wall–Ozawa (FWO) iso-conversional methods. Average activation energy values of around 231, 230, and 223 kJ/mol were obtained for the EPDM host composite using the FM, FWO, and KAS methods, respectively. For a sample loaded with 100 phr lead, the average activation energy values obtained via the same three methods were 150, 159, and 155 kJ/mole, respectively. The results obtained from the three methods were compared with results obtained using the Kissinger and Augis–Bennett/Boswell methods, and strong convergence was found among the results of the five methods. A significant change in the entropy of the sample was detected with the addition of lead powder. For the KAS method, the change in entropy, ΔS, was −3.7 for EPDM host rubber and −90 for a sample loaded with 100 phr lead, α = 0.5.

1. Introduction

EPDM is one of the most popular rubbers due to its many advantages and properties. This type of rubber is utilized in the preparation of composite materials with additives and ingredients. It is also used for the preparation of blended rubber composites with another type of rubber or polymer such as NBR, PVC, etc., to acquire new characteristics and open wide horizons for technological applications. This type of rubber is famous for its radiation shielding properties, especially against gamma rays and neutrons [1,2,3,4]. This is due to its hydrogen abundance. EPDM composites are also used in thermal insulation, leakage prevention, and other relevant applications. Additionally, composites based on EPDM can be made for electromagnetic shielding using additives with magnetic and electrical properties [5,6].
The thermal stability and decomposition mechanisms for polymeric and rubber composites are among the key properties of materials exposed to high temperatures. In view of the exposure of many products based on EPDM rubber to a heating environment, such as reactors, the automobile industry, etc., it is important to study the effect of any additives on the thermal stability of the composite, as well as to predict the most common degradation mechanism [7,8,9,10,11].
Many studies have examined the effectiveness of blending EPDM rubber with other types of rubber to improve its physical or chemical properties; for instance, blending EPDM with natural rubber has been shown to improve its mechanical properties [12]. Scientists have studied the impact of additives such as accelerators [13], vulcanization agents [14], and other additives (during the preparation process) on the properties of EPDM [15,16,17]. Much research has been conducted on the thermal stability and decomposition mechanisms of EPDM composites or blends. The researchers examined how additives and other treatments affect the composite decomposition mechanisms and thermal stability.
Neelesh Ashok et al. [18] investigated the mechanical, viscoelastic, thermal, and transport properties of EPDM-CIIR blends reinforced with bis(3-triethoxysilylpropyl)tetrasulfide (TESPT)-grafted nano-silica after exposure to different cumulative radiation doses. The nanocomposites had a degradation onset temperature of 503 °C, compared to 485 °C for unfilled blends. The thermal degradation activation energy was evaluated using the Coats–Redfern and Horowitz–Metzger models. In addition, this research group performed research related to the use of EPDM in radiation protection, and they investigated the thermal decomposition of different EPDM composites [19,20,21]. A Güngör et al. [22] investigated the effect of surface-modified waste walnut shells (as a filler) on the mechanical, thermal, and rheological properties of EPDM/waste walnut shell filler composite. The TGA results showed similar T50 temperatures for the CB30 and CB30 + S30 samples. This could be attributed to the walnut shells not significantly affecting thermal stability. The thermal degradation of a blend of NR/EPDM (70:30 wt%) under anaerobic and aerobic conditions was investigated using mercapto- and thioacetate-modified EPDM (EPDMSH) [23]. The kinetic parameters of thermal degradation were evaluated from non-isothermal TGA experiments undertaken at different heating rates. The results indicated that the presence of functionalized polymers (in low proportions) led to an increase in the calculated decomposition activation energy, which indicates an increase in resistance to thermal degradation. Using thermal kinetic analysis, the kinetic parameters and decomposition mechanisms of ethylene propylene copolymers (EPM) and ethylene propylene diene terpolymers (EPDM) were also studied by Antonio Perejón et al. [24].
The thermal degradation of all compositions has been described using the random scission model, and it has been shown that the activation energy for the thermal decomposition of EPM decreases with an increase in the proportion of propylene. The EPDM chemical composition is not strongly correlated with the activation energy for pyrolysis.
The thermal decomposition mechanisms of the host EPDM rubber composite and of the same composite loaded with lead powder have not been studied in a detailed manner, despite the importance of this mixture in many applications [25,26]. Therefore, the focus of this research was to study the effect of adding lead powder at different concentrations to the basic composite and to identify the decomposition mechanisms through the study of thermal stability and kinetic thermal analysis. This research also aimed to compare the different models with each other and with the iso-conversion methods, in order to maximize the knowledge outcomes from this work.

2. Materials and Methods

2.1. Preparation of EPDM Composites

The EPDM host composite and the composites loaded with varying phr of lead were synthesized by using a two-roll mill. Carbon black (N-220) and other ingredients (ZnO, stearic acid, paraffin oil, etc.) were mixed on an open rubber roll mill (ASTM D3182–16) at different concentrations of lead (50, 100, and 200 phr). After mixing, the composites were vulcanized via compression molding and hot pressed at about 153 °C, under 4 MPa, for 5 min. The final product was in the form of sheets of suitable thickness (4 mm) and diameter (10 cm) for examination.
The sample ingredients and components are outlined in Table 1. The values are expressed in phr (parts per hundred).

2.2. Thermogravimetric Analysis

Thermogravimetric analysis (TGA) was performed using a TA Instruments Q50 thermogravimetric analyzer, under a nitrogen (N2) flow of 20 mL/min. The measurements were carried out at different heating rates (5, 10, 20, and 30 °C/min) with constant mass of about 10 mg.

3. Theoretical Approach

3.1. Models for Thermal Kinetic Analysis

The decomposition or conversion rate (alpha) is used in analyzing thermal decomposition mechanisms and can be expressed as follows [27]:
α = w 0 w t w 0 w f
where w0, wt, and wf are the initial, residual (at time t), and final weights of the material, respectively. The rate of conversion (/dt) is directly proportional to the conversion function f(α):
d α d t = k   f α = 1 α n
where d α d t is the rate of conversion in s−1, k is the rate constant, α is the degree of conversion, and n is the order of reaction.
In the Arrhenius equation, the activation energy is quantitatively related to the reaction rate. The equation thus provides a quantitative basis for understanding how the reaction rates depend on the activation energy. The rate constant as a function of temperature is given by
k = A   exp   ( E R T )
where E is the activation energy (J.mo1−1), A is the pre-exponential factor, R is the gas constant (8.31 J.mo1−1.K−1), and T is the absolute temperature (K). Using these equations, one can obtain
d α 1 α n = A β   exp   ( E R T )
Let   g α = 0 α d α 1 α n = A β   T 0 T exp E R T d T
The above integral cannot be solved exactly, but it can be expressed in numerical form or by some approximations. For estimating the kinetic parameters, model-free iso-conversional methods such as Friedman’s (FM), Kissinger–Akahira–Sunose (KAS) [28], and Flynn–Wall–Ozawa (FWO) [29,30] were chosen.
Using the KAS method and according to (6), the relationship between ln[β/T2] and 1/T was plotted for each value of conversion α. From the slope and the intersection, Ea and A can be calculated.
ln β T 2 = ln A R g α E R T
Based on Doyle’s approximation (Doyle 1962), the integral in (5) yields the following Equation (7) via simplification, which forms the basis of further iso-conversional methods developed by Flynn and Wall (Chan and Balke 1997) and Ozawa [29,30].
ln β = ln A E R   g α 5.331 1.052   E R T
Thus, the activation energy and pre-exponential factor can be determined from the slope and intercept of the linear relation ln β vs. 1/T for constant values of α.
According to Friedman (1964) [31], the Friedman method (a differential iso-conversional method) provides correct values of the kinetic parameters:
ln d a d t = ln β d α d T = ln   A . f a E R T
where ln(β da/dT) is plotted against 1/T at constant values of α, and from the slope and intercept, both Ea and A can be determined.

3.2. Thermodynamic Analysis

By applying the following equations with the obtained activation energy, the thermodynamic parameters, including the pre-exponential factor (A), changes in enthalpy (ΔH), Gibbs free energy (ΔG), and entropy (ΔS), can be determined [32,33]:
A = β   E exp E R   T m / R   T m 2
Δ H = E R   T α
Δ G = E + R   T m   ln k B T m A   h
Δ S = Δ H Δ G T m
where kB is the Boltzmann constant, h is the Plank constant, Tm is the temperature of the maximum decomposition rate, and Tα is the temperature corresponding to α.

4. Results and Discussion

4.1. Thermogravimetric Analysis

The TGA and DTGA curves of the EPDM host sample and samples loaded with 50, 100, and 200 phr lead are presented in Figure 1. The composite samples underwent two main decomposition stages.
As illustrated in Figure 2, by using the Origin−Lab program (version 8.6), the main interaction stage (in DTGA) can be resolved to two main decomposition process. The first decomposition process (Region I) starts at about 50 °C and theoretically (according to the peak profile) ends at about 600 °C, while the second process (Region II) starts at about 350 °C and ends at about 580 °C. There are no references indicating that EPDM rubber decomposes at temperatures below 350 °C [34,35,36]. Decomposition in region I can thus be attributed to the volatile components used during the preparation of the host rubber composite, such as paraffin oil. Comparing the peak position of the maximum decomposition rate for region II with the literature confirms that the decomposition in this region is caused by the decomposition of the host rubber, EPDM.
The curves clearly indicate that as the concentration of lead powder in the samples increases, the sample becomes more thermally stable, and the residual mass increases [26]. As a result of the addition of lead, the position of maximum decomposition ( 277.5   to   275   ° C ) did not change significantly, indicating that the influence of lead on thermal behavior is physical rather than chemical.
Table 2 lists some of the parameters that could be extracted from the thermal decomposition curves for the EPDM composite samples.
It can be seen from the above analysis that the introduction of lead powder to the EPDM host composite had an important influence on the thermal stability of the host composite. The results can be summarized as follows:
  • The effect of lead powder on the main decomposition peak (maximum decomposition rate) indicates that there was no chemical interaction between lead and EPDM. The variations in these peak positions are within the experimental error. Through previous research, we note that some additives may not affect position the maximum decomposition rate of EPDM rubber [37,38].
  • Lead powder’s effect is evident between 200 and 400 °C, especially at higher concentrations.
  • The decomposed mass at T > 250 °C decreased with increasing lead concentration. This can be attributed to the following: Lead particles absorb a large amount of thermal energy and, in turn, delay the decomposition of the host composite. During the preparation and vulcanization process, the lead increases the thermal homogeneity of the sample, leading to the formation of a homogenous network of cross-linked chains. There has been some evidence that conductive fillers impact EPDM’s thermal behavior, especially its thermal conductivity [39,40,41].
  • As the lead powder concentration increased, the residual mass increased [42].

4.2. Kinetic Analysis

The decomposition of both the host composite and the composites loaded with lead takes place in two stages: a main one for EPDM rubber [43] and a weak one for the volatile matter and other ingredients. There is an overlap between both decompositions in a range of temperatures (as previously shown in Figure 2). This behavior was not affected by change in the heating rate for both sample types, as shown in Figure 3 and Figure 4. In Figure 3b and Figure 4b, the peak corresponding to the maximum decomposition rate shifted to a higher temperature with the heating rate. Based on these data, the kinetic parameters (Ea and A) were estimated for both the host composite and the composites loaded with lead.
Figure 5 shows the relationship between ln(β/Tp2) and 1000/T, as well as that between ln(β/Tp) and 1000/T, according to the Kissinger [28] and Augis–Bennett/Boswell methods [44,45], respectively, for determination of the activation energy. Both methods gave approximately equal results within the permissible error. The activation energy value calculated using the Kissinger method was about 227 kJ/mol, while the Augis–Bennett/Boswell method gave a value of about 229 kJ/mol. The error was less than 1% between the two methods. These results are largely consistent with the references [46,47].
By applying the same two methods (Kissinger and Augis–Bennett), the activation energy was calculated for the sample loaded with 100 phr lead powder (Figure 6). The activation energies calculated by the Kissinger and Augis–Bennett methods were 141 and 144 kJ/mol, respectively. The addition of lead produced a significant decrease in the activation energy of decomposition, and the two methods gave values that are very close to each other. The decrease of the activation energy with the addition of lead can be attributed to the increase in the thermal conductivity of the composite sample and the decrease of entanglements in the composite matrix.
The Kissinger and Augis–Bennett methods ignore the changes in conversion rates (/dt) and conversion (α). They assume that Ea is independent of α. Three methods that are derived from (5) and account for the conversion function f(α) will be illustrated. These methods are Kissinger–Akahira–Sunose, Flynn–Wall–Ozawa, and Friedman ((6), (7), and (8), respectively).
Figure 7 illustrates the dependence of the conversion α (Figure 7a) and its first derivative (Figure 7b) on the temperature. Using these results, the main decomposition stage of the host EPDM with and without lead was analyzed, and a comparison among the three methods was carried out.
Figure 8, Figure 9 and Figure 10 illustrate the application of the KAS, FWO, and Friedman methods, respectively, for EPDM host composite. A resemblance was noted in the behavior of the activation energy and coefficient A for both the KAS (Figure 8c,d) and FWO (Figure 9c,d) methods. The behavior was different when applying Friedman’s method (Figure 10c,d). This may be attributed to the dependence of Friedman’s method (8) on the first derivative with respect to time of the degree of conversion (dα/dt). Despite this apparent difference in behavior, there are some common points that could be concluded. The total average values of activation energy for the three methods are distinctly close (229, 229.6, and 231 kJ/mole for the KAS, FWO, and Friedman methods, respectively). Additionally, the pre-exponential factors for the three methods are of the same order of magnitude (10+15) (see Figure 8, Figure 9 and Figure 10d).
In accordance with the E(α) relation for the three methods mentioned above, there is a region where the activation energy dependence on the degree of conversion (α) is weak or not significant and can be considered constant after looking at the value of its standard deviation (0.64, 0.53, and 3.31, respectively). The average activation energy for each method is close to those of the Kissinger and Augis–Bennett methods, with a slight difference not exceeding 5%, as shown in Table 3.
Figure 11, Figure 12 and Figure 13 show the application of the KAS, FWO, and Friedman methods, respectively, for EPDM host loaded with 100 phr lead powder. It was observed that the activation energy values are close to each other at alpha greater than 0.4, but the activation energy changes when alpha is less than 0.4. There is no specific behavior for the three samples within this range (α < 0.4).
Table 3 summarizes the average values of activation energy over the entire range of degree of conversion, as well as the mean value of activation energy in the less dependent region (from 0.4 to 0.9) and the standard deviation for each data set.
Table 4, Table 5 and Table 6 illustrate the values of the activation energy and the pre-exponent factor, respectively (for the host composite and that loaded with 100 phr lead). The deviation in the activation energy values between the KAS and FWO methods does not exceed 4%; when compared with the Friedman method, however, the difference may reach more than 10% at some values of α. This difference may be due to the dependence of the Friedman method on the time derivative of α. The pre-exponential factor, A, varies between 1015 and 1016 for the three methods. The behavior of A is similar for the KAS and FWO methods, while it differs in the Friedman method. The fluctuation in the predicted values is larger in the case of Friedman method, as shown in Table 5.
For the three methods (KAS, FWO, and FM), (9)–(12) were used to evaluate the pre-exponential factor and thermodynamic parameters, namely, the enthalpy change (ΔH), Gibbs free energy change (ΔG), and entropy change (ΔS), based on the values of activation energy determined via the iso-conversion free-model methods [48]. Their values are listed in Table 7, Table 8 and Table 9. It is evident from these three tables that the activation energy and pre-exponential factors do not depend significantly on the conversion [49], except for small values; this indicates that the material decomposed simply and not in a complex manner. It is noted that the entropy of the rubber composite sample is much higher than the entropy of the host rubber sample (for the KAS method, ΔS was −3.7 for EPDM host rubber and −90 for a sample loaded with 100 phr lead, with α = 0.5). This result is logical, as the randomness of the system is increased by the introduction of lead particles into the host rubber system [50]. The Gibbs free energy for the EPDM composite sample is still close to the activation energy (for each alpha). This can be attributed to the fact that the nature of the rubber is still the same, and it decomposes in a simple and uncomplicated process.
The values of the pre-exponential factor (A) in Table 6 differ from those values calculated in Table 7, Table 8 and Table 9. This can be attributed to the fact that the former values depend on the value of the activation energy (calculated from the slope of the straight line) and the intercept for each value of α (according to the model or methods (6)–(8)), while the latter values (in Table 7, Table 8 and Table 9) were calculated using the general model (9) that depends on the heating rate and temperature at the maximum decomposition (Tm) of the DTGA curve, besides other constants.
Table 10 summarizes the maximum, minimum, and average values of the basic thermodynamic parameters (enthalpy change, Gibb’s free energy change, and entropy change), activation energies (E), and pre-exponent factors (A).
Figure 14 shows great convergence and agreement among the three methods (KAS, FWO, and FM), as well with the two other methods (the Kissinger and Augis–Bennett/Boswell methods (Figure 5)). Based on the KAS method, for example, Figure 14d shows the scope of convergence between the activation energy and Gibbs free energy change in the host sample. There is a high degree of harmony in this reaction, suggesting that it is an uncomplicated process.

5. Conclusions

The thermal stability of ethylene–propylene–diene monomer (EPDM) composite samples loaded with and without lead powder (50, 100, and 200 phr lead) was examined via thermogravimetric analysis (TGA) at a 10 C/min heating rate. Lead addition did not significantly change the onset temperature or peak position corresponding to the maximum decomposition rate of DTGA (onset at about 455 °C and Tm at about 475 °C). Through peak separation of the DTGA curves, it was confirmed that some degradation processes overlapped with the EPDM decomposition peak. Kinetic thermal analysis was conducted for the host sample and a sample loaded with 100 phr lead. The TGA measurements were carried out at different heating rates (5, 10, 20, and 30 °C/min). Average activation energy values of around 231, 230, and 223 kJ/mol were obtained for the EPDM host composite using the Friedman (FM), Flynn–Wall–Ozawa (FWO), and Kissinger–Akahira–Sunose (KAS) methods, respectively. Using 100 phr lead as a sample, the average activation energy values obtained via the three methods were 150, 159, and 155 kJ/mole, respectively. There was strong convergence between the results obtained via these three methods and the results obtained via the Kissinger and Augis–Bennett/Boswell methods. A significant change in the entropy of the sample was detected with the addition of lead powder.

Author Contributions

Conceptualization, S.S.I. and H.A.; methodology, S.S.I.; Data analysis, S.S.I. and N.A.; investigation, S.S.I. and H.A.; Writing—original draft preparation, N.A., and S.S.I.; writing—review and editing, S.S.I.; supervision, N.A.; project administration, H.A.; funding acquisition, H.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia (Project number INST164).

Data Availability Statement

Detailed data supporting the findings of this study are included in the article, as the authors confirm. In any case, the original collected data are available on request from the corresponding author.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia, for funding this research work through the Project number INST164.

Conflicts of Interest

The co-authors have all approved the contents of this manuscript, and the authors declare that they do not have any conflicts of interest. There is no conflict of interest or personal relationship known to the authors that could have appeared to influence this work.

References

  1. Abdulrahman, S.T.; Patanair, B.; Vasukuttan, V.P.; Thomas, S.; Cadel, E.; Cuvilly, F.; Saiter-Fourcin, A.; Ahmad, Z.; Jaroszewski, M.; Strankowski, M. High-density polyethylene/EPDM rubber blend composites of boron compounds for neutron shielding application. Express Polym. Lett. 2022, 16, 558–572. [Google Scholar] [CrossRef]
  2. Barala, S.S.; Manda, V.; Jodha, A.S.; Meghwal, L.R.; Gopalani, D. Ethylene-propylene diene monomer-based polymer composite for attenuation of high energy radiations. J. Appl. Polym. Sci. 2021, 138, 50334. [Google Scholar] [CrossRef]
  3. Chang, B.; Li, S.; Li, M.; Chen, R. Research Progress of Neutron Shielding Materials. In Proceedings of the 29th International Conference on Nuclear Engineering, Shenzhen, China, 8–12 August 2022; p. V002T002A050. [Google Scholar]
  4. Barala, S.S.; Manda, V.; Jodha, A.S.; Gopalani, D. Thermal stability of gamma irradiated ethylene propylene diene monomer composites for shielding applications. J. Appl. Polym. Sci. 2022, 139, e52975. [Google Scholar] [CrossRef]
  5. Lu, S.; Li, B.; Ma, K.; Wang, S.; Liu, X.; Ma, Z.; Lin, L.; Zhou, G.; Zhang, D. Flexible MXene/EPDM rubber with excellent thermal conductivity and electromagnetic interference performance. Appl. Phys. A 2020, 126, 513. [Google Scholar] [CrossRef]
  6. Li, C.; Guo, J.; Xu, P.; Hu, W.; Lv, J.; Shi, B.; Zhang, Z.; Li, R. Facile preparation of superior compressibility and hydrophobic reduced graphene oxide@ cellulose nanocrystals/EPDM composites for highly efficient oil/organic solvent adsorption and enhanced electromagnetic interference shielding. Sep. Purif. Technol. 2023, 307, 122775. [Google Scholar] [CrossRef]
  7. Stelescu, M.D.; Airinei, A.; Manaila, E.; Fifere, N.; Craciun, G.; Varganici, C.; Doroftei, F. Exploring the effect of electron beam irradiation on the properties of some EPDM-flax fiber composites. Polym. Compos. 2019, 40, 315–327. [Google Scholar] [CrossRef] [Green Version]
  8. Güngör, A.; Akbay, I.; Özdemir, T. Waste walnut shell as an alternative bio-based filler for the EPDM: Mechanical, thermal, and kinetic studies. J. Mater. Cycles Waste Manag. 2019, 21, 145–155. [Google Scholar] [CrossRef]
  9. Güngör, A.; Akbay, İ.; Yaşar, D.; Özdemir, T. Flexible X/Gamma ray shielding composite material of EPDM rubber with bismuth trioxide: Mechanical, thermal investigations and attenuation tests. Prog. Nucl. Energy 2018, 106, 262–269. [Google Scholar] [CrossRef]
  10. Chen, G.; Gupta, A.; Mekonnen, T.H. Silane-modified wood fiber filled EPDM bio-composites with improved thermomechanical properties. Compos. Part A Appl. Sci. Manuf. 2022, 159, 107029. [Google Scholar] [CrossRef]
  11. Yao, W.; Xu, X.; Zhou, J.; Chen, M.; Meng, Q.; Wang, P.; Zhang, B.; Zhou, Y.; Ding, Y. Mechanically robust and flame-retarded EPDM composites with high loading of Mg(OH)2 based on reversible crosslinking network from Diels-Alder reactions. Polym. Degrad. Stab. 2022, 202, 110029. [Google Scholar] [CrossRef]
  12. Jung, J.K.; Lee, C.H.; Baek, U.B.; Choi, M.C.; Bae, J.W. Filler Influence on H2 Permeation Properties in Sulfur-CrossLinked Ethylene Propylene Diene Monomer Polymers Blended with Different Concentrations of Carbon Black and Silica Fillers. Polymers 2022, 14, 592. [Google Scholar] [CrossRef] [PubMed]
  13. Mayasari, H.E.; Wirapraja, A.Y. The Effect of Accelerator on Mechanical and Swelling Properties of NBR/EPDM/BIIR Composite. J. Teknol. Proses Dan Inov. Ind. 2021, 6, 7–11. [Google Scholar]
  14. Rahman Parathodika, A.; Raju, A.T.; Das, M.; Bhattacharya, A.B.; Neethirajan, J.; Naskar, K. Exploring hybrid vulcanization system in high-molecular weight EPDM rubber composites: A statistical approach. J. Appl. Polym. Sci. 2022, 139, e52721. [Google Scholar] [CrossRef]
  15. Rizwan, M.; Chandan, M.R. Mechanistic insights into the ageing of EPDM micro/hybrid composites for high voltage insulation application. Polym. Degrad. Stab. 2022, 204, 110114. [Google Scholar] [CrossRef]
  16. Lee, S.-H.; Park, S.-Y.; Chung, K.-H.; Jang, K.-S. Phlogopite-reinforced natural rubber (NR)/ethylene-propylene-diene monomer rubber (EPDM) composites with aminosilane compatibilizer. Polymers 2021, 13, 2318. [Google Scholar] [CrossRef]
  17. Zhang, X.; Li, J.; Chen, Z.; Pang, C.; He, S.; Lin, J. Study on thermal-oxidative aging properties of ethylene-propylene-diene monomer composites filled with silica and carbon nanotubes. Polymers 2022, 14, 1205. [Google Scholar] [CrossRef]
  18. Ashok, N.; Prakash, K.; Selvakumar, D.; Balachandran, M. Synergistic enhancement of mechanical, viscoelastic, transport, thermal, and radiation aging characteristics through chemically bonded interface in nanosilica reinforced EPDM-CIIR blends. J. Appl. Polym. Sci. 2021, 138, 50082. [Google Scholar] [CrossRef]
  19. Ashok, N.; Balachandran, M.; Lawrence, F.; Sebastian, N. EPDM–chlorobutyl rubber blends in γ-radiation and hydrocarbon environment: Mechanical, transport, and ageing behavior. J. Appl. Polym. Sci. 2017, 134, 45195. [Google Scholar] [CrossRef]
  20. Ashok, N.; Balachandran, M.; Lawrence, F. Organo-modified layered silicate nanocomposites of EPDM–chlorobutyl rubber blends for enhanced performance in γ radiation and hydrocarbon environment. J. Compos. Mater. 2018, 52, 3219–3231. [Google Scholar] [CrossRef]
  21. Ashok, N.; Webert, D.; Suneesh, P.; Balachandran, M. Mechanical and sorption behaviour of organo-modified montmorillonite nanocomposites based on EPDM–NBR Blends. Mater. Today Proc. 2018, 5, 16132–16140. [Google Scholar] [CrossRef]
  22. Güngör, A.; Akbay, I.; Özdemir, T. EPDM rubber with hexagonal boron nitride: A thermal neutron shielding composite. Radiat. Phys. Chem. 2019, 165, 108391. [Google Scholar] [CrossRef]
  23. Da Silva Sirqueira, A.; Guenther Soares, B. Kinetic analysis of thermal degradation of NR/EPDM blends: Effect of the reactive compatibilization. J. Appl. Polym. Sci. 2007, 103, 2669–2675. [Google Scholar] [CrossRef]
  24. Perejón, A.; Sánchez-Jiménez, P.E.; Gil-González, E.; Pérez-Maqueda, L.A.; Criado, J.M. Pyrolysis kinetics of ethylene–propylene (EPM) and ethylene–propylene–diene (EPDM). Polym. Degrad. Stab. 2013, 98, 1571–1577. [Google Scholar] [CrossRef] [Green Version]
  25. Kamath, V.A.; Kumaraswamy, S.; Kiran, K.U.; Somashekarappa, H.M. Optimization of PbO/Fe3O4/EPDM flexible composites for gamma shielding applications. In Proceedings of the International Conference on Physics of Materials and Nanotechnology 2019, Mangalore, India, 19–21 September 2019; p. 040004. [Google Scholar]
  26. El-Khatib, A.M.; Doma, A.; Badawi, M.S.; Abu-Rayan, A.E.; Aly, N.; Alzahrani, J.S.; Abbas, M.I. Conductive natural and waste rubbers composites-loaded with lead powder as environmental flexible gamma radiation shielding material. Mater. Res. Express 2020, 7, 105309. [Google Scholar] [CrossRef]
  27. Brown, M.E. Introduction to Thermal Analysis: Techniques and Applications; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  28. Kissinger, H.E. Reaction kinetics in differential thermal analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  29. Ozawa, T. A new method of analyzing thermogravimetric data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef] [Green Version]
  30. Ozawa, T. Kinetic analysis of derivative curves in thermal analysis. J. Therm. Anal. 1970, 2, 301–324. [Google Scholar] [CrossRef]
  31. Friedman, H.L. Kinetics of thermal degradation of char-forming plastics from thermogravimetry. Application to a phenolic plastic. In Journal of Polymer Science Part C: Polymer Symposia; Wiley Subscription Services, Inc., A Wiley Company: New York, NY, USA, 1964; pp. 183–195. [Google Scholar]
  32. Kim, Y.S.; Kim, Y.S.; Kim, S.H. Investigation of Thermodynamic Parameters in the Thermal Decomposition of Plastic Waste−Waste Lube Oil Compounds. Environ. Sci. Technol. 2010, 44, 5313–5317. [Google Scholar] [CrossRef]
  33. Singh, B.; Singh, S.; Kumar, P. In-depth analyses of kinetics, thermodynamics and solid reaction mechanism for pyrolysis of hazardous petroleum sludge based on isoconversional models for its energy potential. Process Saf. Environ. Prot. 2021, 146, 85–94. [Google Scholar] [CrossRef]
  34. Gamlin, C.; Dutta, N.; Roy Choudhury, N.; Kehoe, D.; Matisons, J. Influence of ethylene–propylene ratio on the thermal degradation behaviour of EPDM elastomers. Thermochim. Acta 2001, 367–368, 185–193. [Google Scholar] [CrossRef]
  35. Praveen, K.M.; Taleb, K.; Pillin, I.; Kervoelen, A.; Grohens, Y.; Thomas, S.; Haponiuk, J.T. Comparative mechanical, morphological, rheological, and thermal properties of polypropylene/ethylene-propylene-diene rubber blends. Polym. Adv. Technol. 2022, 33, 3296–3311. [Google Scholar] [CrossRef]
  36. Kim, G.; Sahu, P.; Oh, J.S. Physical Properties of Slide-Ring Material Reinforced Ethylene Propylene Diene Rubber Composites. Polymers 2022, 14, 2121. [Google Scholar] [CrossRef] [PubMed]
  37. Bianchi, M.; Valentini, F.; Fredi, G.; Dorigato, A.; Pegoretti, A. Thermo-Mechanical Behavior of Novel EPDM Foams Containing a Phase Change Material for Thermal Energy Storage Applications. Polymers 2022, 14, 4058. [Google Scholar] [CrossRef] [PubMed]
  38. Gschwind, L.; Jordan, C.S.; Vennemann, N.; Susoff, M.L. Mechanochemical devulcanization of EPDM rubber waste. Correlation of process parameters with sol–gel analyses and revulcanization properties. J. Appl. Polym. Sci. 2023, 140, e53768. [Google Scholar] [CrossRef]
  39. Ma, A.; Wang, X.; Chen, Y.; Yu, J.; Zheng, W.; Zhao, Y. Largely enhanced thermal conductivity of ethylene-propylene-diene monomer composites by addition of graphene ball. Compos. Commun. 2019, 13, 119–124. [Google Scholar] [CrossRef]
  40. Chen, Y.; Gao, S.; Liu, C.; Situ, Y.; Liu, J.; Huang, H. Preparation of PE-EPDM based phase change materials with great mechanical property, thermal conductivity and photo-thermal performance. Sol. Energy Mater. Sol. Cells 2019, 200, 109988. [Google Scholar] [CrossRef]
  41. Zarrinjooy Alvar, M.; Abdeali, G.; Bahramian, A.R. Influence of graphite nano powder on ethylene propylene diene monomer/paraffin wax phase change material composite: Shape stability and thermal applications. J. Energy Storage 2022, 52, 105065. [Google Scholar] [CrossRef]
  42. Brito, Z.; Sanchez, G. Influence of metallic fillers on the thermal and mechanical behaviour in composites of epoxy matrix. Compos. Struct. 2000, 48, 79–81. [Google Scholar] [CrossRef]
  43. Guo, M.; Li, J.; Xi, K.; Liu, Y.; Ji, J. Effect of multi-walled carbon nanotubes on thermal stability and ablation properties of EPDM insulation materials for solid rocket motors. Acta Astronaut. 2019, 159, 508–516. [Google Scholar] [CrossRef]
  44. Boswell, P. On the calculation of activation energies using a modified Kissinger method. J. Therm. Anal. Calorim. 1980, 18, 353–358. [Google Scholar] [CrossRef]
  45. Zhang, X. Applications of kinetic methods in thermal analysis: A review. Eng. Sci. 2020, 14, 1–13. [Google Scholar] [CrossRef]
  46. Arcenegui-Troya, J.; Sánchez-Jiménez, P.E.; Perejón, A.; Pérez-Maqueda, L.A. Determination of the activation energy under isothermal conditions: Revisited. J. Therm. Anal. Calorim. 2022, 148, 1679–1686. [Google Scholar] [CrossRef]
  47. Zaharescu, T.; Catauro, M.; Blanco, I. Kinetics of the Non-isothermal Degradation of POSS/EPDM Hybrids Composites. Macromol. Symp. 2021, 395, 2000221. [Google Scholar] [CrossRef]
  48. Açıkalın, K. Determination of kinetic triplet, thermal degradation behaviour and thermodynamic properties for pyrolysis of a lignocellulosic biomass. Bioresour. Technol. 2021, 337, 125438. [Google Scholar] [CrossRef]
  49. Marinović-Cincović, M.; Janković, B.; Jovanović, V.; Samaržija-Jovanović, S.; Marković, G. The kinetic and thermodynamic analyses of non-isothermal degradation process of acrylonitrile–butadiene and ethylene–propylene–diene rubbers. Compos. Part B Eng. 2013, 45, 321–332. [Google Scholar] [CrossRef]
  50. Prasopdee, T.; Smitthipong, W. Effect of Fillers on the Recovery of Rubber Foam: From Theory to Applications. Polymers 2020, 12, 2745. [Google Scholar] [CrossRef]
Figure 1. The TGA and DTGA curves of EPDM host composite (0 phr lead) and samples loaded with 50, 100, 200 phr lead. Scanning rate, 10 °C min−1. (a) TGA curves; (b) DTGA curves.
Figure 1. The TGA and DTGA curves of EPDM host composite (0 phr lead) and samples loaded with 50, 100, 200 phr lead. Scanning rate, 10 °C min−1. (a) TGA curves; (b) DTGA curves.
Polymers 15 01890 g001
Figure 2. Peak separation of the DTGA curves for EPDM host samples (a) without lead and (b) with 100 phr lead.
Figure 2. Peak separation of the DTGA curves for EPDM host samples (a) without lead and (b) with 100 phr lead.
Polymers 15 01890 g002
Figure 3. (a) TGA and (b) DTGA curves of EPDM host composite at different heating rates of 5, 10, 20, and 30 °C/min.
Figure 3. (a) TGA and (b) DTGA curves of EPDM host composite at different heating rates of 5, 10, 20, and 30 °C/min.
Polymers 15 01890 g003
Figure 4. (a) TGA and (b) DTGA curves of EPDM host composite loaded with 100 phr lead powder at different heating rates of 5, 10, 20, and 30 °C/min.
Figure 4. (a) TGA and (b) DTGA curves of EPDM host composite loaded with 100 phr lead powder at different heating rates of 5, 10, 20, and 30 °C/min.
Polymers 15 01890 g004
Figure 5. Determination of the activation energy (Ea) in accordance with (a) the Kissinger method and (b) the Augis–Bennett/Boswell method.
Figure 5. Determination of the activation energy (Ea) in accordance with (a) the Kissinger method and (b) the Augis–Bennett/Boswell method.
Polymers 15 01890 g005
Figure 6. Determination of the activation energy (Ea) for EPDM/100 phr lead composite in accordance with (a) the Kissinger method and (b) the Augis–Bennett/Boswell method.
Figure 6. Determination of the activation energy (Ea) for EPDM/100 phr lead composite in accordance with (a) the Kissinger method and (b) the Augis–Bennett/Boswell method.
Polymers 15 01890 g006
Figure 7. Temperature dependence of (a) the conversion, α, and (b) the conversion rate, da/dt, for the EPDM host composite sample.
Figure 7. Temperature dependence of (a) the conversion, α, and (b) the conversion rate, da/dt, for the EPDM host composite sample.
Polymers 15 01890 g007aPolymers 15 01890 g007b
Figure 8. The Kissinger–Akahira–Sunose method for the EPDM host composite: (a,b) ln(β/T2) versus temperature; (c,d) the dependence of Ea and A, respectively, on conversion α.
Figure 8. The Kissinger–Akahira–Sunose method for the EPDM host composite: (a,b) ln(β/T2) versus temperature; (c,d) the dependence of Ea and A, respectively, on conversion α.
Polymers 15 01890 g008aPolymers 15 01890 g008b
Figure 9. The Flynn−Wall−Ozawa (FWO) method for the EPDM host composite: (a,b) ln(β/T2) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on conversion α.
Figure 9. The Flynn−Wall−Ozawa (FWO) method for the EPDM host composite: (a,b) ln(β/T2) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on conversion α.
Polymers 15 01890 g009
Figure 10. The Friedman method for EPDM host composite: (a,b) ln(βdα/dT) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on α.
Figure 10. The Friedman method for EPDM host composite: (a,b) ln(βdα/dT) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on α.
Polymers 15 01890 g010
Figure 11. The Kissinger−Akahira−Sunose method for EPDM-100 phr lead composite: (a,b) ln(β/T2) versus temperature; (c,d) the dependence of Ea and A, respectively, on conversion α.
Figure 11. The Kissinger−Akahira−Sunose method for EPDM-100 phr lead composite: (a,b) ln(β/T2) versus temperature; (c,d) the dependence of Ea and A, respectively, on conversion α.
Polymers 15 01890 g011
Figure 12. The Flynn–Wall−Ozawa (FWO) method for EPDM−100 phr lead composite: (a,b) ln(β/T2) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on conversion α.
Figure 12. The Flynn–Wall−Ozawa (FWO) method for EPDM−100 phr lead composite: (a,b) ln(β/T2) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on conversion α.
Polymers 15 01890 g012aPolymers 15 01890 g012b
Figure 13. The Friedman method for EPDM−100 phr lead composite: (a,b) ln(βdα/dT) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on conversion α.
Figure 13. The Friedman method for EPDM−100 phr lead composite: (a,b) ln(βdα/dT) versus 1000/T; (c,d) the dependence of Ea and A, respectively, on conversion α.
Polymers 15 01890 g013
Figure 14. A comparison among the activation energy values of the host and lead-loaded EPDM rubber samples using the three methods (a) KAS, (b) FWO, and (c) FM, and (d) a comparison between the values of the activation energy, Ea, and the change in the Gibbs free energy, ΔG (for KAS, using the host rubber sample as an example).
Figure 14. A comparison among the activation energy values of the host and lead-loaded EPDM rubber samples using the three methods (a) KAS, (b) FWO, and (c) FM, and (d) a comparison between the values of the activation energy, Ea, and the change in the Gibbs free energy, ΔG (for KAS, using the host rubber sample as an example).
Polymers 15 01890 g014
Table 1. The chemical formulations and ingredients used for sample preparation (in phr).
Table 1. The chemical formulations and ingredients used for sample preparation (in phr).
MaterialRoleAmount (phr)
EPDMRubber base100 phr
Carbon Black
(N–220-(Degussa))
Reinforcement (filler)50 phr
Paraffin Oil (ZMTH-Egypt) Softener (plasticizer)50 phr
Stearic acid (Behn Mayer)Activator2 phr
ZnO (Zink Egypt)Accelerator5 phr
ZDMC * (Behn Mayer)Accelerator3 phr
Sulfur (Polychem-egypt)Curing agent3 phr
Lead powder (Pb-Indiamart) additive
reinforcing   filler
Variable
(50, 100, 200)
* ZDMC (Zinc dimethyldithiocarbamate).
Table 2. Thermogravimetric analysis results of EPDM composites.
Table 2. Thermogravimetric analysis results of EPDM composites.
SampleTonset1
°C
Tonset2
°C
Tmax2
°C
M65%m@450 °C
(%)
m@550 °C
(%)
Mres
EPDM0L239.4455.8477.5458.069.927.024.6
EPDM50L237.4455.2475.346172.339.038.6
EPDM100L235.4453.7475.0467.876.547.147.0
EPDM200L252.8454.8475.9479.782.560.560.3
Table 3. The average activation energy values for the EPDM host composite.
Table 3. The average activation energy values for the EPDM host composite.
Eavr (α from 0.05 to 0.95)Eavr (α from 0.2 to 0.8)
Eavr (kJ/mol)STDEavr (kJ/mol)STD
Kissinger226.7-226.7-
Augis–Bennett/Boswell228.7-228.7-
KAS229.25.26231.00.64
FWO229.65.15231.40.55
Friedman231.08.20227.33.31
Table 4. The average activation energy values for EPDM host composite loaded with 100 phr lead.
Table 4. The average activation energy values for EPDM host composite loaded with 100 phr lead.
Eavr (α from 0.05 to 0.95)Eavr (α from 0.4 to 0.95)
Eavr (kJ/mol)STDEavr (kJ/mol)STD
Kissinger141.3-141.3-
Augis–Bennett/Boswell144.0-144.0-
KAS155.03.6152.70.5
FWO159.03.3156.80.7
Friedman149.24.5151.62.8
Table 5. The activation energy values calculated using three methods for the host composite (0 phr lead) and composite loaded with 100 phr lead.
Table 5. The activation energy values calculated using three methods for the host composite (0 phr lead) and composite loaded with 100 phr lead.
αEPDM HostEPDM/100 phr Lead
AlphaKASFWOFreidmanKASFWOFreidman
0.05209.07210.00245.46161.54164.46147.10
0.10222.91223.29240.06161.46165.98150.00
0.15228.12228.34237.75161.27164.16149.08
0.20231.00231.15235.96160.56163.57145.56
0.25231.71231.88231.38158.89161.64142.11
0.30231.29231.52229.67156.50160.13139.29
0.35231.30231.57226.89154.53158.85142.74
0.40231.56231.86227.00154.18157.87146.55
0.45231.14231.49226.23152.99157.15149.24
0.50231.72232.08225.46152.91156.97149.86
0.55231.33231.73225.24152.51156.87149.97
0.60231.21231.65225.28152.51156.09151.14
0.65230.49230.99224.04152.84156.80150.56
0.70229.87230.43224.87152.36156.02151.65
0.75230.06230.63225.21152.85156.55152.28
0.80230.08230.67227.12152.49156.41152.28
0.85229.86230.48226.39152.31156.45152.81
0.90230.89231.50230.01152.90157.23155.64
0.95230.98231.62254.20152.88158.24157.10
Table 6. The pre-exponential factor calculated using three methods for the host composite (0 phr lead) and composite loaded with 100 phr lead.
Table 6. The pre-exponential factor calculated using three methods for the host composite (0 phr lead) and composite loaded with 100 phr lead.
αEPDM HostEPDM/100 phr Lead
AlphaKASFWOFreidmanKASFWOFreidman
0.051.20 × 10151.09 × 10156.73 × 10164.52 × 10119.53 × 10113.43 × 109
0.108.14 × 10156.64 × 10152.54 × 10162.58 × 10116.89 × 10116.74 × 109
0.151.35 × 10161.07 × 10161.77 × 10161.70 × 10113.37 × 10116.69 × 109
0.201.62 × 10161.27 × 10161.39 × 10161.08 × 10112.17 × 10113.16 × 109
0.251.27 × 10161.10 × 10167.05 × 10156.08 × 10101.17 × 10111.87 × 109
0.301.04 × 10168.22 × 10155.80 × 10153.16 × 10107.06 × 10101.40 × 109
0.358.40 × 10156.72 × 10154.00 × 10151.85 × 10104.60 × 10103.08 × 109
0.407.14 × 10155.73 × 10153.77 × 10151.45 × 10103.24 × 10107.06 × 109
0.455.48 × 10154.44 × 10154.03 × 10151.00 × 10102.43 × 10101.31 × 1010
0.504.97 × 10154.02 × 10154.77 × 10158.39 × 1091.99 × 10101.67 × 1010
0.554.66 × 10153.14 × 10154.27 × 10156.62 × 1091.65 × 10101.76 × 1010
0.604.37 × 10152.51 × 10154.58 × 10155.55 × 1091.22 × 10102.56 × 1010
0.654.06 × 10151.83 × 10154.05 × 10154.86 × 1091.13 × 10102.42 × 1010
0.704.12 × 10151.32 × 10155.53 × 10153.63 × 1098.05 × 1093.46 × 1010
0.753.43 × 10151.06 × 10156.65 × 10153.08 × 1096.85 × 1094.19 × 1010
0.803.29 × 10157.92 × 10141.03 × 10162.18 × 1091.18E × 1064.17 × 1010
0.852.73 × 10155.35 × 10141.07 × 10161.46 × 1093.48 × 1094.75 × 1010
0.902.74 × 10153.85 × 10142.25 × 10169.71 × 1082.37 × 1096.27 × 1010
0.952.30 × 10151.77 × 10142.29 × 10164.12 × 1081.18 × 1097.46 × 1010
Table 7. The activation energy, pre-exponential factor, and thermodynamic parameters obtained using the KAS method at a 10 °C/min heating rate for host and EPDM rubber loaded with 100 phr lead.
Table 7. The activation energy, pre-exponential factor, and thermodynamic parameters obtained using the KAS method at a 10 °C/min heating rate for host and EPDM rubber loaded with 100 phr lead.
EPDM 100 Lead at Rate 10 °C/minEPDM 0 phr Lead at Rate 10 °C/min
αEAΔHΔGΔSEAΔHΔGΔS
(kJ/mol)(sec−1)(kJ/mol)(kJ/mol)(J/mol)(kJ/mol)(sec−1)(kJ/mol)(kJ/mol)(J/mol)
0.052061.47 × 1012199.7220.4−27.61623.51 × 109155.8211.6−77.5
0.12171.01 × 1013211.3220.0−11.61613.13 × 109155.1211.6−78.6
0.152212.01 × 1013215.5219.9−5.91613.20 × 109155.1211.6−78.5
0.22242.84 × 1013217.5219.9−3.11602.60 × 109153.9211.6−80.4
0.252253.59 × 1013218.9219.8−1.21592.10 × 109152.6211.7−82.2
0.32253.49 × 1013218.7219.8−1.51561.33 × 109149.9211.8−86.0
0.352253.33 × 1013218.4219.8−1.91551.01 × 109148.3211.8−88.4
0.42243.12 × 1013218.0219.8−2.51549.41 × 108147.9211.8−89.0
0.452242.82 × 1013217.4219.9−3.31537.72 × 108146.7211.9−90.7
0.52232.72 × 1013217.1219.9−3.71537.60 × 108146.6211.9−90.8
0.552243.23 × 1013218.2219.8−2.21537.62 × 108146.6211.9−90.8
0.62232.43 × 1013216.4219.9−4.61537.15 × 108146.2211.9−91.4
0.652242.80 × 1013217.3219.9−3.41537.20 × 108146.3211.9−91.3
0.72232.61 × 1013216.9219.9−4.01526.85 × 108146.0211.9−91.8
0.752222.28 × 1013216.0219.9−5.21537.50 × 108146.5211.9−91.0
0.82232.68 × 1013217.0219.9−3.81527.00 × 108146.1211.9−91.6
0.852232.45 × 1013216.4219.9−4.61526.77 × 108145.8211.9−91.9
0.92222.26 × 1013215.9219.9−5.31526.94 × 108146.0211.9−91.8
0.952232.73 × 1013217.0219.9−3.81537.22 × 108155.8211.6−77.5
Table 8. The activation energy, pre-exponential factor, and thermodynamic parameters obtained using the FWO method at a 10 °C/min heating rate for host and 100 phr loaded EPDM rubber.
Table 8. The activation energy, pre-exponential factor, and thermodynamic parameters obtained using the FWO method at a 10 °C/min heating rate for host and 100 phr loaded EPDM rubber.
FWOEPDM 0 phr Lead at Rate 10 °C/minEPDM 100 phr Lead at Rate 10 °C/min
αEAΔHΔGΔSEAΔHΔGΔS
(kj/mol)(sec−1)(kJ/mol)(kJ/mol)(J/mol)(kj/mol)(sec−1)(kJ/mol)(kJ/mol)(j/mol)
0.05210.03.0 × 1012204.1220.3−21.5164.55.74 × 109158.7211.4−73.5
0.1223.32.7 × 1013217.3219.9−3.4166.07.48 × 109160.1211.4−71.4
0.15228.36.2 × 1013222.3219.73.4164.25.46 × 109158.2211.5−74.1
0.2231.19.9 × 1013225.1219.77.2163.64.92 × 109157.6211.5−75.0
0.25231.91.1 × 1014225.8219.68.2161.63.52 × 109155.6211.6−77.9
0.3231.51.0 × 1014225.4219.67.7160.12.71 × 109154.0211.6−80.1
0.35231.61.1 × 1014225.4219.67.7158.82.17 × 109152.7211.7−82.0
0.4231.91.1 × 1014225.7219.68.1157.91.83 × 109151.7211.7−83.5
0.45231.51.0 × 1014225.3219.67.5157.11.61 × 109151.0211.7−84.5
0.5232.11.1 × 1014225.9219.68.3157.01.57 × 109150.8211.7−84.8
0.55231.71.1 × 1014225.5219.67.8156.91.54 × 109150.7211.7−85.0
0.6231.71.1 × 1014225.4219.67.7156.11.34 × 109149.9211.8−86.1
0.65231.09.6 × 1013224.8219.76.8156.81.52 × 109150.6211.7−85.1
0.7230.48.8 × 1013224.2219.76.0156.01.33 × 109149.8211.8−86.3
0.75230.69.1 × 1013224.4219.76.3156.51.45 × 109150.3211.7−85.5
0.8230.79.1 × 1013224.4219.76.3156.61.46 × 109150.3211.7−85.5
0.85230.58.8 × 1013224.2219.76.0156.51.43 × 109150.2211.7−85.7
0.9231.51.0 × 1014225.2219.67.4157.21.64 × 109150.9211.7−84.6
0.95231.61.1 × 1014225.3219.67.6158.21.95 × 109151.9211.7−83.2
Table 9. The activation energy, pre-exponential factor, and thermodynamic parameters obtained using the Friedman method at a 10 °C/min heating rate for host and 100phr loaded EPDM rubber.
Table 9. The activation energy, pre-exponential factor, and thermodynamic parameters obtained using the Friedman method at a 10 °C/min heating rate for host and 100phr loaded EPDM rubber.
EPDM 0 phr Lead at Rate 10 °C/minEPDM 100 phr Lead at Rate 10 °C/min
αEAΔHΔGΔSEAΔHΔGΔS
(kJ/mol)(sec−1)(kJ/mol)(kJ/mol)(J/mol)(kJ/mol)(sec−1)(kJ/mol)(kJ/mol)(J/mol)
0.05245.51.0 × 1015239.5219.327.0147.12.8 × 108141.3212.1−98.5
0.1240.14.3 × 1014234.1219.419.5150.04.7 × 108144.1212.0−94.5
0.15237.82.9 × 1014231.7219.516.3149.14.0 × 108143.1212.0−95.9
0.2236.02.2 × 1014229.9219.513.8145.62.2 × 108139.6212.2−101.1
0.25231.41.0 × 1014225.3219.67.5142.11.2 × 108136.0212.3−106.1
0.3229.77.7 × 1013223.6219.75.1139.37.2 × 107133.2212.4−110.3
0.35226.94.9 × 1013220.7219.81.3142.71.3 × 108136.6212.3−105.3
0.4228.56.4 × 1013222.4219.73.5146.62.6 × 108140.4212.1−99.8
0.45227.45.3 × 1013221.2219.82.0149.24.1 × 108143.1212.0−96.0
0.5229.87.9 × 1013223.6219.75.2149.94.5 × 108143.7212.0−95.1
0.55228.05.9 × 1013221.8219.72.8150.04.6 × 108143.8212.0−95.0
0.6225.33.8 × 1013219.1219.8−1.0151.15.7 × 108144.9212.0−93.3
0.65223.82.9 × 1013217.5219.9−3.1150.65.1 × 108144.3212.0−94.1
0.7224.93.5 × 1013218.6219.8−1.6151.76.2 × 108145.4211.9−92.6
0.75225.23.7 × 1013219.0219.8−1.1152.36.9 × 108146.0211.9−91.7
0.8227.15.1 × 1013220.9219.81.5152.36.9 × 108146.0211.9−91.7
0.85226.44.5 × 1013220.1219.80.4152.87.6 × 108146.5211.9−91.0
0.9230.08.2 × 1013223.7219.75.4155.61.2 × 109149.3211.8−86.9
0.95254.24.3 × 1015247.9219.138.4157.11.6 × 109150.7211.7−84.9
Table 10. Summary of the maximum, minimum, and average values of thermodynamic parameters (enthalpy change, Gibbs free energy change, and entropy change), activation energies (E), and pre-exponential factors (A).
Table 10. Summary of the maximum, minimum, and average values of thermodynamic parameters (enthalpy change, Gibbs free energy change, and entropy change), activation energies (E), and pre-exponential factors (A).
EPDM 0 phr Lead at Rate 10 °C/minEPDM 100 phr Lead at Rate 10 °C/min
KASMax225.03.6 × 101³218.9220.4−1.2162.03.5 × 109155.8211.9−77.5
Min206.01.5 × 1012199.7219.8−27.6152.06.8 × 108145.8211.6−91.9
Average222.22.6 × 101³216.0219.9−5.2155.11.4 × 109149.3211.8−86.9
FWOMax232.11.1 × 1014225.9220.38.3166.07.5 × 109160.1211.8−71.4
Min210.03.0 × 1012204.1219.6−21.5156.01.3 × 109149.8211.4−86.3
Average229.69.0 × 101³223.5219.75.0159.02.7 × 109152.9211.6−81.8
FMMax254.24.3 × 1015247.9219.938.4157.11.6 × 109150.7212.4−84.9
Min223.82.9 × 101³217.5219.1−3.1139.37.2 × 109133.2211.7−110.3
Average231.53.7 × 1014225.3219.77.5149.25.2 × 109143.1212.0−96.0
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

Alfannakh, H.; Alnaim, N.; Ibrahim, S.S. Thermal Stability and Non-Isothermal Kinetic Analysis of Ethylene–Propylene–Diene Rubber Composite. Polymers 2023, 15, 1890. https://doi.org/10.3390/polym15081890

AMA Style

Alfannakh H, Alnaim N, Ibrahim SS. Thermal Stability and Non-Isothermal Kinetic Analysis of Ethylene–Propylene–Diene Rubber Composite. Polymers. 2023; 15(8):1890. https://doi.org/10.3390/polym15081890

Chicago/Turabian Style

Alfannakh, Huda, Nisrin Alnaim, and Sobhy S. Ibrahim. 2023. "Thermal Stability and Non-Isothermal Kinetic Analysis of Ethylene–Propylene–Diene Rubber Composite" Polymers 15, no. 8: 1890. https://doi.org/10.3390/polym15081890

APA Style

Alfannakh, H., Alnaim, N., & Ibrahim, S. S. (2023). Thermal Stability and Non-Isothermal Kinetic Analysis of Ethylene–Propylene–Diene Rubber Composite. Polymers, 15(8), 1890. https://doi.org/10.3390/polym15081890

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