5.2
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Letter to Editor
Original Article
Research article
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
5.3
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
Corrigendum
Current Issue
Editorial
Erratum
Full Length Article
Full lenth article
Letter to Editor
Original Article
Research article
Retraction notice
Review
Review Article
SPECIAL ISSUE: ENVIRONMENTAL CHEMISTRY
View/Download PDF

Translate this page into:

Original Article
ARTICLE IN PRESS
doi:
10.25259/AJC_1267_2025

Thermogravimetric analysis and pyrolysis kinetics of typical crop straws

Institute of Molecular Engineering and Applied Chemistry, Anhui University of Technology, Hudong Road, Ma’anshan, Anhui, P. R. China

* Corresponding author: E-mail address: zhangqf@ahut.edu.cn (Q.F. Zhang)

Licence
This is an open-access article distributed under the terms of the Creative Commons Attribution-Non Commercial-Share Alike 4.0 License, which allows others to remix, transform, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms.

Abstract

Using corn stalks, rapeseed stalks, and cotton stalks as raw materials, the pyrolysis characteristics and kinetic parameters in a nitrogen atmosphere were analyzed from a kinetic perspective through thermogravimetric (TG) and differential thermogravimetric (DTG) curves at four different heating rates of 10, 20, 30, and 40°C/min. Three model-free methods—Flynn-Wall-Ozawa (FWO), Kissinger-Akahira-Sunose (KAS), and Starink —were employed to calculate the apparent activation energy and pre-exponential factor. Combining the Coats-Redfern (CR) model method with 17 mechanistic models, the optimal reaction mechanism model was ultimately obtained. Research findings indicate that the pyrolytic weight loss process of straw can be divided into three stages: drying and dehydration, pyrolysis of organic compounds such as cellulose/hemicellulose/lignin, and carbonization. The apparent activation energies calculated for corn straw using the FWO method, KAS method, and Starink method were 85.16 kJ mol-1, 75 kJ mol-1, and 75.41 kJ mol-1, respectively; for rapeseed straw: 99.7 kJ mol-1, 89.71 kJ mol-1, and 90 kJ mol-1; for cotton straw: 90.17 kJ mol-1, 79.62 kJ mol-1, and 80.04 kJ mol-1. The CR method revealed that the mechanism equations for pyrolysis of all three straw types conform to the F3 model. Additionally, under N₂ atmosphere, the thermodynamic properties of the reaction—specifically ΔH, ΔG, and ΔS—were significantly higher for rapeseed straw than for corn and cotton straw, confirming the greater complexity of its reaction. These findings contribute to a deeper understanding of straw pyrolysis mechanisms and provide a theoretical foundation for the resource utilization of solid waste.

Keywords

Activation energy
Kinetics
Pre-exponential factor
Reaction mechanism model
Straw
Thermogravimetric analysis

1. Introduction

Driven by global carbon neutrality strategies, biomass energy has emerged as a key renewable resource to replace fossil fuels due to its zero-carbon cycle characteristics [1,2]. As an agricultural powerhouse, China produced over 290 million tons of corn stover in 2024, according to National Bureau of Statistics data. Open-air burning of this crop residue results in PM2.5 emissions reaching millions of tons annually. With its high cellulose content and low ash characteristics, corn stover is considered an ideal feedstock for pyrolysis-derived bio-oil and syngas. Pyrolysis technology enables the conversion of corn stover into high-value bio-oil and syngas, forming a core link in the agricultural waste-to-biofuel pathway [3-5]. Ultimately, efficient energy conversion achieves dual objectives of pollution control and clean energy supply, offering broad application prospects.

Currently, thermogravimetric analysis serves as the core experimental method for deciphering biomass pyrolysis kinetics. Wang et al. employed thermogravimetric-differential thermogravimetric-differential scanning calorimetry (TG-DTG-DSC) to calculate the combustible properties, burnout characteristics, and comprehensive combustion properties of corn stover and corn cobs, thereby analyzing their combustion behavior [6]. Kai et al. investigated the desorption kinetics and evolution of gaseous species during the co-pyrolysis of corn stover and high-density polyethylene using thermogravimetric-fourier transform infrared spectrometry-mass spectrometry (TG-FTIR-MS). Additionally, they explored potential interactions between corn stover and high-density polyethylene during co-pyrolysis [7]. Ding et al. compared differences between the pyrolysis processes of broadleaf wood and coniferous wood via TG analysis. Results indicate that coniferous wood exhibits a higher overall energy barrier during pyrolysis compared to broadleaf wood [8].

Over the past few years, research on biomass pyrolysis has primarily focused on examining how various processing conditions and material types influence product yields and characteristics. These variables include different pretreatment methods, temperature [9-12], heating rate [13], and reaction atmosphere [14]. Pyrolysis temperature significantly influences product distribution and performance. As heating rates increase, corn stover char yield decreases while gas selectivity improves [15]. This arises because heating rates alter heat/mass transfer processes, introducing method-dependent biases in apparent activation energy calculations. Hu et al. conducted a systematic study on the pyrolysis characteristics of corn stover using an innovative infrared-heated pyrolysis reactor. Specifically, they analyzed thermal degradation behavior at heating rates of 10, 20, 30, and 40°C min-1 within a thermogravimetric analyzer, calculating pyrolysis kinetic and thermodynamic parameters [16]. Biomass pyrolysis is typically conducted under an inert atmosphere, though other gases may be introduced to modify the pyrolysis process.

Kinetics analysis techniques primarily encompass model-free fitting methods [17-19] and model-based fitting methods [5,20]. At the methodological level of kinetic analysis, existing research faces threefold limitations: (1) Model-free approaches, while avoiding mechanistic assumptions, cannot directly obtain the reaction mechanism function g(α); (2) Model-based approaches, while relying on predefined g(α), lack objective criteria for selecting solid-phase reaction mechanism models; (3) Parameter fragmentation optimization, where the pre-exponential factor is often obtained through extrapolation from Arrhenius plots, neglecting the kinetic compensation effect between E and A, leading to inaccurate predictions. Therefore, it is essential to thoroughly investigate the pyrolysis kinetics, thermodynamics, and characteristics of corn stover.

To address these challenges, this study proposes a research framework combining model-free methods for parameter estimation and model-based methods for model development: TG/DTG curves for corn stover, rapeseed straw, and cotton stalks were obtained under four heating rates (10, 20, 30, 40°C/min) in a nitrogen atmosphere, and pyrolysis stages were delineated. Three model-free methods—Flynn-Wall-Ozawa (FWO), Kissinger-Akahira-Sunose (KAS), and Starink—were employed to calculate activation energy. Concurrently, the pre-exponential factor was back-calculated using the compensation effect equation, overcoming the limitations of traditional single-point extrapolation. Subsequently, the optimal g(α) was determined using the Coats-Redfern model combined with 17 mechanism models, based on the correlation coefficient R2. Due to variations in computational methods and pyrolysis temperature ranges chosen by different researchers, activation energies for the same material exhibit discrepancies. These findings contribute to understanding the pyrolysis process of straw and hold significant implications for practical applications in straw pyrolysis technology.

The core uniqueness of this study lies in achieving “multidimensional systematic integration and deep structure-property correlation”: Firstly, it encompasses three typical straw feedstocks—corn, rapeseed, and cotton—combined with four heating rate gradients to establish a more universal pyrolysis kinetics research framework; Secondly, it cross-validates kinetic parameter reliability by integrating three non-model methods (FWO, KAS, Starink) with 17 Coats-Redfern mechanistic models, precisely identifying optimal reaction mechanisms; Thirdly, the direct correlation between the changes of the microstructure (pore structure, surface morphology) and the changes of the thermogravimetric data of biomass before and after pyrolysis was established by using the scanning electron microscopy (SEM) characterization technique. Compared to previous studies, this research not only quantifies the kinetic and thermodynamic differences among various straw types but also reveals the fundamental structural causes of pyrolysis variations. It establishes a complete chain of evidence linking “macroscopic thermal behavior–microstructure–reaction mechanism,” providing more precise theoretical support for optimizing processes in the resource utilization of straw-based solid waste.

2. Materials and Methods

2.1. Biomass sample

