Next Article in Journal
Grafting of Anionic Decahydro-Closo-Decaborate Clusters on Keggin and Dawson-Type Polyoxometalates: Syntheses, Studies in Solution, DFT Calculations and Electrochemical Properties
Next Article in Special Issue
Enhanced Adsorption of Sulfonamides by Attapulgite-Doped Biochar Prepared with Calcination
Previous Article in Journal
Structural Dynamics of Chloromethanes through Computational Spectroscopy: Combining INS and DFT
Previous Article in Special Issue
Oak Biomass in the Form of Wood, Bark, Brushwood, Leaves and Acorns in the Production Process of Multifunctional Biochar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermogravimetric Analysis and Kinetic Modeling of the AAEM-Catalyzed Pyrolysis of Woody Biomass

1
Department of Mechanical Engineering, École de Technologie Supérieure, Montreal, QC H3C 1K3, Canada
2
Centre de Recherche de Royallieu, Université de Technologie de Compiègne, EA 4297-TIMR, BP20529, 60205 Compiègne, France
3
École Supérieure de Chimie Organique et Minérale, 1 Rue du Réseau Jean-Marie Buckmaster, 60200 Compiègne, France
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(22), 7662; https://doi.org/10.3390/molecules27227662
Submission received: 23 October 2022 / Revised: 1 November 2022 / Accepted: 1 November 2022 / Published: 8 November 2022
(This article belongs to the Collection Recycling of Biomass Resources: Biofuels and Biochemicals)

Abstract

:
This work analyzes the catalytic effects induced by alkali and alkaline earth metals (AAEMs) on pyrolysis kinetics. To this end, thermogravimetric analyses (TGA) were carried out with raw beech wood and samples impregnated with NaCl, KCl and MgCl2 at four heating rates (5, 10, 15 and 30 °C/min). Obtained results showed that AAEM compounds promote the decomposition of biomass by reducing the initial and peak pyrolysis temperatures. More specifically, the catalytic effect of the alkaline earth metal was shown to be stronger than that of alkali metals. To further interpret the obtained trends, a kinetic modeling of measured data was realized using two isoconversional methods (the Ozawa–Flynn–Wall (OFW) and Kissinger–Akahira–Sunose (KAS) models). With a view to identifying a suitable reaction model, model fitting and master plot methods were considered to be coupled with the isoconversional modeling approaches. The 3-D diffusion reaction model has been identified as being well suited to properly simulate the evolution of the conversion degree of each sample as a function of the temperature. Furthermore, the kinetic parameters derived from the present modeling work highlighted significant decreases of the activation energies when impregnating wood with AAEM chlorides, thus corroborating the existence of catalytic effects shifting the decomposition process to lower temperatures. A survey of the speculated pathways allowing to account for the impact of AAEMs on the thermal degradation of woody biomass is eventually proposed to better interpret the trends identified in this work.

1. Introduction

Alongside increasing concerns about global warming, biomass, a neutral-carbon fuel, is attracting more and more attention, not only for heat and electricity generation, but also for the production of valuable biochemicals and biomaterials. When being pyrolyzed under an inert atmosphere, the biopolymers composing the lignocellulosic biomass typically decompose into biochar, bio-oil and incondensable biogas [1]. However, the yields and selectivity towards target products (e.g., bio-oil) remain limited. Furthermore, and despite the major progress that has been achieved in the field of biomass pyrolysis, bio-oils issued from the thermal conversion of raw biomass generally suffer from a high oxygen content, viscosity and corrosiveness, together with a low heating value [2]. To overcome these limitations, the implementation of a catalytic treatment represents an interesting option. It indeed allows upgrading obtained pyrolysis products by reducing the yields of undesirable components while improving the properties of the target products [3,4].
Among the catalysts currently considered in the framework of biomass pyrolysis, specific attention has recently been devoted to alkali and alkaline earth metals (AAEMs) due to their low toxicity, affordable price, high catalytic efficiency and non-negligible quantity in raw biomass [5]. AAEM compounds can be directly mixed with biomass through wet impregnation (in situ configuration) or brought in downstream in order to react with the vapors produced during pyrolysis (ex situ configuration) [3]. These convenient catalysts can, moreover, be implemented in a wide range of reactor configurations, including fixed bed [6,7] and fluidized bed [8,9] pyrolyzers. To date, most of the studies undertaken to elucidate the impact of AAEMs on the catalytic conversion of biomass have focused on the analysis of the reaction products to evaluate the efficiency of the different alkali and alkaline earth metals in promoting the formation of some target molecules. For instance, Peng et al. investigated the effect of different alkaline additives (NaOH, KOH, Na2CO3 and K2CO3) on the production of phenols during the pyrolysis of lignin [10]. They notably observed that all alkalis promote decarboxylation or decarbonylation reactions as well as the removal of unsaturated alkyl-branched chains from aromatic molecules [10]. Furthermore, they noted that the strong hydroxide alkalis of NaOH and KOH favor the deoxygenation of methoxy groups, thus leading to phenols free of methoxy groups. Zhang et al. then studied the impact of impregnating camphor branch, corn cob and walnut shell with potassium nitrate on the distribution of bio-oils produced from fast pyrolysis in a fluidized bed [9]. They reported that potassium tends to increase the yields of furans and phenols while decreasing the yields of aldehydes esters and sugars. Reaction pathways involving dehydration, ring scission, depolymerization and cracking processes were especially proposed to account for the influence of potassium on biomass pyrolysis. As far as alkaline earth metals are concerned, Veses et al. also noted that calcium-based materials (such as CaO) tend to promote dehydration reactions, thereby decreasing the acidity and oxygen content of the bio-oils issued from the pyrolysis of forest pine wood while increasing their pH and calorific value [11]. As for divalent magnesium cations contained in MgCl2, they have proven to enhance the degradation of hemicellulose to form oxygenate molecules and furans [12] while promoting repolymerization reactions, leading to increased char formation and higher molecular weight compounds in bio-oils [13].
Notwithstanding the findings reported in the above literature survey, kinetic analyses dealing with the AAEM-catalyzed pyrolysis of biomass remain relatively rare, as highlighted in [5]. This lack is all the more critical since pyrolysis, as the first step in the thermochemical conversion of biomass, directly influences the nature and the distribution of the released products, and hence, the subsequent reaction stages, as well as the overall fuel conversion rate [1]. Characterizing pyrolysis kinetics is therefore important for the proper design of reactors and for the optimization of industrial facilities. Among existing models, one can cite the model-free methods [14], the distributed activation energy and the chemical percolation devolatilization models (see [15] and references therein) in addition to simulation tools relying on the density functional theory for instance [16]. These models can be roughly classified into two categories, depending on whether they aim at simulating either the mass loss rate of the fuel, the distribution of the pyrolytic products or both. In the context of TGA-based kinetic analyses, model-free methods remain extensively used in assessing kinetic parameters as they represent an accurate route to directly infer activation energies without the need for any initial assumption regarding the reaction model [17].
With this in mind, the present work aims at analyzing and modeling the kinetics related to the AAEM-catalyzed pyrolysis of a woody biomass. To this end, pyrolysis experiments were conducted with raw beechwood and samples catalyzed with three AAEM compounds (NaCl, KCl and MgCl2) added by wet impregnation. Non-isothermal thermogravimetric measurements were performed with heating rates ranging from 5 to 30 °C/min. Obtained results were then modeled by means of two isoconversional models (namely, the Ozawa–Flynn–Wall (OFW) and Kissinger–Akahira–Sunose (KAS) models) to ease the interpretation of obtained trends.
After introducing the experimental methodology together with the theoretical background underlying the OFW and KAS modeling approaches in Section 2, measured data will be presented and commented on in Section 3. Results issued from the kinetic modeling calculations will notably be detailed therein before being compared with experimental results. In a bid to better interpret observed trends as far as the impact of AAEMs is concerned, Section 4 will finally propose a summary of the main mechanisms at play during the AAEM-catalyzed pyrolysis of biomass.

2. Methodology

2.1. Feedstocks and Sample Preparation

The proximate and ultimate analyses of the beech wood used herein are provided in Table 1.
Three AAEM compounds, namely, NaCl, KCl and MgCl2, were selected to be added to beech wood. Samples were prepared by wet impregnation to favor the interaction between the cations and the biomass (see Figure 1). Beech wood samples were thoroughly impregnated in solutions containing exactly the same amount of cation for each AAEM chloride. 4 g of biomass was typically put into prepared NaCl, KCl and MgCl2 solutions whose concentration of metal cations had been set to 5 g/L for each catalyst. The mixtures were stirred by means of a magnetic stirrer for two hours, as in [18]. Impregnated samples were then filtered to eliminate the extra cations and anions remaining in the solution. As for the control sample consisting of the raw biomass, it was suspended in deionized water for the same duration and within the same stirring conditions to exclude the effects of water washing. Finally, all the samples were dried in an oven at 105 °C [19,20,21] for 24 h to get rid of the extra free water.

