A new drag model for TFM simulation of gas–solid bubbling fluidized beds with Geldart-B particles

In this work, a new drag model for TFM simulation in gas–solid bubbling fluidized beds was proposed, and a set of equations was derived to determine the meso-scale structural parameters to calculate the drag characteristics of Geldart-B particles under low gas velocities. In the new model, the meso-scale structure was characterized while accounting for the bubble and meso-scale structure effects on the drag coefficient. The Fluent software, incorporating the new drag model, was used to simulate the fluidization behavior. Experiments were performed in a Plexiglas cylindrical fluidized bed consisting of quartz sand as the solid phase and ambient air as the gas phase. Comparisons based on the solids hold-up inside the fluidized bed at different superficial gas velocities, were made between the 2D Cartesian simulations, and the experimental data, showing that the results of the new drag model reached much better agreement with experimental data than those of the Gidaspow drag model did.     

鼓泡流化床中流动特性的多尺度数值模拟

鼓泡流化床因其较高的传热特性以及较好的相间接触已经被广泛应用于工业生产中,而对鼓泡流态化气固流动特性的充分认知是鼓泡流化床设计的关键.在鼓泡流化床中,气泡相和乳化相的同时存在使得床中呈现非均匀流动结构,而这种非均匀结构给鼓泡流化床的数值模拟造成了很大的误差.基于此,以气泡作为介尺度结构,建立了多尺度曳力消耗能量最小的稳定性条件,构建了适用于鼓泡流化床的多尺度气固相间曳力模型.结合双流体模型,对A类和B类颗粒的鼓泡流化床中气固流动特性进行了模拟研究,分析了气泡速度、气泡直径等参数的变化规律.研究表明,与传统的曳力模型相比,考虑气泡影响的多尺度气固相间曳力模型给出的曳力系数与颗粒浓度的关系是一条分布带,建立了控制体内曳力系数与局部结构参数之间的关系.通过模拟得到的颗粒浓度和速度与实验的比较可以发现,考虑气泡影响的多尺度曳力模型可以更好地再现实验结果.通过A类和B类颗粒的鼓泡床模拟研究发现,A类颗粒的鼓泡床模拟受多尺度曳力模型的影响更为显著.      

Two-fluid modeling of Geldart A particles in gas-fluidized beds

We have investigated the effect of cohesion and drag models on the bed hydrodynamics of Geldart A particles based on the two-fluid (TF) model.For a high gas velocity Uo=0.03m/s, we found a transition from the homogeneous fluidization to bubbling fluidization with an increase of the coefficient C1, which is used to account for the contribution of cohesion to the excess compressibility. Thus cohesion can play a role in the bed expansion of Geldart A particles. Apart from cohesion, we have also investigated the influence of the drag models. When using the Wen and Yu drag correlation with an exponent n=4.65, we find an under-prediction of the bed expansion at low gas velocities (Uo=0.009 m/s). When using a larger exponent (n=9.6), as reported in experimental studies of gas-fluidization,a much better agreement with the experimental bed expansion is obtained. These findings suggest that at low gas velocity,a scale-down of the commonly used drag model is required. On the other hand, a scale-up of the commonly used drag model is necessary at high gas velocity (Uo=0.2 and 0.06 m/s). We therefore conclude that scaling the drag force represent only an ad hoc way of repairing the deficiencies of the TF model, and that a far more detailed study is required into the origin of the failure of the TF model for simulating fluidized beds of fine powders.     

Development of a filtered drag model considering effect of the solid shear rate

This study investigates a 2D gas–solid fast fluidized bed of typical Geldart A particles using highly resolved simulations with two-fluid model. The results show that the solid shear rate has a considerable impact on the orientation of the meso-scale structures and hence on the filtered drag force. On the basis of the correlation for the filtered drag force established in the literature using the traditional markers (such as filtered solid volume fraction, filtered slip velocity and filter scale), a correction correlated with the solid shear rate in the direction of gravity is proposed for better prediction of the filtered drag force. The corrected model is shown to produce improved results in posterior tests of flows in different fluidization regimes including bubbling, turbulent and fast fluidization.     