Corn stalk samples collected in Zhengzhou, Henan Province, China; rapeseed stalk samples collected in Zhengzhou, Henan Province, China; and cotton stalk samples collected in Lianyungang, Jiangsu Province, China were all pre-processed using a shredder. After grinding, the material was passed through a 50-mesh standard sieve (<0.25 mm) to obtain the initial sample.

2.2. Instrument and experimental conditions

Thermogravimetric analysis was conducted on a Shimadzu DTG-60H analyzer in Japan. Approximately 3–5 mg of straw samples were heated from room temperature to 800°C at heating rates of 10, 20, 30, and 40°C min-1. High-purity N₂ (99.999%) was used as the purge gas at a flow rate of 30 mL/min to maintain an inert atmosphere. To ensure the authenticity of the curves, each experimental group was designed with at least three replicates.

According to GB/T 28731-2012 Industrial Analysis Methods for Solid Biomass Fuels, the moisture, volatile matter, ash, and fixed carbon content in corn stalks, rapeseed stalks, and cotton stalks were estimated. The fixed carbon content was calculated using the mathematical formula: [fixed carbon% = 100% - (moisture content + volatile matter + ash content)%]. Elemental analysis was performed using an organic elemental analyzer (VARIO EL III Organic Elemental Analyzer) in CHNS mode and O mode.

2.3. Kinetic model

2.3.1. FWO model

The apparent activation energy (E, kJ mol-1) of the straw pyrolysis reaction was calculated using the FWO method, whose expression is shown in equation (1).

(1)
ln β = ln A E R g ( α ) 1.052 E R T 5.331

In equation (1): β is the heating rate, °C/min; A is the pre-exponential factor, min⁻1; g(α) is the integral form of the mechanism function; E is the activation energy, kJ mol-1; R is the ideal gas constant, 8.314 J/mol·K); T is the absolute temperature, K. α is the conversion rate of the reaction; the conversion rate (α) at a given moment during biomass pyrolysis is defined as: α is the conversion rate of the reaction; the conversion rate (α) at a given moment during biomass pyrolysis is defined as:

(2)
α = m 0 - m t m 0 - m f ( 0 < α < 1 )

In equation (2): m₀ is the initial mass of the sample; mf is the final mass of the sample; mt is the instantaneous mass of the sample at time t.

Therefore, when α is a constant value, plot a scatter plot of lnβ against 1/T as the independent variable for different heating rates. The slope of the fitted curve can be used to determine E.

2.3.2. KAS model

The apparent activation energy for the pyrolysis reaction of straw was calculated using the KAS method, whose expression is shown in equation (3).

(3)
ln β T 2 = ln A R E g ( α ) E R T

Therefore, when α is a constant value, plot a scatter diagram of ln β T 2 versus 1/T as the independent variable for different heating rates. The slope of the fitted curve can then be used to determine E.

2.3.3. Starink model

The apparent activation energy for the pyrolysis reaction of straw was calculated using the Starink method, whose expression is shown in equation (4).

(4)
ln β T 1.92 = c 1.0008 E R T

Therefore, when α is a constant value, plot a scatter diagram of ln β T 1.92 against 1/T as the independent variable for different heating rates. By fitting the curve, the slope can be used to determine E.

2.3.4. Pre-exponential factor (A) forecast

After obtaining the activation energy E in the model-free approach, the Kissinger method can predict the pre-exponential factor by analyzing the displacement of pyrolysis peak temperatures across multiple heating rates. The expression is shown in equation (5):

(5)
A = β E exp E R T m R T m 2

In equation (5): Tm (K) is the temperature at the maximum value of DTG.

2.3.5. Thermodynamic analysis

Thermodynamic parameters can be predicted from kinetic data. Thermodynamic parameters (ΔH, ΔG, and ΔS) are calculated using the following formula:

(6)
Δ H = E R T

(7)
Δ G = E + R T m ln K B T m h A

(8)
Δ S = Δ H Δ G T m

In equations (6)−(8): ΔH denotes enthalpy change, kJ mol-1; ΔG denotes Gibbs free energy change, kJ mol-1; ΔS denotes entropy change, kJ/K; R denotes the ideal gas constant, 8.314×10⁻3 kJ/(mol·K); T is the characteristic reaction temperature, typically the average or peak reaction temperature, K; Tm is the maximum peak temperature of the DTG curve, K; KB is the Boltzmann constant, 1.381×10⁻23 J/K; h is Planck’s constant, 6.626×10⁻3⁴ J/s.

2.3.6. Coats-Redfern (C-R) model

The CR method was employed to calculate the apparent activation energy and pre-exponential factor for the pyrolysis reaction of straw. The expression for the CR method is shown in equation (9). The g(α) functions for 17 common pyrolysis kinetic mechanisms are listed in Table 1.

(9)
ln g α T 2 = ln A R β E E R T

Table 1. Common functions for pyrolysis kinetic mechanisms.
Mechanism of reaction Symbol f(α) g(α)
First-order reaction F1 1-α -ln (1-α)
Second-order reaction F2 (1-α)2 (1-α)-1-1
Third-order reaction F3 (1-α)3 [(1-α)-2-1]/2
One-dimensional diffusion D1 0.5α-1 α2
Two-dimensional diffusion D2 [-ln (1-α)]-1 α+(1-α) ln(1-α)
Three-dimensional diffusion D3 [3/2(1-α)2/3]/[1-(1-α)1/3] [1-(1-α)1/3]2
Nucleation and growth A2 2(1-α) [-ln(1-α)]1/2 [-ln(1-α)]1/2
Nucleation and growth A3 3(1-α) [-ln(1-α)]2/3 [-ln(1-α)]1/3
Nucleation and growth A4 4(1-α) [-ln(1-α)]3/4 [-ln(1-α)]1/4
Phase boundary reaction R1 1 α
Phase boundary reaction R2 2(1-α)1/2 1-(1-α)1/2
Phase boundary reaction R3 3(1-α)2/3 1-(1-α)1/3
Power law P2/3 2/3α-1/2 α3/2
Power law P2 1/2 α1/2
Power law P3 2/3 α1/3
Power law P4 3/4 α1/4

Therefore, when α is a constant value, plot scatter plots of ln g ( α ) T 2 versus 1/T as the independent variable for different heating rates. The slope of the fitted curve can be used to determine E.

3. Results and Discussion

3.1. Industrial analysis and elemental analysis

The industrial analysis results for straw are shown in Table 2. As indicated in Table 2, the moisture, ash, volatile matter, and fixed carbon content of corn straw are 7.72%, 7.71%, 69.38%, and 15.19%, respectively. Rapeseed straw exhibits moisture, ash, volatile matter, and fixed carbon content of 9.08%, 9.75%, 67.02%, and 14.15%, respectively. Cotton straw contains 6.52%, 15.67%, 63.56%, and 14.25% moisture, ash, volatile matter, and fixed carbon, respectively. Significant differences exist among these three typical crop straws. Corn and rapeseed straw exhibit low ash content, classifying them as low-ash fuels. Cotton straw contains markedly higher ash content than corn and rapeseed straw, attributed to its higher silica content in the outer husk. All three crop residues exhibit high volatile matter content, indicating they ignite more readily and burn more completely than coal. This demonstrates the significant potential of crop residues as a biofuel source [21]. The higher fixed carbon content in corn stalks compared to rapeseed and cotton stalks stems from corn stalks’ greater lignification, resulting in relatively higher fixed carbon levels [22].

Table 2. Industrial analysis of straw (wt.%, ad).
Sample Moisture Ash content Volatile matter Fixed carbon
Corn straw 7.72 7.71 69.38 15.19
Rapeseed straw 9.08 9.75 67.02 14.15
Cotton straw 6.52 15.67 63.56 14.25

The elemental analysis results for straw are shown in Table 3. As indicated in Table 3, corn straw contains 49.72% C, 7.16% H, 41.47% O, 1.49% N, and 0.16% S. Rapeseed straw contains 44.24% C, 6.20% H, 48.26% O, 0.92% N, and 0.38% S. Cotton stalks contain 44.13% C, 6.05% H, 48.77% O, 0.70% N, and 0.35% S. All three crop stalks exhibit relatively high C and O content. The H/C molar ratio reflects the content of light hydrocarbons in the straw. A higher H/C ratio indicates a higher proportion of volatile components. The low sulfur content suggests that crop straw contributes minimally to atmospheric SO2 emissions during energy utilization. Compared to coal, straw qualifies as a clean energy source.

