Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Study on the Law of Diesel Oil Carrying Water in Lanzhou–Chengdu–Chongqing Product Oil Pipeline Based on Large Eddy Simulation
Next Article in Special Issue
Scale-Up Cultivation of Phaeodactylum tricornutum to Produce Biocrude by Hydrothermal Liquefaction
Previous Article in Journal
Exploring E-Waste Resources Recovery in Household Solid Waste Recycling
Previous Article in Special Issue
Environmental Assessment of Olive Mill Solid Waste Valorization via Anaerobic Digestion Versus Olive Pomace Oil Extraction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Hemicellulose, Cellulose, and Lignin Content in Different Types of Biomasses by Thermogravimetric Analysis and Pseudocomponent Kinetic Model (TGA-PKM Method)

1
CARTIF Centre of Technology, Parque Tecnológico de Boecillo, 205, Boecillo, 47151 Valladolid, Spain
2
ITAP Institute, University of Valladolid, Paseodel Cauce 59, 47011 Valladolid, Spain
3
TECNALIA, Basque Research and Technology Alliance (BRTA), Área Anardi 5, E-20730 Azpeitia, Spain
4
VTT-Technical Research Centre of Finland, P.O. Box 1000, VTT, FI-02044 Espoo, Finland
*
Author to whom correspondence should be addressed.
Processes 2020, 8(9), 1048; https://doi.org/10.3390/pr8091048
Submission received: 30 July 2020 / Revised: 20 August 2020 / Accepted: 21 August 2020 / Published: 27 August 2020

Abstract

:
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 efficient, low-cost and fast method for the determination of the main lignocellulosic fraction contents of different 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 purification chemical treatments. In addition, the good quality of the fit obtained between the model and experimental DTG curves (R2Adj = 0.99) allows to obtain the characteristic kinetic parameters of each fraction.

Graphical Abstract

1. 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.

2. Materials and Methods

2.1. 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.

2.2. 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).

2.3. 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 (m0) 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).
X i = m i m m 0 m
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.

2.4. DTG Curves

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

2.5. 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:
d X d t = d X H d t + d X C d t + d X L d t
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:
d X H d t = j = 1 m H d X H j d t = j = 1 m H K H j e x p ( E H j R T ) X H j
d X C d t = j = 1 m C d X C j d t = j = 1 m C K C j e x p ( E C j R T ) X C j
d X L d t = j = 1 m L d X L j d t = j = 1 m L K L j e x p ( E L j R T ) X L j
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); KHj, KCj and KLj: pre-exponential factors of the pseudocomponents of the hemicellulose, cellulose and lignin fractions, expressed in s−1 and EHj, ECj and ELj: 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
d X F j X F j = K F j β e x p ( E F j R T ) d T
The integral of the second term can be resolved by using the exponential integral, defined as follows:
u e u u d u ,   u = E R
Thus, Equation (7), integrated between To and T, can be expressed in the form
X F j , i = X F j , 0 . e x p { K F j β [ T i . e x p ( E F j R T i ) E F j / R T i e x p ( E F j R T ) T d T ] }
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 (XFj,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.
X i = X H i + X C i + X L i = j = 1 m H X H j , i + j = 1 m C X C j , i + j = 1 m L X L j , i
Combining Equations (9) and (10) for each instant of discrete time i gives a system of equations with 3 × (mH + mC + mL) − 1 unknowns, which needs to be solved.

2.6. 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.
O . F . = i = 1 n [ ( d X d t ) i , e x p ( d X d t ) i , m o d e l ] 2
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 (ΔOF5 < 0.01%).
The obtained quality of fit (QOF) between the simulated and experimental curves was evaluated with the expression (12).
Q O F   ( % ) = 100   x   i = 1 n [ ( d X d t ) i , e x p ( d X d t ) i , m o d e l ] 2 / n max [ ( d X d t ) i , e x p ]
where n is the number of experimental points employed (967).
Additionally, the goodness of fit was evaluated by the adjusted R-squared, R2Adj, 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]:
R a d j 2 = 1 ( n 1 ) x S S E ( n ( k + 1 ) ) x S S T = 1 ( n 1 ) x i = 1 n [ ( d X d t ) i , e x p ( d X d t ) i , m o d e l ] 2 ( n ( k + 1 ) ) x i = 1 n [ ( d X d t ) i , e x p ( d X d t ) i , e x p ¯ ] 2
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:

