Experimental and computational study of the bed dynamics of semi‐cylindrical gas–solid fluidized bed

With computational fluid dynamics (CFD) it is possible to get a detailed view of the flow behaviour of the fluidized beds. A profound and fundamental understanding of bed dynamics such as bed pressure drop, bed expansion ratio, bed fluctuation ratio, and minimum fluidization velocity of homogeneous binary mixtures has been made in a semi‐cylindrical fluidized column for gas–solid systems, resulting in a predictive model for fluidized beds. In the present work attempt has been made to study the effect of different system parameters (viz., size and density of the bed materials and initial static bed height) on the bed dynamics. The correlations for the bed expansion and bed fluctuations have been developed on the basis of dimensional analysis using these system parameters. Computational study has also been carried out using a commercial CFD package Fluent (Fluent, Inc.). A multifluid Eulerian model incorporating the kinetic theory for solid particles was applied in order to simulate the gas–solid flow. CFD simulated bed pressure drop has been compared with the experimental bed pressure drops under different conditions for which the results show good agreements.     

Particle‐scale simulation of the flow and heat transfer behaviors in fluidized bed with immersed tube

A kind of new modified computational fluid dynamics‐discrete element method (CFD‐DEM) method was founded by combining CFD based on unstructured mesh and DEM. The turbulent dense gas–solid two phase flow and the heat transfer in the equipment with complex geometry can be simulated by the programs based on the new method when the k‐ε turbulence model and the multiway coupling heat transfer model among particles, walls and gas were employed. The new CFD‐DEM coupling method that combining k‐ε turbulence model and heat transfer model, was employed to simulate the flow and the heat transfer behaviors in the fluidized bed with an immersed tube. The microscale mechanism of heat transfer in the fluidized bed was explored by the simulation results and the critical factors that influence the heat transfer between the tube and the bed were discussed. The profiles of average solids fraction and heat transfer coefficient between gas‐tube and particle‐tube around the tube were obtained and the influences of fluidization parameters such as gas velocity and particle diameter on the transfer coefficient were explored by simulations. The computational results agree well with the experiment, which shows that the new CFD‐DEM method is feasible and accurate for the simulation of complex gas–solid flow with heat transfer. And this will improve the farther simulation study of the gas–solid two phase flow with chemical reactions in the fluidized bed. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

CFD simulation study of the effect of baffles on the fluidized bed for hydrogenation of silicon tetrachloride

In this study, the effect of channel baffles and louver baffles on the flow pattern in the large-scale industrial fluidized beds was studied by computational fluid dynamics (CFD) methods. Then, the effect of flow pattern on the chemical reaction performance was studied for the first time. Simulation results showed that the gas velocity distributed more uniformly, solid particles dispersed more homogeneously and aggregation scarcely occurred in the fluidized bed with louver baffles than that with channel baffles. The residence time distribution indicated that louver baffles remarkably suppressed gas back-mixing in comparison with channel baffles. The reasonable agreements of pressure distribution and reaction results between the simulation in the bed with channel baffles and the data on a large-scale industrial apparatus demonstrated the accuracy of the CFD model. The predicted conversion of SiCl₄ in the bed with louver baffles (27.44%) was higher than that with channel baffles (22.69%), indicating that louver baffles markedly improved the performance of the fluidized bed. This study could provide useful information for future structural improvements of baffles in large-scale fluidized beds.     

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     

A mechanistic study of agglomeration in fluidised beds at elevated pressures

Effect of operating pressure on the hydrodynamics of agglomerating gas–solid fluidised bed was investigated using a combination of discrete element method (DEM) for describing the movement of particles and computational fluid dynamic (CFD) for describing the flow of the gas phase. The inter‐particle cohesive force was calculated based on a time dependent model developed for solid bridging by the viscous flow. Motion of agglomerates was described by the multi‐sphere method. Fluidisation behaviour of an agglomerating bed was successfully simulated in terms of increasing the size of agglomerates. The results showed that increasing the operating pressure postpones de‐fluidisation of the bed. Since the DEM approach is a particle level simulation and study about particle–particle interactions is possible, a micro‐scale investigation in terms of cohesive force and repulsive force during agglomeration at elevated pressures was done. The micro‐scale results showed that although the number of contacts between particles was decreased by increasing operating pressure, stronger solid bridge formed between colliding particles at higher pressures. © 2012 Canadian Society for Chemical Engineering     

Particle‐resolved direct numerical simulation of gas–solid dynamics in experimental fluidized beds

