Next Article in Journal
Multi-Objective Multi-Stage Optimize Scheduling Algorithm for Nonlinear Virtual Work-Flow Based on Pareto
Next Article in Special Issue
Numerical Simulation of Dynamic Variation Characteristics of Particles in a Rolling Fluidized Bed
Previous Article in Journal
Study and Validation of a Novel Grouting Clamp Type Deepwater Oilfield Pipeline Repair Method Based on Computational Fluid Dynamics
Previous Article in Special Issue
Uniaxial Creep Test Analysis on Creep Characteristics of Fully Weathered Sandy Shale
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mesoscale Simulation Approach to Study the Flow Field in an Axial Granular Bed Filter

1
School of Chemical Engineering, Hebei University of Technology, Tianjin 300132, China
2
School of Mechanical Engineering, Hebei University of Technology, Tianjin 300132, China
*
Authors to whom correspondence should be addressed.
Processes 2023, 11(4), 1146; https://doi.org/10.3390/pr11041146
Submission received: 17 February 2023 / Revised: 25 March 2023 / Accepted: 31 March 2023 / Published: 7 April 2023
(This article belongs to the Special Issue Multiphase Flows and Particle Technology)

Abstract

:
In an axial granular bed filter (GBF), a new mesoscale simulation approach is obtained by combining the macroscopic calculation models, i.e., the equations of the total pressure drop and dust-removal efficiency into the porous media model and the source term of the conservation equations. After grid-dependent tests and experimental validation, the effects of different conditions, i.e., granular bed height L, superficial gas velocity ug, dust diameter dp, dust concentration cp, granular diameter dg, initial bed voidage ε0, and filtration time t, on the pressure drop and dust-removal efficiency are investigated. The results show that the pressure drop is related to the inertial and viscous resistance terms, which increase with increasing L, ug, cp and t and decreasing ε and dg. The dust-removal efficiency is related to the Reynolds number, effective Stokes number, and equivalent granular diameter ratio. It increases with increasing L, ug, dp and t (small values), and decreasing cp, ε, and dg. Moreover, the influence of different conditions coincides well with dust-removal efficiency in relevant studies, which further demonstrates the accuracy of the mesoscale simulation approach. With the application of this method, the flow field can also be obtained easily and quickly, which is expected to provide a reference for the simulation study of GBF.

Graphical Abstract

1. Introduction

Large amounts of fine-dust emissions affect air quality and visibility, which also endangers human health. With the development of the economy and the improvement of people’s living standards, people in various countries are becoming more and more concerned about air pollution and have set strict emission standards: the World Health Organization (WHO) air quality guidelines of 2021 [1] stipulate that the annual average emission concentration of PM2.5 should not be higher than 5 μg/m3. A granular bed filter (GBF) is a device that uses granular media as a filter layer to purify dust-containing gases. Granules are accumulated inside the GBF, and the dust inside is captured by the granules through inertial collision, diffusion, interception, gravitational settling, and electrostatic forces [2,3]. Compared with inertial, cyclone, wet, and electrostatic de-dusting devices, GBF has the advantages of a high dust-removal efficiency (especially <10 μm), high-temperature and high-pressure resistance, simultaneous desulfurization, denitrification, the removal of volatile organic compounds (VOCs) and drying of materials, etc. It is widely used in thermal power, steel, petrochemicals, and other industries [4,5].
GBF contains filter media, gas, and dust. The dust can flow with the gas or accumulate on the surface of the filter media by inertial collision and diffusion, making the flow field more complex. At present, a lot of research has been conducted to reveal the flow-field distribution law and to propose macroscopic calculation models for pressure drop and dust-removal efficiency. The pressure-drop model is generally modified based on the Ergun equation by introducing parameters such as the amount of dust deposited on the surface of the filter media. For example, Xiao et al. [6], Robert et al. [7], Denis [8], and Yin et al. [9] have systematically studied and analyzed the pressure drop in the GBF and found that the pressure drop increases with increasing dust deposition. Moreover, the macroscopic calculation model is established for the pressure drop during dust loading by modifying the traditional Ergun equation and considering the effect of dust deposition. Based on different assumptions, different macroscopic models of dust-removal efficiency have been developed. For example, Iwasaki et al. [10] and Ruben et al. [11] assumed that the concentration of dust deposited on the surface of media granules has an exponential distribution, and then a macroscopic model of dust-removal efficiency is proposed. Tardos et al. [12] and Jung et al. [13] developed a single-granule dust-removal efficiency model based on the interaction forces between dust, gas, and filter granules, and modified it to obtain a granular dust-removal-efficiency model. Yu et al. [14] and Bruno et al. [15] constructed a dust-removal-efficiency model by assuming that the bed consists of a series of unit bed elements (UBE) connected in series. The dust-removal efficiency of each UBE is the same under the initial conditions. The dust-removal-efficiency model constructed by UBE theory includes the Reynolds number, Stokes number, equivalent granular diameter ratio, bed height, and other numbers, and takes into account the influence of different conditions and dust deposition, which is chosen in this study.
The macroscopic calculation models of pressure drop and dust-removal efficiency can reflect the bed performance of GBF. However, most of the fluid flows in real industrial installations are non-plug flows and the calculation results may have large errors. Numerical simulation is a much easier way to obtain detailed information about the flow field inside the bed than experimental studies; this method is also suitable for GBF with different structural parameters. If the macroscopic models can be combined into numerical simulations, the flow-field distribution in the granular bed filter can be obtained more quickly and accurately.
Different simulation methods are used in the modeling of the flow field in GBF, e.g., the Eulerian–Lagrangian two-phase flow (CFD-DEM) method and the Eulerian single-phase flow method. For the CFD-DEM method, filter granules are calculated by the Lagrangian method; gases and dusts can be calculated by the combined models of Eulerian single-phase flow, multi-component transport, and DPM. For example, Chen et al. [16] and Li et al. [17] used CFD-DEM to obtain the dust-removal characteristics of a fixed GBF before the formation of filter cake. Wang et al. [18], Chou et al. [19], Kloss et al. [20], and You et al. [21] used the combined models of an Eulerian single-phase flow and DPM to investigate the flow-field characteristics in a randomly filled GBF and obtain the variation in initial pressure drop and dust-removal efficiency. In addition, the flow-field space can be directly simulated after the volume deletion of filtered granules. For example, Chen et al. [22], Guan et al. [23], and Wang et al. [24] used the DEM method to determine the location of filter granules under random accumulation, and then numerically simulated the fluid space after the removal of filter granules in GBF. The above methods allow us to investigate the interaction forces between dusts, gas, and filter granules. The dust-filtration mechanisms are also analyzed according to the simulation results such as dust collision, gravitational settling, diffusion, and interception. However, large installations contain tens of millions or even hundreds of millions of granules. The interaction forces between the gas, the dust, and the filter granules also increase the computational difficulty and workload.
In contrast, the use of a porous media model can significantly reduce the computational effort of the simulation. For example, Ahmed et al. [25], Wojciech et al. [26], Cheng et al. [27], and Rodrigo et al. [28] used the porous media model to simulate the flow-field information in the granules and obtained accurate pressure drop data. However, the variation in dust-removal efficiency and flow field with dust deposition in the granules cannot be obtained by a porous media model.
Therefore, this study combines the macroscopic models, i.e., the pressure drop and the dust-removal efficiency equations based on unit bed elements (UBE) theory [15], into the porous media model and the source term of the conservation equations. A new mesoscale simulation method is then obtained. A three-dimensional full-scale model is developed based on the structural parameters of the axial GBF described in reference [29]. Grid-dependent tests and model validation are performed. The effects of different operating conditions and physical parameters, e.g., bed height, superficial gas velocity, dust diameter, dust concentration, granular diameter, bed voidage, and filtration time, on the flow field are investigated. The feasibility of this mesoscale simulation method can be further verified when the variation laws of pressure drop and dust-removal efficiency obtained by the simulation results coincide well with the experimental results.

