Determination of hemicellulose, cellulose, and lignin content in different types of biomasses by thermogravimetric analysis and pseudocomponent kinetic model (TGA-PKM Method)

: The standard method for determining the biomass composition, in terms of main lignocellulosic fraction (hemicellulose, cellulose and lignin) contents, is by chemical method; however, it is a slow and expensive methodology, which requires complex techniques and the use of multiple chemical reagents. The main objective of this article is to provide a new e ﬃ cient, low-cost and fast method for the determination of the main lignocellulosic fraction contents of di ﬀ erent types of biomasses from agricultural by-products to softwoods and hardwoods. The method is based on applying deconvolution techniques on the derivative thermogravimetric (DTG) pyrolysis curves obtained by thermogravimetric analysis (TGA) through a kinetic approach based on a pseudocomponent kinetic model (PKM). As a result, the new method (TGA-PKM) provides additional information regarding the ease of carrying out their degradation in comparison with other biomasses. The results obtained show a good agreement between experimental data from analytical procedures and the TGA-PKM method ( ± 7%). This indicates that the TGA-PKM method can be used to have a good estimation of the content of the main lignocellulosic fractions without the need to carry out complex extraction and puriﬁcation chemical treatments. In addition, the good quality of the ﬁt obtained between the model and experimental DTG curves (R 2Adj = 0.99) allows to obtain the characteristic kinetic parameters of each fraction.


Introduction
The use of biomass resources for energy generation has been of considerable importance in recent years [1].The global increase in energy demand has been one of the main reasons for their use.Added to this situation, there is also a need for dealing with certain problems, such as the depletion of fossil fuel reserves and the increase in environmental pollution from the use of these energy sources [2].
In this context, biomass has the advantage of being the only renewable resource that can be used in solid, liquid and gaseous forms [3].Furthermore, biomass has the great capacity of producing by-products of high interest, such as catalytic carbons [4] and bioplastics [5].However, biomass has a number of features that make it difficult to use, including its moisture content, low-energy density and complex structure.Lignocellulosic biomass is made up of a structure that includes mainly cellulose, hemicellulose and lignin [6,7].The proportions and distribution of these components in the biomass physical structure is complex and depends on the type of species.The knowledge of this composition is very important for its use in different industrial applications.
Up to now, the determination of biomass composition, in terms of hemicellulose, cellulose and lignin contents, has been made by the chemical method.However, it is a slow and expensive methodology, which requires complex techniques and the use of multiple chemical reagents [8].This means that it is not a suitable method for use in industrial applications.
Thermogravimetric analysis and, especially, the derivative thermogravimetric (DTG) curve is often used for the preliminary study of various thermochemical processes with biomass, since it allows the determination of the different stages of biomass devolatilization.In general, the process of devolatilization of the biomass in the absence of oxygen usually differentiates four stages corresponding to the loss of moisture and the three lignocellulosic components (hemicellulose, cellulose and lignin) [3].Numerous articles have been published in which the thermal decomposition intervals of these lignocellulosic components are presented based on the deconvolution of the DTG curves [9][10][11][12][13][14][15].It has been observed that, after moisture removal that takes place up to 150 • C, the decomposition of the three biomass lignocellulosic components takes place: hemicellulose is the first component to decompose between 200-300 • C, followed by cellulose between 250-380 • C. Regarding the thermal decomposition of lignin, it is the component with the most complex structure, and its decomposition range is the widest [16], occurring from 200 • C up to high temperatures such as 1000 • C [17,18].
There are different studies based on determining the lignocellulosic composition by analyzing DTG curves.However, most of these studies are based exclusively on applying deconvolution methods without taking into account their kinetic interpretation of the process [3,19].
On the other hand, kinetic studies on the thermal decomposition of biomass are extensive, in which the use of different kinetic models is analyzed [20,21], providing the kinetic parameters that best fit the experimental data.However, these studies do not focus on finding a method that allows the quantification of the three main lignocellulosic fractions of the biomass.
The use of kinetic analysis to the quantification of the main lignocellulosic fractions allows to include restrictions for a more precise quantification, while a physical interpretation is added to the deconvolution process.
The main objective of this work is to provide a new efficient, low-cost and fast method for the determination of the hemicellulose, cellulose and lignin contents of different types of biomasses, from agricultural by-products to wood.The method is based on applying deconvolution techniques on DTG pyrolysis curves based on a kinetic analysis of the process, and the kinetic model used is based on the assumption that the degradation of each lignocellulosic fraction can be represented by the evolution of a certain number of pseudocomponents.

