Numerical Simulation of the Aerosol Particle Motion in Granular Filters with Solid and Porous Granules

: In this work, a study was carried out to compare the ﬁltering and hydrodynamic properties of granular ﬁlters with solid spherical granules and spherical granules with modiﬁcations in the form of micropores. We used the discrete element method (DEM) to construct the geometry of the ﬁlters. Models of granular ﬁlters with spherical granules with diameters of 3, 4, and 5 mm, and with porosity values of 0.439, 0.466, and 0.477, respectively, were created. The results of the numerical simulation are in good agreement with the experimental data of other authors. We created models of granular ﬁlters containing micropores with different porosity values (0.158–0.366) in order to study the micropores’ effect on the aerosol motion. The study showed that micropores contribute to a decrease in hydrodynamic resistance and an increase in particle deposition efﬁciency. There is also a maximum limiting value of the granule microporosity for a given aerosol particle diameter when a further increase in microporosity leads to a decrease in the deposition efﬁciency.


Introduction
There are many types of filters, each of which performs its tasks and works under certain conditions [1,2]. Granular filters are the most widely used due to their versatility [3]. Depending on the granule size, shape, and the material from which it is made, it can be used to filter gases and liquids of various degrees of pollution and temperatures. Chemicallyresistant materials in the chemical industry can also be used as catalysts [4][5][6][7]. Granular filters are cheap and simple to manufacture, and have a low flow resistance and large surface area compared to other types of filters.
Granular filters can be divided into several types; these are filters with moving, fixed, and fluidized beds. Fixed bed filters have a low particle deposition efficiency relative to other filters, and cannot operate without stopping. Fluidized bed filters have a low efficiency, but can run continuously. Moving bed filters have the highest particle deposition efficiency, and can be operated continuously [8]. Despite the low efficiency and the lack of the ability to work non-stop, the most common of these three types is the fixed bed filter; this is due to its greater reliability and availability, and the simplicity of modeling and predicting the filter parameters in it.
The calculation and prediction of filtration involve the various mechanisms due to which the deposition of particles occurs; these are interception, inertial deposition, diffusion, and gravitational and electrostatic attraction. Each of the mechanisms has a region of particle sizes within which it makes the main contribution to filtration. For example, diffusion and the gravity mechanism have noticeable effects if the particle diameter is less than 1 µm. In the range of particle sizes 1-50 µm, it is assumed that the main contribution is made by inertial deposition [9,10]. The mechanism of particle capture is closely related to the parameters of the filters themselves-more precisely, to the granules, which are the filtering material-such as the method and porosity of the packing, the size, and shape of the granules; the hydrodynamic resistance of the filter depends on these parameters. The liquid or gas flow rate also makes a significant contribution to the filtration process [11][12][13][14][15][16][17][18][19][20][21].
Combining different granules in terms of particle deposition efficiency and pressure drop allows us to achieve better efficiency than the use of only one type of granules [22][23][24]. The authors of experimental work [25] substantiated the expediency of the use of a multilayer granular filter, which reduces the pressure drop by almost 50% while increasing the particle deposition efficiency by 3.23%. The experimental curves of the change in efficiency and pressure drop are given as a function of the filtration time for granule sizes of 3, 5, 8, 10 mm, and a filtration rate for a range of layers depths of 30, 60, 90, and 120 mm. All of these curves are also shown for the proposed multilayer filter model, consisting of a set of 3, 5, and 8 mm granules filled with equal layers.
Complex studies were carried out by Sahar Bakhshian and Muhammad Sahimi [26] to assess the effect of deformation on the morphology of porous media and the flow properties of packed beds. A wide range of parameters were analyzed before and after the deformation. Depending on the morphology of the porous structure, the opening and closing of the pores during deformation can lead to a significant change in permeability while changing parameters such as pore lengths, pore size distribution, and the porosity of the medium. The deformation and swelling of a porous medium when interacts with a liquid were studied in detail in work [27]. These studies are very important when considering the filtration of liquids in porous materials.
The authors of [28] numerically and experimentally investigated the effect of the depth of the filtering layer on pressure drop and particle deposition efficiency. They showed that an increase in the gas velocity increases the deposition efficiency for particles larger than 7 µm, which can be explained by their increasing inertia. It was also shown that the pressure drop changes linearly with the increasing layer depth. The deposition positions of a granular filter with random packing were determined in the work of the authors [29] for particles with sizes from 1 to 21 µm; the maximum area during deposition on granules is occupied by particles from 1 to 7 µm, and with an increase in particle size, the area of deposition decreases. The authors of the paper [30] investigated the filtration process in models of ordered body-and face-centered cubic packing granules. Parameters such as the particle deposition fraction-depending on the number of rows for two models of granule packing and pressure drop versus flow rate-were estimated. The behavior of the efficiency curves, represented as a function of the particle diameter, changes when the flow rate varies from 0.2 m/s to 1 m/s. In the indicated velocity range, the efficiency curves have a minimum at a particle size of 2 µm for the model's face-centered cubic packing, and at a particle size of 4 µm for the body-centered cubic packing model (except for the velocity value of 0.2 m/s). Later, the same authors [31] investigated a random packing of granules, which did not reveal an obvious minimum in the efficiency of the particle deposition in a similar range of calculated velocities for particles of 2 or 4 µm, which can be explained by the closeness of the last model to the real filter structure, and the absence of the influence of boundary conditions. In [32], heat transfer characteristics studies were added to the study of the deposition efficiency versus the particle diameter, and the pressure drop versus the bed depth.
An experimental study of various flow regimes within pores in a dense granular layer by the electrochemical method was carried out [33]. The Reynolds numbers in the study varied from 20 to 2200; the transition from a laminar to a turbulent flow regime was considered. The results demonstrate the effect of the tube diameter on the flow regime.
Studies of granular filters can be divided into several types. These are empirical, semi-empirical, and Computational Fluid Dynamics (CFD) modeling. Most filters in use today were invented in the early 2000s, and have not undergone any major changes. Their development was carried out based on semi-empirical dependencies, which were developed for a simple prediction. This method of application is universally impossible, as each case requires its correlations, considering the peculiarities of the shape and packing of granules [34]. Experimental studies give accurate results, but this type of research is expensive, and does not allow us to see how the gas moves in the space between the granules, or how the deposition of particles of different sizes occurs. The use of numerical simulation technology allows us to consider in detail each parameter, predict the filter behavior under various conditions, and consider the movement of the particles inside the filter. The main problem of numerical simulation is the creation of a reliable model that will correspond to a real filter; compliance with this condition allows us to obtain data that correlate well with the experimental results. In most studies, in order to create a geometric model of a filter with randomly-packed granules, the Monte Carlo method is combined with the DEM (discrete element method) [35][36][37][38]. Computed tomography technologies are also used to create a geometric model from the images obtained. This method is well suited for making ready-made filter models, but it does not change its parameters.
Introduced by Cundall and Strack [39], the discrete element method (DEM) was widely used to simulate bulk media behavior and, in recent years, was actively used by researchers to simulate particle motion in granular filters. This method offers the possibility of deeply understanding the motion of particles of different shapes in granular filter models. Shire and O'Sullivan, in their work [40], analyzed the constriction size distribution on models of granular filters that were created using DEM. The work [41] is devoted to the study of the motion of a wet granular layer. The breakage phenomenon of particles with sharp edges is considered in the article [42] through the combination of the discrete element method (DEM) and the finite element method (FEM). The work [43] is devoted to the study of the particle's shape's influence on the behavior of the granular layer. Experiments were carried out with rice grains and glass beads on an unbaffled rotary drum using a high-speed camera, the results of which were repeated through numerical simulation using DEM. The authors of [44] estimated the change in permeability during the deposition of large-and small-sized particles in a porous medium consisting of a disordered packing of spheres by combining DEM and the network model.
In their investigations, the authors of work [45] used the DEM method to determine the particle collision frequency on the surface of particles on any plane perpendicular to the flow direction. Although He et al. [46] argued that this method is computationallydemanding in terms of calculations, Guo et al. [47] observed that some particle-scale quantities-such as the position, orientation, and contact force-could be obtained by implementing the DEM simulation technique. Moreover, the authors of work [48] applied a 3D geometric modeling assumption to illustrate a decrease in the computational efforts based on specific particle dispositions in granular filters.
In recent years, the formation of particle piles was studied based on diverse granular shapes. The wedge sequence and the layered sequence were commonly used to assemble a pile. While, in the first method, a pile is formed from a particle point source, in the second, a pile is formed by dropping a layer of particles from a specific height of less than one particle diameter [49][50][51]. Many authors have investigated the relationship between stress distribution and particle sizes in granular models. Matuttis et al. [50] analyzed the particle deposition of a circular and a polygonal shape. On the other hand, Zhu et al. [52] observed that particle elongation is possible with a pressure drop at the bottom of a pile. Furthermore, the authors of the paper [53] focused on the study of the pile formation of ellipsoidal granules.
Coupled with the computational fluid dynamics solver (CFD), the discrete element method (DEM) has become particularly attractive for the simulation of large-scale industrial processes due to excellent mass and heat transfer attributes that result from its application [54].
This paper presents a CFD simulation of the motion of aerosol particles in granular filter models created by DEM for spherical granules with a diameter of 3, 4, and 5 mm. Parameters such as the pressure drop, particle deposition efficiency, and filter quality factor were investigated. It was suggested that it is possible to use filters with porous granules to reduce the pressure drop and increase the particle deposition efficiency. Studies of the influence of the granule micro-porosity on the above-mentioned parameters in the granular layer were carried out. In general, porous granules can be widely used in practice.