Table 3. Elemental analysis of straw.
Sample C H O N S H/C
Corn straw 49.72 7.16 41.47 1.49 0.16 0.144
Rapeseed straw 44.24 6.20 48.26 0.92 0.38 0.140
Cotton straw 44.13 6.05 48.77 0.70 0.35 0.137

3.2. TG analysis

The thermal decomposition behavior of corn stover, rapeseed straw, and cotton straw was investigated using thermogravimetric analysis (TGA) over a temperature range from room temperature to 800°C. Figures 1-3 illustrate the decomposition behavior of corn stover, rapeseed straw, and cotton stalks, respectively, as a function of temperature. The TG and DTG curves reveal that: At all heating rates, the TG and DTG curves exhibit similar trends. As the heating rate increases, the TG and DTG curves shift toward higher temperatures [23,24].

Thermogravimetric curves of corn straw at different heating rates: (a) TG curve; (b) DTG curve.
Figure 1.
Thermogravimetric curves of corn straw at different heating rates: (a) TG curve; (b) DTG curve.
Thermogravimetric curves of rapeseed straw at different heating rates: (a) TG curve; (b) DTG curve.
Figure 2.
Thermogravimetric curves of rapeseed straw at different heating rates: (a) TG curve; (b) DTG curve.
Thermogravimetric curves of cotton straw at different heating rates: (a) TG curve; (b) DTG curve.
Figure 3.
Thermogravimetric curves of cotton straw at different heating rates: (a) TG curve; (b) DTG curve.

Simultaneously, the pyrolytic weight loss process of straw can be divided into three stages: drying and dehydration, pyrolysis of organic compounds such as cellulose/hemicellulose/lignin, and carbonization [25]. During the first stage, the feedstock does not decompose with increasing temperature but undergoes drying and dehydration. This stage spans a temperature range of approximately room temperature to 149.92°C (corn stover), room temperature to 165.89°C (rapeseed straw), and room temperature to 162.8°C (cotton straw). The weight loss percentages for corn stalks, rapeseed stalks, and cotton stalks were approximately 10.9%, 12.4%, and 9.37%, respectively. The rate of weight loss gradually decreased toward the end of this stage. According to related research reports, other crop straws exhibit similar trends. In the second stage, cellulose, hemicellulose, and lignin decompose. The pyrolysis rate of the raw material accelerates from the start to a maximum value before slowing down. The TG curve declines rapidly, and the DTG curve exhibits a peak. During this process, mass decreases rapidly with a relatively fast pyrolysis rate. The primary pyrolysis temperature range for this stage is approximately 149.92°C to 708.68°C (corn stover), 165.89°C to 646.20°C (rapeseed straw), and 162.8°C to 765.16°C (cotton stalk). During the third stage, residual material undergoes slow decomposition and carbonization until mass changes cease. The temperature range for this stage is approximately 708.68°C to 800°C (corn stover), 646.20°C to 800°C (rapeseed straw), and 765.16°C to 800°C (cotton stalk). The residual components after combustion are primarily ash. Industrial analysis indicates that cotton stalks exhibit the highest ash content at 15.67%, resulting in the greatest residue after thermogravimetric analysis. Corn stalks and rapeseed stalks show similar ash content levels, implying comparable thermogravimetric residues.

Taking corn stover as an example, we also observe that when the heating rate increases from 10 to 40°C min-1, the temperature of maximum weight loss rate rises from 300.24°C to 352.70°C, while the maximum mass loss rate increases from 39.057% to 44.177%. The results indicate that the reaction rate is closely related to the heating rate. This phenomenon, known as thermal lag, has been reported in numerous studies on biomass pyrolysis [26]. As the heating rate increased from 10 to 40°C/min, the residual mass range expanded from 4.439% to 10.199%. The highest residual mass was observed at a heating rate of 40°C min-1, while the lowest residual mass occurred at 10°C min-1. This aligns with the volatiles observed in Table 2. The variation in residual mass primarily stems from slower heating rates favoring solid product formation. Results indicate that heating rate has no significant effect on residual mass. Similar conclusions apply to rapeseed straw and cotton stalks.

3.3. Model-free method

Evaluating kinetic factors of biomass, such as activation energy, aids in understanding the reactivity of specific feedstocks. Figures 4-6 show the T-α curves for the pyrolysis of corn stover, rapeseed straw, and cotton stalks, respectively. We observe consistent patterns among the three straw types concerning heating rate, pyrolysis temperature, and α. At a constant heating rate, α increases with rising pyrolysis temperature; at a constant pyrolysis temperature, α decreases with increasing heating rate.

T-α curve of corn stalk pyrolysis.
Figure 4.
T-α curve of corn stalk pyrolysis.
T-α Curve of rapeseed stalk pyrolysis.
Figure 5.
T-α Curve of rapeseed stalk pyrolysis.
T-α curve of cotton stalk pyrolysis.
Figure 6.
T-α curve of cotton stalk pyrolysis.

The apparent activation energies of the three types of straw measured in this study are consistent with previous studies and reflect the differences in raw material characteristics. The FWO, KAS, and Starink methods were employed to predict the thermodynamics of corn stover, rapeseed straw, and cotton straw. Figure 7 shows the linear fitting curve for corn stover at conversion rates ranging from 0.1 to 0.9. Figure 8 shows the linear fitting plot for rapeseed straw at conversion rates ranging from 0.1 to 0.9. Figure 9 presents the linear fitting plot for cotton straw at conversion rates ranging from 0.1 to 0.9. The fitting plots indicate that E changes as the straw pyrolysis conversion rate α increases from 0.1 to 0.9. For corn stover, the FWO method yielded an E range of 56.15–102.31 kJ mol-1 with an average of 85.16 kJ mol-1; The KAS method yielded E values ranging from 43.55 to 93.08 kJ mol-1, with an average of 75 kJ mol-1; The Starink method yielded E values ranging from 44.06 to 93.45 kJ mol-1, with an average of 75.41 kJ mol-1. This aligns with the literature-reported activation energy of corn stalks. The slightly higher value stems from our study’s cross-validation using four heating rates, yielding more precise data [27,28]. For rapeseed straw, the FWO method yielded an E range of 69.45–110.59 kJ mol-1 with an average of 99.7 kJ mol-1; the KAS method yielded an E range of 57.09–102.07 kJ mol-1 with an average of 89.71 kJ mol-1; Calculations using the Starink method yielded an E range of 57.59–102.39 kJ mol-1, with an average value of 90 kJ mol-1. The activation energy of rapeseed straw is higher than that reported in the literature, which is related to the higher lignin content in this study. The stable structure of lignin leads to an increased demand for bond-breaking energy [29-31]. For cotton straw, the FWO method yielded an E range of 47.17–110.77 kJ mol-1, with an average value of 90.17 kJ mol-1; Calculations using the KAS method yielded activation energy values ranging from 33.87 to 100.80 kJ mol-1, with an average of 79.62 kJ mol-1; calculations using the Starink method yielded activation energy values ranging from 34.40 to 101.196 kJ mol-1, with an average of 80.04 kJ mol-1. The activation energy of cotton stalks is highly consistent with the literature, which confirms the reliability of the combined application of the model method in this study [32]. Compared to two other typical crops, rapeseed straw exhibited slightly higher activation energy estimates, while corn stover showed slightly lower estimates. The study revealed that corn stover requires less energy for thermal decomposition than rapeseed straw and cotton stover.