Biomass Samples
Five raw materials representing different types of biomass have been selected, including agricultural biomass (wheat straw) and forest biomass, both as softwood barks (spruce bark and pine bark) and hardwoods (poplar and willow).
The pine bark originated from Sweden, while the other biomasses (wheat straw, poplar, spruce bark and willow) came from the South of France.

Experimental Method
Each sample was crushed in a mill (Model A 10 basic, IKA-Werke GmbH & Co. KG, Staufen, Germany) and then sieved.The sample sizes were all less than 100 µm in order to minimize the heat transfer resistances and mass transfer diffusion effects.
The TG (thermogravimetric) analysis was performed on a TG-DTA analyzer (Model DTG-60H, SHIMADZU Co. Ltd., Kyoto, Japan).The analyses were carried out using a nitrogen atmosphere with a flow rate of 50 mL min −1 .The heating rate used was 5 • C min −1 , from room temperature to a final temperature of 1000 • C. The sample weight was c.a. 10 mg.
To reduce temperature-related errors, the equipment used was calibrated across the entire temperature range.In addition, the actual sample temperature was used directly to solve the kinetic equations and to calculate the actual sample heating rate [22].
The information obtained in these analyses was the weight loss as the temperature and time of analysis increase (TG curve).

Data Treatment
The TG analysis provides the weight loss as a function of temperature over time.The analysis can be used to determine the different fractions of volatiles released as a function of temperature, as well as the solid residue remaining after heat treatment.However, for the determination of kinetics, it is more useful to use the derivative thermogravimetric (DTG) of weight loss as a function of time, because this signal is much more sensitive to small changes.
Before proceeding with its calculation, it is necessary to preprocess the data in order to obtain a curve that depends exclusively on the process variables.
The first step is the normalization of the TG signal.The normalization has been carried out in relation to the initial weight of the sample (m 0 ) and the final weight (m ∞ ) of the sample.To do this, the weight fraction of the volatiles remaining in the sample has been calculated for each instant of discrete time i, as indicated in Equation (1).
In this case, m ∞ represents the mass of char obtained at the end of each TG analysis and includes the mass of ash and fixed carbon at the final temperature of the analysis.

DTG Curves
The DTG curve is obtained from the weight over time derivative for each experimental point, i.e., where ∆ is the interval of the experimental data taken into account.In this case, ∆ = 1 has been used.

Kinetic Model
The thermochemical decomposition of the biomass can be represented by three main kinetics that correspond to the degradation of hemicellulose, cellulose and lignin.The most commonly used model consists of assuming that the process can be represented by the decomposition reactions of each of these compounds [23,24].In addition, the decomposition of these compounds can be represented by a number of parallel and independent first-order Arrhenius-type reactions, named pseudocomponents.
Thus, for the adjustment of the DTG curve of each biomass, it has been assumed that the process follows the model that consists of the decomposition of hemicellulose, cellulose and lignin independently, so that the overall kinetics can then be expressed as follows: where H, C and L represent the mass fraction of hemicellulose, cellulose and lignin, respectively.At the same time, the kinetics of each of these fractions can be represented by a set of parallel reactions, expressed in the form: where T: temperature, in K; R: ideal gas constant, 8.314 × 10 −3 kJ (K mol) −1 ; j: number of pseudocomponents of the fractions of hemicellulose, cellulose and lignin, which take the values from 1 to the total number of pseudocomponents of each fraction of hemicellulose; cellulose and lignin (mH, mC and mL); K Hj , K Cj and K Lj : pre-exponential factors of the pseudocomponents of the hemicellulose, cellulose and lignin fractions, expressed in s −1 and E Hj , E Cj and E Lj : activation energies of the pseudocomponents of the hemicellulose, cellulose and lignin fractions, expressed in kJ mol −1 .
In general, the kinetic equation of each pseudocomponent j, corresponding to fraction F (F = H, C, L), in a nonisothermal process at constant heating rate β = dT/dt, is given by The integral of the second term can be resolved by using the exponential integral, defined as follows: Thus, Equation (7), integrated between To and T, can be expressed in the form Therefore, the kinetics of each pseudocomponent depends on three variables: the pre-exponential factor, the activation energy and the initial concentration of the pseudocomponent in the biomass (X Fj,0 ).
A restriction that the system must satisfy is that the sum of the mass fractions of all the pseudocomponents must be equal to the mass fraction of all volatiles generated for each instant of time t = i.
Combining Equations ( 9) and (10) for each instant of discrete time i gives a system of equations with 3 × (m H + m C + m L ) − 1 unknowns, which needs to be solved.