Particle‐resolved direct numerical simulations (PR‐DNS) of a simplified experimental shallow fluidized bed and a laboratory bubbling fluidized bed are performed by using immersed boundary method coupled with a soft‐sphere model. Detailed information on gas flow and individual particles’ motion are obtained and analyzed to study the gas–solid dynamics. For the shallow bed, the successful predictions of particle coherent oscillation and bed expansion and contraction indicate all scales of motion in the flow are well captured by the PD‐DNS. For the bubbling bed, the PR‐DNS predicted time averaged particle velocities show a better agreement with experimental measurements than those of the computational fluid dynamics coupled with discrete element models (CFD‐DEM), which further validates the predictive capability of the developed PR‐DNS. Analysis of the PR‐DNS drag force shows that the prevailing CFD‐DEM drag correlations underestimate the particle drag force in fluidized beds. The particle mobility effect on drag correlation needs further investigation. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1917–1932, 2016     

Fluid mechanics of liquid–solid fluidization in the centrifugal field

The flow pattern in the bed and in the freeboard of a liquid–solid fluidized bed operated in the centrifugal field was investigated. Centrifugal accelerations between 20 and 500 times the gravity were acting on the particles. Both the radius‐dependent centrifugal acceleration and the Coriolis force were found to strongly influence the fluid dynamics. A CFD simulation of the liquid flow in absence of particles showed high tangential velocities which were two orders of magnitude higher than the radially oriented fluidizing velocity. Measurements with a 1 m diameter rotor confirmed the calculations and revealed that also the liquid–solid fluidized bed is moving in the tangential direction although with lower speed than the pure liquid. The bed expansion under these conditions can be described by a Richardson–Zaki type correlation.     

Measurements and numerical predictions of gas vortices formed by single bubble eruptions in the freeboard of a fluidised bed

Gas vortices generated in the freeboard of a bubbling fluidised bed have become the centre of increasingly more research due to the advances in experimental technology. The behaviour of gas flow in the freeboard of a bubbling fluidised bed is of interest for applications such as the gasification of coal where reactions of gas mixtures, as well as gas–particle heat and mass transfer take place. Knowledge of the hydrodynamics of the gas within the freeboard can be hard to characterise, especially the detailed behaviour of gases escaping from bubbles that erupt at the bed surface. In the present study, experiments were conducted on a rectangular three-dimensional gas–solid fluidised bed. The experiments used a particle imaging velocimetry (PIV) measurement technique to visualise and measure the gas flow within the freeboard after a single bubble eruption. A computational study was carried out using Eulerian–Eulerian, kinetic theory of granular flow approach with a quasi-static flow model and with LES used to account for gas turbulence. Results from a three dimensional simulation of the experimental fluidised bed were compared with experimental velocity profiles of gas flow in the freeboard of the gas–solid fluidised bed after a bubble eruption. The CFD simulations showed a qualitative agreement with the formation of the gas vortices as the bubble erupted. Consistent with experimental findings the CFD simulations showed the generation of a pair of vortices. However, the simulations were unable to demonstrate downward flow at the centre of the freeboard due to particles in free fall after a bubble eruption event was observed in the experiments. Velocity profiles from the CFD data are in reasonably good agreement with the characteristic trends observed in the experiments, whereas the CFD model was able to predict the gas vortices phenomena and the velocity magnitudes were over-predicted.     

提升管内气固流动行为的数值模拟

应用计算流体力学软件Fluent,对空气为连续相、固相为催化裂化反应催化剂的循环流化床提升管内的气固流动行为进行模拟。采用用户自定义函数引入颗粒与壁面的恢复系数和颗粒的镜面反射系数,对颗粒在边壁处的部分滑移运动进行描述。采用不同的计算动力学模型及参数,数值模拟了径向颗粒浓度、轴向床层压降的空间分布,以及用以描述颗粒脉动动能的颗粒温度与固含率的关系,并与文献报道的实验和数值模拟结果进行对比分析。结果表明,选取的颗粒动力学理论模型及参数、颗粒部分滑移边界条件及气固曳力模型,可计算得到合理的颗粒轴向及径向分布,验证了提升管中存在典型的径向环核流动结构和轴向压降分布。进一步分析表明固含率显著影响颗粒温度,当固含率为0.05~0.1,颗粒温度存在转折区。     

Computational Fluid Dynamics Study of a Styrene Polymerization Reactor

