双流体模型中曳力及恢复系数对气固流动的影响

应用双流体模型CFD模拟的方法,从恢复系数和曳力两方面,研究了气固密相流化床中颗粒之间和气固相之间的相互作用对床内非均匀流动结构形成与变化的影响.计算结果表明颗粒间非弹性碰撞和气固间曳力的增大均使气固两相流动的非均匀性增大.通过比较二者对非均匀流动结构的影响,发现气固间曳力是形成非均匀流动结构的决定因素.从碰撞耗散、颗粒动能和颗粒势能的角度分析了二者的作用机理,发现恢复系数和曳力对流动结构的作用主要区别在于对颗粒团聚和床层膨胀的影响程度不同.     

Gas-particle interactions in dense gas-fluidized beds

The occurrence of heterogeneous flow structures in gas-particle flows seriously affects gas–solid contacting and transport processes in dense gas-fluidized beds. A computational study, using a discrete particle method based on Molecular Dynamics techniques, has been carried out to explore the mechanisms underlying the formation of heterogeneous flow structures. Based on energy budget analyses, the impact of non-linear drag force on the flow structure formation in gas-fluidized beds has been examined for both ideal particles (elastic collision, without inter-particle friction) and non-ideal particles (inelastic collision, with inter-particle friction). Meanwhile, the separate role of inter-particle inelastic collisions, accounted for in the model via the restitution coefficient (e) and friction coefficient (μ), has also been studied.It is demonstrated that heterogeneous flow structures exist in systems with both non-ideal particle-particle interaction and ideal particle-particle interaction. The heterogeneous structure in an ideal system, featured with looser packing, is purely caused by the non-linearity of the gas drag: the stronger the non-linearity of the gas drag force with respect to the voidage, the more heterogeneous flow structures develop. A weak dependence of drag on the voidage produces a homogenous flow structure. Collisional dissipation dramatically intensifies the formation of heterogeneous flow structures after the system equilibrium breaks. Quantitative comparisons of flow structures obtained by using various drag correlations in literature will also be reported.     

不同曳力模型及颗粒碰撞恢复系数对短接触旋流反应器内气固流场的影响

为考察曳力模型和颗粒碰撞恢复系数对短接触旋流反应器内流动特性的影响,基于双流体模型, 结合颗粒动力学理论,对反应器内气固两相流场进行模拟研究。分别采用Gidaspow、Wen & Yu和Syamlal-O’Brien 3种曳力模型, 考察颗粒速度特性以及固含率径向分布。对比分析不同曳力模型的计算结果表明,Syamlal-O’Brien模型计算结果与实验结果误差较大,Wen & Yu模型在反应器边壁附近区域的计算结果误差较大,Gidaspow模型计算结果与实验结果最为吻合。此外,颗粒碰撞恢复系数较小时,所得计算值小于实验测量值,当恢复系数为0.95时颗粒扩散效果最好,计算结果与实验数据吻合度最高。     

Study of cluster formation in dense two-phase flow using a multi-lattice deterministic model

To simulate cluster formation in dense two-phase flow, a multi-lattice deterministic trajectory (MLDT) model is developed. The actual inter-particle collision and particle motion are treated by a Lagrangian trajectory model with three sets of lattices to reduce computation time. Initially, the MLDT model is used to simulate dense gas-particle flow in a vertical channel for validation with available simulated and experimental data. Using the MLDT model, effects of particle properties and flow parameters on cluster formation are studied. It is shown that cluster formation is enhanced by the reduction of particle restitution coefficient and increase of original particle volume fraction.     

Numerical predictions of flow behavior and cluster size of particles in riser with particle rotation model and cluster-based approach

Flow behavior of particles in a circulating fluidized bed (CFB) riser is numerically simulated using a two-fluid model incorporating with the kinetic theory for particle rotation and friction stress models. The particle rotations resulting from slightly friction particle-particle collisions was considered by introducing an effective coefficient of restitution based on the kinetic theory for granular flow derived by Jenkins and Zhang [2002. Kinetic theory for identical, frictional, nearly elastic spheres. Physics of Fluids 14, 1228-1235]. The normal friction stress model proposed by Johnson et al. [1990. Frictional-collisional equations of motion for particles flows and their application to chutes. Journal of Fluid Mechanics 210, 501-535] and a modified frictional shear viscosity model proposed by Syamlal et al. [1993. MFIX Documentation and Theory Guide, DOE/METC94/1004, NTIS/DE94000087] were used as the particle frictional stress model. The drag force between gas and particle phases was modified with cluster-based approach (CBA). The flow behavior of particles and the cluster size in a riser of Wei et al. [1998. Profiles of particle velocity and solids fraction in a high-density riser. Powder Technology 100, 183-189] and Issangya et al. [2000. Further measurements of flow dynamics in a high-density circulating fluidized bed riser. Powder Technology 111, 104-113] experiments are predicted. Effects of the rotation and friction stress models on the computed results are analyzed. It is concluded that particle rotations reduce the cluster size and alter the particle flows and distributions through more particle fluctuation energy dissipations. Effects of frictional stress on flow behavior and cluster size are not significant because the particle phase in the CFB riser is not dense enough to take into account for the particle-particle contact interactions.     

颗粒特性对撞击分离器性能影响的实验与数值研究

颗粒-壁面撞击行为和气固相间作用对撞击分离器的性能具有重要影响。基于玻璃珠及煤粉的单颗粒撞击实验数据建立平均撞击恢复系数模型,采用非球形颗粒曳力模型对平板式撞击分离器的分离性能和气固流动开展数值研究。结果表明,采用基于实验的平均恢复系数模型以及考虑颗粒形状的曳力模型,能够准确地预测撞击式分离器的总分离效率和分级分离效率。颗粒分离过程中,Stokes数较大的颗粒对颗粒-壁面撞击模型比较敏感,Stokes数较小的颗粒对气固曳力模型比较敏感。     

Granular and gas-particle flows in a channel with a bidisperse particle mixture

Steady state solutions of granular and gas-particle flows in a channel with a bimodal particle mixture have been computed using kinetic theory. For granular channel flows we find granular energy equipartition breaks down with an increase in the system inelasticity and the mass ratio of particles. The effect of the particle size ratio on breakdown of energy equipartition is very small if the two particle species have the same mass. The species segregation in the solid phase is enhanced with a decrease in the system inelasticity, an increase in the average solid fraction or an increase in the size ratio, due to the competition of three diffusion forces: the thermal diffusion force, the ordinary diffusion force, and the pressure diffusion force. In addition, we find a competition mechanism exists in the equal density case (particles with equal density but different sizes) since in the equal mass case (particles with equal mass but different sizes) small particles have a higher concentration in low granular energy regions, whereas in the equal size case (particles with equal size but different masses) heavy particles have a higher concentration in low granular energy regions. These findings are in agreement with the results for granular Couette flows. For equal density particles, the segregation of large particles has a transition from the walls to the center when the restitution coefficient (e_p) decreases from 1 to 0.99. This sensitivity is reduced when the system becomes more inelastic. For a given monodisperse granular system, we show that if larger particles are mixed in the system the sensitivity of the total solid distribution to the restitution coefficient is suppressed, while if smaller particles are added in the system the situation reverses. Lastly, we extend our work to gas-particle flows in a channel where particles are fluidized by gas flowing upwards, and find that for the kinetic theory models used in the present study, the solid fraction, the species segregation and the granular energy profiles are quite similar between the granular flows and the gas-particle flows.     

Settling velocity of droplets in turbulent flows

Two parameters of particle or droplet dynamics which are of importance in describing their behaviour in turbulent pipe flows are their settling velocity and eddy diffusivity. It is usually assumed that the settling velocity in turbulent flow is equal to that in still fluid and on the basis of this assumption the eddy diffusivity is usually determined experimentally from the distribution of droplets or particles in horizontal turbulent flow. Since the settling velocity has a strong influence on the resultant value of the eddy diffusivity, the influence of turbulence on settling velocity is investigated in this work. A simple stochastic model of settling in turbulent flow is developed and it is shown that a considerable retardation of still fluid settling velocity is possible for a wide range of conditions.     

基于颗粒流运动论的喷动流化床气固流动行为三维模拟

1 INTRODUCTION Spout-fluid beds have been of increasing interest in the petrochemical, chemical and metallurgic indus-tries since spout-fluid beds can reduce some of the limitations of both spouting and fluidization by su-perimposing the two type of systems[1―4]. In recent years, spout-fluid beds have become an alternative for gas/solid contactors in coal gasification. Spout-fluid bed coal gasifiers have been adopted for APFBC-CC (advanced pressurized fluidized bed combus-tion-combined…     

Hydrodynamic modelling of dense gas-fluidised beds: comparison and validation of 3D discrete particle and continuum models

A critical comparison of a hard-sphere discrete particle model, a two-fluid model with kinetic theory closure equations and experiments performed in a pseudo-two-dimensional gas-fluidised bed is made. Bubble patterns, time-averaged particle distributions and bed expansion dynamics measured with a nonintrusive digital image analysis technique are compared to simulation results obtained at three different fluidisation velocities. For both CFD models, the simulated flow fields and granular temperature profiles are compared. The effects of grid refinement, particle-wall interaction, long-term particle contacts, particle rotation and gas-particle drag are studied. The mechanical energy balance for the suspended particles is introduced, and the energy household for both CFD models is compared. The most critical comparison between experiments and model results is given by analysis of the bed expansion dynamics. Though both models predict the right fluidisation regime and trends in bubble sizes and bed expansion, the predicted bed expansion dynamics differ significantly from the experimental results. Alternative gas-particle drag models result in significantly different bed dynamics, but the gap between model and experimental results cannot be closed. In comparison with the experimental results, the discrete particle model gives superior resemblance. The main difference between both CFD models is caused by the neglect of particle rotation in the kinetic theory closure equations embedded in the two-fluid model. Energy balance analysis demonstrates that over 80% of the total energy is dissipated by sliding friction. Introduction of an effective restitution coefficient that incorporates the additional dissipation due to frictional interactions significantly improves the agreement between both models.     

Analysis of gas-particle flow characteristics in impinging streams

A three dimensional Euler–Lagrange model for the gas-particle two-phase impinging streams (GPIS) is developed based on the direct simulation Monte Carlo (DSMC) method with consideration of particle rotation and collision. The gas-particle flow characteristics involved in GPIS as well as the effects of inlet gas velocity and particle rotation are analyzed. The results indicate that two pairs of counter-rotating gas vortices are developed at two sides of the opposite jet flows, which is able to entrain the particles and thus greatly weaken the deposition of particles. Interparticle collisions in the impingement zone produce two effects on the particle behaviors: the direct escaping of particles from impingement zone and the progressive accumulation of particles in impingement zone. Under the same inlet particle mass flow rate, the particle concentration in the impingement zone decreases with increasing inlet velocity of gas due to the increasing impinging reaction of interparticle collisions and growing entrainment of gas vortices. In addition, the rotation of particle provides an additional driving force to push the particles away from the impingement zone, leading to the higher speed of escaping particles and smaller maximum particle concentration at the center of impingement zone than those without particle rotation.     

Using the discrete element method to predict collision-scale behavior: A sensitivity analysis

Discrete element method (DEM) simulations have recently been used to investigate collision-scale measurements such as collision frequency and impact velocity distributions. These simulations are typically validated against particle velocity fields using experimental techniques such as particle image velocimetry or positron emission particle tracking. An important question that has not been addressed is whether validation of a macroscopic velocity field or solid fraction field also implies a validation of collision-scale measurements such as collision frequency. In this study, DEM measurements of solid fraction, shear rate, collision frequency, and impact velocity are made in a small region just beneath the free surface in a rotating drum. The effects of periodic drum length, particle stiffness, coefficient of restitution, and particle size are investigated. The solid fraction and shear rate do not vary with particle stiffness or coefficient of restitution over the range of values studied. However, the collision rate increases with increasing particle stiffness and coefficient of restitution. In addition, the average collision speed decreases as particles become stiffer or less elastic. The shear rate varies with particle size, but the average collision velocity remains constant. These findings indicate that validation against particle velocity and solid fraction fields does not necessarily imply validation of collision frequency and impact velocity. Indeed, the velocity and solid fraction fields were found to be relatively insensitive to a range of DEM contact stiffnesses and coefficients of restitution while the collision distributions were sensitive.     

The aspiration efficiency of a tube sampler operating in a slow moving air stream facing vertically upwards

The aspiration efficiency of a thin-walled, cylindrical aerosol sampler facing vertically upwards in a slow moving vertical air stream is numerically investigated using both a potential flow model and a viscous flow model. In order to predict the air flow around the sampler, for the potential flow model we use a boundary element method whereas for the viscous flow model we use a control volume method. The motion of the particles is then predicted by considering both the drag and gravitational forces. We have found that both numerical models produce similar predictions for the aspiration efficiency and the predictions reveal a more complex sampling behaviour when the sampler is operated in a slow moving air environment, where the air velocity is comparable with the magnitude of the particle settling velocity, than for faster moving air flows. The comparison of the numerical predictions with the only available experimental data indicates that the aspiration efficiency is in qualitative agreement but further investigations are required in order to fully reveal all the sampling characteristics in slow air flows.     

表面活性剂虫状胶束流体中颗粒沉降负尾迹模拟

颗粒在剪切稀释黏弹性表面活性剂形成的蠕虫状胶束流体中沉降时会产生负尾迹,负尾迹的形成对该种复杂流体与固体颗粒之间的相互作用具有重要影响。基于Giesekus本构方程,采用POLYFLOW软件模拟了黏弹性表面活性剂(Viscoelastic Surfactant, VES)蠕虫状胶束流体中单颗粒的沉降过程,分析了流体松弛时间和迁移因子对颗粒周围速度场及应力场的影响,重点研究了颗粒尾部速度负尾迹的产生原因及其对颗粒曳力的影响。结果表明,Giesekus本构方程能够描述VES流体的非线性剪切变稀行为和弹性导致的拉伸变形。流体弹性导致颗粒尾部产生较大的拉伸变形,剪切稀化和流体弹性的共同作用使颗粒尾部产生拉伸变形,导致负尾迹出现。表征流体弹性的De(黛博拉数)越大,流体拉伸黏度的Tr(特劳顿数)越小,负尾迹越长。负尾迹的出现使VES流体中颗粒所受曳力减小,沉降速度开始增加。模拟结果为此种流体的进一步应用提供了一定的研究基础。     

颗粒团绕流曳力系数的LBM计算

采用格子Boltzmann方法(Lattice Boltzmann method, LBM) 中的LBGK(Lattice Bhatnagar- Gross-Krook)模型和二阶精度的曲线边界条件处理方法对二维颗粒团绕流现象进行了数值模拟,并同时使用动量交换法计算了两种颗粒团构型中不同颗粒的曳力系数。结果表明：颗粒团曳力系数与颗粒聚团的构型有着密切联系,颗粒聚团的形成将导致颗粒团曳力系数大幅度减小。除颗粒团构型因素外,颗粒间距和流动Reynolds数也是导致颗粒团曳力系数发生改变的主要因素。当颗粒聚团存在时,颗粒团中不同颗粒的受力有较大差异,若忽略颗粒聚团效应,则颗粒团曳力系数的计算必然将产生偏差。     

Motion of a magnetic flow follower in two‐phase flow — application to the study of airlift reactor hydrodynamics

A low‐cost and simple magnetic particle tracer method was adapted to characterize the hydrodynamic behavior of an internal‐ and an external‐loop airlift reactor (ALR). The residence time distribution of three magnetic particles differing in diameter (5.5, 11.0 and 21.2 mm) and with a density very close to that of water was measured in individual reactor sections. The measured data were analyzed and used to determine the velocity of the liquid phase. Validation of the experimental results for liquid velocity was done by means of the data obtained by an independent reference method. Furthermore, analysis of the differences found in the settling velocity of the particle in single‐liquid and gas‐liquid phases was carried out, using a simplified 3D momentum transfer model. The model considering particle‐bubble interaction forces resulting from changes in the liquid velocity field due to bubble motion was able to predict satisfactorily the increase in the particle settling velocity in the homogeneous bubbly regime. The effective drag coefficient in two‐phase flow was found to be directly dependent on particle Reynolds number to the power of ? 2 but independent of gas flow‐rate for all particle diameters studied. Based on the experimental and theoretical investigations, the valid exact formulation of the effective buoyancy force necessary for the calculation of the correct particle settling velocity in two‐phase flow was done. In addition, recommendations concerning the use of flow‐following particles in internal‐loop ALRs for liquid velocity measurements are presented. Copyright © 2006 Society of Chemical Industry     

Numerical study of cluster formation in a gas–particle circulating fluidized bed

Simulations with two-way coupling are performed for two-dimensional gas–solid flow in a circulating fluidized bed with a total solids concentration of 3% in the riser. The motion of particles is treated by a Lagrangian approach, and particles are assumed to interact through binary, instantaneous, non-frontal, and inelastic collisions with friction. The model for the interstitial gas phase is based on the Navier–Stokes equations for two-phase flow with fluid turbulence calculated by using LES. Several porosity functions exist in the literature relating the drag force for a particle in a cloud to the drag force on an isolated particle. We have studied the influences of this porosity function, observing large differences in the local flow structure. The fluctuating gas–solid motion has been investigated showing a strong anisotropic flow behaviour, which is similar to experimental findings. The instabilities in these flows are strongly linked to the non-linear drag function due to the group effect of particles in a cloud. The collision parameters have been found to have an important influence on the cluster structures.     

Cluster-based drag coefficient model for simulating gas-solid flow in a fast-fluidized bed

Drag coefficient is of essential importance for simulation of heterogeneous gas-solid flows in fast-fluidized beds, which is greatly affected by their clustering nature. In this paper, a cluster-based drag coefficient model is developed using a hydrodynamic equivalent cluster diameter for calculating Reynolds number of the particle phase. Numerical simulation is carried out in a gas-solid fast-fluidized bed with an Eulerian-Lagrangian approach and the gaseous turbulent flow is simulated using large eddy simulation (LES). A Lagrange approach is used to predict the properties of particle phase from the equation of motion. The collisions between particles are taken into account by means of direct simulation Monte Carlo (DSMC) method. Compared with the drag coefficient model proposed by Wen and Yu, results predicted by the cluster-based drag coefficient model are in good agreement with experimental results, indicating that the cluster-based drag coefficient model is suitable to describe various statuses in fast-fluidized beds.     

THE EFFECT OF STATIC CHARGES ON DRAG REDUCTION IN DILUTE GAS-PARTICLE SUSPENSION FLOW

This is an experimental study on the drag reduction phenomenon in a gas-particle suspension two-phase flow system. The purpose of the research was to examine the effect of static charges carried by the conveyed particles on the drag reduction characteristics of a flowing suspension in horizontal and vertical tubes. A specifically designed experimental loop was used to perform simultaneous pressure drop and particle charge measurements. The charge on the particles flowing in the system was artificially increased using high voltage rod-in-cylinder Corona chargers mounted in the closed loop. The particle charge was measured in a Farraday cage sampling device and the effective current generated by particle impact on a specially designed spherical nickel probe mounted in the two-phase flow was also recorded.

It was found during the measurements that the pressure drop in the vertical test section was not significantly influenced by the charge build-up on the flowing particles, while the results obtained in the horizontal test section showed a significant effect on the drag due to charged particles. The electric charge on the particles was shown to have a negative contribution on the drag reduction: the higher the value of the effective charge, the higher the drag increase, i.e., the smaller the drag reduction.