Filter Model
When choosing filters, the question arises as to which filter is the best for a given technological cycle. Current granular filters can have complex internal structures, where gas dynamics is the main contributor to particle deposition efficiency. The different geometries of the internal structure of filters affect the flow field and, consequently, the efficiency of particle deposition. This work aimed to determine the granular filter parameters (granule size, porosity, microporosity) that make the most significant contribution to the change in the deposition efficiency of aerosol particles, the resistance, and the filter quality factor.
For submicron-sized particles, the predominant deposition mechanism in porous structures is inertial deposition. For granular filters, the actual granule size is in the range of 1-10 mm. For example, in works [28,29,32], experimental and numerical studies of the operation of granular filters with granule sizes of 5, 6, and 10 mm were carried out. In works [30,31,55], the operation of a filter with granules with a diameter of 5 mm was investigated. In works [56,57], filters with a granule diameter of 1-5 mm were selected. In [25], experimental studies were carried out for filters with granule diameters of 3, 5, 8, and 10 mm. We made a numerical simulation for a set of models of granular filters with granule diameters of d g = 3, d g = 4, and d g = 5 mm. The discrete element method (DEM) was used to create the geometries in YADE [58]. An example of the creation of a granular filter bed is shown in Figure 1. The figure shows four pictures of the process of pouring granules into the cylinder at different time intervals: 0.15 s, 0.21 s, 0.65 s, and 1.0 s. In the calculations, the density of the granules is 1800 kg/m 3 . Due to the random arrangement of the granules in the filter package, the calculations were averaged for five different geometries, which were created with similar physical parameters to eliminate errors. For the calculations, a cylindrical area with a diameter of 20 mm and a length of 20 mm, filled with spherical granules, was selected ( Figure 2). The number of granules was from 221 to 235 units for the case of d g = 3 mm, from 98 to 103 units for the case of d g = 4 mm, and from 45 to 47 units for the case of d g = 5 mm. For a small size mm, the number of granules is significantly larger. In this case, the deviation in the number of granules relative to the average value is from 4.4% (d g = 5 mm) to 6.3% (d g = 3 mm). An inlet section, 10 mm long, was located before packing the granules to form a gas flow, and after the granules, there was an outlet part 40 mm long. The parameters of the models created for the study are presented in Table 1. We verified the constructed models by comparing the porosity parameters of the bed of spherical granules created by the DEM method with the following formula [59]: where d t is the tube diameter. Calculations using the Formula (1) are also presented in Table 1. In this work, we calculate the porosity of the granule layer, limiting the considered volume from below and from above by the boundary point of the presence of the granule. Therefore, the calculated values of porosity ε are less than those calculated by the formula ε (1) . From Formula (1), it follows that the porosity depends on the diameter of the area into which the granules are poured. Regardless of the granule size, with a large increase in the size of the region (d t /d g > 15), the porosity tends to a value of 0.39 for the dense packing of the spheres. A strong influence in the value of the porosity of the packing of the spheres is observed at d t /d g < 10. This paper does not consider the effect of the tube size on the efficiency of the filtering layer of granules. We will investigate the effect of the granule diameter, as well as the presence of micropores, on the performance of a granular filter for a fixed area.