Calculation Procedure
For the calculation of unknown variables, an optimization method based on the minimization by least squares has been used.As an objective function (OF), the square of the errors between the values of the experimental curve and the model has been used for each instant of time i, in which the model has been evaluated.
The solution has been made with MATLAB using the lsqcurvefit command to find the constants that best fit the system of equations.The final solution was obtained when the percentage variation of the OF was less than 0.01% during five consecutive cycles of 200 iterations each (∆OF 5 < 0.01%).
The obtained quality of fit (QOF) between the simulated and experimental curves was evaluated with the expression (12).
where n is the number of experimental points employed (967).
Additionally, the goodness of fit was evaluated by the adjusted R-squared, R 2 Adj , which represents the response that is explained by the model and was calculated as the ratio between the sum of square of the residuals (SSE) and the total sum of squares (SST) as follows [25]: where k is the number of variables.
The initial values of the constants were taken after an initial analysis of the kinetics, using as initial seed values the restrictions on the concentrations of the hemicellulose, cellulose and lignin fractions obtained from the literature review (Table 1).The decision tree of the calculation process is as Figure 1: Processes 2020, 8, x FOR PEER REVIEW 6 of 22

Analytical Method
The lignocellulosic biomass wt.% composition was determined by chemical methods by the VTT and TECNALIA laboratories; the detailed procedure was described in [32].Biomasses were previously sampled and prepared through TAPPI T257 and then conditioned through TAPPI 264.Table 2 includes the analytical results obtained.The results obtained in Table2 are in-line with the results obtained by other researchers [33].According to the literature, the softwood bark composition corresponds to a cellulose content of 18-38%, the hemicellulose content is 15-33% and the lignin content is 30-60%.For hardwood biomasses, the cellulose content is 43-47%, the hemicellulose content is 25-35% and the lignin content is 16-24%.Finally, the composition of herbaceous biomass, such as cereal straw, is 33-38% cellulose, 26-32% hemicellulose and 17-19% lignin.

Analytical Method
The lignocellulosic biomass wt.% composition was determined by chemical methods by the VTT and TECNALIA laboratories; the detailed procedure was described in [31].Biomasses were previously sampled and prepared through TAPPI T257 and then conditioned through TAPPI 264.Table 2 includes the analytical results obtained.The results obtained in Table 2 are in-line with the results obtained by other researchers [32].According to the literature, the softwood bark composition corresponds to a cellulose content of 18-38%, the hemicellulose content is 15-33% and the lignin content is 30-60%.For hardwood biomasses, the cellulose content is 43-47%, the hemicellulose content is 25-35% and the lignin content is 16-24%.Finally, the composition of herbaceous biomass, such as cereal straw, is 33-38% cellulose, 26-32% hemicellulose and 17-19% lignin.
Therefore, according to the literature review [32,33] and the analyses carried out (Table 2), softwood bark has higher lignin content than hardwood and agricultural biomasses.On the other hand, hardwood has a higher cellulose content than the rest of the biomasses analyzed.
It also should be noted that, during the thermogravimetric analysis (TGA), it is possible to differentiate the biomass into its three main lignocellulosic fractions, but it is not possible to distinguish the extractives from the other fractions.Extractives are a group of compounds that can be obtained from the biomass using organic solvents, such as benzene, alcohol or water [34].The main components of the lipophilic extracts are triglycerides, fatty acids, resin acids, sterile esters and sterols and of hydrophilic extracts are lignin [35].The extractives thermally degrade in the temperature range of 200-400 • C, which falls within the range in which hemicellulose and cellulose and, also, lignin is degraded.For this reason, in order to get comparable results with those obtained by the TGA method, the analytical data are expressed in weight % on a dry and ash and extractives-free basis.

