The combustion of biomass not only falls in energy production, but also in the recovery of waste. The treatment method most used for the recovery of waste is incineration because this method of treatment can minimize the volume of waste. In this work, it comes to realize a numerical modeling of the combustion of biomass in a fixed grate furnace. A literature review allowed us to describe the stages of combustion in terms of mathematical equations. Taking into account the results of elemental analysis and immediate analysis, solid and gaseous species used to simulate their transport equations are: Dry fuel (biomass), char, CH 4, O 2, CO, H 2O, CO 2, and N 2. From equations of energy transportation, we deducted the TS temperature of the solid fuel bed and Tg of gas. Subsequently, we simulated the resolution 1-D transport equations using a computer code written by us and this on the basis of mathematical modeling of the transport equations. This 1-D unstationnary model takes into account the different stages of load transformation. In this calculation code, we used the explicit Euler method for space discretization, and for the time resolution, we used an implicit method which solves stiff problems of differential equations to ordinary derivatives. The results are satisfactory because the calculated numerical profiles follow the experimental profiles, such as, the temperature profiles, the loss of mass of the fuel bed and the speed of propagation of the flame front.
Nowadays, the energies have an important place in the life of man and our current level of development requires an energy consumption increasingly growing. To meet this great need for energy, man has long had use of fossil fuels, rejecting into the atmosphere huge amounts of pollutants, including CO2. Given the need to protect the environment by reducing the emission of greenhouse gases (GHGs), uncertainty about the uncertain future of fossil fuels, and especially our strong energy demand research have turned to new energy sources, which one of the main is combustion of biomass.
Renewable energies are numerous (biomass, solar, geothermal, waste, wind, hydraulic) and energy research grants them special attention. Here we focus especially on the use of biomass as solid fuel that can be burned in a suitable oven to harness the energy released. The technologies used for burning biomass are numerous and the choice of a burner may depend on such factors as the nature of the fuel and its moisture content.
The combustion of biomass not only falls in energy production, but also in the recovery of waste. The treatment method most used for the recovery of waste is incineration because this method of treatment can minimize the volume of waste.
The multiphase combustion is a phenomenon that is part of our daily life: Cooking food using firewood or sawdust furnaces, incineration of household waste (HW), to name a few. This phenomenon, which seems simple, actually falls a physicochemical complexity. The multiphase combustion involves the phases solid, liquid and gaseous, noted, however, that obtaining a good combustion quality depends on the quality of fuel distribution on the grid, on the turning efficiency ensuring the best touch fuel/oxidizer to a sufficient temperature in drying zone-pyrolyse and on the nature of the fuel used. Several studies have already been conducted in this area: This is the case of Shin and Choi [
This work is in search of a numerical method to solve equations, and for approaching the experimental results. So we proposed writing a numerical computation code that allows to numerically modeling the combustion of solid biomass particles. This numerical modeling is a combustion control to optimize the conversion of fuel into heat and reduce the production of pollutants. Controlling combustion involves controlling its various stages, namely drying, pyrolysis and combustion of the gases and the carbon residue released during the pyrolysis phase. One of the peculiarities of our numerical model is to use a stiff method for the resolution of the temporal part, which makes the numerical model more stable and more robust.
The combustion of a load is characterized by a complex set of exothermic chemical reactions due to rapid oxidation of the fuel elements. The kinetics of the combustion of a load involves several stages from the drying to the heterogeneous and homogeneous combustion.
On a smoldering wood particle (
Combustion is a complex process that involves several related equations together and acting on two levels: In solid phase and gas phase.
There is a significant interaction between the combustion in the solid phase and in the gas phase. This interaction is shown in
The main phenomena that take place in the solid phase are: drying or evaporation, pyrolysis and coal combustion. The heat transfer between solid particles here is by conduction and radiation. In this phase, the moisture produced during the drying as well as the volatiles produced during the pyrolysis and the carbon monoxide produced during the combustion of the coal will be returned to the gas phase. In the gas phase, there is therefore combustion of gaseous products from the solid phase. As another interaction, it is the gaseous phase that provides the solid phase with the oxygen necessary for coal combustion, and the heat transfer between these two phases is by convection.
The combustion process thus leads to the destruction of the solid fuel. During its degradation in an incinerator, or a burner, the solid fuel degrades in the presence of oxygen. This combustion results in a gradual decrease of the height of the solid fuel bed into the incinerator, as shown in
To implement all these steps, here we have conducted a 1-D modeling of the combustion of the wood particles having a constant diameter and assumed the same for all particles.
The fuel is made of wooden particles having a density ρ = 600 kg/m3, of radius of between 1 and 3 cm and its PCS is 19.8975 MJ/kg. The height of the fuel bed is 45 cm. The air enters under the grate at a speed of 0.05 m/s.
The conservation equations are written for solid and gaseous phases. These conservation equations applied to a volume element can describe the instantaneous evolution of different variables anywhere in the field. One thus obtains transport equations whose general form is:
∂ ρ ∅ ∂ t ︸ unstationnary term + ∇ ⋅ ρ u ∅ ︸ convection term = ∇ ⋅ ( J ∅ ) ︸ diffusion term + S ∅ ︸ source term (1)
where ∅ represents the variables of mass, the total enthalpy (energy) or species; J ∅ represents the broadcast stream, S ∅ source terms, ρ the densities of the solid or gaseous species and u the speeds of evolution of species.
View test results in
H2O: ∂ ( ε b ρ g Y H 2 O ) ∂ t + ∂ ( ρ g ν g Y H 2 O ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y H 2 O ∂ x ) + R Y H 2 O (2)
O2: ∂ ( ε b ρ g Y O 2 ) ∂ t + ∂ ( ρ g ν g Y O 2 ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y O 2 ∂ x ) + R Y O 2 (3)
CH4: ∂ ( ε b ρ g Y CH 4 ) ∂ t + ∂ ( ρ g ν g Y CH 4 ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y CH 4 ∂ x ) + R Y CH 4 (4)
CO2: ∂ ( ε b ρ g Y CO 2 ) ∂ t + ∂ ( ρ g ν g Y CO 2 ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y CO 2 ∂ x ) + R Y CO 2 (5)
CO: ∂ ( ε b ρ g Y CO ) ∂ t + ∂ ( ρ g ν g Y CO ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y CO ∂ x ) + R Y CO (6)
| Immediate analysis (wt%) | Elemental analysis (wt%) | ||
|---|---|---|---|
| Humidity | 8.2% gross | C | 46.9% dry 49.1% gross |
| Volatile rate | 77.1% dry 69.7% gross | H | 6.4% dry 6.68% gross |
| Ash content | 1.2% dry 1.3% gross | O (by difference) | 38.3% dry 43.86% gross |
| Fixed carbon content | 13.5% dry 29% gross | N | 0.59% dry 0.36% gross |
N2: ∂ ( ε b ρ g Y N 2 ) ∂ t + ∂ ( ρ g ν g Y N 2 ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y N 2 ∂ x ) (7)
Energy Transport equation in:
Gas phase: ∂ ( ε b ρ g c p g T g ) ∂ t + ∂ ( ρ g ν g c p g T g ε b ) ∂ x = ∂ ∂ x ( λ g ∂ T g ∂ x ) + h S ( T S − T g ) + Q g + Q l g (8)
Solid phase: ∂ ( ( 1 − ε b ) ρ s c p s T s ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s c p s T s ) ∂ x = ∂ ∂ x ( λ e f f ∂ T s ∂ x ) + h S ( T g − T s ) + Q s + Q l s (9)
Humidity: ∂ ( ( 1 − ε b ) ρ s Y h u m ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s Y h u m ) ∂ x = ∂ ∂ x ( ( 1 − ε b ) D s ρ s ∂ Y h u m ∂ x ) + R Y h u m (10)
Dry solid fuel: ∂ ( ( 1 − ε b ) ρ s Y F u e l ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s Y F u e l ) ∂ x = ∂ ∂ x ( ( 1 − ε b ) D s ρ s ∂ Y F u e l ∂ x ) + R Y F u e l (11)
Char: ∂ ( ( 1 − ε b ) ρ s Y c h a r ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s Y c h a r ) ∂ x = ∂ ∂ x ( ( 1 − ε b ) D s ρ s ∂ Y c h a r ∂ x ) + R Y c h a r (12)
Ash: ∂ ( ( 1 − ε b ) ρ s Y A s h ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s Y A s h ) ∂ x = ∂ ∂ x ( ( 1 − ε b ) D s ρ s ∂ Y A s h ∂ x ) + R Y A s h (13)
Continuity equation in the gas phase: ε b ∂ ( ρ g ) ∂ t + ε b ∂ ( ρ g ν g ) ∂ x = R g (14)
Continuity equation in the solid phase: ∂ ( ( 1 − ε b ) ρ s ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s ) ∂ x = R s (15)
During combustion of the wood particles, the phenomena of drying, pyrolysis, homogeneous and heterogeneous combustion act to reduce fuel in combustion products. The following equations will model these different physico-chemical phenomena.
Drying is the evaporation process of the water contained in the fuel. It corresponds to the reaction:
Wetbiomass → Drybiomass + Moisture (16)
Drying the reaction rate is given by the Arrhenius equation [
R v a p = A v a p e − E a v a p R T s Y h ρ a p (17)
where A v a p , E a v a p , TS, Yh and ρ a p are respectively the pre-exponential factor ( A v a p = 5.13 × 10 10 [
Pyrolysis is the degradation of the organic portion of the dry fuel under the effect of heat and leading to the formation of carbonaceous residue (coal) and the release of volatile gases. We have considered here for the process of thermal degradation of the wood particles in a two-step mechanism (blasi [
ü The first step corresponds to the production of volatile by the following reaction:
Solidfuel → Volatile (18)
The speed of this devolatilization reaction is given by (blasi [
R d e v = 5.6 × 10 6 exp ( − 10700 T S ) M S _ F u e l (19)
where M S _ F u e l is the mass of dry fuel present in each moment in the combustion chamber and TS the temperature of the solid fuel.
Here the volatile is supposed to be CH4 [
ü The second step is the production of coal according to the reaction:
Solidfuel → Char (20)
The char production rate is given by [
R c h a r = 2.66 × 10 10 exp ( − 12800 T S ) M S _ F u e l (21)
Taking account of the elemental analysis of this solid fuel, the gaseous species that undergo homogeneous combustion (CH4, CO) burn according to the reactions:
CH 4 + 3 2 O 2 → CO + 2 H 2 O [
With the reaction rate [
R g 1 = A 2 [ CH 4 ] a [ O 2 ] b exp ( − E a 2 T g ) (23)
CO + 1 2 O 2 → CO 2 [
With the rate of reaction [
R g 2 = A 3 [ CO ] [ H 2 O ] a [ O 2 ] b exp ( − E a 3 T g ) (25)
α = 0.5 is a weighting factor, the average temperature T e to which the reaction rate is calculated is given by [
T e = α T g + ( 1 − α ) T s si T g ≤ T s ; T e = T g si T g > T s (26)
The constants Ai, E a i / R (i = 2, 3, 4), a and b are given in
Char combustion may be described by the following two processes:
C + O 2 → CO 2 and 2 C + O 2 → 2 CO (27)
By combining these two equations, we obtain: [
C + φ O 2 → 2 ( 1 − φ ) CO + ( 2 φ − 1 ) CO 2 (28)
The constant φ is given by:
φ = 2 + f 2 ( 1 + f ) (29)
where f is the CO/CO2 ratio ( f = CO CO 2 ). It is given by [
f = 12 e − 3300 T s (30)
The speed of the carbon combustion reaction is written as [
R C ( s ) = π d p 2 M C k C ( s ) P O 2 φ (31)
where d p is the diameter of the particle, Mc the molar mass of carbon, P O 2 the partial pressure of oxygen and k C ( s ) the rate constant which is given by [
k C ( s ) = 1 1 k c Ω a c t + 1 k d (32)
This rate constant has three terms, namely:
® The ash factor which takes into account the effect of the presence of ash on the oxidation kinetics and whose expression can be written [
Ω a c t = Y C ( s ) Y C ( s ) + Y a s h (33)
| Reaction | Ai | Eai/R[K] | a | b | Refs. |
|---|---|---|---|---|---|
| Rg1 | 1.6 × 1010 | 24,157 | 0.7 | 0.8 | [ |
| Rg2 | 1.3 × 1011 | 15,105 | 0.5 | 0.5 | [ |
where Y C ( s ) and Y a s h are respectively the mass fractions of carbon residue (char) and ash.
® The rate constant due to chemical kinetics which is given by [
k c = 1.175 T s e − 9000 T s (34)
where TS is the temperature of the solid.
® The mass transfer coefficient k d which is defined by [
k d = S h D O 2 , m d p (35)
where D O 2 , m is the mass diffusion coefficient of oxygen in air and Sh the Sherwood number defined by [
S h = 2 + 0.6 R e 0.5 S c 1 / 3 (36)
where Sc is the Schmidt number given by: S c = μ g ρ g D O 2 , m , and Re the Reynolds number whose expression is: R e = ρ g u g d p / μ g [
μ g = 1.98 × 10 − 5 ( T g 300 ) 2 / 3 (37)
And ρ g the density of the gas calculated at the temperature of the gas phase at a specific point of the bed; it is calculated according to the equation:
ρ g = P R T g ( Y O 2 M O 2 + Y N 2 M N 2 ) (38)
M O 2 and M N 2 are respectively the molar masses of dioxygen and dinitrogen, and P, the gas pressure. Taking air as the inlet gas, we have: Y O 2 = 23.3 % , Y N 2 = 76.7 % .
The Sherwood number can also be written [
S h = J S c 1 / 3 R e (39)
where J which allows taking into account the diffusion of the gaseous mixture flowing through the fuel bed is the Colburn factor whose expression is [
J = 1 ε b ( 0.765 R e 0.82 + 0.365 R e 0.368 ) (40)
With ε b , the porosity of the fuel bed and the expression of which is:
ε b = 1 − ρ a p p , 0 ρ s , 0 (41)
where ρ a p p , 0 is the apparent density determined experimentally ( ρ a p p , 0 = 550 kg / m 3 for wood particles), and ρ s , 0 the initial density of solid biomass particles estimated from the formula:
ρ s , 0 = 1 H b ρ w a t e r + 1 − H b ρ S _ F u e l d r y (42)
With H b , the moisture on crude supplied by immediate analysis, ρ w a t e r the density of water, and ρ S _ F u e l d r y , the density of dry solid biomass particles.
The density of the solid during the thermal conversion is calculated from the following relationship:
ρ s = 1 Y h ρ w a t e r + Y v ρ M V + Y C ( s ) ρ C ( s ) + Y A s h ρ A s h (43)
where Y h , Y v , Y C ( s ) , Y A s h are respectively the mass fractions of moisture, volatile, char and ash, and ρ w a t e r , ρ M V , ρ C ( s ) , ρ A s h , respectively the densities of water, volatiles, char and ash.
The continuity equation for the gas phase is given by [
ε b ∂ ( ρ g ) ∂ t + ε b ∂ ( ρ g ν g ) ∂ x = R g (44)
The transport equation of the gaseous species is given by [
∂ ( ε b ρ g Y i g ) ∂ t + ∂ ( ρ g ν g Y i g ε b ) ∂ x = ∂ ∂ x ( ε b D a , e f f ρ g ∂ Y i g ∂ x ) + R Y i g (45)
With ig = H2O, O2, CH4, CO2, CO, and N2.
The transport equation of the energy in the gas phase is given by [
∂ ( ε b ρ g c p g T g ) ∂ t + ∂ ( ρ g ν g c p g T g ε b ) ∂ x = ∂ ∂ x ( λ g ∂ T g ∂ x ) + h S ( T S − T g ) + Q g + Q l g (46)
where h is the heat transfer coefficient of gas-balance (W/m2K), S the specific exchange surface of the bed of particles (here S = 54 m2/m3), Qg the heat gain due to the combustion in gaseous phase ( Q g = ∑ i r i ( − Δ H i ) , r i and Δ H i the speed and the heat of reaction in the gas phase), and Qlg the loss of energy (enthalpy) along the bed wall. Qlg is given by [
Q l g = ε b ( k / L ) ( T g − T 0 ) / Δ l (47)
where k is the bed conductivity coefficient, L the length/height of the bed, T0 the room temperature, and Δ l the wall thickness.
The continuity equation for the solid phase is given by [
∂ ( ( 1 − ε b ) ρ s ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s ) ∂ x = R s (48)
The general form of the transport equation of solid species is given by [
∂ ( ( 1 − ε b ) ρ s Y i s ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s Y i s ) ∂ x = ∂ ∂ x ( ( 1 − ε b ) D s ρ s ∂ Y i s ∂ x ) + R Y i s (49)
where is = Humidity, dry fuel (solid biomass), char, ash.
The transport equation of the energy in the solid phase is given by [
∂ ( ( 1 − ε b ) ρ s c p s T s ) ∂ t + ∂ ( ( 1 − ε b ) ρ s ν s c p s T s ) ∂ x = ∂ ∂ x ( λ e f f ∂ T s ∂ x ) + h S ( T g − T s ) + Q s + Q l s (50)
Qs is the heat generation due to char combustion (W/m3) given by: Q s = r c h a r x Δ H c h a r where r c h a r and Δ H c h a r represent speed and heat of reaction of char combustion reaction in the solid phase. Qls is the loss of energy (enthalpy) along the bed wall and is given by [
Q l s = ( 1 − ε b ) ( k / L ) ( T s − T 0 ) / Δ l (51)
where k is the bed conductivity coefficient, L the length/height of the bed, T0 the room temperature, and ∆l the wall thickness.
In these transport equations, D s represents the particle mixture coefficient due to random movements of the particles in the bed, D a , e f f the axial dispersion coefficient which is given by [
D a , e f f = D O 2 , m + 0.5 v g d p (52)
And λ e f f the effective thermal conductivity coefficient given by [
λ e f f = λ e f f , 0 + 0.5 λ g P r R e / ε b (53)
D O 2 , m is the mass diffusion coefficient of oxygen in the air, v g the gas velocity, and λ e f f , 0 the conduction coefficient in the absence of gas flow.
The heat transfer coefficient by convection between the particles and the gas is:
h = N u λ g / d p (54)
where λ g is the thermal conductivity of the gas and is [
λ g = 4.8 × 10 4 T g 0.717 (55)
And Nu is the Nusselt number defined by:
N u = 2 ε b + 0.295 R e 0.625 P r 0.33 (56)
where Pr is the Prandtl number which is written:
P r = μ g C p g / λ g (57)
where c p g is the specific heat of the gas.
The reaction heats are derived from the formula:
Q r e a c t i o n = Δ H r e a c t i o n d M i g , s d t (58)
where Δ H r é a c t i o n and d M i g , s d t are respectively the enthalpy and the rate of the reaction.
To solve these equations, we propose a number of assumptions including:
- The diffusion coefficients of temperature and species are assumed to be constant;
- Particles of solid biomass used as fuel are supposed to be dry, which allows us to neglect the drying step during the combustion modeling;
- The particles are assumed to average size, which allows us to accept the hypothesis of the diffusion regime [
The new system obtained after taking into account the above assumptions is:
ü Transport equations of gas species:
ε b ρ g ∂ ( Y i g ) ∂ t + ρ g ν g ε b ∂ ( Y i g ) ∂ x = ε b D a , e f f ρ g ∂ ∂ x ( ∂ Y i g ∂ x ) + R Y i g (59)
where ig = {H2O, O2, CH4, CO2, CO, N2}.
ü Transport equations of temperatures (energy):
Gas Temperature:
ε b ρ g c p g ∂ T g ∂ t + ρ g c p g ν g ε b ∂ T g ∂ x = λ g ∂ ∂ x ( ∂ T g ∂ x ) + W t g (60)
Temperature of the solid bed:
( 1 − ε b ) ρ s c p s ∂ T s ∂ t + ( 1 − ε b ) ρ s ν s c p s ∂ T s ∂ x = λ e f f ∂ ∂ x ( ∂ T s ∂ x ) + W t s (61)
| Source terms | Expressions |
|---|---|
| R Y F u e l | − R d e v − R c h a r |
| R Y c h a r | R c h a r − R C ( s ) |
| R Y v o l a t i l e ( CH 4 ) | R d e v − R g 1 |
| R Y O 2 | − 0.5 R C ( s ) − 1.5 R g 1 − 0.5 R g 2 |
| R Y CO | R C ( s ) + R g 1 − R g 2 |
| R Y H 2 O | 2 R g 1 |
| R Y CO 2 | R g 2 |
| W t g | h CH 4 ∂ ( Y CH 4 ) ∂ t + h CO ∂ ( Y CO ) ∂ t + h H 2 O ∂ ( Y H 2 O ) ∂ t + h CO 2 ∂ ( Y CO 2 ) ∂ t + h S ( T S − T g ) + Q l g |
| W t s | h r _ d e v ∂ ( Y c o m b sec , s ) ∂ t + h r _ c h a r ∂ ( Y c o m b sec , s ) ∂ t + h r _ C ( s ) ∂ ( Y r e s c a r b , s ) ∂ t + h S ( T g − T s ) + Q l s |
As a reminder, we have seen in the theoretical tools that all the equations we are trying to solve are of the form:
a 0 i ∂ Y i ∂ t = − a 1 i ∂ Y i ∂ x + b 1 i ∂ ∂ x ( ∂ Y i ∂ x ) + W 1 i ( Y i , t ) (62)
where Y i = {species (solid Y i s or gaseous Y i g ), temperatures (of gaz T g or of solid bed T s )}.
As numerical methods we use:
ü An explicit method for space discretization and based on the explicit Euler:
The Euler method is based on an explicit finite difference method. It allows approximating the derived differential equations. This method is widely used in programming because it provides great writing simplicity.
ü An implicit method for the time resolution:
For solving the unstationnary term, we used an implicit method based on numerical differentiation formulas (NDFs) that solve stiff problems of differential equations to ordinary derivatives. For this, we put the equation in the form:
∂ ( Y i ) ∂ t = F ( Y i ) (63)
With F ( Y i ) as: F ( Y i ) = A 2 i Y i ( j + 1 ) + B 2 i Y i ( j ) + C 2 i Y i ( j − 1 ) + W 2 i ( Y i ) .
The initial conditions are as follows: The combustion process is initiated above the fuel bed. For this, the upper surface of the bed is exposed to a temperature of about 850˚C. This temperature is that of the source of ignition. The exposure of
| Reaction (Reference) | Mechanism | Reaction rate | |
|---|---|---|---|
| Wood pyrolysis (Blasi [ | Wood → Volatile | d M v o l d t = 5.16 × 10 6 exp ( − 10700 / T b ) M w o o d | (64) |
| Wood → Char ( C ) | (65) | ||
| Char combustion (Smith [ | Char ( C ) + 1 2 O 2 → CO | d M c h a r d t ) 2 = − 60 S C O 2 1 1 / k r ζ + 1 / k m k r = 2.3 T p exp ( − 11100 / T b ) | (66) |
| Volatile combustion (Siminski et al. [ | C x H y + ( 1 2 x + 1 4 y ) O 2 → x CO + y 2 H 2 O | d C v o l d t = − 59.8 T g P 0.3 C v o l 0.5 C O 2 exp ( − 12200 / T g ) | (67) |
| Carbon monoxide combustion (Bartok and Sarofim [ | CO + 1 2 O 2 → CO 2 | d C CO d t = 1.3 × 10 11 C CO C H 2 O 1 / 2 C O 2 1 / 2 exp ( − 15105 / T g ) | (68) |
the top of the bed at this temperature is done during a relatively short time necessary for the appearance of the flame front that gradually spread in the fuel bed.
For solid and gaseous species, we consider here the mass fractions and calculations are initialized at the onset of the flame front on the surface of the fuel bed.
Initial conditions (t = 0):
t = 0 and x = 0 (at the bottom of the fuel bed): Yis,0: YFuel,0 = 0.69; Ychar,0 = 0.1 and Yig,0: Yvolatile,0 = 10−2; Yo2,0 = 0.2; Yco,0 = 10−2; Yh2o,0 = 10−6; Yco2.0 = 10−6 et Yn2.0 = 0.747.
t = 0 and x < H: Tg,0 = TS,0 = T0 = 300 K, where H represents the height of the fuel bed.
t = 0 and x = H (at the top of the fuel bed): Tg,0 = TS,0 = 1123 K.
Boundary conditions:
t > 0 and x = H (at the top of the fuel bed): d Y i s d x = d T s d x = 0 ; is = {dry fuel (Solid biomass), char, ash} and d Y i g d x = d T g d x = 0 ; ig = {H2O, O2, CH4, CO2, CO, N2}.
The combustion process is initiated above the fuel bed. For this, the upper surface of the bed is exposed to a temperature of about 850˚C. This temperature is that of the source of ignition. The exposure of the top of the bed at this temperature is done during a relatively short time necessary for the appearance of the flame front that gradually spread in the fuel bed.
The combustion process is initiated above the fuel bed. For this, the upper surface of the bed is exposed to a temperature of about 850˚C. This temperature is that of the ignition source. The exposure of the surface of the bed at this temperature is for a fairly short time and necessary for the appearance of the flame front which will spread progressively in the fuel bed.
We have in
The combustion process is performed in this burner, after introduction of the fuel. At the very beginning, the fuel (biomass) is introduced into the oven, up to a height of 45 cm. The fuel introduced is supposed to be initially at room temperature.
When the upper surface of the fuel bed is exposed for a short time at a temperature of about 850˚C, this surface immediately undergoes a drying phenomenon with evaporation of moisture. After the drying phase, this high xposure temperature will cause pyrolysis of the fuel particles thus dried, allowing the release of volatiles and coal. These volatile gases will react in the gas phase. It follows a phenomenon of solid phase combustion of coal produced, and a gas phase combustion of volatiles (here mainly CH4) and CO produced during the combustion of coal and CH4.
Due to exposure to high temperatures of the upper surface of the fuel bed, the temperature in the bed increases rapidly, and the progression of the flame in the bed is gradually felt at the thermocouple T2 (see
The numerical temperature curves of thermocouples T2 to T9 have substantially the same ignition delays as the respective experimental curves G2 to G9 (see
For the position of thermocouple T2, we observe a first plateau up to about 10 minutes, at which time the mass fraction of dry fuel decreases sharply until it vanishes, proof that the pyrolysis process is complete. After this time, the mass fraction of dry fuel is kept at 0, thus forming a plateau again at this value. For the position of the thermocouple T3, the same phenomenon is observed with the only difference that the pyrolysis is carried out a little later, or about 18 minutes.This proves that the heat has progressed in the furnace, and that the rise in temperature due to this heat supply is already sufficient to cause the pyrolysis of our fuel.
This proves that the heat has progressed in the furnace, and that the rise in temperature due to this heat supply is already sufficient to cause the pyrolysis of our fuel. We then observe and successively the same phenomena for the points of observation of T4 to T9, and this at different times. These curves therefore have the same profile and differ only in the ignition delay. The lags between the ignition delays of these different profiles are due to the convection and diffusion phenomena that create a transport in the bed.
During the propagation of the flame front along the bed, there is degradation of the fuel with production of ash and energy, which causes a loss of mass. This loss of mass is characterized by a gradual decline in the volume or height of the bed during combustion. This decrease in bed height modeled in
The significant decrease in bed mass observed at each stage correspond to the temperature peaks observed previously, and these levels are due to the limited number of measuring or observation points. The loss of mass in the fuel bed is normal and can be explained by the transformation of the starting fuel into energy, ashes and burnt gases.
This figure presents the experimental results [
the particle size increases, the speed of propagation of the flame decreases. To explain this decrease in the propagation velocity of the flame front, we can say that when the particle size increases, the surface per unit mass decreases, and consequently the cooling effect of the air also decreases.In order to make the combustion stable, therefore, a higher maximum air ratio is required than in the case of smaller particles. Increasing the particle size therefore increases the drying time, and therefore the burning time of the fuel particles, and so the flame front evolves more slowly, which explains the decrease in flame front velocity when particle size increases, as shown in
The study of biomass combustion was carried out here on a fixed bed grid. For this study, we have modeled the different stages of solid combustion in terms of mathematical equations of evolution of species and energy. Multiphasic combustion is a physicochemical complexity. We used a 1-D model here. This 1-D modeling of the transport equations allowed us to propose a numerical model of multiphase combustion of biomass.
This multiphasic combustion is initiated by exposing the upper surface of the fuel bed to a temperature of about 850˚C. The flame front that emerges from this exposure to high temperatures will progress and spread gradually throughout the bed. The ignition of the combustion at the measuring points is justified by the exponential increase of the temperatures of the different thermocouples, the temperature peaks exceeding the 1000˚C. The degradation of the solid fuel during combustion is justified by the decrease of the volume or the loss of mass, which is observed by the decrease of the height of the combustible bed.
The distribution of the fuel on the grid, the size of the particle, the ignition temperature and also the supply air flow are important parameters for obtaining good combustion.
Here, the energy produced during combustion in the form of heat can be used for example in households for cooking meals, drying fruit, etc.
In view of the results obtained, and with some differences, we can say that our model reproduces the physico-chemical and heat transfer phenomena that occur during combustion in this fixed bed.
The authors declare no conflicts of interest regarding the publication of this paper.
Armand, C.T., Akong, O. and Bonoma, B. (2019) Numerical Study of Burning of Biomass in Fixed Bed. Energy and Power Engineering, 11, 35-57. https://doi.org/10.4236/epe.2019.112003
Y i g : Mass fraction of the gas species i
Y i s : Mass fraction of solid species i
u s : Lowering speed of the solid fuel [m/s]
λ s : Thermal conductivity of solid
λ g : Thermal conductivity of gas
μ g : Dynamic viscosity of gas
k: Bed conductivity coefficient
M i : Molar mass of species i [g/mol]
u g : Superficial gas velocity [m/s]
P O 2 : Partial pressure of oxygen [Pa]
h f i ° : Enthalpy of formation of species i [J/mol]
Lvap: Latent heat of vaporization of water [J/kg]
S: Exchange surface of the particles of the bed [m2/m3]
D i : Mass diffusion coefficient of species i
R Y i : Source term of the transport equation of species i
W t g : Term source of the transport equation of T g
W t s : Term source of the transport equation of T S
t: Time [s]
Y i g , 0 : Initial value of Y i g
Y i s , 0 : Initial value of Y i s
T S : Temperature of solid [K]
T g : Gas Temperature [K]
T S , 0 : Initial value of T S [K]
T g , 0 : Initial value of T g [K]
T a : Ambient temperature [K]
Q r i : Heat of the reaction ri [W]
ε b : Porosity of fuel bed
d p : Particle diameter [m]
dx: Space step [m]
ρ i : Density of body i [kg/m3]
h: Convective heat transfer coefficient [W/m2K]
H s : Humidity on dry
H b : Humidity on crude
E a : Activation energy [J]
h r _ r i : Heat of the reaction ri [J/kg]
Ω a c t : Ash Factor
Greek letters
λ s : Thermal conductivity of solid
λ g : Thermal conductivity of gas
μ g : Dynamic viscosity of gas
ξ : Thermal conductivity of solid
ε b : Porosity of fuel bed
ρ i : Density of body i
Ω a c t : Ash Factor
Abbreviations
GHGs: Greenhouse gases
HW: Household waste