Next Article in Journal
Support Vector Machine-Based Global Classification Model of the Toxicity of Organic Compounds to Vibrio fischeri
Next Article in Special Issue
Study on Flame Retardancy Behavior of Epoxy Resin with Phosphaphenanthrene Triazine Compound and Organic Zinc Complexes Based on Phosphonitrile
Previous Article in Journal
Jumping in the Chiral Pool: Asymmetric Hydroaminations with Early Metals
Previous Article in Special Issue
Comprehensive Review of Recent Research Advances on Flame-Retardant Coatings for Building Materials: Chemical Ingredients, Micromorphology, and Processing Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pyrolysis Kinetics of Lignin-Based Flame Retardants Containing MOFs Structure for Epoxy Resins

Jiangsu Co-Innovation Center of Efficient Processing and Utilization of Forest Resources, Jiangsu Provincial Key Lab for the Chemistry and Utilization of Agro-Forest Biomass, College of Chemical Engineering, Nanjing Forestry University, Nanjing 210037, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2023, 28(6), 2699; https://doi.org/10.3390/molecules28062699
Submission received: 27 February 2023 / Revised: 11 March 2023 / Accepted: 14 March 2023 / Published: 16 March 2023
(This article belongs to the Special Issue Flame-Resistant Materials)

Abstract

:
This study describes the preparation of a lignin-based expandable flame retardant (Lignin-N-DOPO) using grafting melamine and covering 9,10-dihydro-9-oxa-10-phosphaphenanthrene-10-oxide (DOPO) using the Mannich reaction. Then, through in situ growth, a metal-organic framework (MOF) HKUST-1 (e.g., Cu3(BTC)2, BTC = benzene-1,3,5-tricarboxylate)/lignin-based expandable flame retardant (F-lignin@HKUST-1) was created. Before that, lignin epoxy resin containing phosphorus (P) and nitrogen (N) components had been created by combining epoxy resin (EP) with F-lignin@HKUST-1. Thermogravimetric analysis was used to examine the thermal characteristics of epoxy resin (EP) composite. The findings indicate that the thermal stability of EP is significantly affected by the presence of F-lignin@HKUST-1. Last but not least, the activation energy (E) of EP/15% F-lignin@HKUST-1 was examined using four different techniques, including the Kissinger-SY iteration method, the Ozawa-SY iteration method, the Lee-Beck approximation-iteration method, and the Gorbatchev approximation-iteration method. It was discovered that the activation energy was significantly higher than that of lignin. Higher activation energy suggests that F-lignin@HKUST-1 pyrolysis requires more energy from the environment, which will be significant about the application of lignin-based flame retardants.

1. Introduction

Epoxy resin (EP), as one of the important thermosetting plastics, has high mechanical strength, excellent adhesion, chemical resistance, and electrical insulation, etc., and has been widely used in construction, automotive, electronics, and aerospace fields [1,2,3]. However, the inherent flammability and smoke emission of EP greatly limit its application in some specific areas [4,5,6]. Therefore, the flame-retardant modification of EP is of great importance.
Materials made of biomass offer various benefits, including being renewable [7], relatively non-polluting [8], widely available, and carbon neutral. Animals, plants, and microbes make up biomass, which is employed in a variety of industries such as the pharmaceutical [9], industrial [10], and packaging sectors [11]. Because biomass produces a lot of carbon after burning, lignin has been employed in the field of flame retardants as a carbon-forming agent [12,13,14,15,16].
Mandlekar et al. [17] investigated the patterns of thermal stability and flame retardancy of bio-based polyamide 11 (PA) having different lignin contents. Microcomposites based on PA and containing 5, 10, 15, and 20 wt.% lignin were prepared using a twin-screw extruder. Morphological analysis showed that the resulting microcomposites had good interfacial interactions and uniform distribution of lignin particles. In addition, thermogravimetric analysis in an inert atmosphere showed that sulfated lignin was able to produce lower coke residues (41–48 wt.% at 600 °C), and unlike sulfated lignin, sulfonated lignin was able to provide higher thermal stability as well as higher coke residues (55–58 wt.%). In addition, vertical flame-spread tests clearly showed that 15 wt.% was the optimum value for kraft or sulfonated lignin loading to achieve improved flame retardancy and V1 rating. In addition, cone-calorimetric tests were used to investigate the forced combustion behavior; in particular, the microcomposites containing sulfonated lignin showed significantly lower peak heat-release rates (−51%) and total heat release (−23%), lower flat-mass loss rates, and a significant mass increase in the final residue (~9 wt.%). In contrast, microcomposites containing kraft lignin showed the opposite effect, as the heat release rate (HRR) and total heat release (THR) values increased in the presence of kraft lignin.
Despite the fact that lignin may be applied directly to materials to act as a flame retardant, the result frequently falls short of expectations, hence it is typically modified chemically. Due to their abundance and strong reactivity, the hydroxyl groups (phenolic or aliphatic hydroxyl groups) in the lignin structure permit chemical change. The use of PN-Lignin as a flame retardant for many polymeric materials, including polypropylene (PP), polylactic acid (PLA), EP, and polyurethane (PU) materials, is frequently investigated. In flame-retardant systems, elemental nitrogen serves as a gas source and expands the material. In flame retardants, elemental phosphorus serves as an acid source and can also be employed as a dehydrating agent, most frequently in the form of POCI3, polyphosphoric acid, DOPO (9,10-dihydro-9-oxa-10-phosphaphenanthrene-10-oxide), and its derivatives. These phosphorus-containing compounds can generate acids that, when heated to a high temperature, dehydrate the carbon layer and esterify polyols, causing the carbon source to create a phosphorus-rich charred layer.
The flame-retardants effect of lignin modified by phosphorus and nitrogen has been improved, but its flame-retardant efficiency and smoke-suppression performance sometimes cannot meet the demand, for example not effectively reducing the release of CO and the emission of carbon particles, which will still have a considerable impact on the environment. According to the research in recent years, metal-organic frameworks(MOFs) materials have good adsorption effect on smoke generation and good catalytic effect on CO oxidation due to their structural characteristics. Zhang et al. [18] used MOFs as precursors to link two different LDHs into composite nanomaterials, and the specimens passed the V-0 rating in the UL-94 test. Xu et al. [19] prepared functionalized reduced graphene oxide (RGO) with Co-ZIF (zeolite imidazolium framework-67) adsorbed borate ions (ZIF-67/RGO-B). The peak heat-release rate (PHRR), total heat release rate (THR), and maximum smoke concentration (Ds, max) of the composites with ZIF-67/RGO-B doping of 2 wt.% were reduced by 65.1%, 41.1% and 66.0%, respectively, compared to the pure EP. Hou et al. [20] successfully synthesized iron- and cobalt-based MOFs and incorporated them into polystyrene (PS) as flame retardants. Thermogravimetric and cone-calorimetric analyses showed that the MOFs had good flame-retardant effects. Hou et al. [21] synthesized a Co-based metal organic backbone (P-MOF) with a phosphorus-based structure using a facile hydrothermal reaction and added the P-MOF to an epoxy resin (EP) to enhance its flame retardancy. Results from the conical calorimeter and steady-state tube-furnace tests showed a significant reduction in total flue-gas yield and total CO.
The purpose of this study was to develop an EP flame retardant that is synthetically simple, utilizing alkaline lignin, melamine, DOPO, and HKUST-1 (a metal-organic framework (MOF), Cu3(BTC)2, BTC = benzene-1,3,5-tricarboxylate). To create lignin-based flame retardants, alkaline lignin was altered in this study to enhance the reaction sites. Then, a phosphor-nitrogen lignin-based flame retardant was created by reacting the phenolized modified lignin with DOPO and melamine complex (Lignin-N-DOPO), and then HKUST-1/lignin-based flame retardant (F-lignin@HKUST-1) was prepared by introducing the MOFs structure as its basis. The flame-retardant properties were significantly enhanced by the addition of EP, and with the use of the Kissinger-SY iterative approach, the Ozawa-SY iterative method, the Lee-Beck approximation, and the Gorbatchev approximation, the kinetics of its pyrolysis were examined.