Insights in hydrodynamics of bubbling fluidized beds at elevated pressure by DEM–CFD approach

A numerical simulation was conducted to study the effect of pressure on bubble dynamics in a gas–solid fluidized bed. The gas flow was modeled using the continuum theory and the solid phase, by the discrete element method (DEM). To validate the simulation results, calculated local pressure fluctuations were compared with corresponding experimental data of 1-mm polyethylene particles. It was shown that the model successfully predicts the hydrodynamic features of the fluidized bed as observed in the experiments. Influence of pressure on bubble rise characteristics such as bubble rise path, bubble stability, average bubbles diameter and bubble velocity through the bed was investigated. The simulation results are in conformity with current hydrodynamic theories and concepts for fluidized beds at high pressures. The results show further that elevated pressure reduces bubble growth, velocity and stability and enhances bubble gyration through the bed, leading to change in bed flow structure.     

Counter-current three-phase fluidization in a turbulent contact absorber: A CFD simulation

A computational fluid dynamics study of three-phase counter-current fluidization occurring in a turbulent contact absorber was performed. A two-dimensional, transient Eulerian multi-fluid model was used, in which the dispersed solid phase was modeled employing a kinetic theory of granular flow. The grid independence of the model, the effect of wall boundary conditions, the choice of granular temperature model, the effects of order of discretization scheme and drag models were studied for a base case setting. The results of simulations were validated against experimental results obtained from the literature. Once the model settings were finalized, simulations were performed for different gas and liquid velocities to predict the hydrodynamics of the absorber. Computed bed expansion and pressure drop were compared with experimental data. Good agreement between the two was found for low velocities of gas and liquid.     

Two-fluid modeling of Geldart A particles in gas–solid micro-fluidized beds

The fluidization behavior of Geldart A particles in a gas–solid micro-fluidized bed was investigated by Eulerian–Eulerian numerical simulation. The commonly used Gidaspow drag model was tested first. The simulation showed that the predicted minimum bubbling velocities were significantly lower than the experimental data even when an extremely fine grid size (of approximately one particle diameter) was used. The modified Gibilaro drag model was therefore tested next. The predicted minimum bubbling velocity and bed voidage were in reasonable agreement with the experimental data available in literature. The experimentally observed regime transition phenomena from bubbling to slugging were also reproduced successfully in the simulations. Parametric studies indicated that the solid-wall boundary conditions had a significant impact on the predicted gas and solid flow behavior.     

STUDY ON THE OVERALL PRESSURE BALANCE OF A DOWNFLOW CIRCULATING FLUIDIZED BED SYSTEM

A pressure balance model for a circulating fluidized bed unit that incorporates a downer has been proposed. The model predictions were validated with the experimental data obtained from a special cold-model circulating fluidized bed. Comparison of the operation stability between a CFB downer and a CFB riser has been carried out. Only one critical gas velocity exists in the CFB-riser for a given riser solids flux, while there can be many critical gas velocities for the operation of a CFB downer. Therefore, it is possible to achieve high solids concentration in a CFB downer if appropriate operating conditions are used.     

STUDY ON THE OVERALL PRESSURE BALANCE OF A DOWNFLOW CIRCULATING FLUIDIZED BED SYSTEM

A pressure balance model for a circulating fluidized bed unit that incorporates a downer has been proposed. The model predictions were validated with the experimental data obtained from a special cold-model circulating fluidized bed. Comparison of the operation stability between a CFB downer and a CFB riser has been carried out. Only one critical gas velocity exists in the CFB-riser for a given riser solids flux, while there can be many critical gas velocities for the operation of a CFB downer. Therefore, it is possible to achieve high solids concentration in a CFB downer if appropriate operating conditions are used.     

