Numerical simulation of gas-particle flow with a second-order moment method in bubbling fluidized beds

Flow behavior of gas and particles is performed by means of gas-solid two-fluid model with the second-order moment model of particles in the bubbling fluidized bed. The distributions of velocity and moments of particles are predicted in the bubbling fluidized beds. Predictions are compared with experimental data measured by Jung et al. (2005) in a bubbling fluidized bed and Patil et al. (2005) experiments in a bubbling fluidized bed with a jet. The simulated second-order moment in the vertical direction is on average 1.5-2.3 times larger than that in the lateral direction in the bubbling fluidized bed (Jung et al., 2005). For a bubbling fluidized bed with a jet, the ratio of normal second-order moment in the vertical direction to in the lateral direction is in the range of 0.5-2.5 (Patil et al., 2005). The bubblelike Reynolds normal stresses per unit bulk density used by Gidaspow et al. (2002) are computed from the simulated hydrodynamic velocities. The simulated bubblelike Reynolds normal stresses in the vertical direction is on average 4.5-6.0 times larger than that in the lateral direction in the bubbling fluidized bed (Jung et al., 2005). The predictions are in agreement with experimental second-order moments measured by Jung et al. (2005) and porosity measured by Patil et al. (2005).     

The motion mechanism and characteristic of bubble in a pseudo-2D tapered fluidized bed

The different carbon nanotube (CNT) particles (^@A and ^@V) were bed materials in the pseudo-2D tapered fluidized bed (TFB) with/without a distributor. A detailed investigation of the motion mechanism of bubbles was carried out. The high-speed photography and image analysis techniques were used to study bubble characteristic and mixing behavior in the tapered angle of TFB without a distributor. The fractal analysis method was used to analyze the degree of particles movement. Results showed that an S-shaped motion trajectory of bubbles was captured in the bed of ^@V particles. The population of observational bubbles in the bed of ^@V particles was more than that of ^@A particles, and the bubble size was smaller in the bed of ^@V particles than that of ^@A particles. The motion mechanism of bubbles had been shown to be related to bed materials and initial bed height in terms of analysis and comparison of bubble diameter, bubble aspect ratio and bubble shape factor. Importantly, compared to the TFB with a distributor, the TFB without a distributor had been proved to be beneficial to the CNT fluidization according to the study of bubble characteristic and the degree of the particle movement. Additionally, it was found that the mixing behavior of ^@V particles was better than ^@A particles in the tapered angle of TFB without a distributor.     

采用改进的曳力模型模拟2D鼓泡流化床的流化特性

针对特定颗粒浓度范围颗粒团聚作用导致的曳力下降问题,基于传统曳力模型在不同颗粒浓度段的特性分析,选择适用于不同颗粒浓度区间的曳力模型,通过光滑函数得到改进的曳力模型,并耦合欧拉双流体模型对2D鼓泡流化床进行数值模拟. 结果表明,与Gidaspow和Syamlal-O'Brien模型相比,改进的曳力模型对床层局部压降的预测结果更好;随表观气速增加,改进的曳力模型能更准确地预测床层膨胀;当表观气速Ug=0.46 m/s时,改进的曳力模型对径向颗粒浓度分布的模拟结果明显好于Syamlal-O'Brien模型.     

Three-dimensional numerical simulation of a horizontal rotating fluidized bed

Three-dimensional numerical simulations of a horizontal rotating fluidized bed (RFB) containing glass bead particles (d_s = 82 μm, ρ_s = 2450 kg/m³) and washed alumina (d_s = 89 μm, ρ_s = 1550 kg/m³) were performed. FLUENT 6.1 software was used to carry out our simulation. The numerical results were compared with the experimental data of Qian and Pfeffer et al. [G.H. Qian, I. Bagyi, I.W. Burdick, R. Pfeffer, H. Shaw, Gas-Solid Fluidization in a Centrifugal Field.” AIChE J. 47 (5) (2001) 1022-1034]. The rotating speed of the RFB was set at 325 rpm (34 rad/s), which is equivalent to a centrifugal acceleration of 7 g.The flow behavior of the solid particles was analyzed; the bed thickness and the calculated pressure drop were compared with the experimental results. Our calculated pressure drop agreed very well with the experimental results.     