2. Results

2.1. Thermogravimetric Analysis

The thermal-stability and thermal-degradation processes of EP and its composites involved are investigated using thermodynamic analysis curves (from room temperature to 700 °C). The TG and DTG curves of EP and its composites are depicted in Figure 1, and Table 1 contains information on the starting thermal decomposition temperature (Td5%), the maximum thermal decomposition rate temperature (Tdmax), and the quantity of residual carbon. Figure 1a illustrates the one-step thermal deterioration of the pure EP sample in N2. The decomposition begins at 333.7 °C and reaches its maximum rate at 405.3 °C, leaving only 16.21 wt.% of carbon behind. The thermogravimetric curves of EP/15% F-Lignin@HKUST-1 are similar to those of pure EP, but exhibit different thermal stability.
Figure 1 shows that the EP/15%F-lignin@HKUST-1 started to degrade before the EP, as may be inferred. Tdmax showed a trend toward decrease, as shown in Table 1, and this finding can be related to the early phosphorus breakdown, the quick DOPO breakdown, and the catalytic action of the breakdown products of MOFs, which makes it easier to create fast-catalytic carbon residues. Additionally, it can be seen from Table 1 that the degradation temperature range of EP/F-lignin@HKUST-1 is smaller, indicating that it can complete the decomposition more quickly and act as a catalyst and carbon promoter [22]. These results show that the addition of a MOFs structure can improve the residual carbon amount of the composite as shown by the higher residual carbon amount at 700 °C compared to pure epoxy resin. These findings suggest that the synergistic interaction between DOPO and MOFs can facilitate the carbon-production process in thermoset composites based on epoxy.

2.2. Kinetic Analysis

