Abstract

The combustion kinetic characteristics of wood powder and pellet were investigated within thermogravimetric analyser (TGA) and tube furnace system. The kinetic parameters of these two different forms of woody fuel were measured and derived by double-step-and-double-equal and isothermal method, respectively. The results showed that the combustion mechanisms of wood powder kept consistent through the whole process, while the combustion mechanisms of wood pellet differed significantly between the volatile and char combustion stages. The most probable mechanism functions of the two different forms of woody fuel were not the same due to the differences in internal heat and mass transfer properties. In addition, activation energy values varied from 92.33 kJ·mol−1 for wood powder to 71.20 kJ·mol−1 for wood pellet, while the preexponential factor value of wood powder (2.55×108 s−1) was far greater than the one of the wood pellet (78.55 s−1) by seven orders of magnitude.

1. Introduction

Suffering the consequences of excessive use of fossil fuel, the countries all over the world contribute more and more to renewable energy development and utilization [13]. Biomass energy is the most potential renewable energy in near future, because it is easy to store and convenient to transport and its application technology is ordinary and mature. Furthermore, from production to utilization, the cycle chain of CO2 emissions from biomass fuel closes to zero [4, 5]. The annual use of wood, widely in construction, furniture manufacturing, packaging, and handicraft, is more than 400 million m3, yet the processing rate is only about 65% [6]. The remaining compounds such as leftover sawdust should be utilized vigorously to promote energy conservation and emission reduction [7, 8]. The most effective way to use the remains is to extrude them into granules with uniform shape and density and then directly burn them for heating and power generation [911].

The previous researches of biomass granular fuel mainly focus on the macro combustion characteristics. Ge et al. [12] calculated the comprehensive combustion characteristics of four kinds of the man-made board in air atmosphere by the thermogravimetric curve and compared the characteristics of synthetic combustion performances. Cuiping et al. [13] piecewise compared the combustion characteristics of three types of granular fuels, namely, corn straw, wood sawdust, and mixed sawdust, through the apparent activation energy values and the preexponential factor values under different combustion stages. Qing et al. [14] calculated the activation energy values and the preexponential factor values of four kinds of biomass materials by the Coats-Redfern method and the first order model, and they came to a conclusion that the wood sawdust is less active than the remaining three species of herbaceous biomass. The calculation compared well with the experimental results. In addition, some scholars explored the granular fuel pyrolysis and combustion rules by mathematical simulation, such as Haseli [15], establishing a mathematical model to describe the combustion process of a single biomass particle, based on a one-dimensional model of quality, momentum, and energy conservation during combustion, with the combination of combustion, pyrolysis kinetic equations.

Owing to the advantages of quick transmission of mass and heat and to acquire more accurate kinetic features of the experimental materials, the samples in kinetic experiments were mostly in fine particles (diameter <0.20 mm) and in small amount (<20 mg) [1620]. However, the internal heat and mass transfer should not be ignored because of the objective existence while researching on the combustion kinetics of the whole granule [17, 21]. For the best of our knowledge, there was no combustion kinetic research of wood pellet based on the isothermal method in previous studies. To further explore the differences in the combustion kinetic characteristics between wood powder and pellet, some thermal kinetic analyses were conducted in the respective status of the same kind of wood. The results are expected to be the supplement to the existing research of wood pellet combustion kinetics.

2. Materials and Methods

2.1. Materials

The experimental materials were wood-based granular fuel produced by a contract-energy-management enterprise in Foshan, Guangdong. The sawdust and residues produced during furniture manufacturing were crushed and squeezed into densified cylindrical granules with approximately 8 mm in diameter, 40-60 mm in length, and about 950 kg/m3 in density, including eucalyptus, fir, beech, and walnut, collectively referred to as miscellaneous wood. The samples could be represented by sampling in dispersive points from the same product batch. According to the national standard DL/T 568 and GB/T 28731, the proximate analysis and ultimate analysis results are listed in Table 1. Before the experiment, some pellets were cut to both ends of the flush, with the mass of 1000 ± 5 mg, and the others were ground and sieved to powders with the diameter less than 0.18 mm. All samples were stored in a desiccator standby.

2.2. Apparatus and Instrument
2.2.1. Wood Powder Analysis

The experiments were conducted in an STA-449F3 thermogravimetric analyser (NETZSCH, Germany), with powder samples weighing approximately 10.0 mg, using air as the atmosphere gas at the flow rate of 60 mL/min. The reaction temperature rose to 1000°C from the environmental temperature in heating rates of 10, 20, 30, and 40°C/min, respectively.