Discrete particle modeling of granular temperature distribution in a bubbling fluidized bed

The discrete hard sphere particle model (DPM) is applied in this work to study numerically the distributions of particle and bubble granular temperatures in a bubbling fluidized bed. The dimensions of the bed and other parameters are set to correspond to those of Müller et al. (2008). Various drag models and operational parameters are investigated to find their influence on particle and bubble granular temperatures. Various inlet superficial gas velocities are used in this work to obtain their effect on flow characteristics. It is found that the superficial gas velocity has the most important effect on granular temperatures including bubble granular temperature, particle translational granular temperature and particle rotational granular temperature. The drag force model affects more seriously the large scale variables such as the bubble granular temperature. Restitution coefficient influences all granular temperatures to some degree. Simulation results are compared with experimental results by Müller et al. (2008) showing reasonable agreement.     

Two-fluid modeling of Geldart A particles in gas-solid micro-fluidized beds

The fluidizadon behavior of Geldart A particles in a gas-solid micro-fluidized bed was investigated by Eulerian-Eulerian numerical simulation.The commonly used Gidaspow drag model was tested first.The simulation showed that the predicted minimum bubbling velocities were significantly lower than the experimental data even when an extremely fine grid size(of approximately one particle diameter) was used.The modified Gibilaro drag model was therefore tested next.The predicted minimum bubbling velocity and bed voidage were in reasonable agreement with the experimental data available in literature.The experimentally observed regime transition phenomena from bubbling to slugging were also reproduced successfully in the simulations.Parametric studies indicated that the solid-wall boundary conditions had a significant impact on the predicted gas and solid flow behavior.     

Numerical investigations on gas–solid flow in circulating fluidized bed risers using a new cluster-based drag model

A cluster-based drag model is proposed for the gas–solid circulating fluidized bed (CFB) riser by including the cluster information collected from image processing and wavelet analysis into the calculation of system drag. The performance of the proposed drag model is compared with some commonly used drag models. A good agreement with the experimental data is achieved by the proposed cluster-based drag model. Error analysis of the proposed cluster-based drag model based on the local distributions of solids holdup and particle velocity is conducted. The clustering phenomenon in the low-density and high-density CFB risers and the effect of the cluster size on the simulation accuracy are also numerically studied by the proposed drag model.     

Important relationship between meso-scale structure and transfer coefficients in fluidized beds

Along with the fast development of computer technology and measurement techniques, fundamental research on fluidization is faced with both new challenges and opportunities. Among others, great attention should be focused on the meso-scale structure of fluidized beds, to study the quantitative prediction theory and optimum control method for the meso-scale structure of fluidized beds, and to establish the modeling of the relationship between meso-scale structure and momentum transfer, heat transfer, mass transfer, and chemical reaction. These efforts, combined with advanced computer simulation, are expected to solve difficult problems in optimum control and scale-up of fluidized bed processes and equipment.     

Statistical and frequency analysis of the pressure fluctuation in a fluidized bed of non-spherical particles

In this paper, the pressure fluctuation in a fluidized bed was measured and processed via standard devia- tion and power spectrum analysis to investigate the dynamic behavior of the transition from the bubbling to turbulent regime. Two types （Geldart B and D） of non-spherical particles, screened from real bed materials, and their mixture were used as the bed materials. The experiments were conducted in a semi- industrial testing apparatus. The experimental results indicated that the fluidization characteristics of the non-spherical Geldart D particles differed from that of the spherical particles at gas velocities beyond the transition velocity Uo The standard deviation of the pressure fluctuation measured in the bed increased with the gas velocity, while that measured in the plenum remained constant. Compared to the coarse particles, the fine particles exerted a stronger influence on the dynamic behavior of the fluidized bed and promoted the fluidization regime transition from bubbling toward turbulent. The power spectrum of the pressure fluctuation was calculated using the auto-regressive （AR） model; the hydrodynamics of the flu- idized bed were characterized by the major frequency of the power spectrum of the pressure fluctuation. By combining the standard deviation analysis, a new method was proposed to determine the transition velocity Uk via the analysis of the change in the major frequency. The first major frequency was observed to vary within the range of 1.5 to 3 Hz.     