2. Numerical Simulation

2.1. Geometry Model and Mesh

Based on the structural parameters of the axial granular bed filter described in reference [29], a cylindrical flow field with a diameter of 58 mm and a height of 170 mm is created. A section of the granular bed is in the middle, as shown in Figure 1. One side of the cylindrical bed is the gas inflow, where the randomly generated dust granules are mixed with air to form a dusty gas. When the dusty gas flows through the granular bed, the dust is removed by the filter granules and the cleaned gas is discharged from the outflow side. The constants in the equations are dependent on the gas, dust, and filter parameters.
The gas is air at room temperature and pressure with a density of 1.205 kg/m3 and a kinetic viscosity of 1.81 × 10−5 Pa·s. The filtered granules and dust are glass spheres and silica with a bulk density of 2500 and 1400 kg/m3, respectively.
Different parameter values are varied in the simulation to examine their effects on the flow field. The parameters examined include bed height L = 20–100 mm, superficial gas velocity ug = 0.15–0.55 m/s, dust concentration cp = 100–100,000 mg/m3, dust diameter dp = 1–10 μm, granular diameter dg = 1–10 mm, bed voidage ε = 0.40–0.55, and filtration time t = 0–120 s. In order to facilitate the analysis, a baseline working condition is chosen: L = 20 mm, ug = 0.345 m/s, cp = 6985 mg/m3, dp = 2.24 μm, dg = 10 mm, ε = 0.49, t = 0 s.
Gambit 2.4.6 software was used to draw the geometry model and mesh. Three fluid calculation areas were set up: left fluid, middle granules, and right fluid. The structured hexahedral mesh was created by combining O-planar partitioning and Cooper’s methods. The maximum distortion of the mesh was less than 0.5, as shown in Figure 1.

2.2. Mathematical Model

In this study, a distinction is made between the dust flowing with the gas and the dust deposited on the surface of the filter granules, which are assumed to belong to different phases. The Eulerian two-phase flow and porous media models in Fluent 6.3.26 software are used to the simulation studies: the gas and the carried dust are considered as one phase and simulated by the multi-component transport model; the dust deposited on the surface of the filter media is the other phase, and the influence of the filter granules on the flow field is simulated by the porous media model. In addition, the macroscopic model of pressure drop and dust-removal efficiency is combined into the porous media model with the source term of the continuity equations in the form of a user-defined function (UDF). The flow-field information in GBF is then obtained by this mesoscale simulation method. It refers to the fact that the variations in the parameters are calculated from a mesoscopic view (each mesh), rather than a macroscopic view (the whole bed) or a microscopic view (each dust and granules).

2.2.1. Continuity and Momentum Conservation Equations

The Eulerian two-phase flow model is applied in this study. The primary phase is the mixture of air and carried dust granules, which are simulated by a multi-component transport model. The local average continuity Equation (1) and the momentum Equation (2) are used. The mass sources of the gas and dust are zero and -dc, respectively. The interaction forces and turbulence are calculated by Equation (3) and k-ε RNG models, respectively, which are not described here [30].
The secondary phase is the dust phase accumulated on the surface of filter granules, which is calculated by using the local average continuity Equation (3). As the dust phase accumulates on the granule surface, the velocity of the dust phase is set to 0 m/s in this study. It is no longer calculated by the momentum conservation equation.
( ρ i α g ε 0 Y i ) t + ( ρ i ν α g ε 0 Y i ) = α g ε 0 S i , i = g   o r   s ,   S g = 0 ,   S s = d c , α g = ε / ε 0
( α g ε 0 ρ i Y i v i ) t + ( α g ε 0 Y i ρ i v i v i ) = α g ε 0 p i + α g ε 0 Y i ρ i g + α g ε 0 τ ¯ i + α g ε 0 ( - ρ i Y i v i v i ¯ ) + α g S ε 0 S v
t ( α s ε 0 ρ s ) + ( α s ε 0 ρ s v s ) = α s ε 0 S s s ,                           S s s = d c , α s = 1 α g

2.2.2. Porous Media Model