拟泡-乳曳力模型在大型高温费托流化床反应器中的CFD模拟研究

以带冷却盘管的大型高温费托流化床反应器为研究对象,开展三维计算流体力学模拟研究。传统双流体模型基于局部平均的假设,认为单位控制体内气固两相均匀分布,网格尺寸必须足够小才能正确揭示局部非均匀结构的所有细节。采用双流体模型模拟大型工业化流化床装置时,将导致网格数量过于庞大,远超现有计算能力。为提高计算效率的同时不损失模拟精度,提出了基于局部非均匀假设、适用于粗网格的拟泡-乳三相非均匀曳力(PBTD)模型。该模型将流化床分为乳化相气体、乳化相颗粒以及气泡三相,分别建立守恒方程,体现气泡的非均匀特性对气固曳力的影响。乳化相内气固曳力以及气泡相与乳化相内颗粒的曳力分开考虑。采用PBTD模型耦合传质和反应模型,建立基于局部非均匀假设的高温费托合成反应器三维流动-传递-反应模型,包括各相守恒控制方程、气泡尺寸模型、相间物质和动量交换模型、高温费托合成反应动力学模型以及初始和边界条件,预测反应器内的流场和组分浓度分布。研究结果表明:在粗网格条件下,非均匀曳力模型可以预测床层内相含率的分布情况,预测的床层膨胀高度与经验公式计算值接近,偏差为1.2%。反应器出口气体组分的质量分数与试验测量值相近,偏差在1.5%~16.0%。模拟结果证实,基于非均匀假设的PBTD模型适用于模拟工业规模的鼓泡流化床反应器,对其设计开发和工业运行具有指导价值。     

Numerical modeling of particle fluidization behavior in a rotating fluidized bed

In this study, numerical modeling of particle fluidization behaviors in a rotating fluidized bed (RFB) was conducted. The proposed numerical model was based on a DEM (Discrete Element Method)-CFD (Computational Fluid Dynamics) coupling model. Fluid motion was calculated two-dimensionally by solving the local averaged basic equations. Particle motion was calculated two-dimensionally by the DEM. Calculation of fluid motion by the CFD and particle motion by the DEM were simultaneously conducted in the present model. Geldart group B particles (diameter and particle density were 0.5 mm and 918 kg/m³, respectively) were used for both calculation and experiment. First of all, visualization of particle fluidization behaviors in a RFB was conducted. The calculated particle fluidization behaviors by our proposed numerical model, such as the formation, growth and eruption of bubble and particle circulation, showed good agreement with the actual fluidization behaviors, which were observed by a high-speed video camera. The estimated results of the minimum fluidization velocity (U_mf) and the bed pressure drop at fluidization condition (ΔP_f) by our proposed model and other available analytical models in literatures were also compared with the experimental results. It was found that our proposed model based on the DEM-CFD coupling model could predict the U_mf and ΔP_f with a high accuracy because our model precisely considered the local downward gravitational effect, while the other analytical models overpredicted the ΔP_f due to ignoring the gravitational effect.     

网格尺寸对循环流化床内气固两相流动的影响

将基于能量最小多尺度方法(EMMS)的曳力模型耦合到双流体模型中,并针对循环流化床内的气固两流动进行了模拟研究。采用全滑移壁面边界条件处理颗粒相,考察了3种网格尺度对轴向空隙率和出口颗粒循环量等气固流动特性的影响。计算结果表明,应用EMMS曳力模型处理相间作用力,同时在采用全滑移壁面边界条件处理颗粒相时,双流体模型能够正确预测轴向空隙率分布。采用网格尺寸为2.325 mm×20 mm时,模拟结果和实测数据吻合较好,表明在循环流化床的数值模拟中选择恰当的网格尺度是极为重要的。     