Mathematical Formulation
We calculated the gas flow (airflow) based on the solution of the equations for the conservation of mass and the conservation of momentum for a stationary problem.
Mass conservation equation: where → V is the gas velocity. The momentum conservation equation: where ρ is the gas density, P is the pressure, and µ is the gas dynamic viscosity coefficient. The trajectories of aerosol particle motion were calculated in the found field of the gas velocities based on Newton's law for single-particle motion: where → R p is the vector of the coordinates of the aerosol particle, → V p is the velocity of the aerosol particle, and τ is the relaxation time of the particle, calculated as in [60]: where ρ p is the particle density used in calculations 1500 kg/m 3 , and d p is the aerosol particle diameter. It is assumed that the aerosol particles that reach the filter surface are deposited on it. In this case, the captured particles are removed from the computational domain [61,62].
At the inlet boundary of the computational domain is set the value of the mass gas flow rate; at the outlet boundary is set the value of atmospheric pressure; other boundaries were taken by default as impermeable walls. The solution to the problem of stationary gas flow (2) and (3) was obtained by the finite volume method in the CFD package ANSYS Fluent (v. 19.2). The geometry of the porous region is such that, even at low flow rates, the active mixing of the flow occurs, and a turbulent flow regime is observed. The average flow rate accepted in the calculation of the particle trajectories was 1 m/s, which corresponds to a turbulent flow regime. We decided to use direct numerical simulation (DNS) in the calculations. The number of elements of the grid division of the region in the calculations is 12,000,000.
After calculating the gas flow field, we calculated the particle motion using Formulas (4) and (5) in ANSYS Fluent (v. 19.2). Spherical particles start at the computational domain's entrance and their direction is normal to the boundary. Then, the particles move against the background of the calculated gas field; for the presented calculations, we started with 30,000 particles. Only the trajectory of motion of the particle center was considered; the density and diameter of the particle affect the inertia and contribute to the deposition, with a sharp deviation of the flow near a solid surface. The particles that reached the surface of the filter granules were considered settled, and were excluded from the further calculation of the trajectories' motion. The particles that reached the surface of the tube were reflected, and their further movement was calculated.