The resistance effect of the granular bed on the flow field is realized by the application of a porous media model (4). In the momentum conservation equation, the macroscopic model of pressure drop given in [15] is used to change the coefficients of the viscous and inertial resistance terms (5 and 6) by means of UDF.
S v = ( μ α v i + C 2 2 ρ g | v i | v i )
1 α = ( 1 + b 1 σ ¯ m ρ p ε 0 ) b 2 ( 1 σ ¯ m ρ p ε 0 ) b 3 150 d g 2 ( 1 ε 0 ) 2 ε 0 3
C 2 = ( 1 + b 1 σ ¯ m ρ p ε 0 ) b 2 ( 1 σ ¯ m ρ p ε 0 ) b 3 1.75 d g ( 1 ε 0 ) ε 0 3
σ ¯ m = ρ p ( ε 0 ε )
where ρp is the dust density, dg is the filter granular diameter, and ε0 is the initial bed voidage. σ ¯ m is the dust deposition quantity and can be calculated by Equation (7). b1, b2, and b3 are constants, which are fitted by the experimental results of pressure drop (taken as 28.92, 1.254, and 4.2171, 159 [15], respectively, in this study) as shown in Equations (5) and (6).

2.2.3. Source Term in Continuity Equations

By changing the source term of the continuity equations, the dust is transferred from one phase to another in the granular bed region. The amount of transfer affects the dust-removal efficiency. According to the UBE theory [15], the conservation and macroscopic equations are connected by Equation (8). In each cell, the gas velocity remains unchanged with the positions. Thus, when “dt” is sufficiently small, the gas flows a distance of “vgdt” where vg is the real gas velocity in each cell. According to the UBE theory, the mass sources of dust “dc” can be calculated by Equation (8) when the distance of gas flow is “vgdt”. However, the gas velocity may vary widely in different cells, which allows for the attainment of the flow-field distribution in an axial granular bed filter by this approach. The thickness l of each layer can be calculated by Equation (9). In addition, the actual bed voidage ε within the granular bed will, in turn, affect l.
d c = ρ g Y g [ 1 ( 1 E ) v g d t l ] 1 d t
l = [ π 6 ( 1 ε ) ] 1 / 3 d g
For the UBE theory, the granular bed is considered as a series of unit beds. According to the work of [15], the total dust-removal efficiency E is related to the initial dust-removal efficiency e0 and can be calculated using Equation (10); e0 is related to the dust-removal efficiency of the single unit bed η0, as shown in Equation (11); F ¯ is the average relative filtration coefficient, which is an important parameter to characterize the deep-filtration performance. The constants of a1, a2, and a3 are fitted by means of the experimental results of the filtration efficiency (taken as 290,600, 0.1374, and 4.338 [15], respectively, in this study), as shown in Equation (12).
E = 1 ( 1 e 0 ) F ¯
e 0 = 1.209 η 0
F ¯ = ( 1 + a 1 σ ¯ m ρ p ε 0 ) a 2 ( 1 σ ¯ m ρ p ε 0 ) a 3
According to Equations (13)–(19) [15], η0 is related to the effective Stokes number Nsteff, dust diameter dp, granular diameter dg, and adhesion probability γ. The constants are obtained directly from the work in [15], which summarizes the relevant research. Nsteff is a modified value of the Stokes number Nst [13], which characterizes the ratio of inertial to diffusion forces. The smaller the value, the smaller the inertia of the dust, the easier it is to follow the fluid motion, and the more pronounced the diffusion effect. The dust-removal efficiency tends to decrease. In order to better predict the traction effect on dust of the gas, the Cunningham correction factor Cs is introduced by taking into account the discontinuous medium effect of dust transition between continuous and free molecular flow. The Reynolds number Re is the ratio of inertial and viscous forces. In Equation (18), ρg, vg, and μg are the density, the superficial gas velocity, and the viscosity of the dusty gas, respectively; dg is the diameter of the filter granules. As Re increases, the inertial effect increases, and the probability of inertial collision between dust and filter granules increases significantly.
η 0 = a 4 N s t e f f 1.3437 ( d p d g ) 0.23 γ
N s t e f f = [ A s + 1.14 R e 0.5 ε 2 / 3 ] N s t 2
γ = 0.00318 N s t 1.248 , N s t 0.01
N s t = ρ p d p 2 v g C s 9 μ g d p
C s = 1 + ( 6.635 × 10 8 d p / 2 ) ( 1.252 + 0.399 e 1.10 ( d p / 2 6.635 × 10 8 ) )
Re = ρ g v g d g μ g
A s = 2 [ 1 ( 1 ε ) 5 / 3 ] 2 3 ( 1 ε ) 1 / 3 + 3 ( 1 ε ) 5 / 3 2 ( 1 ε ) 2

2.3. Parameter Settings

The Eulerian two-phase flow, multi-component transport, and porous media models in Fluent 6.3.24 software are used for the simulation study. The inlet and outlet boundary conditions are set as the velocity inlet and pressure outlet, respectively; the wall surface is selected as a no-slip boundary condition. The simulation type is chosen as transient simulation, and the algorithm is chosen as SIMPLE. The time step is chosen as 0.001 s with a maximum iteration step of 20.
In this study, a single-factor variable method, or the one-factor-at-a-time method, is used to investigate the effect of different conditions on the flow field by changing a parameter value in a baseline working condition (Section 2.1).

2.4. Grid-Dependent Test and Experimental Validation

According to the relevant numerical simulation studies of GBF [21,24,29], the grid sizes of 1, 2, 3, 4, and 5 mm are chosen for study here. Figure 2a indicates that both the pressure drop and dust-removal efficiency vary less for the grid sizes of 1, 2, 3, 4, and 5 mm. In large installations, large grid sizes such as 5 mm can be chosen for the simulation. As the minimum bed height is 20 mm in this study, a 2 mm grid is used for subsequent calculations to obtain more information in the granular bed.
In order to simplify the computational complexity, some assumptions are made in the simulation, e.g., the bed height remains unchanged with radius and the Sauter mean diameter of the granules and dust are used. According to the experimental conditions [29], the superficial gas velocity is set as 0.345 m/s in the simulation; the dust diameter and concentration are 2.24 μm and 6985 mg/m3, respectively; the granular diameter is 10 mm; and the initial bed voidage is 0.49. The simulation results of pressure drop and dust-removal efficiency are in good agreement with the experimental results under different bed heights [29], for which the relative error is less than ±10.0%, as shown in Figure 2b. It indicates that the mesoscale simulation method can calculate the flow field in GBF in a simple and accurate way.

