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

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

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     

外加声场微小流化床中A类与C类颗粒流化的数值模拟

文章基于多相流欧拉模型,应用计算流体动力学（CFD）软件Fluent 6.3,在有无外加声场的情况下,分别对A类与C类颗粒在微小流化床（内径10mm,床高60 mm）中的颗粒分布进行了数值模拟,得到了一系列的固含率云图.模拟结果发现,在通入气体的初始阶段,C类颗粒流化出现了活塞流现象;加上声场后,上述现象消失.无声场时A类颗粒床层有较大气泡;引入声场后气泡变小.引入声场使得C类颗粒与A类颗粒有相似的固含率云图.     

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     

Comparison and evaluation of interphase momentum exchange models for simulation of the solids volume fraction in tapered fluidised beds

An Eulerian computational fluid dynamics (CFD) model with granular flow extension was used to simulate a gas–solid fluidised bed in a tapered reactor. Various drag coefficient models were evaluated, which are used to calculate the drag force, describing the momentum transfer between the gas and solid phases. Comparison and evaluation between time-averaged solids volume fractions obtained from experiments and from simulations with several drag coefficient models were made. The predicted results obtained by the different drag models were verified using experimental data of Depypere et al. (2009). Initial results using a 2-phase Eulerian model showed poor agreement with experimental results. However, extending the Eulerian model to include 3 solid phases—with different mean particle diameter per phase in order to account for the particle size distribution of the fluidised solid material—yielded good agreement with experimental results. Furthermore, quantitative analyses showed that the modified Gidaspow drag model gave the best agreement between CFD simulations and experimental data.     

Experimental validation of CFD simulations of a lab‐scale fluidized‐bed reactor with and without side‐gas injection

Fluidized‐bed reactors are widely used in the biofuel industry for combustion, pyrolysis, and gasification processes. In this work, a lab‐scale fluidized‐bed reactor without and with side‐gas injection and filled with 500–600 μm glass beads is simulated using the computational fluid dynamics (CFD) code Fluent 6.3, and the results are compared to experimental data obtained using pressure measurements and 3D X‐ray computed tomography. An initial grid‐dependence CFD study is carried out using 2D simulations, and it is shown that a 4‐mm grid resolution is sufficient to capture the time‐ and spatial‐averaged local gas holdup in the lab‐scale reactor. Full 3D simulations are then compared with the experimental data on 2D vertical slices through the fluidized bed. Both the experiments and CFD simulations without side‐gas injection show that in the cross section of the fluidized bed there are two large off‐center symmetric regions in which the gas holdup is larger than in the center of the fluidized bed. The 3D simulations using the Syamlal‐O'Brien and Gidaspow drag models predict well the local gas holdup variation throughout the entire fluidized bed when compared to the experimental data. In comparison, simulations with the Wen‐Yu drag model generally over predict the local gas holdup. The agreement between experiments and simulations with side‐gas injection is generally good, where the side‐gas injection simulates the immediate volatilization of biomass. However, the effect of the side‐gas injection extends further into the fluidized bed in the experiments as compared to the simulations. Overall the simulations under predict the gas dispersion rate above the side‐gas injector. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

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     

65 t/h高低差速循环流化床流动特性模拟

对一台65 t/h高低差速循环流化床炉内流动特性进行二维数值模拟。采用基于颗粒动力学理论的欧拉双流体模型来描述气固流动,湍流模型、气固曳力模型和不同粒径颗粒间曳力模型分别采用RNG k-ε per phase模型、Gidaspow模型和Schiller-naumann模型,并应用商业计算流体力学软件Fluent进行数值计算,得到炉内颗粒速度分布、压力分布和颗粒浓度分布,并将压力分布与实测值进行对比。在欧拉双流体模型中分别采用单粒径固相模型和多粒径固相模型,并对模拟结果进行对比分析。结果表明,单粒径固相模型能够较好预测高低差速循环流化床炉内流动特性,为其优化设计、运行及大型化提供了理论依据。     

