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

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

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 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 exoerimental data than those of the Gidasoow dra~ model did.     

A review of multiscale CFD for gas–solid CFB modeling

Meso-scale structure is of critical importance to circulating fluidized bed (CFB) applications. Computational fluid dynamics (CFD) with consideration of meso-scale structures can help understand the structure-oriented coupling between flow, heat/mass transfer and reactions. This article is to review our recent progress on the so-called multiscale CFD (MSCFD), which characterizes the sub-grid meso-scale structure with stability criteria in addition to conservation equations. It is found that the mesh-independent solution of fine-grid two-fluid model (TFM) without sub-grid structures is inexact, in the sense that it overestimates the drag coefficient and fails to capture the characteristic S-shaped axial profile of voidage in a CFB riser. By comparison, MSCFD approach in terms of EMMS/matrix seems to reach a mesh-independent solution of the sub-grid structure, and succeeds in predicting the axial profile and flow regime transitions. Further application of MSCFD finds that neglect of geometric factors is one of the major reasons that cause disputes in understanding the flow regime transitions in a CFB. The operating diagram should, accordingly, include geometric factors besides commonly believed operating parameters for the intrinsic flow regime diagram. Recent extension of MSCFD to mass transfer finds that Reynolds number is insufficient for correlating the overall Sherwood number in a CFB. This is believed the main reason why the conventional correlations of Sherwood number scatter by several orders of magnitude. Certain jump change of state of motion around Reynolds number of 50–100 can be expected to clarify the abrupt decay of Sherwood number in both classical- and circulating-fluidized beds. Finally, we expect that the real-size, 3-D, full-loop, time-dependent multiscale simulation of CFB is an emerging paradigm that will realize virtual experiment of CFBs.     

Effect of acoustic field on heat transfer in a sound assisted fluidized bed of fine powders

The fluidized beds are widely used in a variety of industries where heat transfer properties of the fluidized system become important for successful operation. Fluidized are preferred in heat recovery processes because of their unique ability of rapid heat transfer and uniform temperature. Fine powders handling and processing technologies have received widespread attention due to increased use of fine powders in the manufacture of drugs, cosmetics, plastics, catalysts, energetics and other advanced materials. A better understanding of fluidization behavior of fine powders is of great importance in applications involving heat transfer, mass transfer, mixing, transporting and modifying surface properties etc. The difficulty in putting the fine powders in suspension with the fluidizing gas is related to the cohesive structure and to the physical forces between the primary particles. The sound waves agitate bubbling and this results in improving solids mixing in the fluidized bed. The improved solids mixing results in uniform and smooth fluidization, which leads to better heat transfer rates in the fluidized bed.     

Heat transfer characterization and improvement in an external catalyst cooler fluidized bed

Gas-solid fluidized beds have been widely used in heat transfer processes,and so there have been many studies focused on increasing heat transfer in such units....     

RADIAL PROFILE OF THE SOLID FRACTION IN A LIQUID—SOLID CIRCULATING FLUIDIZED BED

In a liquid-solid circulating fluidized bed, lateral forces acting on the particles determine the movement of the particles in the radial direction, and this creates a radial profile for the solid fraction. This work proposes a model to calculate the radial profile of the solid fraction in a liquid-solid circulating fluidized bed based on the balance of the lateral force and the turbulent dispersion force.     

Discrete particle simulations of bubble-to-emulsion phase mass transfer in single-bubble fluidized beds

A classical Euler–Lagrangian model for gas–solid flows was extended with gas component mass conservation equations and used to obtain fundamental insights into bubble-to-emulsion phase mass transfer in bubbling gas–solid fluidized beds. Simulations of injected single rising bubbles under incipient fluidization conditions were carried out, using Geldart-A and -B particles. Phenomena observed in the simulations and those of various theoretical models used to derive phenomenological models were compared to challenge the assumptions underlying the phenomenological models. The bubble-to-emulsion phase mass transfer coefficients calculated for the simulations using Geldart-B particles were in a good agreement with predictions made using the Davidson and Harrison (1963) model. The bubble-to-emulsion phase mass transfer coefficients for Geldart-A particles were, however, much smaller than the predictions obtained from theoretical models (e.g. Chiba and Kobayashi (1970)). The newly developed model allows a detailed analysis of various hydrodynamic aspects and their effects on the mass transfer characteristics in and around rising bubbles in fluidized beds.     

Heat transfer and hydrodynamic studies on two- and three-phase fluidized beds of floating bubble breakers