A CFD‐PBM‐PMLM integrated model for the gas–solid flow fields in fluidized bed polymerization reactors

Although the use of computational fluid dynamics (CFD) model coupled with population balance (CFD‐PBM) is becoming a common approach for simulating gas–solid flows in polydisperse fluidized bed polymerization reactors, a number of issues still remain. One major issue is the absence of modeling the growth of a single polymeric particle. In this work a polymeric multilayer model (PMLM) was applied to describe the growth of a single particle under the intraparticle transfer limitations. The PMLM was solved together with a PBM (i.e. PBM‐PMLM) to predict the dynamic evolution of particle size distribution (PSD). In addition, a CFD model based on the Eulerian‐Eulerian two‐fluid model, coupled with PBM‐PMLM (CFD‐PBM‐PMLM), has been implemented to describe the gas–solid flow field in fluidized bed polymerization reactors. The CFD‐PBM‐PMLM model has been validated by comparing simulation results with some classical experimental data. Five cases including fluid dynamics coupled purely continuous PSD, pure particle growth, pure particle aggregation, pure particle breakage, and flow dynamics coupled with all the above factors were carried out to examine the model. The results showed that the CFD‐PBM‐PMLM model describes well the behavior of the gas–solid flow fields in polydisperse fluidized bed polymerization reactors. The results also showed that the intraparticle mass transfer limitation is an important factor in affecting the reactor flow fields. © 2011 American Institute of Chemical Engineers AIChE J, 58: 1717–1732, 2012     

Numerical simulations and comparative analysis of two- and three-dimensional circulating fluidized bed reactors for CO2 capture

Carbon dioxide (CO₂), the main gas emitted from fossil burning, is the primary contributor to global warming. Circulating fluidized bed reactor (CFBR) is proved as an energy-efficient method for post-combustion CO₂ capture. The numerical simulation by computational fluid dynamics (CFD) is believed as a promising tool to study CO₂ adsorption process in CFBR. Although three-dimensional (3D) simulations were proved to have better predicting performance with the experimental results, two-dimensional (2D) simulations have been widely reported for qualitative and quantitative studies on gas–solid behavior in CFBR for its higher computational efficiency recently. However, the discrepancies between 2D and 3D simulations have rarely been evaluated by detailed study. Considering that the differences between the 2D and 3D simulations will vary substantially with the changes of independent operating conditions, it is beneficial to lower computational costs to clarify the effects of dimensionality on the numerical CO₂ adsorption runs under various operating conditions. In this work, the comparative analysis for CO₂ adsorption in 2D and 3D simulations was conducted to enlighten the effects of dimensionality on the hydrodynamics and reaction behaviors, in which the separation rate, species distribution and hydrodynamic characteristics were comparatively studied for both model frames. With both accuracy and computational costs considered, the viable suggestions were provided in selecting appropriate model frame for the studies on optimization of operating conditions, which directly affect the capture and energy efficiencies of cyclic CO₂ capture process in CFBR.     

A computational fluid dynamics design of a carbon dioxide sorption circulating fluidized bed

A kinetic theory based hydrodynamic model with experimentally determined sorption rates for reaction of CO₂ with K₂CO₃ solid sorbent is used to design a compact circulating fluidized bed sorption‐regeneration system for CO₂ removal from flue gases. Because of high solids fluxes, the sorber does not require internal or external cooling. The output is verified by computing the granular temperatures, particle viscosities, dispersion, and mass transfer coefficients. These properties agree with reported measurement values except the radial dispersion coefficients, which are much higher due to the larger bed diameter. With the solid sorbent prepared according to published information, the CO₂ removal percentage at the riser top is 69.16%. To improve the CO₂ removal, an effort is needed to develop a better sorbent or to simply lower the inlet gas velocity to operate in a denser mode, leading to a larger system. Also, the effect of temperature rise on the removal efficiency is investigated. © 2010 American Institute of Chemical Engineers AIChE J, 2010     

Insights in distributed secondary gas injection in a bubbling fluidized bed via discrete particle simulations