Development and test of CFD–DEM model for complex geometry: A coupling algorithm for Fluent and DEM

CFD–Discrete Element Method (DEM) model is an effective approach for studying dense gas–solid flow in fluidized beds. In this study, a CFD–DEM model for complex geometries is developed, where DEM code is coupled with ANSYS/Fluent software through its User Defined Function. The Fluent Eulerian multiphase model is employed to couple with DEM, whose secondary phase acts as a ghost phase but just an image copy of DEM field. The proposed procedure preserves phase conservation and ensures the Fluent phase-coupled SIMPLE solver work stable. The model is used to simulate four typical fluidization cases, respectively, a single pulsed jet fluidized bed, fluidized bed with an immersed tube, fluidization regime transition from bubbling to fast, and a simplified two-dimensional circulating fluidized bed loop. The simulation results are satisfactory. The present approach provides an easily implemented and reliable method for CFD–DEM model on complex geometries.     

湍动流化床过渡段固含率分布特征的实验及数值模拟

在湍动流化床中,过渡段对于包括甲醇制烯烃在内的气固催化快反应有着重要的作用。采用PV6D反射型光纤探针对内径95 mm的湍动流化床内过渡段的固含率分布和脉动参数进行了测量,分别考察了表观气速和静床高的影响,并采用修正的基于颗粒动力学的三段曳力双流体模型进行模拟。实验表明,湍动流化床过渡段中固含率的轴向分布呈现S型和指数型两种类型,固含率轴向与径向分布都在过渡段内出现最大梯度,表明过渡段中固体浓度分布比稀相段和密相段更不均匀。表观气速和静床高的变化将导致S型和指数型分布的相互转变,并且对过渡段底部与壁面附近的固体高浓度区影响最为显著。局部固含率脉动概率密度分布表明,在静床高较小时,随着气速的增大,床层下部气含率最大值位置将从中心区移动至环隙区,呈现气含率的双峰型分布。本文提出的修正三段曳力模型考虑了颗粒团聚的影响,对过渡段中分布板影响区之外的固含率分布均能较好地模拟。     

Profiles of solid fraction and heterogeneous phase structure in a gas–solid airlift loop reactor

The time-averaged and transient local solid fractions in a gas–solid airlift loop reactor (ALR) were investigated systematically by experiments and CFD simulations. To demonstrate the macro-flow pattern, the time-averaged local solid fractions in four regions of the ALR were measured by optical fiber probe under the conditions of different superficial gas velocities and particle circulation fluxes. The experimental results show that the lateral distribution of time-averaged local solid fraction is a core-annulus or heterogeneous structure in the three regions (draft tube, bottom region, particle diffluence region), but a uniform lateral distribution in the annulus. The operating conditions have different effects on the lateral distribution of time-averaged local solid fraction in each region. In the CFD simulation, a modified Gidaspow drag model considering the formation of particle clusters was incorporated into the Eulerian–Eulerian CFD model with particulate phase kinetic theory to simulate and analyze the transient local solid fraction and the two-phase micro-structures in the gas–solid ALR. The predicted values of solid fraction were compared with the experimental results, validating the drag model. The contours of transient flow field indicate that the flow field of the ALR should be divided into five flow regions, i.e., draft tube, annulus, bottom region, particle diffluence region and constrained back-mixing region, which further improves the understanding of the airlift reactor where only four divisions were determined from the experiments. The transient local solid fraction and its probability density function profoundly reveal the two-phase micro-structures (dilute phase and emulsion phase or cluster phase in the constrained back-mixing region) and explain the heterogeneous phenomenon of solid fraction in the ALR. The dilute phase tends to exist in the center of bed, while the emulsion phase mainly appears in the wall region. The results also indicate that the gas–solid ALR has the common characteristic of aggregative fluidization similar to that in normal fluidized beds. The simulated two-phase transient micro-structures provide the appropriate explanations for the experimental core-annulus macro-structures of time-averaged local solid fraction.     