The apparent activation energy at the appropriate conversion is frequently determined by linear-regression analysis utilizing the iso-conversion method using the pyrolysis curves produced by the pyrolysis process of solid-state reactants. The most crucial step is picking an appropriate approximation for the temperature integral. With different conversions, the Lee-Beck-approximation and Gorbatchev-approximation iterative methods are employed for P(u), while the third-order expression of Senum-Yang is selected for the temperature-integral function P(u) and iterated by the Kissinger-SY and Ozawa-SY iterative methods. The apparent activation energy Eα at different conversions can be obtained by iteratively calculating the pyrolysis conversion curves of samples Al-lignin, EP/15% F-Lignin@HKUST-1 in 0.05 steps from 0.15–0.90, etc. Based on the aforementioned four methods of calculation, Table 2 and Table 3, respectively, present the apparent activation energy and R2 (regression coefficient) at various conversions.
Table 2 and Table 3 show that the regression coefficients R2 for the apparent activation energy of the Al-lignin pyrolysis process at various conversions. These values may be impacted by experimental errors and result in a mediocre fit, whereas the R2 of the regression coefficients at different conversions. With the increased conversion (0.20 < α < 0.80), the regression coefficients’ R2 values ranged from 0.95 to 0.99. This indicates that the experimental data of sample F-Lignin@HKUST-1 is good on the one hand, while on the other hand the pyrolysis process of this sample is probably a single-step reaction kinetic. The apparent activation-energy values acquired by the four iterative approaches are quite near to each other, and the fit of the data obtained from the simultaneous analysis of the three distinct materials using the four iterative methods is high, suggesting that the findings computed by the four iterative methods are more accurate. Linear-regression analysis plots of iso-conversions of samples Al-lignin and F-Lignin@HKUST-1 are shown below in Figure 2 and Figure 3, respectively.
Table 2 and Table 3 show that the apparent activation energies Eα obtained for the two samples using the four different calculation methods are nearly identical. One of the four calculation methods can be used to determine the apparent activation-energy (Eα) value of lignin, which can then be used to study the activation-energy law during lignin pyrolysis. According to the Ozawa-iterative kinetic model, Figure 4 depicts the curve of conversion versus activation energy for the samples Al-lignin and EP/15% F-Lignin@HKUST-1, respectively. The energy barrier needed for chemical reactions to take place is known as the activation energy; the higher the activation energy, the more challenging the reaction. It establishes the reaction rate’s responsiveness and sensitivity [23]. As a result, the various activation energies at various conversions highlight both the complexity of the reactions and the multi-stage nature of the thermal degradation of solid materials. It also demonstrates the intricacy of the two materials’ pyrolysis and chemical changes, which should have involved many reactions at various phases.
Figure 4 shows that the activation energy needed for the pyrolysis of the sample Al-lignin is low at the early stages of decomposition (0.05 < α < 0.10), which is because the initial process primarily involves the evaporation of water as well as some small molecules, and the energy needed for the evaporation of water is typically low. The local extreme points of EP/15% F-lignin@HKUST-1’s energy versus conversion curve demonstrate that its water content was lower, and its water-evaporation phases were shorter, than in the aforementioned TG and DTG studies; all of them were roughly in the α < 0.10 interval. Similar changes in activation energy were seen in the sample Al-lignin during the initial stages of decomposition.
Both samples begin to reach the severe-pyrolysis phase, which breaks numerous chemical bonds and demands more energy. Lignin pyrolysis is a heat-absorbing process. Al-lignin, one of the samples, had a quick apparent activation energy in the interval of conversion 0.1 < α < 0.4. It eventually rose from around 20 kJ/mol to a range of 100–120 kJ/mol and stayed steady. The apparent activation energy of sample EP/15% F-lignin@HKUST-1 increased rapidly from about 70 kJ/mol to about 130 kJ/mol in the range of conversion 0 < α < 0.1 and remained stable.
The value of the sample Al-lignin activation energy fluctuated slightly as the conversion increased (0.50 < α < 0.70), and the reaction-activation energy showed a downward trend. This was explained by the fact that the formal pyrolysis process of lignin reached its final stage and the chemical bond-breaking process in the lignin molecule slowed down. The activation energy needed for lignin pyrolysis again exhibits a steep rising trend as the conversion rises (0.70 < α < 0.80). This is because the inorganic salts in the lignin molecule start to break down, which consumes a lot of energy. The inorganic salt degradation process was reaching its conclusion when the conversion kept rising (0.80 < α < 0.90), and as a result, the apparent activation energy of lignin began to decline.
The conversion for the sample EP/15% F-lignin@HKUST-1, however, very slightly increased over the range of 0.15 to 0.85 and stayed constant at 130 kJ/mol. The inclusion of flame-retardant components to the modified materials raised the energy barrier during the pyrolysis of lignin, resulting in greater activation energies for the main pyrolysis phases of those materials than for Al-lignin. Also, the difference of the added flame retardant resulted in the activation energy of the main pyrolysis stage of sample EP/15% Lignin-N-DOPO being greater than that of the main pyrolysis stage of sample EP/15% F-lignin@HKUST-1. It can be shown that sample EP/15% F-lignin@HKUST-1 had a shorter main pyrolysis phase, a bigger pyrolysis rate than the former, and a higher fluctuation in conversion when combined with the TG and DTG assays on the modified lignin samples above.
The modified lignin samples’ activation energy versus conversion curves showed some degree of stability throughout a range of conversions, suggesting that the pyrolysis process may be adequately captured by a single-step reaction kinetic-mechanism equation. A multi-step reaction kinetic equation could be required since the activation energy of Al-lignin was increasingly complex with the conversion, and since the pyrolysis process could not be well characterized by a single-step reaction kinetic equation.

2.3. Mechanic Equations

The phases of conversion α 0.1 to 0.9 were replaced into the universal integral Equation (1) and differential Equation (2), respectively, in four sets of lignin pyrolysis data with varying heating rates, and the Avrami-Erofeev equation was chosen as the mechanism function (AEn).
ln G α T T 0 = ln A β E R T
Avrami-Erofeev equation (AEn):
ln d α d T f α E T T 0 R T 2 + 1 = ln A β E R T
The lignin pyrolysis data were calculated using the mechanism functions as equations: AE4, AE5, AE6.
The universal integral equation and differential equation were used for the four sets of EP/15% F-lignin@HKUST-1 pyrolysis data with various heating rates. The mechanism functions were chosen as equations AE2, AE3, and AE4, respectively, and the apparent activation energy produced from the four iterative calculation techniques was utilized to compare and filter the levels of the mechanism functions. The output of the computation is displayed in Table 4 and Table 5 below.
The calculation results in Table 4 and Table 5 demonstrate that the kinetic parameters acquired by the differential-equation approach and the universal-integral method are quite near to one another, demonstrating the relative reliability of the findings produced by the two methods. Table 4 shows the apparent pyrolysis activation energies of Al-lignin, which vary from 65 to 76 kJ/mol, 84 to 98 kJ/mol, and 103 to 120 kJ/mol, respectively. R2, the fitting coefficient, is also more than 0.95. The average apparent activation energy Eα of lignin pyrolysis, as shown in Table 2, is 99.82 kJ/mol, making AE6 the most likely candidate for the mechanism function. Table 5 Mechanism functions AE2, AE3 and AE4 are respectively applied. The apparent activation energies of EP/15% F-lignin@HKUST-1 pyrolysis of samples are 73–113 kJ/mol, 115–175 kJ/mol and 158–237 kJ/mol, respectively, and the fitting coefficients R2 are all higher than 0.9. The apparent activation energy of pyrolysis of the sample obtained in Table 3 above is mostly between 128 kJ/mol and 145 kJ/mol, so the most likely mechanism function is determined as AE3.
The most probable pyrolysis-mechanism function of Al-lignin is AE6, and that of EP/15% F-lignin@HKUST-1 is AE3 according to the universal integral and differential equation methods.

2.4. Compensation Effect