2.2. Experimental Apparatus and Procedure

Non-isothermal pyrolysis tests were performed using a SETARAM SETSYS Evolution thermogravimetric analyzer (TGA). Four heating rates (5, 10, 15 and 30 °C/min) were used to heat the samples from room temperature up to 950 °C, with a plateau at 105 °C, for 20 min to ensure the elimination of free water. To perform all the TGA tests, 10 mg of sample was put in alumina crucibles. An inert atmosphere was continuously maintained around heated samples by means of a 100 mL/min nitrogen flow. Three tests were performed for each sample and operating condition. It should therefore be considered that all the mass loss profiles presented in the following and used for the modeling work presented in Section 3 are based on averaged data. Following [22,23], only the section below 700 °C (which represents the main part of the mass loss process) was taken into account for the calculations presented hereafter to get rid of the possible measurement noise encompassing the data recorded at the highest temperatures. To conclude, the conversion degree ( α ) at any given time (t) was calculated from the initial (i) and final (f) residual masses (noted ‘TG’ and expressed in wt%) based on Equation (1):
α = TG i TG t TG i TG f
where the conversion degrees are set to 0% and 100% for the initial (106 °C) and final (700 °C) measurement points, respectively.

2.3. Kinetic Modeling

2.3.1. Isoconversional Models

The variation of the fuel conversion degree α as a function of the temperature T (expressed in K) during TGA experiments can be expressed using an Arrhenius equation of the type:
d α dT = A β × exp E a R × T × f α
where A is the pre-exponential factor (expressed in s−1), β is the heating rate, E a stands for the activation energy (in kJ/mol), R is the gas constant (equal to 8.314 J/(mol·K)), while f α denotes the reaction model. By integrating both sides of Equation (2) while assuming that the initial temperature and conversion degree are equal to zero, the integral form of the reaction model g α can be expressed as follows:
g α = 0 α d α f α = A β × T 0 T exp E a R × T dT A β × 0 T exp E a R × T dT
The term on the right-hand side of Equation (3) corresponds to an exponential integral which has no exact solution (in closed form). An algebraic approximation of this term must therefore be used to enable the calculation [24,25]. By replacing dT by du in the integral (with u = E a R × R ), the right-hand term of Equation (3) can thus be converted into the form:
g α = A β × 0 T exp E a R × T dT = A β × u ( e u × E a R × u 2 ) du = A × E a β × R × u + ( e u u 2 ) du = A × E a β × R × p u
such that:
ln g α = ln ( A × E a β × R ) + ln p u
where p u is an exponential integral that does not have an analytical solution. Based on Doyle’s assumption [26], the function p u can be converted into a series allowing deriving the expressions of the OFW and KAS models detailed below.
As for the OFW integral isoconversional model [27,28], it consists of plotting the evolution of ln β i as a function of 1 T for different heating rates β i to obtain linearized straight lines whose slopes allow deriving activation energy values for each conversion degree α (see Equation (6)).
ln β i = ln A α × E a , α R × g α 5.331 1.052 × E a , α R × T
Regarding the values of A α , they can be inferred using Equation (7) once the reaction model g α and the activation energies are determined [27]:
A α = R × exp b + 5.331 × g α E a , α
where b is the intercept of the linearized straight lines.
As far as the KAS model is concerned, it consists of plotting straight lines depicting the evolution of ln β i T 2 as a function of 1 T for different heating rates β i . One can then infer the values of the activation energy ( E a , α ) and pre-exponential factor ( A α ) for any conversion degree α by means of the slopes and intercepts of the so-obtained straight lines (see Equations (8) and (9), respectively).
ln β i T 2 = ln A α × R E a , α × g α E a , α R × T
A α = E a , α × exp b × g α R

2.3.2. Methodology Allowing to Identify a Suitable Reaction Model

With a view to selecting a proper reaction model, two possible routes are summarized in the review by Wang et al. [29]. The first solution relies on the selection of the reaction model by comparing the values of the activation energy derived from model-free (i.e., isoconversional) and model fitting methods (see the ‘Model Fitting Approach’ subsection below). Assuming that the isoconversional methods allow assessing more accurate activation energies, these E a values should therefore be considered in the model fitting procedure in order to select an appropriate reaction model allowing to derive similar activation energies [30,31,32,33]. The second route concerns the use of the master plot method (see the ‘Master Plot’ subsection), which allows identifying the most suitable reaction model by comparing experimental curves to some pre-established theoretical ones [33,34,35,36].

Model Fitting Approach

The data fitting method selected herein consists of the Šatava–Šesták model [37], which is commonly used in studies dealing with the thermal decomposition of solid materials [38,39,40,41,42,43]. By plotting the evolution of log g α as a function of 1 T for different g α formulations, straight lines can be obtained (see Equation (10)) with their slopes allowing to derive activation energy values. These latter can then be used to identify a proper reaction model that allows obtaining E a values matching those assessed by means of model-free methods. Finally, the log A × E a β × R term enables inferring the value of the pre-exponential factor A :
log g α = log A × E a β × R 2.315 0.4567 × E a R × T
In the present work, 10 different reaction models were tested, namely, F1, F2, F3, D2, D3, R2, R3, A2, A3 and A4 [44,45].

Master Plot

To ease the identification of a suitable reaction model, Sánchez-Jiménez et al. proposed a calculation procedure consisting of transforming measured data into an experimental master plot that does not depend on experimental conditions [46]. This master plot is then compared with theoretical master plots obtained by using different reaction models. By means of a simple graphical procedure, the most suited kinetic model can be identified. In this process, one must first introduce a so-called generalized time, θ , whose derivative over t can be written as:
d θ dt = exp E a R × T
By combining the general expression of d α /dt (i.e.,   d α dt = A × exp E a R × T × f α ) with d θ /dt (see Equation (11)), d α /d θ can be expressed in the form:
d α d θ = A × exp E a R × T × f α exp E a R × T = A × f α
For a single-step process, the expression of the reaction model is invariable. Using a reference point at α = 50%, one obtains:
d α d θ α d α d θ 50 % theo = f α f 50 %
where the generalized reaction rate can be expressed as follows:
d α d θ = d α dt × exp E a R × T
so that the relationship between the generalized reaction rate and the experimental data are expressed as:
d α d θ α d α d θ 50 % exp = d α dt α d α dt 50 % × exp E a R × T exp E a R × T 50 %
By comparing theoretical and experimental curves (whose plotting requires the estimation of E a through isoconversional methods under non-isothermal conditions), the f α formulation leading to the best match between calculated and measured data can be identified.

3. Results and Discussion

Thermogravimetric analysis results will first be presented and discussed in Section 3.1 before being modeled by means of the OFW and KAS isoconversional approaches in Section 3.2.

3.1. Experimental Characterization of the Pyrolysis Behavior of Raw and Impregnated Beech Wood Samples

