喷动床内气固两相流体动力行为的数值模拟

引　言喷动床被广泛应用于不同工业领域中 ,如石油裂解反应 -再生器、煤和农业废弃物气化和燃烧 ,喷动床还被应用于粮食和药品的干燥等[1] .因此 ,喷动床设计应满足不同应用的要求 .喷动床内气相反应物的反应时间和停留时间依赖于床体几何结构和运行参数 .尽管已有许多的实验对喷动床内气固两相流动进行了研究 ,得到了喷射区、环形区和喷泉区内的气固两相流动流体动力特性 ,然而由于喷动床内气固两相流动的复杂性 ,人们对床体几何结构和运行参数对喷动床动力学的影响至今并不清楚 .因此 ,床体几何结构和运行参数等对喷动床动力学的影响成为…     

CFD simulation of gas–solid flow in a spouted bed with a non-porous draft tube

A multi-fluid Eulerian–Eulerian approach incorporating the kinetic theory of granular flow was used to simulate a spouted bed containing non-porous draft tube. Drag function and coefficient of restitution were investigated. Solid and gas velocity vector, gas flow rate in annulus and spout regions and longitudinal pressure distribution were evaluated. In addition, the effects of the entrainment height and the draft tube diameter were studied. Simulation indicates the formation of three regions namely, annulus, spout and fountain; similar to a conventional spouted bed. Current model predicts acceptable results in both spout and annulus regions. Simulation results indicate that the model can be employed for both mono-size and multi-size particles reasonably. This paper provides useful basis for further works on understanding gas–solid flow mechanism in spouted beds containing a non-porous draft tube.     

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

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

Numerical Simulations of the Effect of Conical Dimension on the Hydrodynamic Behaviour in Spouted Beds

The flow behaviours of gas‐solids were predicted by means of a hydrodynamic model of dense gas‐solid flow in spouted beds. Constitutive equations describing the particulate solids pressure and viscosity were implemented into a hydrodynamic simulation computer program. The effect of operating conditions (inclined angle and gas spouting velocity) on particle velocity and concentration in the spout, annulus and fountain regions were numerical studied. Both vertical and horizontal particle velocities increased with increasing spouting gas velocity. The diameter of the spout increases with decreasing the inclination angle. As the inclination angle is set greater than 60°, the spout cross‐section starts becoming bottlenecked, limiting the upwards flow of solids.     

Vapour phase chemical reaction in spouted beds: Verification of theory

The performance of a 0.15 m diameter spouted bed reactor has been studied for decomposition of ozone on iron oxide catalyst. The variables include bed depth, catalyst particle size and spouting gas velocity.The spouted bed conversion data are compared against predictions from theoretical models proposed previously by Mathur and Lim who assumed vertically upward flow of annular gas in their first, one-dimensional model, and flow along curved streamlines in a subsequent, more rigorous, streamtube model. Experimental results show good agreement with conversions predicted by both models, but the latter model is expected to be superior for design of larger reactors. Additional experiments to discriminate between the two models and to establish their range of validity are planned.     

Motion of solids in spouted beds

Solids movement in a spouted bed has been investigated by means of a radioactively marked particle technique, using a scintillation counter as a detector. The average particle velocity in the spout was found to depend on gas flow rate, gas inlet geometry, axial elevation, bed height, and column diameter. The solids circulation rate is shown to be proportional to the gas flow rate and, further, a function of column and gas inlet geometry, density of fluidising gas and particle diameter. Empirical correlations for both the particle velocity in the spout and the solids circulation rate as a function of axial elevation, are presented.     

Effects of Draft Tubes on Particle Velocity Profiles in Spouted Beds

The vertical particle velocity profiles in a full‐column cylindrical conical spouted bed, with or without a draft tube, are measured using a fibre optic probe system. The profiles have different characteristics for a draft tube spouted bed (DTSB) than for a conventional spouted bed (CSB). The spout of a CSB consists of a central flow where particle velocities fit exponential distributions, and a boundary layer where particle velocities are nearly uniform. The spout of a DTSB has no boundary layer and its radial particle velocity profiles are approximately linear. The particle velocities in the spout of a DTSB increase when superficial gas velocity increases, draft tube diameter decreases, or when entrainment height decreases. A kinematic model has been used to simulate the granular flow in the annulus of a CSB and DTSB, and they are compared with the experiments. The particle velocities in the annulus of a DTSB are much lower than that of a CSB. Their radial profiles are also different with a CSB. The dependence of particle velocities in the annulus of a DTSB on superficial gas velocity, draft tube diameter, and entrainment height are also discussed. One concludes that the draft tube diameter and entrainment height are two key factors for the solid circulation rate of a DTSB.     