Linear fitting curves of corn stover at different conversion rates ranging from 0.1 to 0.9 (a) FWO method; (b) KAS method; (c) Starink method.
Figure 7.
Linear fitting curves of corn stover at different conversion rates ranging from 0.1 to 0.9 (a) FWO method; (b) KAS method; (c) Starink method.
Linear fitting curves of rapeseed stover at different conversion rates ranging from 0.1 to 0.9 (a) FWO method; (b) KAS method; (c) Starink method.
Figure 8.
Linear fitting curves of rapeseed stover at different conversion rates ranging from 0.1 to 0.9 (a) FWO method; (b) KAS method; (c) Starink method.
Linear fitting curves of cotton stover at different conversion rates ranging from 0.1 to 0.9 (a) FWO method; (b) KAS method; (c) Starink method.
Figure 9.
Linear fitting curves of cotton stover at different conversion rates ranging from 0.1 to 0.9 (a) FWO method; (b) KAS method; (c) Starink method.

At a heating rate of 10°C min-1, the values of E, A, and R2 calculated from the FWO method, KAS method, and Starink method at different conversion rates are shown in Tables 4-6. The activation energy was calculated from the slope and intercept of the fitted lines. The pre-exponential factor was predicted based on the activation energy, and R2 was used to assess the accuracy of the kinetic parameters. As shown in Table 4, the R2 values for all corn stover curves fell within the narrow range of 0.9593–0.9999, indicating high precision and good data point fitting. At all heating rates, all R2 values exceeded 0.95, confirming the reliability of the pyrolysis kinetic parameters for straw. Furthermore, the activation energy and frequency factor were not constant but varied with increasing conversion rate, indicating the presence of multiple single-reaction mechanisms in straw pyrolysis. As shown in Table 5, the R2 values for all rapeseed straw curves fall within the narrow range of 0.9831–0.9999. As shown in Table 6, the R2 values for all cotton straw curves fall within the narrow range of 0.9695–0.9999.

Table 4. Kinetic parameters of corn stover under different methods.
Conversion rate Method FWO
KAS
Starink
α E/(kJ/mol) A/s-1 R2 E/(kJ/mol) A/s-1 R2 E/(kJ/mol) A/s-1 R2
0.10 91.20 1.13E+06 0.9989 82.55 1.67E+05 0.9987 82.90 1.80E+05 0.9987
0.20 100.79 9.34E+06 0.9993 91.82 1.30E+06 0.9991 92.18 1.40E+06 0.9991
0.30 102.31 1.31E+07 0.9997 93.08 1.71E+06 0.9997 93.45 1.86E+06 0.9997
0.40 99.54 7.11E+06 0.9999 90.07 8.82E+05 0.9998 90.45 9.59E+05 0.9998
0.50 98.73 5.95E+06 0.9998 89.01 6.98E+05 0.9998 89.40 7.61E+05 0.9998
0.60 89.18 7.23E+05 0.9991 79.09 7.73E+04 0.9989 79.49 8.46E+04 0.9989
0.70 67.09 5.30E+03 0.9991 56.18 4.50E+02 0.9985 56.62 4.97E+02 0.9985
0.80 61.46 1.49E+03 0.9928 49.68 1.02E+02 0.9880 50.15 1.13E+02 0.9882
0.90 56.15 4.47E+02 0.9767 43.55 2.46E+01 0.9593 44.06 2.77E+01 0.9603
Average 85.16 4.15E+06 75.00 5.37E+05 75.41 5.83E+05
results calculated at 10 °C/min.
Table 5. Kinetic parameters of rapeseed stover under different methods.
Conversion rate Method FWO
KAS
Starink
α E/(kJ/mol) A/s-1 R2 E/(kJ/mol) A/s-1 R2 E/(kJ/mol) A/s-1 R2
0.10 110.59 1.28E+08 0.9857 102.07 1.91E+07 0.9831 102.39 2.05E+07 0.9832
0.20 107.81 6.88E+07 0.9918 99.01 9.63E+06 0.9901 99.35 1.04E+07 0.9902
0.30 103.18 2.45E+07 0.9973 94.08 3.19E+06 0.9966 94.48 3.49E+06 0.9967
0.40 102.80 2.25E+07 0.9994 93.58 2.85E+06 0.9993 93.91 3.07E+06 0.9993
0.50 107.17 5.96E+07 0.9991 97.73 7.23E+06 0.9989 98.08 7.82E+06 0.9989
0.60 103.99 2.93E+07 0.9914 94.15 3.24E+06 0.9897 93.68 2.92E+06 0.9890
0.70 101.11 1.54E+07 0.9937 90.13 1.32E+06 0.9923 90.57 1.45E+06 0.9924
0.80 91.24 1.69E+06 0.9932 79.51 1.20E+05 0.9907 79.98 1.33E+05 0.9908
0.90 69.45 1.22E+04 0.9901 57.09 7.16E+02 0.9846 57.59 8.04E+02 0.9848
Average 99.70 3.89E+07 89.71 5.19E+06 90.00 5.53E+06
results calculated at 10 °C/min.
Table 6. Kinetic parameters of cotton stover under different methods.
Conversion rate Method FWO
KAS
Starink
α E/(kJ/mol) A/s-1 R2 E/(kJ/mol) A/s-1 R2 E/(kJ/mol) A/s-1 R2
0.10 99.93 7.18E+06 0.9884 91.029 1.02E+06 0.9859 91.38 1.10E+06 0.9860
0.20 105.06 2.21E+07 0.9934 95.855 2.94E+06 0.9920 96.22 3.18E+06 0.9921
0.30 107.47 3.74E+07 0.9960 98.023 4.73E+06 0.9951 98.40 5.13E+06 0.9951
0.40 110.51 7.25E+07 0.9975 100.805 8.70E+06 0.9969 101.19 9.46E+06 0.9969
0.50 110.77 7.68E+07 0.9993 100.798 8.68E+06 0.9992 101.20 9.48E+06 0.9992
0.60 104.51 1.96E+07 0.9928 93.987 1.95E+06 0.9910 94.41 2.14E+06 0.9911
0.70 70.64 1.11E+04 0.9946 59.140 8.40E+02 0.9917 59.60 9.32E+02 0.9918
0.80 55.48 3.67E+02 0.9891 43.063 2.12E+01 0.9804 43.56 2.38E+01 0.9809
0.90 47.17 5.49E+01 0.9858 33.865 2.44E+00 0.9695 34.40 2.77E+00 0.9706
Average 90.17 2.62E+07 79.62 3.11E+06 80.04 3.39E+06
results calculated at 10 °C/min.

Low pre-exponential factor estimates indicate that thermal degradation processes exhibit lower reactivity compared to surface reactions. High pre-exponential factor estimates suggest high reactivity associated with complex polyhedral reaction processes. The estimated value of A increases with decomposition progression, attributed to heightened collision frequency between molecules at elevated conversion levels and heating rates.

Figure 10 shows the relationship between E and A of corn stover and α at different conversion rates. Figure 11 shows the relationship between E and A of rapeseed stover and α at different conversion rates. Figure 12 shows the relationship between E and A of cotton stover and α at different conversion rates. Regarding the relationship between E and α, the trends in activation energy under different methods were largely consistent for biomass of the same type. For corn stover, activation energy increased as α rose from 0.10 to 0.30: from 91.20 to 102.31 kJ mol-1 (FWO), from 82.55 to 93.08 kJ mol-1 (KAS), and from 82.90 to 93.45 kJ mol-1 (Starink). This indicates the endothermic nature of pyrolysis during the pyrolysis stage. When α increased from 0.30 to 0.90, the activation energy showed a decreasing trend. It decreased from 102.31 to 56.15 kJ mol-1 (FWO), from 93.08 to 43.55 kJ mol-1 (KAS), and from 93.45 to 44.06 kJ mol-1 (Starink). It peaks at α = 0.30. For rapeseed straw, the activation energy decreases as α increases from 0.10 to 0.40. From 110.59 to 102.8 kJ mol-1 (FWO), 102.07 to 93.58 kJ mol-1 (KAS), 102.39 to 93.91 kJ mol-1 (Starink). Activation energy increased as α increased from 0.40 to 0.50. From 102.80 to 107.17 kJ mol-1 (FWO), 93.58 to 97.73 kJ mol-1 (KAS), 93.91 to 98.08 kJ mol-1 (Starink). When α increases from 0.50 to 0.90, the activation energy shows a decreasing trend. From 107.17 to 69.45 kJ mol-1 (FWO), 97.73 to 57.09 kJ mol-1 (KAS), and 98.08 to 57.59 kJ mol-1 (Starink). This fluctuation suggests a complex reaction mechanism during rapeseed straw pyrolysis. For cotton stalks, the activation energy increases as α rises from 0.10 to 0.50. Values range from 99.93 to 110.77 kJ mol-1 (FWO), 91.029 to 100.798 kJ mol-1 (KAS), and 91.38 to 101.20 kJ mol-1 (Starink). When α increased from 0.50 to 0.90, the activation energy showed a decreasing trend. It decreased from 110.77 to 47.17 kJ mol-1 (FWO), from 100.798 to 33.865 kJ mol-1 (KAS), and from 101.20 to 34.40 kJ mol-1 (Starink). It reached its peak at α = 0.50.

