Next Article in Journal
Impact of Formation Dip Angle and Wellbore Azimuth on Fracture Propagation for Shale Reservoir
Next Article in Special Issue
A Modelling-Based Framework for Carbon Emissions Calculation in Additive Manufacturing: A Stereolithography Case Study
Previous Article in Journal
Multi-Objective Optimization of Drilling GFRP Composites Using ANN Enhanced by Particle Swarm Algorithm
Previous Article in Special Issue
Effect of Organic Powders on Surface Quality in Abrasive Blasting Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CFD-DEM Simulation of Fast Fluidization of Fine Particles in a Micro Riser

1
School of Mathematics and Statistics, Chongqing Three Gorges University, Chongqing 404020, China
2
School of Mathematical Sciences, Yangzhou University, Yangzhou 225002, China
*
Author to whom correspondence should be addressed.
Processes 2023, 11(8), 2417; https://doi.org/10.3390/pr11082417
Submission received: 18 July 2023 / Revised: 29 July 2023 / Accepted: 4 August 2023 / Published: 11 August 2023
(This article belongs to the Special Issue 10th Anniversary of Processes: Women's Special Issue Series)

Abstract

:
In recent years, the discrete element method (DEM) has gradually been applied to the traditional fluidization simulation of fine particles in a micro fluidized bed (MFB). The application of DEM in the simulating fast fluidization of fine particles in MFB has not yet received attention. This article presents a drag model that relies on the surrounding environment of particles, namely the particle circumstance-dependent drag model or PCDD model. Fast fluidization in an MFB of fine particles is simulated using DEM based on the PCDD model. Simulations indicate that the local structure in an MFB exhibits particle aggregation, which is a natural property of fast fluidization, forming a structure where a continuous dilute phase and dispersed concentrated phase coexist. There exists a strong effect of solid back-mixing in an MFB, leading to relatively low outlet solid flux. The gas back-mixing effect is, however, not so distinct. The axial porosity shows a monotonically increasing distribution with the bed height but does not strictly follow the single exponential distribution. The solid volume fraction at the bottom of the bed is significantly lower than the correlated value in CFB. The axial heterogeneous distribution of the cross-sectional average porosity in the lower half of the bed is also weakened. The radial porosity shows a higher distribution pattern in the central region and a lower one in the sidewall region.

1. Introduction

Gas–solid fluidized beds are widely used in environmental, energy, chemical and other process engineering fields [1]. Numerical simulation has become an indispensable means of studying gas–solid two-phase fluid dynamics for the cognition of a fluidized bed. The discrete element method [2,3,4,5], or DEM, providing particle-level information has been a tool needed to guide process design and operation in the field of fluidization engineering. The particle and device sizes simulated with DEM show a wide range. With the development of computer technology, high-performance algorithms and modeling capabilities, DEM has been able to simulate laboratory-scale coarse particle reactors. This method is gradually recognized by the industry and will inevitably be widely used [6]. On the other hand, with the improvement of the model and the calculation accuracy, the particle size used in DEM simulation gradually turns refined. Yu’s research group once simulated and analyzed different fluidization types of particle systems with particles as small as 30 μm [7].
In order to distinguish the particles in the gas–solid system with practical significance, Geldart [8] divided the particles into four classes: A, B, C and D. This is a widely accepted particle classification method in academia. Class C particles or ultra-fine particles are usually considered difficult to be fluidized; Class D particles or ultra-coarse particles have been used in numerous DEM simulation studies but are limited to fewer operations such as spout and slugging in practice; Class B particles are fluidized materials with a relatively small size in previous DEM simulations but still classified as coarse particles; and Class A particles or fine particles have a more smaller size (30–100 μm) and density (less than 1400 kg/m3). Fine particles are rarely used in DEM simulation, but they are widely used in engineering. For example, the catalytic cracking catalyst (FCC) particles, as a typical fine fluidized material, show advantages and characteristics different from other types of particles, such as a high gas–solid mass transfer rate, high bed expansion ratio and high heat transfer rate; they significantly expand before bubbling. Moreover, their cohesion is stronger, and the van der Waals force should not be ignored in many cases.
Ye et al. [9] modeled the fluidization of class A particles using DEM in 2004. The next year, Potic et al. [10] from the same study group proposed the concept of a micro fluidized bed (MFB), which is far small than an experimental circulating fluidized bed (CFB). They also studied the micro fluidized bed experimentally and theoretically. Han et al. [11] indicated that MFB has the advantages of convenience, safety, low cost and high efficiency. As an ideal isothermal reactor, it can be used as a reaction analyzer. As a small high selectivity module reactor, it can be used to produce high-value fine chemicals. Since only a few coarse particles side by side can fill the whole bed diameter, it is inevitable to abandon the use of class D particles or even class B particles. One has to choose Class A particles as fluidized materials in MFB experiments and DEM simulations. Therefore, the DEM simulation of class A particle fluidization in an MFB may become a hot research field.
Due to the lack of experimental data and simulation research, it is difficult to find reference literature. The existing DEM simulation results are mainly limited to fixed, uniform expansion, bubbling, slugging and turbulence fluidized beds [12,13,14,15,16,17]. Taking additional research in the past two years as examples, Guo et al. [18] presented a study of the solid-like and fluid-like states in the homogeneous fluidization regime of Class A particles. Li et al. [19] also performed a detailed analysis of bed hydrodynamics for polydisperse gas–solid flow of Class A particles. Both of the recent works implemented 3-D simulations. It is currently imperative to carry out the DEM simulation study of unconventional fast fluidization in MFB, even if it is a 2-D simulation study.
One difficult problem is how to reasonably simulate the flow structures, especially the macroscopic structure in an MFB. The local flow structure of particle agglomeration could be well simulated in [20]. However, the gas–solid back-mixing behavior could not be reasonably simulated, leading to the macroscopic structure seriously not insistent with the well-known core-annulus theoretical model. Can DEM reasonably simulate the fast fluidization of fine particles in MFB? Looking at the problem, this paper proposes a drag force model considering the surrounding particle environment and tries to solve the problem by reasonably calculating the drag force.

2. Drag Model

The drag model formally adopts the drag formula of a single particle in the particle group
F d i = 1 8 π d p C d i ρ g | u i v i | ( u i v i ) ε i 2 .
where d p is the particle diameter, C d i is the apparent drag coefficient of the particle group around particle i , ρ g is the gas density, u i is the local gas velocity, v i is the local gas velocity and ε i is the local porosity. A gird-independent calculation method for ε i is given in the following.
A neighborhood of the target particle is firstly determined to consider the influence of each surrounding particle in the neighborhood on the target particle. The influence of each surrounding particle is distinguished according to the different distances from the surrounding particle. This grid-independent calculation method is similar to the kernel approximation method in smooth particle dynamics (SPH) [21]. In this way, Xu et al. [22] calculated ε i as
ε i = 1 j = 1 N i W ( r i r j , h ) 1 6 π d p 2
where N i is the number of particles in the neighborhood, h is the characteristic radius or smooth length of the neighborhood, r i r j is the distance between particle i and j and W ( d , h ) is the normalized kernel function which is calculated as
W ( d , h ) = 7 478 π h 2 ( 3 d / h ) 5 6 ( 2 d / h ) 5 + 15 ( 1 d / h ) 5 , 0 d / h < 1 , ( 3 d / h ) 5 6 ( 2 d / h ) 5 ,   1 d / h < 2 , ( 3 d / h ) 5 , 2 d / h < 3 , 0 ,   d / h 3 ,
However, the anomalies close to the boundary have been addressed for the above Equations (2) and (3). The local porosity close to the boundary is effectively improved by adding boundary virtual particles in a previous study [23]. Figure 1 shows the diagram of the neighborhood, real particles and virtual particles in the kernel approximation method.
In ref. [23], h is set to be 2.5 d p . In fact, the accurate tune is difficult even in SPH simulations. On the other hand, in Equation (2), the contribution of each particle j for the local solid volume fraction of particle i is W ( r i r j , h ) π d p 2 / 6 . However, this contribution is formally not three-dimensional but two-dimensional. Moreover, the contribution should not be fixed for different particle systems. Thirdly, the basis for the correction is that the global mean values of local porosity cannot deviate from the global mean values of grid porosity. This problem is solved by introducing a solid volume fraction multiplier as
ε i = 1 λ j = 1 N i W ( r i r j , h ) 1 6 π d p 2
In order to avoid parameter testing, λ is determined by the following equation as
i = 1 N j = 1 N i W ( r i r j , h ) λ 6 N π d p 2 = ε ts
where ε ts is the particle filling ratio of the bed, i.e., the total solid volume fraction and N is the total number of real particles that are randomly distributed throughout the bed.
Combined with the apparent drag coefficient proposed by Wen and Yu [24] and the single-particle drag coefficient proposed by Schiller and Naumann [25], the apparent drag coefficient is calculated as
C d i = 24 μ g ε i 3.7 ρ g | u g - v p | d p + 3.6 μ g 0.313 ε i 4.387 ( ρ g | u g - v p | d p ) 0.313
where ε i is calculated according to Equation (4). Thus, the particle circumstance-dependent drag (PCDD) model is constructed according to Equations (1) and (6).

3. Simulation Method

The Navier–Stokes equations of the gas phase motion are expressed as Equations (7) and (8),
( ε g ρ g ) t + ( ε g ρ g u ) = 0
( ε g ρ g u ) t + ( ε g ρ g u u ) = ε g p S p ( ε g τ g ) + ε g ρ g g
where ε g is the mean porosity, p is the pressure, u is the gas velocity, t is the time, τ g is the stress tensor and S p is the momentum exchange source term. S p is calculated as
S p = i = 1 N k A i F d i A V
where N k is the number of particles that overlap with grid k , A is the particle disc area, A i is the overlap area of particle i with grid k and V is the grid volume, which is calculated as if the grid had a thickness of d p .
The exact area fraction model is used [26] to calculate the 2-D porosity, and the following conversion formula [3] is used to convert ε 2 D to the 3-D porosity ε 3 D as
ε 3 D = 1 2 π 3 ( 1 ε 2 D ) 3 2
The finite volume method is used for discretizing the Navier–Stokes equations, using the consistent velocity inlet, pressure outlet and wall impenetrable conditions as boundary conditions and using the SIMPLER method [27] to solve the discretized algebraic equations.
The particle translation is described according to the Newton’s Second Law as
ρ p V i d v i d t = ρ p V p g + F d i + F v i + F c i V i Δ p i
where F c i is the sum of contact force on particle i , F v i is the sum of van der Waals force on particle i , V i is the volume of the particle, p i is the local pressure. The soft sphere method is used to handle the particle–particle and particle–wall collision. The mass of the wall is infinite. The drag force F d i is calculated using the PCDD model. The van der Waals force between particles is calculated as
F v i j = H a d p e i j 24 ( d i j d p ) 2
where H a is the Hamaker constant, e i j is the unit vector from particle i to particle j and d i j is the distance between particle i and particle j . The van der Waals force of the wall applied to particle i is calculated as
F v i w = H a d p e i w 12 ( d i w d p ) 2
where e i w is the unit vector from the particle to the wall, and d i w is the distance between the particle and the wall. Moreover, to calculate the non-contact van der Waals force according to Equations (12) and (13), one needs to set the cutoff distance H 0 , i.e., the lower bound of d i j d p and d i w d p . The current simulations take H 0 = 0.4 nm.
The particle rotation is described as
I d ω i d t = T c i
where ω i is the particle angular velocity, I is the inertia of the particle and T c i is the torque of the collision.
In the present two-dimensional simulations, the size of the micro riser is D × H = 2.5 mm × 40 mm. The number of grids discretizing the corresponding flow field is 10 × 160. The maximal particle number used is 11,452. Initially, 8230 real particles are set static and randomly distributed throughout the whole bed, which determines λ = 0.948 according to Equation (5). The particle feed way follows the import-and-export balance. Table 1 shows the other fixed parameters of the gas–solid phases in the simulations.

4. Results and Discussion

4.1. Particle Agglomeration

The transient particle and porosity distribution are given in Figure 2. The coexisting structures of the dilute phase as disperse particles and the dense phase as granular aggregates are consistently presented in the MFB. Kuipers et al. [28] recommended using the porosity of 0.85 as a standard to distinguish the emulsion phase and the bubble phase in a bubbling fluidized bed. Figure 2 also uses this standard to distinguish the dense phase and the dilute phase. As can be seen from the porosity distribution, the dilute phase is continuous while the dense phase is dispersed. In different areas of the bed, the boundaries of the two phases are sometimes blurred and sometimes clear. With the change of time and space, clusters show the dynamic behaviors of continuous formation and fragmentation, and they vary in shape and size.

4.2. Gas–Solid Back-Mixing

Figure 3 presents the instantaneous distribution of the particle velocity within the 0.01 m height range in the bed. As shown in the figure, there is a serious back-mixing effect of the particles. The large cluster forms near the wall and then moves down along the wall. The radial velocity of the particles inside the cluster is small, and the axial velocity is large. Where the cluster passes, a strong tail vortex is formed to attract the tail particles to follow down. At the same time, a part of the peripheral particles begins to break away from the cluster due to the large radial velocity. Eventually, the cluster disintegrates due to the penetration of the fluid. At the bottom of the bed, the newly entering particles have not been fully accelerated, and the upper particles are in frequent contact. This results in a low radial and axial velocity of the particles in this area and also in a disorder orientation. In addition to the concentrated phase area, the particle velocity near the wall is relatively small, while the velocity in the central area is relatively large. Moreover, for the particles in the central area, whether they are the dispersed particles or particle groups, the axial velocity is mainly upward. As also indicated in Figure 3, the positive and negative velocity might be detected in almost all radial and axial positions.
Figure 4 shows the instantaneous distribution of the gas velocity in the 0.01 m height range, which corresponds to time in Figure 3. The gas tends to bypass the large cluster, around which the gas has a high radial velocity. In the interior of the cluster, both the radial and axial gas velocity are greatly reduced. Although the orientation of the radial velocity becomes chaotic, the larger part of the axial velocity is negative. Only in this area, the gas shows a sign of back-mixing. This proves that the back-mixing of the gas phase is mainly caused by the back-mixing of particles, i.e., the solid phase flow dominates. Except for the dense phase area, the gas velocity tends to be evenly distributed. Gas flow is mainly axially upward, and the radial velocity is small. That is to say, the dilute phase area is generally dominated by the upward gas flow. Thus, most of the dispersed particles can be continuously transported upward.
Figure 5 shows the outlet solid flux over time. After about 0.2 s, the two-phase flow within the MFB reaches a relatively stable state. In the relatively stable state, the outlet solid flux still has large fluctuations. The average value is about 7.8 kg/(m2·s), and it is much smaller than in the literature [20], which is 90~110 kg/(m2·s). This significant difference indicates that the overall drag force calculated with the PCDD model is largely reduced in the current simulations.
The saturated entrainment rate of the particle G s can be calculated using the following correlation formula of Bai and Kato [29] as
G s d p μ g = 0.125 F r 1.85 A r 0.63 ( ρ p ρ g ρ g ) 0.44
where F r is the Frodes number and A r is the Archimedes number. When the operating gas velocity is 1.52 m/s, the correlated G s is 16.18 kg/(m2·s). The experimental result given by Li and Kwauk [30] and the simulation result presented by Yang et al. [31] are both approximately 14.3 kg/(m2·s), indicating that the correlation formula has good predictive power. The correlation value at the current gas operating velocity 1.7 m/s is 19.9 kg/(m G s ·S). Therefore, due to the failure to reasonably simulate the back-mixing behavior of particles, the outlet solid flux is overestimated in [20], which is even much higher than the saturation entrainment rate of particles. The currently selected gas–solid properties are basically the same as those in Yang et al. In this work the operating gas velocity is slightly higher, and the particle filling ratio is slightly lower. However, the simulated average outlet solid flux is much lower than the simulated and correlated values in CFB. The authors speculate that the most direct reason might be the geometry effect. In MFB compared with CFB, the bed diameter is very small, and thus for the wall, the relative area of contact with particles significantly enlarges [11]. The strong wall friction force causes the particles near the wall to be unable to be transported upward, which through particle collisions, gradually passed to the surrounding particles and even the central area. And the gas phase needs to consume a lot of energy or pressure drop for suspended particles, which seriously hinders the transport of the particles.

4.3. Axial Structure

Figure 6 shows the distribution of the mean porosity at the bottom and top sections of the bed. It shows that the top of the bed is a dilute phase area, and the fluctuation range of the mean porosity is small or even close to 1. At the bottom of the bed, the porosity is greatly reduced, and the fluctuation range is large. The lowest value can descend to below 0.8, which indicates that the dense phase can continuously form in different radial positions of the bottom.
Figure 7 shows the instantaneous and time averaging cross-sectional porosity between 0.2–0.5 s, with the height 0.000125 m, 0.005125 m, 0.015125 m, 0.025125 m, 0.035125 m and 0.039875 m, respectively. In CFB, the axial porosity distribution of fast fluidized bed is usually a monotonically increasing exponential distribution unless the outlet solid flux reaches the saturated entrainment rate of the particles. The current simulated axial distribution, although a monotonically increasing function of height, does not strictly obey the simple exponential distribution. It appears to have an inflection point between 0.015125 m and 0.025125 m. Because the parameters of the simple exponential function given by the vast majority of experimental studies are not uniform, it is difficult to compare the current simulated distribution of axial porosity with the results of the literature. Consider the solid phase volume fraction at the bottom of the bed ε sd . As G s < G s ; the correlation formula of Bai and Kato [29] is used as
ε sd ρ p ( u g u t ) G s = 1 + 0.00614 ( ρ p u g G s ) 0.23 ( ρ p ρ g ρ g ) 1.21 ( u g g D ) 0.383
where u t represents the terminal velocity of the particles. In the moderate Reynolds number condition, the formula of Lewis et al. [32] can be used to calculate u t . From (16) the correlated ε sd is 0.22, while the presently simulated value is 0.136. The simulated value is far below the correlated value. The possible reasons for such a deviation include the following. The two-dimensional simulation cannot reach the reality of the three-dimensional simulation; the simple set import and export boundary conditions cannot simulate the real complex import and export effects. Additionally, it should be more inclined to think that the relative area of wall-particle contacts increase significantly in the MFB, resulting in a large decrease of the outlet solid flux. It is not difficult to analyze according to Formula (16), and ε sd decreases with the decrease in the outlet solid flux. Due to the low level of ε sd , the heterogeneous distribution of cross-sectional porosity in the bottom area of the bed is weakened. This may be the main reason why the simulated axial voidage does not strictly obey the simple exponential distribution.

4.4. Radial Structure

The radial flow law of gas–solid phases in the fast fluidized bed is extremely important. The results of extensive experimental analysis show that as long as the cross-sectional porosity is determined, the radial distribution of the porosity is only a function of the radial position, i.e., the dimensionless radius, regardless of the operating conditions. The cross-sectional porosity at the height 0.015125 m is 0.912. The correlation formula of Patience and Chaouki [33] can be used as the following
ε ¯ 0.4 ε ε ¯ 0.4 ε = 4 r 6 .
The correlation result of the radial porosity in CFB is then obtained.
Figure 8 shows the radial distribution of the porosity at the height 0.015125 m. r can be the dimensionless radius of the bed, and only the results from the right of the MFB can be selected. The simulation results in Figure 8 and the correlation results both show the high porosity in the central and the low porosity in close the wall, which follows the typical core-annulus distribution in the trend. The difference between the correlated and the simulated lies in the relatively low porosity in the central area and the relatively high porosity in the close wall in the simulated results, indicating a weak core-annulus effect. This significant difference lies in the amplification effect of wall friction in the MFB. The narrow bed diameter itself is easy to cause material slug, which can easily fill most of the bed diameter. As has also been stated, there exists the trace of clusters in the central area of the bed. Therefore, the heterogeneity of the radial distribution is also weakened.

5. Conclusions

In this paper, we give a drag model considering the particle environment, which is heterogeneous and complex. This model is implanted into CFD-DEM for simulation of fast fluidization in the MFB of fine particles. Simulations show that gas–solid two-phase flow in the MFB satisfies the common law of fast fluidization to some extent. Because the geometric size of the container is much smaller than that of a CFB, the relative area of the wall-particle contact is large. The wall friction factor becomes very significant, and thus, simulations exhibit a special rule of a large difference with those in CFB. The specific conclusions are drawn as follows:
(1)
The local structure in the MFB satisfies the natural property of fast fluidized particle agglomeration, forming a disperse dilute phase and continuous dense phase;
(2)
There is serious gas–solid back-mixing in the MFB, and the dense phase is the main area of gas–solid back-mixing. The wall friction factor aggravates the particle remixing effect, resulting in a relatively low outlet solid flux;
(3)
The axial porosity presents an increasing distribution with the bed height but does not strictly satisfy the monotonic exponential distribution. The solid volume fraction at the bottom of the bed is much lower than the correlated results for a CFB;
(4)
The radial porosity exhibits a weak core-annulus structure with a higher central region and lower side region. Compared with the correlated results for a CFB, the central region in the MFB is relatively dense, while the side region is relatively dilute;
(5)
All distinct variation results of the present simulation can be successfully explained using the relatively strong friction factor. Thus, the present drag model and simulation are both validated and effective, at least in the sense of a qualitative trend.

Author Contributions

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

Funding

This research received fundings supported by the Natural Science Foundation of Chongqing (CSTB2023NSCQ-MSX0974 and CSTB2022NSCQ-MSX0290) and the National Natural Science Foundation of China (12226352).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Nomenclature

Aarea, m2
ArArchimedes number
Cdrag coefficient
Dbed wide, m
dparticle diameter or distance between particles, m
eunit vector
Fforce on particle, and N
FrFred number
Goutlet solid flux, kg·m−2·s−1
ggravity acceleration, m·s−2
hsmooth length, m
Hheight of bed, m
HaHamaker constant, N·m
H0Truncation distance, m
Iinertia moment of particle as spherical, kg·m2
i, j, kparticle or grid index
Nnumber of particles
ppressure, Pa
rparticle position vector
rdimensionless radius
Spmomentum exchange source term
Ttorque, N·m
ttime, s
u0inlet gas velocity, m·s−1
ugas velocity, m·s−1
utparticle terminal speed
Vvolume, m3
vparticle velocity, m·s−1
ε porosity
ε ¯ cross-sectional porosity
ε sd solid volume fraction at bottom of bed
κ stiffness coefficient, N·m−1
λ solid volume fraction multiplier
μ viscosity, N·s·m−2
ρ density, kg·m−3
τ viscous stress tensor, Pa
ω particle angular velocity, s−1
ξ restitution coefficient
subscript
2Dtwo-dimension
3Dthree-dimension
ccontact
ddrag
ggas
i, j, kparticle or grid index
mfminimal fluidized state
pparticle
ssolid
ttotal
vvan der Waals
wbed wall

References

  1. Jin, Y.; Zhu, J.X.; Wang, Z.W.; Yu, Z.Q. Fluidization Engineering Principles; Tsinghua University Press: Beijing, China, 2001. [Google Scholar]
  2. Tsuji, Y.; Kawaguchi, T.; Tanake, T. Discrete particle simulation of two-dimensional fluidized bed. Powder Technol. 1993, 77, 79–87. [Google Scholar] [CrossRef]
  3. Hoomans, B.P.B.; Kuipers, J.A.M.; Briels, W.J.; Van Swaaij, W.P.M. Discrete particle simulation of bubble and slug formation in a two-dimensional gas-fluidised bed: A hard-sphere approach. Chem. Eng. Sci. 1996, 51, 99–108. [Google Scholar] [CrossRef] [Green Version]
  4. Xu, B.H.; Yu, A.B. Numerical simulation of the gas-solid flow in a fluidized bed by combing discrete particle method with computational fluid dynamics. Chem. Eng. Sci. 1997, 52, 2785–2809. [Google Scholar] [CrossRef]
  5. Ouyang, J.; Li, J.H. Particle-motion-resolved discrete model for simulating gas-solid fluidization. Chem. Eng. Sci. 1999, 54, 2077–2083. [Google Scholar] [CrossRef]
  6. Zhou, C.L.; Zhao, Y.Z. Discrete element method and its applications in fluidizatio. CIESC J. 2014, 65, 2520–2534. [Google Scholar]
  7. Yu, A.B.; Xu, B.H. Particle-scale modelling of gas–solid flow in fluidisation. J. Chem. Technol. Biotechnol. 2003, 78, 111–121. [Google Scholar] [CrossRef]
  8. Geldart, D. Types of gas fluidization. Powder Technol. 1973, 7, 285–292. [Google Scholar] [CrossRef]
  9. Ye, M.; van der Hoef, M.; Kuipers, J. A numerical study of fluidization behavior of Geldart A particles using a discrete particle model. Powder Technol. 2004, 139, 129–139. [Google Scholar] [CrossRef]
  10. Potic, B.; Kersten, S.R.A.; Ye, M.; van der Hoef, M.A.; Kuipers, J.A.M.; van Swaaij, W.P.M. Fluidization with hot compressed water in micro-reactors. Chem. Eng. Sci. 2005, 54, 5982–5990. [Google Scholar] [CrossRef]
  11. Han, Z.; Yue, J.; Geng, S.; Hu, D.; Liu, X.; Bello Suleiman, S.; Cui, Y.; Bai, D.; Xu, G. State-of-the-art hydrodynamics of gas-Solid micro fluidized beds. Chem. Eng. Sci. 2021, 232, 116345. [Google Scholar] [CrossRef]
  12. Ye, M.; van der Hoef, M.A.; Kuipers, J.A.M. The effects of particle and gas properties on the fluidization of Geldart A particles. Chem. Eng. Sci. 2005, 60, 4567–4580. [Google Scholar] [CrossRef]
  13. Liu, P.; LaMarche, C.Q.; Kellogg, K.M.; Hrenya, C.M. Fine-particle defluidization: Interaction between cohesion, Young’s modulus and static bed height. Chem. Eng. Sci. 2016, 145, 266–278. [Google Scholar] [CrossRef] [Green Version]
  14. Hou, Q.F.; Zhou, Z.Y.; Yu, A.B. Micromechanical modeling and analysis of different flow regimes in gas fluidization. Chem. Eng. Sci. 2012, 84, 449–468. [Google Scholar] [CrossRef]
  15. Kobayashi, T.; Tanaka, T.; Shimada, N.; Kawaguchi, T. DEM-CFD analysis of fluidization behavior of Geldart Group A particles using a dynamic adhesion force model. Powder Technol. 2013, 248, 143–152. [Google Scholar] [CrossRef]
  16. Li, T.; Rabha, S.; Verma, V.; Dietiker, J.F.; Xu, Y.; Lu, L.; Rogers, W.; Gopalan, B.; Breault, G.; Tucker, J.; et al. Experimental study and discrete element method simulation of Geldart Group A particles in a small-scale fluidized bed. Adv. Powder Technol. 2017, 28, 2961–2973. [Google Scholar] [CrossRef]
  17. Li, S.; Zhao, P.; Xu, J.; Zhang, L.; Wang, J. Direct comparison of CFD-DEM simulation and experimental measurement of Geldart A particles in a micro-fluidized bed. Chem. Eng. Sci. 2021, 64, 622–625. [Google Scholar] [CrossRef]
  18. Guo, Q.; Bordbar, A.; Ma, L.; Yu, Y.; Xu, S. A CFD-DEM study of the solid-like and fluid-like states in the homogeneous fluidization regime of Geldart A particles. AIChE J. 2022, 68, e17420. [Google Scholar] [CrossRef]
  19. Li, S.; Zhao, P.; Xu, J.; Zhang, L.; Wang, J. CFD-DEM simulation of polydisperse gas-solid flow of Geldart A particles in bubbling micro-fluidized beds. Chem. Eng. Sci. 2022, 253, 117551. [Google Scholar] [CrossRef]
  20. Wu, G.R.; Zuo, Z.F.; Li, Y.G.; Israr, M. Improvement of relative DEM time step range in fast fluidization simulation of Type-A FCC particles. Processes 2023, 11, 1155. [Google Scholar] [CrossRef]
  21. Lucy, L.B. Numerical approach to testing the fission hypothesis. Astron. J. 1977, 82, 1013–1024. [Google Scholar] [CrossRef]
  22. Xu, M.; Ge, W.; Li, J.H. A discrete particle model for particle–fluid flow with considerations of sub-grid structures. Chem. Eng. Sci. 2007, 62, 2302–2308. [Google Scholar] [CrossRef]
  23. Wu, G.R.; Ouyang, J.; Li, Q. Revised drag calculation method for coarse grid Lagrangian-Eulerian simulation of gas-solid bubbling fluidized bed. Powder Technol. 2013, 235, 959–967. [Google Scholar] [CrossRef]
  24. Wen, C.Y.; Yu, Y.H. Mechanics of fluidization. Chem. Eng. Progr. Symp. Ser. 1966, 62, 100–111. [Google Scholar]
  25. Schiller, V.L.; Naumann, A. Uber die grundlegenden berechnungen bei der schwerkraftaufbereitung. Z. Ver. Dtsch. Ing. 1993, 77, 318–320. [Google Scholar]
  26. Wu, G.R.; Ouyang, J. Fine grid DEM simulation of bed layer height in bubbling fluidized-bed. CIESC J. 2014, 65, 2092–2097. [Google Scholar]
  27. Patankar, T.V. Numerical Heat Transfer and Fluid Flow; Hemisphere Publishing Corporation: New York, NY, USA, 1980. [Google Scholar]
  28. Kuipers, J.A.M.; van Uuin, K.J.; van Beckum, F.P.H.; van Swaaij, W.P.M. A numerical model of gas-fluidized beds. Chem. Eng. Sci. 1992, 47, 1913–1924. [Google Scholar] [CrossRef] [Green Version]
  29. Bai, D.; Kato, K. Generalized correlations of solids holdups at dense and dilute regions of circulating fluidized beds. In Proceedings 7th SCEJ Symposium on CFB; Society of Chemical Engineers: Tokyo, Japan, 1994. [Google Scholar]
  30. Li, J.; Kwauk, M. Particle-Fluid Two-Phase Flow: The Energy-Minimization Multi-Scale Method; Metallurgical Industry Press: Beijing, China, 1994. [Google Scholar]
  31. Yang, N.; Wang, W.; Ge, W.; Li, J.H. CFD simulation of concurrent-up gas-fluid flow in circulating fluidized beds with structure-dependent drag coefficient. Chem. Eng. Sci. 2003, 96, 71–80. [Google Scholar] [CrossRef]
  32. Lewis, W.K.; Gilliland, E.R.; Bauer, W.C. Characteristics of fluidized particles. Ind. Eng. Chem. 1949, 41, 1104–1117. [Google Scholar] [CrossRef]
  33. Patience, G.S.; Chaouki, J. Solids Hydrodynamics in the Fully Developed Region of CFB Riser, 4th ed.; Large, J.F., Laguerie, C., Eds.; Engineering Foundation: New York, NY, USA, 1996. [Google Scholar]
Figure 1. Neighboring area, real particle and virtual particle in kernel estimator method.
Figure 1. Neighboring area, real particle and virtual particle in kernel estimator method.
Processes 11 02417 g001
Figure 2. Particle and porosity distribution.
Figure 2. Particle and porosity distribution.
Processes 11 02417 g002
Figure 3. Particle velocity distribution below 0.01 m.
Figure 3. Particle velocity distribution below 0.01 m.
Processes 11 02417 g003
Figure 4. Gas velocity distribution below 0.01 m.
Figure 4. Gas velocity distribution below 0.01 m.
Processes 11 02417 g004
Figure 5. Outlet solid flux.
Figure 5. Outlet solid flux.
Processes 11 02417 g005
Figure 6. Cross-sectional porosity at bottom and top.
Figure 6. Cross-sectional porosity at bottom and top.
Processes 11 02417 g006
Figure 7. Axial distribution of cross-sectional porosity.
Figure 7. Axial distribution of cross-sectional porosity.
Processes 11 02417 g007
Figure 8. Radial porosity distribution.
Figure 8. Radial porosity distribution.
Processes 11 02417 g008
Table 1. Fixed parameters for particle and gas.
Table 1. Fixed parameters for particle and gas.
ParticleGas
Density ρp = 930 kg·m−3Viscosity μg = 1.7 × 10−5 N·s·m−2
Particle diameter dp = 54 μmDensity ρg = 1.28 kg·m−3
Porosity at minimum fluidization εmf = 0.45Inlet gas velocity u0 = 1.7 m·s−1
Stiffness Coef. κ = 10 N·m−1CFD time step Δ tg = 2 × 10−6 s
Restitution Coef. ξ = 0.9
DEM time step Δ tp = 2.5 × 10−7 s
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

Wu, G.; Li, Q.; Zuo, Z. CFD-DEM Simulation of Fast Fluidization of Fine Particles in a Micro Riser. Processes 2023, 11, 2417. https://doi.org/10.3390/pr11082417

AMA Style

Wu G, Li Q, Zuo Z. CFD-DEM Simulation of Fast Fluidization of Fine Particles in a Micro Riser. Processes. 2023; 11(8):2417. https://doi.org/10.3390/pr11082417

Chicago/Turabian Style

Wu, Guorong, Qiang Li, and Zhanfei Zuo. 2023. "CFD-DEM Simulation of Fast Fluidization of Fine Particles in a Micro Riser" Processes 11, no. 8: 2417. https://doi.org/10.3390/pr11082417

APA Style

Wu, G., Li, Q., & Zuo, Z. (2023). CFD-DEM Simulation of Fast Fluidization of Fine Particles in a Micro Riser. Processes, 11(8), 2417. https://doi.org/10.3390/pr11082417

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