Analysis of the dynamics of paste drying in a spouted bed

Although there are some models available in the literature for paste drying in spouted beds, few of them have focused on transient analysis of dynamical systems. Our objective was to integrate experiments and simulations of a dynamic model to investigate the transient response to disturbances and interruptions in the feed flow during paste drying in a spouted bed with inert particles. The spouted bed consisted of a cylindrical column with 50.0?cm of height and 20.0?cm of diameter. Drying tests were performed at inlet gas temperatures of 70?°C and 100?°C and inlet air flow 30% above the minimum spouting velocity. A 5% w/w suspension of calcium carbonate was used as paste material, and glass spheres of 2.2?mm were used as inert materials. Different patterns of step function changes were tested in the paste feed flow rate. A lumped parameter model was used to predict mass and heat transfer during the drying. Experiments and simulations were in good agreement.     

Numerical simulation of two-dimensional spouted bed with draft plates by discrete element method

A discrete element method (DEM)-computational fluid dynamics (CFD) two-way coupling method was employed to simulate the hydrodynamics in a two-dimensional spouted bed with draft plates. The motion of particles was modeled by the DEM and the gas flow was modeled by the Navier-Stokes equation. The interactions between gas and particles were considered using a two-way coupling method. The motion of particles in the spouted bed with complex geometry was solved by combining DEM and boundary element method (BEM). The minimal spouted velocity was obtained by the BEM-DEM-CFD simulation and the variation of the flow pattern in the bed with different superficial gas velocity was studied. The relationship between the pressure drop of the spouted bed and the superficial gas velocity was achieved from the simulations. The radial profile of the averaged vertical velocities of particles and the profile of the averaged void fraction in the spout and the annulus were statistically analyzed. The flow characteristics of the gas-solid system in the two-dimensional spouted bed were clearly described by the simulation results.     

三维喷动床内异径干湿颗粒混合特性数值模拟

基于计算流体力学-离散单元法,建立了三维喷动床内气固两相流数学模型,采用Fortran语言编制了并行数值模拟程序。对三维喷动床内两种不同直径的干颗粒及湿颗粒的混合特性进行了数值模拟,并从颗粒角度分析了双组分颗粒的运动机制。利用Lacey混合指数对床内整体以及特定区域的混合程度进行了定量分析,并研究了液桥体积、颗粒密度比以及表观气速对异径颗粒混合的影响。结果表明：在单孔射流喷动床内,干湿两种颗粒流动方式相似,湿颗粒无明显的聚团现象;液桥力对小直径的颗粒影响较大,使不同直径湿颗粒速度差减小;环隙区内颗粒的混合是影响整床颗粒混合的关键因素;液桥体积对颗粒混合的影响较大,对颗粒密度比以及表观气速的影响有限。     

Hydrodynamics of a gas-driven gas-liquid-solid spouted bed with a draft tube

A cylindrical gas-liquid-solid spouted bed, driven exclusively by gas flow, has been developed with a high potential for use in biochemical processes, such as a biological wastewater treatment. A plexiglass column with a 152 mm inner diameter was used in combination with a 53 mm inner diameter plexiglass draft tube. Three particle types were studied with densities ranging from 1044 kg/m³-1485 kg/m³ and average particle sizes ranging from 0.7-2.5 mm. Four flow regimes were observed when increasing the gas velocity, including fixed bed, semispouted bed, full spouted bed, and internal circulating fluidized bed. The transition gas velocities between those regimes were experimentally measured and termed as minimum spouting velocity, full spouting velocity, and minimum circulating velocity, respectively. A measurement of the downward particle flux in the annulus was used to identify the minimum spouting velocity, while the particle velocity and dense phase retraction in the annulus were monitored for the full spouting and minimum circulating velocities. All regime transition velocities increased with more dense particles and longer draft tubes. The minimum spouting velocity and full spouting velocity were not affected when varying the nozzle-tube gap, while the minimum circulating velocity increased with longer nozzle-tube gaps. Experiments without a draft tube were carried, though the spouting stability was significantly reduced without the draft tube.     