Corn stalks. (a) Relationship between E and α; (b) Relationship between A and α.
Figure 10.
Corn stalks. (a) Relationship between E and α; (b) Relationship between A and α.
Rapeseed stalks. (a) Relationship between E and α; (b) Relationship between A and α.
Figure 11.
Rapeseed stalks. (a) Relationship between E and α; (b) Relationship between A and α.
Cotton stalks. (a) Relationship between E and α; (b) Relationship between A and α.
Figure 12.
Cotton stalks. (a) Relationship between E and α; (b) Relationship between A and α.

For the relationship between A and α, the trends in the indexing factor under the KAS and Starink methods are very similar for biomass of the same type, whereas the FWO method differs from them. For corn stover, the indexing factor shows an upward trend when the conversion rate ranges from 0.10 to 0.30. From 1.13E+06 to 1.31E+07 kJ mol-1 (FWO), 1.67E+05 to 1.71E+06 kJ mol-1 (KAS), and 1.80E+05 to 1.86E+06 kJ mol-1 (Starink). This is primarily attributed to the thermal decomposition of cellulose and non-cellulose polysaccharides. However, when the conversion rate is between 0.30 and 0.90, the prefactor shows a decreasing trend. From 1.31E+07 to 4.47E+02 kJ mol-1 (FWO), 1.71E+06 to 2.46E+01 kJ mol-1 (KAS), 1.86E+06 to 2.77E+01 kJ mol-1 (Starink). It peaks at 0.30. For rapeseed straw, the prefactor exhibits a decreasing trend when the conversion rate ranges from 0.10 to 0.40. From 1.28E+08 to 2.25E+07 kJ mol-1 (FWO), 1.91E+07 to 2.85E+06 kJ mol-1 (KAS), 2.05E+07 to 3.07E+06 kJ mol-1 (Starink). When the conversion rate is between 0.40 and 0.50, the prefactor shows an increasing trend. From 2.25E+07 to 5.96E+07 kJ mol-1 (FWO), 2.85E+06 to 7.23E+06 kJ mol-1 (KAS), 3.07E+06 to 7.82E+06 kJ mol-1 (Starink). When the conversion rate is between 0.50 and 0.90, the prefactor shows a decreasing trend. From 5.96E+07 to 1.22E+04 kJ mol-1 (FWO), 7.23E+06 to 7.16E+02 kJ mol-1 (KAS), 7.82E+06 to 8.04E+02 kJ mol-1 (Starink). Peak reached at 0.10. For cotton stalks, the pre-exponential factor shows an increasing trend when the conversion rate is between 0.10 and 0.50. From 7.18E+06 to 7.68E+07 kJ mol-1 (FWO), 1.02E+06 to 8.70E+06 kJ mol-1 (KAS), 1.10E+06 to 9.48E+06 kJ mol-1 (Starink). However, when the conversion rate is between 0.50 and 0.90, the prefactor shows a decreasing trend. From 7.68E+07 to 5.49E+01 kJ mol-1 (FWO), 8.70E+06 to 2.44E+00 kJ mol-1 (KAS), and 9.48E+06 to 2.77E+00 kJ mol-1 (Starink). It peaks at 0.50.

3.4. Thermodynamic parameter analysis

Thermodynamic parameters are crucial for understanding pyrolysis processes. Knowledge of these parameters provides valuable insights into energetics and kinetics, laying the foundation for subsequent industrial design. We employed equations (6), (7), and (8) to determine thermodynamic parameters for α values ranging from 0.1 to 0.9. The results for corn stover, rapeseed straw, and cotton straw are presented in Tables 7-9, respectively. Notably, ΔH values remain positive across all conversion rates. Rapeseed straw exhibits higher ΔH values compared to other straws. Average ΔH values for corn straw are 81.36 kJ mol-1 (FWO), 71.21 kJ mol-1 (KAS), and 71.61 kJ mol-1 (Starink). The average ΔH values for rapeseed straw were 95.89 kJ mol-1 (FWO), 85.89 kJ mol-1 (KAS), and 86.19 kJ mol -1 (Starink). The average ΔH values for cotton stalks were 86.13 kJ mol-1 (FWO), 75.58 kJ mol-1 (KAS), and 76 kJ mol-1 (Starink). As α increased from 0.1 to 0.3, ΔH gradually rose for all three biomasses, indicating the endothermic nature of the pyrolysis process. After this conversion level, ΔH values progressively decreased. At each conversion level, the difference between ΔH and E for the three biomasses ranged from approximately 3 to 10 kJ mol-1. Consequently, an additional 3 to 10 kJ mol-1 of energy was required for product formation. This minor discrepancy also indicates the feasibility of the pyrolysis process.

Table 7. Thermodynamic parameters of corn stalk under different models.
Conversion rate Method FWO
KAS
Starink
ΔH ΔG ΔS ΔH ΔG ΔS ΔH ΔG ΔS
0.10 87.40 168.30 -0.14 78.75 168.77 -0.16 79.10 168.75 -0.16
0.20 96.99 167.82 -0.12 88.02 168.27 -0.14 88.38 168.25 -0.14
0.30 98.51 167.75 -0.12 89.28 168.20 -0.14 89.65 168.18 -0.14
0.40 95.74 167.88 -0.13 86.27 168.36 -0.14 86.65 168.34 -0.14
0.50 94.93 167.92 -0.13 85.21 168.41 -0.15 85.60 168.39 -0.14
0.60 85.38 168.41 -0.14 75.29 168.98 -0.16 75.69 168.95 -0.16
0.70 63.29 169.76 -0.19 52.38 170.61 -0.21 52.82 170.57 -0.21
0.80 57.66 170.18 -0.20 45.88 171.19 -0.22 46.35 171.15 -0.22
0.90 52.35 170.61 -0.21 39.75 171.82 -0.23 40.26 171.77 -0.23
Average 81.36 168.74 -0.15 71.21 169.40 -0.17 71.61 169.37 -0.17
results calculated at 10 °C/min.
Table 8. Thermodynamic parameters of rapeseed stalk under different models.
Conversion rate Method FWO
KAS
Starink
ΔH ΔG ΔS ΔH ΔG ΔS ΔH ΔG ΔS
0.10 106.78 164.08 -0.10 98.26 164.45 -0.12 98.58 164.44 -0.12
0.20 104.00 164.20 -0.11 95.20 164.60 -0.12 95.54 164.58 -0.12
0.30 99.37 164.40 -0.12 90.27 164.83 -0.13 90.67 164.81 -0.13
0.40 98.99 164.42 -0.12 89.77 164.86 -0.13 90.10 164.84 -0.13
0.50 103.36 164.22 -0.11 93.92 164.66 -0.13 94.27 164.64 -0.13
0.60 100.18 164.37 -0.11 90.34 164.83 -0.13 89.87 164.85 -0.13
0.70 97.30 164.50 -0.12 86.32 165.04 -0.14 86.76 165.01 -0.14
0.80 87.43 164.98 -0.14 75.70 165.62 -0.16 76.17 165.59 -0.16
0.90 65.64 166.26 -0.18 53.28 167.17 -0.20 53.78 167.13 -0.20
Average 95.89 164.60 -0.12 85.89 165.12 -0.14 86.19 165.10 -0.14
results calculated at 10 °C/min.
Table 9. Thermodynamic parameters of cotton stalk under different models.
Conversion rate Method FWO
KAS
Starink
ΔH ΔG ΔS ΔH ΔG ΔS ΔH ΔG ΔS
0.10 95.89 168.46 -0.13 86.99 168.90 -0.14 87.34 168.89 -0.14
0.20 101.02 168.22 -0.12 91.82 168.66 -0.13 92.18 168.64 -0.13
0.30 103.43 168.11 -0.11 93.99 168.55 -0.13 94.36 168.53 -0.13
0.40 106.47 167.98 -0.11 96.77 168.42 -0.12 97.15 168.40 -0.12
0.50 106.73 167.96 -0.11 96.76 168.42 -0.12 97.16 168.40 -0.12
0.60 100.47 168.24 -0.12 89.95 168.75 -0.14 90.37 168.73 -0.14
0.70 66.60 170.12 -0.18 55.10 170.97 -0.20 55.56 170.93 -0.20
0.80 51.44 171.27 -0.21 39.03 172.48 -0.23 39.52 172.43 -0.23
0.90 43.13 172.05 -0.22 29.83 173.63 -0.25 30.36 173.56 -0.25
Average 86.13 169.16 -0.14 75.58 169.86 -0.16 76.00 169.83 -0.16
results calculated at 10 °C/min.