3. Results and Discussion

In order to further validate the accuracy of this mesoscale simulation method, the influences of different conditions are also studied by this method. The variation trend in the pressure drop and the dust-removal efficiency obtained in this paper are compared with previous studies. In particular, a single-variable method is adopted to investigate the effects of different conditions, e.g., bed height, superficial gas velocity, dust concentration, dust diameter, granular diameter, bed voidage, and filtration time, on the flow field. A certain parameter value in the baseline working condition is varied with reference to the related studies of GBF.
As shown in Figure 3 and Figure 4, the pressure and dust-mass fraction vary less along the radial direction and more along the axial direction. Compared to other parts of the bed, the pressure and dust-mass fraction vary more in the granular bed, and the variation varies with different conditions. In order to quantitatively analyze the effects of different conditions on pressure drop and dust-removal efficiency, the average pressure difference between the gas inlet and outlet sections is defined as the pressure drop of the bed; the reduction ratio of dust concentration between the two sections is defined as the dust-removal efficiency of the bed.

3.1. Flow-Field Distribution under Different Operating Conditions

The operating conditions mainly include the initial granular bed height L and superficial gas velocity ug, which will significantly affect the pressure drop and dust-removal efficiency. According to the relevant studies in GBF [29], L and ug change in the range of 20–100 mm and 0.15–0.55 m/s, respectively. In particular, it is better to display and compare the variation in results between the ratios of the results in the working condition and those in the baseline working condition. The results in the baseline working condition are marked as * in Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11.

3.1.1. Flow-Field Distribution under Different Bed Heights

When the bed height L increases from 20 mm to 100 mm gradually, the pressure in the granular bed varies less in the axial direction, while the pressure drop increases linearly, as shown in Figure 2 and Figure 5a. This is because when L increases, it does not affect the gas flow cross-sectional area of the granular bed. The gas–solid slip velocity does not change nor does the gas resistance of the granular bed with per unit thickness. Thus, the inertial and viscous resistance terms do not change (Figure 5a), and the pressure drop is approximately linear with L. The conditions and methods in the relevant references are summarized in Appendix A (Table A1).
The dust-removal efficiency increases from 17.55% to 62.01% gradually with the increase in L. The dust concentration in the gas phase within the granular bed, near the gas inlet side, changes little with the increase in L, as shown in Figure 4 and Figure 5b. The dust-removal efficiency is related to inertial collision, gravitational sedimentation, and adsorption, etc. As L increases, the probability of inertial collision between dust and filter granules increases significantly. As indicated by Figure 5b, when L grows, Re remains unchanged, as do the Nst and effective Stokes number Nsteff. Meanwhile, the number of granular units N increases, which effectively improves the dust-removal efficiency.
E = 1 ( 1 E 0 ) N , N = L / l

3.1.2. Different Superficial Gas Velocities

When the superficial gas velocity ug increases from 0.15 m/s to 0.55 m/s gradually, the pressure inside the granular bed increases in the axial direction. The pressure drop increases with increasing ug, whose change trend is consistent with the work described in [21,22,24,29,31], as shown in Figure 3 and Figure 6a. The conditions and methods in the relevant references are summarized in the Appendix (Table A1). This is due to the fact that as ug increases, the gas–solid slip velocity increases in the granular bed. The inertial force of the gas grows, as well as the gas resistance of filtered granules in granular bed, which causes the pressure drop to increase. According to Equations (5)–(7) and Figure 6a, the viscous resistance of the filtered granules to the gas is proportional to ug, and the inertial resistance is proportional to ug2. The pressure drop is then approximately proportional to 1–2 times ug.
With the increase in ug, the dust-removal efficiency increases; while the increase rate decreases, as shown in Figure 4 and Figure 6b. This is due to the fact that the dust-removal efficiency is related to inertial collision, gravitational settling, diffusion, and interception, etc. The superficial gas velocity is generally recommended to be 0–0.55 m/s in GBF [21,22,24,30]. When ug grows, the probability of inertial collision between dust and filter granules increases significantly. Moreover, as shown in Figure 6b, Re, Nst, and Nsteff increase, which effectively improves the dust-removal efficiency.

3.2. Flow-Field Distribution under Different Dust and Filter Granular Parameters

The characteristics of the flow field in a GBF are closely related to the forces between the dust and filter granules. For this reason, the effect of different dust concentrations, particle diameters of dust, granular diameters, and bed voidages on the flow field must be investigated.

3.2.1. Different Dust Concentrations

With the increase in dust concentration cp, the pressure variation and pressure drop in the granular bed along the axial direction increase gradually, which is consistent with the work in [32,33], as shown in Figure 3 and Figure 7. The conditions and methods in the relevant references are summarized in the Appendix (Table A1). As indicated in Figure 7a, with the increase in cp, the overall density of dusty gas increases, as well as the inertial and viscous resistance of filter granules to the gas. The pressure drop grows accordingly.
The dust-removal efficiency decreases gradually with increasing cp, as shown in Figure 4 and Figure 7b. As the decrease in gas velocity during dust filtration is not taken into account, the dimensionless numbers of Re, Nst, and Nsteff, which are computed by theoretical methods, vary less. However, during the actual filtration process, the total amount of dusty gas decreases for the decrease in dust and the gas velocity decreases. According to the analysis in Section 3.1.2, the probability of inertial collision between the dust and filter granules decreases at low ug, which causes the dust-removal efficiency to decrease.

3.2.2. Different Dust Diameters

When the dust diameter dp increases from 1 to 10 μm gradually, the pressure variation in the granular bed are basically unchanged along the axial direction as well as the pressure drop, which is consistent with the work in [21,22,24,29], as shown in Figure 3 and Figure 8a. The conditions and methods in the relevant references are summarized in the Appendix (Table A1). This is because the dust concentration is low, and thus the effect of dust diameter on the pressure drop seems to be small.
With the increase in dp, the dust-removal efficiency increases gradually, while the increase rate decreases, as shown in Figure 4 and Figure 8b. With the increase in dp, the inertial collision and interception probability of dust and filter granules decrease significantly. As shown in Figure 8b, when dp grows, Nst increases, as well as Nsteff and the equivalent granular diameter ratio Nr. The dust-removal efficiency increases afterwards.

