Ajay Sharma* and
Bikash Mohanty
Department of Chemical Engineering, Indian Institute of Technology, Roorkee, India 247667. E-mail: asharma@ch.iitr.ac.in
First published on 12th April 2021
For better utilization of 11 million tonnes of Mangifera indica wood (MIW) sawdust produced annually in India, the present study was planned for its characterization followed by determination of its pyrolysis kinetics from TGA data under a N2 atmosphere. The characterization process included proximate-, ultimate-, biopolimeric components-, and heating value-analysis, as well as TG/DTG analysis. The distributed activation energy (DAE)- and Starink-methods were implemented on non-isothermal thermograms to compute the isoconversional values of activation energy for the pyrolysis of MIW. Further, the reaction mechanism for the pyrolysis of MIW was predicted using the Coats–Redfern (C–R) model-fitting method. Two distinct pyrolysis regions, region-I from 0.05–0.5 and region-II from 0.51–0.7, were observed in the complete conversion ranges. The estimated activation energy for region-I ranged from 143.03 to 176.46 kJ mol−1 with an average value of 157.12–157.97 kJ mol−1 and that of region-II varied between 143.03 and 161.68 kJ mol−1 with an average of 151.51–152.45 kJ mol−1. The one-dimensional diffusion model (D1) followed by the five and a half reaction order model (F5.5) were recommended to describe the pyrolysis reaction mechanism of MIW for the two above regions, respectively. Further, the activation energies obtained via the DAE and Starink methods were used for the computation of thermodynamic parameters such as frequency factor, and change in-enthalpy, -entropy, and -Gibbs free energy.
Parikka estimated the total worldwide energy potential from biomass, on a sustainable basis, to be 104 EJ per year, of which woody biomass, energy crops, and straw constituted 40.1%, 36%, and 16.6%, respectively.1 Only about 40% of potential biomass energy is currently utilized. The present global energy demand is about 470 EJ per year.1 Forest residues, which remain largely unutilized, present solid waste disposal problems. Among these residues, wood sawdust, which is created as byproducts of the wood processing at wood cutting and furniture making industry, are generally generated in large quantity every year.2 India is producing around 149 million cubic meters (≈100 million tonnes) of Mangifera indica wood per year,3 which when converted into planks will produce 11 million tonnes of sawdust. This sawdust is typically used in making of particle board and road laying. Still, a large amount of its waste is dumped around roadside, water bodies, and/or burn directly in an open environment, which creates anthropogenic hazard for all the living species and thus are an unacceptable solution. Mangifera indica wood is increasingly being used in furniture making, due to its pros as hardwood furniture material as well as its sustainability quotient. Mangifera indica species are found all around the world (major contributors are Bangladesh, Burma, Cambodia, China, Hawaii (US), India, Indonesia, Laos, Malaysia, Papua New Guinea, Philippines, Thailand, and Vietnam). Mangifera indica wood is readily available and is affordable compared to other timber species.
Reaction kinetics, which includes the computations of kinetic triplets (the activation energy, pre-exponential factor and reaction models), is essentially required in order to understand the chemistry behind the thermal degradation of material and further to transform laboratory-scale studies to large-scale reactor conditions.4 The procedures used for the computation of kinetic parameters (A, E, and n) of material degradation during solid–gas heterogeneous reactions include model-fitting and model-free methods, which have been extensively used in the thermal decomposition of solids, thermal degradation of polymers and crystallization of glasses.5 TG-DTG coupled FTIR and/or GC-MS is widely used to extract information about material degradation and further to interpret the products released from the pyrolyzer unit.6 The reaction kinetics provide deeper insights into the mechanism of solid thermal degradation which is important for the proper selection-, optimization- and operation – of a reactor.7 Moreover, kinetic analysis of pyrolysis reaction further opens the path for the estimation of various thermodynamic parameters enthalpy (ΔH), Gibbs free energy (ΔG), and entropy (ΔS) of the pyrolysis process.8 The thermodynamic parameters are essentially needed to design an efficient pyrolysis system and also for choosing the target biofuels.9
Amongst all thermochemical processes such as combustion, gasification, and pyrolysis, pyrolysis is the least energy-consuming process thus, the assessment of thermodynamic parameters on pyrolysis conditions is imperative. In all, kinetic parameters, reaction mechanism and thermodynamic parameters (change in-enthalpy, -Gibbs free energy, and -entropy) are important for the design, optimization, and scaling of process reactor.
Even though India is the largest producer of mango (Mangifera indica) wood in the world, scientific management of waste mango wood created during processing needs special attention. From literature review, it has been observed that there is hardly any research work conducted to study pyrolysis kinetics and its thermodynamic analysis for Mangifera indica wood (MIW). Also, a few research papers are available in an open literature that provides systematic kinetic and thermodynamic data of wood pyrolysis with detailed explanations of phenomena. Hence, the present study was planned to investigate the characterization of MIW sawdust, determination of pyrolysis kinetics and reaction mechanism. In addition, the thermodynamic parameters were computed to assess the feasibility of MIW pyrolysis for effective process designs and operations.
Moisture content (MC) | Volatile matter (VM) | Ash | |
---|---|---|---|
a Where m1 = mass of the sample after moisture removal, m2 = mass of sample after the removal of moisture and volatile matter, and m3 = mass of ash content present in the sample.b All mass should be measured after cooling down the sample in a desiccator up to room temperature. | |||
Standard | ASTME871 | ASTME872 | ASTMD1102 |
Sample mass (initial) | 1 g | m1 g | m2 g |
Temperature | 103 ± 2 °C | 925 ± 5 °C | 580–600 °C |
Holding time | Until constant mass | 7 min | 4 h |
Crucible | Silica spherical | Silica cylindrical | Silica cylindrical |
Sample mass (final) | m1 g | m2 g | m3 g |
Calculation (%) | (1 − m1) × 100 | (m1 − m2)×100 | m3 × 100 |
Fixed carbon (mass%) = 100 − MC − VM − ash = (m2 − m3) × 100 |
Further, the elemental compositions (CHNS/O) of MIW sawdust were determined using an elemental analyzer (Perkin Elmer 2400, USA). The elemental compositions which include the determination of the mass% of carbon, hydrogen, nitrogen, sulfur, and oxygen content were then used for the computation of calorific value. The calorific value of material can be categorized as higher and lower heating values (HHV and LHV). The HHV and LHV measure the amount of energy which can be obtained from the unit quantity of material. The HHV and LHV were computed using the analytical expressions given by Channiwala and Parikh10 and Basu,11 respectively. In addition to it, mass% of extractives, hemicellulose, cellulose, and lignin were determined using Li et al.12 The detailed experimental procedure is given elsewhere.13
![]() | (1) |
Therefore,
![]() | (2) |
Further, a constant term β (=dT/dt) called heating rate is introduced as a conversion factor, used whenever required for the transformation of dynamic form into non-isothermal form. As in the present case, for non-isothermal TG kinetics, eqn (2) converts as below:
![]() | (3) |
Integrating above on both sides:
The second term in R.H.S is comparatively small. So, neglecting low-temperature end of the interval, the general form of all the integral isoconversional kinetic methods,
![]() | (4) |
![]() | (5) |
From eqn (4),
![]() | (6) |
Assuming (1 − 2RT/Ea) ≈ 1 and taking loge of eqn (6), the final expression becomes:
ln[g(α)/Tα2] = −(Ea/RTα) + ln(AR/βEa) | (7) |
Based on mechanistic assumptions, reaction models are divided into nucleation, geometrical contraction, diffusion, and reaction order.18 Under nucleation models, the formation and growth of nuclei are considered to be rate-limiting. In geometrical contraction, the progress of product layer from the surface to the inner boundary is considered to be rate-limiting. In diffusion, the diffusion of reactant or product is considered to be rate-limiting. And finally, in reaction order based model, the rate law is based on homogeneous kinetics. Based on the above mechanisms, the integral form of the reaction model [g(α)] was decided. Table 2 lists different models as a function of the extent of conversion. Here, g(α) is approximated to defined linear or non-linear expressions based on conversion (α) and order of reaction (n). Corresponding to each kinetic mechanism, a global value of activation energy (Ea) and pre-exponential factor (A) was obtained.
![]() | (8) |
ΔHα = Eα − RTα | (9) |
![]() | (10) |
![]() | (11) |
Moreover, the energy barriers (EB) for the formulation of activated complex from the initial reactants were also computed.
Mathematically,
EB = Eα − ΔHα. | (12) |
MIW | |
---|---|
a O2 content (mass%) = 100 − C (mass%) − H (mass%) − N (mass%) − S (mass%) − ash (mass%) − moisture (mass%).b Empirical formula of MIW was developed using its ultimate analysis data. | |
Proximate analysis (mass%) | |
Moisture content | 7.06 |
Volatile matter | 73.38 |
Ash | 4.47 |
Fixed carbon | 15.09 |
![]() |
|
Ultimate analysis (mass%) | |
C | 38.98 |
H | 5.065 |
N | 0.08 |
S | 0.00 |
Oa | 44.345 |
HHV (MJ kg−1) | 14.89 |
LHV (MJ kg−1) | 13.64 |
H/C molar ratio | 1.559 |
O/C molar ratio | 0.853 |
Empirical formulab | CH1.559O0.853N0.002 |
![]() |
|
Biopolimeric components (mass%) | |
Extractives | 5 |
Hemicellulose | 25 |
Cellulose | 46.25 |
Lignin | 23.75 |
The proximate analysis of MIW shows volatiles and fixed carbon as 73.38 mass% and 15.09 mass% with moisture and ash content as 7.06 mass% and 4.47 mass%, respectively. The adequate amount of volatiles (73.38 mass%) in MIW facilitates the thermal conversion process, as they lead to more combustible gases during the conversion, which further converts to bio-oil by the condensation process. The fixed carbon in MIW (15.09 mass%) is also an important element as it gets converted into biochar during the pyrolysis.21 Further, the moisture content of 7.06 mass%, in MIW can be considered as low,22 meaning that it can be removed with a small amount of additional energy. The low ash content in MIW (4.47 mass%) indicates its fast thermal degradation since with high ash content biomass will tend to form slag and may cause fouling problems in many thermal conversion conveying parts. The moisture and ash content found in proximate analysis refers to the same moisture and ash as reported in the ultimate analysis. However, in the proximate analysis, obtained fixed carbon as 15.09 mass% differs from the reported carbon in the ultimate analysis which is 38.98 mass%. The plausible reason for this difference is that %C measured in ultimate analysis is the total mass% of carbon in MIW, including the carbon in volatile matter. Also, the sulfur content in MIW is 0 mass%. This indicates its favorability as there are no concerns related to SO2 emission & equipment corrosion when using it as a fuel. The HHV and LHV of MIW were computed as 14.89 MJ kg−1 and 13.64 MJ kg−1, respectively. The empirical formula of MIW was obtained as CH1.559O0.853N0.002. Moreover, the biopolimeric components of MIW such as extractives, hemicellulose, cellulose and lignin (in mass%) were found as 5, 25, 46.25 and 23.75, respectively.
![]() | ||
Fig. 1 (a) Mass loss (%) of MIW with rise in temperature, (b) mass loss rate (mg min−1) of MIW with rise in temperature. |
Fig. 2a and b depict the regression line plots of DAE and Starink methods, respectively. Fig. 2a and b shows that for conversions between 0.05 and 0.7, the correlation coefficients (R2) for above linear fit methods are close to unity. In the present study, depending upon the heating rates [β(1–4) = 10, 20, 30, and 40 °C min−1], the kinetic parameters were computed in the conversion range 0.05 to 0.7 which corresponds to a temperature range of 200 °C to 395 ± 15 °C. After a fractional conversion of 0.7, the computed activation energies (Eα) were not taken into consideration as these were unreliable due to low value of correlation coefficient (R2). The average activation energy (in kJ mol−1) computed using DAE and Starink methods were found to be 154.86 and 155.74, respectively. Table 4 enlists the comparison of average activation energy of MIW with other biomass reported in the literature using isoconversional methods. The matching of average E-values of MIW with other biomass such as poplar wood,26 eucalyptus wood,27 European beech wood,28 and sal wood29 makes it suitable for the co-pyrolysis with several other biomass feedstocks.
Biomass feedstock | Heating rates (°C min−1) | Average activation energy (kJ mol−1) | References |
---|---|---|---|
MIW | 10, 20, 30, 40 | 154.86–155.74 | Present work |
Dalbergia sissoo wood | 5, 10, 15, 20, 30 | 136.06–147.75 | 30 |
Pine wood | 5, 10, 20, 40 | 181.00–181.28 | 31 |
Babool wood | 10, 15, 20 | 91.6 ± 4.7 to 102.7 ± 6.06 | 32 |
Wood sawdust | 5, 10, 20 | 164.24–173.41 | 33 |
Poplar wood | 2, 5, 10, 15 | 157.27–158.58 | 26 |
Eucalyptus wood | 10, 40 | 108.39–192.93 | 27 |
European beech wood | 5, 10, 15, 25, 35, 50 | 157.20–185.40 | 28 |
Pine sawdust | 5, 10, 15, 20, 25 | 168.58–206.62 | 29 |
Sal wood | 148.44–181.53 | ||
Areca nut husk | 160.45–184.61 | ||
Olive wood | 5, 10, 20, 40, 100 | 54.05–116.78 | 34 |
Invasive Reed Canary | 10, 20, 30, 40 | 159.61–161.29 | 35 |
In reference to Fig. 3, For 0.05 ≤ α ≤ 0.1 (avg. temperature from 270 to 288 °C), the decrease in activation energy was noticed which indicates first breaking of the strong hemicellulose–lignin linkage then the comparatively weaker hemicellulose–cellulose and lignin–cellulose linkage.36,37 Then, the activation energy remains almost constant for 0.1<α ≤ 0.2 (289 to 306 °C), indicating no severe bond cracking in this conversion range (Fig. 1b). After that, for 0.2 < α ≤ 0.4 (307 to 334 °C), an increase in activation energies attributes to the thermal cracking of strong bonds (aromatic rings, CO, C(
O)OH, C–C bonds within and between the alkyl chains) inside the biopolimeric components.38 Further, for 0.4 < α ≤ 0.65 (335 to 379 °C), the decrement in activation energies might be due to the breaking of C–OH and aliphatic groups, which requires comparatively less energy to break.39 After that, for 0.65 < α ≤ 0.7 (380 to 395 °C), the increment in activation energy indicates requirement of additional energy to benzene ring, aliphatic groups, methyl groups, and methoxy groups at ortho position of hydroxyl group in MIW.38 In the present study, the pattern of change in activation energies with conversion values is more consistent with the other biowastes such as Prosopis juliflora,40 pine, sal and areca nut sawdust.29 The change in activation energy with respect to conversion (Fig. 3) was due to the complex nature of pyrolysis reaction which shows the involvement of parallel, complex and competitive reactions.41
Region | Reaction model | g(α) | 10 °C min−1 | 20 °C min−1 | 30 °C min−1 | 40 °C min−1 | Average E-values | |||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
E, kJ mol−1 | R2 | E, kJ mol−1 | R2 | E, kJ mol−1 | R2 | E, kJ mol−1 | R2 | E, kJ mol−1 | R2 | |||
I | A2 | [−ln(1 − α)]1/2 | 37.35 | 0.99 | 37.46 | 0.99 | 40.15 | 0.99 | 39.28 | 0.99 | 38.56 | 0.99 |
A3 | [−ln(1 − α)]1/3 | 21.76 | 0.98 | 21.76 | 0.99 | 23.52 | 0.99 | 22.91 | 0.99 | 22.49 | 0.99 | |
A4 | [−ln(1 − α)]1/4 | 13.96 | 0.98 | 13.91 | 0.99 | 15.20 | 0.99 | 14.72 | 0.99 | 14.45 | 0.98 | |
F0 | α | 73.09 | 0.98 | 73.49 | 0.98 | 78.30 | 0.98 | 76.83 | 0.99 | 75.43 | 0.98 | |
F1 | −ln(1 − α) | 84.15 | 0.99 | 84.56 | 0.99 | 90.04 | 0.99 | 88.38 | 0.99 | 86.78 | 0.99 | |
F2 | [1 − 1/(1 − α)]/−1 | 96.56 | 1.00 | 96.96 | 1.00 | 103.19 | 1.00 | 101.32 | 1.00 | 99.51 | 1.00 | |
F3 | [1 − 1/(1 − α)2]/−2 | 110.27 | 1.00 | 110.67 | 1.00 | 117.72 | 1.00 | 115.62 | 1.00 | 113.57 | 1.00 | |
F4 | [1 − 1/(1 − α)3]/−3 | 125.22 | 1.00 | 125.61 | 1.00 | 133.56 | 1.00 | 131.22 | 1.00 | 128.90 | 1.00 | |
F5 | [1 − 1/(1 − α)4]/−4 | 141.30 | 0.99 | 141.68 | 0.99 | 150.60 | 0.99 | 148.00 | 0.99 | 145.39 | 0.99 | |
F5.5 | [1 − 1/(1 − α)4.5]/−4.5 | 149.73 | 0.99 | 150.10 | 0.99 | 159.53 | 0.99 | 156.79 | 0.99 | 154.04 | 0.99 | |
D1 | α2 | 155.61 | 0.98 | 156.62 | 0.99 | 166.35 | 0.99 | 163.49 | 0.99 | 160.52 | 0.99 | |
D2 | (1 − α)ln(1 − α) + α | 162.56 | 0.99 | 163.57 | 0.99 | 173.71 | 0.99 | 170.74 | 0.99 | 167.65 | 0.99 | |
D3 | [1 − (1 − α)1/3]2 | 170.08 | 0.99 | 171.09 | 0.99 | 181.68 | 0.99 | 178.58 | 0.99 | 175.36 | 0.99 | |
D4 | [1 − (2/3)α] − (1 − α)2/3 | 165.07 | 0.99 | 166.07 | 0.99 | 176.36 | 0.99 | 173.36 | 0.99 | 170.21 | 0.99 | |
R2 | 1 − (1 − α)1/2 | 78.45 | 0.99 | 78.86 | 0.99 | 84.00 | 0.99 | 82.43 | 0.99 | 80.94 | 0.99 | |
R3 | 1 − (1 − α)1/3 | 80.32 | 0.99 | 80.72 | 0.99 | 85.97 | 0.99 | 84.38 | 0.99 | 82.85 | 0.99 | |
II | A2 | [−ln(1 − α)]1/2 | 8.76 | 0.91 | 8.98 | 0.96 | 8.90 | 0.97 | 8.00 | 0.95 | 8.66 | 0.95 |
A3 | [−ln(1 − α)]1/3 | 2.37 | 0.61 | 2.43 | 0.79 | 2.35 | 0.84 | 1.70 | 0.67 | 2.21 | 0.73 | |
A4 | [−ln(1 − α)]1/4 | — | — | — | — | — | — | — | — | — | — | |
F0 | α | 12.82 | 0.89 | 13.20 | 0.94 | 13.15 | 0.96 | 12.07 | 0.94 | 12.81 | 0.93 | |
F1 | −ln(1 − α) | 27.93 | 0.96 | 28.61 | 0.98 | 28.56 | 0.99 | 26.89 | 0.98 | 28.00 | 0.98 | |
F2 | [1 − 1/(1 − α)]/−1 | 48.64 | 0.98 | 49.70 | 0.99 | 49.64 | 1.00 | 47.18 | 0.99 | 48.79 | 0.99 | |
F3 | [1 − 1/(1 − α)2]/−2 | 74.47 | 0.99 | 75.99 | 1.00 | 75.91 | 1.00 | 72.46 | 1.00 | 74.71 | 1.00 | |
F4 | [1 − 1/(1 − α)3]/−3 | 104.22 | 0.99 | 106.27 | 1.00 | 106.19 | 1.00 | 101.59 | 1.00 | 104.57 | 1.00 | |
F5 | [1 − 1/(1 − α)4]/−4 | 136.62 | 0.99 | 139.24 | 1.00 | 139.14 | 1.00 | 133.31 | 1.00 | 137.08 | 1.00 | |
F5.5 | [1 − 1/(1 − α)4.5]/−4.5 | 153.46 | 1.00 | 156.39 | 1.00 | 156.28 | 1.00 | 149.80 | 1.00 | 153.98 | 1.00 | |
D1 | α2 | 36.07 | 0.94 | 37.07 | 0.97 | 43.18 | 0.95 | 35.02 | 0.97 | 37.83 | 0.96 | |
D2 | (1 − α)ln(1 − α) + α | 44.41 | 0.96 | 45.58 | 0.98 | 52.11 | 0.96 | 43.22 | 0.98 | 46.33 | 0.97 | |
D3 | [1 − (1 − α)1/3]2 | 54.97 | 0.97 | 56.35 | 0.99 | 63.24 | 0.97 | 53.58 | 0.99 | 57.04 | 0.98 | |
D4 | [1 − (2/3)α] − (1 − α)2/3 | 47.90 | 0.96 | 49.14 | 0.98 | 55.79 | 0.97 | 46.64 | 0.98 | 49.87 | 0.97 | |
R2 | 1 − (1 − α)1/2 | 19.69 | 0.94 | 20.21 | 0.97 | 23.54 | 0.95 | 18.80 | 0.97 | 20.56 | 0.96 | |
R3 | 1 − (1 − α)1/3 | 22.28 | 0.95 | 22.85 | 0.98 | 26.27 | 0.96 | 21.34 | 0.98 | 23.19 | 0.97 |
For better linear trend fitting (R2 > 0.97) of data between ln[g(α)/Tα2] against 1/Tα for all four heating rates, the conversion (α) range can be divided into two regions (Fig. 1a), region-I (0.05 ≤ a ≤ 0.5) and region-II (0.51 ≤ α ≤ 0.7). These ranges were based on primary analysis.
Depending on the degradation mechanism, the functional relationship of the function g(α) changes as can be seen from Table 5. As per C–R method, plots between ln[g(α)/Tα2] against 1/Tα were created for various degradation mechanisms to determine the value of E (slope of the linear trend line in Fig. 4a and b) for each of these degradation mechanisms at different heating rates as reported in Table 5. Then the E-values of a given degradation mechanism for all the four heating rates were added to produce an average E-value to be compared with average E-values determined using DAE and Starink methods. As per C–R method, if for a given degradation mechanism, these two average values match well, then the given degradation mechanism is considered to be the prevailing one.
![]() | ||
Fig. 4 C–R plots for the pyrolysis of MIW; (a) region-I (0.05 ≤ α ≤ 0.5) and (b) region-II (0.51 ≤ α ≤ 0.7). |
It is shown from Table 5 that in region-I, the average E-values are different for different reaction mechanisms characterized by different functional relationship of g(α), under all the four different heating rates. For random nucleation (A2, A3, and A4) and contraction geometry (R2 and R3) mechanisms, the average E-values ranged from 14.45 to 38.56 kJ mol−1 and 80.94 to 82.85 kJ mol−1 respectively. Further, for the reaction order model (F0, F1, F2, F3, F4, F5, and F5.5), the average E-values of four heating rates ranged from 75.43 to 154.04 kJ mol−1, whereas in the diffusion model (D1, D2, D3, and D4) these values ranged from 160.52 to 175.36 kJ mol−1 which were closest to the average E-values obtained from DAE (157.12 kJ mol−1) and Starink (157.97 kJ mol−1) methods. Thus it can be safely concluded that diffusion appears to be the rate-determining step for the pyrolysis reaction of MIW at the initial stage (0.05 ≤ α ≤ 0.5), denoted by region-I.
Diffusion based models assume that the reaction rate is higher than reaction front propagation throughout the pyrolysis medium. Accordingly, the pyrolytic degradation reaction occurs at the boundary of the two phases like solid and gas. Therefore, the gas must diffuse into the formed product layer around the sample to reach to the unreacted solid for a reaction to progress and thus makes diffusion the rate-determining step.
Similarly, for region-II (0.51 ≤ α ≤ 0.7) the average E-values for different reaction order models ranged from 12.81 for zeroth order (F0) to 153.98 kJ mol−1 for five and a half order reaction (F5.5) model. However, it can be seen that the E-values obtained from DAE (=151.51 kJ mol−1) and Starink methods (=152.45 kJ mol−1) were close to the five and half order reaction model indicating that region-II is a five and a half order reaction controlled region.
Therefore, the mechanism for the pyrolysis reaction of MIW was one-dimensional diffusion model (D1; region-I: 0.05 ≤ α ≤ 0.5) followed by five and a half order (F5.5; region-II: 0.51 ≤ α ≤ 0.7) reaction, which is proportional to the concentration of reactant(s) in a particular reaction. The reaction mechanism obtained in present work for pyrolysis of MIW is consistent with the earlier research30,42,43 (i.e. for Dalbergia sissoo wood30 it was governed by D3 followed by F3 controlled, for waste tea42 D3 followed by F2, and for four different wood species43 D1,2,3 followed by F1).
α | DAE method | Starink method | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Eα, kJ mol−1 | Aα, min−1 | k, min−1 | ΔHα, kJ mol−1 | Energy barrier, kJ mol−1 | ΔGα, kJ mol−1 | ΔSα, J mol−1 K−1 | Eα, kJ mol−1 | Aα, min−1 | k, min−1 | ΔHα, kJ mol−1 | Energy barrier, kJ mol−1 | ΔGα, kJ mol−1 | ΔSα, J mol−1 K−1 | |
0.05 | 149.54 | 8.8 × 1012 | 0.05 | 144.98 | 4.56 | 151.43 | −10.57 | 150.33 | 1.0 × 1013 | 0.05 | 145.77 | 4.56 | 151.40 | −9.23 |
0.1 | 143.03 | 2.3 × 1012 | 0.15 | 138.32 | 4.71 | 151.66 | −21.83 | 143.84 | 2.8 × 1012 | 0.15 | 139.13 | 4.71 | 151.63 | −20.45 |
0.15 | 145.43 | 3.8 × 1012 | 0.25 | 140.64 | 4.79 | 151.57 | −17.90 | 146.25 | 4.5 × 1012 | 0.25 | 141.46 | 4.79 | 151.54 | −16.51 |
0.2 | 144.20 | 3.0 × 1012 | 0.38 | 139.34 | 4.86 | 151.62 | −20.09 | 145.04 | 3.5 × 1012 | 0.38 | 140.18 | 4.86 | 151.59 | −18.66 |
0.25 | 153.03 | 1.8 × 1013 | 0.56 | 148.11 | 4.92 | 151.31 | −5.25 | 153.88 | 2.1 × 1013 | 0.56 | 148.96 | 4.92 | 151.29 | −3.81 |
0.3 | 160.68 | 8.5 × 1013 | 0.83 | 155.70 | 4.98 | 151.07 | 7.58 | 161.53 | 1.0 × 1014 | 0.83 | 156.55 | 4.98 | 151.04 | 9.02 |
0.35 | 168.79 | 4.4 × 1014 | 1.26 | 163.75 | 5.04 | 150.82 | 21.17 | 169.65 | 5.2 × 1014 | 1.27 | 164.61 | 5.04 | 150.79 | 22.62 |
0.4 | 175.59 | 1.7 × 1015 | 1.72 | 170.51 | 5.08 | 150.62 | 32.56 | 176.46 | 2.1 × 1015 | 1.73 | 171.38 | 5.08 | 150.59 | 34.02 |
0.45 | 170.17 | 5.8 × 1014 | 2.28 | 165.04 | 5.13 | 150.77 | 23.35 | 171.05 | 7.0 × 1014 | 2.30 | 165.92 | 5.13 | 150.75 | 24.83 |
0.5 | 160.78 | 8.7 × 1013 | 2.88 | 155.60 | 5.18 | 151.06 | 7.42 | 161.68 | 1.0 × 1014 | 2.91 | 156.50 | 5.18 | 151.03 | 8.94 |
0.55 | 151.72 | 1.4 × 1013 | 3.96 | 146.47 | 5.25 | 151.36 | −8.01 | 152.64 | 1.7 × 1013 | 4.01 | 147.39 | 5.25 | 151.33 | −6.45 |
0.6 | 144.99 | 3.5 × 1012 | 6.05 | 139.64 | 5.35 | 151.59 | −19.56 | 145.93 | 4.2 × 1012 | 6.15 | 140.58 | 5.35 | 151.56 | −17.97 |
0.65 | 143.03 | 2.3 × 1012 | 9.95 | 137.57 | 5.46 | 151.66 | −23.06 | 144.00 | 2.8 × 1012 | 10.15 | 138.54 | 5.46 | 151.63 | −21.42 |
0.7 | 157.04 | 4.0 × 1013 | 24.32 | 151.46 | 5.58 | 151.18 | 0.45 | 158.02 | 4.9 × 1013 | 24.90 | 152.44 | 5.58 | 151.15 | 2.11 |
The frequency factor (Aα) for each value of a were computed using eqn (8) and found to vary from 2.3 × 1012 min−1 to 1.7 × 1015 min−1 for DAE and 2.8 × 1012 min−1 to 2.1 × 1015 min−1 for Starink methods for α value varying from 0.1 (around 293 °C) to 0.4 (around 338 °C). It is a fact that if Aα < 109 min−1, it is considered as low and often indicates less surface reactivity of the system.44 In present case, higher Aα values were noticed for 0.3 ≤ α ≤ 0.5 (from 8.5 × 1013 min−1 to 2.1 × 1015 min−1), which corresponds to the temperature range of 326–350 °C, indicating high surface reactivity of MIW particle along with formation of a simple activated complex.45 From Fig. 5a and Table 6, the value of Aα was relatively lower at hemicellulose shoulder (at α = 0.25, T = 320 °C, Aα = 1.8 × 1013 or 2.1 × 1013 min−1) than the main peak of cellulose (at α = 0.4, T = 338 °C, Aα = 1.7 × 1015 or 2.1 × 1015 min−1) for the pyrolysis of MIW, indicating a comparative lower rate of molecular collisions for the hemicellulose component.
Fig. 5b shows that all values of ΔHα are positive. This indicates the endothermic nature of the pyrolysis reaction. The highest and lowest values of ΔHα recorded were 171.38 kJ mol−1 at α = 0.4 (338 °C) and 137.57 kJ mol−1 at α = 0.65 (384 °C), respectively. The highest and lowest positive value of ΔHα attributes to the maximum and minimum dissociation of MIW during its pyrolysis at 338 °C and 384 °C, respectively and it is also evident from DTG thermogram of MIW (Fig. 1b). The ΔHα value increased for α from 0.2 to 0.4 and 0.65 to 0.7 indicating an increase in endothermicity of MIW pyrolysis in these ranges. However, for other ranges of α such as 0.05 to 0.1 and 0.4 to 0.65, it gradually decreased indicating a comparative decrease in endothermicity of reaction in these ranges. For α ranging from 0.05 to 0.7, the average ΔH values determined by DAE and Starink methods were 149.79 and 150.67 kJ mol−1. For the present case, values of ΔHα and activation energy (Eα) values followed the same trend. This trend is similar to the trends observed for pyrolysis process of lignocellulosic biomass, such as castor residue46 and garlic stalk.47 The difference between ΔHα and Eα corresponds to the energy barriers (EB) for the formation of activated complex for the reactants.48 A reaction having a higher energy barrier requires a higher temperature to be conducted. In the present case, the energy barrier for reaction increases (from 4.56 to 5.58 kJ mol−1) as evident from Fig. 5c when conversion increases from 0.05 to 0.7. During this, the temperature also increases from 276 to 398 °C and thus helps to overcome the energy barrier of required for the reactions to take place. The small value of EB (≈5.06 kJ mol−1) during complete reaction period indicates that small amount of energy is required to form the activated complex, indicating that the product formation would be easier to achieve and thus bioenergy production through pyrolysis of MIW will be easier.
Further, the degree of stability of a system may be estimated by the value of change in Gibbs free energy (ΔG). The higher values indicate the harder reactivity and thus more energy will be absorbed by the system during the whole reaction process, and lower values indicate that the product can be produced with lower energy supply because of its reactive nature. Fig. 5d depicts the change in Gibbs free energy (ΔGα) computed by DAE and Starink methods. In the present case, all ΔGα values are positive indicating a non-spontaneous behavior of the pyrolytic process which consumes a considerable amount of energy to occur. The ΔGα values varied with the conversion rate (0.05 ≤ α ≤ 0.7), indicating that the system passed through ups and downs of the degree of difficulty [Fig. 5d; 0.05–0.1 (up), 0.1–0.15 (down), 0.15–0.2 (up), 0.2–0.4 (down), 0.4–0.65 (up), 0.65–0.7 (down)] for the MIW pyrolysis. The average values of ΔG were computed as 151.27 and 151.24 kJ mol−1 for DAE and Starink methods.
The entropy, being a state function, is an index to indicate degree of randomness or disorderness associated with the reaction system. The change of reaction entropy (ΔSα) reflects how near the system is to its own thermodynamic equilibrium. The lower value of entropy (ΔSα) signifies that the reactive system is approaching to its own thermodynamic equilibrium state. In this case, the substance shows low reactivity and thus, takes more reaction time to form the activated complex.49 On the other hand, a large value of ΔSα means the system is far from its thermodynamic equilibrium. In this case, the reactivity of material is high, and the reaction system will react faster to produce activated complex, and consequently, takes short reaction time to form the activated complex. For MIW pyrolysis, the ΔSα values were found both positive as well as negative. The ΔSα values were negative for each value of conversion (α) except α ranging from 0.3 to 0.5 (Fig. 5e).
Pyrolysis is a complex process and is governed by a set of parallel, complex and competitive reactions taking place simultaneously which involve the disintegration of natural polymers (hemicellulose, cellulose, and lignin) within the biomass material as well as integration of smaller molecules to form bigger molecules with the increase in temperature. The entropy of the system increases when bulky polymer molecules break to form a large number of volatiles/gas molecules and decrease when the smaller molecules combine to create a larger molecule during the pyrolysis process. The increment in ΔSα values for conversion from 0.2 to 0.4 (311 to 338 °C) was largely due to the formation of a large number of molecules, owing to the simultaneous dissociation of hemicellulose, cellulose and lignin molecules.50 The increment in ΔSα values for above conversion range shows increase in reactivity which is supported by the fact that in this range ΔGα values decrease monotonously indicating decrease in endothermicity of reaction. For, above conversion range, the increase in activation energy (Eα) and ΔHα values show an increase in EB values as well as endothermicity of the reactive system. It is indicative of the fact that the system absorbs large amount of energy for dissociation of larger molecules in this range for which the entropy of the system has increased. After the above conversion range, the system entropy decreases up to a conversion value of 0.65 (384 °C), which is indicative of the association of small molecules to form product molecules within the reactive system.44 This is consistent with the decrease in ΔHα values that show a decrement in endothermicity of the reactive system. This supports the fact that in this region the entropy of the system decreases due to the dominance of reaction in which association of molecules takes place. In this range, EB slightly increases which is taken care of by the rising temperature.
For the conversion range (0.05 ≤ α ≤ 0.7), the maximum negative value of ΔSα was −23.06 J mol−1 K−1 for DAE method at α = 0.65 and maximum positive value of 34.02 J mol−1 K−1 for Starink method at a = 0.4 (Table 6). This indicates that the system was in the highest state of stability at conversion value 0.65 (T = 384 °C) and most reactive around conversion value of 0.4 (T = 338 °C). It is interesting to see that the value of ΔGα around conversion value 0.65 is also highest (Fig. 5d and Table 6) indicating that the endothermicity is very high and thus conducting reaction is difficult. Also around α = 0.4, the Aα vs. α plot (Fig. 5a) showed the highest peaks, confirming that the system was most reactive around the conversion value at 0.4 which is also confirmed by DTG plot shown in Fig. 1b. For the rest of the conversion points, the negative value of ΔSα reflects that activated complex needs comparatively more thermal energy to propagate the reaction in the forward direction or we can say that more or less the thermal equilibrium has been reached and a thermally stable product has been produced.51 This statement are in accordance of the lower values of Aα as demonstrated in Aα vs. α plot (Fig. 5a), confirming the less reactivity of the system in the conversion range from 0.05 ≤ α < 0.3 and further for 0.5 < α ≤ 0.7.
The reaction rate constant (k) was computed using the Arrhenius correlation and found to be continuously increasing with conversion value up to 0.7 (Fig. 5f). The lowest and highest value of the rate constant for the pyrolysis of MIW was estimated to be ranged from 0.05 to 24.90 min−1. It was expected as the temperature continuously increased in the range 0.05 ≤ α ≤ 0.7 and reaction rate constant is a strong function of temperature.
1. The characterization for as-received MIW revealed high volatile content, fixed carbon, and higher heating value as 73.38 mass%, 15.09 mass%, and 14.89 MJ kg−1, respectively. Further, no sulfur was detected in MIW sawdust.
2. There was no major discrepancy in apparent activation energy values computed by DAE (154.86 kJ mol−1) and Starink (155.74 kJ mol−1) methods.
3. The reaction mechanism for the pyrolysis of MIW can be classified as one-dimensional diffusion type (D1; g(α) = α2) for 5–50% conversion followed by five and a half order (F5.5; g(α) = [1 − 1/(1 − α)4.5]/−4.5) heterogeneous rate expression up to the rest conversion value of 70%.
4. The average values of A, ΔH and ΔG obtained were 2.2 × 1014 min−1, 149.79 kJ mol−1 and 151.27 kJ mol−1 for DAE method and 2.6 × 1014 min−1, 150.67 kJ mol−1 and 151.24 kJ mol−1 for Starink method, respectively.
α | Conversion or reaction progress, dimensionless |
β | Heating rate, °C min−1 |
Ea | Apparent activation energy, kJ mol−1 |
Eα | Activation energy at conversion α, kJ mol−1 |
A | Frequency factor, min−1 |
k | Reaction rate constant, min−1 |
n | Order of reaction, dimensionless |
R | Universal gas constant, kJ mol−1 K−1 |
T | Temperature, K in all equations and °C elsewhere |
p(x) | Exponential/temperature integral approximation |
f(α) | Differential form of reaction model |
g(α) | Integral form of reaction model |
R2 | Correlation coefficient, dimensionless |
MIW | Mangifera indica wood |
TG | Thermogravimetric analysis/thermogravimetry |
DTG | Derivative thermogravimetry |
DAE | Distributed activation energy |
C–R | Coats–Redfern |
ASTM | American Society for Testing and Materials |
HHV | Higher heating value or gross calorific value |
LHV | Lower heating value or net calorific value |
EB | Energy barrier |
This journal is © The Royal Society of Chemistry 2021 |