An essential component of the research of thermal-analysis kinetics is the kinetic compensatory effect. By using the equal-conversion approach, the appropriate pre-exponential factor lnA may be determined once the mechanism function of the pyrolysis process of solid reactants has been established. In Table 6 and Table 7, the apparent activation energy Eα and preexponential factor lnA of the lignin samples EP/15% F-lignin@HKUST-1 and Al-lignin are shown, respectively, based on the Kissinger-SY iteration approach and the Ozawa-SY iteration methods, where the mechanism function is the most plausible mechanism function corresponding to it.
When the mechanism function is the most likely, Figure 5 and Figure 6 illustrate the related compensating-effect curves for lignin Al-lignin and sample EP/15%F-lignin@HKUST-1 based on Kissinger-SY iteration technique and Ozawa-SY iteration method, respectively. Figure 5 shows that in the pyrolysis of the three materials, the apparent activation energies Eα and lnA exhibit some linear connection.
Figure 5 illustrates this relationship between the apparent activation energy, E, during lignin pyrolysis and the preexponential component lnA. This relationship occurs under the assumption that mechanism function AE6, lnA = 0.1589E − 1.9599, with an R2 of 0.88, is the equation that fits using the Kissinger-SY technique. The fitted equation using the iterative Ozawa-SY approach is lnA = 0.1589E − 1.9595 with an R2 of 0.88. According to the Kissinger-SY iteration approach, the fitted equation for sample EP/15%F-lignin@HKUST-1 is lnA = 0.3574E − 25.0022, with R2 equal to 0.88. According to the Ozawa-SY iteration method, the fitted equation is lnA = 0.3574E − 25.0021, with R2 equal to 0.88. The two samples’ calculation outcomes using the two approaches are almost identical, demonstrating the validity of the calculation outcomes. All of the R2 coefficients for linear fitting were larger than 0.88, although there was no evidence of substantial linear association. This suggests that the two samples’ pyrolysis processes are rather complicated, and it is challenging to characterize the reaction-mechanism function of pyrolysis by a mechanism function with matching physical relevance. This is due to the fact that several various chemical processes take place during the pyrolysis of materials, and the power series of the mechanism function may be rather complicated.

3. Experiments

3.1. Materials

Alkaline lignin (Al-lignin) was purchased from TCI Shanghai Co., Ltd., Shanghai, China. DOPO, carboxymethyl cellulose sodium (CMC) and 4,4-diaminodiphenylmethane (DDM) were provided by Shanghai Macklin Co., Ltd., Shanghai, China. Diglycidyl ether of bisphenol A type epoxy resin (commercial name: E 51) was purchased from Xingchen Epoxy Resins Factory (Nantong, China). Melamine (MEL) was obtained from Shanghai Lingfeng Co., Ltd., Shanghai, China. Formaldehyde, N,N-dimethylformamide (DMF), anhydrous copper acetate (Cu(OAc)2), trimesic acid (H3BTC,C9H6O6), glacial acetic acid (CH3COOH), diethyl ether, ethanol, and phenol were purchased from Nanjing Chemical Reagent Co., Ltd., Nanjing, China.

3.2. Phenolation of Alkaline Lignin(Ph-Lignin)

We prepared 2 mol/L sulfuric acid (H2SO4) solution, took 80 mL with a three-neck flask and set the temperature at 80 °C. Then we added 20 g of lignin (Al-lignin) and stirred at condensation reflux for 1.5 h. Then we slowly added 18 g of phenol, dropwise, during the temperature rise to 95 °C and stirred at condensation reflux for 1.5 h at 95 °C. After the reaction was completed, it was cooled to room temperature, the solid was rinsed three times using ether, and then dried at 70 °C for 12 h.

3.3. Preparation of MEL and DOPO Mixture(MEL-DOPO)

We combined 280 mL of ethanol aqueous solution (70 wt.%) and 12.96 g of DOPO (0.06 mol), and the mixture was agitated for 30 minutes to thoroughly dissolve the DOPO. When the temperature reached 70 °C, 10 g of MEL (0.08 mol) was added and agitated for 6 h. Following completion, the mixture was cooled to ambient temperature, three times washed with ethanol, and then dried for 12 h in a vacuum oven at 70 °C.

3.4. Preparation of PN-Lignin

We added 4 g of phenolized lignin (Ph-lignin) and 20 g of MEL-DOPO to a 50 mL single-mouth flask, then 300 mL of N,N dimethylformamide (DMF) solution was added and the temperature was set to 75 °C. We added 7.2 g (0.024 mol) of formaldehyde to react for 3 h. After being cooled to room temperature, the solid was repeatedly washed several times by adding an appropriate amount of distilled water, and then the solid was dried in a drying oven at 70 °C for 24 h.

3.5. Preparation of PN-Lignin@HKUST-1

First of all, exactly 1.467 g of CMC (6.057 mol) with 2 g of PN-lignin was added into 70 mL of deionized water and mechanically stirred at room temperature for 1 h. Secondly, 1.326 g of Cu(OAc)2 (6.642 mmol) dissolved in 30 mL of deionized water and 1 mL of HAc were added dropwise and stirred for 5 min. Thirdly, 1.155 g of H3BTC (5.496 mmol) dissolved in 15 mL ethanol was sequentially added and the mixture was continuously stirred in a closed system for 4 h. Finally, the end product was obtained by centrifugation and dried in the vacuum oven at 70 °C for 24 h.

3.6. Preparation of Flame-Retardant EP Composites

A traditional curing method was applied to fabricate epoxy composites with some modifications. Typically, 3 g of PN-lignin@HKUST-1 and 15 g of EP were mechanically stirred at 1000 rpm for 1.5 h to obtain a uniform mixture system. Next, 3 g of curing agent DDM (mass ratio, DDM/EP = 1:5) were added and stirred for 1.5 h. The uniformly-mixed thermoset was poured into the mould and cured at 100 °C for 2 h after bubbles were removed in the vacuum oven. For further curing, the temperature was raised to 150 °C for 2 h. To evaluate the effect of modification, the neat EP thermoset and EP composites containing PN-lignin (referred to as EP/PN-lignin) were also prepared according to the similar procedure.