3.2.3. Different Granular Diameters

When the granular diameter dg increases from 1 to 20 mm gradually, the pressure inside the granular bed changes along the axial direction, and the pressure drop of GBF decreases. Moreover, the decrease amplitude of pressure drop becomes smaller and smaller, which is consistent with the references [24,29], as shown in Figure 3 and Figure 9a. The conditions and methods in relevant references are summarized in the Appendix (Table A1). With the increase in dg, the contact area between the gas and filter granules in the granular bed decreases. According to Equation (21), it is inversely proportional to dg. At this time, Figure 9a shows that the viscous and inertial resistance of the filter granules to the gas phase decreases, which causes the pressure drop to decrease.
S = S V = n S 0 V = n π d g 2 V = ( 1 ε ) V 1 / 6 π d g 3 π d g 2 V = 6 1 ε d g
When dg increases, the dust-removal efficiency decreases gradually, as shown in Figure 4 and Figure 9b. When dg grows, the probability of inertial collision and interception of dust and filter granules decreases significantly. Meanwhile, as indicated by Figure 9b, the decrease in Nst, Nsteff, and Nr is greater than the increase in Re, which causes the dust-removal efficiency to decrease.

3.2.4. Different Initial Bed Voidages

When the initial bed voidage ε0 increases from 0.40 to 0.55 gradually, the pressure drop decreases and the rate of decrease becomes slower and slower, which is consistent with the work in [34], which can be seen in Figure 3 and Figure 10a. The conditions and methods in the relevant references are summarized in the Appendix (Table A1). With the increase in ε0, the gaps between filter granules become larger, the contact area between the gas and filter granules in the granular bed decreases (Equation (21)). Thus, the viscous and inertial resistance of the filter granules to gas decreases, as shown in Figure 10a. According to Equations (4)–(6), as the viscous and inertial resistances are inversely proportional to 1–2 and 2–3 times ε0, respectively, the pressure drop is approximately inversely proportional to 1–3 times ε0.
With the increase in ε0, the dust-removal efficiency decreases gradually, and the decrease rate becomes smaller and smaller, which can be seen in Figure 4 and Figure 10b. As ε0 increases, the contact area between the gas and filter granules decreases. Meanwhile, as Figure 10b indicates, the probability of inertial collision and interception between the dust and filter granules decreases, as well as Nsteff, and the dust-removal efficiency decreases.

3.3. Flow-Field Distribution under Different Filtration Times

The spatial distribution of dust deposition varies with filtration time t, which affects the distribution of the flow field in GBF. In this study, the average dust deposition in the granular bed is chosen in the theoretical calculation, which is used to compute the viscous and inertial resistance, Re, and other dimensionless numbers. The effect of t on pressure drop and dust-removal efficiency is then analyzed.
When t increases from 0 to 120 s gradually, the pressure changes in the axial direction and pressure drop grow in the granular bed, which is consistent with the work in [21], as shown in Figure 3 and Figure 11. The conditions and methods in the relevant references are summarized in the Appendix (Table A1). This is due to the fact that the amount of dust deposited in the granular bed increases with increases in t. The actual bed voidage decreases, which causes the viscosity and inertial resistances of the filter granules to the gas phase to increase. The pressure drop increases accordingly.
With the increase in t, the dust-removal efficiency increases gradually, and the increase rate becomes smaller and smaller, as shown in Figure 4 and Figure 11. With the increase in dust deposition, the actual bed voidage of the granular bed decreases. The probability of inertial collision and interception of the dust and filter granules also increases. Moreover, Nsteff grows, as well as the number of filter unit layers N and the dust-removal efficiency.

4. Conclusions

In this study, a mesoscale simulation method is proposed to simulate the flow field in a granular bed filter (GBF). After grid-independent and experimental validation, the effects of different conditions, i.e., operating conditions, dust and filter granular parameters, and filtration time, on the flow field inside the filter are investigated. The influence of different conditions coincides well with flow field in relevant studies, which further demonstrates the accuracy of the mesoscale simulation method. The following main conclusions are also obtained:
(1)
The mesoscale simulation method incorporates the macroscopic calculation models of pressure drop and dust-removal efficiency into the porous media model and the source term of the continuity equations. The accurate simulation results can then be obtained with a large grid size and small computational effort. The trends of pressure drop and dust-removal efficiency with different conditions are consistent with the results reported.
(2)
Different operating conditions have a large impact on the flow field in GBF. With the increase in bed height and superficial gas velocity, the gas residence time in the granular bed grows; the viscous and inertial resistance increases, as well as the pressure drop of the bed. Meanwhile, the inertial collision between dust and filter granules are enhanced. As the Reynolds number Re increases, as well as the number of granular unit N and effective Stokes number Nsteff, the dust-removal efficiency increases.
(3)
Dust and filter granular parameters also affect the flow-field distribution in GBF. With the increase in dust concentration or the decrease in granular diameter and bed voidage, the density of dusty gas increases, as well as the contact area between the gas and the filter granules. The pressure drop increases afterwards. At a low dust concentration, the pressure drop has little correlation with the dust diameter. With the increase in dust diameter or the decrease in dust concentration, granular diameter, and bed voidage, the inertial force increases, as well as the actual gas velocity and the contact area S’. Moreover, many dimensional numbers, e.g., Re, Nsteff, and Nr change. The dust-removal efficiency increases accordingly.
(4)
The dust deposition in the fixed GBF increases over time, which causes the pressure drop, and the gas viscous and inertial resistances to increase. In addition, as the actual bed voidage decreases with the increase in dust deposition, the probability of inertial collision and interception between the dust and filter granules increases, which causes the Nsteff, N, and dust-removal efficiency to grow.

Author Contributions

