CFD-DEM simulation of fluidized bed with an immersed tube using a coarse-grain model
Graphical abstract
Introduction
Fluidized beds are frequently used in various industries due to their excellent heat and mass exchange efficiency (Asegehegn et al., 2011b, Goldschmidt et al., 2004, Liu et al., 2020a, Taghipour et al., 2005, Tsuji et al., 2008, Yue et al., 2019, Zhang et al., 2019). In many applications, immersed tubes are introduced into the fluidized beds to enhance the transfer and conversion rates, control the bed temperature and facilitate the particle mixing (Asegehegn et al., 2011a, Hou et al., 2016, Li et al., 2011, Schreiber et al., 2011, Wahyudi et al., 2016, Yurong et al., 2004). However, understanding the gas–solid flow properties is difficult ascribed to the complex underlying mechanism in the systems and shortcomings in research methods, which greatly impede their design and scale-up. The presence of tubes further complicates the gas-flow dynamics of the systems and then increases the difficulty of investigation (Asegehegn et al., 2011b). Therefore, an in-depth study of gas–solid behavior in the fluidized bed with immersed tubes has aroused the interest of researchers. In recent years, many experiments have been conducted to gain further insight into the various characteristics of fluidized bed with immersed tubes involving heat transfer (Kamble et al., 2014), bubble behavior (Asegehegn et al., 2011a, Asegehegn et al., 2011b), and erosion around the tube surface (Wiman and Almstedt, 1997). However, experimental measurements have some limitations, e.g., the harsh operational conditions, high expense, and restricted optical access (Ostermeier et al., 2019, Qi et al., 2019). Thanks to the rapid progress in available computational resources and numerical models, numerical simulation gains particular importance as a powerful tool for the prediction of gas–solid flow (Almuttahar and Taghipour, 2008, Liu et al., 2020b), e.g., CFD-DEM (computational fluid dynamic - discrete element method) and CFD-DPM (computational fluid dynamic - discrete phase method). As typical Eulerian-Lagrangian methods, CFD-DEM can track every particle individually taking into account the particle–particle and particle–wall interaction, and it has been widely used in the simulation of fluidized beds (Gan et al., 2016, Lu et al., 2014, Moliner et al., 2019, Su et al., 2011, Zhao et al., 2017); while CFD-DPM neglects the interaction between particles, and it is suitable for the simulation of dilute flow. To date, the industrial application of CFD-DEM with billions of particles is still limited because of excessive calculation cost even though great advances have been made in the computing capacity of computers. The calculation time, therefore, becomes extremely long to the extent that the simulation will last for months by a single personal computer. Hence, coarse-grain model is proposed to overcome exhaustive computational cost existing in the present DEM method because it can notably reduce the number of particles tracked in the numerical domain. The coarse-grain model combining the CFD-DEM method (CG CFD-DEM) has been proved to be effective by many researchers and applied to a variety of situations. Sakai and coworkers simulated the macroscopic behavior in three-dimensional fluidized beds numerically and experimentally using the CG model (Sakai et al., 2014). They simulated a group of original particles with large-sized particles, referred to as the CG particles. In their CG model, the total energy was estimated to agree between the CG particle system and the original particle system. They verified the effectiveness of the CG model through the test simulation, where the calculated particle size was the same as the CG particle without using the CG model. Chu et al. developed a comprehensive CG CFD-DEM to simulate multiphase flow in a cyclone, and they analyzed the error caused by this model (Chu et al., 2016). In their model, original particles with the same properties, e.g., size and density, can be grouped into a CG particle. They emphasized that particles included in the CG particle had equal translational and rotational movement. The conclusion demonstrated that the CG CFD-DEM model was more accurate in the scenario where the flow was dilute. Lu et al. extended the CG method to the simulation of liquid–solid reacting flows (Lu et al., 2016). They determined the diameter of the CG particle with the energy minimization multiscale (EMMS) model and calculated the collision force using the collision diameter. Based on CG CFD-DEM, Lan et al., 2020, Lu et al., 2020 studied the residence time distribution of particles. Hu et al. (2019) investigated the effect of different operating parameters on the coal gasification process with CG model. About the coarse ratio, Mu et al. (2020) emphasized that it possessed an upper limit which was related to the length scales at the mesoscale that can still be resolved.
Wear occurs when a surface is subjected to the intensive impingement of particles for being exposed to the dense gas–solid flow (Finnie, 1960). This will notably shorten the service life of the device and affect the continuous operation and the economics of fluidized bed plants. Therefore, researches on the wear mechanism should be paid attention to. Investigating erosion properties with numerical methods has become a focus in this field, and many erosion models have been proposed such as the kinetic theory erosion model and kinetic energy dissipation erosion model based on the Eulerian-Eulerian method. These models have low demand for computing resources as opposed to the Eulerian-Lagrangian models whose detailed information of particles is obtained by the Lagrangian method. However, due to the latent defect of those methods that only mean values on a small computational volume can be given, developing a new erosion model to accurately predict the wear of the immersed tubes is an urgent need. Zhao et al. (2017) first proposed a quantitative erosion model called shear impact energy model (SIEM) based on Ashrafizadeh’s work (Ashrafizadeh and Ashrafizadeh, 2012) and proved its accuracy through comparing the simulation results with previous experimental data. In the current work, SIEM is also adopted to study the erosion of the immersed tube in the fluidized bed.
Up to now, the CG model is seldom applied to the fluidized bed with an immersed tube, especially the wear occurring on it. In this study, a CG CFD-DEM is proposed and verified to be applicable and accurate to investigate the fluidized bed with an immersed tube through comparing the macroscopic results obtained from CG CFD-DEM and that predicted by original CFD-DEM as well as experiment results, i.e., the flow pattern, the solid velocity, the pressure drop, the bed height, and the void fraction around the tube. In the CG model developed by us, the energy dissipation rate of the CG system agrees with that of the original system through modifying key collision parameters, referred to as the coefficient of the sliding friction. For the first time, the CG CFD-DEM combining SIEM is applied to the prediction of erosion around the tube. Besides, the simulation time consumed by CG CFD-DEM and original CFD-DEM for a given time size is recorded to prove the capacity of reducing computing cost for the CG model.
Section snippets
Model of discrete phase
According to Newton’s laws of motion, the translational motion and rotational motion of a particle can be described by:where m, v, I, ω, Fc, Fd, Fb, g, Tc indicate solid mass, translational velocity, inertia tensor, rotational velocity, contact force, drag force, buoyancy, gravitational acceleration, contact torque, respectively. The contact force Fc includes the normal and tangential component, which is given by:where Fc,n is normal contact force, Fc,t is
Experiment and simulation conditions
The set-up used in the experiments and simulations is the same and schematically shown in Fig. 1(a), which has a section of 200 mm × 30 mm (width × depth) and a height of 1000 mm, and a circular tube with a diameter of 60 mm is inserted into the bed 120 mm above the gas inlet. The properties of the material are given in Table 1. A total of 310,000 glass beads with an average diameter of 1.5 mm and a density of 2576 kg/m3 are used. In the experiment, the material of the wall is transparent PMMA
Particle flow patterns
Two typical flow patterns are distinguished according to the minimum fluidization velocity (Ku et al., 2013): the packed bed and the fluidized bed. When the superficial velocity is less than the minimum fluidization velocity, particles heap at the bottom of the bed and the packed bed is formed. The initial packed states of experiments and simulations with different coarse ratio are presented in Fig. 2. As illustrated in this figure, the initial bed heights agree well between the experiments and
Conclusions
In this study, a CFD-DEM method combining the CG model is developed to investigate the fluidized bed with an immersed tube and the following conclusions are obtained:
- (1)
The applicability and the adequacy of the new-developed CG CFD-DEM are validated by comparison with the original CFD-DEM and the experiments qualitatively and quantitatively. The flow patterns of original particles at different superficial velocity can be reasonably simulated by CG particles but sacrificing a little micro-scale
CRediT authorship contribution statement
Lianyong Zhou: Validation, Data curation, Visualization, Writing - original draft, Writing - review & editing. Yongzhi Zhao: Conceptualization, Supervision, Methodology, Software.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Acknowledgement
This research is financially supported by the National Key Research and Development Program of China (2019YFC1805605).
References (57)
- et al.
Computational fluid dynamics of high density circulating fluidized bed riser: Study of modeling parameters
Powder Technol.
(2008) - et al.
Investigation of bubble behavior in fluidized beds with and without immersed horizontal tubes using a digital image analysis technique
Powder Technol.
(2011) - et al.
Numerical simulation and experimental validation of bubble behavior in 2D gas–solid fluidized beds with immersed horizontal tubes
Chem. Eng. Sci.
(2011) - et al.
A numerical 3D simulation for prediction of wear caused by solid particle impact
Wear.
(2012) - et al.
An experimentally validated coarse-grain DEM study of monodisperse granular mixing
Powder Technol.
(2020) - et al.
Applicability of a coarse-grained CFD–DEM model on dense medium cyclone
Miner. Eng.
(2016) The voidage function for fluid-particle interaction systems
Int. J. Multiphase Flow.
(1994)Erosion of surfaces by solid particles
Wear.
(1960)- et al.
Particle scale study of heat transfer in packed and fluidized beds of ellipsoidal particles
Chem. Eng. Sci.
(2016) - et al.
Flow patterns near a solid obstacle in a fluidized bed
Chem. Eng. Sci.
(1964)