二维导流管喷动床离散颗粒动力学模拟

采用离散单元法(DEM)-计算流体力学(CFD)双向耦合数值方法对二维导流管喷动床进行了模拟,颗粒的运动通过DEM模型描述,而气体的运动用Navier-Stokes方程进行求解,气体和固体颗粒之间的相互作用通过曳力形式传递。文中将DEM和边界元方法(BEM)结合起来解决颗粒在具有复杂边界设备内的运动。通过采用BEM+DEM-CFD相结合的方法进行模拟计算,得到了喷动床的最小喷动速度,研究了不同表观气速下床内的流型,得到了二维导流管喷动床的床层压降与表观气速的关系,统计分析了喷射区、环隙区内颗粒的运动速度和相应的空隙率,全面地描述了二维导流管喷动床内的气固流动特征。     

Discrete element method simulation of three-dimensional conical-base spouted beds

A discrete element method (DEM) simulation of three-dimensional conical-base spouted beds is presented. The overall height and diameter of the vessel are 0.5 and 0.15 m, respectively, and the nozzle diameter is 0.02 m. The inclined angle of the conical section varies from 0 to 60 degrees. The gas flow is described by the continuity and Navier-Stokes equations and solved by a finite difference method of second order accuracy in space and time. For gas-particle interaction, the Ergun equation (for void fraction smaller than 0.8) and the Wen-Yu model (for void fraction of 0.8 and above) are employed. A new method for treatment of the boundary condition for 3-D gas flow along the cone surface is proposed. This boundary condition satisfies both the continuity and momentum-balance requirements for the gas phase. Usefulness of the present simulation for studying gas flow pattern and particle motion in conical-base spouted beds is demonstrated. The effects of the inclined angle and draft tube on gas and particle flow in spouted beds are discussed.     

喷动流化床流动特性数值模拟

为了捕捉喷动流化床中微观层次上的颗粒运动信息,建立了基于CFD的二维非稳态喷动流化床欧拉-欧拉两相流模型。分析了不同流化气速对喷动流化床气固流动特性的影响,即不同工况下的炉内压力降、颗粒浓度、床内空隙率分布、气体速度分布和固体颗粒速度分布。数值模拟研究结果表明:随流化气速的增大,压降和炉内平均空隙率逐渐增大,密相床层高度逐渐增加,沿着轴向方向的气体流量增大,喷动气的射流深度逐渐增加,同时射流半径也逐渐增加。     

基于动态双重网格下喷动床滞止区流动特性CFD-DEM数值模拟

基于Fortran语言自行开发了基于动态双重网格方法下的喷动床内气固两相流动的CFD-DEM方法,同时开展了喷动床内径向混合实验与模拟研究,又结合单网格方法对喷动床内0~2.0 s内的滞止区特性进行对比分析,验证了动态双网格方法计算结果的准确性。然后利用动态双网格方法对不同进口气速下和不同初始堆积高度下的喷动床进行数值模拟研究,对滞止区颗粒流动过程进行追踪,结果表明：径向混合实验结果与数值模拟结果有很好一致性;在喷动床内存在一定的滞止区,滞止区内的颗粒流动性较差;初始堆积高度不变,随着进口速度的增加,滞止区高度下降速率和向喷口延伸速度无明显变化;进口速度不变,随着初始堆积高度的增加,滞止区颗粒下降速度随之增加,但其向喷口延伸速度逐渐变慢。     

Numerical simulation of two-dimensional spouted bed with draft plates by discrete element method

A discrete element method (DEM)-computational fluid dynamics (CFD) two-way coupling method was employed to simulate the hydrodynamics in a two-dimensional spouted bed with draft plates. The motion of particles was modeled by the DEM and the gas flow was modeled by the Navier-Stokes equation. The interactions between gas and particles were considered using a twoway coupling method. The motion of particles in the spouted bed with complex geometry was solved by combining DEM and boundary element method (BEM). The minimal spouted velocity was obtained by the BEMDEM-CFD simulation and the variation of the flow pattern in the bed with different superficial gas velocity was studied. The relationship between the pressure drop of the spouted bed and the superficial gas velocity was achieved from the simulations. The radial profile of the averaged vertical velocities of particles and the profile of the averaged void fraction in the spout and the annulus were statistically analyzed. The flow characteristics of the gas-solid system in the two-dimensional spouted bed were clearly described by the simulation results. __________ Translated from Chemical Engineering (China), 2007, 35(6): 24–28 [译自: 化学工程]     

Hydrodynamic simulations of gas–solid spouted bed with a draft tube

Flow behavior of particles in a two-dimensional spouted bed with a draft tube is studied using a continuous kinetic-friction stresses model. The kinetic stress of particles is predicted from kinetic theory of granular flow, while the friction stress is computed from a model proposed by Johnson et al. (1990). A stitching function is used to smooth from the rapid shearing viscous regime to the slow shearing plastic regime. The distributions of concentration and velocities of particles are predicted in the spouted bed with a draft tube. Simulated results compare with the vertical velocity of particles (Zhao et al., 2008) measured and in the spout bed with draft plates and solid circulation rate (Ishikura et al., 2003) measured in the spouted bed with a draft tube. The impact of the friction stress of particles on the spout, annulus, fountain and entrancement regions is analyzed in gas–solid spouted bed with a draft tube. Numerical results show that the gas flow rate through the annulus increases with the increase of the entrainment zone. The solids circulation rate decreases with the decrease of inlet gas velocity and the height of the entrainment zone. The effect of spouting gas velocity on distributions of concentration, velocity and particle circulation is discussed.     

Numerical simulations of flow behavior of gas and particles in spouted beds using frictional-kinetic stresses model

Flow behavior of gas and particles is simulated in the spouted beds using a Eulerian-Eulerian two-fluid model on the basis of kinetic theory of granular flow. The kinetic-frictional constitutive model for dense assemblies of solids is incorporated. The kinetic stress is modeled using the kinetic theory of granular flow, while the friction stress is from the combination of the normal frictional stress model proposed by Johnson and Jackson (1987) and the frictional shear viscosity model proposed by Schaeffer (1987) to account for strain rate fluctuations and slow relaxation of the assembly to the yield surface. An inverse tangent function is used to provide a smooth transitioning from the plastic and viscous regimes. The distributions of concentration, velocity and granular temperature of particles are obtained in the spouted bed. Calculated particle velocities and concentrations in spouted beds are in agreement with the experimental data obtained by He et al. (1994a, b). Simulated results indicate that flow behavior of particles is affected by the concentration of the transition point in spouted beds.     

Liquid spouting of pulp fibers in a conical vessel

Pulp fibers can be spouted in water in a conical vessel. The entities which are spouted are fiber flocs rather than individual fibers. Synthetic fibers, which do not flocculate, cannot be spouted. For comparison, rigid spherical particles were spouted with water in the same conical vessel. Liquid spouting of rigid particles was similar to gaseous spouting. For pulp spouting, the minimum spouting velocity is proportional to the mass of fibers in the bed and inversely proportional to the diameter of the inlet. For rigid particles, the minimum spouting velocity is proportional to the height of the bed and inversely proportional to the square of the diameter of the inlet. A model for the minimum spouting velocity was developed for pulp spouting.     

Particle velocity and solids circulation rate in a jet-spouted bed

Radial distribution of particle velocities, solids circulation rate and radial variation in bed voidage in a jet-spouted bed were investigated in a 0.2 m diameter conical-cylindrical column. The upward particle velocity is strongly dependent on gas velocity and its radial profile is the Gaussian. The solids circulation rates are of the same order of magnitude as those in a conventional spouted bed and are proportional to the gas velocity. Estimation of the variation in bed voidage in the radial direction revealed that the annulus voidage of about 0.90 was greater than the spout voidage of around 0.70. This tendency is contrary to that in a conventional spouted bed.