Next Article in Journal
Drying Hot Red Chilies: A Comparative Study of Solar-Gas-Fired, Tunnel, and Conventional Dryers
Previous Article in Journal
Challenges Associated with the Production of Nanofibers
Previous Article in Special Issue
Energy Recovery Efficiency of Integrating Anaerobic Co-Digestion of Pig Slurry and Feedlot Cattle Manure and Hydrothermal Carbonization of Anaerobic Sludge Cake
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Systematic Approach to Determining the Kinetics of the Combustion of Biomass Char in a Fluidised Bed Reactor

Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, UK
*
Author to whom correspondence should be addressed.
Processes 2024, 12(10), 2103; https://doi.org/10.3390/pr12102103
Submission received: 5 August 2024 / Revised: 15 September 2024 / Accepted: 23 September 2024 / Published: 27 September 2024
(This article belongs to the Special Issue Biomass Pretreatment for Thermochemical Conversion)

Abstract

:
The aim of this work was to investigate the combustion of biochar in a fluidised bed and determine the intrinsic kinetic parameters for combustion: pre-exponential constant A i and activation energy E i . When analysing the rates of reaction, Regimes I, II and III were demonstrated, with values for the activation energy of 155, 57 and 9 kJ/mol, respectively, when combustion was limited by different factors: intrinsic kinetics, intraparticle and external mass transport phenomena. These mass transport phenomena were decoupled from a set of ‘apparent’ kinetics incorporating effectiveness factors, which we used as a starting point in the determination of the intrinsic kinetic parameters. We also investigated a simple approach to model the evolution of the char structure over the course of oxidation using an empirical function, f X , fitted with an O(7) polynomial. We then reassessed the division into three combustion regimes by exploring the changes in f X and the intraparticle effectiveness factor that occurred upon increasing the combustion temperature. Overall, we demonstrate that experiments in a fluidised bed can be used to determine biochar kinetics in a simplified but trustworthy way.

1. Introduction

Biomass is the most accessible and affordable fuel worldwide. The primary benefit of biomass compared to fossil-derived fuels is the potential to achieve net-negative CO2 emission when utilised for bioenergy with carbon capture and storage (BECCS) [1]. Among the available technologies, fluidised bed combustors are particularly well fitted for the thermal conversion of biomass because of their operational flexibility. They offer fast heat transfer, resulting in a uniform temperature in the bed, while long residence times ensure the high conversion of solid fuels [2]. Combustion in fluidised beds has been already studied at length, with the majority of work focused on fossil fuels rather than biomass.
A key parameter in designing a combustor is information on the kinetics of combustion, especially useful when aiming for efficient operation at relatively low operating temperatures, e.g., in heating applications. To determine the kinetics of combustion, various methods have been presented in the literature. Conveniently, thermogravimetric analysers (TGA) inform us about the rates of mass loss and provide means for quick experimentation; thus, they are often employed in combustion research. However, the measured reactivity of chars in TGA differs from that in fluidised beds, arising from differences in gas delivery, heating rates and possible char restructuring [3]. Difficult to eliminate are the effects from various transport phenomena; thus, the results reported in the literature predominantly provide the apparent kinetics of char combustion without isolating the intrinsic reaction kinetics. To assess the true kinetics, careful experimental design is required to ensure purely kinetic limitations. For example, Morin et al. [4] investigated beech combustion at low temperatures (603 to 643 K), yet mass transfer remained important because of the large particles of char that were used (cylinders of 4 mm diameter and 9 mm length).
As recently reviewed by Kwong, reaching a set of trustworthy kinetic parameters for a very reactive fuel, such as biochar, is not a trivial task [1]. This paper aims to present a systematic approach to determining the intrinsic kinetics of the combustion of biomass char directly from fluidised beds.
To include the effects of different transport phenomena in a fluidised bed, the apparent rate of reaction—defined here as the observed rate of char consumption per unit mass of the fed char and thus in g s−1—can be expressed as
r = d X d t = η i p η e η p k i P O 2 n f X
where X   and P O 2 are the char conversion and partial pressure of O2, respectively, while f X is a function connecting the evolution of the pore structure to the char conversion. Then, η i p , η e , η p are effectiveness factors that account for interphase, external and intraparticle mass transfer phenomena, respectively. The interphase mass transfer refers to the transport of O2 from gas bubbles to the particulate phase in the fluidised bed. The mass transfer of O2 from the bulk of the particulate phase to the reactive surface of the char particle is referred to here as external mass transfer, whereas intraparticle mass transfer concerns the transport of O2 from the external surface of the char particle through its pores to the site on the particle surface where combustion occurs.
The parameter k i in Equation (1) represents the reaction constant, describing the intrinsic kinetics of combustion—here, in bar−n s−1. Here, the dependence of the apparent rate, r , on the partial pressure of O2 is assumed in the form of a power law, where n in Equation (1) is the reaction order with respect to oxygen—a similar approach to Morin et al. [4]. The evolution of the reactive surface of the char particle during combustion is represented by the function f X , which depends on the conversion of the particle, X .
For combustion mechanisms in different operating conditions and regimes, we provided an overview in our recent review paper [1]. Multiple mechanistic models were revised by Hurt and Calo for coal chars. McMurchie et al. from our group looked into the combustion of char from oak wood, assessing the order of reaction n to infer about the global mechanism and the rate-determining elementary steps [5]. A good level of discussion on the Langmuir–Hinshelwood approach to explain combustion was also provided by Morin et al. [4]. The mechanism of combustion by itself is a substantial and complex topic. In contrast, our goal in this paper is to focus on the intricacies of determining global kinetic parameters of biochar combustion in fluidised beds, presenting a step-by-step approach to decoupling the effect of transport phenomena from the experimentally observed rates. This aims to provide guidance for the extraction of information about the intrinsic kinetics of combustion directly from fluidised bed experiments.

