1. Introduction
In the present energy and environmental crisis, biomass energy could be an interesting alternative solution. It presents many economic and environmental benefits such as reducing waste, worldwide availability, and a neutral carbon balance [
1]. Biomass combustion already provides a significant amount of world energy consumption and is the first supplier of renewable energy [
2], but its use should be increased in the future. Furthermore, biomass combustion is known to cause major emission of multiple gaseous pollutants, such as carbon monoxide (CO), nitrogen oxides (NOx), and polycyclic aromatic hydrocarbons (PAHs) [
3], resulting in incomplete combustion. In addition, biomass combustion produces high amounts of particulate matter (PM). These particles could be divided into three ranges: coarse (d > 1 µm), fine (0.1 µm < d < 1 µm) [
4], and ultrafine (d < 0.1 µm). The last two are known to be emitted by biomass combustion and are studied in this work. Moreover, there are three kinds of fine particles, with their own way of forming [
5]:
Soot, formed in the flame by complex reactions of surface, coagulation, and agglomeration of gaseous hydrocarbons, with PAHs being precursor for soot;
Fly ashes, which are volatile metals (Zn, Pb, Ni, Cd, Cu, As, etc.) encased and unburnt in wood;
Organic materials, issued from the decomposition of biomass during the heating process.
All these pollutant emissions have a significant effect on air quality and human health [
6], with biomass combustion being one of the main sources of air pollution in Europe [
7]. Many investigations have shown that the boiler characteristics, the nature of biomass, and the combustion operating conditions have a great influence on the formation of pollutants [
8,
9,
10].
Some numerical studies have been carried out on particle-laden flows, such as boiler chimneys. Accurate simulation of the deposition and distribution of fine particles in turbulent-laden flows requires correct modeling of both the carrier flow and the dispersion of micro and nanoparticles. One of the key challenges lies in capturing the complex interplay between turbulent fluid motion and particle dispersion.
Modeling of particle dispersion in particle-laden flows typically falls into two main categories: Eulerian–Eulerian and Eulerian–Lagrangian approaches [
11,
12,
13]. In the Eulerian–Eulerian approach, particles are considered to be a continuous phase, and conservation equations, such as continuity, momentum, and energy, are solved similarly to the carrier gas flow field. Hence, the motions of individual particles are not considered, and the formulation of boundary conditions for the second phase is challenging [
14]. This approach is well-suited for denser suspensions where particle–particle interactions are significant and the feedback from particles on the flow cannot be ignored. Nevertheless, accurately specifying inter-phase exchange rates presents a challenge. On the other hand, the Eulerian–Lagrangian approach treats particles as discrete entities dispersed in the continuous carrier flow phase. It is preferred over the Eulerian–Eulerian method because it accounts for the discrete nature of particles and thus enables a more accurate management of turbulent diffusion. Particle trajectories are determined by integrating Newton’s second law, accounting for relevant forces such as gravity, drag, lift, etc. In the present work, Lagrangian particle tracking methodology is employed, assuming the dispersed phase does not significantly affect the continuous flow field (one-way coupling).
Presently, three primary approaches exist for simulating turbulent flows. Direct numerical simulation (DNS) resolves all spatial and temporal turbulent scales, which makes it the most accurate approach, although its application in industrial contexts is very challenging due to its prohibitive computational cost. On the other hand, large eddy simulation (LES) resolves turbulent scales larger than the grid cell size, while subgrid-scale fluctuations are modeled. Despite being less computationally intensive than DNS, LES still poses challenges for simulating complex flows in high Reynolds number environments due to the drastic increase in computational costs with the Reynolds number. In the context of industrial applications, the most commonly used alternative to DNS/LES is the Reynolds-averaged Navier–Stokes (RANS) approach, which models the entire turbulence spectrum. In this study, the
turbulence closure model is used for its ability to represent flow anisotropy near the wall [
14].
In contrast to the DNS/LES methods, the RANS approach evaluates only the ensemble-averaged flow fields. Consequently, the RANS equations lack the instantaneous fluid velocity vector at the particle position. It becomes essential to model instantaneous fluctuations in order to consider the complex coupling between turbulent flow motion and particle dispersion. This approach, inspired by stochastic models, aims to simulate turbulence fluctuations by employing additional equations in order to calculate instantaneous turbulent velocities using local computed quantities such as distance to the wall, mean turbulent kinetic energy, and Eulerian time scale. Random walk models are significant examples of this approach, known for their moderate computational costs and ease of implementation. The DRW model [
15], also known as the Eddy interaction model (ETM), emerges as a robust approach, originally proposed for homogeneous turbulent flows.
The aim of this work is to experimentally investigate the gaseous and PM emissions in the chimney of a domestic-scale boiler with a nominal power of 10 kW, during the combustion of natural wood pellets. This boiler is equipped with an underfeed burner, instrumented with several sensors (thermocouples, flowmeters, etc.) located in strategic areas and several points of measurement. These sensors allow for a complete characterization of the laden flow in terms of temperature and velocity, a measurement of gaseous pollutants, and a study of the evolution of the emitted PM characteristics along the chimney and at the outlet. Additionally, a numerical model of the carrier flow and particle transportation was developed and validated by experimental data.
3. Experimental Results
All the experiments were conducted during the steady phase of the boiler. To ensure the reproducibility of the system, four tests of about 1 h of recording were performed to monitor the evolution of the operating parameters of the boiler.
In
Table 3, the results show that the boiler functioning is repeatable, with maximum variation to the mean, over the four tests, around 5%.
3.1. Measurement of Gaseous Pollutant
For the gaseous emissions, eight tests were recorded at position P2. The results are shown in
Table 4.
The main result of these measurements is that, even if the boiler functioning is stable and reproducible, this is not the case for pollutant emissions, especially for gaseous pollutants, with a maximal variation to the mean of 41% around a mean value of 465.3 ppm for the CO mole fraction, and a variation of 52% around 32.6 ppm for the THC mole fraction. The CO2 mole fraction also shows a significant variation but around a low value of concentration. Only the O2 levels stay rather stable over the different tests.
Comparisons are made with French standards, and to other similar studies [
34,
35] in
Table 5, showing that the results are in the order of magnitude of the literature review for other studies using wood pellets.
Moreover,
Figure 4 shows an example of the temporal evolution of gaseous concentration levels. It is important to note that CO and THC emissions are unstable over time but are temporally coherent.
These temporal fluctuations could be explained by the fluctuation of the flame, visible through the optical window of the boiler, even during the steady phase, due to the geometry of the burner: more than 20 air intakes, distributed around the burner, leading to an asymmetric distribution of the air through the flame.
3.2. Evolution of PM Levels along the Chimney
One of the aims of this study was to highlight the evolution of PM granulometry and concentration along the chimney of the boiler. Six tests were performed at several positions (P1 to P4 in
Figure 2), to check reproducibility. The results are presented in
Figure 5 and
Table 6.
Except for position P2, where measurements could be disturbed by the T-shape of the chimney (noticeable on
Figure 2), the data show, in one side, a decrease in the concentration, from 5.1 × 10
7 particles/cm
3 at P1 to 3.9 × 10
7 at P4, and in another side a slight increase in the median diameter, from 50.7 to 58.0 nm.
These results endorse the theory of a deposition of the emitted PM, which contribute to the decrease in the concentration, but also of an agglomeration of these particles which could explain the increase in the median diameter.
Two remarks can be made from these results:
There is only a small decrease in the PM concentration between P1 and P3, with a great increase in the median diameter. Moreover,
Figure 5a shows a decrease in the concentration of the smaller particles (with diameter < 80 nm) and an increase in the concentration for particles with higher diameter. It may be presumed that the predominant phenomenon between positions P1 and P3 is the agglomeration of particles.
Between the positions P3 and P4, the median diameter shows only a slight evolution, but there is an important decrease in the particles’ concentration. It could be assumed that the predominant phenomenon between position P3 and P4 is the particles’ deposition.
To investigate furthermore these experimental results, a CFD modeling of the particles’ transportation in the chimney was performed, particularly to confirm the hypothesis of the deposition.
4. Numerical Results and Comparison with Experimental Results
Before analyzing the particle distribution, validation work was performed on the carrying flow. The numerical and experimental sampling of the velocity profiles, in
Figure 6, shows a good agreement, at position P3, where the flow is strongly disturbed by the pipe bend.
The temperature is obtained along the chimney as shown in
Figure 7. The temperature of the fluid is imposed at location P1 and is compared to experimental data at locations P2 and P5.
Table 7 gives the comparison of the measured values and the simulated temperatures. At position P2, there is a difference of approximately 4% between the mean measured temperature and the numerical ones. At position P5, the temperature is measured with a PT100 sensor with a higher accuracy. The difference between the numerical value and the experimental value is less than 1 °C. These results give confidence in the prediction of the temperature field, then on the thermophoresis effect, considered in the particle transportation model.
In
Figure 8, the particle concentration obtained with the CFD simulation are plotted according to the particle’s diameter for positions P1, P3, and P4. The concentration imposed in the numerical model at the location P1 is taken from the experimental measurement. A Gaussian distribution is used on the most encountered diameters [
5]. The range of diameter, [6–124 nm], has been chosen according to experimental measurement. Each parcel of the DPM model contains 1 million of particles.
For location P3, the particle concentration is very similar to the imposed condition at position P1. A significant difference is observed at location P4, with a decrease in the particle number on all the ranges of diameters.
In
Figure 9, the experimental and numerical results are compared at locations P3 and P4. For both positions, the numerical results predict quite well the measurements. The maximum concentration is in the same range around 6 × 10
6 particles/cm
3 (P3) and 5 × 10
6 particles/cm
3 (P4). However, the model overestimates the particle numbers in the range 20–50 nm and underestimates the particle numbers in the range 60–120 nm. One consistent assumption to explain these differences can be the agglomeration of particles along the chimney. The current CFD model does not consider the agglomeration of ultrafine particles. Nevertheless, soot particles are known to agglomerate. The size of the particles may change along the chimney. For the experimental results, it has been highlighted that the median diameter is growing along the chimney.
Finally, the total number of particles are compared in
Figure 10. The total numbers of particles at position P1 are equal to the mean measured value, as this is the imposed condition in the chimney inlet. At locations P3 and P4, the same tendency is observed regarding the decrease. At location P3, the simulated total number decreases about 3.8% compared to 6.8% for the experimental measurement. At this point, there is a non-negligible dispersion of the experimental data. At location P4, the simulated and experimental total numbers of particles are closer. A decrease of 20.9% is observed for the numerical simulation between P1 and P4. For the experimental data, the decrease is about 23.3%.
5. Conclusions
The paper deals with the characterization of gaseous and solid pollutants emitted from a 10 kW wood pellet boiler and with the prediction of particle matter transport in the chimney. One of the aims of this work was to predict the modification of the granulometry of the particles along the chimney and at the outlet. In this work, advanced measurement systems (EEPS for PM, spectroscopy and FID for gases) are applied on a domestic-scale installation, instrumented with several sensors. Moreover, numerical CFD simulations were applied on the geometry of the chimney, considering heat transfer phenomena, using a turbulence closure model and a DRW turbulent particles dispersion model. These simulations take into account the operating conditions of the boiler. The novelty of this work is to study experimentally and numerically a domestic biomass boiler in real operating conditions and to compare the measured and modeled PM emissions.
The obtained results show that the following:
The emissions of the main gaseous pollutants (CO and THC) are variable over the different tests, even if the boiler operating parameters stay stable. However, the concentration levels of CO and THC versus time show a really similar evolution.
The PM level evolves significantly on the different measured locations: the concentration decreases along the chimney while the median diameter increases. The evolution of the median diameter can be associated with a particle agglomeration, as biomass soot are known to form agglomeration of nanoparticles.
Additionally, the decrease in the total number of particles along the chimney can be explained by the deposition. Some deposits are observed inside, especially in the horizontal part.
A comparison with a CFD model gives similar results on the deposition regarding the decrease in the total number of particles along the chimney.
However, the particle–particle interaction is not taken into account in the used model. The particles’ adhesion and agglomeration are not simulated, and some differences are observed for the granulometric distribution. Integrating these behaviors will be the next aim to improve the forecast. This will bring richer information regarding the numerical simulation to understand the causes of the evolution of the particle’s concentration along the chimney.