Experiments have shown that distributed secondary gas injection via a fractal injector in fluidized beds can significantly reduce the bubble size, and may also decrease the bubble fraction. In order to gain insight into the distribution of the gas between the phases and the mechanisms behind these effects simulations of small bubbling fluidized beds with one or two secondary gas injection points were carried out using a discrete particle model. Although the systems are very small, so that wall effects cannot be excluded, the model predicts that the bubble size and bubble fraction both decrease with secondary gas injection, while the gas flow through the dense phase increases. The secondary gas tends to stay in the dense phase, which limits the amount of gas available to form bubbles and is the main contributor to the decrease in the bubble size and fraction. The gas-solid contact improves as a result.     

CFD modeling of tapered circulating fluidized bed reactor risers: Hydrodynamic descriptions and chemical reaction responses

A comprehensive two-dimensional transient Eulerian model combined with the kinetic theory of granular flow was developed to obtain the hydrodynamic and chemical reaction behaviors in tapered circulating fluidized bed reactor risers. In this study, the focus was on the chemical reactions and its behaviors inside three different riser geometries. The model was verified by using an experimental dataset from the literature, and was then used for both predicting the hydrodynamic behaviors and computing the system turbulent properties. The tapered-out riser improves the system turbulence or mixing which can be explained by the dispersion coefficients. On the other hand, the tapered-in riser increases the solid particle residence time and gives a more uniform temperature distribution, because it does not have sufficient force to support the weight of the particles. The same riser geometries but with the addition of the chemical reaction were then used for evaluating the previously proposed criteria that the riser geometry should be chosen with respect to the characteristics of the reactions. Reactions with a medium reaction rate were best suited to the typical riser, whilst reactions with a fast and slow reaction rate best fitted the tapered-out and tapered-in risers, respectively.     

Dynamics of gas-particle flow in circulating fluidized beds

The flow of gas—particle mixtures in a circulating fluidized bed has been studied, probing the flow behavior under both stable and unstable operating conditions. A novel feature of our work is the use of electrical capacitance tomography to image particle distribution over the cross-section at one elevation in the standpipe. In addition, we have also obtained data on pressure profile and aeration rate in the standpipe, particle circulation rate in the circulating fluidized bed and riser gas flow rate under various operating conditions. Here, we report experimental results obtained for two different particles, both belonging to Geldart type A. At low aeration rates, a stable dense phase flow is obtained in the standpipe. At high aeration rates, the flow became unstable, manifesting low frequency oscillations in the flow characteristics. Our results suggest that, under conditions explored in the present study, this instability originates in the standpipe and that any attempt to model it should consider the interaction between the various components of the circulating fluidized bed system.     

Computation of system turbulences and dispersion coefficients in circulating fluidized bed downer using CFD simulation

In this study, the Eulerian computational fluid dynamics model with the kinetic theory of granular flow model was effectively used to compute the system turbulences and dispersion coefficients in a circulating fluidized bed (CFB) downer. In addition, the obtained model was used to simulate all the system velocities.     

Particle-scale simulation of fluidized bed with immersed tubes

In order to simulate gas-solids flows with complex geometry, the boundary element method was incorporated into the implementation of a combined model of computational fluid dynamics and discrete element method. The resulting method was employed to simulate hydrodynamics in a fluidized bed with immersed tubes. The transient simulation results showed particle and bubble dynamics. The bubble coalescence and break-up behavior when passing the immersed tubes was successfully predicted. The gassolid flow pattern in the fluidized bed is changed greatly because of the immersed tubes. As particles and gas are come in contact with the immersed tubes, the gas bubbles will be deformed. The collisions between particles and tubes will make the tubes surrounded by air pockets most of the time and this is unfavorable for the heat transfer between particles and tubes. __________ Translated from Chemical Engineering, 2007, 35(11): 21–24 [译自:化学工程]     

沉浸管式流化床的颗粒尺度模拟