3.7. Measurements

Using the DTG-60AH, the thermogravimetric analysis (TGA) was achieved (SHIMADZU, Kyoto, Japan). The temperature range was from 30 °C to 900 °C, the sample value range was 5–10 mg, nitrogen was used as the protective gas, and the gas flow rate was set at 20 mL/min.
Al-lignin was tested at 8 °C/min, 12 °C/min, 16 °C/min and 20 °C/min. The ramp rates were set at 5 °C/min, 10 °C/min, 15 °C/min and 20 °C/min for EP/15% Flignin@HKUST-1.

3.8. Kinetic Modeling

The non-isothermal constant-rate-of-warming approach is frequently employed in thermal analysis to analyze the kinetic behavior of pyrolysis of solid-state reactants. As demonstrated in Equation (3), the thermal decomposition rate for the pyrolytic behavior of solid-state reactants can be stated as a function of temperature k(T) and conversion f(α):
d α d t = β d α d T = k T f α
where
α is the conversion of the sample (%)
β is the constant temperature rise rate (°C/min)
T is the reaction temperature (Kelvin, K)
where the conversion α can be obtained by Equation (4):
α = M 0 M T M 0 M
where
M0 is the initial mass of the sample (%)
MT is the mass of the sample at the decomposition temperature (%)
M is the mass of the sample after decomposition (%)
The function k(T) can generally be obtained from the Arrhenius equation:
k T = A × exp E R T
where
A is the prefactor
E is the apparent activation energy (KJ/mol)
R is the gas constant (8.3145 J/mol/K)
To generate Equation (7) by integrating the left and right sides of the equation from 0 to and T0 to T, respectively, and designating the integral kinetic mode function as G(α), Equation (5) can be substituted into 1 and distorted to produce Equation (6). The initial temperature T0 of the reaction is low, and consequently, the reaction rate is very small.
d α f α = A β exp E R T d T
G α = 0 α d α f α = 0 T A β exp E R T d T
If it is stipulated that:
u = E R T ,   p u = u e u u 2 d u
The class I kinetic equations for thermal analysis can be obtained (Equation (8)).
G α = 0 α d α f α = 0 T A β exp E R T d T = A E β R u e u u 2 d u = A E β R P u
If P(u) = PFK(u) = e − u/u2 or take the Doyle approximation, i.e., P(u) = PD(u) = 0.00484e − 1.0516u, the Kissinger equation (Equation (9)) and the Ozawa equation (Equation (10)) can be obtained by deformation when substituting Equation (8), respectively.
ln β T = ln A R E G α E R T
ln β = ln 0.00484 A E R G α 1.0516 E R T
If the function P(u) takes the Senum-Yang approximation (P(u) = PSY(u) = [(e − u/u2) × Q4(u)]), where
Q 4 u = u 3 + 18 u 2 + 88 u + 96 u 4 + 20 u 3 + 120 u 2 + 240 u + 120
H ( u ) = P S Y u P D u = e u Q 4 u 0.00484 u 2 e 1.0516 u
Immediately combining the aforementioned equations yields the Kissinger-Senum-Yang-iteration equation (Equation (13)) and the Ozawa-Senum-Yang-iteration equation (Equation (14)).
ln β Q 4 u T 2 = ln A R E G α E R T
ln β H ( u ) = ln 0.00484 A E R G α 1.0516 E R T
The Kissinger-iterative equation and the Ozawa-iterative equation, both of which take into account how H(u) and Q(u) fluctuate slowly with u without being constrained by the range of u, are both members of the equal conversion integration method. To locate E, the iterative technique is utilized. To find E, first set H(u) = 1 or Q4(u) = 1. The starting value of E is then determined from the slope using the least-squares approach and the linear relationship between ln(1/T2) − 1/T and ln(−1/T) using the least-squares method. The new value of E is then calculated from the slope by substituting the initial value into H(u) and Q(u) using the linear connection with ln(1/H(u) − 1/T and ln[1/Q(u)T2] − 1/T, respectively. Another corrected value can be derived by iterating once more using this corrected value of E as the initial value. A more logical E-value fulfilling less than 0.1 kJ/mol is thus reached after a number of repetitions.

4. Conclusions

In this study, Al-lignin and EP/15% F-lignin@HKUST-1 samples were analyzed kinetically using the Kissinger-SY, Ozawa-SY, Lee-Beck approximation, and Gorbatchev approximation iteration methods. The differential-equation method and the universal-integral methods were used to calculate the two samples’ most likely mechanism functions.
  • Four kinetic methods were used to determine the activation energies of the Al-lignin and EP/15% F-lignin@HKUST-1 samples. The four procedures were reasonable since they produced outcomes that were almost identical. Al-lignin pyrolysis has a wide range of apparent activation energies and can be separated into three stages. The complexity of the pyrolysis process is reflected in the multiple increases and drops in the apparent activation energy of al-lignin pyrolysis. The change of apparent activation energy during pyrolysis of EP/15% F-lignin@HKUST-1 samples can also be divided into three stages, but different from lignin, the apparent activation energy during pyrolysis of these two samples experienced three stages: rapid increase, slow change, and rapid rise, reflecting that their pyrolysis-rate curves only have a single obvious rate peak.
  • Using universal integral and differential equation approaches, the most likely mechanism functions for the pyrolysis of Al-lignin and EP/15% F-lignin@HKUST-1 were identified as AE6, AE4, and AE3. The regression-intercept of the Kissinger-SY and Ozawa-SY iterative approaches were used to get the relevant preexponential component lnA, and the related apparent activation energy Eα had a strong linear relationship.

Author Contributions

Conceptualization, T.Y., J.C. and X.G.; methodology, T.Y., R.Y., H.Y., J.C. and X.G.; software, C.S.; validation, T.Y., J.C. and X.G.; formal analysis, T.Y., R.Y., C.S., Y.L., R.L., J.C. and X.G.; investigation, T.Y., Y.L. and R.L.; data curation, T.Y., J.C. and X.G.; writing—original draft preparation, T.Y.; writing—review and editing, J.C. and X.G..; visualization, T.Y., J.C. and X.G.; supervision, J.C. and X.G.; funding acquisition, X.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported from the National Natural Science Foundation of China (no. 21774059), the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions and the College Students’ Practice and Innovation Training Project (202210298002Z, 202210298003K).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Analysis and testing facility was provided by the advanced analysis and testing center of Nanjing Forestry University, China.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Kumar, S.; Krishnan, S.; Samal, S.K.; Mohanty, S.; Nayak, S.K. Toughening of petroleum based (DGEBA) epoxy resins with various renewable resources based flexible chains for high performance applications: A review. Ind. Eng. Chem. Res. 2018, 57, 2711–2726. [Google Scholar] [CrossRef]
  2. Wang, X.; Guo, W.; Song, L.; Hu, Y. Intrinsically flame retardant bio-based epoxy thermosets: A review. Compos. Part B Eng. 2019, 179, 107487. [Google Scholar] [CrossRef]
  3. Jin, F.-L.; Li, X.; Park, S.-J. Synthesis and application of epoxy resins: A review. J. Ind. Eng. Chem. 2015, 29, 1–11. [Google Scholar] [CrossRef]
  4. Liu, Q.; Zhao, Y.; Gao, S.; Yang, X.; Fan, R.; Zhi, M.; Fu, M. Recent advances in the flame retardancy role of graphene and its derivatives in epoxy resin materials. Compos. Part A Appl. Sci. Manuf. 2021, 149, 106539. [Google Scholar] [CrossRef]
  5. Luo, H.; Rao, W.; Zhao, P.; Wang, L.; Liu, Y.; Yu, C. An efficient organic/inorganic phosphorus–nitrogen–silicon flame retardant towards low-flammability epoxy resin. Polym. Degrad. Stab. 2020, 178, 109195. [Google Scholar] [CrossRef]
  6. Fang, F.; Ran, S.; Fang, Z.; Song, P.; Wang, H. Improved flame resistance and thermo-mechanical properties of epoxy resin nanocomposites from functionalized graphene oxide via self-assembly in water. Compos. Part B Eng. 2019, 165, 406–416. [Google Scholar] [CrossRef]
  7. Wang, S.; Bai, J.; Innocent, M.T.; Wang, Q.; Xiang, H.; Tang, J.; Zhu, M. Lignin-based carbon fibers: Formation, modification and potential applications. Green Energy Environ. 2021, 7, 578–605. [Google Scholar] [CrossRef]
  8. Liu, M.; Aravind, P. The fate of tars under solid oxide fuel cell conditions: A review. Appl. Therm. Eng. 2014, 70, 687–693. [Google Scholar] [CrossRef]
  9. Liao, J.J.; Abd Latif, N.H.; Trache, D.; Brosse, N.; Hussin, M.H. Current advancement on the isolation, characterization and application of lignin. Int. J. Biol. Macromol. 2020, 162, 985–1024. [Google Scholar] [CrossRef]
  10. Galkin, M.V.; Samec, J.S. Lignin valorization through catalytic lignocellulose fractionation: A fundamental platform for the future biorefinery. ChemSusChem 2016, 9, 1544–1558. [Google Scholar] [CrossRef]
  11. Asgher, M.; Qamar, S.A.; Bilal, M.; Iqbal, H.M. Bio-based active food packaging materials: Sustainable alternative to conventional petrochemical-based packaging materials. Food Res. Int. 2020, 137, 109625. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, H.; Yu, B.; Xu, X.; Bourbigot, S.; Wang, H.; Song, P. Lignin-derived bio-based flame retardants toward high-performance sustainable polymeric materials. Green Chem. 2020, 22, 2129–2161. [Google Scholar] [CrossRef]
  13. Wu, Q.; Ran, F.; Dai, L.; Li, C.; Li, R.; Si, C. A functional lignin-based nanofiller for flame-retardant blend. Int. J. Biol. Macromol. 2021, 190, 390–395. [Google Scholar] [CrossRef] [PubMed]
  14. Liang, D.; Zhu, X.; Dai, P.; Lu, X.; Guo, H.; Que, H.; Wang, D.; He, T.; Xu, C.; Robin, H.M. Preparation of a novel lignin-based flame retardant for epoxy resin. Mater. Chem. Phys. 2021, 259, 124101. [Google Scholar] [CrossRef]
  15. Dai, P.; Liang, M.; Ma, X.; Luo, Y.; He, M.; Gu, X.; Gu, Q.; Hussain, I.; Luo, Z. Highly efficient, environmentally friendly lignin-based flame retardant used in epoxy resin. ACS Omega 2020, 5, 32084–32093. [Google Scholar] [CrossRef]
  16. Chen, S.; Lin, S.; Hu, Y.; Ma, M.; Shi, Y.; Liu, J.; Zhu, F.; Wang, X. A lignin—Based flame retardant for improving fire behavior and biodegradation performance of polybutylene succinate. Polym. Adv. Technol. 2018, 29, 3142–3150. [Google Scholar] [CrossRef]
  17. Mandlekar, N.; Cayla, A.; Rault, F.; Giraud, S.; Salaun, F.; Malucelli, G.; Guan, J. Thermal stability and fire retardant properties of polyamide 11 microcomposites containing different lignins. Ind. Eng. Chem. Res. 2017, 56, 13704–13714. [Google Scholar] [CrossRef]
  18. Zhang, Z.; Qin, J.; Zhang, W.; Pan, Y.-T.; Wang, D.-Y.; Yang, R. Synthesis of a novel dual layered double hydroxide hybrid nanomaterial and its application in epoxy nanocomposites. Chem. Eng. J. 2020, 381, 122777. [Google Scholar] [CrossRef]
  19. Xu, W.; Wang, X.; Wu, Y.; Li, W.; Chen, C. Functionalized graphene with Co-ZIF adsorbed borate ions as an effective flame retardant and smoke suppression agent for epoxy resin. J. Hazard. Mater. 2019, 363, 138–151. [Google Scholar] [CrossRef]
  20. Hou, Y.; Hu, W.; Gui, Z.; Hu, Y. Preparation of metal–organic frameworks and their application as flame retardants for polystyrene. Ind. Eng. Chem. Res. 2017, 56, 2036–2045. [Google Scholar] [CrossRef]
  21. Hou, Y.; Hu, W.; Gui, Z.; Hu, Y. A novel CoSII)–based metal-organic framework with phosphorus-containing structure: Build for enhancing fire safety of epoxy. Compos. Sci. Technol. 2017, 152, 231–242. [Google Scholar] [CrossRef]
  22. Liao, S.-H.; Liu, P.-L.; Hsiao, M.-C.; Teng, C.-C.; Wang, C.-A.; Ger, M.-D.; Chiang, C.-L. One-step reduction and functionalization of graphene oxide with phosphorus-based compound to produce flame-retardant epoxy nanocomposite. Ind. Eng. Chem. Res. 2012, 51, 4573–4581. [Google Scholar] [CrossRef]
  23. Ounas, A.; Aboulkas, A.; Bacaoui, A.; Yaacoubi, A. Pyrolysis of olive residue and sugar cane bagasse: Non-isothermal thermogravimetric kinetic analysis. Bioresour. Technol. 2011, 102, 11234–11238. [Google Scholar] [CrossRef] [PubMed]