The computational fluid dynamics (CFD) approach was adopted to simulate benzoyl peroxide (BPO)‐initiated styrene polymerization in a laboratory‐scale continuous stirred‐tank reactor (CSTR). The CFD results revealed the effects of non‐homogeneity and the short‐circuiting of the unreacted styrene and initiator on the reactor performance. The study also investigated the effects of the impeller speed and the residence time on the conversion and the flow behavior of the system. The CFD simulation showed that intense mixing remained confined to a small region near the impeller. With increasing impeller speed, it was found that the perfectly mixed region near the impeller expanded, thus reducing non‐homogeneity. Different contours were generated and exhibited the effect of the mixing parameters on the propagation rate and styrene conversion. The monomer and initiator conversions predicted with the CFD model were compared to those obtained with a CSTR model. The CFD model accounts for the non‐ideality behavior of the polymerization reactor, and hence conversion predictions are more realistic.     

Modeling of Gas‐Particle Turbulent Flow in Spout‐Fluid Bed by Computational Fluid Dynamics with Discrete Element Method

Gas and solid turbulent flow in a cylindrical spout‐fluid bed with conical base were investigated by incorporating various gas‐particle interaction models for two‐way coupling simulation of discrete particle dynamics. The gas flow field was computed by a k‐ϵ two‐equation turbulent model, the motion of solid particles was modeled by the discrete element method. Drag force, contact force, Saffman lift force, Magnus lift force and gravitational force acting on individual particles were considered in the mathematical models. Calculations on the cylindrical spout‐fluid bed with an inside diameter of 152 mm, a height of 700 mm, a conical base of 60° and the ratio of void area of 3.2 % were carried out. Based on the simulation, the gas‐solid flow patterns at various spouting gas velocities are presented. Besides, the changes in particle velocity, particle concentration, collision energy, particle and gas turbulent intensities at different proportions of fluidizing gas to total gas flow are discussed.     

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     

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

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

CFD Simulation of Mixing and Segregation of Binary Solid Mixtures in a Dense Fluidized Bed

The mixing and segregation behaviour of binary solid mixtures has been extensively studied through various experiments, while accurate CFD simulations are difficult to achieve due to process complexity and a lack of reliable constitutive relations. In this study, CFD simulations of a dense fluidized bed with glass and polystyrene particles were performed in order to identify a universal set of simulation parameters and models for simulating binary mixtures with different mixed and segregation behaviour. Through a comparison to experimental data, it was found that the EMMS drag model coupled with the Ma-Ahmadi solid pressure and radial distribution models predicted more a reasonable axial distribution of solid phases than the Syamlal O'Brien drag model coupled with the Lun et al. solid pressure and radial distribution models. The increase in the solid-solid drag further improved the simulation results.     

Simulation of the Hydrodynamics and Drying in a Spouted Bed Dryer

Gas-particle flow behavior in a spouted bed of spherical particles was simulated using the Eulerian-Eulerian two-fluid modeling approach, incorporating a kinetic-frictional constitutive model for dense assemblies of the particulate solid. The interaction between gas and particles was modeled using the Gidaspow drag model and the predicted hydrodynamics is compared with published experimental data. To investigate drying characteristics of particulate solids in axisymmetric spouted beds, a heat and mass transfer model was developed and incorporated into the commercial computational fluid dynamics (CFD) code FLUENT 6.2. The kinetics of drying was described using the classical and diffusional models for surface drying and internal moisture drying, respectively. The overall flow patterns within the spouted bed were predicted well by the model; i.e., a stable spout region, a fountain region, and an annular downcomer region were obtained. Calculated particle velocities and concentrations in the axisymmetric spouted bed were in reasonable agreement with the experimental data of He et al. (Can. J. Chem. Eng. 1994a, 72:229; 1994b, 72:561). Such predictions can provide important information on the flow field, temperature, and species distributions inside the spouted bed for process design and scale-up.     

Computationally efficient incorporation of microkinetics into resolved-particle CFD simulations of fixed-bed reactors

A method is developed to couple microkinetics with fluid flow in a fixed bed and transport inside the catalyst particles using computational fluid dynamics. Initially, the microkinetics model is solved for a wide range of different temperatures, and partial pressures. Next, reaction rates are mapped into quadratic multivariate splines. Splines coefficients are then imported into our user-defined function, and consequently the reaction rates are evaluated at each iteration simultaneously with the CFD simulations. This method has been applied to our solid particle model to investigate the effects of fluid flow, transport and elementary reaction steps on each other for ethylene and methanol partial oxidations. Reaction rates of all elementary steps as well as species surfaces sites and compositions are evaluated inside the particle. The suggested method can couple complex reaction mechanisms with detailed CFD simulations without increasing the computational time compared with global kinetics methods.     

A CFD model for predicting the heat transfer in the industrial scale packed bed

