1. Introduction
Human health depends on air quality. Emissions, especially particulate matter (PM), contained in the air can penetrate the lungs and blood circulation, which can cause various respiratory and cardiovascular diseases [
1,
2]. Therefore, it is important to search for possibilities to reduce the concentration of PM in the air. The air pollution with particulate matter is mainly caused by traffic, industrial activity, domestic fuel combustion, natural sources, and unspecified sources of human origin [
3,
4]. The largest producer of particulate matter in Central Europe, where Slovakia is also located, is domestic fuel combustion [
4], which is also the interest of this article.
There are also the air quality standards for particulate matter based on the European Union (EU) limits and the air quality guidelines (AQG) set by the World Health Organization (WHO). The limit values for particulate matter with a diameter of 10 μm or less (PM
10) are as follows: (I) 50 μm/m
3 per 1 day given by the EU limit value and also the WHO AQG, (II) 40 μm/m
3 per calendar year given by the limit value, and (III) 20 μm/m
3 per calendar year given by the WHO WQG. The limit values for particulate matter with a diameter of 2.5 μm or less (PM
2.5) are as follows: (I) 25 μm/m
3 per 1 day given only by the WHO AQG, (II) 25 μm/m
3 per calendar year given by the EU limit value, and (III) 10 μm/m
3 per calendar year given by the WHO WQG [
5].
It is important to capture these particles before they escape into the air. With this aim, there is an effort to optimize the combustion process, and also various geometric constructions can be proposed, which are inserted into the heat sources. Their installation and maintenance are easier compared to other types of separators placed into the chimney. The produced fly ash from combustion can also be used for making nanoporous materials, e.g., for capturing carbon dioxide (CO
2), as is stated in the work of Miricioiu and Niculescu [
6].
The optimization of a biomass boiler by an air manifold was investigated by Bianco et al. The authors found out that reducing the entropy meant improving the device performances [
7]. The research of Smith et al. was focused on the optimization of biomass combustion by burner designs [
8]. The main interest of the work of Nussbaumer et al. was the secondary air injection with the aim of the optimization of fluid dynamics for combustion improvement [
9]. Athanasios et al. investigated an optimized wood boiler manufactured by a Greek company, which should belong in Class 3 of the EN303-5 standard [
10]. All these works chose the approach of using CFD simulations, some of which were verified experimentally or, in the case of Nussbaumer et al., using the PIV method. However, none of these works dealt with optimization aimed at capturing particulate matter. Therefore, there is little knowledge regarding particulate matter tracking in the combustion chamber or flue gas tract of the heat source.
In the cases of combustion simulation, realizable k−ε models are used, where k is the turbulence kinetic energy and ε is the dissipation rate of the turbulence kinetic energy. This approach was also chosen by Bianco et al., where it was stated that the advantage of the realizable k−ε model is the modified dissipation rate equation and damping function [
7]. A realizable k-ε model was also used by Silva et al. in their research of industrial furnaces. The authors showed results of fields for temperature, velocity, and concentration of carbon monoxide (CO) and carbon dioxide (CO
2), which should be a tool for the optimization [
11]. In the research of Gómez et al., turbulence was modeled by a realizable k–ε model. They simulated a biomass boiler and analyzed the thermal behavior using exhaust gas recirculation [
12]. A realizable k−ε model with enhanced wall functions was chosen in the work of Zadravec et al. They investigated air staged combustion, which led to lower emissions and higher thermal efficiency [
13]. The standard k−ε model was used in the works of Smith et al. and Athanasios et al.
Measurements of particle velocity can be performed using particle imaging velocimetry (PIV). The laser light sheet illuminates the observed flow, and the instantaneous figures are captured. The figure analysis allows obtaining flow fields [
14]. This method has been used for the observation of the flue gas velocity during solid fuel combustion in several works. Mizeraczyk et al. investigated the velocity field in an electrostatic precipitator by using the PIV method. Their work showed the importance of the secondary flows, which influenced the motion and precipitation of small particles [
15]. The flow field of the pulverized coal flame was observed by the PIV method in the work of Sung and Choi. They focused on the heat release region and the internal recirculation zone [
16]. The combined approach of PIV and particle tracking velocimetry (PTV) was used by Schneider et al. in their work. They also investigated the flow field of solid fuel particles in the combustion chamber [
17]. However, most of the works were focused on the observation of velocity in areas outside the flue gas tract.
Based on this knowledge, it is important to search for possibilities to decrease the particulate matter in the environment. One of the possibilities can be capturing particulate matter in the flue gas tract caused by geometric modification of the heat source, which is also the aim of this article. For this case, a wood stove modified by three baffles placed in the flue gas tract was chosen. The effect of baffles was observed using CFD simulation, while geometric parameters of the first baffle were changed, and was verified using the PIV method. In the other work of Nussbaumer et al., it was stated that CFD is a suitable tool for improving the design of flow during combustion. PIV is a complementary tool, which can be used for validation of flow situations [
18].
2. Materials and Methods
2.1. The CFD Model
The created model, which was used in CFD simulation, is shown in
Figure 1. It represents the flue gas tract of the wood stove. However, the flue gas tract was modified by using three baffles, which were designed to capture particulate matter. The position of the first baffle was changed twice, and therefore the flowing was observed during three different settings. The flue gas tract was designed as a prism with a height of 500 mm, a width of 330 mm, and a thickness of 500 mm. The inlet, through which the particles entered the model, was realized as a rectangular area with a length of 500 mm and a width of 80 mm. The outlet was realized by using a rectangular area with a length of 500 mm and a width of 85 mm. The flue gas tract also contained a second baffle, which was located 170 mm from the sidewall of the model, 400 mm in height, and was divided into two parts at the end. The third baffle was placed 260 mm from the sidewall, and perpendicular to the bottom wall of the model in the shape of an extended letter
M.
The first position of the first baffle used was 80 mm from the wall and an unchangeable height of 400 mm. At the end of the straight part of the baffle was a curvature which directs the flow of flue gases. The model settings were modified by changing the following parameters: (I) width X, (II) height Y, (III) angle α, (IV) angle β, and (V) angle γ. The values of changeable parameters are stated in
Table 1.
The meshes of the observed models are shown in
Figure 2. The shape of mesh elements was mostly hexagonal in all proposed settings. The number of elements, and the average values of the orthogonal quality, the skewness, and the aspect ratio are stated in
Table 2. Based on the data stated in
Table 2, the mesh quality was evaluated as acceptable to the criteria in all proposed settings.
Wall boundary conditions were set up for individual parts of baffles. The CFD model was developed in Ansys software (ANSYS, Inc., Canonsburg, PA, USA), where the CFD solver was Ansys Fluent.
Boundary conditions are stated in
Table 3. The velocity of the flue gas was set up for the inlet and the pressure was set up for the outlet. Firstly, the mass flow in the inlet was evaluated by stoichiometric calculations with a consideration of 15 kW heat power. The simulation was considered with the fuel of wood pellets. Their characteristics were as follows: (I) 49.8% of carbon, (II) 6.4% of hydrogen, (III) 0.3% of nitrogen, (IV) 10% of moisture, (V) 0.8% of ash, and (VI) 15.4 MJ.kg
−1 of calorific value [
19]. The inlet velocity was converted from the determined mass flow using the continuity equation. However, after recalculating the velocity into the flue gas tunnel leading to the model, it was not possible to achieve an experimental velocity of 1.97 m∙s
−1 resulting from the stoichiometric calculations. Therefore, the measurements were performed at the minimum velocity achievable in the experimental set-up, which was 2.5 m∙s
−1. The simulations were performed at this velocity converted to the area in the inlet of the model. They were realized at a temperature of 20 °C, the same as during the experimental measurement, without including the combustion process. The air density was set at 1.18 kg∙m
−3 due to a temperature of 20 °C. Outlet pressure was evaluated at −12 Pa according to STN EN 13240.
The simulation of particulate matter motion in the flue gas was performed using a two-phase model. The airflow was simulated using the Navier–Stokes equations. Lagrangian particle tracking was used for particle motion, which is stated in Equation (1).
where
Fx is an additional acceleration [m∙s
−2],
FD is the drag force [s
−1],
u is the fluid velocity [m∙s
−1],
ρ is the fluid density [kg∙m
−3],
up is the particle velocity [m∙s
−1], and
ρp is the particle density [kg∙m
−3] [
20].
The particle size distribution was performed using the Rosin–Rammler exponential equation, which is given in Equation (2). The size range was divided into several discrete intervals, this interval representing the mean average for which the trajectories were calculated.
where
Yd is a volume fraction of particles with a diameter larger than
d,
is the mean parameter, and
n is the spread parameter [
20]. An experimental particle size distribution had to be performed to solve this equation. It was realized using a vibrating sieve shaker with the following results: (I) 10.8% particles with a diameter less than 20 µm, (II) 31.65% particles with a diameter between 20 and 50 µm, (III) 35.61% particles with a diameter between 50 and 100 µm, (IV) 14.83% particles with a diameter between 100 and 150 µm, (V) 6.68% particles with a diameter between 150 and 250 µm, and (VI) 0.43% particles with a diameter between 250 and 500 µm [
21].
The simulation was realized by using the k–ε realizable model with a scalable wall function. The combustion process was not included in this model. The model was simulated by considering the temperature of 20 °C, which was reached during the experimental measurement.
2.2. PIV Measurement
The CFD model was verified by PIV measurement. The experimental set-up consisted of the model with baffles created based on the CFD simulations, a closed air tunnel to ensure the flow conditions, and the equipment needed for the PIV method. The model of the flue gas tract was inserted into the fiberboard. This model was constructed from transparent polymethyl methacrylate (PMMA) material, which was connected using metal elements. The dimensions of the model were the same as when using the CFD, and the settings of the model were also modified in the same way.
The closed air tunnel consisted of a flow fan, a differential pressure sensor for setting the fan operating velocity, a temperature sensor for controlling the flow conditions, an anemometer for measuring and setting the airflow velocity, a droplet generator that produced and distributed a tracer in the form of oil droplets to the air tunnel, and a logger. Oil droplets represented the particulate matter. The experimental set-up is shown in
Figure 3.
PIV equipment included a charge-coupled device (CCD) camera, a neodymium-doped yttrium aluminum garnet (Nd: YAG) laser, and the Dynamic Studio software for the evaluation measurements. Oil droplets were illuminated by the Nd: YAG laser. The main laser beam had a wavelength of 1064 nm, and the second laser beam with a wavelength of 532 nm was found in the spectrum of green light. The laser beam was transformed into a divergent laser light sheet by light sheet optics. The CCD camera recorded the figures at each pulse. This camera was rotated 90° from the laser plane and also equipped with a filter that let through only green particles. PIV equipment is shown in
Figure 4.
The measurement was realized in a dark room without the influence of other light sources. The first measurement area was located at the top of the model, where the first baffle was curved. These measurements were realized with three different settings of the baffle. The second measurement area was realized in the place of the divide of the second baffle. The measurement realized in this area was only for setting 1. Before the start of the measurements, it was necessary to calibrate the camera and select points A and B. Based on these points, the scale was determined. The velocity of flowing oil droplets in the tunnels reached 2.5 m∙s−1. The camera exposure time was set to 300 µs. Depending on the flow velocity, the time between figures was also set to 450–900 µs. The optimal displacement of the particles in the evaluation area should not be greater than 8 pixels. In this case, the maximum displacement was 6 pixels during the measurement between the two figures. The laser intensity should be adjusted so that the figure quality from the first laser and the second laser is the same. Based on the visibility of the particles and the number of shading and obstacle elements, the intensity of the first laser was set to 5.7–6.2 and the intensity of the second laser 6.0–6.3. Firstly, the size of the evaluated area was 16 × 16 pixels. Then, if no particle was found, the area was gradually enlarged to 64 × 64 pixels. The size of the evaluation area was determined so that all particles were recorded.
The Dynamic Studio software (Dantec Dynamics A/S, Skovlunde, Denmark) used two figures from the measurement, the correlation, a mask, and the removal of erroneous vectors with the use of the range, peak, and moving average. The result is the time middle-average vector velocity field, which is created from all validated instantaneous states of the velocity field in a given set of measurements.
Figure records of PIV measurements were evaluated by mathematical statistics, which were correlations between two consecutive images. The measurement uncertainty is expressed by a correlation coefficient with a range from 0 to 1, where 0 means the match is 0% and it was not possible to identify the particle shift between frames, and 1 represents 100% correlation certainty. The value of the correlation coefficient is influenced, among other things, by the quality of the image recording, poorly lighted measured areas, too large a variance of velocities in the monitored space, and the formation of lost pairs. Lost pairs are caused by the incapability to identify a particle in the interrogation area on two consecutive frames, and this can be caused by the incorrect setting of time between laser pulses, when the particle shift is greater than the maximum allowable value of 8 pixels, large velocity variability in the monitored area, significant 3D flow, etc.
During the experimental measurements, a correlation coefficient was achieved in the range of 0.96 to 0.98 in the main monitored areas, which can be considered as a very accurate measurement. A more significant error was in poorly lighted parts of the monitored areas, where the correlation coefficient achieved a value in the range of 0.82 to 0.87. Poorly lighted areas represented shadows behind the baffles.
4. Conclusions
The reduction in particulate matter during combustion can be achieved by optimizing the flue gas flow. This is important in terms of its impact on human health. This article deals with the optimization of the flue gas flow using three separation baffles located in the flue gas tract. The geometric parameters of the first separation baffle were changed and their influence was investigated using CFD simulations and PIV measurements.
Based on CFD simulations, the results showed that the most effective is setting 1. Setting 2 can capture 1.8% less particles and setting 3 can capture 0.6% less particles than setting 1. In terms of stoichiometric calculations, it would be possible to capture up to 62.3% of the particles. Based on PIV measurements, the results showed the greatest vortexes arising at setting 3, but the weakest flue gas stream at setting 1.
The comparison of PIV measurements and CFD simulations confirmed the supposed character of the turbulent flow with vortexes appearing in the flue gas tract. Inaccuracies in PIV measurements that lead to a decrease in the measured velocity could be due to the lost pairs, shading elements, or local insufficient particle saturation. In CFD simulations, higher velocities were achieved. They were idealized, and the flow of flue gas was unaffected by external influences or material properties. The results, despite the mentioned inaccuracies, are sufficient to follow this problem.
Based on the results, it can be concluded that CFD is a suitable tool for optimizing the flue gas flow during combustion, and PIV can be used for validation of this flow, similar to what has also been mentioned in various other works. To obtain more accurate results, it would be appropriate to include the effect of the combustion temperature, which could be a part of future research.