Results issued from the TGA analyses conducted with raw and impregnated beech wood samples are detailed in Figure 2. Curves depicting variations of the mass loss (TG) and mass loss rate (dTG) as a function of the temperature are notably reported therein for heating rates of 5, 10, 15 and 30 °C/min. As explained in Section 2.2, three tests were performed for each sample and operating condition. Error bars plotted in Figure 2 consequently account for the dispersion of the experimental data around the mean values.
Based on obtained TG and dTG curves, it can first be noted that the pyrolysis behaviors of the four tested samples significantly diverge from one another regardless of the heating rate. More particularly, the temperatures at which pyrolysis takes place appear to be higher for the raw biomass than for the impregnated samples.
To better highlight and discuss the observed discrepancies, Table 2 reports the characteristic pyrolysis temperatures. This includes the initial (Ti) and final (Tf) decomposition temperatures estimated for conversion degrees of 10 and 90%, respectively, as well as the peak temperature (Tp) assessed at the peak mass loss rate.
Based on the data gathered in Figure 2 and Table 2, it can be noted that the beech wood starts to decompose at temperatures between 257.5 and 286.4 °C for 5 °C/min < β < 30 °C/min.
As shown in Figure 2b,d,f,h, the dTG curves reach a peak for temperatures between 340 and 361 °C (corresponding to the decomposition of cellulose), with shoulders located on the left for temperatures of around 250–300 °C (denoting the degradation of hemicellulose) [47]. The mass loss rates then decrease to plateau above 412–432 °C, depending on the considered heating rate. As far as impregnated samples are concerned, plotted TG and dTG curves exhibit more or less reduced Ti and Tp values, which is consistent with the role played by alkali and alkaline earth metals on the reduction of the pyrolysis temperatures [13,20,48,49,50,51,52,53]. AAEM cations are indeed likely to cleave chemical bonds during the impregnation process, and hence favor the decomposition of the biomass structure, especially at the beginning of the pyrolysis process (i.e., for relatively low temperatures). It is interesting here to note that the tested alkaline earth metal (magnesium) seems to exhibit a higher catalytic effect than the two alkali metals (sodium and potassium), as exemplified by the greater Ti reductions observed with MgCl2 regardless of the heating rate. Impregnating beech wood with magnesium chloride indeed allows reducing the initial pyrolysis temperatures by 43 °C on average when considering all β values as compared to the Ti values measured for the raw biomass sample. On the other hand, potassium only induces a mean Ti drop of 10 °C against 1.6 °C for sodium, which mainly shows a catalytic effect for temperatures higher than 300 °C.
It is noteworthy that these trends are actually in line with those identified in previous investigations dealing with the impact of AAEMs on the pyrolysis of woody biomass components. Yu et al. indeed reported that magnesium was able to lower the onset temperature of cellulose pyrolysis [51]. Gao et al. and Leng et al. also reported that Ti values were much more reduced when using alkaline earth metals such as calcium and magnesium instead of potassium to catalyze the pyrolysis of rice straw and cellulose, respectively [52,53].
As for the peak temperatures, similar trends are observed for all catalysts. Tp indeed decreases by 22, 29 and 15 °C on average for the samples impregnated with NaCl, KCl and MgCl2, respectively (see Table 2). It is interesting to observe that the effect of sodium and potassium on Tp is much more significant than on Ti. Furthermore, the reduction of Tp induced by the use of magnesium chloride turns out to be the lowest, which contrasts with the fact that such a catalyst was previously shown to induce the highest Ti drop. Looking at the residual masses measured at 700 °C (noted ‘TG700°C‘ in Table 2), it can be noted that these latter increase (from +3.7 for MgCl2 to +8.2 wt% for NaCl on average) when impregnating beech wood with AAEMs. Mean TG700°C values calculated based on all heating rates indeed pass from 30.4 wt% for the raw biomass to 38.6 wt%, 37.7 wt% and 34.1 wt% when using sodium-, potassium- and magnesium-based catalysts, respectively. Such a trend, which here again, is in line with the observations made in previous studies [13,21,49,51,54,55,56,57], can be traced to an enhanced char formation and/or to recombination reactions occurring at the last stage of the pyrolysis process when the temperature exceeds 350 °C. All the results presented within this section will, however, be further discussed in Section 4, which deals with the description of the main mechanisms underlying the AAEM-catalyzed pyrolysis of biomass components.

3.2. Kinetic Analysis

As illustrated in Section 2.3.1, isoconversional models allow inferring E a values without requiring any a priori assumption regarding the selection of the reaction model. As such, these methods generally lead to the derivation of consistent activation energies [29,58], which can be used as inputs in subsequent calculation stages to identify suitable reaction models and assess proper values of pre-exponential factors. This type of calculation procedure implemented in [30] has since been widely used and/or recommended for the modeling of data issued from non-isothermal pyrolysis experiments [29,31,38,41,59,60,61,62]. Following a similar approach, OFW and KAS models were thus implemented herein (see Section 3.2.1) before being coupled with Šatava–Šesták and master plot calculation procedures (see Section 3.2.2).

3.2.1. Isoconversional Modeling of TGA Results

The OFW and KAS isoconversional models were implemented following the calculation procedures described in Section 2.3.1. Doing so led to the obtention of the linearized straight lines depicted in Figure 3 and Figure 4, from which the values of the α-dependent activation energy reported in Table 3 were inferred.
As previously mentioned (see Section 2.2), results issued from three tests were averaged for each sample and operating condition to mitigate the deviations possibly observed from test to test due to measurement noise and uncertainties. As such, it was possible to obtain very good linear correlations, as shown in Figure 3 and Figure 4, and evidenced by the high coefficients of determination reported in Table 3 (R2 ≈ 0.985 on average). As can be seen, for each 10% step of the conversion degree, identical activation energies are obtained with the OFW and KAS models (see Table 3). The mean relative standard deviation between the E a values assessed using both methods is indeed of only 1.51%. Furthermore, and overall, the higher the conversion degree, the higher the activation energy, as also noted in [58,63,64,65,66,67], which is consistent with what may have been expected since the species emitted at high temperatures typically require more energy to be released.
When analyzing the results gathered in Table 3 in greater detail, it can be noted that the activation energies assessed for raw biomass with the OFW and KAS models are quite constant for 20% ≤ α ≤ 80%, with values ranging from 153.0 to 160.5 kJ/mol with the OFW model and from 150.3 to 159 kJ/mol with the KAS one. On the other hand, the activation energy increases significantly, up to around 310 kJ/mol, when the conversion degree reaches 90%, which may be related to the decomposition of lignin that contains more rigid carbon-carbon linkages requiring more energy to be cleaved [29]. As for the samples impregnated with NaCl, KCl and MgCl2, the range of E a values estimated for 20% ≤ α ≤ 70% is also relatively narrow (i.e., between ~132 and ~140 kJ/mol for NaCl and KCl compared to values between ~145 and ~161 kJ/mol for MgCl2). The activation energies then increase for higher conversion degree, thus indicating a pyrolysis behavior relatively similar to that described in the case of beech wood. Another interesting feature emerging from the tests performed with NaCl, KCl and MgCl2 concerns the global reduction of the E a values resulting from the impregnation of beech wood with these AAEM-catalysts. Mean decreases of the activation energies between ~5 and ~19 kJ/mol are indeed observed for 10% ≤ α ≤ 70%. Since the activation energy accounts for the minimum energy required for a reaction to occur, any decrease of E a tends to indicate that less external energy is required to overcome the energy barrier allowing the pyrolysis to take place. As a consequence, the data reported in Table 3 tend to corroborate one of the experimental observations made in Section 3.1 regarding the catalytic effect induced by AAEM compounds, which are able to shift the decomposition of biomass to lower temperatures. Lastly, it is noteworthy that the E a values inferred herein are consistent with those reported in the literature for widely varied biomass sources with values which are generally between 100 and 300 kJ/mol [29]. Nevertheless, interpreting the kinetic features related to the AAEM-catalyzed pyrolysis of beech wood in greater detail requires conducting additional analyses with the view to identifying a suitable reaction model, as detailed in Section 3.2.2.

3.2.2. Identification of a Proper Reaction Model

Selection of a Proper g α Formulation through Model-Fitting Calculations

The Šatava–Šesták model was implemented following the procedure described in Section 2.3.2. Curves depicting the evolution of log g α as a function of 1 T were plotted for 10 reaction models commonly used in the literature. To this end, only the data associated with conversion degrees between 10 and 80% (i.e., those free from potentially important measurement noise) were considered to enable consistent data fitting calculations. Among tested reaction models, the focus was placed on order-based (Fn), diffusion (Dn), geometrical contraction (Rn) and nucleation (An) ones, as listed in Section 2.3.2 [44,45,64]. An example of data fitting based on the results obtained for a heating rate of 10 °C/min with the D3 reaction model is presented in Figure 5. By repeating the same procedure for the 10 reaction models considered herein (i.e., F1, F2, F3, D2, D3, R2, R3, A2, A3 and A4), a series of linearized curves were obtained, with their slopes allowing to infer the values of the apparent activation energy.
Although most of the tested models allow obtaining a high coefficient of correlation (0.93 ≤ R2 ≤ 0.99), identifying a proper model based solely on R2 values may be inadequate, as highlighted by Khawam and Flanagan [30]. According to these authors, the selection of an adapted reaction model should be preferentially achieved based on isoconversional model plots (IMO) of activation energies as a function of α . Following this approach, E a values derived from the OFW and KAS models were compared with those inferred by means of the Šatava–Šesták model integrating the 10 mechanisms listed above. To this end, linear data fittings were performed at each investigated heating rate with all the g α formulations (data not reported herein for brevity). Obtained results then showed that β had no significant influence on the so-inferred E a values. The relative standard deviation between the activation energies assessed at 5, 10, 15 and 30 °C/min was indeed found to be less than 0.6% on average for all the reaction models in the case of the samples impregnated with NaCl, as an example. Mean activation energies were therefore considered, as in [38], and plotted in Figure 6 for comparison purposes. The best suited reaction model can then be identified, and corresponds to that whose implementation in the model fitting procedure allows assessing E a values closest to those derived when using the model-free approaches [30,38].
Obtained plots thus show that order-based (F3) and diffusion (D2 and D3) models allow obtaining activation energies quite close to those inferred using the OFW and KAS models. Nevertheless, the D3 model still appears to be the most suited. Its use indeed leads to the lowest deviations between the E a values inferred by means of the model fitting and isoconversional approaches (~4% compared to 4.7 and 8.8% for the D2 and F3 models, respectively, based on the results reported in Figure 6). Nevertheless, and before drawing a related clear-cut conclusion, the suitability of the D3 model will be verified.