3. Results and Discussion

3.1. 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.

3.2. 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.

3.3. 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 (Table 4 and Table 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 Table 6 and Table 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, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 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 Figure 3, Figure 4, Figure 5, Figure 6 and Figure 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 R2 and QOF% can be observed in Table 8.
In addition, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 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 (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.
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.

3.4. 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 R2Adj 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.
Table 11 shows the kinetic parameters by fitting the model using a single heating rate and three heating rates simultaneously. The obtained results by both datasets are very similar. For example, the activation energy obtained is identical for almost all the pseudocomponents, and only pseudocomponents 3 and 7 have a relative standard deviation of 4%.
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.

4. 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.

Author Contributions

Conceptualization, D.D. and A.U.; methodology, D.D. and A.U.; software, D.D. and A.U.; validation, D.D., A.U., R.P.; formal analysis, D.D., A.U., A.B. and T.T.; investigation, D.D. and A.U.; resources, A.B., T.T. and R.P.; data curation, D.D., A.U.; writing (original draft preparation), D.D., A.U.; writing (review and editing), R.P, A.B. and T.T.; visualization, D.D., A.U., R.P., A.B. and T.T.; supervision, D.D., A.U., R.P., A.B. and T.T.; project administration, A.B; funding acquisition, A.B., T.T. and R.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union’s Horizon 2020 research and innovation programme under grant agreement No 723670, with the title “Systemic approach to reduce energy demand and CO2 emissions of processes that transform agroforestry waste into high added value products (REHAP)”.

Acknowledgments

The authors would like to thank María González Martínez from IMT Mines Albi (Université de Tolousse) for her technical contribution and support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, S.; Dale, B.E. All biomass is local: The cost, volume produced, and global warming impact of cellulosic biofuels depend strongly on logistics and local conditions. Biofuels Bioprod. Biorefining 2015, 9, 422–434. [Google Scholar] [CrossRef]
  2. McKendry, P. Energy production from biomass (part 1): Overview of biomass. Bioresour. Technol. 2002, 83, 37–46. [Google Scholar] [CrossRef]
  3. Rego, F.; Dias, A.P.S.; Casquilho, M.; Rosa, F.C.; Rodrigues, A. Fast determination of lignocellulosic composition of poplar biomass by thermogravimetry. Biomass Bioenergy 2019, 122, 375–380. [Google Scholar] [CrossRef]
  4. Álvarez-Mateos, P.; Alés-Álvarez, F.J.; García-Martín, J.F. Phytoremediation of highly contaminated mining soils by Jatropha curcas L. and production of catalytic carbons from the generated biomass. J. Environ. Manag. 2019, 231, 886–895. [Google Scholar] [CrossRef] [PubMed]
  5. Snell, K.D.; Peoples, O.P. PHA bioplastic: A value-added coproduct for biomass biorefineries. Biofuels Bioprod. Biorefining Innov. Sustain. Econ. 2009, 3, 456–467. [Google Scholar] [CrossRef]
  6. Pereira, B.L.C.; Carneiro, A.D.C.O.; Carvalho, A.M.M.L.; Colodette, J.L.; Oliveira, A.C.; Fontes, M.P.F. Influence of chemical composition of Eucalyptus wood on gravimetric yield and charcoal properties. BioResources 2013, 8, 4574–4592. [Google Scholar] [CrossRef]
  7. Melzer, M.; Blin, J.; Bensakhria, A.; Valette, J.; Broust, F. Pyrolysis of extractive rich agroindustrial residues. J. Anal. Appl. Pyrolysis 2013, 104, 448–460. [Google Scholar] [CrossRef]
  8. Park, J.I.; Liu, L.; Ye, X.P.; Jeong, M.K.; Jeong, Y.S. Improved prediction of biomass composition for switchgrass using reproducing kernel methods with wavelet compressed FT-NIR spectra. Expert Syst. Appl. 2012, 39, 1555–1564. [Google Scholar] [CrossRef]
  9. Yu, J.; Paterson, N.; Blamey, J.; Millan, M. Cellulose, xylan and lignin interactions during pyrolysis of lignocellulosic biomass. Fuel 2017, 191, 140–149. [Google Scholar] [CrossRef]
  10. Shen, D.; Xiao, R.; Gu, S.; Luo, K. The pyrolytic behavior of cellulose in lignocellulosic biomass: A review. RSC Adv. 2011, 1, 1641–1660. [Google Scholar] [CrossRef]
  11. Shen, D.K.; Gu, S.; Bridgwater, A.V. The thermal performance of the polysaccharides extracted from hardwood: Cellulose and hemicellulose. Carbohydr. Polym. 2010, 82, 39–45. [Google Scholar] [CrossRef]
  12. Shen, D.K.; Gu, S. The mechanism for thermal decomposition of cellulose and its main products. Bioresour. Technol. 2009, 100, 6496–6504. [Google Scholar] [CrossRef] [PubMed]
  13. Li, S.; Lyons-Hart, J.; Banyasz, J.; Shafer, K. Real-time evolved gas analysis by FTIR method: An experimental study of cellulose pyrolysis. Fuel 2001, 80, 1809–1817. [Google Scholar] [CrossRef]
  14. Qiao, Y.; Wang, B.; Ji, Y.; Xu, F.; Zong, P.; Zhang, J.; Tian, Y. Thermal decomposition of castor oil, corn starch, soy protein, lignin, xylan, and cellulose during fast pyrolysis. Bioresour. Technol. 2019, 278, 287–295. [Google Scholar] [CrossRef]
  15. Wang, S.; Lin, H.; Ru, B.; Sun, W.; Wang, Y.; Luo, Z. Comparison of the pyrolysis behavior of pyrolytic lignin and milled wood lignin by using TG–FTIR analysis. J. Anal. Appl. Pyrolysis 2014, 108, 78–85. [Google Scholar] [CrossRef]
  16. Di Blasi, C. Modeling chemical and physical processes of wood and biomass pyrolysis. Prog. Energy Combust. Sci. 2008, 34, 47–90. [Google Scholar] [CrossRef]
  17. Zhou, H.; Long, Y.; Meng, A.; Li, Q.; Zhang, Y. The pyrolysis simulation of five biomass species by hemi-cellulose, cellulose and lignin based on thermogravimetric curves. Thermochim. Acta 2013, 566, 36–43. [Google Scholar] [CrossRef]
  18. Skreiberg, A.; Skreiberg, Ø.; Sandquist, J.; Sørum, L. TGA and macro-TGA characterisation of biomass fuels and fuel mixtures. Fuel 2011, 90, 2182–2197. [Google Scholar] [CrossRef]
  19. Carrier, M.; Loppinet-Serani, A.; Denux, D.; Lasnier, J.M.; Ham-Pichavant, F.; Cansell, F.; Aymonier, C. Thermogravimetric analysis as a new method to determine the lignocellulosic composition of biomass. Biomass Bioenergy 2011, 35, 298–307. [Google Scholar] [CrossRef]
  20. Hu, S.; Jess, A.; Xu, M. Kinetic study of Chinese biomass slow pyrolysis: Comparison of different kinetic models. Fuel 2007, 86, 2778–2788. [Google Scholar] [CrossRef]
  21. Aboyade, A.O.; Carrier, M.; Meyer, E.L.; Knoetze, J.H.; Görgens, J.F. Model fitting kinetic analysis and characterisation of the devolatilization of coal blends with corn and sugarcane residues. Thermochim. Acta 2012, 530, 95–106. [Google Scholar] [CrossRef]
  22. 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]
  23. Manya, J.J.; Velo, E.; Puigjaner, L. Kinetics of biomass pyrolysis: A reformulated three-parallel-reactions model. Ind. Eng. Chem. Res. 2003, 42, 434–441. [Google Scholar] [CrossRef]
  24. Burhenne, L.; Messmer, J.; Aicher, T.; Laborie, M.P. The effect of the biomass components lignin, cellulose and hemicellulose on TGA and fixed bed pyrolysis. J. Anal. Appl. Pyrolysis 2013, 101, 177–184. [Google Scholar] [CrossRef]
  25. O’Brien, C.M. Statistical Applications for Environmental Analysis and Risk Assessment by Joseph Ofungwu. Int. Stat. Rev. 2014, 82, 487–488. [Google Scholar] [CrossRef]
  26. Raitanen, J.E.; Järvenpää, E.; Korpinen, R.; Mäkinen, S.; Hellström, J.; Kilpeläinen, P.; Jaana Liimatainen, J.; Ora, A.; Tupasela, T.; Jyske, T. Tannins of Conifer Bark as Nordic Piquancy—Sustainable Preservative and Aroma? Molecules 2020, 25, 567. [Google Scholar] [CrossRef] [Green Version]
  27. Hofbauer, H.; Kaltschmitt, M.; Nussbaumer, T. Energie aus Biomasse–Grundlagen, Techniken, Verfahren; Springer: Berlin, Germany, 2009. [Google Scholar]
  28. Demirbaş, A. Calculation of higher heating values of biomass fuels. Fuel 1997, 76, 431–434. [Google Scholar] [CrossRef]
  29. Rowell, R.M.; Pettersen, R.; Han, J.S.; Rowell, J.S.; Tshabalala, M.A. Cell wall chemistry. In Handbook of Wood Chemistry and Wood Composites; CRC Press: Boca Ratón, FL, USA, 2005; Volume 2. [Google Scholar]
  30. 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]
  31. Pasangulapati, V.; Ramachandriya, K.D.; Kumar, A.; Wilkins, M.R.; Jones, C.L.; Huhnke, R.L. Effects of cellulose, hemicellulose and lignin on thermochemical conversion characteristics of the selected biomass. Bioresour. Technol. 2012, 114, 663–669. [Google Scholar] [CrossRef]
  32. González Martinez, M. Woody and Agricultural Biomass Torrefaction: Experimental Study and Modelling of Solid Conversion and Volatile Species Release Based on Biomass Extracted Macromolecular Components. Ph.D. Thesis, University of Toulouse, Toulouse, France, 2018. [Google Scholar]
  33. Saini, J.K.; Saini, R.; Tewari, L. Lignocellulosic agriculture wastes as biomass feedstocks for second-generation bioethanol production: Concepts and recent developments. 3 Biotech 2015, 5, 337–353. [Google Scholar] [CrossRef] [Green Version]
  34. González Martínez, M.; Dupont, C.; da Silva Perez, D.; Mortha, G.; Thiéry, S.; Meyer, X.M.; Gourdon, C. Understanding the torrefaction of woody and agricultural biomasses through their extracted macromolecular components. Part 1: Experimental thermogravimetric solid mass loss. Energy 2020, 205, 118067. [Google Scholar]
  35. Guo, X.J.; Wang, S.R.; Wang, K.G.; Qian, L.I.U.; Luo, Z.Y. Influence of extractives on mechanism of biomass pyrolysis. J. Fuel Chem. Technol. 2010, 38, 42–46. [Google Scholar] [CrossRef]
  36. Ebringerová, A.; Hromádková, Z.; Heinze, T.; Hemicellulose, T.H. Polysaccharides I. Adv. Polym. Sci. 2005, 186, 67. [Google Scholar]
  37. Yang, H.; Yan, R.; Chen, H.; Lee, D.H.; Zheng, C. Characteristics of hemicellulose, cellulose and lignin pyrolysis. Fuel 2007, 86, 1781–1788. [Google Scholar] [CrossRef]
  38. Chen, W.H.; Wang, C.W.; Ong, H.C.; Show, P.L.; Hsieh, T.H. Torrefaction, pyrolysis and two-stage thermodegradation of hemicellulose, cellulose and lignin. Fuel 2019, 258, 116168. [Google Scholar] [CrossRef]
  39. Yeo, J.Y.; Chin, B.L.F.; Tan, J.K.; Loh, Y.S. Comparative studies on the pyrolysis of cellulose, hemicellulose, and lignin based on combined kinetics. J. Energy Inst. 2019, 92, 27–37. [Google Scholar] [CrossRef]
  40. Korpinen, R.; Kallioinen, M.; Hemming, J.; Pranovich, A.; Mänttäri, M.; Willför, S. Comparative evaluation of various lignin determination methods on hemicellulose-rich fractions of spruce and birch obtained by pressurized hot-water extraction (PHWE) and subsequent ultrafiltration (UF). Holzforschung 2014, 68, 971–979. [Google Scholar] [CrossRef]
  41. Ioelovich, M. Methods for determination of chemical composition of plant biomass. J. SITA 2015, 17, 208–214. [Google Scholar]