The change in Gibbs free energy indicates the spontaneity of a reaction. A lower ΔG value signifies thermodynamically favorable conditions, indicating the feasibility of converting biomass into desired products. A higher ΔG value indicates a reaction prone to extreme heat flows and randomness, suggesting less favorable conditions. Under FWO conditions, the average ΔG values for corn stover, rapeseed stover, and cotton stover were estimated at 168.74 kJ mol-1, 164.6 kJ mol-1, and 169.16 kJ mol-1, respectively. Under KAS conditions, the average ΔG values for corn stover, rapeseed straw, and cotton stalk were estimated at 169.4 kJ mol-1, 165.12 kJ mol-1, and 169.86 kJ mol-1, respectively. Under Starink conditions, the average ΔG values for corn stover, rapeseed straw, and cotton stalk were estimated at 169.37 kJ mol-1, 165.1 kJ mol-1, and 169.83 kJ mol-1, respectively. It can be observed that the ΔG values for the three straws are nearly identical, independent of the model. The positive ΔG values indicate the non-spontaneous nature of pyrolysis [33].

The entropy change represents the degree of disorder within a system. During pyrolysis, ΔS denotes the entropy difference between the products and reactants involved in the thermal decomposition of organic materials. All three types of straw exhibited negative ΔS values under different models [33]. This confirms that the disorder level of pyrolysis products is lower than that of the initial reactants. Simultaneously, these negative values indicate that pyrolysis products in the activated state possess a more ordered structure than pre-pyrolysis materials. This also suggests a low incidence of randomness during the dissociation of complex chemical bonds in actual reactants. The average ΔS values for corn stover were -0.15 kJ K-1 (FWO), -0.17 kJ K-1 (KAS), and -0.17 kJ K-1 (Starink). For rapeseed straw, the average ΔS values were -0.12 kJ K-1 (FWO), -0.14 kJ K-1 (KAS), and -0.14 kJ/K (Starink). The average ΔS values for cotton stalks were -0.14 kJ K-1 (FWO), -0.16 kJ K-1 (KAS), and -0.16 kJ K-1 (Starink).

The changes in ΔH, ΔG, and ΔS during pyrolysis of the three straw types are closely related to the reaction mechanisms and provide key theoretical basis for process optimization in bioenergy production. From the perspective of reaction mechanism correlation, ΔH, as the enthalpy change of reaction, directly reflects the energy demand of the pyrolysis reaction. The ΔH of rapeseed straw is significantly higher than that of corn and cotton straw, consistent with the higher apparent activation energy mentioned earlier. This stems fundamentally from its higher lignin content. The stable benzene ring structure of lignin requires the breaking of more chemical bonds during pyrolysis, increasing the energy input required for the reaction. This aligns with the F3 three-dimensional diffusion mechanism. A high enthalpy change indicates that the reaction system must overcome stronger intermolecular forces to achieve three-dimensional diffusion and conversion of the material. ΔG represents Gibbs free energy; its positive values indicate that pyrolysis of all three straw types is a non-spontaneous reaction requiring external heat input [34]. The lower ΔG value for rapeseed straw further confirms its lower energy barrier and greater spontaneity. ΔS, as entropy change, reflects alterations in system disorder. The larger ΔS value for rapeseed straw indicates more severe molecular structural disruption during pyrolysis, with rapid breakdown of the cross-linked structures of cellulose, hemicellulose, and lignin, yielding more small-molecule products. This aligns with the mechanism in the F3 model where “product diffusion dominates the reaction process.” Greater entropy increase implies stronger diffusion driving force, facilitating deeper conversion of the reaction.

From a bioenergy production perspective, ΔH provides guidance for heating process optimization: Corn and cotton straw exhibit lower ΔH values, allowing for energy savings through moderate-to-low heating intensities. Rapeseed straw requires moderately increased heating intensity to ensure complete reaction. ΔG guides reaction condition control. Optimizing pyrolysis temperature or introducing catalysts can lower the reaction free energy barrier, enhance spontaneity, and improve conversion efficiency. ΔS reflects product distribution tendencies. Rapeseed straw exhibits higher ΔS, favoring the formation of gaseous hydrocarbons (e.g., methane, ethylene), making it suitable for targeted biogas production. Corn and cotton straw show relatively moderate ΔS, yielding primarily biochar and bio-oil. Processes can be specifically optimized to enhance target product yields. The synergistic analysis of these three factors provides precise thermodynamic guidance for large-scale bioenergy production from crop residues, facilitating the minimization of energy consumption and maximization of product value.

3.5. Model method

The pyrolysis kinetics of three straw samples was analyzed using the Coats-Redfern model method to elucidate the pyrolysis reaction mechanism. Figures 13-15 show the linear fitting curves of the C-R model for corn stalks, rapeseed stalks, and cotton stems, respectively, under a heating rate of 10°C/min. Seventeen kinetic models were selected for each temperature range. The optimal model for each sample was determined by comparing the correlation coefficients of different reaction models. A correlation coefficient closer to 1 indicates a more ideal fit.

Linear fitting plot of corn stalk C-R method.
Figure 13.
Linear fitting plot of corn stalk C-R method.
Linear fitting plot of rapeseed stalk C-R method.
Figure 14.
Linear fitting plot of rapeseed stalk C-R method.
Linear fitting plot of cotton stalk C-R method.
Figure 15.
Linear fitting plot of cotton stalk C-R method.

The results indicate that for all three straw types, the optimal ideal kinetic model under a heating rate of 10°C/min is identical: the F3 model. For the main pyrolysis stage, the F3 model was selected as the optimal reaction mechanism. The linear correlation coefficients for corn stover, rapeseed straw, and cotton stalk were 0.98471, 0.94482, and 0.97261, respectively.

The F3 model determined by C-R method is consistent with most of the research conclusions of biomass pyrolysis mechanism. In the pyrolysis of apple residue, the chemical reaction of the constituent elements is the dominant one, which follows the F3 model [35]. In addition, the pyrolysis mechanism of biomass is regulated by the proportion of raw material components, and the higher the proportion of lignin, the more complex diffusion behavior is likely to be presented [36]. The higher values of ΔH, ΔS and E of rapeseed straw in this study are consistent with the conclusion of the literature, which further confirms the influence of component differences on the reaction mechanism.

The core principle of the F3 model is that during the primary reaction phase of straw pyrolysis (the volatilization and decomposition stage), heat and product molecules must diffuse through the three-dimensional pore structure within the straw. The reaction rate is thus controlled by the diffusion process rather than the chemical reaction itself. Experimental results show that all three straw types conform to the F3 model at a heating rate of 10°C/min. This indicates that during pyrolysis, volatile components generated by the decomposition of internal fiber structures (cellulose, hemicellulose) must escape through the three-dimensional porous medium channels. This diffusion process becomes the key rate-limiting step, highly consistent with the diffusion-controlled mechanism described by the F3 model. Further validation indicates that the F3 mechanism serves as the optimal kinetic model for biomass microwave pyrolysis at all heating rates ranging from 10 to 40°C/min.

3.6. Characterization of solid product