为模拟具有复杂几何结构的气固流动系统,文中将计算流体力学和离散单元法与边界元方法结合起来,对沉浸管式流化床内颗粒及气泡的运动行为进行了数值模拟。模拟计算得到的瞬态流型图揭示了气泡绕流沉浸管束时出现的合并和破碎状态及颗粒群的详细运动行为,发现床内气固二相的流动受到沉浸管束存在的显著影响。当颗粒及气相的流动受到沉浸管的阻碍而绕管流动过程中气泡会发生变形,变得扭曲狭长且易被撕碎。同时颗粒与管道壁面碰撞会造成气固二相复杂的动态运动形式,床内的管道大部分时间会被气穴包围,将严重阻碍管道与颗粒之间的传热。     

气固搅拌流化床内压力脉动特性的数值模拟

在传统气固流化床中引入搅拌桨,可减轻聚合物颗粒的黏附并强化流态化过程。采用计算流体力学(CFD)方法对搅拌流化床内的压力脉动特性进行数值模拟,考察流态化过程中的气泡行为。模拟过程采用多重参考坐标系方法解决搅拌桨区域的运动问题,由欧拉双流体模型和颗粒动力学方法模拟气固两相流。床层压力脉动的统计分析和功率谱分析表明,随着搅拌桨转速的增加,流化床内的压力脉动标准偏差和功率谱幅值变小,床层内的平均气泡尺寸减小,床层可由鼓泡流态化向散式流态化转变。     

CFD modeling and validation of the turbulent fluidized bed of FCC particles

An experimental and computational study is presented on the hydrodynamic characteristics of FCC particles in a turbulent fluidized bed. Based on the Eulerian/Eulerian model, a computational fluid dynamics (CFD) model incorporating a modified gas‐solid drag model has been presented, and the model parameters are examined by using a commercial CFD software package (FLUENT 6.2.16). Relative to other drag models, the modified one gives a reasonable hydrodynamic prediction in comparison with experimental data. The hydrodynamics show more sensitive to the coefficient of restitution than to the flow models and kinetics theories. Experimental and numerical results indicate that there exist two different coexisting regions in the turbulent fluidized bed: a bottom dense, bubbling region and a dilute, dispersed flow region. At low‐gas velocity, solid‐volume fractions show high near the wall region, and low in the center of the bed. Increasing gas velocity aggravates the turbulent disorder in the turbulent fluidized bed, resulting in an irregularity of the radial particle concentration profile. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

Effect of bed size on hydrodynamics in 3‐D gas–solid fluidized beds

It is well known that hydrodynamics observed in large scale gas–solid fluidized beds are different from those observed in smaller scale beds. In this article, an efficient two‐fluid model based on kinetic theory of granular flow is applied, with the goal to highlight and investigate hydrodynamics differences between three‐dimensional fluidized beds of diameter 0.10, 0.15, 0.30, 0.60, and 1.0 m, focusing on the bubble and solids flow characteristics in the bubbling regime. Results for the 0.30 m diameter bed are compared with experimental results from the literature. The bubble size evolution closely follows a correlation proposed by Werther for small beds, and a correlation proposed by Darton for sufficiently large beds. The bubble size increases as the bed diameter is increased from 0.10 to 0.30 m, and remains approximately constant for bed diameters from 0.30 to 1.0 m. Concurrently, an increase in bubble rise velocity is observed, with a much high bubble rise velocity in the largest bed of diameter 1.0 m due to gulf stream circulations. The dynamics in shallow and deep beds is predicted to be different, with marked differences in bubble size and solids circulation patterns. © 2015 American Institute of Chemical Engineers AIChE J, 61: 1492–1506, 2015     

催化裂化流化床内构件的研究进展

介绍了近年来催化裂化流化床中内构件的研究现状,主要包括挡板、垂直管束和填料。阐述了气、固两相在装有不同内构件的流化床中的流动特点,指明了各类内构件的优缺点。并详细阐述了挡板式和填料式内构件的改进历程,根据现有内构件存在的优缺点提出了新的内构件的开发方向。