Figure 1. Decision tree of the calculation procedure.
Figure 1. Decision tree of the calculation procedure.
Processes 08 01048 g001
Figure 2. DTG curves comparison.
Figure 2. DTG curves comparison.
Processes 08 01048 g002
Figure 3. Model fitted to the experimental pine bark DTG curve.
Figure 3. Model fitted to the experimental pine bark DTG curve.
Processes 08 01048 g003
Figure 4. Model fitted to the experimental spruce bark DTG curve.
Figure 4. Model fitted to the experimental spruce bark DTG curve.
Processes 08 01048 g004
Figure 5. Model fitted to the experimental poplar DTG curve.
Figure 5. Model fitted to the experimental poplar DTG curve.
Processes 08 01048 g005
Figure 6. Model fitted to the experimental willow DTG curve.
Figure 6. Model fitted to the experimental willow DTG curve.
Processes 08 01048 g006
Figure 7. Model fitted to the experimental wheat straw DTG curve.
Figure 7. Model fitted to the experimental wheat straw DTG curve.
Processes 08 01048 g007
Figure 8. Model fitted to the experimental pine bark TG curve.
Figure 8. Model fitted to the experimental pine bark TG curve.
Processes 08 01048 g008
Figure 9. Model fitted to the experimental spruce bark TG curve.
Figure 9. Model fitted to the experimental spruce bark TG curve.
Processes 08 01048 g009
Figure 10. Model fitted to the experimental poplar TG curve.
Figure 10. Model fitted to the experimental poplar TG curve.
Processes 08 01048 g010
Figure 11. Model fitted to the experimental willow TG curve.
Figure 11. Model fitted to the experimental willow TG curve.
Processes 08 01048 g011
Figure 12. Model fitted to the experimental wheat straw TG curve.
Figure 12. Model fitted to the experimental wheat straw TG curve.
Processes 08 01048 g012
Figure 13. Model fitted to the experimental poplar DTG and TG curves: (a) DTG at 3 °C min−1, (b) 5 °C min−1 and (c) 10 °C min−1. (d) TG at the three heating rates.
Figure 13. Model fitted to the experimental poplar DTG and TG curves: (a) DTG at 3 °C min−1, (b) 5 °C min−1 and (c) 10 °C min−1. (d) TG at the three heating rates.
Processes 08 01048 g013aProcesses 08 01048 g013b
Table 1. Literature references of the main lignocellulosic fraction compositions related to used biomasses.
Table 1. Literature references of the main lignocellulosic fraction compositions related to used biomasses.
BiomassRef.Hemicellulose,
wt.%
Cellulose,
wt.%
Lignin,
wt.%
Extractives,
wt.%
Ash,
wt.%
Pine bark[26] a25.019.038.018.0
Spruce bark[24] a27.042.026.0
[27] a24.341.030.0
[28] a21.250.827.5
[26] a28.022.031.019.0
Poplar[17] b26.050.024.0
[19] a28.043.025.05
[29] a18.0–26.646.5–52.016.0–25.9
[3] b22.049.028.0
[30] a24.049.020.05.91.0
Willow[30] a16.741.729.39.72.5
Wheat straw [30] a,c24.639.217.0
[28] a29.038.015.0
[28] a39.128.818.6
[30] a25.037.520.24.03.7
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).
Biomass
Component
Analysis MethodPine BarkSpruce BarkPoplarWillowWheat Straw
HemicelluloseTAPPI T24918.30 13.90 21.70 22.60 23.80
CelluloseTAPPI T24921.90 29.70 42.70 44.30 37.50
LigninTAPPI T22240.70 45.10 26.90 25.10 20.50
ExtractivesInternal Method15.20 4.40 8.00 15.70
TAPPI 204 4.90
AshXP CEN/TS 147752.80 2.80 2.30 8.30
TAPPI 211 5.22
Table 3. Minimum number of components for each biomass fraction.
Table 3. Minimum number of components for each biomass fraction.
ComponentTemperature Range, °CNumber of Pseudocomponents
Water25–1501
Hemicellulose200–3502
Cellulose250–4001
Lignin150–10003
Table 4. Upper bonds of the pseudocomponents (PC).
Table 4. Upper bonds of the pseudocomponents (PC).
Kinetic
Parameters
PC
2
PC
3
PC
4
PC
5
PC
6
PC
7
K (s−1)1.00 × 1091.50 × 1052.40 × 10155.00 × 1013.001.80
E (kJ mol−1)120.0080.00240.0060.0060.0068.00
Xj,0 (wt.%)50.0050.0060.0060.0020.00-
Table 5. Lower bonds of the pseudocomponents (PC).
Table 5. Lower bonds of the pseudocomponents (PC).
Kinetic
Parameters
PC
2
PC
3
PC
4
PC
5
PC
6
PC
7
K (s−1)7.00 × 1081.40 × 1051.50 × 10154.00 × 1072.301.00 × 10−1
E (kJ mol−1)100.0070.00160.0055.0045.0050.00
Xj,0 (wt.%)0.115.0015.000.10-
Table 6. Kinetic parameters of the pseudocomponents.
Table 6. Kinetic parameters of the pseudocomponents.
WaterHemicelluloseCelluloseLignin
BiomassKinetic ParametersPC
1
PC
2
PC
3
PC
4
PC
5
PC
6
PC
7
Pine barkK (s−1)9.14 × 1046.00 × 1081.50 × 1051.69 × 10155.00 × 102.449.83 × 10−1
E (kJ mol−1)48.58120.0075.60204.7455.0050.2061.31
Xj,0 (wt.%)6.8114.717.0024.6627.6813.165.98
Spruce barkK (s−1)4.52 × 1037.00 × 1081.42 × 1051.51 × 10155.00 × 102.331.11
E (kJ mol−1)40.74119.9974.93205.1755.0048.7960.95
Xj,0 (wt.%)8.2514.046.5624.5324.6114.187.82
PoplarK (s−1)7.22 × 1056.00 × 1081.40 × 1052.30 × 10154.81 × 102.795.07 × 10−1
E (kJ mol−1)52.44120.0080.00207.3955.0053.2259.71
Xj,0 (wt.%)3.8121.721.0051.8515.344.242.04
WillowK (s−1)2.98 × 1056.00 × 1081.47 × 1051.94 × 10155.00 × 102.608.77 × 10−1
E (kJ mol−1)50.69120.0074.88208.0355.0049.4362.61
Xj,0 (wt.%)4.6921.911.9144.3315.008.253.91
Wheat strawK (s−1)8.23 × 1056.00 × 1081.50 × 1051.51 × 10155.00 × 103.001.62 × 10−1
E (kJ mol−1)53.65120.0071.38200.6455.0051.0852.16
Xj,0 (wt.%)5.1824.161.0039.5119.736.903.51
Table 7. Kinetic parameters from other studies.
Table 7. Kinetic parameters from other studies.
ComponentTemperature, °CE, kJ mol−1K, min−1Reference
Hemicellulose200–350127.009.5 × 1010[38]
83.20–96.404.55 × 106–1.57 × 108[39]
Cellulose300–340227.023.36 × 1018[37]
239.70–325.0016.30 × 1019–3.62 × 1026[39]
Lignin220–3807.802.96 × 10−3[37]
25–90047.90–54.506.80 × 102–6.60 × 104[17]
160–68025.204.70 × 102[18]
20.00–29.105.35 × 10–3.18[39]
Table 8. Quality of the fit expressed as R2Adj and QOF%.
Table 8. Quality of the fit expressed as R2Adj and QOF%.
BiomassQOF%R2Adj
Pine bark1.510.9939
Spruce bark1.940.9905
Poplar1.090.9960
Willow1.420.9933
Wheat straw1.790.9921
Table 9. Comparison between the analytical and thermogravimetric analysis-pseudocomponent kinetic model (TGA-PKM) results.
Table 9. Comparison between the analytical and thermogravimetric analysis-pseudocomponent kinetic model (TGA-PKM) results.
BiomassComponent Analytical Method
wt.%,
Dry, Ash and Extractives-Free Basis
TGA-PKM Method
wt.%,
Dry, Ash and Extractives-Free Basis
Error,
wt.%
PoplarHemicellulose23.77 23.62−0.15
Cellulose46.77 53.907.14
Lignin29.46 22.48−6.99
WillowHemicellulose24.57 25.000.43
Cellulose48.15 46.51−1.64
Lignin27.28 28.491.21
Wheat strawHemicellulose29.10 26.54−2.56
Cellulose45.84 41.67−4.18
Lignin25.06 31.806.73
Spruce BarkHemicellulose15.67 22.456.78
Cellulose33.48 26.74−6.74
Lignin50.85 50.81−0.04
Pine barkHemicellulose22.62 23.300.68
Cellulose27.07 26.47−0.61
Lignin50.31 50.24−0.07
Table 10. Quality of the fit for each dataset.
Table 10. Quality of the fit for each dataset.
Quality of the Fit 3 °C min−15 °C min−110 °C min−1Global
QOF%1.161.580.961.35
R2Adj0.99570.99170.99710.9959
Table 11. Kinetic parameters of each pseudocomponents calculated using a single heating rate and three simultaneous heating rates.
Table 11. Kinetic parameters of each pseudocomponents calculated using a single heating rate and three simultaneous heating rates.
HemicelluloseCelluloseLignin
Number of
Heating Rates
Kinetic ParametersPC
2
PC
3
PC
4
PC
5
PC
6
PC
7
Single heating rateK (s−1)6.00 × 1081.40 × 1052.30 × 10154.81 × 1012.79 × 1005.07 × 10−1
E (kJ mol−1)120.0080.00207.3955.0053.2259.71
Xj,0 (wt.%)21.721.0051.8515.344.242.04
Three simultaneous heating ratesK (s−1)6.56 × 1081.50 × 1052.23 × 10155.49 × 1012.79 × 1005.57 × 10−1
E (kJ mol−1)119.9877.18207.4855.0053.9957.48
Xj,0 (wt.%)22.781.4649.9715.083.881.51
Table 12. Comparison between TGA-PKM results using three simultaneous heating rates and the analytical method.
Table 12. Comparison between TGA-PKM results using three simultaneous heating rates and the analytical method.
BiomassComponent Analytical Method
Wt.%,
Dry, Ash and Extractives-Free Basis
TGA-PKM Method
(Three Simultaneous Heating Rates)
wt.%,
Dry, Ash and Extractives-Free Basis
Error,
wt.%
PoplarHemicellulose23.77 25.601.83
Cellulose46.77 52.786.01
Lignin29.46 21.62−7.85

Share and Cite

MDPI and ACS Style

Díez, D.; Urueña, A.; Piñero, R.; Barrio, A.; Tamminen, T. Determination of Hemicellulose, Cellulose, and Lignin Content in Different Types of Biomasses by Thermogravimetric Analysis and Pseudocomponent Kinetic Model (TGA-PKM Method). Processes 2020, 8, 1048. https://doi.org/10.3390/pr8091048

AMA Style

Díez D, Urueña A, Piñero R, Barrio A, Tamminen T. Determination of Hemicellulose, Cellulose, and Lignin Content in Different Types of Biomasses by Thermogravimetric Analysis and Pseudocomponent Kinetic Model (TGA-PKM Method). Processes. 2020; 8(9):1048. https://doi.org/10.3390/pr8091048

Chicago/Turabian Style

Díez, David, Ana Urueña, Raúl Piñero, Aitor Barrio, and Tarja Tamminen. 2020. "Determination of Hemicellulose, Cellulose, and Lignin Content in Different Types of Biomasses by Thermogravimetric Analysis and Pseudocomponent Kinetic Model (TGA-PKM Method)" Processes 8, no. 9: 1048. https://doi.org/10.3390/pr8091048

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop