The Analysis of the Possibility of Feeding a Liquid Catalyst to a Coal Dust Channel

This article presents the concept of a method of improving the dynamics of combustion in boilers operating in power plants, cogeneration plants, and heating plants by introducing a catalyst that is fed with a carrier in the form of droplets. Thanks to the proposed method, a greater degree of fuel burnout can be obtained, which, in turn, results in lower energy consumption in the case of producing the same amount of power. The parameters of the emitted exhaust gases and ash are also improved. The method described in the article involves the adding of a catalyst to the dust pipe of the boiler, which improves the combustion parameters. The catalyst was implemented using a sprayer/nebulizer. In order to obtain the correct flow parameters, the sprayer was modeled using CFD calculations. The calculations include trajectories, velocities and concentrations with regards to various flow parameters. Particular attention should be paid to the model of the evaporation of moving droplets. The results of these calculations enable the parameters that guarantee that the catalyst reaches the dust channel outlet in the desired form to be assessed. The analysis is an introduction to experimental research that is carried out on a medium and large scale.


Introduction
Due to the rising share of energy from renewable sources, both in the case of domestic and foreign production capacities, there is an increasing need to improve the operating parameters of classic power units. It is worth emphasizing that soot from combustion processes is a significant pollutant of the Earth's atmosphere, and it is the second-most important factor influencing the climate, right after carbon dioxide (CO 2 ) [1]. These and other environmental problems have prompted the industry to implement both newer solutions and more expensive and complex exhaust gas treatment technologies in combustion processes [2][3][4]. The general energy policy and global trends are focused on achieving three basic goals, i.e., minimizing the emission of harmful substances to the environment, keeping energy prices at the lowest possible level, and ensuring the security of the energy supply in the short and long term [5,6]. On the other hand, there is also an aim to increase the efficiency of electricity production. These goals are achieved thanks to the possibility of using modern technologies and materials in individual elements of conventional power plants. One such element is the boiler. Increasing requirements for the operation of boilers are forced by the variable loads of energy systems, in which the unstable natures of work, wind, and solar power plants have an increasing share. This results in the necessity to increase the dynamics (flexibility) of coal-fired power plants. In this case, the ability to transfer from basic to regulating work and, also, the speed and range of changes of the loads that are associated with it can be seen to be important.
There are many possibilities of modifying the course of the combustion process, which would allow the energy efficiency of boilers to be increased [7]. One of the most popular solutions is the use of catalysts [8][9][10]. It is not only the choice of the type of modifier that is the subject of extensive considerations by technologists and scientists but, also, the way it is delivered to the right place and in the right form. In the literature, it can be found that the course of the combustion process also depends on factors such as: environmental parameters, the properties of the feed fuel (e.g., gas), or the geometry of the injector and dust channel [2,3,[11][12][13]. Dust channels are not only widely used in combustion processes but, also, in many industries related to, among others, mining and the ventilation of rooms, industrial buildings, and utility buildings [14][15][16][17][18][19]. For example, the authors of paper [2] carried out an analysis of the combustion process with regards to various environmental parameters, the properties of feed fuel, and the various geometric parameters of injectors. As a result, channel fuel injectors have been proven to be effective in reducing soot levels by 35-100% at a wide range of operating conditions and geometries. Moreover, guidance on the geometric parameters that are most important for the purpose of improving the process efficiency were also provided. In turn, the authors of paper [20], in addition to the analysis of the above-mentioned factors, also attempted to better understand the process and the phenomenon of soot formation during combustion processes. These studies aimed to reduce the amount of soot emissions. In the process of coal combustion in power boilers, various additives are used to modify this process. They are designed to prevent the accumulation of soot and other deposits. The appropriate use of combustion process modifiers can significantly contribute to the improvement of the quality parameters of exhaust emissions. In the process of coal combustion in power boilers, various additives are used to modify this process. They are meant to prevent the accumulation of soot and other deposits. The correct use of combustion process modifiers can significantly contribute to improving the quality parameters of exhaust emissions. Therefore, it is important to select the chemical composition of the fuel additive, as well as the method of its administration. Many different fuel additive concepts have been proposed, with the most popular additives being copper compounds and sodium chloride. Fuel additives enable the load associated with coal combustion to be decreased and, therefore, the amount of fuel needed to produce a unit amount of energy to also be reduced [10,21,22].
The precise dosing of catalysts to the flux that is fed to the fuel combustion chamber is still a difficult issue, and many variables that influence the dosing accuracy need to be considered. The parameters of the catalyst transporting flux, the assumed quantity and form of the catalyst, its concentration, and the geometry of the system should definitely be analyzed.
The main goal of the study was to develop a system that would allow for the precise dosing of the catalyst into the combustion chamber and, also, to analyze the effectiveness of the developed concept by performing numerical calculations. The method described in this article involves the intensification of the combustion of solid fuels in a boiler by using alkyl alcohol and Al-Ni catalyst as combustion promoters. Work [23] described in detail the operating principle of the catalyst. There is no chemical reaction of the fuel with the catalyst, but at the appropriate temperature, isopropanol in the presence of the catalyst is hydrogenated, and free radicals are formed, and the combustion reaction is more efficient and faster. It should be noted that the catalyst only works when there is alcohol in the droplets. If it evaporates, the desired effect will not be achieved.

Materials and Methods
The main assumption of the method, i.e., delivering a droplet of alcohol with a suspended catalyst to the pyrolysis zone (the droplets cannot evaporate), was tested using simulations (Ansys Fluent, Canonsburg, PA, USA). The obtained results will be verified on a specially constructed measuring stand. If the droplets in the 1.0-1.5-m section of the pipeline do not evaporate, the proposed method can be successfully applied. It is assumed that the catalyst is fed using a nebulizer. Its system is designed in such a way that the size and quantity of the catalyst can be controlled. In order to achieve the correct parameters of the flow, it was modeled with the use of CFD calculations in Ansys Fluent 2021R1 software. These calculations include trajectories, velocities, and concentrations for various flow parameters. Due to this, the obtained spectrum of results can be adjusted to a specific case. Particular attention should be paid to the evaporation model of moving particles. This information allows the delivery of the appropriate amount of droplets (with an appropriate size) that contain the catalyst to be programmed.

Evaporation-Condensation Model
The evaporation-condensation model is a mechanistic model [24] that has a physical basis and which can be used in the case of mixed and Eulerian multiphase models. In this model, the mass transfer is described by the following formulas: in the case of the evaporation process, i.e., when T > T sat : whereas, in the case of the condensation process, i.e., when T < T sat : .
where . m e→v is the mass transfer flux, i.e., from the liquid phase to the vapor phase, coeff is a coefficient that is different when analyzing different problems and which can be interpreted as a relaxation time, α is a volume fraction, and ρ is the density. The source term for the energy equation can be derived by multiplying the mass transfer flux by the latent heat.
In turn, the Hertz-Knudsen formula describes the evaporation-condensation flux on the basis of the kinetic theory of gases: where M is mass of molecule, P is the pressure, T is the temperature, R is the gas constant, β is the so-called accommodation coefficient (which refers to the amount of vapor molecules that go into the liquid surface, which are then adsorbed by this surface), and P* is the vapor pressure. The abberviation sat refers to the saturation condition. The Clapeyron-Clausius equation shows the relation between the pressure and the temperature in the case of the saturation process: where v g is the inverse of density of the gas, v l is the inverse of density of the liquid (volume per mass unit), and L is the latent heat (J/kg). Based on this differential equation, a variation of temperature with regards to the variation of pressure, which is close to the saturation condition, can be obtained. When P * and T * are close to the saturation condition (Figure 1), the Clapeyron-Clausius equation takes the following form: Taking the above into consideration, the Hertz-Knudsen equation can be written as follows: where β is determined using the accomodation coefficient and the physical characteristics of the gas. In the case of conditions that are close to equlibrium, β is equal to 1.0. In the Eulerian and mixed multiphase models, the flow is assumed to be dispersed. In the case when all vapor bubbles have the same diameter, the interfacial area density is described by the formula where Vcell is the cell volume, and the phase source term takes the following form: In this case, coeff, which is inversed to the relaxation time (1/s), can be determined as follows: This leads to the final expression for the vaporization, which is defined in Equation (2). It can be treated implicitly as a source term in the phase conservation equation. A similar expression can be obtained in the case of the condensation process. However, small droplets in a continuous vapor phase need to be considered, even if the primary phase is a liquid.
It should be highlighted that coefficient coeff should theoretically be different for the condensation and evaporation expressions. Moreover, the flat interface, the dispersed regime of the flux of droplets with a constant diameter, and known β were the main assumptions for the obtained theoretical expression. Taking the above into consideration, the Hertz-Knudsen equation can be written as follows: where β is determined using the accomodation coefficient and the physical characteristics of the gas. In the case of conditions that are close to equlibrium, β is equal to 1.0.
In the Eulerian and mixed multiphase models, the flow is assumed to be dispersed. In the case when all vapor bubbles have the same diameter, the interfacial area density is described by the formula where V cell is the cell volume, and the phase source term takes the following form: In this case, coeff, which is inversed to the relaxation time (1/s), can be determined as follows: This leads to the final expression for the vaporization, which is defined in Equation (2). It can be treated implicitly as a source term in the phase conservation equation. A similar expression can be obtained in the case of the condensation process. However, small droplets in a continuous vapor phase need to be considered, even if the primary phase is a liquid.
It should be highlighted that coefficient coeff should theoretically be different for the condensation and evaporation expressions. Moreover, the flat interface, the dispersed regime of the flux of droplets with a constant diameter, and known β were the main assumptions for the obtained theoretical expression.
Due to the fact that the diameter of a bubble and the accommodation coefficient are usually not fully known, the coefficient coeff can be adjusted in order to correlate with the Energies 2021, 14, 8521 5 of 14 experimental data. By default, this coefficient, in the case of evaporation and condensation, is equal to 0.1.

Geometry of a Dust Channel
Numerical calculations were carried out for a model of geometry corresponding to that of the actual apparatus ( Figure 2) in which the experimental studies were conducted. The stand consists of dust channels that are fixed to the OP-430 boiler No. 15 at the Siekierki CHP Plant (Warsaw) in Poland. The temperature in the dust channel is 110-115 • C (>100 • C). The dust channel is a cylindrical stainless-steel pipe with a diameter of 300 mm and a length of 2.5 and 6.0 m. The tests were carried out for the isopropanol evaporation process (100% of isopropanol as the catalyst's main system) with the following parameters: isopropanol density of 785 kg/m 3 , specific heat of 2570 J/kg × K, heat of vaporization of 700 kJ/kg, and boiling point of 355 K. The geometry and numerical mesh of the dust channel created in the ANSYS Fluent 2021R1 program are shown in Figure 3. The numerical mesh applied in the computations is poly-hexcore, with the boundary layer being applied at all the domain walls. The mesh contains around 324 k numerical cells inside a volume of 0.623 m 3 (i.e., volume of the computational domain), and therefore, the average cell size (treated as a cube) of the whole domain is 12.4 mm. The mesh is presented in Figure 4.
Due to the fact that the diameter of a bubble and the accommodation coefficient ar usually not fully known, the coefficient coeff can be adjusted in order to correlate with th experimental data. By default, this coefficient, in the case of evaporation and condensa tion, is equal to 0.1.

Geometry of a Dust Channel
Numerical calculations were carried out for a model of geometry corresponding t that of the actual apparatus ( Figure 2) in which the experimental studies were conducted The stand consists of dust channels that are fixed to the OP-430 boiler No. 15 at the Siek ierki CHP Plant (Warsaw) in Poland. The temperature in the dust channel is 110-115 °C (>100 °C). The dust channel is a cylindrical stainless-steel pipe with a diameter of 300 mm and a length of 2.5 and 6.0 m. The tests were carried out for the isopropanol evaporatio process (100% of isopropanol as the catalyst's main system) with the following parame ters: isopropanol density of 785 kg/m 3 , specific heat of 2570 J/kg × K, heat of vaporizatio of 700 kJ/kg, and boiling point of 355 K. The geometry and numerical mesh of the dus channel created in the ANSYS Fluent 2021R1 program are shown in Figure 3. The numer ical mesh applied in the computations is poly-hexcore, with the boundary layer being ap plied at all the domain walls. The mesh contains around 324 k numerical cells inside volume of 0.623 m 3 (i.e., volume of the computational domain), and therefore, the averag cell size (treated as a cube) of the whole domain is 12.4 mm. The mesh is presented i Figure 4.           the pulverized coal before the bent section was leveled and had a homogeneous dust concentration of approx. 250 g/m 3 . Downstream of the bent section, after changing its flow direction, the pulverized coal mainly accumulated near the upper part of the dust channel (Figure 5b), where the dust concentration was above 1000 g/m 3 . Then, some of the pulverized coal flowed along the wall of the dust channel, and its flow was generally only observed around the wall along the entire perimeter (high dust concentrations above 1000 g/m 3 ). In turn, the transporting gas flowed through the center of the dust channel (dust concentration below 50 g/m 3 - Figure 5d). The 30-µm-diameter droplets of 100% isopropanol spray delivered to the top of the dust channel completely evaporated before the outlet opening (Figure 5c). before the bent section was leveled and had a homogeneous dust concentration of approx. 250 g/m 3 . Downstream of the bent section, after changing its flow direction, the pulverized coal mainly accumulated near the upper part of the dust channel (Figure 5b), where the dust concentration was above 1000 g/m 3 . Then, some of the pulverized coal flowed along the wall of the dust channel, and its flow was generally only observed around the wall along the entire perimeter (high dust concentrations above 1000 g/m 3 ). In turn, the transporting gas flowed through the center of the dust channel (dust concentration below 50 g/m 3 - Figure 5d). The 30-µm-diameter droplets of 100% isopropanol spray delivered to the top of the dust channel completely evaporated before the outlet opening (Figure 5c).

Droplet Motion and Evaporation Analysis
In the current study, the Euler-Lagrange approach is applied to simulate the evolution of evaporating droplets and pulverized coal, i.e., the gas phase is treated in the Eulerian manner, while the droplets and pulverized coal particles are treated as Lagrange discrete particles moving in the gas phase. Both continuous and discrete phases interact with each other by exchange of mass, momentum, and energy. The mass, momentum, and energy equations that govern the droplet evolution are presented below. The mass conservation equation of the droplet is defined as [25,26]: where is the droplet mass, and is the droplet evaporation rate defined as follows: in which is the droplet surface area, is the gas density, and is the Spalding mass number defined as: where is the vapor mass fraction at the surface (s) and in the bulk gas (∞). The mass transfer coefficient is calculated from the Sherwood number correlation: in which is the droplet diameter, D is the diffusion coefficient of the vapor in the bulk gas, Sc is the Schmidt number, and is the droplet Reynolds number. The droplet momentum equation is defined as follows [27]: where is the droplet velocity, is the fluid velocity, and is the standard acceleration due to gravity, is the droplet relaxation time, is the droplet density, is the fluid viscosity,

Droplet Motion and Evaporation Analysis
In the current study, the Euler-Lagrange approach is applied to simulate the evolution of evaporating droplets and pulverized coal, i.e., the gas phase is treated in the Eulerian manner, while the droplets and pulverized coal particles are treated as Lagrange discrete particles moving in the gas phase. Both continuous and discrete phases interact with each other by exchange of mass, momentum, and energy. The mass, momentum, and energy equations that govern the droplet evolution are presented below. The mass conservation equation of the droplet is defined as [25,26]: where m d is the droplet mass, and r d is the droplet evaporation rate defined as follows: in which A d is the droplet surface area, ρ is the gas density, and B M is the Spalding mass number defined as: where Y is the vapor mass fraction at the surface (s) and in the bulk gas (∞). The mass transfer coefficient k c is calculated from the Sherwood number correlation: in which d p is the droplet diameter, D is the diffusion coefficient of the vapor in the bulk gas, Sc is the Schmidt number, and Re d is the droplet Reynolds number. The droplet momentum equation is defined as follows [27]: where u p is the droplet velocity, u f is the fluid velocity, and g is the standard acceleration due to gravity, is the droplet relaxation time, ρ d is the droplet density, µ f is the fluid viscosity, is the normalized drag coefficient, and C D is the drag coefficient. Finally, the droplet temperature T d is updated according to the heat balance, which takes into account the convective and latent heat transfer between the droplet and the continuous phase [25,26]: where c p is the droplet heat capacity, T f is the gas temperature, h f g is the latent heat of the droplet, and h is the convective heat transfer coefficient that is calculated with the modified Nu number as follows [13]: (18) where k is the thermal conductivity of the gas phase, Pr is the Prandtl number of the gas phase, and B T is the Spalding heat transfer number defined as: where c pv is the specific heat of the droplet vapor, and . q d is the heat transferred to the droplet.
Afterwards, using the model presented above, the process of vaporization of the isopropanol droplets of selected diameters during the flow through the dust channel was modeled for different reagent injection conditions. The boundary conditions applied in the numerical simulations are as follows: The next step was to compare the trajectory of 20-µm droplets for isopropanol, water, and the N-C solution. The following parameters were considered for the water: ρ = 994 kg/m 3 , cp = 4200 J/kgK, b = 2382 J/g, and Tp = 100 • C and, for the N-C solution: ρ = 970 kg/m 3 , cp = 3874 J/kgK, b = 2060 J/g, and Tp = 83 • C. The concentration of the isopropanol in the N-C solution was 20%, with the remaining 80% being water. Figure 6 shows the effect of the droplet diameter and the type of solution on the path, followed by the reagent aerosol until it disappears completely (evaporation). The reagent injection took place at a high speed of 120 m/s. Due to the physical properties, for droplets with a composition of 100% isopropanol, the evaporation process is fast and ends before the dust channel's outlet, even in the case of large droplets with a diameter of 30 µm (Figure 6a). When comparing 20 and 30-µm isopropanol droplets, it can be seen that the 20-µm droplet is able to reach a distance of approx. 0.5 m (Figure 6b), while the 30-µm droplet can reach a distance of approx. 0.85 m. For the pure water aerosol, complete evaporation does not occur for the droplets of 20 µm (Figure 6c). Therefore, it should be assumed that, in the case of a certain concentration of the reagent, the vaporization of the aerosol will occur exactly in the desired area, i.e., near the outlet of the dust channel. An appropriate reagent concentration is provided by the N-C solution. This is due to the fact that, for the N-C solution (for high-speed injection conditions), evaporation of the aerosol occurs near the outlet in the case of the droplet with a diameter of 20 µm (Figure 6d). However, in the case of lower injection velocities, which are usually used in practice, the diameter of the droplets, which ensures evaporation in the area of the dust channel's outlet, will be greater than 20 µm. Therefore, as can be seen from further analysis of the obtained images, in order for the droplet to reach the end of the dust channel, it should have a diameter greater than 30 µm.
Energies 2021, 14, x FOR PEER REVIEW 10 of 14 case of lower injection velocities, which are usually used in practice, the diameter of the droplets, which ensures evaporation in the area of the dust channel's outlet, will be greater than 20 µ m. Therefore, as can be seen from further analysis of the obtained images, in order for the droplet to reach the end of the dust channel, it should have a diameter greater than 30 µ m.
(a) Successive simulations show the trajectories of the N-C solution droplets for the selected diameters (from 6 to 30 μm) for three different inlet velocities (10,20, and 30 m/s). Figure 7 shows the effect of changing the droplet's diameter during the flow in the dust channel under selected conditions. Droplets with a size of 20 μm evaporate before the dust channel's outlet. This is in contrast to the droplets with a size of 30 μm, which, for the tested injection speeds, are completely vaporized at the desired distance from the injection point, i.e., near the dust channel's outlet. A droplet size of 20 μm reaches a distance of about 0.85 m. The results show that droplets smaller than 30 μm evaporate, and therefore, for further experimental work, a droplet diameter equal to 30 μm will be considered.  Figure 7 shows the effect of changing the droplet's diameter during the flow in the dust channel under selected conditions. Droplets with a size of 20 µm evaporate before the dust channel's outlet. This is in contrast to the droplets with a size of 30 µm, which, for the tested injection speeds, are completely vaporized at the desired distance from the injection point, i.e., near the dust channel's outlet. A droplet size of 20 µm reaches a distance of about 0.85 m. The results show that droplets smaller than 30 µm evaporate, and therefore, for further experimental work, a droplet diameter equal to 30 µm will be considered. Increasing the inlet velocity from 10 to 30 m/s contributes to the "injection" of the droplet into the flowing flux to a greater depth. Therefore, from the simulations made using Ansys Fluent, it can be seen that the droplet size that is needed to deliver the catalyst (carried by the liquid droplets) to the dust channel's outlet is, in practice, between 30 and 50 µ m (>30 µ m).
By analyzing the obtained results, the velocities and droplet sizes required for the droplets to reach of the end of the dust channel can be indicated. A very important component of the performed simulations is the modeling of the droplets' evaporation. Droplet evaporation for the same diameters and for two different velocities (30 and 120 m/s) was compared for the setting angle of the inlet pipe of the atomizer, which was equal to 45° ( Figure 8). When compared to an angle of 30°, the droplets, at an angle of 45°, get deep inside the dust channel; this means that the catalyst is distributed in a larger volume of the dust channel. In the case of the setting angle of 30°, the droplets reach a distance of approx. 0.12 m from the axis of the dust channel. Therefore, a change of the angle to 45° benefits the delivery process of the droplets. Figure 8b, which shows the trajectory of a droplet with a diameter of 30 µ m, proves that such a droplet will reach the end of the dust channel, effectively transferring the catalyst. Therefore, from the simulations made using Ansys Fluent, it can be seen that the droplet size that is needed to deliver the catalyst (carried by the liquid droplets) to the dust channel's outlet is, in practice, between 30 and 50 µm (>30 µm).
By analyzing the obtained results, the velocities and droplet sizes required for the droplets to reach of the end of the dust channel can be indicated. A very important component of the performed simulations is the modeling of the droplets' evaporation. Droplet evaporation for the same diameters and for two different velocities (30 and 120 m/s) was compared for the setting angle of the inlet pipe of the atomizer, which was equal to 45 • (Figure 8). When compared to an angle of 30 • , the droplets, at an angle of 45 • , get deep inside the dust channel; this means that the catalyst is distributed in a larger volume of the dust channel. In the case of the setting angle of 30 • , the droplets reach a distance of approx. 0.12 m from the axis of the dust channel. Therefore, a change of the angle to 45 • benefits the delivery process of the droplets. Figure 8b, which shows the trajectory of a droplet with a diameter of 30 µm, proves that such a droplet will reach the end of the dust channel, effectively transferring the catalyst. The obtained results show which droplets will not evaporate before reaching the end of the dust channel. The modeling of the above-described phenomena allows for the change of inlet parameters of the applied dispersion (adaptation to various types of boilers) and, thus, the selection of appropriate velocities, droplet size range, and the angle of inclination of the pipe (applicator) with the catalyst. In future research work, the conducted simulations will be verified with tests on an appropriate laboratory stand.

Conclusions
As a result of using the simulation of the real model, a very wide spectrum of results was obtained. The desired droplet size was determined, i.e., the critical (minimum) diameter was determined for the given parameters of the dust channel's operation. The possibilities of adapting the created model to a specific case were analyzed. The paper presented selected calculation results that confirmed the correctness of the made assumptions. Their analysis gives appropriate guidelines for the construction of a test stand. A system that allows for the precise dosing of the catalyst was presented, and the effectiveness of the developed method was analyzed. The method involved the intensification of the combustion of solid fuels in a boiler with the use of alkyl alcohol and Al-Ni (as combustion promoters). Due to the proposed method, it will be possible to obtain a higher degree of fuel burnout, which, in turn, will result in a lower energy consumption in the case of producing the same amount of power and, also, in an increase of the flame's stability. The parameters of the emitted exhaust fumes and ash will also improve. Based on the conducted process simulations, it was shown that: • in order to supply the catalyst to the system, droplets with a diameter of at least 30 µ m should be generated; • an angle of 45° for the inlet pipe of the isopropanol feeding atomizer is beneficial; this allows the catalyst flux to be introduced into the fuel flux.
The study clearly shows that additional factors must be considered when designing the process of introducing a catalyst into a fuel system. The atomizing device should be carefully matched with the liquid properties and other process parameters. The performed simulations are currently being verified on experimental stands. The research should be continued by extending it with tests on dust channels of various geometries, The obtained results show which droplets will not evaporate before reaching the end of the dust channel. The modeling of the above-described phenomena allows for the change of inlet parameters of the applied dispersion (adaptation to various types of boilers) and, thus, the selection of appropriate velocities, droplet size range, and the angle of inclination of the pipe (applicator) with the catalyst. In future research work, the conducted simulations will be verified with tests on an appropriate laboratory stand.

Conclusions
As a result of using the simulation of the real model, a very wide spectrum of results was obtained. The desired droplet size was determined, i.e., the critical (minimum) diameter was determined for the given parameters of the dust channel's operation. The possibilities of adapting the created model to a specific case were analyzed. The paper presented selected calculation results that confirmed the correctness of the made assumptions. Their analysis gives appropriate guidelines for the construction of a test stand. A system that allows for the precise dosing of the catalyst was presented, and the effectiveness of the developed method was analyzed. The method involved the intensification of the combustion of solid fuels in a boiler with the use of alkyl alcohol and Al-Ni (as combustion promoters). Due to the proposed method, it will be possible to obtain a higher degree of fuel burnout, which, in turn, will result in a lower energy consumption in the case of producing the same amount of power and, also, in an increase of the flame's stability. The parameters of the emitted exhaust fumes and ash will also improve. Based on the conducted process simulations, it was shown that: • in order to supply the catalyst to the system, droplets with a diameter of at least 30 µm should be generated; • an angle of 45 • for the inlet pipe of the isopropanol feeding atomizer is beneficial; this allows the catalyst flux to be introduced into the fuel flux.
The study clearly shows that additional factors must be considered when designing the process of introducing a catalyst into a fuel system. The atomizing device should be carefully matched with the liquid properties and other process parameters. The performed simulations are currently being verified on experimental stands. The research should be continued by extending it with tests on dust channels of various geometries, and with tests conducted for a carrier liquid containing a catalyst with different rheological properties.