Figure 1. TG (a) and DTG (b) curves of pure EP and EP composites under N2 atmosphere.
Figure 1. TG (a) and DTG (b) curves of pure EP and EP composites under N2 atmosphere.
Molecules 28 02699 g001
Figure 2. Linear regression fitting diagram of the equal conversion method of Al-lignin, based on four temperature approximations.
Figure 2. Linear regression fitting diagram of the equal conversion method of Al-lignin, based on four temperature approximations.
Molecules 28 02699 g002
Figure 3. Linear regression fitting diagram of equal conversion method of EP/15% F-Lignin@HKUST-1 based on four temperature approximation.
Figure 3. Linear regression fitting diagram of equal conversion method of EP/15% F-Lignin@HKUST-1 based on four temperature approximation.
Molecules 28 02699 g003
Figure 4. Relationship between activation energy and conversion of Al-lignin (a) and EP/15% F-lignin@HKUST-1.
Figure 4. Relationship between activation energy and conversion of Al-lignin (a) and EP/15% F-lignin@HKUST-1.
Molecules 28 02699 g004
Figure 5. Linear correlation between the activation energy of Al-lignin and the preinhetial factor based on two iterative methods.
Figure 5. Linear correlation between the activation energy of Al-lignin and the preinhetial factor based on two iterative methods.
Molecules 28 02699 g005
Figure 6. Linear correlation between the activation energy of EP/15% F-lignin@HKUST-1 and the preinhetial factor based on two iterative methods.
Figure 6. Linear correlation between the activation energy of EP/15% F-lignin@HKUST-1 and the preinhetial factor based on two iterative methods.
Molecules 28 02699 g006
Table 1. TG data of EP, Al-lignin and EP composites under N2 atmosphere.
Table 1. TG data of EP, Al-lignin and EP composites under N2 atmosphere.
SamplesTd5% (°C)Tdmax (°C)Residues (wt.%)
EP333.7405.316.21
Al-lignin105.5333.341.57
EP/15%F-lignin@HKUST-1334.1360.720.14
Table 2. Activation energy of Lignin at different conversions.
Table 2. Activation energy of Lignin at different conversions.
αKissiger-SY IterationOzawa-SY IterationLee-Beck IterationGorbatchev Iteration
E/(KJ/mol)R2E/(KJ/mol)R2E/(KJ/mol)R2E/(KJ/mol)R2
0.1557.610.9257.610.9357.670.9257.670.92
0.2092.190.9692.190.9792.220.9692.220.96
0.25112.610.98112.620.99112.640.98112.640.98
0.30122.640.99122.640.99122.670.99.122.670.99
0.35126.600.99126.600.99126.620.99126.620.99
0.40127.890.99127.890.99127.920.99127.920.99
0.45127.730.99127.730.99127.760.99127.760.99
0.50123.280.99123.280.99123..310.99123.310.99
0.55111.540.99111.540.99111.580.99111.580.99
0.60100.410.97100.420.97100.470.97100.470.97
0.6564.240.8964.240.9064.370.8964.370.89
0.7049.340.9149.340.9249.560.9249.560.92
0.75129.120.97129.120.97129.190.97129.190.97
0.80148.400.96148.410.97148.470.96148.470.96
0.85124.020.95124.020.95124.120.95124.120.94
0.9071.870.8471.870.8572. 100.8472.100.84
Table 3. Activation energy of EP/15% F-Lignin@HKUST-1 at different conversions.
Table 3. Activation energy of EP/15% F-Lignin@HKUST-1 at different conversions.
αKissiger-SY IterationOzawa-SY IterationLee-Beck IterationGorbatchev Iteration
E/(KJ/mol)R2E/(KJ/mol)R2E/(KJ/mol)R2E/(KJ/mol)R2
0.15133.400.95133.400.94133.420.93133.420.95
0.20129.260.96129.260.96129.290.96129.290.95
0.25129.840.97129.840.97128.910.98129.860.97
0.30128.890.98128.890.99128.220.99128.910.98
0.35128.190.99128.190.98128.970.99128.220.99
0.40128.940.99128.940.99129.460.99128.970.99
0.45129.440.99129.440.991 29.950.99129.460.99
0.50129.930.99129.930.99131.360.99129.950.99
0.55131.330.99131.330.99132.900.99131.360.99
0.60132.870.98132.870.99134.600.99132.900.96
0.65134.570.98134.570.97136.700.98134.600.97
0.70136.670.95136.670.94137.780.96136.700.95
0.75137.760.95137.760.92138.500.96137.780.95
0.80138.480.96138.480.91140.780.94138.500.94
0.85140.750.95140.750.92144.140.95140.780.93
0.90144.110.95144.1 10.91146.390.94144.140.94
Table 4. Kinetic parameters of universal integral equations and differential equations based on the Avrami-Erofeev equation at different temperature rates of Al-lignin.
Table 4. Kinetic parameters of universal integral equations and differential equations based on the Avrami-Erofeev equation at different temperature rates of Al-lignin.
Mechanism Functionβ
/K∙min−1
Universal IntegrationDifferential Equation
E/
KJ∙mol−1
Ln (A/min−1)R2E/
KJ∙mol−1
Ln (A/min−1)R2
AE4865.275.970.9771.626.790.97
1274.167.320.9775.817.380.96
1671.877.190.9876.897.770.97
2075.487.770.9672.497.220.95
AE5884.268.90.9890.619.730.98
1295.3610.490.9897.0410.550.97
1692.5010.250.9897.5410.840.98
2097.0410.930.9794.1010.370.96
AE68103.2511.820.98109.612.670.98
12116.5713.650.98118.2713.720.97
16113.1313.310.98118.1913.910.98
20118.6114.080.97115.713.530.96
Table 5. Kinetic parameters of universal integral equations and differential equations based on the Avrami-Erofeev equation at different temperature rates of EP/15% F-lignin@HKUST-1.
Table 5. Kinetic parameters of universal integral equations and differential equations based on the Avrami-Erofeev equation at different temperature rates of EP/15% F-lignin@HKUST-1.
Mechanism Functionβ
/K∙min−1
Universal IntegrationDifferential Equation
E/
KJ∙mol−1
Ln (A/min−1)R2E/
KJ∙mol−1
Ln (A/min−1)R2
AE4573.968.970.9272.4710.370.90
1095.4111.840.9386.8812.600.91
1581.449.990.9580.0712.080.94
20112.1814.860.94101.8615.730.92
AE55116.8516.310.93115.2217.680.92
10149.7320.740.94140.9321.420.93
15127.6617.570.96126.2419.640.95
20174.2124.810.95163.8125.630.94
AE65159.7423.660.93158.0425.010.93
10204.0529.630.94195.0830.270.94
15173.8925 150.96172 4327 220.96
20236.2434.750.95225.7835.550.94
Table 6. Activation energy and preintential factors of Al-lignin at different conversions under Kissinger-SY iteration and Ozawa-SY iteration method.
Table 6. Activation energy and preintential factors of Al-lignin at different conversions under Kissinger-SY iteration and Ozawa-SY iteration method.
αKissinger-SY IterationOzawa-SY Iteration
E/KJ∙mol−1lnA/min−1E/KJ∙mol−1lnA/min−1
0.1557.610.9257.610.92
0.2092.199.8392.19.9.83
0.25112.6115.02112.6215.02
0.30122.6417.74122.6417.74
0.35126.619.07126.619.07
0.40127.8919.75127.8919.75
0.45127.7320.06127.7320.06
0.50123.2819.43123.2819.43
0.55111.5417.22111.5417.22
0.60100.4114.95100.4214.95
0.6564.248.5364.248.53
0.7049.345.8449.345.84
0.75129.1216.92129.1216.92
0.80148.419.43148.4119.43
0.85124.0216.63124.0216.63
0.9071.8710.8071.8710.80
Table 7. Activation energy and preintential factors of EP/15% F-lignin@HKUST-1 at different conversions under Kissinger-SY iteration and Ozawa-SY iteration method.
Table 7. Activation energy and preintential factors of EP/15% F-lignin@HKUST-1 at different conversions under Kissinger-SY iteration and Ozawa-SY iteration method.
αKissinger-SY IterationOzawa-SY Iteration
E/KJ∙mol−1lnA/min−1E/KJ∙mol−1lnA/min−1
0.15129.8419.38129.8419.38
0.20128.8919.97128.8919.97
0.25128.1920.44128.1920.44
0.30128.9421.08128.9421.08
0.35129.4421.59129.4421.59
0.40129.9322.05129.9322.05
0.45131.3322.64131.3322.64
0.50132.8723.22132.8723.22
0.55134.5723.80134.5723.80
0.60136.6724.42136.6724.42
0.65137.76.24.81137.76.24.81
0.701384825.091384825.09
0.75140.7525.60140.7525.60
0.80144.1126.20144.1126.20
0.85146.3626.42146.3626.42
0.90142.2625.32142.2625.32
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