Heat transfer characteristics in three-phase fluidized beds of floating bubble breakers have been studied in a 0.142 m I.D. x 2.0 m high Plexiglas column fitted with an axially mounted cylindrical heater.Effects of the liquid and gas velocities, the particle size, the volume ratio of floating bubble breaker to particles on phase holdup, the vertical bubble length, and the heat transfer coefficient have been determined.In the bubble-coalescing regime, the heat transfer coefficient in three-phase fluidized beds having the volume ratio V_f/V_s of 10–15% produced a maximum increase in heat transfer coefficient of about 20% in comparison to that in the bed without floating bubble breakers. Also, bubble length and gas-phase holdups exhibited their maximum and minimum values at a volume ratio of 10–15%. The heat transfer coefficient in three-phase fluidized beds of floating bubble breakers can be estimated from the surface renewal model with isotropic turbulence theory.Heat transfer coefficients expressed in terms of the Nusselt number have been correlated with the particle Reynolds number and the volume ratio of floating bubble breakers to particles.     

气固两相流介尺度LBM-DEM模型

LBM-DEM耦合方法通常是指一种颗粒流体系统直接数值模拟算法,即是一种不引入经验曳力模型的计算方法,颗粒尺寸通常比计算网格的长度大一个量级,颗粒的受力通过表面的粘性力与压力积分获得,其优点是能描述每个颗粒周围的详细流场,产生详细的颗粒-流体相互作用的动力学信息,可以探索颗粒流体界面的流动、传递和反应的详细信息及两相相互作用的本构关系,但其缺点是计算量巨大,无法应用于真实流化床过程模拟。本文针对气固流化床中的流体以及固体颗粒间的多相流体力学行为,建立了一种稠密气固两相流的介尺度LBMDEM模型,即LBM-DEM耦合的离散颗粒模型,实现在颗粒尺度上流化床的快速离散模拟。该耦合模型采用格子玻尔兹曼方法(LBM)描述气相的流动和传递行为,离散单元法(DEM)用于描述颗粒相的运动,并利用能量最小多尺度(EMMS)曳力解决气固耦合不成熟问题,以提高其模拟精度。通过经典快速流态化的模拟,验证了介尺度LBM-DEM耦合模型的有效性。模拟结果表明介尺度LBM-DEM模型是一种探索实验室规模气固系统的有力手段。     

Micro/meso simulation of a fluidized bed in a homogeneous bubbling regime

Due to the wide range of spatial scales and the complex features associated to fluid/solid and solid/solid interactions in a dense fluidized bed, the system can be studied at different length scales, namely micro, meso and macro. In this work, we select a flow configuration relevant of a homogeneous liquid/solid fluidization and compare computed results from Particle Resolved Simulation (PRS) with those from locally averaged Euler/Lagrange simulation. PRS at the micro-scale is carried out by a parallel Distributed Lagrange Multiplier (DLM) solver in the framework of fictitious domain methods (Wachs, 2011a, 2015). For meso-scale simulations, the set of mass and momentum conservation equations is averaged in control volumes encompassing few particles and momentum transfer between the two phases is modeled using appropriate drag laws. Both methods are coupled to a Discrete Element Method (DEM) combined with a soft-sphere contact model to solve the Newton–Euler equations with collisions for the particles in a Lagrangian framework (Wachs et al., 2012). A test case of intermediate size with 2000 spheres is chosen as a sensible compromise between size limitations of the meso-scale model for an appropriate averaging process and computational resources required to run micro-scale simulations. These two datasets yield new insight on momentum transfer at different spatial scales in the flow, and question the validity of certain approximations adopted in the meso-scale model. Results demonstrate an acceptable agreement between the micro- and meso-scale predictions on integral measures as pressure drop and bed height. Investigating more detailed features of the flow, it has been shown that particles fluctuations are considerably suppressed in meso-scale simulations and in particular the particles transverse motion is underestimated, regardless of the selected drag law. The origin of these dependencies is carefully investigated by reconstructing the closure laws based on PRS results and comparing them to the closure laws proposed in the literature.     

Effects of magnetic fields on improving mass transfer in flue gas desulfurization using a fluidized bed

The effects of magnetic fields on improving the mass transfer in flue gas desulfurization using a fluidized bed are investigated in the paper. In this research, the magnetically fluidized bed (MFB) is used as the reactor in which ferromagnetic particles are fluidized with simulated flue gas under the influence of an external magnetic field. Lime slurry is continuously sprayed into the reactor. As a consequence, the desulfurization reaction and the slurry drying process take place simultaneously in the MFB. In this paper, the effects of ferromagnetic particles and external magnetic fields on the desulphurization efficiency are studied and compared with that of quartz particles as the fluidized particles. Experimental results show that the ferromagnetic particles not only act as a platform for lime slurry to precipitate on like quartz particles, but also take part in the desulfurization reaction. The results also show that the specific surface area of ferromagnetic particles after reaction is enlarged as the magnetic intensity increases, and the external magnetic field promotes the oxidation of S(IV), improving the mass transfer between sulphur and its sorbent. Hence, the efficiency of desulphurization under the effects of external magnetic fields is higher than that in general fluidized beds.     