Numerical Simulation of the Aerosol Motion in a Granular Filter with Solid Spherical Granules
The numerical calculations were compared with the data of the semi-empirical dependence for the pressure drop ∆P relative to length L of a granular medium [63], in order to assess the correctness of the gas flow model: In [64], a formula was proposed to take into account the tube diameter: Figure 3 shows a comparison of the numerical simulation results and the application of semi-empirical Formulas (6) and (7) for the pressure drop. The graphs show that the pressure drop curves correlate well with each other for all three granule diameters: 3, 4, and 5 mm. Formula (7) was designed to consider the effect of the tube walls near the granules on the pressure drop. In [64], the existing models for the calculation of the pressure drop are considered in detail, and the difference in the value of the pressure drop for an unrestricted layer of granules and a layer of granules enclosed in a narrow tube is described. In our case, despite the presence of a tube, the graphs of the calculated pressure drop are closer to the curve calculated by Formula (6) for an unlimited layer of granules than to the curve calculated by Formula (7) for a limited layer of granules. We believe that this is because a layer of granules of small thickness was considered, and in the filter zone, there was no time to form a near-wall layer sufficient to influence the pressure drop. In [32], experimental and numerical calculations of the pressure drop for a granular layer in a confined tube are performed, which show good agreement with the calculations by Formula (6). Furthermore, we studied the particle deposition efficiency, which is defined as the ratio of the number of particles deposited on the filter surface to the total number of particles that started at the inlet boundary, which is given as where n trap p is the number of deposited particles, and n total p is the total number of starting particles. Figure 4 shows an example of the particle trajectories in a filter for a particle diameter of d p = 1.8 · 10 −6 m. We performed a comparative calculation to verify our research. In studies [28,29], the filtration efficiency curve for a granular filter with a granule diameter of d g = 10 mm is presented. We constructed a cylindrical region with a granular filter, as in [29], calculated the gas movement at a similar velocity of 0.345 m/s, and calculated the motion and deposition of aerosol particles with a diameter of 1 µm to 21 µm. The results are presented in Figure 5. The result of our calculation is in good agreement with the results of the experimental studies. Figure 5. Comparison of the filtration efficiency of a granular filter with the experimental data [29]. Figure 6a shows the results of a study of the particle deposition efficiency for the filter models with spherical granules for various values of the porosity of the medium. Calculations of the efficiency were carried out for all of the constructed geometries, and thus, for each value of the granule diameter, we averaged the results over five calculations. We can see that the filter with d g = 3 mm granules has the best efficiency among the cases of our calculations. In this case, the largest scatter in the results (up to 8.6%) can also be observed for granules with a diameter of 3 mm. For packages of granules with a diameter of 5 mm, the scatter of the results is not large (up to 2.35%). Let us consider the filter quality factor. In experimental studies [65], where the number of particles that are not settled on the filter is estimated, one can find the expression Q f = − ln(C out /C in )/∆P, where C out and C in are the concentrations of particles at the outlet and inlet of the filter, respectively. In our research, we calculated the filter quality factor as the ratio of the particle deposition efficiency to the pressure drop characterizing the hydrodynamic resistance of the medium. Figure 6b shows the calculation results. For small particles (up to d p = 3 · 10 −6 m), a filter with spherical granules with a diameter of d g = 3 mm has the highest quality factor among the cases of our calculations. The largest deviation on the graphs is for deposition efficiency. As the particle diameter increases, its quality factor curve is lower than that of all of the other granule models. The model with granules with a diameter of d g = 5 mm has the highest quality factor for larger particle sizes. This effect is explained by a decrease in the deposition efficiency value deviation with an increase in the diameter of the aerosol particles. Obviously, with a large value of the particle diameter, the deposition efficiency will tend to 1 (provided no transparency and holes in the filter), at which point the pressure drop will completely determine the filter quality.

Numerical Simulation of the Aerosol Motion in a Granular Filter with Porous Spherical Granules
Furthermore, we studied the possibility of improving the filtering properties of granular filters with spherical granules. For this, we modified the filter with porous granules with a diameter of d g = 5 mm (Figure 7). The micropore diameter was d mic = 1.9 mm. By increasing and decreasing the distance between the micropores, we controlled the porosity value. The porosity parameters of the created models are presented in Table 2.  The addition of micropores entails an increase in the mesh elements required for the calculations. As such, if for the case of solid spheres there were 12,000,000 elements in the mesh, then for the case of total porosity ε = 0.635 the number of elements is about 14,000,000, and for the case of total porosity ε = 0.825 the number of elements is about 16,000,000. Figure 8 shows the results of a study of the particle deposition efficiency for filter models with modified granules compared to the model without micropores. It is obvious that the addition of micropores affects the behavior of the movement of the liquid and aerosol particles inside the granular filter. Particles, moving with the gas flow, can penetrate into the granules and settle during their movement in the internally-developed porous structure. Figure 8 shows that with the modification of a granular filter with micropores, an increase in particle deposition efficiency is observed in the entire range of particle sizes, which ranges from 10% to 25%. In this case, with a slight increase in porosity (ε = 0.55, ε = 0.635), the deposition efficiency curves qualitatively repeat the character of the curve, as for a bed of granules without micropores. Then, with an increase in porosity (ε = 0.75), the character of the curve begins to change. Models of filters with different total porosity values show different results for different ranges of values of particle diameters. Thus, a filter with total porosity ε = 0.635 demonstrates an efficiency of particle deposition higher than others by 5-8% in the area of particles up to 2·10 −6 m, and a further increase in porosity leads to a decrease in efficiency. Then, we can observe the intersection of the curves of the aerosol particle deposition efficiency for different porosity values. In the area of large particles, filters with total porosity of ε = 0.811, ε = 0.825 and ε = 0.843 work better. A structure with porosity ε = 0.825 has the most excellent efficiency, and a further increase in porosity also leads to a decrease in efficiency. Thus, small and large particles have their own porosity, which determines the maximum deposition efficiency. Obviously, the addition of multiple pores should increase the total surface area of the bed. However, too many pores can reduce the surface area of the bed. Figure 9a shows the calculated dependence of the bed area on porosity. The largest surface area is observed for the porosity range ε = 0.635 − 0.811. Here the analysis of the averaged efficiency curves is carried out; the analysis of the scatter of the results for specific values of the diameter of the aerosol particles is discussed below. This conclusion is valid for the special cases of geometric characteristics that we have chosen for the calculations in this study.  Figure 9b shows a graph of pressure drop versus flow rate. The structure with the highest porosity has the lowest pressure drop. The addition of micropores helps to reduce the flow resistance of the filter. Figure 10 shows the dependence of the filter quality factor on the particle diameter. The original filter model without micropores has the minimum value of the quality factor. This result is due to the combination of the study results of the particle deposition efficiency and pressure drop; in both indicators, the model without micropores is inferior to the models with micropores. Since the pressure drop makes the main contribution to the quality factor, the structure with the highest porosity value ε = 0.843, which corresponds to the lowest pressure drop, has the highest performance. Figure 10. Changes in the quality factor versus the particle diameter for a filter with microporous granules.
The calculation results show that it is impossible to single out a filter with one porosity value as the most effective in a wide range of aerosol particle sizes. Next, we considered the values of the deposition efficiency of aerosol particles and the filter quality factor for particle size values d p = 1 µm, d p = 4 µm, d p = 10 µm, and d p = 30 µm. Figure 11 shows the results.
For small particle diameters, the most excellent efficiency is achieved for porosity ε = 0.635, when the filter has the largest surface area. For the deposition of small, weaklyinertial particles, the curvature of the rotation of the streamline of the carrier gas flow near the filter surface plays an important role. The porosity ε = 0.635 provides a large surface area due to the hollow spherical inclusions inside the granule, with several connecting channels. Small particles trapped inside the porous granule along with the gas are more likely to settle inside the micropore. A further increase in porosity leads to the destruction of the granule's surface, so the total surface area remains almost unchanged. Many micropores and connecting holes inside the granule do not contribute to the capture of small particles, which deviate slightly from the gas flow lines. At the same time, we found that the maximum value of the filter surface area does not provide the greatest value of the large particle deposition efficiency. The highest efficiency is determined at ε = 0.825. A further increase in porosity leads to the formation of straight, open channels in the filter, and particle motion without contact with the granules. Figure 11a also shows the scatter bars for the calculation of the deposition efficiency of aerosol particles. For each value of the total porosity, we created five sets of granule packs, as for solid spheres. For each constructed filter geometry, the calculation of the movement and deposition of the aerosol particles was carried out. We observed the greatest scatter of results (up to 7%) in our calculations for the porosity value ε = 0.635. The scatter of results decreases both with a decrease in porosity until a solid sphere is reached, and with a strong increase in porosity. This may be because, for ε = 0.635, the number of pores in a separate granule is small (see Figure 7), and their position of the granule in space can have a strong effect on the movement of gas and aerosol particles in the filter. With a high porosity of ε = 0.843, each granule has a highly-developed internal structure, such that the position of the granules in the filter layer does not have a critical effect on the filtering properties, and the filter itself becomes closer to an open cell foam. At the same time, the scatter of the results is less than 1%. Figure 11. Dependence of the filter characteristics for particle size values d p = 1 µm, d p = 4 µm, d p = 10 µm and d p = 30 µm: (a) particle deposition efficiency on porosity; (b) particle deposition efficiency on surface area; (c) filter quality on porosity; (d) filter quality on surface area.
Considering the scatter of the results for the selected particle diameter d p = 1 µm and d p = 30 µm, we can say that, for the filter parameters used, the deposition efficiency of the aerosol particles is almost independent of the porosity of the granules. The main attention here can be paid to the reduction of the aerodynamic drag of the filter. For the case of d p = 4 µm, the porosity value ε = 0.635 provides an advantage over filters with other porosity values in particle deposition efficiency, even taking into account the larger scatter of the calculation results. This conclusion was obtained for the specific case considered, but these results could be useful in studies to optimize the operation of granular filters.