A scanning electron microscope (JSM-6490LV) with 15 kV acceleration voltage was used to analyze the surface morphology of biomass straw before and after pyrolysis, as shown in Figure 16. The analysis revealed that corn straw exhibited a loose fiber bundle stacking structure with a rough surface and irregular gaps between fibers Figure 16(a). The high cellulose and hemicellulose content resulted in a loose fiber arrangement with minimal pores. Rapeseed straw showed the highest lignification degree, featuring a dense and smooth surface where lignin acted as a binder filling fiber gaps, leaving almost no visible voids Figure 16(b). Cotton straw formed a fine, elongated cellulose fiber network with a smooth surface, where narrow and regular gaps formed at fiber intersections Figure 16(c). The high cellulose purity contributed to strong fiber toughness and a well-structured morphology. Under mild pyrolysis conditions, the pore size remained small, likely due to persistent blockage by volatile matter and ash.

SEM images of (a) corn straw, (b) rapeseed straw, (c) cotton straw, (d) corn straw bio-char, (e) rapeseed straw bio-char, (f) cotton straw bio-char.
Figure 16.
SEM images of (a) corn straw, (b) rapeseed straw, (c) cotton straw, (d) corn straw bio-char, (e) rapeseed straw bio-char, (f) cotton straw bio-char.

Pyrolysis-induced decomposition and volatilization of organic components lead to morphological transformations in surface structures. Corn stalk biochar develops a moderately developed uniform mesoporous structure with thick, structurally stable pore walls, reflecting the pyrolysis characteristics of moderate cellulose/hemicellulose decomposition Figure 16(d). Rapeseed stalk biochar exhibits a honeycomb-like interconnected porous structure with dense and highly interconnected pores, where volatilization occurs during lignin pyrolysis while aromatic structures prevent pore collapse Figure 16(e). Cotton stalk biochar retains elongated interconnected channels and traces of residual fibrous frameworks, with its fibrous structure remaining partially intact during pyrolysis, allowing volatilization to occur directionally through fiber gaps to form characteristic pores Figure 16(f). These morphological differences directly correlate with their pyrolysis kinetics and thermodynamic properties, providing microscopic structural evidence for targeted biochar applications.

4. Conclusions

This study investigated the thermodynamic behavior of three typical crop straws (corn stover, rapeseed straw, and cotton straw) during pyrolysis. The main conclusions are as follows:

Straw exhibits the highest carbon content and relatively low sulfur content among the three crops, making it a relatively clean energy source. Studies on straw pyrolysis characteristics indicate that the process comprises three stages: drying and dehydration, main pyrolysis, and carbonization. As the heating rate increases, both TG and DTG curves shift toward higher temperatures. TG analysis indicates that the pyrolysis performance of straw is influenced by the heating rate. Corresponding characteristic temperatures increase with higher heating rates. For the T-α curve of straw, at a constant heating rate, α increases with rising pyrolysis temperature. At a constant pyrolysis temperature, α decreases with increasing heating rate;

The apparent activation energy for corn stalk pyrolysis was calculated using the FWO method, yielding a range of 56.15–102.31 kJ mol-1 with an average value of 85.16 kJ mol-1; Using the KAS method, E ranged from 43.55 to 93.08 kJ mol-1, with an average value of 75 kJ mol-1; using the Starink method, E ranged from 44.06 to 93.45 kJ mol-1, with an average value of 75.41 kJ mol-1; For rapeseed straw, the FWO method yielded an E range of 69.45–110.59 kJ mol-1 with an average of 99.7 kJ mol-1; the KAS method yielded an E range of 57.09–102.07 kJ mol-1 with an average of 89.71 kJ mol-1; Calculations using the Starink method yielded an E range of 57.59–102.39 kJ mol-1, with an average value of 90 kJ mol-1; For cotton straw, the FWO method yielded an E range of 47.17–110.77 kJ mol-1, with an average value of 90.17 kJ mol-1; Calculations using the KAS method yielded an E range of 33.87–100.80 kJ mol-1, with an average value of 79.62 kJ mol-1; calculations using the Starink method yielded an E range of 34.40–101.196 kJ mol-1, with an average value of 80.04 kJ mol-1;

The apparent activation energies and pre-exponential factors of different samples showed significant differences, indicating varying degrees of pyrolysis difficulty. Comparison revealed that rapeseed straw exhibited the highest apparent activation energy, suggesting its pyrolysis process requires higher energy input and is relatively challenging. In contrast, corn stover exhibits the lowest apparent activation energy, suggesting easier pyrolysis. Regarding pre-exponential factors, rapeseed straw has the highest A value, indicating a faster reaction rate. The relatively lower A values for other samples reflect differences in their reaction rates. Overall, these pyrolysis kinetic parameters provide a robust theoretical basis for further research into the pyrolysis characteristics and resource utilization of straw. Finally, based on thermogravimetric data analysis, the optimal models describing the second stage of straw pyrolysis were determined to be the F3 model;

For thermodynamic parameter analysis, the ΔH values for all three straw types remained positive across different models and conversion rates. Compared to other straws, rapeseed straw exhibited higher ΔH values. ΔG values were also positive, indicating the non-spontaneous nature of pyrolysis. Simultaneously, it can be observed that the ΔG values for all three straw types are nearly identical, independent of the model. Conversely, ΔS values are all negative. These negative values indicate that pyrolysis products in the activated state possess a more ordered structure than before pyrolysis. They also suggest that the probability of random events occurring during the dissociation of complex chemical bonds in actual reactants is very low.

For future research, the composition of straw and specific reactions during pyrolysis remains to be investigated. Subsequently, we will incorporate the obtained kinetic parameters into the Aspen Plus reactor module to predict bio-oil yield and component distribution across different scales (from laboratory to ten-thousand-ton level). This will advance the industrial application of straw pyrolysis technology, support agricultural and rural production, and enhance the resource utilization rate of straw.

Acknowledgment

This project was partly supported by National Natural Science Foundation of China (grant no. 42271301).

CRediT authorship contribution statement

Hui Fang: Methodology, investigation, writing – original draft. Rui-Shen Xie: Methodology, investigation. Li-Xia Wu: Methodology, investigation. Mei-Feng Chen: Writing – review & editing. Qian-Feng Zhang: Conceptualization, investigation, supervision, writing – review & editing.

Declaration of competing interest

There are no conflicts of interest.

Data availability

Data will be made available on request.

Declaration of generative AI and AI-assisted technologies in the writing process

The authors confirm that there was no use of artificial intelligence (AI)-assisted technology for assisting in the writing or editing of the manuscript and no images were manipulated using AI.