Validation of the Selected g α Formulation through Direct Calculations

To corroborate the conclusion drawn above, we carried out complementary calculations aimed at stimulating the evolution of the conversion degree of each sample as a function of the temperature when considering the 10 reaction models used herein. Obtaining such curves, however, required integrating both E a and   A in Equation (2). Pre-exponential factors were therefore assessed in addition to the activation energies listed in Table 3. To that end, the ten g α formulations were integrated in Equations (7) and (9) to perform the simulations based on the OFW and KAS models, respectively. Proceeding as such allowed to derive the pre-exponential factors listed in Tables S1–S4, which are provided as Supplementary Materials to this article (see the Supplementary Materials). Introducing these parameters in Equation (2) then led to the obtention of the kinetic profiles depicted in Figure 7, which display the evolution of α as a function of T calculated based on the KAS model for a heating rate of 10 °C/min (data related to the OFW model are reported in Figure S1 of the Supplementary Materials, for their part). To obtain these conversion degree profiles, the temperature values reported in Table 2 for an α of 10% were used as initial reaction temperatures. Furthermore, calculations were carried out by applying the kinetic parameters derived for a given α on a ±5% conversion degree range (e.g., rate constant parameters estimated for α = 20% were kept constant to perform the calculations on the 15% < α < 25% range).
As can be seen by looking at Figure 7, all tested models globally reproduce the overall shape of the measured conversion degree profiles. Some mechanisms (e.g., F2, F3, D2), however, exhibit better predictive capabilities, as is especially the case of the D3 model previously identified as being the most suited. Such a direct calculation procedure thus tends to confirm the relevance of the 3-D diffusion reaction model, as further verified through master plot calculations below.

Validation of the Selected g α Formulation through Master Plot Calculations

Following the procedure described in Section 2.3.2, the generalized master plot method proposed by Sánchez-Jiménez et al. was implemented on a range of α for which the rate constant parameters were quite constant. Indeed, and as reiterated in Section 2.3.2, the master plot method is only applicable to single-step processes for which the activation energy does not vary with the conversion degree. As illustrated in Table 3, such a condition is fulfilled when 10% < α < 80% for the different tested samples, which corresponds to the main pyrolysis stage. Measured data were thus transformed into experimental master plots which were compared with theoretical master plots derived from the use of the F1, F2, F3, D2, D3, R2, R3, A2, A3 and A4 reaction models.
Figure 8 depicts the results obtained with the samples impregnated with NaCl and KCl, as examples for a heating rate of 10 °C/min. Results depicted therein confirm that the D3 reaction model appears to be the most suited to match experimental points (similar trends being obtained with the other samples). This conclusion is actually consistent with the fact that diffusion models are generally recommended in kinetic analyses aimed at investigating biomass pyrolysis [63,64,67,68,69]. This, moreover, confirms the interest of applying the master plot method to select a proper reaction model, as highlighted in a recent review by Wang et al. focusing on biomass pyrolysis [29].

3.2.3. Simulation of Conversion Degree Profiles

By integrating the parameters assessed by means of the OFW and KAS models in Equation (2) while considering the so-identified D3 mechanism, one obtains the conversion degree profiles depicted in Figure 9 for a heating rate of 10 °C/min, as an example (profiles related to β values of 5, 15 and 30 °C/min are reported in Figures S2–S4 of the Supplementary Materials).
As can be seen, simulated data satisfactorily reproduce their experimental counterparts regardless of the considered sample. As such, the results depicted in Figure 9 tend to corroborate the consistency of the rate constant parameters and reaction mechanism assessed through this modeling work. To conclude the analysis of the catalytic effect induced by AAEMs on the pyrolysis kinetics of beechwood, the rate constant ( k = A × exp E a / ( R × T ) ) of each sample was calculated for temperatures between 250 and 400 °C, which corresponds to the main pyrolysis stage. Here again, the D3 model was considered together with mean apparent rate constant parameters assessed by means of the KAS model for 10% < α < 80% (similar trends are obtained for the OFW model). It is then noteworthy that the rate constant values are higher for impregnated samples. As an example, the rate constants calculated for pure wood and samples impregnated with NaCl, KCl and MgCl2 reach values of 1.83 × 10−4 s−1, 2.92 × 10−4 s−1, 3.83 × 10−4 s−1 and 4.26 × 10−3 s−1, respectively, for a temperature of 325 °C (i.e., 598.15 K). Furthermore, the ability of the tested AAEM chlorides to promote the pyrolysis follows the sequence: MgCl2 > KCl > NaCl, which is in line with the experimental observations made in Section 3.1. The modeling work realized herein thus allowed deriving kinetic parameters suitable for experimentally reproducing monitored conversion degree profiles while corroborating, from a kinetic perspective, the higher catalytic efficiency of the tested alkaline earth metal over alkali metals. The mechanisms likely to explain such trends will, however, be discussed further in detail in Section 4.

4. Summary of the Mechanisms Underlying the Catalytic Impact of AAEMs on Biomass Pyrolysis

To better understand the observations made in Section 3 regarding the effects of alkali and alkaline earth metals on the thermal degradation of beech wood, this section aims at reviewing the major mechanisms likely to influence the AAEM-catalyzed pyrolysis of biomass. Actually, there are essentially two routes allowing to account for the catalytic effects induced by AAEMs [1,5].
First, alkali and alkaline earth metals have been demonstrated to promote the cleavage of intra- and inter-molecular bonds (i.e., glycosidic and hydrogen bonds) in biomass. They, moreover, favor the cracking of primary volatiles, thus enhancing the yields of light oxygenated compounds and incondensable gaseous species while contributing to limiting the formation of levoglucosan. In that respect, Zhao and Li notably showed that small sodium cations were likely to pierce through rice husk textures to act on biopolymers (i.e., cellulose, hemicellulose and lignin), thus changing the pyrolysis reaction pathways and favoring biomass degradation [70]. More specifically, Na+ cations have been proven to act on cellulose through hemolytic cleavage in pyranose rings and heterocyclic cleavage of glycosidic linkage while promoting ring scission, isomerization, dehydration, decarbonylation and/or decarboxylation reactions, which promote the formation of smaller furans and alcohols [48,70]. Besides, sodium cations can also act on the branched-chain structure of hemicellulose to ease depolymerization, ring scission, dehydration as well as rearrangement reactions. Finally, NaCl can enhance the decomposition of lignin by facilitating the chemical bonds of the lignin structure to undergo dehydration, protonation, aromatization and rearrangement reactions. All these processes therefore contribute to favor biomass degradation, which is consistent with the experimental observations and kinetic analyses made throughout the present work. As far as potassium chloride is concerned, Leng et al. showed that such a catalyst was likely to promote the cleavage of glycosidic bonds and the scission of pyran rings directly by homolytic reaction to form low molecular weight species during the fast pyrolysis of cellulose [53]. Furthermore, Safar et al. noted that the amount of intermolecular bonds of cellulose was reduced when impregnating woody biomass with potassium cations, hence decreasing the crystallinity of this biopolymer and increasing the biomass reactivity during pyrolysis [21]. As for magnesium, it has been shown to enhance the degradation of hemicellulose, which results in the formation of furans, while contributing to the initial dehydration/decomposition of cellulose [12,13,71]. All these reaction pathways are therefore likely to explain the decreases in the pyrolysis temperatures and activation energies highlighted in Section 3.
Secondly, AAEMs can play an important catalytic role as bridges linking adjacent oxygenated functional groups. As such, AAEMs are thus likely to promote recombination reactions favoring the formation of char. Hwang et al. notably reported in their study on the pyrolysis of MgCl2-impregnated yellow poplar that magnesium could enhance the repolymerization of the volatile molecules released, which would increase the average molecular weight, the viscosity and the solid content of the produced bio-oil despite its higher water content [13]. Such a phenomenon, here again, is well in line with the observations made based on the TG curves related to the AAEM-impregnated samples herein, which exhibited higher residual masses at 700 °C than in the case of raw wood regardless of the considered heating rate.

5. Conclusions