Conceptualization, T.L. and Z.Z.; methodology, T.L. and R.W.; software, T.L. and M.T.; validation, T.L.; formal analysis, T.L.; investigation, T.L.; resources, T.L.; data curation, T.L. and M.T.; writing—original draft preparation, T.L.; writing—review and editing, T.L. and R.W.; visualization, T.L.; supervision, S.Z.; project administration, D.W. and R.W.; funding acquisition, R.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Natural Science Foundation of Hebei Province: (B2021202022, B2022202003). The APC was funded by Natural Science Foundation of Hebei Province.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the Natural Science Foundation of Hebei Province (B2021202022, B2022202003).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a1, a2, a3, a4constantdimensionless
b1, b2, b3constantdimensionless
C s Cunningham correction factordimensionless
1/αcoefficient of viscous resistancedimensionless
C 2 coefficient of inertia resistancedimensionless
dpdust diameter μm
dggranular diametermm
eremoval efficiency in eachunit bed elementdimensionless
Etotal removal efficiencydimensionless
F ¯ average relative filter coefficientdimensionless
gacceleration of gravitym/s2
L granular bed heightmm
cpdust concentration mg/m3
N S t Stokes numberdimensionless
N s t e f f effective Stokes numberdimensionless
ReReynolds numberdimensionless
N R equivalent granular diameter ratiodimensionless
Δppressure dropPa
ppressurePa
u g superficial gas velocitym/s
v velocity vectorm/s
Ymass fractiondimensionless
γ adhesion probabilitydimensionless
ε bed void ratiodimensionless
ε 0 initial bed void ratiodimensionless
η 0 efficiency of the individual collectorsdimensionless
μ dynamic viscosity of the fluidPa·s
ρ densitykg/m3
tfiltration times
σ ¯ m average mass specific depositkg/m3
Subscripts
ispecies in the gas phase
ggas in the gas phase (except dg)
sdust in the gas phase
ssdust accumulated on the surface of filter particles
0indicates the initial state, i.e., t = 0 s

Appendix A

Table A1. Relevant investigations on the flow field in a granular bed.
Table A1. Relevant investigations on the flow field in a granular bed.
AuthorsBed L* (m)ug* (m/s)cp*dp*(μm)dg*(mm)ε*t*(s)Methods
Minghao You et al. [21]Cylindrincal0.40.2–0.810 g/m31–1001---0–3600Sim.
Junlin Chen et al. [22]Axial 0.02–0.050.25–0.550.108 g/min3–50100.40---Sim.
FeiLong Wang et al. [24]Rectangular0.030.2–1.02 g/m31–51–5------Sim.
L. Guan et al. [29]Axial0.02–0.100.15–0.553.281 m3/h1–215–200.485---Exp. and Sim.
Shaowu Yin et al. [31]Rectangular0.03, 0.040.3–0.70.25 g/m334.73, 50.502---Exp.
Yinsheng Yu et al. [32]Rectangular0.030.1, 0.50–3.6 g/m35, 203, 8------Exp. and Sim.
Ming Chang et al. [33]Cylinder2.579.065.95–59.13 g/m3102.070.37---Exp. and Sim.
T.E. Bustnes et al. [34]Axial0.0320.0127120 g/m320, 30, 400.40.917–0.98---Exp. and Sim.
Note: L, ug, cp, dp, dg, ε, and t are the bed height, superficial gas velocity, dust concentration, dust diameter, granular diameter, bed voidage, and filtration time, respectively.