2. Materials and Methods

Particles of biochar, created from mixed sources of hardwood, were purchased from Soilfixer (UK). Granules of biochar (<8 mm diameter) were crushed and sieved to different ranges, e.g., 200–250, 250–300, 300–355, 355–422, 422–500 and 500–600 µm. The fuel was used in our previous study [6], where the results of the proximate and ultimate analysis, BET area and BJH pore size are presented. Silica sand (David Ball Group Ltd., Cambridge, UK, Fraction D) was sieved to 200–250 µm, washed in deionised water and dried for 48 h at 393 K before being used as the bed material in the fluidised bed reactor.
The combustion experiments were performed in a bubbling fluidised bed housed in a quartz reactor (i.d. 30 mm), with a porous disk (grade 1) located 110 mm above the gas inlet. The disk’s function was to distribute the gas and support the bed material. The quartz reactor was filled with ~30 g of silica sand and was heated electrically by a tubular furnace. The temperature was raised to the desired setpoint as measured by a K-type thermocouple, which was inserted into the reactor with the tip located ~1 cm above the porous disk. The outlet gas from the reactor was sampled continuously with a pump (1 L min−1 at NTP) and passed through a drying tube of CaCl2, and then directed to an ABB EL3020 analyser equipped with an NDIR cell to measure the concentrations of CO and CO2. A lid was placed on top of the quartz tube to prevent air ingress to the reactor, as also described elsewhere [6,7].
Gases with different concentrations of O2 (5, 11, 21 vol% O2) in N2 were used to fluidise the bed with U / U m f ~3, where the superficial, minimum fluidisation velocity, U m f , was determined experimentally at each experimental temperature. To study char combustion, batches of biochar—3 to 40 mg—were dropped into the reactor at a temperature between 723 and 1123 K. From preliminary experiments, a minimum temperature of 673 K and mass of 3 mg were required to obtain CO and CO2 profiles that were distinguishable from the background noise. The offset of combustion, i.e., t = 0 , was defined as the moment when the measured concentration of CO2 increased above 0.01 vol%, while the endpoint, i.e., t = t b u r n o u t , was achieved once the CO2 reading decreased below 0.01 vol%.
A summary of the conditions used in the experimental campaigns to investigate different transport phenomena is presented in Table 1. Every experiment was repeated at least once, and only results with a carbon balance above 80% were analysed for kinetics. Consequently, for the experiments at 5 vol% O2 and lower temperatures, only a reduced number of data were obtained, further contributing to the simplified assessment of Regime I for this O2 concentration—see Section 3.3.
Assuming that carbon combusts only to CO and CO2, the mass of the carbon burnt ( m C ) during an experiment and the conversion at a given time, X t , could be calculated from the molar flowrates of the produced CO and CO2, n ˙ C O , o u t and n ˙ C O 2 , o u t :
m C = M C t = 0 t b u r n o u t n ˙ CO , o u t + n ˙ CO 2 , o u t   d t
X t = M C t = 0 t n ˙ CO , o u t + n ˙ CO 2 , o u t   d t m C
where M C is the molecular mass of carbon (12.01 g mol−1) and t is an arbitrary duration of time, aligned with the time of char combustion and the time of measurement of CO and CO2 traces ( t t b u r n o u t ) .
The instantaneous rate of reaction, as inferred by the rate of consumption of carbon for a unit mass of biochar, d X / d t in s−1, is determined as
r = d X d t = M C n ˙ CO , o u t + n ˙ CO 2 , o u t m C

3. Results

3.1. Influence of Interphase Mass Transfer (Campaign 1)

Fluidised beds can be considered as systems of two interacting phases—(1) a particulate phase, with solid particles fluidising in a portion of the provided gas, and (2) a gas phase, where the rest of the provided gas travels through the particulate phase as bubbles. Solid particles of char combust in the particulate phase, consuming, in the first instance, oxygen from the interstitial gas in this phase. The depletion of O2 in the particulate phase drives interphase diffusion. If the mass transfer between the two phases is slower than the rate of reaction in the particulate phase, the transport of gaseous O2 from bubbles to the mixture of particles and interstitial gas (particulate phase) limits the rate of combustion. Consequently, the char particles in the particulate phase are exposed to a different oxygen concentration than the O2 concentration at the entrance to the reactor. To determine the intrinsic kinetics, first, the observed rate should be decoupled from the interphase mass transfer. To demonstrate the significance of the interphase mass transfer, the mass loads of biochar used in the experiments were varied from 6 to 20 mg, and the results of the instantaneous rate of combustion in 11 vol% O2 at 823 K are presented in Figure 1a.
No change in the instantaneous rate of combustion was found when varying the mass of biochar, m b i o c h a r , indicating that the measured rate was independent of the introduced amount of carbon. Thus, we conclude that, as seen in Figure 1a, combustion was not limited by interphase mass transfer. In addition to the experiments carried out at 11 vol% O2 and 823 K, the instantaneous rate of reaction was also measured for 5 vol% and 873 K, comparing the values at a fixed conversion point of X = 0.1. As presented in Figure 1b, the rate changed very little with the input mass—the apparent trend of some decrease in rate is within the experimental error as seen when comparing the results at 5vol% O2, 873 K with those at 21vol% O2, 823 K. Consequently, for the conditions shown in Figure 1, i.e., temperatures of 823–873 K, pO2 and m b i o c h a r , the apparent rate of combustion in Equation (1) can be simplified by taking η ip = 1.
Learning from Campaign 1, experiments that aimed to determine the kinetics of combustion were performed using 10 mg of biochar (see Campaign 3 in Table 1), noting also the rest of the conditions (temperatures, pO2) which helped to mitigate the effect of the interphase mass transfer of O2.