Yao, T.; Yang, R.; Sun, C.; Lin, Y.; Liu, R.; Yang, H.; Chen, J.; Gu, X. Pyrolysis Kinetics of Lignin-Based Flame Retardants Containing MOFs Structure for Epoxy Resins. Molecules 2023, 28, 2699. https://doi.org/10.3390/molecules28062699

AMA Style

Yao T, Yang R, Sun C, Lin Y, Liu R, Yang H, Chen J, Gu X. Pyrolysis Kinetics of Lignin-Based Flame Retardants Containing MOFs Structure for Epoxy Resins. Molecules. 2023; 28(6):2699. https://doi.org/10.3390/molecules28062699

Chicago/Turabian Style

Yao, Tianyu, Ruohan Yang, Cong Sun, Yuzhu Lin, Ruoqi Liu, Hongyu Yang, Jiajia Chen, and Xiaoli Gu. 2023. "Pyrolysis Kinetics of Lignin-Based Flame Retardants Containing MOFs Structure for Epoxy Resins" Molecules 28, no. 6: 2699. https://doi.org/10.3390/molecules28062699

APA Style

Yao, T., Yang, R., Sun, C., Lin, Y., Liu, R., Yang, H., Chen, J., & Gu, X. (2023). Pyrolysis Kinetics of Lignin-Based Flame Retardants Containing MOFs Structure for Epoxy Resins. Molecules, 28(6), 2699. https://doi.org/10.3390/molecules28062699

Article Metrics

Back to TopTop