-
Nitrogen atmosphere was selected to conduct the pyrolysis tests of beech wood using a NETZSCH STA449F3 thermal analyzer. Three different heating rates of 5, 10 and 20 K/min were used to heat the 5−7 mg wood powders from 290 to 1,070 K. Heating rate plays an important role in collecting proper experimental data. Very slow heating rate, such as 1 K/min, allows the reactions to come closer to equilibrium and there is less thermal lag in the sample and apparatus. Contrarily, high heating rates give faster experiments, which are more representative of the heating rates in fires, but deviate more from equilibrium and result in greater thermal lag. Larger heating rates are suitable for finding a wide range of decomposition, while smaller heating rates show better performance in the separation of individual events. Heating rates used in tests should preferably be in the range of 1−20 K/min as recommended by ICTAC (International Confederation for Thermal Analysis and Calorimetry) Kinetics Committee[18]. Consequently, these three representative heating rates, 5, 10 and 20 K/min, were used in this study. A ceramic crucible was employed and a small hole was set in the center of the crucible lid to allow the release of volatiles during pyrolysis. All samples were dried in an oven for at least 72 h before tests to minimize the impact of moisture. Given that thermogravimetric experiments are highly reproducible, only three replicate experiments were performed at each heating rate to estimate the experimental uncertainty.
Model-free methods
-
Utilization of a suitable kinetic parameter calculation method allows a rough estimate of the search range for the subsequent optimization. Model-free methods, including the Kissinger method and isoconversional methods, are introduced in this section and will be utilized to calculate the kinetics of wood. Previous studies showed that the kinetic parameters calculated by the Kissinger method are very close to other isoconversional methods, such as KAS, Tang, DAEM methods.
In a single-step pyrolysis reaction, the degree of conversion of the solid, α, can be expressed as:
$ \alpha = \displaystyle\frac{{\left( {{m_0} - m} \right)}}{{\left( {{m_0} - {m_\infty }} \right)}} $ (1) where
,$m$ and${m_0}$ are the transient, initial and final masses of the sample, respectively. The reaction kinetic equation is:${m_\infty }$ $ \displaystyle\frac{{d\alpha }}{{dt}} = \lambda \left( T \right)f\left( \alpha \right) $ (2) where
is time, T is the absolute temperature,$t$ is the differential form of the reaction model and$f\left( \alpha \right)$ is the reaction rate constant. For a first order reaction,$\lambda $ can be expressed as:$\lambda $ $ \lambda = Aexp\left( { - \displaystyle\frac{{{E_a}}}{{RT}}} \right) $ (3) where
,$A$ and${E_a}$ refer to the pre-exponential factor, activation energy and ideal gas constant, respectively. With constant heating rate$R$ , Eq. (2) can be converted to:$\beta = dT/dt$ $ \displaystyle\frac{{d\alpha }}{{dT}} = \displaystyle\frac{A}{\beta }exp\left( { - \displaystyle\frac{{{E_a}}}{{RT}}} \right)f\left( \alpha \right) $ (4) The integral function of the conversion rate,
, can be expressed as:$g\left( \alpha \right)$ $ g\left( \alpha \right) = \mathop \int \nolimits_0^\alpha \displaystyle\frac{{d\alpha }}{{f\left( \alpha \right)}} = \displaystyle\frac{A}{\beta }\mathop \int \nolimits_{{T_0}}^T exp\left( { - \displaystyle\frac{{{E_a}}}{{RT}}} \right)dT $ (5) where
is the initial temperature.${T_0}$ Kissinger method
-
The Kissinger[19] method uses the peak temperature of the reaction peak for plotting, and the equation can be expressed as:
$ ln\left( {\displaystyle\frac{{{\beta _n}}}{{{T_{p,n}}^2}}} \right) = ln\left( { - \displaystyle\frac{{AR}}{{{E_a}}}f'\left( {{\alpha _p}} \right)} \right) - \displaystyle\frac{{{E_a}}}{{R{T_{p,n}}}} $ (6) where the subscript
and$n$ refer to the n-th heating rate and peak of MLR curve. When the model is a first-order kinetic model, the value of$p$ is −1, Eq. (6) can be further simplified as:$f'\left( {{\alpha _p}} \right)$ $ ln\left( {\displaystyle\frac{{{\beta _n}}}{{{T_{p,n}^2}}}} \right) = ln\left( {\displaystyle\frac{{AR}}{{{E_a}}}} \right) - \displaystyle\frac{{{E_a}}}{{R{T_{p,n}}}} $ (7) A straight line can be obtained by plotting
versus$ln\left( {{\beta _n}/{T_{p,n}^2}} \right)$ . The slope and intercept of the line can be employed to estimate$1/{T_{p,n}}$ and A, respectively.${E_a}$ KAS method
-
The formula for the KAS[20] method is expressed as:
$ ln\left( {\displaystyle\frac{\beta }{{{T^2}}}} \right) = ln \displaystyle\frac{{AR}}{{{E_a}g\left( \alpha \right)}} - \displaystyle\frac{{{E_a}}}{{RT}} $ (8) Using this formula, a straight line can be made by plotting
to$ln\left( {\beta /{T^2}} \right)$ at different heating rates and arbitrary conversion rate. The slope of this line can be used to derive$1/T$ . For first order reactions, the reaction mechanism functions are${E_a}$ and$f\left( \alpha \right) = 1 - \alpha $ , and Eq. (8) becomes:$g\left( \alpha \right) = - ln\left( {1 - \alpha } \right)$ $ ln\left( {\displaystyle\frac{\beta }{{{T^2}}}} \right) = ln\left( {\displaystyle\frac{{AR}}{{ - {E_a}ln\left( {1 - \alpha } \right)}}} \right) - \displaystyle\frac{{{E_a}}}{{RT}} $ (9) Tang method
-
Tang et al.[21] proposed the Arrhenius temperature integral approximation for Eq. (5):
$ - lnP\left( x \right) = 0.377739 + 1.894661ln\left( x \right) + 1.00145x $ (10) $ x = \displaystyle\frac{{{E_a}}}{{RT}} $ (11) $ P(x) = \displaystyle\frac{{g(\alpha )\beta R}}{{A{E_a}}} $ (12) Substituting the approximation into Eq. (5) and taking logarithms on both sides, the expression of the Tang method is obtained:
$ \begin{split} ln\left( {\frac{\beta }{{{T^{1.894661}}}}} \right) =& ln\left( {\frac{{A{E_a}}}{{Rg\left( \alpha \right)}}} \right) + 3.635041 -\\& 1.894661ln{E_a} - 1.001450\frac{{{E_a}}}{{RT}}\end{split} $ (13) At different heating rates and arbitrary conversion rate, plotting
versus$ln\left( {\beta /{T^{1.894661}}} \right)$ makes a straight, and$1/T$ can be estimated based on the slope. Similarly, for a first order reaction, Eq. (13) can be transformed into:${E_a}$ $\begin{split} ln\left( {\frac{\beta }{{{T^{1.894661}}}}} \right) =\;& ln\left( {\frac{{A{E_a}}}{{ - Rln\left( {1 - \alpha } \right)}}} \right) + 3.635041- \\& 1.894661ln{E_a} - 1.001450\frac{{{E_a}}}{{RT}}\end{split} $ (14) The pre-exponential factor,
, can be calculated from the intercept.$A$ DAEM method
-
DAEM[22] is an effective method to study the reaction behavior of complex systems including an infinite number of parallel first order reactions. The simplified DAEM model can be expressed as:
$ ln\left( {\displaystyle\frac{\beta }{{{T^2}}}} \right) = - \displaystyle\frac{E}{{RT}} + ln\left( {\displaystyle\frac{{AR}}{E}} \right) + 0.6075 $ (15) Similarly,
and${E_a}$ can be calculated based on the slope and intercept of the linearly fitted line of$A$ to$ln ( {\beta /{T^2}} )$ .$1/T$ -
To simulate the measured mass and mass loss rate (MLR) of sample in the thermogravimetric analysis experiments, a 0D numerical model for the pyrolysis of a thermally thin solid is developed. The 0D model implies the temperature gradient inside the condensed phase is neglected compared to the traditional 1D heat transfer model. The general forms of pyrolysis reaction and reaction rate are:
$ {\theta _1}Com{p_1} + {\theta _2}Com{p_2} \to {\theta _3}Com{p_3} + {\theta _4}Com{p_4} $ (16) $ {w}_{j}=-{c}_{Com{p}_{1}}^{{n}_{j,1}}{c}_{Com{p}_{2}}^{{n}_{j,2}}{A}_{j}exp\left(-\displaystyle\frac{{E}_{a,j}}{RT}\right)\text{ }\left(j=1,2,3\cdots \right) $ (17) where
denotes component,$Comp$ is the stoichiometric coefficient by mass,$\theta $ and$w$ are reaction rate and the reactant concentration index,$n$ is the transient mass of the reactant after normalization, the subscript$c$ refers to the j-th reaction. The transient mass change rate of each component is calculated as:$j$ $ \frac{{d{c_i}}}{{dt}} = \mathop \sum \limits_{j = 1}^{{N_j}} {\theta _{j,i}}{w_j} $ (18) where the subscript i denotes i-th component,
is the total number of reactions.${N_j}$ is positive or negative when the i-th component serves as a reactant or a product. The total transient mass (m) and MLR can be calculated as:${\theta _{j,i}}$ $ m = \mathop \sum \limits_{s = 1}^{{N_s}} {c_s};\quad MLR = \mathop \sum \limits_{g = 1}^{{N_g}} \mathop \sum \limits_{j = 1}^{{N_j}} {V_{j,g}}{w_j} $ (19) where
and${N_s}$ are the total numbers of solid and gaseous components, respectively. To commence simulation, the initial mass fraction of the starting reactants, the mass stoichiometry coefficients, and the three components of Arrhenius kinetics for each reaction need to be assigned. The initial values of m and MLR for the intermediate and final products are set to be zero.${N_g}$ Optimization algorithms
-
Three most frequently utilized optimization algorithms, GA, PSO, and SCE, are used to compare their performance by determining kinetics of wood combining the numerical model and experimental results. Detailed information of the three algorithms was introduced in our recent publications[23,24]. A four-component reaction scheme is applied to describe the pyrolysis of wood, namely the evaporation of water, pyrolysis of cellulose, hemicellulose and lignin. Each reaction includes four unknown parameters, namely
,${E_a}$ , stoichiometric coefficient of solid product, and the reaction order. To gain comparison purpose, the initial search ranges, population sizes, iteration numbers, and objective functions of the three algorithms when implementing optimizations are set to be identical.$A$ GA (Genetic algorithm)
-
GA simulates the process of population evolution, adopting a series of genetic operations such as selection, crossover and mutation for the current population to create a new generation and gradually progress the population to a state close to the optimal solution. In GA, each set of unknown parameters is referred to as an individual, and a combination of tens to thousands of individuals is defined as a population. Offspring, namely all potential solutions, are continuously produced by the overall population through genetic, crossover and mutation operations. An objective function, also known as fitness function, is essential for the assessment process. The objective function utilized in current study takes both m and MLR into account, as recommended by ICTAC Kinetics Committee[25]:
$ {R^2} = \displaystyle\mathop \sum \limits_{l = 1}^{{N_m}} \displaystyle\frac{{\left( {{m_{l,exp}} - {m_{l,num}}} \right)}}{{{m_{l,exp}} - {{\bar m}_{l,exp}}}} + \mathop \sum \limits_{l = 1}^{{N_{MLR}}} \displaystyle\frac{{\left( {ML{R_{l,exp}} - ML{R_{l,num}}} \right)}}{{ML{R_{l,exp}} - {{\overline {MLR} }_{l,exp}}}} $ (20) where
,${N_m}$ are the total numbers of experimental data points of m and MLR, the subscript exp and num denote experimental and numerical values, respectively,${N_{MLR}}$ and$\bar m$ are the average values of experimental m and MLR, respectively. A lower value of R2 generally represents higher accuracy of the algorithm and better fit between experimental and simulation results. This objective function will also be used in the following PSO and SCE algorithms.$\overline {MLR} $ PSO (particle swarm optimization) algorithm
-
Inspired by the behaviors of bird populations, PSO is developed as an alternative optimization algorithm. PSO algorithm includes velocity and position models, where the velocities of particles are used to update the positions of the particles, and the positions represent the potential solutions in the search ranges. The velocity and position updating processes require each particle to keep in mind the previous optimal position as well as the global optimal position searched by all particles. At the beginning of iteration, the velocities and positions of the particles are randomly assigned according to a specified range. Then they are updated according to the following relationship:
$ v_{ij}^{k + 1} = v_{ij}^k + {q_1}(x_{ij}^{pb} - x_{ij}^k) + {q_2}(x_{ij}^{gb} - x_{ij}^k) $ (21) $ x_{ij}^{k + 1} = x_{ij}^k + v_{ij}^{k + 1} $ (22) where
and$i$ are the numbers of particles and parameters in the particles,$j$ and$ x $ represent the position and velocity of the particles,$ v $ ,$ {x^{pb}} $ are the local and global best positions of the particles, q1 and q2 are two random numbers located within [0,2].$ {x^{gb}} $ SCE (shuffled complex evolution) algorithm
-
SCE is an efficient optimization algorithm due to its excellent global search performance and convergence speed, and it is suitable for solving high-dimensional complex nonlinear problems. The main principle of SCE is that each parameter has its own specific search range, within which the fitting value is calculated for each randomly generated individual as well as the ranking. The probability that an individual is selected is:
$ f\left( {{x_k}} \right) = n + 1 - k \quad (k = 1,2,3 \cdots ,n) $ (23) $ p\left( {{x_k}} \right) = \displaystyle\frac{{f\left( {{x_k}} \right)}}{{\displaystyle\mathop \sum \nolimits_{k = 1}^n f\left( {{x_k}} \right)}} = \displaystyle\frac{{2\left( {n + k - 1} \right)}}{{n\left( {n + 1} \right)}} $ (24) where
represents the k-th individual ranked from the lowest to the highest individual fitness value,${x_k}$ is the total number of individuals,$n$ is the function determining the assignment of individuals,$f\left( {{x_k}} \right)$ is the probability that the k-th individual is selected. Individuals are divided into multiple complexes for evolution based on sorting. During evolution, the better adapted particles would replace the less adapted particles as parents to generate the next generation. Then, all groups are mixed and reordered and the process is repeated until the convergence condition is satisfied.$p\left( {{x_k}} \right)$ -
In this section, we analyze the kinetics of wood pyrolysis based on thermogravimetric experimental data obtained at the three different heating rates (5, 10 and 20 K/min), followed by an analytical discussion of the performance of the three algorithms, GA, PSO and SCE, in terms of both efficiency and accuracy in optimizing the kinetic parameters.
Analyses of thermogravimetric results
-
Figure 1 shows the measured mass and MLR curves during wood pyrolysis at different heating rates. Both mass and MLR data are normalized by the initial mass
. As expected, higher heating rates shift the mass and MLR curves toward higher temperature range. In Fig. 1b, it can be reasonably inferred that wood pyrolysis encompasses four main decomposition reactions: water evaporation, decomposition of hemicellulose, cellulose, and lignin, corresponding to the first minor peak before 400 K, the asymmetrical left shoulder of the main peak, the main peak, and the long tail, respectively. In order to evaluate the contribution of each reaction to the total MLR and estimate the relevant kinetics, Gauss multi-peak fitting method[26] is adopted to resolve the MLR curve at each heating rate into multiple elemental curves, as shown in Fig. 2. Each elemental curve corresponds to an elemental pyrolysis reaction. Gauss multi-peak fitting method needs only one parameter when separating overlapped multiple peaks, which renders its simplicity in shape and easy optimization compared with other separation methods, such Weibull, Gauss, Gamma, and Symmetric logistic, as demonstrated by ICTAC Kinetics Committee[25].${m_0}$ Figure 1.
Dependencies of (a) mass and (b) MLR of beech wood on temperature at different heating rates.
Kinetic analyses
-
Pyrolysis reaction scheme of beech wood, including four elemental reactions, is listed in Table 1. Based on previous studies[24], the reaction of water, hemicellulose and cellulose can be described by first order reactions, while lignin is a high-order reaction. Consequently, there are totally 12 parameters, excluding the stoichiometric coefficient of water evaporation, in each optimization run. Therefore, in this section the kinetic analyses are implemented for these four components using the methods introduced previously.
Table 1. Reaction mechanism of beech wood.
# Reaction 1 Water → Vapor 2 Hemicellulose → θ1 Char + (1−θ1) Gas_H 3 Cellulose → θ2 Char + (1−θ2) Gas_C 4 −Lignin → θ3 Char + (1-θ3) Gas_L Based on the separated curves in Fig. 2, the Kissinger method is first used to estimated kinetics of each reaction. Plotting
~${\text{ln}}( {\beta /{T_p^2}} )$ and executing linear fitting, shown in Fig. 3,$1/{T_p}$ and$A$ of each reaction can be estimated based on the slop and intercept, as listed in Table 2. The relatively good linearity of the fitted lines confirms the reliability of the Kissinger method. In addition,${E_a}$ and$A$ are also estimated by the KAS method (${E_a}$ ~${\text{ln}} ( {\beta /{T^2}} )$ ), Tang method ($1/T$ ~$\ln ( {\beta /{T^{1.894661}}} )$ ), and DAEM method ($1/T$ ~${\text{ln}}( {\beta /{T^2}} )$ ). The derived A and$1/T$ are listed in Table 2, and the linear fittings at different conversion rates are depicted in Figs 4−6, respectively.${E_a}$ Figure 3.
Linear fittings of ${\text{ln}}( {\beta /{T_p}^2} )$ vs $1/{T_p}$ in the Kissinger method.
Table 2. Estimated $A$ (s−1) and ${E_a}$ (kJ/mol) of wood pyrolysis by the Kissinger, KAS, Tang, and DAEM methods.
Component Kissinger KAS Tang DAEM Average A Ea A Ea A Ea A Ea A Ea Water 1.49 × 103 44.8 1.01 × 103 42.6 2.61 × 103 42.9 1.53 × 103 42.6 1.72 × 103 42.7 Hemicellulose 1.12 × 1012 147.8 1.30 × 1013 144.7 3.06 × 1013 154.9 9.31 × 1012 144.7 1.77 × 1013 148.1 Cellulose 4.09 × 1012 166.3 4.70 × 1012 174.8 4.84 × 1012 169.7 4.12 × 1012 174.8 4.55 × 1012 173.1 Lignin 2.34 × 1011 180.5 6.40 × 1011 170.7 1.46 × 1012 171.2 7.96 × 1011 170.7 9.66 × 1011 170.9 Figure 4.
Linear fittings of ${\text{ln}}( {\beta /{T^2}} )$ vs $1/T$ plots in the KAS method: (a) water, (b) hemicellulose, (c) cellulose, (d) lignin.
Figure 5.
Linear fittings of $\ln ( {\beta /{T^{1.894661}}} )$ vs $1/T$ plots in the Tang method: (a) water, (b) hemicellulose, (c) cellulose, (d) lignin.
Figure 6.
Linear fittings of ${\text{ln}}( {\beta /{T^2}} )$ vs $1/T$ plots in the DAEM method: (a) water, (b) hemicellulose, (c) cellulose, (d) lignin.
The calculated A and
by the four analytical methods in Table 2 are similar despite some minor deviations. Figure 7 shows the variation trends of calculated${E_a}$ of the four main components in wood with varying conversion rates. Since in the Kissinger and DAEM methods${E_a}$ linearly depends on${\text{ln}}( {\beta /{T^2}} )$ which is similar to those in the KAS and Tang methods, only the results of the KAS and Tang methods are plotted in Fig. 7. In Fig. 7a & c,$1/T$ of water and cellulose decline linearly with increasing α, whereas in Fig. 7b & d opposite variation trend is observed for hemicellulose and lignin. Vyazovkin et al.[25] suggested a pyrolysis process can be described by a single step reaction only if the difference between the maximum and the minimum values of${E_a}$ is smaller than 20% of the average value. In Fig. 7a−d, the ratios of${E_a}$ are 18.6%, 4.01%, 3.38%, 4.15%, implying all these reactions could be regarded as single step reactions.$\left( {{E_{a,max}} - {E_{a,min}}} \right)/{E_{a,ave}}$ Figure 7.
Dependencies of calculated ${E_a}$ on α using the KAS and Tang methods: (a) water, (b) hemicellulose, (c) cellulose, (d) lignin.
Comparison of accuracy of GA, PSO and SCE
-
To compare the optimization performance of the three focused algorithms, the same optimization settings are used for each algorithm. The average values of the four analytical methods in Table 2 are used in determining the search ranges. More specifically, these average values are employed as the mean values of the search ranges during optimization. The lower bounds of search ranges are set to be 0.1 times of the mean values, while the upper bounds are selected to ensure these mean values are the average values of the search ranges. Figure 8 shows the values of objective function (R2) and computation times (tcom) of the optimization runs using the three algorithms with 200-3000 population sizes and fixed iteration number of 200. R2 of GA changes irregularly as the population size is smaller than 1000, but it declines with further increase of population size. This phenomenon indicates the accuracy of GA strongly depends on population size. R2 of PSO descents quickly for population size smaller than 800 and increase slightly beyond this range. While the R2 of SCE are always very low and changes slightly with varying population size, suggesting the accuracy of SCE is very high and is barely affected by population size. In Fig. 8b, tcom of the three algorithms all increase with population size. The difference is that both tcom and its increasing rate of SCE are much larger than those of GA and PSO, implying the computation efficiencies of GA and PSO are approximately identical to each other and both are higher than that of SCE. As introduced in by Shi et al.[24], each iteration of SCE involves multiple complex systems, and therefore the computational complexity is higher than the other two algorithms, requiring a larger amount of computational and storage resources. Meanwhile, Table 3 lists tcom and R2 of the three algorithms with varying population sizes. Distinctly, the computation efficiencies of the three algorithms are ranked as GA ≈ PSO > SCE, while the accuracies are ranked as SCE > PSO > GA.
Figure 8.
Objective function values and the computation times of the three algorithms when optimizing kinetics of wood pyrolysis with 200-3000 population sizes and 200 iterations.
Table 3. Computation times and ${R^2}$ of GA, PSO and SCE optimizations.
Population size GA PSO SCE tcom × 104 (s) R2 × 10−2 tcom × 104 (s) R2 × 10−2 tcom × 104 (s) R2 × 10−2 200 0.16 8.82 0.15 7.86 1.79 7.43 400 0.34 9.09 0.30 7.70 1.82 7.43 600 0.40 8.26 0.49 7.68 2.57 7.43 800 0.55 8.48 0.58 7.47 2.78 7.43 1,000 0.76 8.83 0.75 7.50 3.31 7.43 2,000 1.27 8.30 1.54 7.52 5.46 7.47 3,000 2.26 8.44 2.66 7.67 5.55 7.56 Comparison of efficiency of GA, PSO and SCE
-
Similarly, the convergence efficiencies of the three algorithms are compared using fixed population size of 3000 and 6000 iterations, and the evolutions of R2 are portrayed in Fig. 9. GA, PSO, and SCE converge at 1000, 500, and 800 iterations, respectively. Meanwhile, the decreasing rate of R2 of PSO is higher than the two others. Consequently, it can be reasonably concluded that the convergence efficiencies of three algorithms can be ranked as PSO > SCE > GA. However, in Fig. 8b, the computation efficiencies are ranked as GA ≈ PSO > SCE. In Fig. 8b and Table 3, SCE consumes much more time than GA and PSO, up to approximately 12 times. Apparently, each algorithm has its inherent merits and limits. Overall, PSO is a more favorable algorithm featuring high accuracy, computation efficiency and convergence efficiency. Nevertheless, particular care should be taken when applying it since it may fall into local optimal solution.
Figure 9.
Objective function value evolutions of the three algorithms when optimizing kinetics of wood pyrolysis with 3,000 population size and 6,000 iterations.
Parameter validation
-
In all the optimization runs in above, the minimum values of R2 for GA, PSO, and SCE are 0.0826, 0.0747, and 0.0743, respectively. These optimizations are conducted at 10 K/min heating rate and the corresponding results are listed in Table 4. Apparently, with identical optimization setting, only minor deviations exit among different algorithms. In order verify the reliability of the optimized parameters, the experimental results of 5 and 20 K/min heating rates are predicted using the numerical model and the kinetics listed in Table 4, as exhibited in Figs 10−12 corresponding to GA, PSO, and SCE, respectively. Existing studies[23] revealed that hemicellulose pyrolysis is mainly responsible for the asymmetric shoulder of the MLR curve, cellulose decomposition is related to the main peak, and lignin decomposition is located at a higher temperature range. In Figs 10 & 11, corresponding to GA and PSO algorithms, the numerical curves fit the experimental data well and the relative locations of the subpeaks agree with the literature. Nevertheless, in Fig. 12 the peak temperature of lignin is lower than that of cellulose despite the overall good agreement, conflicting with the existing conclusion. This divergence is presumably caused by the compensation effect among A,
and reaction order when extracting kinetics from TGA data using optimization algorithms, which is still an unsolved problem[25].${E_a}$ Table 4. Best optimized kinetics of wood pyrolysis by GA, PSO and SCE.
Component Parameter Search range GA PSO SCE Water A (s−1) 1.72 × 102−1.72 × 104 1.47 × 104 1.49 × 104 1.49 × 104 Ea (kJ/mol) 4.28 × 104−4.48 × 104 4.4 × 104 4.4 × 104 4.4 × 104 Hemicellulose A (s−1) 9.41 × 1011−9.41 × 1013 5.53 × 1012 4.49 × 1012 4.18 × 1012 Ea (kJ/mol) 1.38 × 105−1.58 × 105 1.4 × 105 1.38 × 105 1.38 × 105 θ 0−0.5 0.28 0.38 0.37 Cellulose A (s−1) 4.32 × 1011−4.32 × 1013 4.09 × 1012 3.95 × 1012 4.09 × 1012 Ea (kJ/mol) 1.6 × 105−1.8 × 105 1.63 × 105 1.63 × 105 1.64 × 105 θ 0−0.5 0.13 0.12 0.12 Lignin A (s−1) 6.0 × 1010−6.0 × 1012 2.22 × 1011 1.96 × 1011 1.15 × 1011 Ea (kJ/mol) 0−3.0 × 105 1.23 × 105 1.16 × 105 1.13 × 105 θ 0−1 0.14 0.01 0.01 n 0−5 4.67 4.99 5 Figure 10.
Comparison between experimental and numerical MLRs using optimized parameters of GA at 5, 10 and 20 K/min heating rates.
-
Pyrolysis of beech wood was investigated experimentally at three heating rates of 5, 10, and 20 K/min. Based on the measured MLR curves, the overlapping peaks of wood was first separated by Gauss multi-peak fitting method to identify their contributions. Then, four analytical methods were used to determine
and A of each reaction. Subsequently, the accuracy, computation efficiency, and convergence efficiency of GA, PSO and SCE algorithms were compared at 10 K/min heating rate. It was found that in terms of optimization accuracy, SCE was the best followed PSO, and then GA. While for computation efficiency, PSO was the best, then GA and SCE. Whereas considering convergence efficiency, the three algorithms were ranked as PSO > SCE > GA. All these indicated each algorithm had its inherent advantages and limits, and PSO featured better overall performance. Furthermore, the reliability of the optimized kinetics was verified by predicting the remaining experimental data at 5 and 10 K/min which were not used during parametrization.${E_a}$ An interesting conclusion attained is that no any single algorithm excels others in all aspects. A potential solution to this issue may be developing more advanced hybrid algorithms which could better balance the accuracy, computation efficiency, convergence efficiency, storage resource, etc. Meanwhile, these heuristic optimization algorithms can also be coupled with some artificial intelligence (AI) algorithms, such as machine learning, deep learning, support vector machines, decision trees, random forest, and metaheuristics. Even though AI has been successfully applied in many engineering applications, few attempts invoking AI have been made to challenge the complex pyrolysis process of biomass. All these need in-depth exploration in future studies.
-
About this article
Cite this article
Wang H, Gong J. 2023. A comparative study of GA, PSO and SCE algorithms for estimating kinetics of biomass pyrolysis. Emergency Management Science and Technology 3:9 doi: 10.48130/EMST-2023-0009
A comparative study of GA, PSO and SCE algorithms for estimating kinetics of biomass pyrolysis
- Received: 13 June 2023
- Accepted: 21 August 2023
- Published online: 08 September 2023
Abstract: Optimization performances of three most frequently utilized optimization algorithms, GA (Genetic Algorithm), PSO (Particle Swarm Optimization), and SCE (Shuffled Complex Evolution), are compared to examine their accuracy, computation efficiency, and convergence efficiency. Micro scale TGA (thermogravimetric analysis) experiments of wood were conducted at three heating rates to collect the necessary data for analysis. Gauss multi-peak fitting method was first applied to identify the contribution of each component of wood to the mass loss rate (MLR) curves. Then the Kissinger method and three isoconversional methods, including KAS, Tang, and DAEM methods, were employed to extract kinetics of wood pyrolysis. The average values of the four sets of solutions were adopted to determine the search range in the following optimizations. A thermally thin numerical model was developed to inversely model the collected experimental data combining the three algorithms. The results showed that wood pyrolysis can be described by a four-component parallel reaction scheme. The four sets of kinetic parameters derived using different analytical methods are very close to each other. When extracting kinetics from experimental data using numerical model and optimization algorithms, the accuracies of the three algorithms are ranked as SCE > PSO > GA. While the computation efficiencies and convergency efficiencies are ranked as GA ≈ PSO > SCE and PSO > SCE > GA, indicating each algorithm has its inherent advantages and limits. In most optimization applications, PSO is more favorable considering its better overall performance.