Devolatilization Behavior
The performance of the DTG curves shows similar behavior (Figure 2).At first sight, two large peaks can be observed in all of them: the first one appears from room temperature to about 150 • C and corresponds to the loss of moisture.At temperatures exceeding 150 • C, degradation of lignocellulosic compounds begins [30,32,36,37].The second large peak is located in the range of temperature between 250 and 380 • C and corresponds to the degradation of cellulose.Two other peaks, which are more or less perceptible depending on the type of biomass, can be seen overlapping the cellulose peak.Thus, at temperatures between 200 and 300 • C, the degradation of hemicellulose occurs, which proves a deformation of the cellulose peak in that temperature range.Finally, lignin is the component with the most complex structure, and its decomposition range is the widest, occurring from 200 • C to the final temperature of the analysis.The degradation of lignin is more significant near the 400 • C zone, where a small peak can be observed that overlaps with the end of the cellulose degradation.
Processes 2020, 8, x FOR PEER REVIEW 7 of 22 Therefore, according to the literature review [32,33] and the analyses carried out (Table 2), softwood bark has higher lignin content than hardwood and agricultural biomasses.On the other hand, hardwood has a higher cellulose content than the rest of the biomasses analyzed.
It also should be noted that, during the thermogravimetric analysis (TGA), it is possible to differentiate the biomass into its three main lignocellulosic fractions, but it is not possible to distinguish the extractives from the other fractions.Extractives are a group of compounds that can be obtained from the biomass using organic solvents, such as benzene, alcohol or water [34].The main components of the lipophilic extracts are triglycerides, fatty acids, resin acids, sterile esters and sterols and of hydrophilic extracts are lignin [35].The extractives thermally degrade in the temperature range of 200-400 °C, which falls within the range in which hemicellulose and cellulose and, also, lignin is degraded.For this reason, in order to get comparable results with those obtained by the TGA method, the analytical data are expressed in weight % on a dry and ash and extractives-free basis.

Devolatilization Behavior
The performance of the DTG curves shows similar behavior (Figure 2).At first sight, two large peaks can be observed in all of them: the first one appears from room temperature to about 150 °C and corresponds to the loss of moisture.At temperatures exceeding 150 °C, degradation of lignocellulosic compounds begins [30,32,36,37].The second large peak is located in the range of temperature between 250 and 380 °C and corresponds to the degradation of cellulose.Two other peaks, which are more or less perceptible depending on the type of biomass, can be seen overlapping the cellulose peak.Thus, at temperatures between 200 and 300 °C, the degradation of hemicellulose occurs, which proves a deformation of the cellulose peak in that temperature range.Finally, lignin is the component with the most complex structure, and its decomposition range is the widest, occurring from 200 °C to the final temperature of the analysis.The degradation of lignin is more significant near the 400 °C zone, where a small peak can be observed that overlaps with the end of the cellulose degradation.In relation to the development of each type of biomass, it is observed that pine bark and spruce bark have very similar development patterns.Both barks, as compared to the rest of the biomasses (poplar, willow and white straw), have a higher peak near 400 • C corresponding to the degradation of lignin and a lower peak height corresponding to the degradation of cellulose (~350 • C) and hemicellulose (~300 • C).Therefore, these softwood barks have a higher lignin content and lower cellulose and hemicellulose contents, as compared to other biomasses (Table 2).It is also observed that these two biomasses have the lowest DTG area, so they are the ones that release the least amounts of total volatiles.
On the other hand, willow and poplar show very similar behaviors, which indicates that their compositions will be very similar.Both biomasses present a greater generation of volatiles in the cellulose degradation zone.This is in agreement with the fact that both biomasses have higher cellulose contents and lower lignin contents compared to the rest of the biomasses analyzed (Table 2).
Finally, wheat straw presents a single peak in the degradation zone of hemicellulose and cellulose and is slightly displaced to the low temperature zone.This suggests a higher hemicellulose content, while the evolution of the lignin content is very similar to that of poplar and willow.