The pyrolysis of raw beech wood and samples impregnated with three AAEM catalysts (NaCl, KCl and MgCl2) was investigated at four heating rates by means of thermogravimetric analyses. Obtained results showed that the tested alkali and alkaline earth metals promote biomass decomposition by reducing the initial and peak pyrolysis temperatures. The alkaline earth metal (magnesium) was, moreover, found to exhibit a higher catalytic effect than alkali metals (sodium and potassium), as exemplified by its higher ability to reduce initial decomposition temperatures. A comprehensive kinetic analysis of measured data was then realized using two model-free methods (OFW and KAS models). In terms of highlights, the OFW and KAS models led to inferring very similar activation energies, which were, moreover, found to be lower for the samples impregnated with NaCl, KCl and MgCl2, thus corroborating the existence of a catalytic effect when such AAEM compounds are added to beech wood. Coupling the above-mentioned isoconversional models with model fitting and master plot approaches allowed to identify the D3 mechanism as being well suited to reproduce experimentally monitored data. While they allowed to simulate the evolution of the conversion degree of each sample as a function of the temperature regardless of the heating rate, the kinetic parameters inferred within this work also led to an illustration of increases of the pyrolysis rate constants when impregnating biomass with AAEMs.
To conclude, the main mechanisms at play during the catalytic pyrolysis of biomass (i.e., the enhanced cleavage of chemical bonds in biomass and the promotion of char formation reactions) have been discussed herein in a bid to interpret the observations made throughout this study. The trends depicted herein thus allowed to corroborate some speculated pathways proposed to account for the impact of AAEMs on the thermal degradation of woody biomass. Furthermore, obtained results also confirmed the interest of coupling varied modeling approaches to contribute to elucidating the fundamentals of AAEM-catalyzed pyrolysis, noting that complementary simulation tools including refined global kinetic schemes as well as phenomenological models would benefit from being considered in complementary studies, thus paving the way for future works to be undertaken.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27227662/s1, Table S1: Kinetic parameters derived from the use of the OFW and KAS approaches for different reaction models—Wood; Table S2: Kinetic parameters derived from the use of the OFW and KAS approaches for different reactions model—Wood + NaCl; Table S3: Kinetic parameters derived from the use of the OFW and KAS approaches for different reaction models—Wood + KCl; Table S4: Kinetic parameters derived from the use of the OFW and KAS approaches for different reaction models—Wood + MgCl2; Figure S1: Evolution of α as a function of the temperature for a heating rate of 10 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the OFW model integrating 10 different reaction mechanisms; Figure S2: Evolution of α as a function of the temperature for a heating rate of 5 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the OFW, KAS models integrating the D3 reaction mechanism; Figure S3: Evolution of α as a function of the temperature for a heating rate of 15 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the OFW, KAS models integrating the D3 reaction mechanism; Figure S4: Evolution of α as a function of the temperature for a heating rate of 30 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the OFW, KAS models integrating the D3 reaction mechanism.

Author Contributions