Compared to the traditional lumped-parameter model,computational fluid dynamics (CFD) attracted more attentions due to facilitating more accurate reactor design and optimization methods when analyzing the heat transfer in the industrial packed bed.Here,a model was developed based on the CFD theory,in which the heterogeneous fluid flow was resolved by considering the oscillatory behavior of voidage and the effective fluid viscosity.The energy transports in packed bed were calculated by the convection and diffusion incorporated with gaseous dispersion in fluid and the contacting thermal conductivity of packed particles in solids.The heat transfer coefficient between fluid and wall was evaluated by considering the turbulence due to the packed particles adjacent to the wall.Thus,the heat transfer in packed bed can be predicted without using any adjustable semi-empirical effective thermal conductivity coefficient.The experimental results from the literature were employed to validate this model.     

Convective heat transfer and solid conversion of reacting particles in a copper(II) chloride fluidized bed

This paper develops a predictive model of convective heat transfer and conversion of cupric chloride particles in a fluidized bed reactor of a copper–chlorine (Cu–Cl) cycle of thermochemical hydrogen production. The hydrolysis reaction of particles in the fluidized bed is endothermic and it requires excess steam for complete conversion of cupric chloride solid. The excess steam supply may be used for partial heat supply to the endothermic reaction, and also to avoid defluidization in the bed. To avoid defluidization, the change of gas flow in the bed due to the reaction should be minimized at a given operating condition. The model predicts the maximum possible steam inlet temperature, steam conversion, amount of partial heat supply, and also gas flow rates through the bed to avoid defluidization. The new model presents significant new insight by analyzing the hydrodynamic and mass transport processes, considering the equilibrium limitation on the conversion of cupric chloride solid. The model results indicate that the chemical reaction requires a high mole ratio of steam for complete conversion of cupric chloride particles. The maximum steam conversion is limited by temperature, pressure, and the presence of hydrogen chloride gas. The maximum conversion of steam at 400 °C is 3.75% and it requires excess steam of 12.8 moles per unit mole of cupric chloride solid for complete conversion of solid. The heat supply by steam for the reaction, as well as raising the solid feed to the reaction temperature, varies with reaction temperature. The paper also adds significant new insight by analyzing the steam flow requirement in terms of temperature, conversion rate, and quality of fluidization. Additional new results are presented and applications discussed for the Cu–Cl cycle of nuclear hydrogen production.     

Discrete particle simulation of gas fluidization of ellipsoidal particles

Fluidization is widely used in industries and has been extensively studied, either experimentally or theoretically, in the past decades. In recent years, a coupled simulation approach of discrete element method (DEM) and computational fluid dynamics (CFD) has been successfully developed to study the gas–solid flow and heat transfer in fluidization at a particle scale. However, to date, such studies mainly deal with spherical particles. The effect of particle shape on fluidization is recognized but not properly quantified. In this paper, the CFD–DEM approach is extended to consider the fluidization of ellipsoidal particles. In the simulation, particles used are either oblate or prolate, with aspect ratios varying from very flat (aspect ratio=0.25) to elongated (aspect ratio=3.5), representing cylinder-type and disk-type shaped particles, respectively. The commonly used correlations to determine the fluid drag force acting on a non-spherical particle are compared first. Then the model is verified in terms of solid flow patterns. The effect of aspect ratio on the flow pattern, the relationship between pressure drop and gas superficial velocity, and microscopic parameters such as coordination number, particle orientation and force structure are investigated. It is shown that particle shape affects bed permeability and the minimum fluidization velocity significantly. The coordination number generally increases with aspect ratio deviating from 1.0. The analysis of particle orientations shows that the bed structures for ellipsoids are not random as that for spheres. Oblate particles prefer facing upward or downward while prolate particles prefer horizontal orientation. Spheres have the largest particle–particle contact force and fluid drag force under the comparable conditions. With aspect ratio deviating from 1.0, particle–particle interaction and fluid drag become relatively weak. The proposed model shows a promising method in examining the effect of particle shape on different flow behaviour in gas fluidization.     

CFD modeling of the gas-particle flow behavior in spouted beds

Gas-particle flow behavior in a cylindrical spouted bed and a three dimensional spout-fluid bed of spherical particles was simulated using the Eulerian-Eulerian two-fluid modeling approach, incorporating a kinetic-frictional constitutive model for dense assemblies of the particulate solid. The interaction between gas and particles was modeled using the Gidaspow drag model and the predicted hydrodynamic characteristics are compared with published experimental data. The overall flow patterns within the cylindrical spouted bed were predicted well by the model, i.e. a stable spout region, a fountain region and an annular downcomer region were correctly predicted by the model. The flow instabilities which develop in the spout-fluid bed are along with discussion of the mechanisms leading to instabilities. Bubble formation and motion of the bubbles inside the spout-fluid bed are also described. Such predictions can provide important information on the flow field within the spouted beds for process design and scale-up.