2.2.2. Wood Pellet Analysis

The experiments of wood pellet were conducted in a self-made tube furnace system. The main part of the device was a tube reactor, typed SK-4-4-16Q, with an inner diameter of 60 mm and mid-effective heating length of 300 mm, produced by Wuhan Yahua Electric Furnace Ltd. When the experiment began, the influx of air was fed into the tube furnace at the flow rate of 3 L/min. As the furnace atmosphere was stable and the temperature reached the set value, i.e., 800, 1000, 1200, and 1400°C, a corundum crucible carrying a pellet was put into the center of the reactor rapidly, and then the sample was ignited. Afterwards at the selected time moment, i.e., 12, 24, …, 60, 90, …, 210 s, the crucible was taken out and put into a vessel filled with argon immediately, and then the vessel was sealed tightly and put into a water sink to cool down quickly. After that, the sample which kept the morphology, composition, and properties was weighed, and finally a set of isothermal TG (Thermal Gravity Analysis) curves were drawn.

2.3. Method

The combustion process of wood obeys the Arrhenius law; the basic kinetic equation is generally described as follows:where is the conversion of the combustible part of the sample, which is described as follows:where is the mass of the sample; , , and are, respectively, the mass of the sample at initial time, moment, and finally time; is the combustion rate constant; is preexponential factor, s−1; is apparent activation energy, kJ/mol; is the ideal gas law constant, J/(mol·K); and is the reaction model.

2.3.1. Wood Powder Thermogravimetric Analysis

By the double-step-and-double-equal method for wood powder, concrete steps are as follows [22].

(a) Determination of Most Probable Mechanism Function G(α). Equation (3) is integrated and transformed into Flynn-Wall-Ozawa equation aswhere is the heating rate, obtained from the formula as follows:

transforming (3) into

By substituting the at the same temperature on the TG curves corresponding to different heating rates and different mechanism function in (5), a series of correlation coefficients can be achieved with lg G(α) -lg β linear regression. The is closer to 1, the corresponding G(α) describes the wood powder combustion process better, and, hence, the closest one is the most probable mechanism function.

(b) Derivation of E and A. According to (3), a plot of lg β-1/T corresponds to a straight line with the same at different heating rates and the determinate form of , from which the values of and can be derived from the slope and intercept, respectively.

2.3.2. Wood Pellet Thermogravimetric Analysis

(a) Determination of Most Probable Mechanism Function G(). The most probable mechanism function of granular fuel combustion can be deduced from the coordination of the experimental data and the kinetic models, namely, reduced time method [22]. Nine commonly used mechanism functions (Table 2) were tested in this paper. If an experimental curve overlaps a theoretical curve, or the experimental data points almost fall in one theoretical curve by comparing the square of deviance [22], the corresponding to the theoretical curve is the most probable mechanism function, through where t0.5 is the time length while α=0.5; G(α)α=0.5 is the G(α) while α=0.5; and G(α)α=1.0 is the G(α) while α=1.0.

(b) Derivation of E and A. By transforming the determinate most probable mechanism function to plot a straight line including and , and can be achieved according to the slope and intercept. And then substituting the and corresponding to various temperature sections into Arrhenius equation, one has

Plotting the straight line of against , , and can be derived from the slope and intercept, respectively.

3. Results and Discussion

3.1. The Solution of Combustion Kinetic Parameters of Wood Powder

The TG and DTG (Differential Thermal Gravity) curves of wood powder combustion are shown in Figures 1 and 2, respectively.

By substituting the at the same temperature in four TG curves, heating rates, and different mechanism functions in (5), with - linear regression, the fitting results were achieved (shown in Table 3). It was clear that the closest to 1 corresponds to the function G(α)=(1-α)n (n=-1), namely, the most probable mechanism function to describe the combustion process of wood powder.

And then, by substituting the and corresponding to the same at four TG curves of different heating rates and the determinate in (3), with linear regression of -1/T, a series of values were derived from the slopes and values from the intercepts, and the results are shown in Table 4.

The fitting results of values shown in Table 4 were reliable for the correlation coefficient values varying from 0.9030 to 0.9991. The results of values shown in Table 4 were solved according to the values and the determinate . Global average and was 92.33 kJ mol−1 and 2.55×108 s−1, respectively.

3.2. The Solution of Combustion Kinetic Parameters of Wood Pellet
3.2.1. The Determination of Most Probable Mechanism Function G(α)