TGA-PKM Method
The first step was to determine the minimum number of pseudocomponents needed to adequately represent the evolution of each of the three main lignocellulosic fractions and all volatiles generated during the thermal degradation process.
This analysis was carried out by means of an initial kinetic analysis, in which a division of the DTG was established according to the degradation temperatures of the three main constituents of the biomass (hemicellulose, cellulose and lignin), in addition to water.Each of these regions was initially attributed a single pseudocomponent; then, the number of pseudocomponents was gradually increased, until an adequate performance of the evolution of the volatiles was achieved.The minimum numbers of pseudocomponents necessary for the quantifications of each fraction are shown in the Table 3.The use of a larger number of pseudocomponents could induce overfitting.The next step was to determine the minimum number of heating rates needed to achieve the objective of quantifying the main biomass fractions.The use of three or more heating rates while reducing the effect of kinetic compensation and improving the accuracy of kinetic parameters requires the use of significantly different heating rates, which involves, in practice, the use of higher heating rates.However, higher heating rates worsen the separation of lignocellulosic fractions, making their identification more difficult.Additionally, the use of various heating rates for the quantification of the lignocellulosic fractions is more time-consuming.
Therefore, a low heating rate achieves a better separation of the degraded compounds and is less time-consuming.This is the reason why a single heating rate of 5 • C min −1 has been employed in the determination of the main lignocellulosic fractions.However, a validation of the method using three heating rates has been carried out and is reported in Section 3.4.
To improve the accuracy of the kinetic parameters, it was found that the use of upper and lower limits of the kinetic parameters (Tables 4 and 5) was necessary, not only to ensure adequate values of the pre-exponential and activation energy but, also, to provide adequate seed values for the determination of the hemicellulose, cellulose and lignin fractions.Finally, taking into account the above procedure, the values of the kinetic parameters of each pseudocomponent were calculated by the TGA-PKM method and are summarized in Table 6.
Overall, the results obtained in this study are in reasonable ranges when compared to the results corresponding to the kinetics of other biomasses published, as can be seen in Table 7.As shown in Tables 6 and 7, cellulose is the compound with the highest activation energies.This is attributed to the fact that the cellulose is a very long polymer of glucose units without any branches [18], while hemicellulose has a random branched amorphous structure that gives a lower activation energy; this is the reason why hemicellulose decomposes more easily in a lower temperature range [38].
Lignin has a very complex structure composed of three kinds of heavily crosslinked phenylpropane structures [18].Additionally, it is observed that the activation energy is lower than for hemicellulose and cellulose, which indicates that its thermal degradation is easier.However, it presents much lower values of pre-exponential factors that cause a lower reaction rate; this fact is reflected in the wide range of temperatures in which its degradation takes place and in the high temperature required to reach a complete degradation.
In addition, Figures 3-7 show the fit of the model to the DTG experimental data, as well as the contribution of the different pseudocomponents to the model.In all the figures, it can be seen that a good fit is achieved between the global model, obtained as the envelope resulting from the sum of the seven pseudocomponents, and the experimental DTG curve.By comparison, between the kinetic constants in Table 6 and Figures 3-7, it can be seen that low activation energy leads to a reaction in the low temperature zone and vice versa.With respect to the factor, low values cause the reaction rate to be slower and to take place over a wider temperature range, which is characteristic of the lignin pseudocomponents.On the contrary, high values of the pre-exponential factor increase the reaction rate, leading to a narrower temperature range, which is characteristic of cellulose, for example.
On the other hand, at the same activation energy, a higher pre-exponential factor causes the By comparison, between the kinetic constants in Table 6 and Figures 3-7, it can be seen that low activation energy leads to a reaction in the low temperature zone and vice versa.With respect to the pre-exponential factor, low values cause the reaction rate to be slower and to take place over a wider temperature range, which is characteristic of the lignin pseudocomponents.On the contrary, high values of the pre-exponential factor increase the reaction rate, leading to a narrower temperature range, which is characteristic of cellulose, for example.
On the other hand, at the same activation energy, a higher pre-exponential factor causes the reaction to take place in the high temperature zone.For example, there are lignin pseudocomponents with a similar activation energy as hemicellulose pseudocomponents (Table 6) but with much lower pre-exponential factors, which cause the reaction to take place at higher temperatures.
The quality of the fit expressed as R 2 and QOF% can be observed in Table 8.  9) and (10).As can be seen, the TG curve model achieves good results not only with respect to the model fitting to the experimental TG curve along the operating temperature but, also, with respect to the final value.Table 9 shows the comparison between the analytical composition and the data obtained with the TGA-PKM method.As can be seen, there is a good agreement between the data obtained through analytical procedures and the TGA-PKM model.This indicates that the new method can be used to have a good estimation of the content of the main lignocellulosic fractions of the analyzed biomasses without the need to carry out complex extraction and purification chemical treatments.Table 9 shows the comparison between the analytical composition and the data obtained with the TGA-PKM method.As can be seen, there is a good agreement between the data obtained through analytical procedures and the TGA-PKM model.This indicates that the new method can be used to have a good estimation of the content of the main lignocellulosic fractions of the analyzed biomasses without the need to carry out complex extraction and purification chemical treatments.The following error ranges are obtained between the values measured analytically and those measured by the TGA-PKM method for each of the main lignocellulosic fractions: hemicellulose (−2.56-6.78),cellulose (−6.74-7.14)and lignin (−6.99-6.73).The level of accuracy achieved is considered suitable, taking into account that it is within the error range of the chemical methods.For example, Korpinen et al. found that the determination of lignin by different chemical methods can be as high as 10 wt.% [39]; Ioelovich [40] also determined a difference of 4 wt.% between the TAPPI and NERL methods in determination of the cellulose content.In this way, the TGA-PKM method allows to obtain a fast estimation of the contents of the main lignocellulosic fractions within the ranges that would be obtained by a chemical analysis.