References

  1. Hoffmann, B.; Boogaard, H.; de Nazelle, A.; Andersen, Z.J.; Abramson, M.; Brauer, M.; Brunekreef, B.; Forastiere, F.; Huang, W.; Kan, H.; et al. WHO air quality guidelines 2021–aiming for healthier air for all: A joint statement by medical, public health, scientific societies and patient representative organisations. Int. J. Public Health 2021, 66, 1604465. [Google Scholar] [CrossRef]
  2. Chen, Y.-S.; Hsiau, S.-S.; Lai, S.-C.; Chyou, Y.-P.; Li, H.-Y.; Hsu, C.-J. Filtration of dust particulates with a moving granular bed filter. J. Hazard. Mater. 2009, 171, 987–994. [Google Scholar] [CrossRef]
  3. Yu, Y.; Tao, Y.; Wang, F.-L.; Chen, X.; He, Y.-L. Filtration performance of the granular bed filter used for industrial flue gas purification: A review of simulation and experiment. Sep. Purif. Technol. 2020, 251, 117318. [Google Scholar] [CrossRef]
  4. Liu, S.; Tong, L.; Jiang, M.; Feng, X.; Yin, S.; Liu, C.; Wang, L.; Ding, Y. Flow stratification characteristics of binary particles in a moving granular bed. Powder Technol. 2020, 374, 482–491. [Google Scholar] [CrossRef]
  5. Zhao, J.; Huang, J.; Wu, J.; Fang, Y.; Wang, Y. Modeling and optimization of the moving granular bed for combined hot gas desulfurization and dust removal. Powder Technol. 2008, 180, 2–8. [Google Scholar] [CrossRef]
  6. Xiao, G.; Wang, X.; Zhang, J.; Ni, M.; Gao, X.; Luo, Z.; Cen, K. Granular bed filter: A promising technology for hot gas clean-up. Powder Technol. 2013, 244, 93–99. [Google Scholar] [CrossRef]
  7. Brown, R.C.; Shi, H.; Colver, G.; Soo, S.-C. Similitude study of a moving bed granular filter. Powder Technol. 2003, 138, 201–210. [Google Scholar] [CrossRef]
  8. Bémer, D. Granular bed filtration of a liquid aerosol. Powder Technol. 2022, 395, 218–225. [Google Scholar] [CrossRef]
  9. Yin, S.; Wang, X.; Guo, Y.; Wang, L. Filtration characteristics of granular bed with layered drawers for removing dust from gas streams. Particuology 2021, 55, 191–198. [Google Scholar] [CrossRef]
  10. Bradford, S.; Simunek, J.; Bettahar, M.; Van Genuchten, M.T.; Yates, S.R. Modeling colloid attachment, straining, and exclusion in saturated porous media. Environ. Sci. Technol. 2003, 37, 2242–2250. [Google Scholar] [CrossRef]
  11. Kretzschmar, R.; Barmettler, K.; Grolimund, D.; Yan, Y.-D.; Borkovec, M.; Sticher, H. Experimental determination of colloid deposition rates and collision efficiencies in natural porous media. Water Resour. Res. 1997, 33, 1129–1137. [Google Scholar] [CrossRef]
  12. Tardos, G.; Abuaf, N.; Gutfinger, C. Dust Deposition in Granular Bed Filters: Theories and Experiments. J. Air Pollut. Control. Assoc. 1978, 28, 354–363. [Google Scholar] [CrossRef] [Green Version]
  13. Jung, Y.; Walata, S.A.; Tien, C. Experimental Determination of the Initial Collection Efficiency of Granular Beds in the Inertial-Impaction-Dominated Region. Aerosol Sci. Technol. 1989, 11, 168–182. [Google Scholar] [CrossRef]
  14. Yu, Y.; Tao, Y.; Wang, F.-L.; He, Y.-L. Parameter study and optimization on filtration and resistance characteristics of granular bed filter. Adv. Powder Technol. 2018, 29, 3250–3256. [Google Scholar] [CrossRef]
  15. Wenzel, B.M.; Porciúncula, C.B.; Marcilio, N.R.; Menegolla, H.B.; Dornemann, G.M.; Godinho, M.; Martins, C.B. Filtration of dust in an intermittent moving granular bed filter: Performance and modeling. Sep. Purif. Technol. 2014, 133, 108–119. [Google Scholar] [CrossRef]
  16. Chen, J.L.; Li, X.F.; Wang, Y.W.; Huai, X.L. Numerical study on the dust removal characteristics in granular bed filter. J. Eng. Thermophys. 2018, 39, 1997–2004. [Google Scholar]
  17. Li, S.-Q.; Marshall, J. Discrete element simulation of micro-particle deposition on a cylindrical fiber in an array. J. Aerosol Sci. 2007, 38, 1031–1046. [Google Scholar] [CrossRef]
  18. Wang, F.-L.; He, Y.-L.; Tang, S.-Z.; Kulacki, F.A.; Tao, Y.-B. Real-time particle filtration of granular filters for hot gas clean-up. Fuel 2019, 237, 308–319. [Google Scholar] [CrossRef]
  19. Chou, C.; Tseng, C.; Smid, J.; Kuo, J.; Hsiau, S. Numerical simulation of flow patterns of disks in the asymmetric louvered-wall moving granular filter bed. Powder Technol. 2000, 110, 239–245. [Google Scholar] [CrossRef]
  20. Kloss, C.; Goniva, C.; Aichinger, G.; Pirker, S. Comprehensive DEM-CFD Simulations-Model Synthesis, Experimental Validation and Scalability. In Proceedings of the Seventh International Conference on CFD in the Minerals and Process Industries, Melbourne, Australia, 9–11 December 2009. [Google Scholar]
  21. You, M.; Li, Z.; Zhan, M.; Liu, M.; Sun, G.; Chen, J. Flow simulation and performance analysis of a cyclone–granular bed filter. Powder Technol. 2020, 361, 210–219. [Google Scholar] [CrossRef]
  22. Chen, J.; Li, X.; Huai, X.; Wang, Y.; Zhou, J. Numerical study of collection efficiency and heat-transfer characteristics of packed granular filter. Particuology 2019, 46, 75–82. [Google Scholar] [CrossRef]
  23. Guan, L.; Yuan, Z.; Gu, Z.; Yang, L.; Zhong, W.; Wu, Y.; Sun, S.; Gu, C. Numerical simulation of ash particle deposition characteristics on the granular surface of a randomly packed granular filter. Powder Technol. 2017, 314, 78–88. [Google Scholar] [CrossRef]
  24. Wang, F.-L.; Tang, S.-Z.; He, Y.-L.; Kulacki, F.A.; Tao, Y.-B. Parameter study of filtration characteristics of granular filters for hot gas clean-up. Powder Technol. 2019, 353, 267–275. [Google Scholar] [CrossRef]
  25. Elrahmani, A.; Al-Raoush, R.I.; Abugazia, H.; Seers, T. Pore-scale simulation of fine particles migration in porous media using coupled CFD-DEM. Powder Technol. 2022, 398, 117130. [Google Scholar] [CrossRef]
  26. Sobieskia, W.; Zhang, Q. Multi-scale modeling of flow resistance in granular porous media. Math. Comput. Simul. 2017, 132, 159–171. [Google Scholar] [CrossRef]
  27. Cheng, K.; Zhu, J.; Qian, F.; Cao, B.; Lu, J.; Han, Y. CFD–DEM simulation of particle deposition characteristics of pleated air filter media based on porous media model. Particuology 2023, 72, 37–48. [Google Scholar] [CrossRef]
  28. Tobisawa, R.Y.; de Lemos, M.J. Filtration efficiency of particle-laden flows for thermal plug and abandonment of oil wells using turbulence modeling in porous media. Int. Commun. Heat Mass Transf. 2022, 135, 106108. [Google Scholar] [CrossRef]
  29. Guan, L.; Gu, Z.; Yuan, Z.; Yang, L.; Zhong, W.; Wu, Y.; Sun, S. Numerical study on the penetration of ash particles in a three-dimensional randomly packed granular filter. Fuel 2016, 163, 122–128. [Google Scholar] [CrossRef]
  30. Manual, F. Fluent 6.3 User’s Guide; Fluent Inc.: New York, NY, USA, 2006. [Google Scholar]
  31. Yin, S.; Kang, P.; He, Y.; Wang, L.; Liu, C.; Tong, L.; Ding, Y. Experimental study on filtration characteristics of a novel moving granular bed filter. Sep. Purif. Technol. 2021, 267, 118624. [Google Scholar] [CrossRef]
  32. Yu, Y.; Tao, Y.; He, Y.; He, Y.-L. Structure optimization of granular bed filter for industrial flue gas filtration containing coagulative particles: An experimental and numerical study. Adv. Powder Technol. 2020, 31, 2244–2256. [Google Scholar] [CrossRef]
  33. Chang, M.; Fan, Y.; Lu, C. The inner/outer self-cleaning phenomenon in the granular bed filter-cyclone coupled separator. Sep. Purif. Technol. 2022, 287, 120460. [Google Scholar] [CrossRef]
  34. Bustnes, T.; Kaminski, C.; Mackley, M. The capture and release of biomass in a high voidage fibrous microstructure: Mechanisms and shear stress levels. J. Membr. Sci. 2006, 276, 208–220. [Google Scholar] [CrossRef]
