Particle size distribution in CPFD modeling of gas–solid flows in a CFB riser

A computational particle fluid dynamics (CPFD) numerical method to model gas–solid flows in a circulating fluidized bed (CFB) riser was used to assess the effects of particle size distribution (PSD) on solids distribution and flow. We investigated a binary PSD and a polydisperse PSD case. Our simulations were compared with measured solids concentrations and velocity profiles from experiments, as well as with a published Eulerian-Eulerian simulation. Overall flow patterns were similar for both simulation cases, as confirmed by experimental measurements. However, our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser, as seen in other numerical studies. Above this bottom region, distributions of particle volume fraction and particle vertical velocity were consistent with our experiments, and the simulated average particle diameter decreased as a power function with riser height. Interactions between particles and walls also were successfully modeled, with accurate predictions for the lateral profiles of particle vertical velocity. It was easy to implement PSD into the CPFD numerical model, and it required fewer computational resources compared with other models, especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Particle size distribution in CPFD modeling of gas-solid flows in a CFB riser

A computational particle fluid dynamics(CPFD) numerical method to model gas-solid flows in a circulating fluidized bed(CFB) riser was used to assess the effects of particle size distribution(PSD) on solids distribution and flow.We investigated a binary PSD and a polydisperse PSD case.Our simulations were compared with measured solids concentrations and velocity profiles from experiments,as well as with a published Eulerian-Eulerian simulation.Overall flow patterns were similar for both simulation cases,as confirmed by experimental measurements.However,our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser,as seen in other numerical studies.Above this bottom region,distributions of particle volume fraction and particle vertical velocity were consistent with our experiments,and the simulated average particle diameter decreased as a power function with riser height.Interactions between particles and walls also were successfully modeled,with accurate predictions for the lateral profiles of particle vertical velocity.It was easy to implement PSD into the CPFD numerical model,and it required fewer computational resources compared with other models,especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Circulating turbulent fluidization—A new fluidization regime or just a transitional phenomenon

While circulating fluidized bed （CFB） reactor has many advantages over the more conventional turbulent fluidized bed （TFB） reactor, it does at least have one significant shortcoming-the rather dilute solids volume concentration in CFB reactor gives rise to less ideal reaction intensity. On the other hand, while having higher reaction intensity, TFB reactor has one fatal drawback of particle back-mixing, making it not suitable for certain reactions such as catalytic reaction where the catalyst requires frequent regeneration. This paper describes some key issues in the development of a circulating turbulent fluidized bed （CTFB） reactor that combines the advantages of both TFB and CFB, that is, to have the high reaction intensity as in TFB but and also to have a suppressed solids back-mixing as in CFB due to a continuous net upflow of solids flux through the bed. Experimental results show enough evidence to suggest that a new fluidization regime is formed, the characteristics of which appears to be distinct from those observed in a regular TFB and from those in either the bottom or the upper sections of regular CFB and/or high-density CFB （HDCFB）. Fundamentally, the difference is that particle-particle interaction （collision） dominates the motion of particles in CTFB and TFB, while gas-particle interaction （drag force） is the key element that determines the two phase flow in CFB including HDCFB.     

CPFD modeling of hydrodynamics,combustion and NOx emissions in an industrial CFB boiler

The ultra-low NOx emission requirement (50 mg/m³) brings great challenge to CFB boilers in China. To further tap the NOx abatement potential, full understanding the fundamentals behind CFB boilers is needed. To achieve this, a comprehensive CPFD model is established and verified; gas-solid flow, combustion, and NOx emission behavior in an industrial CFB boiler are elaborated; influences of primary air volume and coal particle size on furnace performance are evaluated. Simulation results indicate that there exists a typical core-annular flow structure in the boiler furnace. Furnace temperature is highest in the bottom dense-phase zone (about 950 °C) and decreases gradually along the furnace height. Oxygen-deficient combustion results in high CO concentration and strong reducing atmosphere in the lower furnace. NOx concentration gradually increases in the bottom furnace, reaches maximum at the elevation of secondary air inlet, and then decreases slightly in the upper furnace. Appropriate decreasing the primary air volume and coal particle size would increase the CO concentration and intensify the in-furnace reducing atmosphere, which favors for NOx reduction and low NOx emission from CFB boilers.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. 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 wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.     

Simulation of mass balance behavior in a large-scale circulating fluidized bed reactor

We determine using a compound model the influence of the mass of granular matter on the behavior of a supercritical circulating fluidized bed (CFB) reactor. Population balance enables a stationary-regime modeling of the mass flow of granular matter inside a CFB unit in a large-scale. The simulation includes some important dynamic processes of gas-particle flows in fluidized bed such as attrition, fragmentation, elutriation, and fuel combustion. Numerical calculations with full boiler loading were performed of operational parameters such as furnace temperature, furnace pressure, feeding materials mass flows, and excess air ratio. Furthermore, three bed inventory masses were adopted as experimental variables in the simulation model of mass balance. This approach enables a sensitivity study of mass flows of granular matter inside a CFB facility. Some computational results from this population balance model obtained for a supercritical CFB reactor are presented that show consistency with the operational data for large-scale CFB units.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics （CFD） modeling of the gas-solids two-phase flow in a circulating fluidized bed （CFB） riser is carried out. 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 wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory （Johnson ＆ Jackson, 1987） with specularity coefficient and particle-wall restitution coefficient.The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operatin~ conditions of the CFB riser.     

An image processing method for feature extraction and dynamic tracking of particle clusters in CFBs

A new image processing method based on the high-speed camera is proposed to identify, locate, and track clusters. The instantaneous characteristic parameters of particle clusters in the riser of the circulating fluidized bed (CFB) can be acquired, such as solids holdup, vertical velocity, lateral displacement, aspect ratio and near-circularity. Experiments were carried out with glass bead particles, river sand particles and FCC particles. The time series of images of gas–solid flow in a CFB riser with a 100 mm × 25 mm cross-section and 3.2 m in length were obtained using high-speed cameras. The k-means++ clustering algorithm is utilized to identify the clusters, centroid is applied to locate the clusters, and the cross-correlation algorithm is employed to track the specific clusters and number them to get the instantaneous characteristic parameters. The results illustrate that the shapes of clusters in the center area are closest to circle, moving upwards at a uniform speed, while the clusters in the side-wall area are mostly elongated or long chain-like, moving slowly downwards. In the transition area, the clusters are more complex, moving upwards at a constant speed, and having large lateral displacement. The results show that the image processing method used in this study is successful in acquiring the dynamic and structural parameters of the clusters simultaneously.     

复杂气固两相系统的微观结构

流化床中的气固两相流动是一个高度复杂的非线性混沌系统。本文利用激光粒子动态分析仪（ＰＤＡ）得到的循环流化床中颗粒脉动速度信号,采用ＦＦＴ分析了脉动信号的宽频谱特征,在此基础上应用小波法分析了脉动信号的动态特征,得到了颗粒脉动速度的微观结构,指出颗粒脉动速度的非线性特性是流化床具有混沌特性的根源,且在不同的尺度上颗粒脉动速度表现出各向异性的特征。     

Investigation on lateral dispersion characteristics of fuel particles in dense phase zone of a large-scale circulating fluidized bed boiler under combustion conditions

The dispersion characteristics of fuel particles in the dense phase zone in circulating fluidized bed (CFB) boilers have an important influence on bed temperature distribution and pollutant emissions. However, previous research in literature was mostly on small-scale apparatus, whose results could not be applied directly to large-scale CFB with multiple dispersion sources. To help solve this problem, we proposed a novel method to estimate the lateral dispersion coefficient (D_x) of fuel particles under partial coal cut-off condition in a 350 MW supercritical CFB boiler based on combustion and dispersion models. Meanwhile, we carried out experiments to obtain the D_x in the range of 0.1218–0.1406 m²/s. Numerical simulations were performed and the influence of operating conditions and furnace structure on fuel dispersion characteristics was investigated, the simulation value of D_x was validated against experimental data. Results revealed that the distribution of bed temperature caused by the fuel dispersion was mainly formed by char combustion. Because of the presence of intermediate water-cooled partition wall, the mixing and dispersion of fuel and bed material particles between the left and right sides of the furnace were hindered, increasing the non-uniformity of the bed temperature near furnace front wall.     

气固两相混合层二维涡配对的数值研究

采用双向耦合模型中的速度耦合模型,数值模拟了气固两相混合层中涡的配对、合并过程,文中采用拟谱方法对流场进行了直接数值模拟,用Lagrange模型跟踪固粒,结果发现,在双向耦合过程中Stokes数仍然是重要的参数,但不是唯一影响流场的参数,流场的发展还与固粒的体积浓度、固粒的相对密度以及固粒大小等因素有关。固粒抑制流场的发展,阻碍涡的配对、合并,加快了涡量的扩散;小St数的固粒仍然跟随流体运动,大St数的固粒趋向于低涡量区的趋势减弱。     

Particle-scale and sub-grid drag models coupled CFD for simulating the CO methanation in a CFB riser

The diffusion and chemical reactions inside the catalyst particles and the heterogeneous flow structure in the computational cells are key factors to affect the accuracy of the coarse-grid simulation in circulating fluidized bed (CFB) methanation reactors. In this work, a particle-scale model is developed to calculate the effective reaction rate considering the transient diffusion and chemical reactions in the particle scale, i.e., the scale of the single catalyst particle. A modified sub-grid drag model is proposed to consider the effects of the meso-scale and chemical reactions on the heterogeneous gas-solid interaction, where the meso-scale is between the single particle and the whole reactor and featured with the particle cluster. Subsequently, a coupled model is developed by integrating the particle-scale and modified sub-grid drag models into CFD. Moreover, the coupled model is validated to achieve accurate predictions on the CO methanation process in a CFB riser. Notably, the coupled model can be performed with a coarse grid (∼58 times particle diameter) and a large time step (0.005 s) to accelerate the simulation. By simply changing the reaction kinetics, different gas-solid catalytic reaction systems can be simulated by using the coupled model.     

Particle/Turbulence Interactions, Mass Transfer and Gas/Solid Chemistry in a CFBC Riser

Gas/solid chemistry in the upper, dilute region of a circulating fluidised bed combustor (CFBC) riser is addressed. The limitations of turbulent mixing are illustrated by the example of the heterogeneous NO/CO/char reaction, relevant in CFB combustion of nitrogen-containing solid fuels. The mass transfer of the gaseous reactants to the char surface is determined, and how the conversion is influenced by the degree of mixing of the multiphase system by turbulent dispersion. Particle/turbulence interactions are taken into account by a (Lagrangian) frequency spectrum of velocity fluctuations, which determine the dispersion of the char particles described here with the Tchen–Hinze model. Chars from solid fuels characterised by fuel ratio (FR) ranging from 0.1 (wood) via 0.5 (peat) and 1 (coal A) to 2 (coal B) were considered. The effective rate of the NO/CO/char reaction is determined as a function of the size and type of the char particle, temperature, particle concentration, reactor dimensions and fluidization velocity, at atmospheric pressure. It was found that for this case the effective gas/solid conversion rate in the upper, dilute region of the riser is much lower than the gas/solid chemistry, mainly due to mass transfer limitations for char particles with sizes of typically 300 μm. The concentration of NO at the char particle surface is only a few % of that in the bulk gas phase. Strong influences were found for particle size and temperature, whilst the fluidisation velocity and the reactor size have only a small influence. It is concluded that for a typical CFB riser, for particles larger than approx. 20 μm, mass transfer has a stronger influence on the heterogeneous NO/CO/char reduction mechanism than the ``unmixedness' due to particle eddy dispersion limitations. It is recommended that this or a similar approach to turbulent dispersive mixing is implemented in CFD codes when these are used for boiler and furnace calculations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dual characters of interaction between particles and fluid

The unique characteristics of gas-solids two-phase flow and fluidization in terms of the flow structures and the apparent behavior of particles and fluid-particle interactions are closely linked to physical properties of the particles, operating conditions and bed configurations. Fluidized beds behave quite differently when solid properties, gas velocities or vessel geometries are varied. An understanding of hydrodynamic changes and how they, in turn, influence the transfer and reaction characteristics of chemical and thermal operations by variations in gas-solid contact, residence time, solid circulation and mixing and gas distribution is very important for the proper design and scale-up of fluidized bed reactors. In this paper, rather than attempting a comprehensive survey, we concentrate on examining some important positive and negative impacts of particle sizes, bubbles, clusters and column walls on the physical and chemical aspects of chemical reactor performance from the engineering application point of view with the aim of forming an adequate concept for guiding the design of multiphase fluidized bed chemical reactors.One unique phenomenon associated with particle size is that fluidized bed behavior does not always vary monotonically with changing the average particle size. Different behaviors of particles with difference sizes can be well understood by analyzing the relationship between particle size and various forces. For both fine and coarse particles, too narrow a distribution is generally not favorable for smooth fluidization. A too wide size distribution, on the other hand, may lead to particle segregation and high particle elutriation. Good fluidization performance can be established with a proper size distribution in which inter-particle cohesive forces are reduced by the lubricating effect of fine particles on coarse particles for Type A, B and D particles or by the spacing effect of coarse particles or aggregates for Type C powders.Much emphasis has been paid to the negative impacts of bubbles, such as gas bypassing through bubbles, poor bubble-to-dense phase heat & mass transfer, bubble-induced large pressure fluctuations, process instabilities, catalyst attrition and equipment erosion, and high entrainment of particles induced by erupting bubbles at the bed surface. However, it should be noted that bubble motion and gas circulation through bubbles, together with the motion of particles in bubble wakes and clouds, contribute to good gas and solids mixing. The formation of clusters can be attributed to the movement of trailing particles into the low-pressure wake region of leading particles or clusters. On one hand, the existence of down-flowing clusters induces strong solid back-mixing and non-uniform radial distributions of particle velocities and holdups, which is undesirable for chemical reactions. On the other hand, the formation of clusters creates high solids holdups in the riser by inducing internal solids circulations, which are usually beneficial for increasing concentrations of solid catalysts or solid reactants.Wall effects have widely been blamed for complicating the scale-up and design of fluidized-bed reactors. The decrease in wall friction with increasing the column diameter can significantly change the flow patterns and other important characteristics even under identical operating conditions with the same gas and particles. However, internals, which can be considered as a special wall, have been used to improve the fluidized bed reactor performance.Generally, desirable and undesirable dual characteristics of interaction between particles and fluid are one of the important natures of multiphase flow. It is shown that there exists a critical balance between those positive and negative impacts. Good fluidization quality can always be achieved with a proper choice of right combinations of particle size and size distribution, bubble size and wall design to alleviate the negative impacts.     

Direct numerical simulation of particle dispersion in a turbulent jet considering inter-particle collisions

To investigate the behaviour of inter-particle collision and its effects on particle dispersion, direct numerical simulation of a three-dimensional two-phase turbulent jet was conducted. The finite volume method and the fractional-step projection algorithm were used to solve the governing equations of the gas phase fluid and the Lagrangian method was applied to trace the particles. The deterministic hard-sphere model was used to describe the inter-particle collision. In order to allow an analysis of inter-particle collisions independent of the effect of particles on the flow, two-way coupling was neglected. The inter-particle collision occurs frequently in the local regions with higher particle concentration of the flow field. Under the influence of the local accumulation and the turbulent transport effects, the variation of the average inter-particle collision number with the Stokes number takes on a complex non-linear relationship. The particle distribution is more uniform as a result of inter-particle collisions, and the lateral and the spanwise dispersion of the particles considering inter-particle collision also increase. Furthermore, for the case of particles with the Rosin–Rammler distribution (the medial particle size is set d₅₀ = 36.7 μm), the collision number is significantly larger than that of the particles at the Stokes number of 10, and their effects on calculated results are also more significant.     

Direct numerical simulation of particle behaviour in the wall region of turbulent flows in horizontal channels

The motion of small particles in the wall region of turbulent channel flows has been investigated using direct numerical simulation. It is assumed that the particle concentration is low enough to allow the use of one-way coupling in the calculations, i.e. the fluid moves the particles but there is no feedback from the particles on the fluid motion. The velocity of the fluid is calculated by using a pseudospectral, direct solution of the Navier-Stokes equations. The calculations indicate that particles tend to segregate into the low-speed regions of the fluid motion near the wall. The segregation tendency depends on the time constant of the particle made non-dimensional with the wall shear velocity and kinematic viscosity. For very small and very large time constants, the particles are distributed more uniformly. For intermediate time constants (of the order 3), the segregation into the low-speed fluid regions is the highest. The finding that segregation occurs for a range of particle time constants is supported by experimental results. The findings regarding the more uniform distributions, however, still remain to be verified against experimental data which is not yet available. For horizontal channel flows, it is also found that particles are resuspended by ejections (of portions of the low-speed streaks) from the wall and are, therefore, primarily associated with low-speed fluid. The smaller particles are flung further upwards and, as they fall back towards the wall, they tend to be accelerated close to the fluid velocity. The larger particles have greater inertia and, consequently, accelerate to lower velocities giving higher relative velocities. This velocity difference, as a function of wall-normal distance, follows the same trend as in experiments but is always somewhat smaller in the calculations. This appears to be due to the Reynolds number for the numerical simulation being smaller than that in the experiment. It is concluded that the average particle velocity depends not only on the wall variables for scaling, but also on outer variables associated with the mean fluid velocity and fluid depth in the channel. This is because fluid depth in combination with the wall shear velocity determines how much time a particle, of a given size and density, spends in the outer flow and, hence, how close it gets to the local fluid velocity.     

Hydrodynamics in a new liquid–solid circulating conventional fluidized bed

A new type of liquid–solid fluidized bed, named circulating conventional fluidized bed (CCFB) which operates below particle terminal velocity was proposed and experimentally studied. The hydrodynamic behavior was systematically studied in a liquid–solid CCFB of 0.032 m I.D. and 4.5 m in height with five different types of particles. Liquid–solid fluidization with external particle circulation was experimentally realized below the particle terminal velocity. The axial distribution of local solids holdup was obtained and found to be fairly uniform in a wide range of liquid velocities and solids circulation rates. The average solids holdup is found to be significantly increased compared with conventional fluidization at similar conditions. The effect of particle properties and operating conditions on bed behavior was investigated as well. Results show that particles with higher terminal velocity have higher average solids holdup.     

Modeling of particle transport and combustion phenomena in a large-scale circulating fluidized bed boiler using a hybrid Euler-Lagrange approach

The constantly developing fiuidized combustion technology has become competitive with a conventional pulverized coal （PC） combustion. Circulating fluidized bed （CFB） boilers can be a good alternative to PC boilers due to their robustness and lower sensitivity to the fuel quality. However, appropriate engineering tools that can be used to model and optimize the construction and operating parameters of a CFB boiler still require development. This paper presents the application of a relatively novel hybrid Euler-Lagrange approach to model the dense gas-solid flow combined with a combustion process in a large-scale indus- trial CFB boiler. In this work, this complex flow has been resolved by applying the ANSYS FLUENT 14.0 commercial computational fluid dynamics （CFD） code. To accurately resolve the multiphase flow, the original CFD code has been extended by additional user-defined functions. These functions were used to control the boiler mass load, particle recirculation process （simplified boiler geometry）, and interphase hydrodynamic properties. This work was split into two parts. In the first part, which is referred to as pseudo combustion, the combustion process was not directly simulated. Instead, the effect of the chemi- cal reactions was simulated by modifying the density of the continuous phase so that it corresponded to the mean temperature and composition of the flue gases, In this stage, the particle transport was simu- lated using the standard Euler-Euler and novel hybrid Euler-Lagrange approaches, The obtained results were compared against measured data, and both models were compared to each other. In the second part, the numerical model was enhanced by including the chemistry and physics of combustion. To the best of the authors＇ knowledge, the use of the hybrid Euler-Lagrange approach to model combustion is a new engineering application of this model, In this work, the combustion process was modeled for air-fuel combustion. The simulation results were compared with experimental data.     

Flow regime diagrams for gas-solid fluidization and upward transport

Flow regime maps are presented for gas-solids fluidized beds and gas-solids upward transport lines. For conventional gas solids fluidization, the flow regimes include the fixed bed, bubbling fluidization, slugging fluidization and turbulent fluidization. For gas solids vertical transport operation, solids flux must be incorporated in the flow regime diagrams. The flow regimes then include dilute-phase transport, fast fluidization or turbulent flow, slug/bubbly flow, bubble-free dense-phase flow and packed bed flow. In practical circulating fluidized beds and transport risers, operation below the fast fluidization regime is commonly impossible due to equipment limitations. Practical flow regime maps are proposed with the flow regimes, including homogeneous dilute-phase flow, core-annular dilute-phase flow (where there are appreciable lateral gradients but small axial gradients) and fast fluidization (where there are both lateral and axial gradients). The boundary between fast fluidization and dilute-phase pneumatic transport is set by the type A choking velocity, at which the uniform suspension collapses and particles start to accumulate in the bottom region of the transport line, while the mechanism of transition from fast fluidization to dense-phase flow depends on the column and particle diameters.     

Constraint strength and axial/radial particle velocity profiles for an integrated riser outlet

To study axial/radial profiles of particle velocity in the affected region of an integrated riser outlet,a cold model was developed for the integrated riser reactor combining the gas-solid distributor with the fluidized bed.Constraints,related to the gas-solid distributor and the upper fluidized bed,imposed on the particle flow in the riser outlet region,were investigated experimentally.The experimental results showed that with increasing superficial gas velocity,these constraints have strong influences on particle flow behavior,the particle circulation flux in the riser,and the height of the static bed material of the upper fluidized bed.When the constraints have greater prominence,the axial profile of the cross-sectionally averaged particle velocity in the outlet region initially increases and then decreases,the rate of decrease being proportional to the constraint strength.Along the radial direction of the outlet section,the region where the local particle velocity profile tends to decrease appears near the dimensionless radius r/R = 0.30 initially and then,with increasing constraint strength,gradually extends to the whole section from the inner wall.Based on the experimental data,an empirical model describing the constraint strength was established.The average relative error of the model is within 7.69%.