Drag of a blocked fluidized bed

The characteristics of a blocked fluidized bed are presented. The drag of the first and second layers of disks forming the solid phase of the bed and the static pressure at any depth of the structurally inhomogeneous bed as a function of the characteristic parameters were obtained from experimental data. The relation between the density of arrangement of the drag elements in the bed and the stability of its specific drag depthwise is examined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1977.In conclusion, the author thanks A. A. Barmin for editing the articel and valuable methodological advice.     

DEM simulation of polydisperse systems of particles in a fluidized bed

Numerical simulations based on three-dimensional discrete element model (DEM) are conducted for mono-disperse, binary and ternary systems of particles in a fluidized bed. Fluid drag force acting on each particle depending on its size and relative velocity is assigned. The drag coefficient corresponding to Ergun’s correlation is applied to the system of fluidized bed with particle size ratios of 1:1 for the mono-disperse system, 1:1.2, 1:1.4 and 1:2 for the binary system and 1:1.33:2 for the ternary system b...     

Effects of initial static bed height on fractional conversion and bed pressure drop in tapered-in and tapered-out fluidized bed reactors

In this article, a standard 2D Two-Fluid Model (TFM) closed by the kinetic theory of granular flow (KTGF) has been applied to simulate the behavior of tapered-in and tapered-out fluidized bed reactors. In this regard, two types of chemical reactions with gas volume reduction and increase were considered to investigate the effects of initial static bed height on the fractional conversion and bed pressure drop. To validate the CFD model predictions, the results of hydrodynamic simulations concerning bed pressure drop and bed expansion ratio were compared against experimental data reported in the literature and excellent agreement was observed. The obtained simulation results clearly indicate that there is an appropriate static bed height in a tapered-in reactor in which the fractional conversion becomes maximum at this height; whereas variations of static bed height in a tapered-out reactor have insignificant influences on the fractional conversion. Moreover, it was found that the residence time, temperature, and intensity of turbulence of the gas phase are three important factors affecting the fractional conversion in tapered fluidized bed reactors. In addition, it was observed that increasing the static bed height increases the bed pressure drop for both the tapered-in and tapered-out fluidized bed reactors.     

Euler–Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method (DEM) using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model. The model takes into account detailed gas phase chemistry, zero-dimensional modeling of the pyrolysis and gasification of each individual particle, particle shrinkage, and heat and mass transfer between the gas phase and the particulate phase. We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.     

Particle-scale and sub-grid drag models coupled CFD for simulating the CO methanation in a CFB riser

The diffusion and chemical reactions inside the catalyst particles and the heterogeneous flow structure in the computational cells are key factors to affect the accuracy of the coarse-grid simulation in circulating fluidized bed (CFB) methanation reactors. In this work, a particle-scale model is developed to calculate the effective reaction rate considering the transient diffusion and chemical reactions in the particle scale, i.e., the scale of the single catalyst particle. A modified sub-grid drag model is proposed to consider the effects of the meso-scale and chemical reactions on the heterogeneous gas-solid interaction, where the meso-scale is between the single particle and the whole reactor and featured with the particle cluster. Subsequently, a coupled model is developed by integrating the particle-scale and modified sub-grid drag models into CFD. Moreover, the coupled model is validated to achieve accurate predictions on the CO methanation process in a CFB riser. Notably, the coupled model can be performed with a coarse grid (∼58 times particle diameter) and a large time step (0.005 s) to accelerate the simulation. By simply changing the reaction kinetics, different gas-solid catalytic reaction systems can be simulated by using the coupled model.