Figure 1. Geometry model and mesh in an axial granular bed filter [29].
Figure 1. Geometry model and mesh in an axial granular bed filter [29].
Processes 11 01146 g001
Figure 2. Pressure drop and dust-removal efficiency under different meshing sizes, as well as bed heights, when filtration time t = 0 s [29].
Figure 2. Pressure drop and dust-removal efficiency under different meshing sizes, as well as bed heights, when filtration time t = 0 s [29].
Processes 11 01146 g002
Figure 3. Pressure drop under different conditions in an axial granular bed filter. (a) The baseline working condition: L = 20 mm, ug = 0.345 m/s, cp = 6985 mg/m3, dp = 2.24 μm, dg = 10 mm, ε = 0.49, t = 0 s; (b) bed height L = 40 mm; (c) superficial gas velocity ug = 0.55 m/s; (d) dust concentration cp = 10,000 mg/m3; (e) dust diameter dp = 10 μm; (f) granular diameter dg = 5 mm; (g) bed voidage ε = 0.55; (h) filtration time t = 120 s. Unspecified values of other parameters of (bh) are the same as those in (a) the baseline working condition.
Figure 3. Pressure drop under different conditions in an axial granular bed filter. (a) The baseline working condition: L = 20 mm, ug = 0.345 m/s, cp = 6985 mg/m3, dp = 2.24 μm, dg = 10 mm, ε = 0.49, t = 0 s; (b) bed height L = 40 mm; (c) superficial gas velocity ug = 0.55 m/s; (d) dust concentration cp = 10,000 mg/m3; (e) dust diameter dp = 10 μm; (f) granular diameter dg = 5 mm; (g) bed voidage ε = 0.55; (h) filtration time t = 120 s. Unspecified values of other parameters of (bh) are the same as those in (a) the baseline working condition.
Processes 11 01146 g003
Figure 4. Mass fraction of dust under different conditions in an axial granular bed filter. (a) The baseline working condition: L = 20 mm, ug = 0.345 m/s, cp = 6985 mg/m3, dp = 2.24 μm, dg = 10 mm, ε = 0.49, t = 0 s; (b) bed height L = 40 mm; (c) superficial gas velocity ug = 0.55 m/s; (d) dust concentration cp = 10,000 mg/m3; (e) dust diameter dp = 10 μm; (f) Granular diameter dg = 5 mm; (g) bed voidage ε = 0.55; (h) filtration time t= 120 s. Unspecified values of other parameters of (bh) are the same as those in (a) baseline working condition.
Figure 4. Mass fraction of dust under different conditions in an axial granular bed filter. (a) The baseline working condition: L = 20 mm, ug = 0.345 m/s, cp = 6985 mg/m3, dp = 2.24 μm, dg = 10 mm, ε = 0.49, t = 0 s; (b) bed height L = 40 mm; (c) superficial gas velocity ug = 0.55 m/s; (d) dust concentration cp = 10,000 mg/m3; (e) dust diameter dp = 10 μm; (f) Granular diameter dg = 5 mm; (g) bed voidage ε = 0.55; (h) filtration time t= 120 s. Unspecified values of other parameters of (bh) are the same as those in (a) baseline working condition.
Processes 11 01146 g004
Figure 5. Dimensionless numbers under different bed heights, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Figure 5. Dimensionless numbers under different bed heights, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Processes 11 01146 g005
Figure 6. Dimensionless numbers under different superficial gas velocities as well as pressure drop and dust-removal efficiency when filtration time t = 0 s.
Figure 6. Dimensionless numbers under different superficial gas velocities as well as pressure drop and dust-removal efficiency when filtration time t = 0 s.
Processes 11 01146 g006
Figure 7. Dimensionless numbers under different dust concentrations, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Figure 7. Dimensionless numbers under different dust concentrations, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Processes 11 01146 g007
Figure 8. Dimensionless numbers under different particle diameters of dust, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Figure 8. Dimensionless numbers under different particle diameters of dust, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Processes 11 01146 g008
Figure 9. Dimensionless numbers under different granular diameters, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Figure 9. Dimensionless numbers under different granular diameters, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Processes 11 01146 g009
Figure 10. Dimensionless numbers under different initial bed voidages, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Figure 10. Dimensionless numbers under different initial bed voidages, as well as pressure drop and dust-removal efficiency, when filtration time t = 0 s.
Processes 11 01146 g010
Figure 11. Dimensionless numbers under different filtration times, as well as pressure drop and dust-removal efficiency.
Figure 11. Dimensionless numbers under different filtration times, as well as pressure drop and dust-removal efficiency.
Processes 11 01146 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, T.; Zhao, Z.; Wang, R.; Tang, M.; Wang, D.; Zhang, S. A Mesoscale Simulation Approach to Study the Flow Field in an Axial Granular Bed Filter. Processes 2023, 11, 1146. https://doi.org/10.3390/pr11041146

AMA Style

Liu T, Zhao Z, Wang R, Tang M, Wang D, Zhang S. A Mesoscale Simulation Approach to Study the Flow Field in an Axial Granular Bed Filter. Processes. 2023; 11(4):1146. https://doi.org/10.3390/pr11041146

Chicago/Turabian Style

Liu, Tao, Zhifeng Zhao, Ruojin Wang, Meng Tang, Dewu Wang, and Shaofeng Zhang. 2023. "A Mesoscale Simulation Approach to Study the Flow Field in an Axial Granular Bed Filter" Processes 11, no. 4: 1146. https://doi.org/10.3390/pr11041146

APA Style

Liu, T., Zhao, Z., Wang, R., Tang, M., Wang, D., & Zhang, S. (2023). A Mesoscale Simulation Approach to Study the Flow Field in an Axial Granular Bed Filter. Processes, 11(4), 1146. https://doi.org/10.3390/pr11041146

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop