Experimental and Numerical Study of Energy Losses in a Barbecue Oven in Burkina Faso ()
1. Introduction
In developing countries, the common energy source consists of biomass, which is becoming increasingly rare [1]. Indeed, several studies in the energy and environment sector are sounding the alarm about the degraded state of forest resources and its consequences on climate change [2] [3] [4] [5] [6]. The World Food Program (WFP) has estimated at 105,000 ha the regression rate of forest formations in Burkina while the Ministry of Environment and Sustainable Development stands at 107,626 ha/year [7].
In order to reduce biomass consumption, several energy-saving technologies have been developed. Indeed, from three (03) stones fireplace up to its improved version, efforts have been made especially for domestic cooking [8]. At the national level, the catalog of innovations includes several types of improved fireplaces with wood, charcoal and butane gas [9] [10]. As energy wood is increasingly scarce, researchers have tested fuels such as animal excrement [11], ethanol [12], rice crosses [13], paraffin [14] and so on. In addition, researchers have long prioritized sheet metal because of their low thermal inertia and their ease of access [15] [16]. For domestic use, this is practical because they are shock-resistant and easily transportable, but the efficiency of households between 15% and 17% is low [17]. On the other hand, for grilling activities, it is very important to take into account energy losses as well as the service life of the oven. In fact, the use of iron sheet increases energy losses and exposes the users to very high temperatures. In addition, the ovens, subjected to high temperatures, degrade very quickly [18]. As a result, energy, sanitary and hygienic constraints have been identified on the different technologies [19] [20]. These ovens consume huge amounts of charcoal and operate for a long time. Unfortunately, the techniques of charcoal production at the national level have very low performance (around 20%) [21].
In Burkina Faso, the low energy efficiency of the ovens has a negative impact on the grilling sector because the energy losses represent net losses of revenue for the actors.
Unfortunately, there are no scientific studies which highlight this problem.
Therefore, the main objective of this work is to quantify the energy losses generated by these equipments through an experimental and numerical study. The results obtained should later motivate work to improve the energy efficiency of grilling equipment in Burkina Faso.
2. Experimental Study
2.1. Materials and Methods
2.1.1. Materials
A set of experimental equipment was used to characterize the thermal behavior of the most used oven model. We have among others:
· An oven made of iron sheets identical to the most common model. The dimensioning is done according to the needs of the user. In these ovens the charcoal is directly in contact with all the walls which have a thickness of 2 mm. The charcoal pan of our oven has the following dimensions: 1.2 m × 0.5 m × 0.2 m.
Figure 1 shows the oven and Figure 2 its diagram.
· Temperature recorders and thermocouples of type K (precision: 1.5˚C), a precision balance (precision: 1 g), a RAYCAM thermal imaging camera with measurement uncertainty of 2˚C and a combustion analyzer.
2.1.2. Methods
Four kilograms of charcoal were introduced into the oven corresponding for a grilling cycle. Before starting the charcoal burning, thermocouples were placed on the outer walls of the oven, on the grate and inside the combustion chamber.
The infrared captures of the oven, made at intervals of five (05) minutes, allows to follow temperature distribution on the external walls of the oven. The measurements were repeated sufficiently under the same conditions in order to have average values. The experiment has a duration of two (02) hours and at the end of each experiment, the charcoal has been extinguished, dried and weighed to determine the mass of fuel burned. In addition, we also monitor the production of gases such as carbon monoxide (CO), carbon dioxide (CO2), methane (CH4), dihydrogen ratio on the dinitrogen (H2:N2) to analyze the quality of the combustion. Figure 3 shows the experimental device and Figure 4 the diagram of the experimental device.
2.2. Results and Discussion
2.2.1. Evolution of the Temperatures of the Different Parts of the Oven
The profiles represented are averages of the temporal evolutions of the temperatures of the various places of the oven. Standard deviations on measurements were also represented.
Figures 5-8 show that the maximum standard deviations of the temperatures of the different places of the oven do not have the same values. Indeed, that of the combustion temperature is 34˚C and a maximum relative uncertainty of 7.34%.
Figure 5. Temporal evolution of average combustion temperatures.
Figure 6. Temporal evolution of average temperatures of the bottom wall.
Figure 7. Temporal evolution of the average temperatures of the grid.
Figure 8. Temporal evolution of lateral mean parietal temperatures.
For grid temperatures the maximum standard deviation is 16.5˚C and a maximum relative uncertainty of 7.77%. The temperatures below the oven have a maximum standard deviation of 25.5˚C and a maximum relative uncertainty of 9.98%. Sidewall temperatures have a standard deviation of 18.75˚C and a maximum relative uncertainty of 9.08%. These values of the errors were obtained at the beginning of the experiments, which is to say for the low values of the temperatures.
The standard deviations are a little big for the cave and grid temperatures because the experiments on the oven are very sensitive to external conditions. Indeed, the direction and the speed of the wind, its humidity are elements which influence not only on the combustion but also on the temperatures of the walls. These factors are very random depending on the seasons, days, hours and even minutes. A change of speed or wind direction can cause large temperature variations. In addition, these variations are also due to the heterogeneity of charcoal which comes from different woods in most cases. These variations show the difficulty that lies in the experiments using ovens in the conditions of the broiler. Indeed, several studies show that many factors can influence enormously on the stove experimental studies [22] [23] [24].
2.2.2. Global Report
To make an interpretation of the overall operation of the oven we have represented the evolution of the temperatures at different places of the oven.
From Figure 9, it can be seen that the temperatures of the different parts of the oven have the same behavior than that of the fuel. However, these temperatures are at about 200˚C lower than that of the oven. The temperature under the oven remains greater than those of the side walls and those of the grid. The side
Figure 9. Temporal evolutions of the profiles of the average temperatures of the various places of the oven.
wall temperatures increase from 60.7˚C to 248.5˚C, and then decrease more slowly to 165.5˚C after two hours of experiments. Temperatures under the charcoal tank increased from 86.8˚C to 269˚C during 45 minutes and decreased to 188.45˚C. Note that the temperatures of the side walls remain higher than those of the bottom wall during the experiments. This difference in temperature is explained by the fact that the fuel rests on the bottom wall during the whole experiment.
The temperature of the grate also increases with the heat generated by the charcoal and it decreases as the temperature of the fuel decreases. It increases from 84.5˚C to 257˚C during the first 45 minutes and then decreases to 150.6˚C. The temperature oscillations of the grid are due to the gusts of wind which arrive at random. It is on this grid that the products to be grilled must be placed. The temperatures of the different walls as well as that of the grid are close; therefore, indicating that the heat lost by the different walls is the same available at the grid. The high thermal conductivity of the iron sheet facilitates thermal losses, hence the low energy efficiency of the oven. The temperature of the gate is close to that of the side walls and remains even below that of the wall from below the oven while it should be the hottest part after the fuel temperature. This clearly shows that this type of oven, although it is more used, is not energy efficient and therefore contributes to the waste of the country's energy resources, which are becoming increasingly rare. The heat produced by the fuel must be confined and directed towards the grid but is lost through the walls of the oven.
2.2.3. IFR Captures of the Oven
In order to support our analysis in the first part, we will view a few shots of the thermal camera shown through Figures 10-15 and make brief interpretations.
It is noted that at the beginning of the combustion of the charcoal, the heat generated is distributed on the grid and on the outer wall of the oven. Obviously,
Figure 11. 20 minutes after the start of the experiment.
Figure 12. 40 minutes after the start of the experiment.
Figure 13. 60 minutes after the start of the experiment.
the thermal energy seen on the wall will be lost by convective and radiative exchanges. The other thermograms show that this very important flow increases with time until about 40 minutes and then decreases slightly until the end of the
Figure 14. 90 min after the start of the experiment.
Figure 15. 120 min after the start of the experiment.
experiment. This heat is a loss of energy for grilling. Another consequence will be the thermal discomfort felt by the users. In fact, users exposed to high temperatures are subject to many diseases but to injuries due to contact with the walls of the oven which are very hot. These thermograms highlight the poor energy efficiency of this type of grill oven.
2.2.4. Analysis of Combustion
The quality of the combustion can be determined from the gases produced. Figure 16 shows the different proportions of the gases produced by combustion as a function of time.
The maximum standard deviations of the measurements are 0.2% for CO2, 0.1% for CO and 0.01% for CH4 and H2:N2. The presence of carbon monoxide indicates incomplete combustion although the grate is not loaded. This is due to insufficient air intake because for good combustion it takes an excess of air. In addition, the air entering from the top to the countercurrent with the burnt gas reduces its diffusion in the tray and causes incomplete combustion with a combustion bed essentially above.
3. Numerical Study
3.1. Physical Model
The oven consists essentially of a casing sheet iron. We considered the same dimensions as those of the study of the experimental study. The combustion chamber shown in Figure 17 will be modeled by likening it to a semi-open enclosure on its upper wall. It has a parallelepiped shape and the six (06) faces have been noted from A to F. the air inlets and the outputs of the burned gases are from the top.
3.2. Mathematical Model
3.2.1. Numerical Method
The use of this method requires the definition of a common parameter in solid and fluid regions: thermal conductance. The description of the main steps of the realization of this modelization is the following one [25]:
· Cutting the system into elementary blocks. Much of the value of the model is based on the quality of the block division.
· Determination of the different internal heat transfer modes associated with the blocks of the studied system. The blocks are defined by the nature of the materials or an equivalence and/or heat exchange modes. The analysis of these exchanges leads to installing between the nodes (centers of the blocks) connections called thermal conductors.
Figure 16. Proportions of gases produced by combustion as a function of time.
· Establishment for each block equations of the balance of heat flows and boundary conditions: definition of elementary models.
· Determination of the coefficients intervening in these different relations.
The nodal method consists of defining isothermal zones. Thus, each zone constitutes a node. We will start by making a sample after the discretization (the geometry chosen for each node is cubic).
3.2.2. Simplifying Hypotheses
We consider the following assumptions:
· Conduction transfers are unidirectional.
· The thermo-physical properties of the materials are constant.
· The ambient temperature is the same everywhere on the outside of the oven.
· The distribution of heat is homogeneous inside the oven.
· The properties of the air depend on the speed outside and inside they depend on the speed and the temperature.
The determination of the overall thermal behavior of the charcoal tank consists in determining the thermal behavior of the following points:
We have a total of 12 nodes to which we will determine the thermal behavior because the temperature at node M is known. The nodes are numbered in Figure 18. Thus, we obtain a system composed of square matrices of order 12. The nodes D1 and D2 are the front plane and the nodes E1 and E2 are in the rear plane.
In general, in each node we have the following equation [26]:
(1)
Some places of the oven exchange heat by convection and by radiation so in writing the equations in these nodes we will take it into account.
According to the heat transfer mode, the following equations are adopted:
· The conductive conductance:
(2)
· The convective conductance:
(3)
Figure 18. Characteristic points of the oven.
· The radiative conductance:
(4)
All internal walls exchange heat by conduction, convection and radiation with the fuel. We get the following relation:
(5)
With
for all the walls;
for the side walls; and
for the bottom wall of the oven.
The outer walls exchange heat by convection and radiation with the external environment. We get the following relation:
(6)
With
and
and
.
The discretization of these equations gives for the external walls:
(7)
With
;
and
.
Likewise at the internal walls, we have:
(8)
We generally have:
.
3.2.3. Heat Transfer Coefficients
· Convective heat transfer
In the characterization of thermal transfers between any fluid and a wall, the number of Nusselt (Nu) is used. Indeed:
(9)
With h the heat exchange coefficient, the characteristic length and the thermal conductivity of the fluid.
Empirical correlations make it possible to determine the Nusselt number as a function of the type and the transfer regime. In our study the type of convection inside the chamber is forced convection because of the chimney and the effect of the draft. The following relationship will then be used [27]:
(10)
and
;
Outside the combustion chamber, we have a natural convection with the external environment. Indeed, Bernard Eyglunent, in his manual estimates that in natural convection around a vertical flat plate [28]:
for
(11)
and
for
(12)
The outer and inner faces of the combustion chamber exchange heat by convection with the surrounding environment. The air velocity being not negligible, we chose the Mc Adam correlation to determine the heat exchange coefficient of these faces:
(13)
In a charcoal fireplace the correlations used for the characteristics of the air inside are [29]:
Massive heat:
(14)
Thermal conductivity:
(15)
Kinematic viscosity:
(16)
Dynamic viscosity:
(17)
The density:
(18)
(19)
· Radiant heat transfer
The radiative transfer coefficient between the outer wall of the pressure cooker and the celestial vault is determined by the expression [26].
(20)
In the case of a vertical wall with the celestial vault, the radiative Fciel form factor is deduced by the following expression [30] [31]:
(21)
where b represents a parameter, which is a function of the anisotropy of the sky. For an isotropic sky (b = 0) the radiative form factor corresponds to the value 0.5.
Among the several reported correlations to determine the temperature of the celestial vault, the one proposed by Swinbank is chosen in this work [32].
(22)
where Tamb is the ambient temperature (K).
The radiation of the coal is quantified by the following relation [29]:
(23)
3.3. Comparison of Numerical and Experimental Results and Validation of the Theoretical Model
The duration of both numerical simulation and experimental study is two (02) hours.
3.3.1. Comparison of Theoretical and Experimental Temperature Profiles of Combustion Temperatures
In order to compare our theoretical results with the experimental results, we interpolated the fuel temperature in order to integrate it into our simulation code.
Figure 19 shows that the 5th order polynomial interpolation curve makes a good approximation of the combustion temperature. In fact, apart from the first two points whose relative errors are 10.4% and 9.3%, the other values have relative errors of less than 5%. This curve will be used to determine fuel temperatures regardless of the time in the numerical study.
Figure 19. Experimental fuel temperatures and interpolation curve.
3.3.2. Comparison of Theoretical and Experimental Temperature Profiles of External Sidewall Temperatures
Figure 20 shows the theoretical and experimental temperature profiles of the sidewalls.
Both profiles have a first phase of rapid growth and a second phase of decay a little slow. During these two phases, the temperature profiles remain close, indicating that the theoretical model manages to describe the overall behavior of the side walls. The relative errors are large during the first thirty (30 min) and about 6% after. This is explained by the fact that in the experimental study, the early charcoal burning zone is centered. The combustion front has not yet reached the side walls. After thirty minutes (30 min) of combustion, all the walls are in contact with the incandescent fuel.
3.3.3. Comparison of Theoretical and Experimental Temperature Profiles of External Wall Temperatures of the Oven
The theoretical and experimental temperature profiles of the bottom wall of the oven are shown in Figure 21.
The temperature profiles also have a first phase of growth and a second phase of decay. During these two phases, the temperature profiles remain similar, therefore indicating that the theoretical model describes the overall behavior of the bottom wall of the coal. The relative error between the theoretical results and the experimental results is 9.68%.
3.3.4. Comparison of the Theoretical and Experimental Temperature Profiles of the Grid
Figure 22 depicts the theoretical and experimental temperature profiles of the grid on which the products to be grilled are placed.
Figure 20. Experimental lateral parietal temperature profiles.
Figure 21. Profiles of external bottom wall theoretical and experimental temperature.
Figure 22. Temporal evolution of the theoretical and experimental temperatures of the grid.
Both profiles also have a first phase of rapid growth and a second phase of decay a little slow. During these two phases, the temperature profiles remain close, which allows us to say that our theoretical model manages to describe the overall thermal behavior of the oven grate. Since the combustion zone is centered at the beginning, the walls are not in contact with the incandescent coal, which explains the small variation of the experimental temperature. The relative error between the experimental and theoretical temperatures is therefore large at first but after 20 minutes the values decrease to about 7.19%.
3.3.5. Comparison of Theoretical and Experimental Temperature Profiles of oven Temperatures
Through Figure 23 and Figure 24, we can observe the experimental and theoretical temperature profiles of the oven.
These two figures show us that all the temperatures follow that of the fuel. With the increase of the fuel’s temperature, the temperatures of the other parts keep increasing as well. The theoretical temperature under the oven remains greater than those of the side walls and those of the grid as in the experimental study. The theoretical temperatures of the different walls as well as that of the grid are similar. These results allowed us to validate our theoretical model of the oven.
Figure 23. Experimental temperature profiles.
3.4. Quantification of the Energy Lost by the Walls and the Corresponding Charcoal Mass
From the temporal profile of the temperatures, the energy lost by the oven walls can be calculated. In this study, energy stored in the walls and energy lost to outside environment were considered lost energy. The amount of energy lost per time unit is represented by the following figure.
Figure 25 shows that the amount of energy lost is very large from the beginning. This is due to the high conductivity of the iron sheet. Charcoal heating value is estimated by P. Khummongkol et al. [33]:
. The amount of energy lost with the oven walls during the two hours of experience is
, which gives about
on 3 kg of burned charcoal. This shows that almost half of the burned charcoal is lost through the barbecue walls.
4. Conclusion
The experimental and numerical study carried out in this work made it possible to understand the thermal behavior of a model of oven commonly used in the sector of grilling in Burkina Faso. From this study, it appears that the oven has low energy efficiency due to the thermophysical properties of the materials used and the poor combustion of charcoal. As a consequence, this results in wasted energy and the discomfort of users exposed to high temperatures (over 200 ˚ C). In addition, the amount of energy lost by the oven walls represents about half of the energy produced by the combustion of charcoal. This study clearly shows that these ovens need to be improved in terms of energy efficiency, which will be the subject of our next work.
Figure 25. Temporal heat lost with barbecue walls.
Acknowledgements
The authors express their deep gratitude to the International Science Program (ISP) for their financial support.
Nomenclature
: Emissivity of the wall
: Thermal conductivity (W∙m−1∙K−1)
: Density (kg∙m−3)
: Kinematic viscosity (m2∙s−1)
: Stefan Boltzmann constant (5.67 × 10−8 m−2∙K−4)
: Dynamic viscosity
e:Characteristic thickness (m)
Ci: Massive heat (kJ/kg∙K)
D: Hydraulic diameter (m)
Fciel: Sky form factor
: Thermal conductance
Gr: Grashof
hc: Convection heat transfer coefficient (W∙m−2∙K−1)
hr: Radiative conductance (W∙m−2∙K−1)
mi: Massa of node I (kg)
Nu: Nusselt number
P: Perimeter
Pr: Prandlt number
Re: Neynold number
S: Surface (m2)
Tamb: Ambient temperature (K)
Tciel: Sky temperature (K)
V: Air velocity (m/s)