Validation of the TGA-PKM Method
In order to check the validity of the method, an additional fit of the poplar biomass devolatilization was performed using, simultaneously, three heating rates: 3, 5 and 10 • C min −1 datasets.
Figure 13 shows the graphical results by fitting the model to the DTG and TG curves for each heating rate.
Additionally, the quality of the fit achieved for each heating rate and for the global fit are summarized in Table 10, where QOF% and R 2 Adj are shown for each heating rate dataset and for the three heating rates simultaneously.The results obtained (Figure 13 and Table 10) indicate that the quality of the fit obtained is very satisfactory, since the model is capable of representing the evolution of the devolatilization process when different heating rates are used.
NERL methods in determination of the cellulose content.In this way, the TGA-PKM method allows to obtain a fast estimation of the contents of the main lignocellulosic fractions within the ranges that would be obtained by a chemical analysis.

Validation of the TGA-PKM Method
In order to check the validity of the method, an additional fit of the poplar biomass devolatilization was performed using, simultaneously, three heating rates: 3, 5 and 10 °C min   Finally, the results obtained in the determination of the lignocellulosic fractions are shown in Table 12.The results obtained with a single heating rate are comparable to those obtained with three heating rates, because the deviation between the results calculated by the model and by the analytical method are of the same order when a single heating rate or three heating rates are considered.However, slightly better results are achieved if a single heating rate of 5 • C min −1 is used, but mainly, it requires considerably less analysis time, which justifies the use of a single heating rate.