3.2. Influence of External Mass Transfer (Campaign 2)

External mass transfer, the transport of O2 from the bulk concentration in the particulate phase to the reactive surface of the char particle, limits the rate of combustion when large char particles are used, especially severely at higher temperatures, when the intrinsic kinetics are rapid. Thus, first, combustion experiments at varied particle sizes were conducted. The results of the instantaneous rates of reaction at 773 K, fluidised by air, with varied sizes of biochar particles, dp, are presented in Figure 2. As seen in Figure 1a, the maximum rates of combustion of the studied biochar align with the early stage of combustion at conversion X below 0.2. Here, instantaneous rates at X = 0.1 are used for the analysis, representing the fastest stages of combustion but also the initial rates when the char diameter is still close to the nominal particle size. When the particle burns in a regime limited by O2 delivery to the external surface, its diameter decreases with conversion, as in a shrinking-core situation. Thus, taking X for the largest particle size in the burnout process corresponds to the “worst-case scenario”, i.e., a situation most severely affected by external mass transport.
For a single particle, the rate of external mass transfer is proportional to the external surface of the char particle, i.e., d p 2 . As the mass of the particle is proportional to d p 3 , the rate of combustion per unit mass would be inversely proportional to d p , if the process were solely limited by external mass transfer. However, as shown in Figure 2, the experimental rate of combustion was found to be independent of the particle diameter in the investigated region, indicating strongly that the effect of the external mass transfer had little influence on the rate of combustion.
To further investigate the importance of external mass transfer, we estimated the range of mass transfer rates expected in the presented experiments. The rate of mass transfer of O2 (molar flowrate) to a single particle of char is 4 π d p k g c B c S , where c B is the bulk concentration of O2 and c S is the concentration of O2 at the external surface of the char. Assuming a fast reaction and slow O2 delivery, the surface concentration, c S , can be taken as 0, i.e., O2 becoming fully depleted. The external mass transfer coefficient, k g , connects to the Sherwood number, S h = k g d p / D G , where D G is the molecular diffusivity of O2 in N2. The Sherwood number, S h , can be estimated from a suitable correlation—for a bubbling fluidised bed, this is provided by Scala [8]:
S h = 2 ε m f + 0.70 R e m f / ε m f 1 / 2 S c 1 / 3
where ε m f is the voidage of the particulate phase at the onset of fluidisation—here, taken as 0.51, as estimated by fitting the Carman–Kozeny equation to the experimental results. Then, R e m f is the Reynolds number based on U m f , and S c is the Schmidt number. The resulting rate of external mass transfer (per mass of biochar) to particles of char differing in size is presented in Figure 2.
As expected, the theoretical rate of mass transfer decreases with the particle size. In the considered conditions, i.e., pO2 = 0.21 bar and 773 K, the calculated rate of external mass transfer is at least an order of magnitude larger than the rates observed experimentally. Thus, we conclude that the external mass transport in our experiments was not a limiting factor. Consequently, the particle size of the fraction 250–300 μm was deemed suitable for further work in Campaign 3, where we aimed to study the intrinsic kinetics of combustion, as shown in Table 1. Deciding on a particle size and experimental conditions in which the effects of external mass transfer are mitigated enables the setting of η e = 1 in Equation (1).

3.3. Influence of Combustion Temperature (Campaign 3)

The combustion reaction of biochar could proceed under three possible reaction regimes depending on the temperature: Regime I, where combustion is limited by the intrinsic kinetics; Regime II, where combustion is influenced by both the kinetics and mass transfer; and Regime III, where only the mass transfer limits the combustion. Therefore, the intrinsic kinetics can only be extracted directly when the combustion is in Regime I. Alternatively, the assessment necessitates the decoupling of the kinetics from transport phenomena.
To identify the extent of the three regimes, batches of particles of biochar ( d p = 250–300 µm) were combusted, varying the combustion temperature. Results taken as logarithms of the instantaneous rate of reaction (per mass of biochar added) at conversion of 0.1 are plotted in Figure 3.
Three straight lines, each with different gradients because of the difference in the apparent activation energy, are depicted in Figure 3a,b. The decision about which points to classify under each regime was made by maximising the sum of the R2 (correlation coefficient) values for the three lines.
As shown in Figure 3a, the transition temperature between Regimes I and II is ~773 K and that between Regimes II and III is ~883 K during combustion in 21 vol% O2. In Regime II, the gradient is nearly half the gradient of Regime I, as expected from the theoretical findings in the literature [9]. The gradient of Regime III is close to zero.
The transition between different regimes was also observed in experiments carried out in 5 vol% O2, as shown in Figure 3b. Here, we made a deliberate choice to limit Regime I to the same temperature range as for 21 vol% O2. This then necessitated the description of Regime I with two temperature points only, as we did not have further results for temperatures between 750 and 798 K. Consequently, our analysis is a showcase and only a first approximation of Regime I at 5 vol% O2, as the fitted line—needed for parameter estimations—has no associated standard error.
The transition between Regimes II and III, judged solely on the changes in the gradient, occurred at a similar temperature for both the 5 vol% and 21 vol% O2 experiments. This is again an approximated result. With a lower concentration of O2, the driving force for external mass transfer is lower; hence, the reaction would have been more readily limited by external mass transfer in the 5 vol% case, with the transition to Regime III occurring at a lower temperature—not clearly observed in Figure 3b.
The kinetic constant of a chemical reaction can be expressed as a function of the temperature, following the Arrhenius equation, i.e., k i = A i exp E i R T . Substituting k i in Equation (1) leads to the dependency of the rate of combustion on the temperature. After taking a log on Equation (1) and asserting from the previous results that η e = η i p = 1 (i.e., negligible external and interphase mass transfer), Equation (1) now becomes
ln d X d t = E i R 1 T + ln η p A i P O 2 n f X
The gradient and the y-intercept of the lines of best fit in Figure 3a,b provide information on the apparent activation energy, E i , and pre-exponential factor, A i , from the transformation of Equation (1). From Equation (6), the activation energies at different reaction regimes are determined, as presented in Table 2. Under a kinetically controlled regime, the activation energy for 21 vol% O2 is slightly higher than in 5 vol%, but both remain within the standard error of the results at 21 vol% O2. Worth noticing is that the obtained values of Ei are similar to that for coal chars given by Smith, 179.4 kJ mol−1, and often used by the combustion community, also when studying biochar combustion [1].
As mentioned, the values of E i can be calculated directly from the slopes of the linear fits in Figure 3. In contrast, the ln term in Equation (6), which denotes the intercepts, demonstrates that the pre-exponential factors, A i , depend on η p (effectiveness factor for intraparticle mass transfer). As we discuss in Section 4.2, η p contains information about the reaction kinetics; thus, the ln term is rich in information that is highly coupled, including the coupling between A i and E i .   Additionally, because E i   is already present in the first term of Equation (6), ( E i / R T ) , the two terms (and, in the linear fits, the intercept and slope) are also coupled. Finally, the structure function f X is also coupled with intra-particle mass transfer, η p , a relation that is indirectly explored in Section 4.1.

4. Discussion

4.1. The Conversion Function, f X

To determine the value of A i in Equation (1), the effect of the conversion function, f X , has to be decoupled. In Section 3.3, the rate of combustion was taken at the conversion of 0.1, often close to the maximum rate in a given experiment. Often, a small but arbitrarily selected value of X is used to identify the initial rate of combustion. The goal is to find a rate corresponding to a situation where only a small amount of char has combusted so that the porous structure of the char can be assumed relatively unchanged, so as to then approximate the value of the conversion function as f X = 0 = 1 . Rates at very low values of X are avoided to mitigate the influence of oxygen adsorbed on the surface of the char, reacting with the fuel upon introduction to the hot fluidised bed [10]. Additionally, the experimental noise is often more pronounced at the start and end of the experiment; thus, unless a good resolution in the detection of the combustion products is ensured, the starting sections of experiments should be avoided.
Taking X = 0.1 in Section 3.3 to establish the combustion regimes was an arbitrary choice. However, from Figure 1a, in a typical experiment, an increase in the rates was observed with X up to a value of ~0.15, followed by a decrease in the rate until complete burnout. The rise in the rate during the initial stage might be attributed to an increase in the char surface area from the gradual improvements in the access to micropores [11], whereas the drop in the rate after the maximum is because of the decrease in the net surface area caused by the enlargement and coalescence of pores [12].
To account for these effects in our work, an f X function was fitted empirically to the experiment at which the particle combusted in Regime I, i.e., in the absence of an intra-particle gradient in the gas concentration. A similar approach was adopted in the literature earlier [13,14,15] when discussing the gasification of coal chars in CO2. The conversion function, f X , was determined empirically as the ratio of the instantaneous rate at conversion X to the initial rate, taken here at X = 0.01 :
f X = r X / r X = 0.01
Equation (7) can be understood as a “normalised rate”. Examples of such normalised rates from experiments at temperatures between 723 and 923 K in an oxygen concentration of 21 vol% are given in Figure 4a. The normalised rates describe the structural evolution in char particles only if experiments are controlled by the reaction kinetics [13]. Changes in shape indicate a change in the regime [16]. Such changes are expected when oxygen delivery to the reactive sites in the char particle becomes too slow. Then, the effect from the opening of pores and from the evolution of the internal porosity becomes obscured by the diffusional effects, indicating the transition from Regime I (internal control) to Regime II (intraparticle control).
In Figure 4a, the normalise trends show no significant variation in trends between 723 and 773 K, but, with even higher temperatures, the maximum values in the normalised curves shift towards higher values of X. We can, thus, indirectly confirm the temperature points for the transition from Regimes I to II, demonstrating a simple approach to assessing the controlling phenomena in combustion. The temperature range for Regime I in Figure 4a agrees with that in Figure 3a.
To represent f X , we performed a multi-order regression in Excel (LINEST function), fitting the polynomials to the normalised rate from an experiment at 723 K, and the results are presented in Figure 4a. Clearly, the higher the order of the model, the more accurate the representation and fit of the f X curve to the shape of the normalised rate. Still, all three models reach high correlation coefficients, R2, with, the seventh-order polynomial giving the best overall fit. Taking the experimental noise into account, all three models are deemed satisfactory. The full regression equations for f X are given in the Appendix A.

4.2. Effect of Intraparticle Mass Transfer

Thus far, our simplistic assessment of the effect of intraparticle mass transfer is embodied by the assumption of η p inherently being a unity, corresponding to Regime I, where combustion is limited by the intrinsic kinetics. However, this assumption is worth revisiting, as η p = 1 might not be valid, especially close to the transition temperature from Regimes I to II. To decouple the effect of the intraparticle mass transfer, the effectiveness factor, η p , is calculated from the Thiele modulus, ϕ , a dimensionless number based on the ratio of the maximum rate of reaction under a purely kinetic limitation to the rate of diffusional transport through the porous structure of the char particle [17]:
ϕ = d p n + 1 k i 2 D p , E ρ C m r , C R T   P O 2 n 1
where D p , E is the effective diffusion coefficient in the pore, defined as the product of the molecular diffusion coefficient with the porosity of the char, D ε . Then, k i is the intrinsic rate constant to represent the chemical reaction. For a spherical particle, the effectiveness factor for intraparticle mass transfer, η p , can be calculated as
η p = f c ϕ 1 tanh 3 ϕ 1 3 ϕ
where f c is the correction factor to account for the effect of the reaction order being different from 1, and it is defined as [17]
f c = 1 + 1 / 2 2 ϕ 2 + 1 2 ϕ 2 0.5 1 n 2
The effectiveness factor, as described in Equation (9), is an implicit function of the intrinsic kinetics; hence, a numerical or iterative approach is needed to decouple the effect of intraparticle diffusion. When evaluating η p and k i , the activation energy E i is taken to be constant, using the values from Table 2 obtained for Regime I; hence, only the pre-exponential factor A i in k i is iterated. However, another unknown that adds uncertainty is the order of the reaction with respect to oxygen, n . From the literature, the values of n usually range from 0.6 to 1.0 [18]. Here, we take n = 0.8 and represent the connected uncertainty as bounds for the resulting value of η p , as presented in Figure 5.
The results in Figure 5 consider particles of char (250–300 μm) combusting in 21 vol% O2, with η p describing the influence of intraparticle diffusion at conversion of X = 0.1, i.e., close to the maximum observed rate for experiments in the low-temperature range. Our intention in taking X = 0.1 is to assess the section of the experiment most prone to being limited by intraparticle mass transfer (worst-case scenario), pinpointing the transition between Regimes I and II. In experiments at higher temperatures, the maximum rates shift beyond X of 0.1, but this is less important as, at these conditions, the limitations arise from the external mass transfer, transitioning to or within Regime III. Finally, when calculating η p from the intercepts in the Arhenius plots, we also account for the structure evolution function, f(X), taken as the O(7) polynomial from Figure 4.
As shown in Figure 5, η p decreases with the combustion temperature because the intrinsic rate of reaction increases more rapidly than the rate of intraparticle mass transport; hence, the latter phenomenon becomes rate-limiting. A non-uniform profile of O2 within the particles of char is formed at a higher temperature; a lower value of η p , therefore, is expected. As previously discussed in Section 4.1, the polynomial fitted f X is normalised at X = 0.01 ; hence, at X = 0.1 , the value of f X is greater than 1, and, here, f(0.1) = 1.27.

4.3. Intrinsic Kinetics of Combustion

With the effect of the intraparticle mass transfer quantified, the intrinsic kinetics of the combustion are evaluated at Regime I using the proposed O(7) model for f X and the kinetic parameters (pre-exponential factor and activation energy) are presented in Table 3.
The activation energy values E i found here, as shown in Table 3, are in the range of 140–155 kJ mol−1. For the combustion of lignocellulosic char, Di Blasi [19] summarised that the E i is in the range between 76 and 229 kJ mol−1, where the activation energy in this work falls. In more recent work, Recari et al. [20] extracted the kinetics of combustion of spruce chars at 723 K in a 5 vol% O2 environment using a TGA, assuming a shrinking-core model, and found that the activation energy was ~97 kJ mol−1, lower than the value in Table 3. In a fluidised bed, assuming a volumetric model, i.e., setting f X = 1 X , Morin et al. [4] found that the pre-exponential factor and the activation energy were 8.0 ± 1.3 × 10 8 bar−n s−1 and 144.3 kJ mol−1, respectively—similar to the values presented here, regardless of the differences in the fuel type. Similar to the work here, but using a TGA instead of a fluidised bed reactor, Peterson and Brown [21] studied the combustion of biomass char in a wide range of temperatures to establish the reaction regime. In Regime I, the pre-exponential factor and activation energy from the combustion of hardwood were found to be in the range of 3.12 ×   10 6 to 1.08 × 10 8 s−1 and 109 to 129 kJ mol−1, respectively, which are also consistent with the values reported in Table 3 in this work.

4.4. Effect of Heat Transfer

When a char particle is combusted in a bed fluidised by air, the temperature of the char particle will increase above the bed temperature. The temperature of the char particle studied here can be estimated by a heat balance at steady state, assuming that (a) the temperature within the char particle is uniform (the validity of this assumption was demonstrated in our previous work [22]); (b) the solid matrix and the gas within the char pores are always in thermal equilibrium; (c) the heat evolved by combustion to the bed temperature is negligible due to the small mass of char, as measured by the thermocouple in the reactor; (d) the contribution of transpiration to the overall heat transfer is small; and (e) CO2 is the only product from the combustion of char.
As the enthalpy of forming CO2, Δ H C O 2 , is greater than the enthalpy of forming CO, Δ H C O , the assumptions above would result in the maximum temperature deviation between the bed T b and particle temperature T p . The energy balance can be written as
r Δ H C O 2 = π d p 2 h T p T b + π d p 2 ϵ σ T p 4 T b 4
where r is the rate of combustion for a single particle of char; ϵ is the radiative emissivity of char, taken as 0.9 [23]; and σ is the Stefan–Boltzmann constant. The heat transfer coefficient, h , is calculated from the Nusselt number, N u = h d p / k a i r , where k a i r is the thermal conductivity of air. As the particle sizes of the bed materials and char are approximately the same, the N u number was calculated using the correlation by Baskakov et al. [24]:
N u = 2 1 1 ε m f 1 / 3 + 0.117 A r 0.39 P r 0.33
where A r and P r are the Archimedes number and Prandtl number, respectively. For combustion at a bed temperature of 773 K, using the rate of combustion at the conversion of 0.1, the particle temperature was estimated to be ~774 K, only a degree greater than the bed temperature. Therefore, the temperature of the bed particle and char particle could be assumed as the same temperature, and the effect of heat transfer could be neglected.