References

  1. , , , , . Biomass pyrolysis: Past, present, and future. Environment, Development and Sustainability. 2020;22:17-32. https://doi.org/10.1007/s10668-018-0200-5
    [Google Scholar]
  2. , , , , , , , . A review of recent advances in biomass pyrolysis. Energy & Fuels. 2020;34:15557-15578. https://doi.org/10.1021/acs.energyfuels.0c03107
    [Google Scholar]
  3. , , , , , . Comparative study of pyrolysis behaviors of corn stalk and its three components. Journal of Analytical and Applied Pyrolysis. 2013;104:185-193. https://doi.org/10.1016/j.jaap.2013.08.005
    [Google Scholar]
  4. , , , , , . Investigating the potential for energy, fuel, materials and chemicals production from corn residues (cobs and stalks) by non-catalytic and catalytic pyrolysis in two reactor configurations. Renewable and Sustainable Energy Reviews. 2009;13:750-762. https://doi.org/10.1016/j.rser.2008.01.004
    [Google Scholar]
  5. , , , , , . Study on the co-pyrolysis of oil shale and corn stalk: Pyrolysis characteristics, kinetic and gaseous product analysis. Journal of Analytical and Applied Pyrolysis. 2022;163:105456. https://doi.org/10.1016/j.jaap.2022.105456
    [Google Scholar]
  6. , , , , , , . The study of combustion characteristics of corn stalks and cobs via TGA-DTG-DSC analysis. IOP Conference Series: Earth and Environmental Science. 2019;354:012130. https://doi.org/10.1088/1755-1315/354/1/012130
    [Google Scholar]
  7. , , , . TG-FTIR-MS study of synergistic effects during co-pyrolysis of corn stalk and high-density polyethylene (HDPE) Energy Conversion and Management. 2019;181:202-213. https://doi.org/10.1016/j.enconman.2018.11.065
    [Google Scholar]
  8. , , , , . Comparative pyrolysis behaviors and reaction mechanisms of hardwood and softwood. Energy Conversion and Management. 2017;132:102-109. https://doi.org/10.1016/j.enconman.2016.11.016
    [Google Scholar]
  9. , , , , , , , , . In-depth understanding of high temperature and low residence time on the corn straw rapid pyrolysis char structure evolution. Fuel. 2025;379:133052. https://doi.org/10.1016/j.fuel.2024.133052
    [Google Scholar]
  10. , , , , , . Influence of temperature and particle size on structural characteristics of chars from beechwood pyrolysis. Journal of Analytical and Applied Pyrolysis. 2018;130:127-134. https://doi.org/10.1016/j.jaap.2018.01.018
    [Google Scholar]
  11. , , , . Influence of pyrolysis temperature on physicochemical properties of biochar obtained from the fast pyrolysis of pitch pine (Pinus rigida) Bioresource Technology. 2012;118:158-162. https://doi.org/10.1016/j.biortech.2012.04.094
    [Google Scholar]
  12. , , , , . Effect of pyrolysis temperature on char structure and chemical speciation of alkali and alkaline earth metallic species in biochar. Fuel Processing Technology. 2016;141:54-60. https://doi.org/10.1016/j.fuproc.2015.06.029
    [Google Scholar]
  13. , , . Study on the pyrolysis properties of corn straw by TG–FTIR and TG–GC/MS. Journal of Thermal Analysis and Calorimetry. 2021;143:3783-3791. https://doi.org/10.1007/s10973-020-09778-6
    [Google Scholar]
  14. , , , , , . Syngas composition and ash characteristics of corn straw under a CO2 atmosphere. Biomass and Bioenergy. 2022;166:106630. https://doi.org/10.1016/j.biombioe.2022.106630
    [Google Scholar]
  15. , , . Devolatilization non-isothermal kinetic analysis of agricultural stalks and application of TG-FT/IR analysis. Journal of Thermal Analysis and Calorimetry. 2012;107:1177-1189. https://doi.org/10.1007/s10973-011-1959-x
    [Google Scholar]
  16. , , , , , , , . Infrared heated pyrolysis of corn stover: Determination of kinetic and thermodynamic parameters. Journal of Analytical and Applied Pyrolysis. 2021;158:105273. https://doi.org/10.1016/j.jaap.2021.105273
    [Google Scholar]
  17. , . Pyrolysis of corn cob: Physico-chemical characterization, thermal decomposition behavior and kinetic analysis. Chemical Product and Process Modeling. 2021;16:117-127. https://doi.org/10.1515/cppm-2020-0048
    [Google Scholar]
  18. , , , , , . Thermodynamic and kinetic analyses of corn crop residue. In: Proceedings of the 2nd International Conference on Water Energy Food and Sustainability (ICoWEFS 2022). Cham: Springer International Publishing; . p. :206-217. https://doi.org/10.1007/978-3-031-26849-6_22
    [Google Scholar]
  19. , . Insight of biopolymeric degradation and kinetic triplets of corn cob and its biochar using thermogravimetric analysis. Thermochimica Acta. 2025;748:179990. https://doi.org/10.1016/j.tca.2025.179990
    [Google Scholar]
  20. , , , , . Thermal and kinetic analyzing of pyrolysis and combustion of self-heating biomass particles. Process Safety and Environmental Protection. 2021;151:39-50. https://doi.org/10.1016/j.psep.2021.05.011
    [Google Scholar]
  21. , , , , , , , , , . Chemical characterization and oxidative potential of particles emitted from open burning of cereal straws and rice husk under flaming and smoldering conditions. Atmospheric Environment. 2017;163:118-127. https://doi.org/10.1016/j.atmosenv.2017.05.037
    [Google Scholar]
  22. , , , , , . Combustion of pellets produced from corn straw and cashew nut shell in residential and industrial scenarios: Chemical, thermal and emission analyses. Waste Management & Research : The Journal of the International Solid Wastes and Public Cleansing Association, ISWA. 2023;41:1486-1495. https://doi.org/10.1177/0734242X221135518
    [Google Scholar]
  23. , , , . Comparative studies on the pyrolysis of cellulose, hemicellulose, and lignin based on combined kinetics. Journal of the Energy Institute. 2019;92:27-37. https://doi.org/10.1016/j.joei.2017.12.003
    [Google Scholar]
  24. , , , . Kinetic analysis of rice husk pyrolysis using kissinger-akahira-sunose (KAS) method. Procedia Engineering. 2016;148:1247-1251. https://doi.org/10.1016/j.proeng.2016.06.486
    [Google Scholar]
  25. , , , , . Kinetic study of regional agro-industrial wastes pyrolysis using non-isothermal TGA analysis. Applied Thermal Engineering. 2016;106:1157-1164. https://doi.org/10.1016/j.applthermaleng.2016.06.084
    [Google Scholar]
  26. , , , , , , . Characteristics of biomass fast pyrolysis in a wire-mesh reactor. Fuel. 2017;200:225-235. https://doi.org/10.1016/j.fuel.2017.03.070
    [Google Scholar]
  27. , , , . Effect of torrefaction on combustion kinetic parameters of groundnut shell and maize stalk. Biomass Conversion and Biorefinery. 2025;15:18311-18329. https://doi.org/10.1007/s13399-025-06595-5
    [Google Scholar]
  28. , , , , , , , . Co-hydrothermal carbonization of corn stalk and swine manure: Combustion behavior of hydrochar by thermogravimetric analysis. Bioresource Technology. 2019;271:75-83. https://doi.org/10.1016/j.biortech.2018.09.100
    [Google Scholar]
  29. , , , , . Pyrolysis and co-pyrolysis of cattle manure, rape straw, and their blend: Physicochemical characterization, kinetic triplets, reaction mechanism, and thermodynamic analysis. Energy. 2024;292:130520. https://doi.org/10.1016/j.energy.2024.130520
    [Google Scholar]
  30. , , , , . A study on the characteristics and kinetic of co-catalytic pyrolysis with rape straw and ABS waste plastics. Renewable Energy. 2025;242:122450. https://doi.org/10.1016/j.renene.2025.122450
    [Google Scholar]
  31. , , , , , , , , , , , . Synergistic effects and kinetics in co-pyrolysis of waste tire with five agricultural residues using thermogravimetric analysis. Journal of Energy Resources Technology. 2023;145:121501. https://doi.org/10.1115/1.4062826
    [Google Scholar]
  32. , , , , . Understanding pyrolysis mechanisms of corn and cotton stalks via kinetics and thermodynamics. Journal of Analytical and Applied Pyrolysis. 2022;164:105521. https://doi.org/10.1016/j.jaap.2022.105521
    [Google Scholar]
  33. , , , , , , , , . Co-pyrolysis of corn stalk and high-density polyethylene with emphasis on the fibrous tissue difference on thermal behavior and kinetics. The Science of the Total Environment. 2024;957:177847. https://doi.org/10.1016/j.scitotenv.2024.177847
    [Google Scholar]
  34. , . Investigation of thermodynamic parameters in the pyrolysis conversion of biomass and manure to biochars using thermogravimetric analysis. Bioresource Technology. 2013;146:485-493. https://doi.org/10.1016/j.biortech.2013.07.086
    [Google Scholar]
  35. , . Kinetic, thermodynamic, and reaction mechanism analysis of apple pomace (Malus Domestica Borkh) for bioenergy potential using single and multi-step (Deconvolution) approaches. Bioresource Technology Reports. 2025;32:102401. https://doi.org/10.1016/j.biteb.2025.102401
    [Google Scholar]
  36. , , , . Impact of feedstock composition on pyrolysis of low-cost, protein- and lignin-rich biomass: A review. Journal of Analytical and Applied Pyrolysis. 2020;147:104780. https://doi.org/10.1016/j.jaap.2020.104780
    [Google Scholar]
Show Sections