Conceptualization, R.L. and A.B.; methodology, W.W., R.L. and A.B.; software, W.W. and R.L.; validation, R.L.; formal analysis, W.W. and R.L.; investigation, W.W. and R.L.; resources, R.L. and A.B.; writing—original draft preparation, W.W. and R.L.; writing—review and editing, R.L.; visualization, W.W. and R.L.; supervision, R.L., A.B. and D.L.; project administration, R.L., A.B. and D.L.; funding acquisition, R.L. and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research has received the support of the French Ministry of Higher Education, Research and Innovation (Ministère de l’Enseignement supérieur, de la Recherche et de l’Innovation).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Collard, F.-X.; Blin, J. A review on pyrolysis of biomass constituents: Mechanisms and composition of the products obtained from the conversion of cellulose, hemicelluloses and lignin. Renew. Sustain. Energy Rev. 2014, 38, 594–608. [Google Scholar] [CrossRef]
  2. Broumand, M.; Albert-Green, S.; Yun, S.; Hong, Z.; Thomsson, M.J. Spray combustion of fast pyrolysis bio-oils: Applications, challenges, and potential solutions. Prog. Energy Combust. Sci. 2020, 79, 100834. [Google Scholar] [CrossRef]
  3. Chen, X.; Che, Q.; Li, S.; Liu, Z.; Yang, H.; Chen, Y.; Wang, X.; Shao, J.; Chen, H. Recent developments in lignocellulosic biomass catalytic fast pyrolysis: Strategies for the optimization of bio-oil quality and yield. Fuel Process. Technol. 2019, 196, 106180. [Google Scholar] [CrossRef]
  4. Qiu, B.; Tao, X.; Wang, J.; Liu, Y.; Li, S.; Chu, H. Research progress in the preparation of high-quality liquid fuels and chemicals by catalytic pyrolysis of biomass: A review. Energy Convers. Manag. 2022, 261, 115647. [Google Scholar] [CrossRef]
  5. Wang, W.; Lemaire, R.; Bensakhria, A.; Luart, D. Review on the catalytic effects of alkali and alkaline earth metals (AAEMs) including sodium, potassium, calcium and magnesium on the pyrolysis of lignocellulosic biomass and on the co-pyrolysis of coal with biomass. J. Anal. Appl. Pyrolysis 2022, 163, 105479. [Google Scholar] [CrossRef]
  6. Guo, F.; Liu, Y.; Wang, Y.; Li, X.; Li, T.; Guo, C. Pyrolysis kinetics and behavior of potassium-impregnated pine wood in TGA and a fixed-bed reactor. Energy Convers. Manag. 2016, 130, 184–191. [Google Scholar] [CrossRef]
  7. Shen, Y.; Zhang, N.; Zhang, S. Catalytic pyrolysis of biomass with potassium compounds for co-production of high-quality biofuels and porous carbons. Energy 2020, 190, 116431. [Google Scholar] [CrossRef]
  8. Hwang, H.; Oh, S.; Cho, T.-S.; Choi, I.-G.; Choi, J.W. Fast pyrolysis of potassium impregnated poplar wood and characterization of its influence on the formation as well as properties of pyrolytic products. Bioresour. Technol. 2013, 150, 359–366. [Google Scholar] [CrossRef]
  9. Zhang, H.; Ma, Y.; Shao, S.; Xiao, R. The effects of potassium on distributions of bio-oils obtained from fast pyrolysis of agricultural and forest biomass in a fluidized bed. Appl. Energy 2017, 208, 867–877. [Google Scholar] [CrossRef]
  10. Peng, C.; Zhang, G.; Yue, J.; Xu, G. Pyrolysis of lignin for phenols with alkaline additive. Fuel Process. Technol. 2014, 124, 212–221. [Google Scholar] [CrossRef]
  11. Veses, A.; Aznar, M.; Martínez, I.; Martínez, J.D.; López, J.M.; Navarro, M.V.; Callén, M.S.; Murillo, R.; García, T. Catalytic pyrolysis of wood biomass in an auger reactor using calcium-based catalysts. Bioresour. Technol. 2014, 162, 250–258. [Google Scholar] [CrossRef]
  12. Santana, J.A.; Sousa, N.G.; Cardoso, C.R.; Carvalho, W.S.; Ataíde, C.H. Sodium, zinc and magnesium chlorides as additives for soybean hulls pyrolysis: Influence on the temperature range of reactions and product selectivity. J. Therm. Anal. Calorim. 2016, 125, 471–481. [Google Scholar] [CrossRef]
  13. Hwang, H.; Oh, S.; Choi, I.-G.; Choi, J.W. Catalytic effects of magnesium on the characteristics of fast pyrolysis products—Bio-oil, bio-char, and non-condensed pyrolytic gas fractions. J. Anal. Appl. Pyrolysis 2015, 113, 27–34. [Google Scholar] [CrossRef]
  14. Vyazovkin, S.; Wight, C.A. Model-free and model-fitting approaches to kinetic analysis of isothermal and nonisothermal data. Thermochim. Acta 1999, 340–341, 53–68. [Google Scholar] [CrossRef]
  15. Lemaire, R.; Menage, D.; Seers, P. Study of the high heating rate devolatilization of bituminous and subbituminous coals-Comparison of experimentally monitored devolatilization profiles with predictions issued from single rate, two-competing rate, distributed activation energy and chemical percolation devolatilization models. J. Anal. Appl. Pyrolysis 2017, 123, 255–268. [Google Scholar] [CrossRef]
  16. Wang, J.; Zhao, B.; Liu, S.; Zhu, D.; Huang, F.; Yang, H.; Guan, H.; Song, A.; Xu, D.; Sun, L.; et al. Catalytic pyrolysis of biomass with Ni/Fe-CaO-based catalysts for hydrogen-rich gas: DFT and experimental study. Energy Convers. Manag. 2022, 254, 115246. [Google Scholar] [CrossRef]
  17. Sharifzadeh, M.; Sadeqzadeh, M.; Guo, M.; Borhani, T.N.; Murthy Konda, N.V.S.N.; Garcia, M.C.; Wang, L.; Hallett, J.; Shah, N. The multi-scale challenges of biomass fast pyrolysis and bio-oil upgrading: Review of the state of art and future research directions. Prog. Energy Combust. Sci. 2019, 71, 1–80. [Google Scholar] [CrossRef]
  18. Zhou, L.; Jia, Y.; Nguyen, T.-H.; Adesina, A.A.; Liu, Z. Hydropyrolysis characteristics and kinetics of potassium-impregnated pine wood. Fuel Process. Technol. 2013, 116, 149–157. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, M.; Wang, J.; Zhang, M.; Chen, M.; Zhu, X.; Min, F.; Tan, Z. Catalytic effects of eight inorganic additives on pyrolysis of pine wood sawdust by microwave heating. J. Anal. Appl. Pyrolysis 2008, 82, 145–150. [Google Scholar] [CrossRef]
  20. Wang, W.-L.; Ren, X.-Y.; Li, L.-F.; Chang, J.-M.; Cai, L.-P.; Geng, J. Catalytic effect of metal chlorides on analytical pyrolysis of alkali lignin. Fuel Process. Technol. 2015, 134, 345–351. [Google Scholar] [CrossRef]
  21. Safar, M.; Lin, B.-J.; Chen, W.-H.; Langauer, D.; Chang, J.-S.; Raclavska, H.; Pétrissans, A.; Rousset, P.; Pétrissans, M. Catalytic effects of potassium on biomass pyrolysis, combustion and torrefaction. Appl. Energy 2019, 235, 346–355. [Google Scholar] [CrossRef]
  22. Balasundram, V.; Ibrahim, N.; Md Kasmani, R.; Hamid, M.K.A.; Isha, R.; Hasbullah, H.; Ali, R.R. Thermogravimetric catalytic pyrolysis and kinetic studies of coconut copra and rice husk for possible maximum production of pyrolysis oil. J. Clean. Prod. 2017, 167, 218–228. [Google Scholar] [CrossRef]
  23. Wu, Y.; Zhu, J.; Wang, Y.; Yang, H.; Jin, L.; Hu, H. Insight into co-pyrolysis interactions of Pingshuo coal and high-density polyethylene via in-situ Py-TOF-MS and EPR. Fuel 2021, 303, 121199. [Google Scholar] [CrossRef]
  24. Senum, G.I.; Yang, R.T. Rational approximations of the integral of the Arrhenius function. J. Therm. Anal. 1977, 11, 445–447. [Google Scholar] [CrossRef]
  25. Flynn, J.H. The isoconversional method for determination of energy of activation at constant heating rates: Corrections for the Doyle approximation. J. Therm. Anal. 1983, 27, 95–102. [Google Scholar] [CrossRef]
  26. Doyle, C.D. Kinetic analysis of thermogravimetric data. J. Appl. Polym. Sci. 1961, 5, 285–292. [Google Scholar] [CrossRef]
  27. Ozawa, T. A new method of analyzing thermogravimetric data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef] [Green Version]
  28. Flynn, J.H.; Wall, L.A. General treatment of the thermogravimetry of polymers. J. Res. Natl. Bur. Stand. Sect. A 1966, 70, 487–523. [Google Scholar] [CrossRef]
  29. Wang, S.; Dai, G.; Yang, H.; Luo, Z. Lignocellulosic biomass pyrolysis mechanism: A state-of-the-art review. Prog. Energy Combust. Sci. 2017, 62, 33–86. [Google Scholar] [CrossRef]
  30. Khawam, A.; Flanagan, D.R. Complementary use of model-free and modelistic methods in the analysis of solid-state kinetics. J. Phys. Chem. B 2005, 109, 10073–10080. [Google Scholar] [CrossRef]
  31. Moriana, R.; Zhang, Y.; Mischnick, P.; Li, J.; Ek, M. Thermal degradation behavior and kinetic analysis of spruce glucomannan and its methylated derivatives. Carbohydr. Polym. 2014, 106, 60–70. [Google Scholar] [CrossRef]
  32. Wang, S.; Lin, H.; Ru, B.; Dai, G.; Wang, X.; Xiao, G.; Luo, Z. Kinetic modeling of biomass components pyrolysis using a sequential and coupling method. Fuel 2016, 185, 763–771. [Google Scholar] [CrossRef]
  33. Prasad Nath, H.; Kumar Dutta, B.; Kalita, D.; Saikia, B.K.; Saikia, N. Evaluation of the effect of high sulfur subbituminous coal on the devolatilization of biomass residue by using model free, model fitting and combined kinetic methods. Fuel 2022, 310, 122235. [Google Scholar] [CrossRef]
  34. Chen, C.; Ma, X.; He, Y. Co-pyrolysis characteristics of microalgae Chlorella vulgaris and coal through TGA. Bioresour. Technol. 2012, 117, 264–273. [Google Scholar] [CrossRef]
  35. He, Y.; Chang, C.; Li, P.; Han, X.; Li, H.; Fang, S.; Chen, J.; Ma, X. Thermal decomposition and kinetics of coal and fermented cornstalk using thermogravimetric analysis. Bioresour. Technol. 2018, 259, 294–303. [Google Scholar] [CrossRef]
  36. Florentino-Madiedo, L.; Vega, M.F.; Diaz-Faes, E.; Barriocanal, C. Evaluation of synergy during co-pyrolysis of torrefied sawdust, coal and paraffin. A kinetic and thermodynamic dataset. Data Brief 2021, 37, 107170. [Google Scholar] [CrossRef]
  37. Škvára, F.; Šesták, J. Computer calculation of the mechanism and associated kinetic data using a non-isothermal integral method. J. Therm. Anal. 1975, 8, 477–489. [Google Scholar] [CrossRef]
  38. Chen, F.; Fu, L.; Feng, L.; Liu, C.; Ren, B. Non-isothermal decomposition kinetics of diosgenin. Russ. J. Phys. Chem. 2013, 87, 1611–1614. [Google Scholar] [CrossRef]
  39. Zhao, N.; Li, J.; Gong, H.; An, T.; Zhao, F.; Yang, A.; Hu, R.; Ma, H. Effects of α-Fe2O3 nanoparticles on the thermal behavior and non-isothermal decomposition kinetics of nitrocellulose. J. Anal. Appl. Pyrolysis 2016, 120, 165–173. [Google Scholar] [CrossRef]
  40. Zhao, B.; Xu, X.; Li, H.; Chen, X.; Zeng, F. Kinetics evaluation and thermal decomposition characteristics of co-pyrolysis of municipal sewage sludge and hazelnut shell. Bioresour. Technol. 2018, 247, 21–29. [Google Scholar] [CrossRef]
  41. Ma, Z.; Wang, J.; Yang, Y.; Zhang, Y.; Zhao, C.; Yu, Y.; Wang, S. Comparison of the thermal degradation behaviors and kinetics of palm oil waste under nitrogen and air atmosphere in TGA-FTIR with a complementary use of model-free and model-fitting approaches. J. Anal. Appl. Pyrolysis 2018, 134, 12–24. [Google Scholar] [CrossRef]
  42. Lv, X.; Huang, R.; Wu, Q.; Xu, B.; Zhang, J. Non-isothermal reduction kinetics during vacuum carbothermal reduction of ilmenite concentrate. Vacuum 2019, 160, 139–145. [Google Scholar] [CrossRef]
  43. Liu, Y.; Yang, L.; Zhang, Y. Thermal behavior and kinetic decomposition of sweet potato starch by non-isothermal procedures. Arch. Thermodyn. 2019, 40, 67–82. [Google Scholar] [CrossRef]
  44. Vyazovkin, S.; Burnham, A.K.; Criado, J.M.; Pérez-Maqueda, L.A.; Popescu, C.; Sbirrazzuoli, N. ICTAC Kinetics Committee recommendations for performing kinetic computations on thermal analysis data. Thermochim. Acta 2011, 520, 1–19. [Google Scholar] [CrossRef]
  45. Vyazovkin, S.; Chrissafis, K.; Di Lorenzo, M.L.; Koga, N.; Pijolat, M.; Roduit, B.; Sbirrazzuoli, N.; Sunol, J.J. ICTAC Kinetics Committee recommendations for collecting experimental thermal analysis data for kinetic computations. Thermochim. Acta 2014, 590, 1–23. [Google Scholar] [CrossRef]
  46. Sánchez-Jiménez, P.E.; Pérez-Maqueda, L.A.; Perejón, A.; Criado, J.M. Generalized master plots as a straightforward approach for determining the kinetic model: The case of cellulose pyrolysis. Thermochim. Acta 2013, 552, 54–59. [Google Scholar] [CrossRef] [Green Version]
  47. Helsen, L.; Van den Bulck, E. Kinetics of the low-temperature pyrolysis of chromated copper arsenate-treated wood. J. Anal. Appl. Pyrolysis 2000, 53, 51–79. [Google Scholar] [CrossRef]
  48. Wang, J.; Zhang, M.; Chen, M.; Min, F.; Zhang, S.; Ren, Z.; Yan, Y. Catalytic effects of six inorganic compounds on pyrolysis of three kinds of biomass. Thermochim. Acta 2006, 444, 110–114. [Google Scholar] [CrossRef]
  49. Nowakowski, D.J.; Jones, J.M.; Brydson, R.M.D.; Ross, A.B. Potassium catalysis in the pyrolysis behaviour of short rotation willow coppice. Fuel 2007, 86, 2389–2402. [Google Scholar] [CrossRef]
  50. Shimada, N.; Kawamoto, H.; Saka, S. Different action of alkali/alkaline earth metal chlorides on cellulose pyrolysis. J. Anal. Appl. Pyrolysis 2008, 81, 80–87. [Google Scholar] [CrossRef]
  51. Yu, Y.; Liu, D.; Wu, H. Formation and characteristics of reaction intermediates from the fast pyrolysis of NaCl- and MgCl2- loaded celluloses. Energy Fuels 2014, 28, 245–253. [Google Scholar] [CrossRef]
  52. Gao, P.; Xue, L.; Lu, Q.; Dong, C. Effects of alkali and alkaline earth metals on N-containing species release during rice straw pyrolysis. Energies 2015, 8, 13021–13032. [Google Scholar] [CrossRef] [Green Version]
  53. Leng, E.; Wang, Y.; Gong, X.; Zhang, B.; Zhang, Y.; Xu, M. Effect of KCl and CaCl2 loading on the formation of reaction intermediates during cellulose fast pyrolysis. Proc. Combust. Inst. 2017, 36, 2263–2270. [Google Scholar] [CrossRef]
  54. Jensen, A.; Dam-Johansen, K.; Wójtowicz, M.A.; Serio, M.A. TG-FTIR study of the influence of potassium chloride on wheat straw pyrolysis. Energy Fuels 1998, 12, 929–938. [Google Scholar] [CrossRef]
  55. Shah, M.H.; Deng, L.; Bennadji, H.; Fisher, E.M. Pyrolysis of potassium-doped wood at the centimeter and submillimeter scales. Energy Fuels 2015, 29, 7350–7357. [Google Scholar] [CrossRef] [Green Version]
  56. Zhu, C.; Maduskar, S.; Paulsen, A.D.; Dauenhauer, P.J. Alkaline-earth-metal-catalyzed thin-film pyrolysis of cellulose. ChemCatChem 2016, 8, 818–829. [Google Scholar] [CrossRef]
  57. Dalluge, D.L.; Kim, K.H.; Brown, R.C. The influence of alkali and alkaline earth metals on char and volatile aromatics from fast pyrolysis of lignin. J. Anal. Appl. Pyrolysis 2017, 127, 385–393. [Google Scholar] [CrossRef]
  58. Cai, J.; Xu, D.; Dong, Z.; Yu, X.; Yang, Y.; Banks, S.W.; Bridgwater, A.V. Processing thermogravimetric analysis data for isoconversional kinetic analysis of lignocellulosic biomass pyrolysis: Case study of corn stalk. Renew. Sustain. Energy Rev. 2018, 82, 2705–2715. [Google Scholar] [CrossRef] [Green Version]
  59. Tang, W.; Chen, D.; Sheng, Y. Study on non-isothermal decomposition kinetics of ephedrini hydrochloridum. Wuhan Univ. J. Nat. Sci. 2003, 8, 443–446. [Google Scholar] [CrossRef]
  60. Aboulkas, A.; El harfi, K.; El Bouadili, A. Thermal degradation behaviors of polyethylene and polypropylene. Part I: Pyrolysis kinetics and mechanisms. Energy Convers. Manag. 2010, 51, 1363–1369. [Google Scholar] [CrossRef]
  61. Huang, M.-X.; Zhou, C.-R.; Han, X.-W. Investigation of thermal decomposition kinetics of taurine. J. Therm. Anal. Calorim. 2013, 113, 589–593. [Google Scholar] [CrossRef]
  62. Anca-Couce, A.; Berger, A.; Zobel, N. How to determine consistent biomass pyrolysis kinetics in a parallel reaction scheme. Fuel 2014, 123, 230–240. [Google Scholar] [CrossRef]
  63. Mishra, G.; Bhaskar, T. Non isothermal model free kinetics for pyrolysis of rice straw. Bioresour. Technol. 2014, 169, 614–621. [Google Scholar] [CrossRef] [PubMed]
  64. Hu, M.; Chen, Z.; Wang, S.; Guo, D.; Ma, C.; Zhou, Y.; Chen, J.; Laghari, M.; Fazal, S.; Xiao, B.; et al. Thermogravimetric kinetics of lignocellulosic biomass slow pyrolysis using distributed activation energy model, Fraser–Suzuki deconvolution, and iso-conversional method. Energy Convers. Manag. 2016, 118, 1–11. [Google Scholar] [CrossRef]
  65. Rony, A.H.; Kong, L.; Lu, W.; Dejam, M.; Adidharma, H.; Gasem, K.A.M.; Zheng, Y.; Norton, U.; Fan, M. Kinetics, thermodynamics, and physical characterization of corn stover (Zea mays) for solar biomass pyrolysis potential analysis. Bioresour. Technol. 2019, 284, 466–473. [Google Scholar] [CrossRef]
  66. Masawat, N.; Atong, D.; Sricharoenchaikul, V. Thermo-kinetics and product analysis of the catalytic pyrolysis of Pongamia residual cake. J. Environ. Manag. 2019, 242, 238–245. [Google Scholar] [CrossRef]
  67. Hu, L.; Wei, X.-Y.; Guo, X.-H.; Lv, H.-P.; Wang, G.-H. Investigation on the kinetic behavior, thermodynamic and volatile products analysis of chili straw waste pyrolysis. J. Environ. Chem. Eng. 2021, 9, 105859. [Google Scholar] [CrossRef]
  68. Masnadi, M.S.; Habibi, R.; Kopyscinski, J.; Hill, J.M.; Bi, X.; Lim, C.J.; Ellis, N.; Grace, J.R. Fuel characterization and co-pyrolysis kinetics of biomass and fossil fuels. Fuel 2014, 117, 1204–1214. [Google Scholar] [CrossRef]
  69. Ding, Y.; Ezekoye, O.A.; Lu, S.; Wang, C.; Zhou, R. Comparative pyrolysis behaviors and reaction mechanisms of hardwood and softwood. Energy Convers. Manag. 2017, 132, 102–109. [Google Scholar] [CrossRef]
  70. Zhao, N.; Li, B.-X. The effect of sodium chloride on the pyrolysis of rice husk. Appl. Energy 2016, 178, 346–352. [Google Scholar] [CrossRef]
  71. Xu, Q.; Ma, X.; Yu, Z.; Cai, Z. A kinetic study on the effects of alkaline earth and alkali metal compounds for catalytic pyrolysis of microalgae using thermogravimetry. Appl. Therm. Eng. 2014, 73, 357–361. [Google Scholar] [CrossRef]
Figure 1. Diagram of sample preparation and analysis procedures.
Figure 1. Diagram of sample preparation and analysis procedures.
Molecules 27 07662 g001
Figure 2. Evolution of mass loss noted ‘TG’ (a,c,e,g) and derivate mass loss rate noted ‘dTG’ (b,d,f,h) as a function of the temperature for heating rates of 5 °C/min (a,b), 10 °C/min (c,d), 15 °C/min (e,f) and 30 °C/min (g,h).
Figure 2. Evolution of mass loss noted ‘TG’ (a,c,e,g) and derivate mass loss rate noted ‘dTG’ (b,d,f,h) as a function of the temperature for heating rates of 5 °C/min (a,b), 10 °C/min (c,d), 15 °C/min (e,f) and 30 °C/min (g,h).
Molecules 27 07662 g002
Figure 3. Linear fitting obtained with the OFW model (i.e., plots of ln β as a function of 1000 R × T for conversion degrees between 10 and 90%) for (a) wood, (b) wood + NaCl, (c) wood + KCl and (d) wood + MgCl2.
Figure 3. Linear fitting obtained with the OFW model (i.e., plots of ln β as a function of 1000 R × T for conversion degrees between 10 and 90%) for (a) wood, (b) wood + NaCl, (c) wood + KCl and (d) wood + MgCl2.
Molecules 27 07662 g003
Figure 4. Linear fitting obtained with the KAS model (i.e., plots of ln β T 2 as a function of 1000 R × T for conversion degrees between 10 and 90%) for (a) wood, (b) wood + NaCl, (c) wood + KCl and (d) wood + MgCl2.
Figure 4. Linear fitting obtained with the KAS model (i.e., plots of ln β T 2 as a function of 1000 R × T for conversion degrees between 10 and 90%) for (a) wood, (b) wood + NaCl, (c) wood + KCl and (d) wood + MgCl2.
Molecules 27 07662 g004
Figure 5. Example of linear data fitting issued from the implementation of the Šatava–Šesták model while considering the D3 reaction model for a heating rate of 10 °C/min.
Figure 5. Example of linear data fitting issued from the implementation of the Šatava–Šesták model while considering the D3 reaction model for a heating rate of 10 °C/min.
Molecules 27 07662 g005
Figure 6. Isoconversional model plot of E a obtained with the OFW, KAS and Šatava–Šesták models in the case of samples impregnated with (a) NaCl and (b) KCl, as examples.
Figure 6. Isoconversional model plot of E a obtained with the OFW, KAS and Šatava–Šesták models in the case of samples impregnated with (a) NaCl and (b) KCl, as examples.
Molecules 27 07662 g006
Figure 7. Evolution of α as a function of the temperature for a heating rate of 10 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the KAS model integrating 10 different reaction mechanisms.
Figure 7. Evolution of α as a function of the temperature for a heating rate of 10 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the KAS model integrating 10 different reaction mechanisms.
Molecules 27 07662 g007
Figure 8. Comparison of experimental ( β = 10 °C/min) and theoretical master plots for samples impregnated with (a) NaCl and (b) KCl.
Figure 8. Comparison of experimental ( β = 10 °C/min) and theoretical master plots for samples impregnated with (a) NaCl and (b) KCl.
Molecules 27 07662 g008
Figure 9. Evolution of α as a function of the temperature for a heating rate of 10 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the OFW and KAS models integrating the D3 reaction mechanism.
Figure 9. Evolution of α as a function of the temperature for a heating rate of 10 °C/min in the case of (a) wood and samples impregnated with (b) NaCl, (c) KCl and (d) MgCl2; comparison of experimental data (noted ‘Exp’) with predicted ones obtained from the use of the OFW and KAS models integrating the D3 reaction mechanism.
Molecules 27 07662 g009
Table 1. Proximate and ultimate analyses of the tested beech wood sample.
Table 1. Proximate and ultimate analyses of the tested beech wood sample.
SampleProximate AnalysisUltimate Analysis
Fixed Carbon (wt%, db )Volatiles (wt%, db) Ash (wt%, db)C (wt%, daf )H (wt%, daf)O * (wt%, daf)N (wt%, daf)S (wt%, daf)
Beech wood14.3684.481.1651.65.443.0--
db: dry basis; daf: dry ash-free basis; * calculated by difference.
Table 2. Characteristic decomposition temperatures, maximal mass loss rates and residual masses at 700 °C for raw and impregnated beech wood samples for β values of 5, 10, 15 and 30 °C/min.
Table 2. Characteristic decomposition temperatures, maximal mass loss rates and residual masses at 700 °C for raw and impregnated beech wood samples for β values of 5, 10, 15 and 30 °C/min.
β 5 °C/min10 °C/min
Sample Ti (°C)Tf (°C)Tp (°C)dTGmax (%/min)TG700°C (wt%)Ti (°C)Tf (°C)Tp (°C)dTGmax (%/min)TG700°C (wt%)
Wood257.5 (±0.9)411.9 (±37.1)340 (±2)−3.74 (±0.23)30.9 (±2.1)267.5 (±0.2)414.3 (±6.5)350 (±1)−6.55 (±0.29)30.4 (±2.8)
Wood + NaCl254.7 (±0.6)408.8 (±11.0)315 (±1)−3.49 (±0.04)40.0 (±2.0)264.7 (±1.1)425.5 (±18.0)325 (±2)−5.88 (±0.31)41.2 (±2.3)
Wood + KCl246.7 (±0.8)411.1 (±5.0)309 (±1)−3.51 (±0.21)36.9 (±2.8)256.8 (±0.5)421.4 (±5.1)320 (±3)−5.92 (±0.19)39.3 (±1.7)
Wood + MgCl2216.0 (±5.1)451.6 (±3.4)329(±0)−2.59 (±0.08)34.9 (±1.0)225.4 (±1.6)458.2 (±3.1)340 (±1)−4.72 (±0.11)34.5 (±1.4)
β 15 °C/min30 °C/min
SampleTi(°C)Tf(°C)Tp(°C)dTGmax (%/min)TG700°C (wt%)Ti(°C)Tf(°C)Tp(°C)dTGmax (%/min)TG700°C (wt%)
Wood275.5 (±2.4)419.7 (±10.3)353 (±1)−8.32 (±0.32)32.7 (±3.1)286.4 (±2.6)431.6 (±3.3)361 (±3)−13.9 (±0.4)27.6 (±2.1)
Wood + NaCl273.7 (±1.4)420.1 (±5.8)329 (±2)−8.61 (±0.39)36.6 (±2.9)287.6 (±1.5)435.0 (±5.1)344 (±1)−12.9 (±0.5)36.5 (±2.7)
Wood + KCl264.4 (±1.0)433.1 (±11.6)323 (±1)−8.15 (±0.45)37.1 (±3.3)278.0 (±3.2)434.2 (±4.9)337 (±3)−12.5 (±0.3)37.5 (±1.6)
Wood + MgCl2231.4 (±1.1)470.2 (±21.7)342 (±3)−6.18 (±0.44)33.4 (±2.6)242.6 (±13.4)484.4 (±55.6)331 (±5)−9.7 (±0.5)33.6 (±3.3)
Table 3. Activation energies issued from the implementation of the OFW and KAS models.
Table 3. Activation energies issued from the implementation of the OFW and KAS models.
α WoodWood + NaClWood + KClWood + MgCl2
E a   ( kJ / mol ) R2 E a   ( kJ / mol ) R2 E a   ( kJ / mol ) R2 E a   ( kJ / mol ) R2
OFW
10%143.70.9978125.40.9938128.60.9961134.40.9991
20%155.80.9982134.00.9926135.10.9967161.00.9983
30%160.50.9963140.00.9953138.20.9974147.60.9946
40%157.90.9954139.00.9984137.10.9980146.90.9937
50%156.50.9966137.30.9993136.30.9983152.30.9986
60%154.40.9971136.80.9994136.30.9983153.20.9978
70%153.00.9971137.20.9993137.30.9982152.70.9955
80%155.00.9967155.00.9949162.10.9972171.20.9744
90%309.50.9096239.00.8348248.50.8853224.70.9647
KAS
10%142.10.9976122.90.9927126.30.9954133.00.9990
20%154.40.9979131.60.9914132.80.9962160.50.9981
30%159.00.9958137.60.9945135.80.9969145.90.9939
40%156.00.9947136.30.9981134.30.9977144.70.9929
50%154.30.9961134.30.9991133.40.9980150.10.9984
60%151.80.9966133.70.9993133.30.9979150.80.9974
70%150.30.9966134.10.9992134.30.9979150.10.9947
80%152.20.9962152.60.9941160.10.9967169.20.9709
90%314.20.9034239.90.8211249.90.8757223.90.9607
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, W.; Lemaire, R.; Bensakhria, A.; Luart, D. Thermogravimetric Analysis and Kinetic Modeling of the AAEM-Catalyzed Pyrolysis of Woody Biomass. Molecules 2022, 27, 7662. https://doi.org/10.3390/molecules27227662

AMA Style

Wang W, Lemaire R, Bensakhria A, Luart D. Thermogravimetric Analysis and Kinetic Modeling of the AAEM-Catalyzed Pyrolysis of Woody Biomass. Molecules. 2022; 27(22):7662. https://doi.org/10.3390/molecules27227662

Chicago/Turabian Style

Wang, Wei, Romain Lemaire, Ammar Bensakhria, and Denis Luart. 2022. "Thermogravimetric Analysis and Kinetic Modeling of the AAEM-Catalyzed Pyrolysis of Woody Biomass" Molecules 27, no. 22: 7662. https://doi.org/10.3390/molecules27227662

APA Style

Wang, W., Lemaire, R., Bensakhria, A., & Luart, D. (2022). Thermogravimetric Analysis and Kinetic Modeling of the AAEM-Catalyzed Pyrolysis of Woody Biomass. Molecules, 27(22), 7662. https://doi.org/10.3390/molecules27227662

Article Metrics

Back to TopTop