The weight loss curves of wood pellet during isothermal combustion are shown in Figure 3. The fitting results of α-t/t0.5 theory and experimental curves are shown in Figure 4. It was clear that the experimental data fell on the theoretical curve of function G(α)=A2(α) or nearby when the conversion rate of varied in the range of [0.10,0.80).

(a) The Derivation of E and A. The most probable mechanism function was transformed from

transforming (8) into

By plotting the straight line of ln[-ln(1-α)]-ln t as shown in Figure 5, the n and k could be derived from the slope and intercept, respectively.

Substitute the ki and Ti of the four temperature sections into Arrhenius equation as

Plotting the straight line of against as shown in Figure 6, and were derived from the slope and the intercept to be 71.20 kJ·mol−1 and 78.55 s−1, respectively, correlation coefficient to be 0.9844.

3.2.2. Comparison of Combustion Kinetic Characteristics between Wood Powder and Pellet

Wood pellet combustion is a complex reaction of multiple materials, multiple phases, and multiple steps [15, 24, 25]; thus the kinetics triplet, which might be the apparent result of multistep reactions, usually varies with the combustion condition and the reaction progress. Designed towards precision and automation, the TGA (thermogravimetric analyser) currently is most suitable for thermal measurement on small amounts and small size of the sample. Hence, there is a certain distortion of wood pellet thermal measurement due to the intense internal heat and mass transfer. The comparison between wood powder and pellet combustion kinetic parameters is shown in Table 5.

From Table 5, it is clear that the mechanism functions of wood powder and pellet are different. The discrepancies indicated that sample size, shape, amount, and heating rate could lead to the difference in combustion mechanism of the same kind of fuel [19]. For wood pellet, when α<0.80, the combustion process could be described by the mechanism function G(α)=[-ln(1-α)]1/2, but when α≥0.80, no function that had been selected could fit the corresponding combustion process, whereas for wood powder, there was no piecewise feature; the whole combustion process could be described by the mechanism function G(α)=(1-α)]−1. Referring to the proximate analysis result shown in Table 1, was the dividing point where volatile combustion transits to charcoal combustion, illustrating the fact that the complicated heat and mass transfer effects on combustion of the wood pellet were much greater than those on wood powder due to the complexity of internal structure. Therefore, origin pellet should be used as far as possible in the research targeting kinetic parameters determination.

The values of apparent activation energy , characterizing the molecular activity of wood pellet and powder, were of the same order with the difference 21.13 kJ·mol−1, indicating the consistency of the two forms of materials. As for the characterization of combustible molecule and oxygen molecule collision frequency [26], the preexponential factor values of the two materials differed greatly, with the latter being seven orders of magnitude higher than the former, showing that the wood powder was much more probable to collide with O2 than wood pellet in the combustion process. All the calculations were consistent with the apparent phenomenon. In engineering, the former was mainly applied to the layered combustion, and hence good ventilation should be ensured to make the fuel surface fully contacted with the fresh air [27, 28], while the latter was often used in suspension combustion; despite the relatively static state between the fuel and air [29], the fuel still showed good burnout performance for the high molecule collision frequency.

4. Conclusions

By the kinetic analysis and results comparison of the combustion process of wood powder and pellet of the same kind of wood, the following conclusions are drawn:

(1) For wood powder, the combustion mechanism of the whole combustion process was basically unchanged, and the most probable mechanism function had no obvious segmental characteristics; for wood pellet, there was a great discrepancy between volatile combustion stage and char combustion stage, respectively, described by different most probable mechanism functions.

(2) The activation energy values were kJ mol−1 for wood powder and kJ mol−1 for wood pellet, respectively, with a difference less than 30%. The preexponential factor value of the former (2.55 x 108 s−1) was far greater than that of the latter (78.55 s−1), by seven orders of magnitude.

(3) Currently, most of the TGA is only suitable for combustion kinetic analysis in wood powder with less sample quantity and small particle size; hence it is difficult to truly reveal the combustion kinetic characteristics of bulk granular biomass fuel. In academic researches and engineering applications for the granular fuel, the TGA that is suitable for bulk fuel should be adopted to reduce the experimental errors induced by the sample.

Data Availability

The datasets generated and analysed during this study are available in the FAIRsharing repository: https://figshare.com/s/d51ecef2f1992ae0b7c3.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors thank Shunchun Yao and Saihua Jiang for comments and insights leading to the development of this manuscript. The authors gratefully acknowledge support by South China University of Technology and Guangdong Institute of Special Equipment Inspection and Research, Shunde Branch. This work was supported by Administration of Quality and Technology Supervision of Guangdong Province Research Funds (Grant no. 2017PB08).