Evaluation of wall boundary condition parameters for gas–solids fluidized bed simulations

Wall boundary conditions for the solids phase have significant effects on numerical predictions of various gas–solids fluidized beds. Several models for the granular flow wall boundary condition are available in the open literature for numerical modeling of gas–solids flow. A model for specularity coefficient used in Johnson and Jackson boundary conditions by Li and Benyahia (Li and Benyahia, AIChE J. 2012;58:2058–2068) is implemented in the open‐source CFD code‐MFIX. The variable specularity coefficient model provides a physical way to calculate the specularity coefficient needed by the partial‐slip boundary conditions for the solids phase. Through a series of two‐dimensional numerical simulations of bubbling fluidized bed and circulating fluidized bed riser, the model predicts qualitatively consistent trends to the previous studies. Furthermore, a quantitative comparison is conducted between numerical results of variable and constant specularity coefficients to investigate the effect of spatial and temporal variations in specularity coefficient. Published 2013 American Institute of Chemical Engineers AIChE J, 59: 3624–3632, 2013     

Microstructural aspects of the flow behaviour in a liquid—solids circulating fluidized bed

The local instantaneous and time‐average suspension densities were determined in a 76 mm diameter by 3 m tall liquid‐solids circulating fluidized bed riser using a fibre‐optic probe. Attempts were made to qualify the microflow structure through statistical analysis of the local bed voidage fluctuations obtained under different operating conditions for the first time. The results show that local microflow structure is uniform in the axial direction but non‐uniform in the radial direction with more flow fluctuation near the wall than in the core of the column for a given axial position. The standard deviation and intermittency index tend to increase with increasing solids circulating rates. Comparing with the gas—solids CFB, the liquid—solids CFB shows much more homogeneous flow structure in both the axial and radial microscopic flow behaviours. The microflow behaviours in the conventional liquid—solids fluidization, liquid—solids circulating fluidization and dilute‐phase liquid transport regimes are also characterized by examining the probability distribution and the intermittency index of the solids holdup.     

Annular wall layer thickness in circulating fluidized bed risers

The thickness of downward-flowing annular wall layers in circulating fluidized bed risers has been determined in the literature based on measured radial profiles of both local particle velocity and solids flux. The thickness of the wall layer is shown to be larger based on solids flux profiles than when based on particle velocity profiles, because fluctuations in local instantaneous particle velocity are correlated with fluctuations in local solids concentration. A new correlation is developed to predict the time-average thickness of the downflowing particle streamer layer based on solids flux measurements as a function of the cross-sectional average voidage. It successfully accounts for the variation of the wall layer thickness with axial location and solids circulation rate.     

Mixing of a lateral gas stream in a two-dimensional riser of a circulating fluidized bed

A cold model of a circulating fluidized bed having a two-dimensional riser, with a 12 × 120 mm section and a 6.4 m height, was equipped with a device to inject a lateral gas stream along the riser. The apparatus was operated under conditions ranging from those characteristic of combustors to those of gas-conversion processes. Flow structures in the interaction region between the rising gas—solids suspension and the lateral gas stream were studied by means of a motion analysis system. Three main configurations were identified. The ratio between the momentum of the lateral gas stream and that of the rising suspension was found to be the parameter able to discriminate among the three configurations. A satisfactory agreement was found with mixing data obtained using a cylindrical riser.     

Flow structure and thickness of annular downflow layer in a circulating fluidized bed riser

Flow structures were determined in a circulating fluidized bed (CFB) riser (0.203 m i.d.×5.9 m high) of FCC particles (d_p=70 μm, ρ_s=1700 kg/m³). A momentum probe was used to measure radial momentum flux profiles at several levels and to distinguish between upward and downward flow regions. Time-mean dynamic pressure (ΔP_m) decreases towards the wall in the range U_g=5-8 m/s, G_s=10-340 kg/m² s. The thickness of the annular downflow layer based on ΔP_m=0 reaches a maximum with increasing height. The annular downflow layer disappears locally with increasing solids mass flux (G_s) at a constant gas velocity, with achievement of the dense suspension upflow (DSU) regime. A new correlation is developed to predict the time-mean thickness of solids down-flowing layer based on solids mass flux and momentum flux. It successfully accounts for the variation of the annular layer thickness with height and G_s, and covers a wide G_s range right up to near the onset of the DSU regime.     

Cocurrent downflow circulating fluidized bed (downer) reactors — A state of the art review

In the last decade, cocurrent downflow circulating fluidized bed reactors, called “downer” reactors, have been proposed as an alternative to upflow circulating fluidized bed, or “riser”, reactors. In this paper, published results on downer studies are summarized and future directions of research are recommended. Downer reactors are shown to have several distinct advantages over upflow circulating fluidized bed reactors and can potentially be used in many industrial processes, mainly due to a more uniform gas and solids flow structure compared with risers.     

Modelling and simulation of a gas–solids dispersion flow in a high-flux circulating fluidized bed (HFCFB) riser

Radial solids velocity profiles were computed on seven axial levels in the riser of a high-flux circulating fluidized bed (HFCFB) using a two-phase 3-D computational fluid dynamics model. The computed solids velocities were compared with experimental data on a riser with an internal diameter of 76 mm and a height of 10 m, at a high solids flux of 300 kg m⁻² s⁻¹ and a superficial velocity of 8 m s⁻¹. Several hundreds of experimental and numerical studies on CFBs have been carried out at low fluxes of less than 200 kg m⁻² s⁻¹, whereas only a few limited useful studies have dealt with high solids flux. The k– two-phase turbulence model was used to describe the gas–solids flow in an HFCFB. The model predicts a core–annulus flow in the dilute and developed flow regions similar to that found experimentally, but in the region of highest solids concentration it is somewhat overpredicted at the level close to the inlet.     

Phase distributions in a gas–liquid–solid circulating fluidized bed riser

The distributions of the three phases in gas–liquid–solid circulating fluidized beds (GLSCFB) were studied using a novel measurement technique that combines electrical resistance tomography (ERT) and optical fibre probe. The introduction of gas into a liquid–solid circulating fluidized bed (LSCFB), thus forming a GLSCFB, caused the increase of solids holdup due to the significantly decreased available buoyancy with the lower density of the gas, even with a somewhat increased liquid velocity due to the decreased liquid holdup giving space for the gas holdup. The gas passed through the riser in the form of bubbles, which tended to flow more through the central region of the riser, leading to more radial non‐uniformity in radial holdup of the phases. The gas velocity has the most significant effect on the gas phase holdup. While the gas velocity also has an obvious effect to the solids holdups, the liquid flow rate had a much more considerable effect on the phase holdups. The solids circulation rate also had a significant effect on the phase holdups, with increasing solids circulation rate causing much more increased solids holdup in the central region than close to the wall. A correlation was developed for the relative radial distributions of solids holdup in GLSCFB, as such radial profiles were found similar over a wide range of operating conditions, like those in a typical gas–solid circulating fluidized beds (GSCFB). Finally, the axial solids profiles in a GLSCFB was found to be much closer to those in an LSCFB which are very uniform, than those found in a GSCFB which are less uniform and sometime having a S shape. Water was used as the continuous and conductive phase, air was the gas phase and glass bead and lava rock particles were used as the solid and non‐conductive phase.     

Influence of ring-type internals on axial pressure distribution in circulating fluidized bed

In order to improve non-uniform radial and axial flow structure of a circulating fluidized bed, the influence of ring-type internals on the axial pressure distribution and gas—solids flow structure in a riser of 7.6 cm in diameter and 3 m in height was investigated experimentally. Four different opening areas, 70%, 80%, 90% and 95%, were used and the superficial gas velocity and solids circulation rate were in the ranges of 5 to 10 m/s and 20 to 233 kg/m²·s, respectively. With the presence of internals, the axial pressure gradient distribution shows the formation of a zigzag type profile instead of the regular exponential or S-shape profile and the bottom acceleration region is shortened. The opening ratio of the rings plays an important role in affecting the flow structure. The optimal opening ratio is tightly related to the operating conditions. In the circulation fluidized bed used in this study, 90% open area was found to be most suitable for obtaining a more uniform gas—solids flow structure.     

Hydrodynamics and reactor performance evaluation of a high flux gas‐solids circulating fluidized bed downer: Experimental study

Reactor performance of a high flux circulating fluidized bed (CFB) downer is studied under superficial gas velocities of 3–7 m/s with solids circulation rate up to 300 kg/m²s using ozone decomposition reaction. Results show that the reactant conversion in the downer is closely related to the hydrodynamics, with solids holdup being the most influential parameter on ozone decomposition. High degree of conversion is achieved at the downer entrance region due to strong gas‐solids interaction as well as higher solids holdup and reactant concentration. Ozone conversion increases with the increase of solids circulation rate and/or the decrease of superficial gas velocity. Overall conversion in the CFB downer is less than but very close to that in an ideal plug flow reactor indicating a good reactor performance in the downer because of the nearly “ideal” hydrodynamics in downer reactors. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3412–3423, 2014     

Modelling of circulating fluidized bed combustion of a char

The combustion of a char in the 41 mm ID riser of a laboratory circulating fluidized bed combustor has been investigated at different air excesses and rates of solids (char and sand) circulating in the loop. Riser performance was characterized by an axial oxygen concentration profile as well as by the overall carbon content and particle size distribution. The proposed model accounts for carbon surface reaction, intraparticle and external diffusion, and attrition. External diffusion effects were relevant in the riser dense region where char was potentially entrapped in large clusters of inert solids. Experimental data and results of the model calculations are in satisfactory agreement.     

The impacts of standard wall function and drag model on the turbulent modelling of gas–particle flow in a circulating fluidised bed riser

A transient turbulence model was applied to simulate the gas–particle system in a circulating fluidised bed riser. The k–epsilon turbulent equations coupled with the fluctuating energy equation were used to simulate the gas–particle system in a riser. The simulation results were validated by the experimental data of a CFB system. A grid study was implemented to examine the impact of grid discretisation. A comparison between the conventional drag models and the EMMS model was also conducted. Other factors, like the restitution coefficient particle to particle, was also found to have a significant impact on the turbulence model. © 2013 Canadian Society for Chemical Engineering     

Axial hydrodynamic studies in a gas-liquid-solid circulating fluidized bed riser

Axial distribution of phase holdups was studied in the riser of a gas-liquid-solid circulating fluidized bed (GLSCFB). The effects of gas and liquid superficial velocities as well as solids circulation rate on radial distribution of phase holdups at different axial locations were investigated. Electrical resistance tomography (ERT) and optical fiber probe were employed online in the experiments for a precise determination of phase holdups. An empirical model was developed for the determination of gas bubbles in analysis of data obtained by fiber optic sensor. Gas holdup was higher at the central region of the riser and increased axially due to coalescence of small bubbles and decrease of hydrostatic pressure at higher levels in the riser. This led to an increase in solids holdup in regions close to the wall which was slightly higher than the solids holdup at the wall. Both solids and liquid holdups were lower in the central region and increased radially towards the wall. Gas holdup decreased with increasing solids circulation rate but opposite trend was observed for solids holdup. Solids circulation rate had negligible effect on liquid holdup at lower axial locations compared to top of the riser. Cross-sectional average of solids, gas and liquid holdups did not change significantly at higher liquid superficial velocities.     

Effect of temperature on gas bypassing and solids circulation rate in an internally circulating fluidized bed

The effect of bed temperature in the annulus region on gas bypassing from annulus to riser and solids circulation rate has been determined in an internally circulating fluidized bed. The bed voidage in the annulus region increases with increasing temperature. The average fractions of gas bypassing of annulus gas are 59.4% and 87% at room temperature and high temperature (400°C to 800°C), respectively. Solids circulation rate increases with increasing annulus gas velocity, but is nearly independent of riser velocity. Solids circulation rate increases with increasing temperature.     

Exit effect on hydrodynamics of the internal circulating fluidized bed riser

This is the first time an extensive investigation has been carried out regarding the effects of riser exit geometry on pressure drop and solid behaviour inside the Internal Circulating Fluidized Bed (ICFB) riser, using different riser exit geometries at several operating conditions.The Radioactive Particle-Tracking (RPT) technique was used for solid concentration measurements and solid residence time distribution at the exit zone. Experiments were conducted using Geldart B particles, in the gas superficial velocity range of 4 to 10 m/s. Axial solid hold-up, solid residence time distribution in the exit zone, and the reflux ratio factor k_m, (defined earlier by [E.H. Van der Meer, R.B. Thorpe, J.F. Davidson, Flow patterns in the square cross-section riser of a circulating fluidized bed and the effect of riser exit design, Chem. Eng. Sc. 55 (19) (2000) 4079-4099]), were the main criteria used to investigate the impact of gas-solid separator devices implemented at the ICFB riser exit.Solid residence time distribution results and axial solid hold-up profiles provided clear evidence that the separator device at the riser exit strongly influences the hydrodynamic structure of the ICFB riser. The V-shaped riser exit geometry was found to be the optimum of all the configurations studied.     

循环流化床提升管压力瞬时波动的功率谱分析

在φ0.4 m×9.1 m循环流化床提升管中,采集了分布板以上不同轴向高度的压力瞬时波动信号,并采用功率谱进行分析。分析发现,压力瞬时波动功率谱谱图中存在一个振幅最大点,即为主频。增大固体颗粒循环速率或减小提升管操作气速,主频对应振幅增大,而功率谱主频减小。相同操作条件下,随轴向位置的增高,压力瞬时波动的主频基本不变,而主频对应振幅减小。     

Characterization on hydrodynamic behavior in liquid‐containing gas‐solid fluidized bed reactor

In many industrial processes involving gas–solid fluidized bed rectors, the addition of a liquid phase significantly alters the hydrodynamics. To fully characterize the hydrodynamics in the fluidized bed, pressure and acoustic measuring techniques were applied to study the behavior of gas bubbles and particles. A camera was used to take pictures to verify the pressure and acoustic results. During the liquid‐addition process, the pressure technique captured the bubble size variation and bubble motion while the acoustic technique reflected particle motion and particle size growth. Hurst and V‐statistics analyses of acoustic emission were used for the first time to detect periodic behavior during the injection process. The new break formation and change trend of V_max were used as the criteria to judge occurrence of abnormal fluidization states, such as agglomeration and gas channeling formation. These measurement techniques are beneficial in the elimination of adverse effects caused by the addition of liquid. © 2012 American Institute of Chemical Engineers AIChE J, 59: 1056–1065, 2013     

Hydrodynamic modelling of the axial density profile in the riser of a low‐density circulating fluidized bed

The axial density profile is an important characteristic of the CFB riser and a key parameter in the CFB design. A simple, but reliable model is needed to predict the density profiles. Elutriation‐based models treat the dense phase at the bottom of the riser as a dense bubbling bed whereas the dilute phase higher up can be looked upon as the entrainment zone above the dense bed. The elutriation model, as originally presented by Rhodes et al. (1986) and based on Wen et al. (1982) is extensively studied and modified. In spite of the modifications, the use of entrainment models has certain clear limitations due to a wide range of predictions as evident from Table 1. Elutriation rates are calculated based on the hydrodynamic phenomena in the dense bed and a fitting procedure for the entrainment decay constant (σ) was performed.     

Experimental and numerical research for fluidization behaviors in a gas–solid acoustic fluidized bed

The effects of sound assistance on fluidization behaviors were systematically investigated in a gas–solid acoustic fluidized bed. A model modified from Syamlal–O'Brien drag model was established. The original solid momentum equation was developed and an acoustic model was also proposed. The radial particle volume fraction, axial root‐mean‐square of bed pressure drop, granular temperature, and particle velocity in gas–solid acoustic fluidized bed were simulated using computational fluid dynamics (CFD) code Fluent 6.2. The results showed that radial particle volume fraction increased using modified drag model compared with that using the original one. Radial particle volume fraction was revealed as a parabolic concentration profile. Axial particle volume fraction decreased with the increasing bed height. The granular temperature increased with increasing sound pressure level. It showed that simulation values using CFD code Fluent 6.2 were in agreement with the experimental data. © 2009 American Institute of Chemical Engineers AIChE J, 2010