Geldart A类颗粒节涌床气固流动特性的实验及模拟

针对气固节涌床,在实验基础上,基于欧拉?欧拉双流体模型结合颗粒动力学理论,考虑Geldart A类颗粒聚团对气固间曳力的影响,采用修正后的Gidaspow曳力模型对气固节涌床进行数值模拟。结果表明,通过与实验结果及经验公式进行对比,修正的模型可准确合理地模拟流化床内节涌特性。表观气速0.09 m/s≤Ug≤0.39 m/s时,床层内部压力脉动标准偏差随表观气速增加而增加,流型由鼓泡转变为节涌直至节涌程度最大,床内气固流动主要受轴对称栓运动特性影响,床内压降、床层膨胀比、气栓平均上升速度、最大轴对称栓长度随表观气速增加而增加,最大轴对称栓产生位置随表观气速增加而降低;Ug>0.39 m/s后,床内压力脉动标准偏差随表观气速增加而降低,节涌程度降低至向湍动流态化流型转变,床内气固流动主要受壁面栓运动特性影响,增加表观气速,节涌床内压降变化幅度较小,气栓平均上升速度增加幅度加大,床层膨胀比及最大轴对称栓长度降低,最大轴对称栓产生的位置略有升高。     

CFD Modeling of Microwave-Assisted Fluidized Bed Drying of Moist Particles Using Two-Fluid Model

The drying behavior of moist spherical particles in a microwave-assisted fluidized bed dryer was simulated. The two-fluid Eulerian model incorporating the kinetic theory of granular flow was applied to simulate the gas–solid flow. The simulations were carried out using the commercial computational fluid dynamics (CFD) package Fluent 6.3.26. The effects of different levels of microwave power densities as well as initial gas temperature on the prediction of solids moisture content, gas temperature, and gas absolute humidity were investigated. The effect of microwaves was incorporated into calculations using a concatenated user-defined function (UDF). The simulation results were compared with experimental data obtained from drying of soybeans in a pilot-scale microwave-assisted fluidized bed dryer and reasonable agreement was found. The mean relative deviation for prediction of solids moisture content, gas temperature, and gas absolute humidity were less than 3, 10, and 5%, respectively. Further work is needed to validate the proposed model for large-scale plants.     

Modelling of binary fluidization of coal and ash and validation using radioactive particle tracking and densitometry

Binary fluidization finds wide application in a variety of gas–solid catalytic and non-catalytic industrial fluidization systems. In the present study, a three-dimensional transient computational fluid dynamics (CFD) model was used to model the binary fluidization of coal and ash in a laboratory-scale cold flow fluidized bed. In parallel, phase velocity measurements using radioactive particle tracking (RPT) and gamma-ray densitometry were performed, which provided a rich database for validation of the CFD model. RPT being a time-resolved Lagrangian technique, it was possible to extract velocity fluctuations and their correlations in addition to the mean velocity profiles. The latter provided additional validation for the CFD model, in addition to the typical validation that is done with time-averaged profiles of phase velocity and volume fraction. The robust validation procedure opens up the possibility of expanding this model to a pilot plant-scale fluidized bed.     

Direct numerical simulations of dynamic gas‐solid suspensions

Direct numerical simulation results for gas flow through dynamic suspensions of spherical particles is reported. The simulations are performed using an immersed boundary method, with careful correction for the grid resolution effect. The flow systems we have studied vary with mean flow Reynolds number, solids volume fraction, as well as particle/gas density ratio. On the basis of the simulation results, the effect of particle mobility on the gas‐solid drag force is analyzed and introduced into the existing drag correlation that was derived from simulations of stationary particles. This mobility effect is characterized by the granular temperature, which is a result of the particle velocity fluctuation. The modified drag correlation is considered so‐far the most accurate expression for the interphase momentum exchange in computational fluid dynamics models, in which the gas‐solid interactions are not directly resolved. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1958–1969, 2016     