5. Conclusions

In summary, this study demonstrates a step-by-step approach towards the decoupling of the mass transfer phenomena from the apparent kinetics in order to find the intrinsic kinetic parameters of biochar combustion. The obtained activation energy for combustion in air, Ei = 154 ± 19 kJ mol−1, is similar to the values commonly quoted for coal chars. The value of E i drops quickly when moving to Regime II (internal transfer limitation), with E i of 57 ± 8 kJ/mol, and then to Regime III (external transfer limitations), with E i as low as 9 ± 5 kJ/mol. The obtained intrinsic pre-exponential constant Ai depends on the reaction model and thus is more difficult to compare. Here, we show that Ai is a strong function of the intraparticle effectiveness factor and f(X) and is also strongly coupled to E i . As expected, the combustion of biochar is fast and thus the intrapacticle mass transfer becomes quickly prominent—in our study, already at 798 K. As we propose here, a useful approach when investigating the kinetics of biochar combustion in fluidised beds is to first determine the temperature range for the three combustion regimes, I, II and III, in order to minimise the effects of mass transfer deliberately.

Author Contributions

S.G.N. and K.Y.K., experiments, data analysis, writing; E.J.M., methodology, resources, data analysis, writing—original draft preparation, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

A set of all collected experimental and supporting data for this work can be found at https://www.repository.cam.ac.uk.

Acknowledgments

D. Balan is thanked for the support in the experimental work.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The sixth-, seventh-, and eighth-order polynomial fitting in Figure 4 is given by Equations (A1)–(A3), respectively. Coefficients have been rounded to the nearest whole number, expect for the zeroth-order coefficient, which has been rounded to two decimal places.
f X = 82 X 6 + 280 X 5 374 X 4 + 248 X 3 83 X 2 + 11 X + 0.74
f X = 410 X 7 1567 X 6 + 2414 X 5 1917 X 4 + 833 X 3 192 X 2 + 19 X + 0.62
f X = 1076 X 8 + 4846 X 7 9047 X 6 + 9051 X 5 5227 X 4 + 1747 X 3 320 X 2 + 27 X + 0.55