Conclusions
We carried out a numerical simulation of aerosol motion in models of granular filters with granules of various sizes, i.e., 3, 4, and 5 mm in diameter. We investigates the particle deposition efficiency over a wide range of sizes and aerodynamic drag as a pressure drop.
The d g = 3 mm granule filter has the highest particle deposition efficiency, but also the highest resistance. The analysis of the filter quality factor parameter made it possible to divide the scale of the aerosol particle diameters into areas in which, for small particles, it is preferable to use a filter with a denser packing of granules (d g = 3 mm), and for large particles, the deposition efficiency values are close, and its aerodynamic resistance determines the filter quality.
The analysis of the calculation results for a filter modified by micropores with spherical granules with a diameter of d g = 5 mm shows that, for almost all of the models with any microporosity of the granules, there is an increase in the particle deposition relative to the original model. The particle deposition efficiency does not change much with the addition of microporosity to the granules for very small and very large particles. For medium-sized particles, the optimum porosity value can be determined, i.e., the value at which the highest deposition efficiency is achieved. In our calculations for the case of d p = 4 µm, we obtained the optimal value of porosity ε = 0.635, for the case of d p = 30 µm-ε = 0.825. At the same time, by adding micropores, we can significantly reduce the aerodynamic drag of the filter. For all of the considered cases, the value of the filter quality factor reaches the highest value at the maximum value of the total porosity of the granule layer.
In this work, the presented results and conclusions refer to a specific filter geometry. These results may be of interest for the design of new granular filters, and the optimization of the operation of existing granular filters. We foresee further research in this direction in the tasks of the analysis of the effect of micropore size, granule size, and filter size on the performance of granular filters.