The local structure of gas fluidized beds —II. The spatial distribution of bubbles

The spatial distribution of bubbles in gas fluidized beds has been investigated with the measuring system described in Part I of this paper in beds of 0.10, 0.20, 0.45 and 1.0 m dia. The results indicate that in gas fluidized beds a characteristic flow profile of the bubble phase exists such that near the distributor a zone of increased bubble development exists in an annulus near the wall. This zone moves towards the vessel centre-line with increasing height above the distributor. The merging of the annular zone in the vessel centre-line marks the beginning of the transition of the fluidized bed to the state of slugging. The spatial flow profile of the bubble phase is shown to be responsible for the existence of characteristic solids circulation patterns in gas fluidized beds.     

Multi-scale analysis of flow structures in fluidized beds with immersed tubes

Multi-scale analysis and non-linear analysis were combined to investigate the hydrodynamics of fluidized beds with and without horizontal tubes.Pressure fluctuations were measured and analyzed employing discrete wavelet analysis,recurrence plot analysis,and recurrence quantification analysis.A systematic procedure was followed to determine wavelet parameters.At low gas velocities,the energy of macrostructures reduces with the addition of the first tube and then increases with the addition of a second tube.However,there is no notable difference at high gas velocities.Determinism is high for the bed without tubes,which can be attributed to the periodic behavior of bubbles.Determinism decreases with the addition of tubes because the breakage of bubbles results in less periodic behavior.The three methods of analysis used in this study captured the effects of immersed tubes on the hydrodynamics of fluidized beds.Recurrence quantitative analysis was found to be a powerful and easy-to-use method that can capture the nonlinear characteristics of fluidized beds much more quickly than conventional methods of nonlinear analysis.This method can thus be effectively used for the online monitoring of hydrodynamic changes in fluidized beds.     

Multi-scale analysis of flow structures in fluidized beds with immersed tubes

Multi-scale analysis and non-linear analysis were combined to investigate the hydrodynamics of fluidized beds with and without horizontal tubes. Pressure fluctuations were measured and analyzed employing discrete wavelet analysis, recurrence plot analysis, and recurrence quantification analysis. A systematic procedure was followed to determine wavelet parameters. At low gas velocities, the energy of macro-structures reduces with the addition of the first tube and then increases with the addition of a second tube. However, there is no notable difference at high gas velocities. Determinism is high for the bed without tubes, which can be attributed to the periodic behavior of bubbles. Determinism decreases with the addition of tubes because the breakage of bubbles results in less periodic behavior. The three methods of analysis used in this study captured the effects of immersed tubes on the hydrodynamics of fluidized beds. Recurrence quantitative analysis was found to be a powerful and easy-to-use method that can capture the nonlinear characteristics of fluidized beds much more quickly than conventional methods of nonlinear analysis. This method can thus be effectively used for the online monitoring of hydrodynamic changes in fluidized beds.     

CFD–DEM study of effect of bed thickness for bubbling fluidized beds

The effect of bed thickness in rectangular fluidized beds is investigated through the CFD–DEM simulations of small-scale systems. Numerical results are compared for bubbling fluidized beds of various bed thicknesses with respect to particle packing, bed expansion, bubble behavior, solids velocities, and particle kinetic energy. Good two-dimensional (2D) flow behavior is observed in the bed having a thickness of up to 20 particle diameters. However, a strong three-dimensional (3D) flow behavior is observed in beds with a thickness of 40 particle diameters, indicating the transition from 2D flow to 3D flow within the range of 20–40 particle diameters. Comparison of velocity profiles near the walls and at the center of the bed shows significant impact of the front and back walls on the flow hydrodynamics of pseudo-2D fluidized beds. Hence, for quantitative comparison with experiments in pseudo-2D columns, the effect of walls has to be accounted for in numerical simulations.     

Phase holdups in three-phase fluidized beds in the presence of disc promoter

Three-phase fluidized beds are found to have wide applications in process industries. The present investigation essentially comprises of the studies on gas holdup, liquid holdup and bed porosity in three-phase fluidized beds with coaxially placed disc promoter. Holdup data were obtained from bed expansion and pressure drop measurements. Analysis of the data was done to elucidate the effects of dynamic and geometric parameters on gas holdup, liquid holdup and bed porosity. Data were correlated and useful equations were obtained from empirical modeling.