References

  1. Kwong, K.Y.; Marek, E.J. Combustion of Biomass in Fluidized Beds: A Review of Key Phenomena and Future Perspectives. Energy Fuels 2021, 35, 16303–16334. [Google Scholar] [CrossRef]
  2. Spliethoff, H. Power Generation from Solid Fuels; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  3. Mueller, A.; Haustein, H.D.; Stoesser, P.; Kreitzberg, T.; Kneer, R.; Kolb, T. Gasification Kinetics of Biomass- and Fossil-Based Fuels: Comparison Study Using Fluidized Bed and Thermogravimetric Analysis. Energy Fuels 2015, 29, 6717–6723. [Google Scholar] [CrossRef]
  4. Morin, M.; Pécate, S.; Hémati, M. Kinetic study of biomass char combustion in a low temperature fluidized bed reactor. Chem. Eng. J. 2018, 331, 265–277. [Google Scholar] [CrossRef]
  5. McMurchie, E.J.; Olatunji, J.; Kwong, K.Y.; Marek, E.J.; Hayhurst, A.N. The kinetics and mechanism of the combustion of the carbon in biochar from oak, as studied in an electrically heated, fluidised bed of sand. Fuel 2024, 378, 132677. [Google Scholar] [CrossRef]
  6. Kwong, K.Y.; Kwong, K.Y.; Harrison, A.R.P.; Harrison, A.R.P.; Gebers, J.C.; Gebers, J.C.; Dennis, J.S.; Dennis, J.S.; Marek, E.J.; Marek, E.J. Chemical Looping Combustion of a Biomass Char in Fe2O3-, CuO-, and SrFeO3−δ-Based Oxygen Carriers. Energy Fuels 2022, 36, 9437–9449. [Google Scholar] [CrossRef]
  7. Kwong, K.; Mao, R.; Scott, S.; Dennis, J.; Marek, E. Analysis of the rate of combustion of biomass char in a fluidised bed of CLOU particles. Chem. Eng. J. 2021, 417, 127942. [Google Scholar] [CrossRef]
  8. Scala, F. Mass transfer around freely moving active particles in the dense phase of a gas fluidized bed of inert particles. Chem. Eng. Sci. 2007, 62, 4159–4176. [Google Scholar] [CrossRef]
  9. Levenspiel, O. Chemical Reaction Engineering, 2nd ed.; Wiley: New York, NY, USA, 1972. [Google Scholar]
  10. Fennell, P.; Dennis, J.; Hayhurst, A. The order with respect to oxygen and the activation energy for the burning of an anthracitic char in O2 in a fluidised bed, as measured using a rapid analyser for CO and CO2. Proc. Combust. Inst. 2009, 32, 2051–2058. [Google Scholar] [CrossRef]
  11. Bhatia, S.K.; Perlmutter, D.D. A random pore model for fluid-solid reactions: I. Isothermal, kinetic control. AIChE J. 1980, 26, 379–386. [Google Scholar] [CrossRef]
  12. Ballal, G.; Zygourakis, K. Evolution of pore surface area during noncatalytic gas-solid reactions. 2. Experimental results and model validation. Ind. Eng. Chem. Res. 1987, 26, 1787–1796. [Google Scholar] [CrossRef]
  13. Saucedo, M.A.; Lim, J.Y.; Dennis, J.S.; Scott, S.A. CO2-gasification of a lignite coal in the presence of an iron-based oxygen carrier for chemical-looping combustion. Fuel 2014, 127, 186–201. [Google Scholar] [CrossRef]
  14. Dai, P.; González, B.; Dennis, J.S. Using an experimentally-determined model of the evolution of pore structure for the gasification of chars by CO2. Fuel 2016, 171, 29–43. [Google Scholar] [CrossRef]
  15. Marek, E.J.; Zheng, Y.; Scott, S.A. Enhancement of char gasification in CO2 during chemical looping combustion. Chem. Eng. J. 2018, 354, 137–148. [Google Scholar] [CrossRef]
  16. Zheng, Y.; Marek, E.; Scott, S.A. The effect of CO on CO2-char gasification. Proc. Combust. Inst. 2021, 38, 4261–4269. [Google Scholar] [CrossRef]
  17. Hong, J.; Hecker, W.C.; Fletcher, T.H. Improving the Accuracy of Predicting Effectiveness Factors for mth Order and Langmuir Rate Equations in Spherical Coordinates. Energy Fuels 2000, 14, 663–670. [Google Scholar] [CrossRef]
  18. Haugen, N.E.L.; Loong, B.K.Y.; Mitchell, R.E. Numerical approaches for thermochemical conversion of char. Prog. Energy Combust. Sci. 2022, 91, 100993. [Google Scholar] [CrossRef]
  19. Di Blasi, C. Combustion and gasification rates of lignocellulosic chars. Prog. Energy Combust. Sci. 2009, 35, 121–140. [Google Scholar] [CrossRef]
  20. Recari, J.; Berrueco, C.; Abelló, S.; Montané, D.; Farriol, X. Effect of temperature and pressure on characteristics and reactivity of biomass-derived chars. Bioresour. Technol. 2014, 170, 204–210. [Google Scholar] [CrossRef]
  21. Peterson, C.A.; Brown, R.C. Oxidation kinetics of biochar from woody and herbaceous biomass. Chem. Eng. J. 2020, 401, 126043. [Google Scholar] [CrossRef]
  22. Kwong, K.Y.; Dennis, J.S.; Marek, E.J. Application of particle-scale modelling to the combustion of a char particle in a fluidised bed of CLOU particles. Fuel 2023, 342, 127902. [Google Scholar] [CrossRef]
  23. von Berg, L.; Soria-Verdugo, A.; Hochenauer, C.; Scharler, R.; Anca-Couce, A. Evaluation of heat transfer models at various fluidization velocities for biomass pyrolysis conducted in a bubbling fluidized bed. Int. J. Heat Mass Transf. 2020, 160, 120175. [Google Scholar] [CrossRef]
  24. Leckner, B. Heat and mass transfer to/from active particles in a fluidized bed—An analysis of the Baskakov-Palchonok correlation. Int. J. Heat Mass Transf. 2021, 168, 120860. [Google Scholar] [CrossRef]