Catalytic Ozone Decomposition in a Gas‐Solids Circulating Fluidized‐Bed Riser

Computational fluid dynamics (CFD) modeling of the catalytic ozone decomposition reaction in a circulating fluidized‐bed (CFB) riser, using iron‐impregnated fluid catalytic cracking particles as catalyst, is carried out. The catalytic reaction is defined as a one‐step reaction, and the reaction equation is modified by with respect to the particle surface area, A_p, and an empirical coefficient. The Eularian‐Eularian method with the kinetic theory of granular flow is used to solve the gas‐solids two‐phase flow in the CFB riser. The simulation results are compared with experimental data, and the reaction rate is modified by using an empirical coefficient, to provide better simulation results than the original reaction rate. Moreover, the particle size has great effects on the reaction rate. The generality of the CFD model is further validated under different operating conditions of the riser.     

Euler multiphase-CFD simulation on a bubble-driven gas–liquid–solid fluidized bed

An integrated flow model was developed to simulate the fluidization hydrodynamics in a new bubble-driven gas–liquid–solid fluidized bed using the computational fluid dynamic (CFD) method. The results showed that axial solids holdup is affected by grid size, bubble diameter, and the interphase drag models used in the simulation. Good agreements with experimental data could be obtained by adopting the following parameters: 5 mm grid, 1.2 mm bubble diameter, the Tomiyama gas–liquid model, the Schiller–Naumann liquid–solid model, and the Gidaspow gas–solid model. At full fluidization state, an internal circulation of particles flowing upward near the wall and downward in the centre is observed, which is in the opposite direction compared with the traditional core-annular flow structure in a gas–solid fluidized bed. The simulated results are very sensitive to bubble diameters. Using smaller bubble diameters would lead to excessive liquid bed expansions and more solid accumulated at the bottom due to a bigger gas–liquid drag force, while bigger bubble diameters would result in a higher solid bed height caused by a smaller gas–solid drag force. Considering the actual bubble distribution, population balance model (PBM) is employed to characterize the coalescence and break up of bubbles. The calculated bubble diameters grow up from 2–4 mm at the bottom to 5–10 mm at the upper section of the bed, which are comparable to those observed in experiments. The simulation results could provide valuable information for the design and optimization of this new type of fluidized system.     

Numerical and experimental investigation of a fluidized bed chamber hydrodynamics with heat transfer

The hydrodynamics and heat transfer of a gas-solid fluidized bed chamber was investigated by computational fluid dynamic (CFD) techniques. A multifluid Eulerian model incorporating the kinetic theory for solid particles was applied to simulate the unsteady state behavior of this chamber. For momentum exchange coefficients, Syamlal-O’Brien drag functions were used. A suitable numerical method that employed finite volume method was applied to discretize the equations. The simulation results also indicated that small bubbles were produced at the bottom of the bed. These bubbles collided with each other as they moved upwards forming larger bubbles. Also, the solid particle temperature effect on heat transfer and hydrodynamics was studied. Simulation results were compared with the experimental data in order to validate the CFD model. Pressure drops and mean gas temperature predicted by the simulations at different positions in the chamber were in good agreement with experimental measurements at gas velocities higher than the minimum fluidization velocity. Furthermore, this comparison showed that the model could predict hydrodynamics and heat transfer behaviors of gas solid fluidized bed reasonably well.     

喷动床内气固二相流动行为的计算模拟

采用双流体模型结合颗粒动理学理论对喷动床内气固二相流体流动行为进行了计算模拟研究。模型中运用颗粒动理学理论描述颗粒相应力封闭流体控制方程,使用Gidaspow曳力模型描述气固相间作用。喷动床内颗粒在浓相区的体积分数很大,采用Schaeffer′s模型描述颗粒间的摩擦应力。模拟计算结果表明,喷动床内分喷射区、喷泉区、环隙区3个区域,在射流入口处形成一个瓶颈。模拟计算得到的颗粒速度和空隙度分布与实验数据进行比较,计算结果与实验结果吻合较好。