Conclusions
Five lignocellulosic samples have been characterized by the TGA-PKM experimental protocol, covering different types of woody and herbaceous biomasses from both forest and agricultural origins (spruce bark, pine bark, poplar, willow and wheat straw).
The TGA-PKM method developed allows the determination of the main lignocellulosic fractions of biomasses without the need to use long and complex chemical methods; e.g., TAPPI methods T222 and T249 require several long successive steps (hydrolysis, extraction, filtration, neutralization, reduction, etc.) [41], which may require several days of work in the laboratory, while the new method may be performed in a few hours.Thus, it would be possible to reduce the cost of analysis and processing time by 80-90%.
The accuracy of the TGA-PKM method was tested and proved to be significantly good and consistent within the order of magnitude of the standard analytical methods to determine the contents of the main lignocellulosic fractions.

Figure 1 .
Figure 1.Decision tree of the calculation procedure..

Figure 1 .
Figure 1.Decision tree of the calculation procedure.

Figure 3 .
Figure 3. Model fitted to the experimental pine bark DTG curve.

Figure 3 .
Figure 3. Model fitted to the experimental pine bark DTG curve.

Figure 3 .
Figure 3. Model fitted to the experimental pine bark DTG curve.

Figure 4 .
Figure 4. Model fitted to the experimental spruce bark DTG curve.

Figure 4 .
Figure 4. Model fitted to the experimental spruce bark DTG curve.Processes 2020, 8, x FOR PEER REVIEW 12 of 22

Figure 5 .
Figure 5. Model fitted to the experimental poplar DTG curve.

Figure 5 .
Figure 5. Model fitted to the experimental poplar DTG curve.

Figure 5 .
Figure 5. Model fitted to the experimental poplar DTG curve.

Figure 7 .
Figure 7. Model fitted to the experimental wheat straw DTG curve.

Figure 7 .
Figure 7. Model fitted to the experimental wheat straw DTG curve.

Figure 8 .
Figure 8. Model fitted to the experimental pine bark TG curve.

Figure 8 .
Figure 8. Model fitted to the experimental pine bark TG curve.

Figure 8 .
Figure 8. Model fitted to the experimental pine bark TG curve.

Figure 9 .
Figure 9. Model fitted to the experimental spruce bark TG curve.

Figure 9 .
Figure 9. Model fitted to the experimental spruce bark TG curve.Processes 2020, 8, x FOR PEER REVIEW 15 of 22

Figure 10 .
Figure 10.Model fitted to the experimental poplar TG curve.

9 Figure 10 .
Figure 10.Model fitted to the experimental poplar TG curve.

Figure 10 .
Figure 10.Model fitted to the experimental poplar TG curve.

Figure 11 .
Figure 11.Model fitted to the experimental willow TG curve.

Figure 11 .
Figure 11.Model fitted to the experimental willow TG curve.Processes 2020, 8, x FOR PEER REVIEW 16 of 22

Figure 12 .
Figure 12.Model fitted to the experimental wheat straw TG curve.

Figure 12 .
Figure 12.Model fitted to the experimental wheat straw TG curve.

Table 1 .
Literature references of the main lignocellulosic fraction compositions related to used biomasses.
a By chemical methods, b by thermogravimetric analysis (TGA) and c cellulose as glucan and hemicellulose as xylan.

Table 2 .
Composition by chemical methods for the raw biomasses (wt.%, dry basis).

Table 2 .
Composition by chemical methods for the raw biomasses (wt.%, dry basis).

Table 3 .
Minimum number of components for each biomass fraction.

Table 6 .
Kinetic parameters of the pseudocomponents.

Table 7 .
Kinetic parameters from other studies.

Table 8 .
Quality of the fit expressed as R 2 Adj and QOF%.Figures 8-12 show the fit of the global model to the TG experimental data.The TG curve model has been obtained simultaneously with the DTG curve model by solving Equations (

Table 9 .
Comparison between the analytical and thermogravimetric analysis-pseudocomponent

Table 10 .
Quality of the fit for each dataset.

Table 11 .
Kinetic parameters of each pseudocomponents calculated using a single heating rate and three simultaneous heating rates.

Table 12 .
Comparison between TGA-PKM results using three simultaneous heating rates and the analytical method.