Figure 1. Instantaneous rate of combustion, per unit mass of biochar, when varying the mass of biochar (6–20 mg and particle size of 200–250 µm) added to a fluidised bed: (a) when 11 vol% O2/N2 is used as a fluidising gas in a bed of 823 K, tracing the instantaneous rate against char conversion; (b) when different concentrations of O2 (vol%) and temperatures are used, and the presented instantaneous rate refers to the rate at X = 0.1 conversion.
Figure 1. Instantaneous rate of combustion, per unit mass of biochar, when varying the mass of biochar (6–20 mg and particle size of 200–250 µm) added to a fluidised bed: (a) when 11 vol% O2/N2 is used as a fluidising gas in a bed of 823 K, tracing the instantaneous rate against char conversion; (b) when different concentrations of O2 (vol%) and temperatures are used, and the presented instantaneous rate refers to the rate at X = 0.1 conversion.
Processes 12 02103 g001
Figure 2. Instantaneous rate of combustion per unit mass of biochar at X = 0.1 and T = 773 K for varied diameters of biochar particles, with bed fluidised by air. The green markers denote the experimental results and the orange ones the rate of external mass transfer, with S h calculated by Equation (5).
Figure 2. Instantaneous rate of combustion per unit mass of biochar at X = 0.1 and T = 773 K for varied diameters of biochar particles, with bed fluidised by air. The green markers denote the experimental results and the orange ones the rate of external mass transfer, with S h calculated by Equation (5).
Processes 12 02103 g002
Figure 3. Logarithm of the instantaneous rate of combustion per unit mass of biochar ( d p = 250–300 µm), at X = 0.1, at different temperatures, with O2 concentrations of (a) 21 vol% O2 and (b) 5 vol% O2. The combustion reaction was classified into different reaction regimes based on the gradients in the plot. The dashed lines represent the linear lines of best fit, to guide the eye.
Figure 3. Logarithm of the instantaneous rate of combustion per unit mass of biochar ( d p = 250–300 µm), at X = 0.1, at different temperatures, with O2 concentrations of (a) 21 vol% O2 and (b) 5 vol% O2. The combustion reaction was classified into different reaction regimes based on the gradients in the plot. The dashed lines represent the linear lines of best fit, to guide the eye.
Processes 12 02103 g003
Figure 4. The rate of reaction normalised to the initial rate, i.e., X = 0.01, at different conversion levels for biochar (250–300 µm) in 21 vol% O2: (a) at various temperatures, (b) at 723 K and compared to f X models obtained by fitting 6th- to 8th-order polynomial regressions.
Figure 4. The rate of reaction normalised to the initial rate, i.e., X = 0.01, at different conversion levels for biochar (250–300 µm) in 21 vol% O2: (a) at various temperatures, (b) at 723 K and compared to f X models obtained by fitting 6th- to 8th-order polynomial regressions.
Processes 12 02103 g004
Figure 5. Intraparticle effectiveness factor, η p , evaluated during combustion of char particles (250–300 μm) at different temperatures in 5 and 21 vol% O2, at conversion X = 0.1. The error bars represent the effectiveness factor calculated using the order of reaction n = 0.6 and 1.0. The regimes of the reaction are approximated as labelled.
Figure 5. Intraparticle effectiveness factor, η p , evaluated during combustion of char particles (250–300 μm) at different temperatures in 5 and 21 vol% O2, at conversion X = 0.1. The error bars represent the effectiveness factor calculated using the order of reaction n = 0.6 and 1.0. The regimes of the reaction are approximated as labelled.
Processes 12 02103 g005
Table 1. Experimental campaigns conducted to discern mass transfer and reaction kinetics.
Table 1. Experimental campaigns conducted to discern mass transfer and reaction kinetics.
CampaignTemperature (K)Mass of Char (mg)Particle Size
(µm)
Concentration
of O2 (vol%)
Purpose
1823, 8736, 10, 14, 19, 21200–2505, 11, 21assessment of the influence of interphase mass transfer on the apparent kinetics
27736, 12200–255,
250–300,
300–355,
355–422,
422–500,
500–600
21influence of interparticle (external) mass transfer on the apparent kinetics
3723, 748, 773, 798, 823, 848, 863, 873, 883, 898, 923, 948, 973, 1023, 1073, 112310 1250–3005, 21identification of the combustion regimes,
assessment of the intrinsic kinetics of combustion
1 Doubled and halved masses were also used to study the influence of interphase mass transfer.
Table 2. Extracted activation energies and apparent pre-exponential factors for different regimes of combustion in different concentrations of O2, according to Equation (6). The values in parentheses give the estimated parameters within the standard error, assessed from the experimental results. Regime I describes combustion limited by intrinsic kinetics, Regime II corresponds to limitations from intraparticle mass transfer, and Regime III from interparticle mass transfer.
Table 2. Extracted activation energies and apparent pre-exponential factors for different regimes of combustion in different concentrations of O2, according to Equation (6). The values in parentheses give the estimated parameters within the standard error, assessed from the experimental results. Regime I describes combustion limited by intrinsic kinetics, Regime II corresponds to limitations from intraparticle mass transfer, and Regime III from interparticle mass transfer.
ParameterConcentration of O2 (vol%)Regime
IIIIII
Temperature range (K)21<798798–873>883
5<798798–883>898
E i (kJ mol−1)21155 (126, 183)57.0 (49.0, 65.0)9.3 (4.2, 14.3)
5142 (-) *57.0 (55.8, 58.1)29.2 (24.6, 33.8)
A i η p P O 2 n f X (s−1)214.7 × 108 (2.56 × 106, 1.61 × 1010)1.26 × 102 (4.02 × 101, 3.97 × 102)1.78 × 10−1 (9.24 × 10−2, 3.44 × 10−1)
53.11 × 107 (-) *5.36 × 101 (4.53 × 101, 6.34 × 101)1.46 × 100 (8.09 × 10−1, 2.64 × 100)
* Assessed based on only two temperature points; thus, the linear fit gives R2 = 1. See the text for a discussion.
Table 3. Intrinsic kinetic parameters (activation energy and pre-exponential factors) from Regime I at different concentrations of O2, using different f X .
Table 3. Intrinsic kinetic parameters (activation energy and pre-exponential factors) from Regime I at different concentrations of O2, using different f X .
f X Modelled as O(7) Polynomial; f(0.1) = 1.27
P O 2 = 0.21 bar P O 2 = 0.05 bar
E i (kJ mol−1)154.6141.8
A i (bar−n s−1)(2.18 ± 0.06) × 109(4.33 ± 0.23) × 108
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Newman, S.G.; Kwong, K.Y.; Marek, E.J. A Systematic Approach to Determining the Kinetics of the Combustion of Biomass Char in a Fluidised Bed Reactor. Processes 2024, 12, 2103. https://doi.org/10.3390/pr12102103

AMA Style

Newman SG, Kwong KY, Marek EJ. A Systematic Approach to Determining the Kinetics of the Combustion of Biomass Char in a Fluidised Bed Reactor. Processes. 2024; 12(10):2103. https://doi.org/10.3390/pr12102103

Chicago/Turabian Style

Newman, S. G., K. Y. Kwong, and E. J. Marek. 2024. "A Systematic Approach to Determining the Kinetics of the Combustion of Biomass Char in a Fluidised Bed Reactor" Processes 12, no. 10: 2103. https://doi.org/10.3390/pr12102103

APA Style

Newman, S. G., Kwong, K. Y., & Marek, E. J. (2024). A Systematic Approach to Determining the Kinetics of the Combustion of Biomass Char in a Fluidised Bed Reactor. Processes, 12(10), 2103. https://doi.org/10.3390/pr12102103

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

Article Metrics

Back to TopTop