Abstract
In this article, the sideentry agitator commonly used in the wet desulfurization absorption tower of a 300 MW thermal power unit was taken as the research object, and the gas–liquid twophase flow field in the stirring tank was numerically simulated by using the twofluid model based on Euler–Euler method. The vertical and horizontal velocity field distribution and gas phase distribution in the stirring tank were studied. The results showed that the sideentry agitator formed a “fountain” circulating flow in the stirring tank. The highspeed region of the liquid flow was mainly concentrated near the mixing blade and the bottom banded region. The range of large changes in the fluid velocity was mainly concentrated from the bottom of the mixing tank to 2–3 m from the bottom. The distribution of the oxidizing air was concentrated in the middle of the stirring tank. The mixing effect of two phases in the stirred tank was good, and the oxidizing air was fully diffused. The numerical study in this article has reference significance for indepth understanding of the flow field and mixing effect of the sideentry agitator.
1 Introduction
The rapid development of the electric power industry has also brought about problems such as environmental pollution. Thermal power plants have always been a major emitter of greenhouse gases [1], especially the emission of SO_{2} has caused a bad impact on the atmospheric environment. In order to effectively control the emission of air pollutants in the thermal power industry, China has carried out ultralow emission transformation for the thermal power industry, among which the desulfurization of thermal power plants has become the key to control air pollution [2]. Limestone wet flue gas desulfurization is a relatively mature and widely used desulfurization technology [3,4]. The limestone slurry is used for countercurrent washing of the flue gas in the absorption tower to absorb sulfur dioxide in the flue gas, and undergo a series of reactions and dehydration treatments. A powdery gypsum with water content of about 10% is produced. In the absorption tower, the sideentry agitator located in the slurry tank is an important component to maintain the normal operation of the wet flue gas desulfurization system, and it is also the most commonly used stirring method [5]. It has the functions of precipitation, strengthening oxidation air diffusion, and limestone dissolution.
The sideentry agitator of the absorption tower is a kind of mechanical agitation. The liquid is propelled by the stirring blades to generate three basic flows: axial flow, swirling flow, and turbulent flow. Decompose these three flows into an axial part (axial flow) and a nonaxial part (rotational or helical flow), as shown in Figure 1. The spiral highspeed axial flow enables the liquid in the lowspeed flow area to be carried into the highspeed liquid flow area, so as to play a mixing role, and the nonaxial flow formed by the normal thrust will cause the liquid to swell up and down, so that the slurry pool can achieve uniform mixing and prevent the accumulation of sediment at the bottom of the tank.
The stirring performance and flow characteristics of sideentering agitators have been studied by some researchers. For the study of twophase flow, many scholars use computational fluid dynamics (CFD) method for numerical simulation to explore the characteristics of internal twophase flow field. Gómez et al. [7] investigated the flow in a reduced size rectangular tank equipped with a sideentering axial flow impeller in the laminar regime (
In order to further study the flow characteristics of the internal flow field of the sideentry agitator, numerical simulation of the flow field in the stirring tank with sideentry agitator was carried out, and the internal flow field velocity distribution and gas phase distribution were analyzed.
2 Numerical simulation method of sideentering agitator
2.1 Simulation working condition of sideentering stirring device
The numerical calculation model in this article was simplified from the actual structure size of the slurry pool at the bottom of the absorption tower for wet desulfurization of 300 MW thermal power unit, as shown in Figure 2. The stirring tank had a flat bottom with a diameter of T = 11.63 m and a liquid level of H = 8.79 m. At the bottom of the stirring tank, four sideentry agitators were evenly arranged along the circumference of the tank. The agitator blades were installed in the slurry pool of stirring tank, and the stirring shaft formed a certain inclination angle with the horizontal line and the center line of the tank, respectively, which is marked as β = 10° and θ = 7°.
The paddle type is a four narrow blade detachable rotary paddle agitator [17], as shown in Figure 3. The agitator consists of four blades spaced 90° apart. In order to enhance its axial flow performance, the blade placement angle was β = 45°. The diameter of the paddle is D = 1 m, the width of the paddle is W = 0.2 m, and the thickness of the paddle is δ = 0.012 m. The ratio of tank diameter to paddle diameter is T:D ≈ 9:1. The stirring shaft diameter is d = 0.1 m, and the extension length is l = 0.8 m. The lateral installation height of the agitator is h = 1.41 m and the speed of the agitator paddle is 194 rpm. The stirring effect requires that when one agitator is not working, the inlet of the pump cannot be blocked, the oxidized air diffuses evenly, and the slurry is stirred evenly and suspended far from the bottom of the tank. The stirring slurry is limestone slurry with a suspended solid volume fraction of 15%, density of ρ = 1,150 kg/m^{3}, viscosity of μ = 2 × 10^{−3} Pa s, particle diameter of 4 × 10^{−5} m, and pH between –5 and –6.
2.2 Governing equation of sideentering stirring device
In this article, a twofluid model based on Euler–Euler method was used to simulate the gas–liquid twophase flow field in stirring tank. The Euler–Euler method regards the dispersed phase and fluid such as bubbles and particles as continuous medium. The dispersed phase and the continuous phase coexist and penetrate each other. Both phases are processed in the Euler coordinate system, so they have the same form of governing equations. This model requires less calculation and is the main simulation method for multiphase systems at present [18–20]. The basic equations of the model mainly include the continuity equation and the momentum conservation equation. This article focused on the flow characteristics in the stirring tank, so the temperature change of gas–liquid twophase flow is not considered. The governing equation of the model is as follows:
The continuity equation [21,22]:
where
The momentum equation:
where p is the pressure acting on the two phases, τ _{eff,k } is the Reynolds stress under the average velocity gradient, R _{ k } is the momentum exchange term, and the interphase momentum exchange is mainly due to the action of the interphase force, which includes drag force, additional mass force, etc. F _{ k } is the volume force.
The main purpose of using the gas and liquid turbulence models is to solve the turbulent viscosity, and then obtain the effective viscosity of the fluid in the twofluid model, so that the twofluid model is closed. Generally, the k–ε turbulent model is used to calculate the continuousphase turbulent viscosity [23]. The standard k–ε twoequation turbulence model is a model close to the steady state of turbulent flow in the Reynolds timeaveraged method of turbulence calculation method. Considering the calculation cost and calculation accuracy, the standard k–ε turbulence model was used in this article [24], the transport equation of turbulent kinetic energy k and dissipation rate ε of turbulent kinetic energy is:
where u
_{
i
} denotes the component of fluid velocity along i direction, k denotes the turbulent kinetic energy, ε denotes the turbulent dissipation rate, μ
_{
t
} is the turbulent viscosity coefficient, and
2.3 Computing domain and grid division
Due to the geometric symmetry of the stirring tank, a quarter of the stirring tank was selected as the computational domain in order to simplify the calculation. For the sideentering stirring device, the computational domain was divided into two parts, as shown in Figure 4(a), due to the existence of rotating dynamic boundary and stationary static boundary. One part (a) contained the moving blades and the other part (b) contained the stationary baffle and tank. Two reference frames were used for the calculation of two regions. The rotating reference frame was used in the area where the blades were located, and the stationary reference frame was used for the other part. At the interface of two reference frames, the flow was assumed to be steady, i.e., the velocity at the interface must be the same for both reference frames. However, the cell node was static throughout the calculation. When using the above method to divide the dynamic and static areas, it must be ensured that the normal velocity of the nodes on the interface is zero, i.e., the top view of the interface must be circular rather than quadrilateral or other shapes.
In this article, the preprocessing software Gambit was used to generate grids. The computational domain was large and irregular, so the method of combining structured and unstructured grid was used to divide the grid into blocks, and the computational domain was divided into three parts: the rotating part of blades, the stationary region at the bottom of the stirring tank, and the stationary region at the upper part of the stirring tank. The stationary region at the upper part of the stirring tank was regular in shape, so the structured grid was divided by hexahedral mesh. The structure of the blade rotation area and the static area at the bottom of the stirring tank was more complex, so the adaptable tetrahedral unstructured mesh was used to divide the geometry. The method of combining unstructured grids of different sizes was used to divide the rotating part of the blade and the static part at the bottom of the stirring tank into grids, respectively. The rotating part of the blade adopted a relatively fine tetrahedral grid, while the stationary part at the bottom of the stirred tank adopted a slightly larger tetrahedral grid. The rotating part of the blade, agitator shaft, and the stationary part at the bottom of stirring tank were meshed encrypted to increase the calculation accuracy. The total number of grids was about 970,000, and the distribution of grid nodes is shown in Figure 4(b).
2.4 Numerical calculation method of sideentering stirring device
Fluent software was used to simulate the twophase flow field in the stirring tank. The multireference frame (MRF) method was used to deal with the interaction between the moving blade and the stationary slot wall. The rotational coordinate system was adopted for the blade and its nearby fluid region, and the stationary coordinate system was adopted for other regions. The momentum and energy exchange in the two different regions were realized through the interface conversion. The steadystate implicit separation algorithm was used in the solution process.
The fluid in the moving region (i.e., blade region) was set to rotate at the same speed as the stirring blade due to the MRF method used for simulation, while the fluid in the stationary region (i.e., other regions in the tank) was static. The wall of stirring tank was defined as the boundary condition of static wall, while the agitator shaft and agitator paddle were defined as the boundary condition of moving wall. The stirring shaft was in the static fluid region and moved relative to the fluid in the region, the paddle was in the moving fluid region and moved at the same speed as the surrounding fluid, so it was stationary with respect to the fluid in the region. Because the agitator blades were far from the free liquid surface, there is almost no liquid movement at the free liquid surface. Therefore, the free liquid surface can be defined as a symmetric boundary condition, so as to ignore its influence on the whole mixing process. The rotating periodic boundary without pressure drop was used for a pair of quarter crosssections divided in the tank body, while the interface between the blade area and other areas in the tank was defined as Interface after the surface was bonded in the grid division module.
The finite volume method was used for calculation, the delay modified QUICK difference was used to discrete the continuous phase, and the secondorder upwind difference scheme was used for the diffusion phase. The pressure relaxation factor was set as 0.3, the momentum relaxation factor as 0.7, the kinetic energy relaxation factor as 0.8, and the dissipation rate relaxation factor as 0.8. SIMPLE algorithm was used for pressure–velocity coupling. Due to the use of unstructured grid, a secondorder upwind difference scheme was used to improve the calculation accuracy.
3 Analysis of numerical simulation results of sideentering agitator
3.1 Analysis of macroflow field simulation results
From the simulated flow field distribution, it can be seen that the highspeed zone of liquid flow was mainly concentrated near the agitator blade and the banded region formed below it, and the maximum velocity appeared near the agitator blade end, which was 10 m/s. The velocity near the blade area was large and decreased rapidly with the increase of the distance from the blade area. Figures 5 and 6 show the simulated macroscopic velocity fields in different planes. As shown in Figure 5, the rotating impeller produced an axial flow with a certain downward angle, which flowed to the center of the tank. Part of the fluid turned to flow along the tank bottom to the tank center after hitting the tank bottom. Because four agitators were evenly arranged in the tank, the four streams of fluid collided with each other in the center of the tank and formed a strong axial liquid flow, rising up to the liquid surface. Then they turned into a radial flow to the tank wall, and flowed down to the bottom of the tank along the tank wall, forming an obvious up and down circulating flow in the tank.
The horizontal crosssection flow field at the axial height z = 1.4 m is shown in Figure 6. Due to the 7° inclination angle to the right when the agitator was arranged, the axial flow generated by the four blades formed a counterclockwise circular flow around the tank axis on the horizontal plane. The flow in the center of the tank was mainly upward along the axial direction of the stirring tank, so the horizontal flow was not obvious. A small vortex in a clockwise direction was formed near the wall of the tank between each two blades, which was generated by the partial liquid flow outside the large circulation under the action of the adsorption force at the rear end of the blade. The flow pattern of the upper part of the tank was relatively simple, and the flow field was like a fountain. The fluid in the center of the tank flowed upward from the bottom of the tank to the liquid surface, while the fluid around the tank flowed down from the liquid surface to the bottom of the tank. The fluid flow velocity in the center of the tank was slightly larger than that in the surrounding tank wall, but the horizontal flow velocity was relatively slow overall after reaching a certain height, and the velocity was about 0.2–0.4 m/s.
3.2 Velocity distribution on the central axis of the stirring tank
The fluid velocity components in the stirring tank included axial velocity, radial velocity, and tangential velocity. In this article, the speed direction was specified as follows: the reference coordinate system was established with the center of the circle at the bottom of the stirring tank as the origin. The radial velocity was positive in the positive direction of the Xaxis, and negative in the opposite direction. The axial velocity was positive in the positive direction of the Zaxis, and negative in the opposite direction. Tangential velocity was positive in the positive direction of the Yaxis, and negative in the opposite direction. Figure 7 shows the combined velocity on the central axis of the stirring tank and the distribution of each velocity component. It can be seen from the figure that the main velocity on the central axis of the stirring tank was axial velocity, and its change trend was consistent with that of the combined velocity. The axial velocity increased from the bottom of the stirring tank with the increase in the height of liquid level, and reached the maximum value about 2 m from the bottom, then gradually decreased to 0 at the liquid level. There was a certain tangential velocity component near the center of the tank bottom, but it decreased rapidly with the increase of z value, indicating that there was vortex in the fluid near the center of the tank bottom, which was caused by the large cycle formed in the lower part of the tank.
3.3 Velocity distribution in the vertical plane of the stirring tank
The purpose of stirring is to make the material have speed, so the stirring effect was evaluated by analyzing the velocity distribution of the material flow field. Different x and z values were taken in the y = 0 plane to analyze the axial and radial velocity distributions, respectively. In the bottom and middle of the stirring tank where the flow field velocity changed sharply, the value interval was small. In the part where the flow field changed slightly, the value interval was large. The values of x were 0.3, 2.0, 3.5, and 5.0 m, and the values of z were 0.5, 1.5, 4.0, and 8.5 m. The axial velocity distribution in the plane y = 0 is shown in Figure 8, and the radial velocity distribution is shown in Figure 9.
It can be seen from Figure 8 that the range that fluid flow velocity changed largely and was mainly concentrated at the bottom of the stirring tank to 2–3 m away from the bottom, and the change trend was consistent. The velocity near the tank wall increased sharply at first, and then decreased after maintaining a section of highspeed zone. This kind of distribution was mainly caused by the axial flow with a certain downward inclination angle generated by the rotation of the stirring blade, which rushed to the bottom of the stirring tank and then developed in the radial negative direction of the stirring tank. The fluid velocity in the middle section of the stirring tank changed little, and it increased slightly near the liquid level. With the increase of x value, the flow velocity at the bottom of the stirring tank increased until it reached the maximum near the end of stirring blade. However, the flow velocity decreased gradually with the increase of x value from the middle of the stirring tank to the liquid level and increased slightly near the wall. The velocity near the axis of the stirring tank changed gently, which indicated that the fluid developed toward the liquid level after interacting with the other three axial flows in the center of the stirring tank and then flowed toward the wall of the stirring tank after passing through the liquid level. The fluid distribution presented a “fountain shape.” A circulating flow was formed on the vertical surface of the fluid stirring tank to achieve the purpose of stirring.
It can be seen from Figure 9 that at the bottom of the stirring tank, the axial velocity of the fluid increased sharply due to the rotation of the stirring blades. As the fluid flowed toward the center of the stirring tank, the axial velocity changed slightly when it decreased to a certain value, especially in the highspeed axial flow area, which also proved that the agitator blade mainly generated axial flow. The fluid in the middle of the stirring tank flowed upward in the center of the stirring tank and downward around the tank wall, and the flow velocity in the center of the stirring tank was greater than the flow velocity near the tank wall. The flow velocity from the tank wall to the center of the stirring tank decreased first and then increased. The fluid formed a cycle with high peripheral flow velocity and low internal flow velocity. The axial velocity near the liquid level changed little and the fluid near the liquid level was mainly a radial flow from the center of the stirring tank to the tank wall.
3.4 Gas phase distribution in stirring tank
Through the intermediate transition reaction of oxygen and sulfite in the flue gas, part of the calcium sulfite is converted to gypsum, which is chemically called calcium sulfate dihydrate. The remaining calcium sulfite in the slurry of the absorption tower is oxidized by the air pumped in by the oxidation fan to produce calcium sulfate. The reaction process is mainly completed by the reaction of hydrogen sulfite and oxygen. The second purpose of agitation is to enhance the diffusion of oxidized air, promote the oxidation of calcium sulfite, the growth of gypsum crystals, and the dissolution of limestone. Therefore, the more evenly the air diffuses, the better the oxidation reaction will be. The ventilation device was located 1.4 m away from the bottom of the stirring tank, placed 1.8 m away from the tank wall horizontally, with a diameter of 200 mm and a ventilation speed of 7 m/s. The volume distribution of gas phase on the horizontal plane of the stirring tank z = 1.4 m is shown in Figure 10, and the volume distribution of gas phase on the vertical plane is shown in Figure 11.
It can be seen from Figures 10 and 11 that the oxidized air flowed to the bottom of the stirring tank together with the fluid under the action of axial flow with a certain downward inclination angle after being ejected from the nozzle. It flowed upward after collision with the bottom of the stirring tank. At this time, it collided with other fluid flowing to the center of the stirring tank again. Part of the oxidized air moved upward along the central axis of the stirring tank under the dual action of buoyancy and fluid, and part of the oxidized air changed direction on the horizontal plane to form a bottom circulation. The oxidized air speed decreased gradually when it flowed from the nozzle to the center of the stirring tank, and part of the energy was transferred to the surrounding fluid, which plays a stirring effect at the same time. After the oxidation air collided in the center of the stirred tank, it flowed to the liquid surface and diffused along the center of the stirred tank to the tank wall after obtaining a radial velocity.
4 Conclusion
In this article, numerical simulation was carried out on the sideentering agitator used in desulfurization tower of power plant. The flow field distribution obtained from simulation shows that the highspeed zone of liquid flow is mainly concentrated near the agitator blade and the belt region formed below it, and the maximum flow rate is near the end of the agitator blade, which is 10 m/s. The four streams of fluid generated by the four stirring blades collide with each other in the center of the tank and form a strong axial liquid flow, which rises up to the liquid level, then turns into a radial flow to the wall of the surrounding tank, and flows down to the bottom of the tank along the wall. An obvious “fountain” up and down circulating flow is formed in the tank body to achieve the purpose of stirring. The intense agitation is beneficial to the formation of grain, improve the utilization rate of oxygen, and promote the mixing of slurry and prevent scaling.

Funding information: The authors acknowledge the support of the National Natural Science Foundation of China (No. 52176040, 52006126).

Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

Conflict of interest: The authors state no conflict of interest.

Data availability statement: The data that support the findings of this study are available from the corresponding author upon reasonable request.
References
[1] Liu HZ, Guo YX. Comprehensive discussion on control technology of sulfur dioxide pollution from coal burning in China. Energy Environ Prot. 2005;19(1):13–4.Search in Google Scholar
[2] Hao YR, Meng C, Liu KK. Analysis on the development of industrial air pollution control technology. China Resour Compr Util. 2020;38(6):132–4.Search in Google Scholar
[3] Sai JC, Wu SH, Wang HT, Qin YK. Present situation of Chinese FGD technology and localization of manufacture issue. Power Syst Eng. 2003;19(1):53–4.Search in Google Scholar
[4] Hao JM, Wang SX, Lu YQ. Technical manual of sulfur dioxide pollution control from coal burning. Beijing: Chemical Industry Press; 2001.Search in Google Scholar
[5] Cai XJ, Wu LZ. The construct of sideinserted stirrer. PetroChem Eq. 2001;30(B05):45–6.Search in Google Scholar
[6] National Chemical Equipment Design Technology Center Station Mixing Engineering Technology Committee. Mixing equipment. Beijing: Chemical Industry Press; 2019.Search in Google Scholar
[7] Gómez C, Bennington CPJ, Taghipour F. Investigation of the flow field in a rectangular vessel equipped with a sideentering agitator. J Fluid Eng – Trans ASME. 2010;132(5):051106.10.1115/1.4001575Search in Google Scholar
[8] Havryliv R, Kostiv I, Maystruk V. Using the computational fluid dynamic software to mixing process modeling in the industrial scale vessel with sidemounted agitator. 2020 10th International Conference on Advanced Computer Information Technologies (ACIT). Deggendorf, Germany: IEEE; 2020. p. 192–5.10.1109/ACIT49673.2020.9208986Search in Google Scholar
[9] Fathonah NN, Madhania S, Nurtono T, Winardi S. Simulation of scaleup criteria for sideentry stirred tank in laboratory scale. IOP Conference Series: Materials Science and Engineering. Vol. 1053. Issue 1; 2021. p. 012106.10.1088/1757899X/1053/1/012106Search in Google Scholar
[10] Fathonah NN, Nurtono T, Kusdianto K, Madhania S, Wahyudiono W, Winardi S. Numerical study of the hydrodynamic characteristics in an agitated tank with sideentry mixer: the effect of stirrer entry angle. Int J Technol. 2019;10(3):521–30.10.14716/ijtech.v10i3.2925Search in Google Scholar
[11] AlKouz W, Bendrer BAI, Aissa A, Almuhtady A, Jamshed W, Nisar KS, et al. Galerkin finite element analysis of magneto twophase nanofluid flowing in double wavy enclosure comprehending an adiabatic rotating cylinder. Sci RepUK. 2021;11(1):1–15.10.1038/s41598021958462Search in Google Scholar
[12] Manzoor R, Khalid A, Khan I, Islam SU, Baleanu D, Nisear KS. Numerical simulation of drag reduction on a square rod detached with two control rods at various gap spacing via lattice Boltzmann method. Symmetry. 2020;12(3):475.10.3390/sym12030475Search in Google Scholar
[13] Kehn RO. Comparing top entry versus side entry agitator performance in low viscosity blending. Can J Chem Eng. 2011;89(5):1059–67.10.1002/cjce.20580Search in Google Scholar
[14] Grenville RK, Giacomelli JJ, VanOmmeren GJ, Hastings CF, Walters MJ. Blending in above ground storage tanks with sideentering agitators. Chem Eng Res Des. 2018;137:395–402.10.1016/j.cherd.2018.07.017Search in Google Scholar
[15] Wesselingh JA. Mixing of liquids in cylindrical storage tanks with sideentering propellers. Chem Eng Sci. 1975;30(8):973–81.10.1016/00092509(75)800649Search in Google Scholar
[16] Fang J, Sang ZF, Yang BQ. Numerical simulation of threedimensional flow field in side entry agitator. Petro Mach. 2009;37(1):30–4.Search in Google Scholar
[17] Chen ZP, Zhang XW, Lin XH. Design and selection manual of stirring and mixing equipment. Beijing: Chemical Industry Press; 2004. p. 141–2Search in Google Scholar
[18] Lahey RT, Drew DA. On the development of multidimensional twofluid models for vapor/liquid twophase flows. Chem Eng Commun. 1992;118:125–39.10.1080/00986449208936090Search in Google Scholar
[19] Drew DA. Mathematical modeling of twophase flow. Annu Rev Fluid Mech. 1983;15:261–91.10.21236/ADA114535Search in Google Scholar
[20] Hua JS, Wang CH. Numerical simulation of bubble driven liquid flows. Chem Eng Sci. 2000;55(19):4159–73.10.1016/S00092509(00)000750Search in Google Scholar
[21] Wang JJ, Li LC, Gu XP, Feng LF. Progress on CFD simulation of gas–liquid twophase flow in stirred tank reactor. Process Equip Pip. 2012;49(1):1–4.Search in Google Scholar
[22] Chen J, Xiao WD. Gas–liquid flow dynamics simulation in sideentering stirred tank. CIESC J. 2013;64(7):2344–52.Search in Google Scholar
[23] Song YL. Numerical simulation of gas–liquid twophase flow in multilayer paddle stirring tank. M.Sc. thesis. Beijing, China: Beijing University of Chemical Technology; 2006.Search in Google Scholar
[24] Wang FJ. Computational fluid dynamics analysis principle and application of CFD software. Beijing: Tsinghua University Press; 2004.Search in Google Scholar
[25] Launder BE, Spalding DB. Lectures in mathematical models of turbulence. London: Academic Press; 1972.Search in Google Scholar
© 2022 Ruonan Fan et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.