CFD study: Effect of pulsating flow on gas–solid hydrodynamics in FCC riser

Gas–solid flow in a fluid catalytic cracking (FCC) riser exhibits poor mixing in the form of a core–annulus flow pattern and a dense bottom/dilute top distribution of solids. To enhance gas–solid mixing, studies on dense fluidized beds have suggested using a pulsating flow of gas. The present study investigates the effect of pulsating flow on gas–solid hydrodynamics inside the FCC riser employing computational fluid dynamics. Two flow conditions are investigated: a cold flow of air-FCC catalyst in a pilot-scale riser and a reactive flow in an industrial-scale FCC riser. In the cold-flow riser, pulsating flows cause the slug flow of solids and thus increase the average solid accumulation in the flow domain and solid segregation towards the wall. In the industrial FCC riser, pulsating flows produce radial profiles that are more homogeneous. Pulsating flows further improve the conversion and yield in the initial few metres of height. At 7 m, the conversion from pulsating flow is 59%, compared with 44% in without pulsating flow. The results and analysis presented here will help optimize flow conditions in the circulating fluidized bed riser, in not only FCC but also applications such as fast pyrolysis and combustion.     

Modeling the electrostatic effect on the hydrodynamic behavior in FCC risers: From understanding to application

A CFD simulation was proposed to investigate the electrostatic effect on the hydrodynamic behavior of turbulent gas–solid flow in FCC risers. The simulation was first verified using the open experimental data with expected electrostatic effects observed in FCC risers. The influences of several operating parameters on the degree of electrification in FCC risers were analyzed, such as surface charge densities, pressure, gas velocity. It was noted that the gas velocity played a highly significant role compared with solid flux, while the effect of pressure was relatively weak. Further analysis showed that a much stronger electrostatic effect was found in small-scale FCC risers than their large-scale counterparts, and in addition, the major regions affected by the electrostatic charge depend on the scale of the riser. Finally, an external electric field was applied to optimize the flow field distribution in the FCC riser. The results of the electrostatic effects on the hydrodynamic behaviors in FCC risers are of great use in providing a reference for the optimization of FCC risers and their scaling.     

Solids flow diagram of a CFB riser using Geldart B-type powders

Riser operating modes are vital to designing a circulating fluidized bed (CFB) reactor for a required process of either a gas-solid or a gas-catalytic nature. Different operating modes provide different solids’ residence times and mixing behaviors, which define the reactions’ efficiency and yield. The literature demonstrates distinct operating modes resulting from observed differences in slip factors and the range of particle velocities and their associated residence time distribution. The present research uses positron emission particle tracking (PEPT) in a riser of B-type bed material to determine the different operating modes by measuring (i) particle velocities and residence time distribution, (ii) population densities of these particles in the cross-sectional area of the riser, and (iii) solids flow pattern at the bottom of the riser. Data treatment defines four distinct solids hold-up regimes in the riser and proposes a "phase diagram" depicting the existence of the different operating modes (dilute, dense, core-annulus and combined) as a function of the superficial gas velocity and solids circulation flux in the riser. The delineated regimes have good agreement with available literature data and known industrial operations. Comparison with literature data for risers using A-type powders is also fair. The diagram enables CFB designers to better delineate operating characteristics.     

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.     

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.     

3D full-loop simulation and experimental verification of gas-solid flow hydrodynamics in a dense circulating fluidized bed

Because of their advantages of high efficiency and low cost, numerical research methods for large-scale circulating fluidized bed （CFB） apparatus are gaining ever more importance. This article presents a numer- ical study of gas-solid flow dynamics using the Eulerian granular multiphase model with a drag coefficient correction based on the energy-minimization multi-scale （EMMS） model. A three-dimensional, full-loop, time-dependent simulation of the hydrodynamics of a dense CFB apparatus is performed. The process parameters （e.g., operating and initial conditions） are provided in accordance with the real experiment to enhance the accuracy of the simulation. The axial profiles of the averaged solid volume fractions and the solids flux at the outlet of the cyclone are in reasonable agreement with experimental data, thereby verifying the applicability of the mathematical and physical models. As a result, the streamline in the riser and standpipe as well as the solids distribution contours at the cross sections is analyzed. Computational fluid dynamics （CFD） serves as a basis for CFB modeling to help resolve certain issues long in dispute but difficult to address experimentally. The results of this study provide the basis of a general approach to describing dynamic simulations of gas-solid flows.     

Extending EMMS-based models to CFB boiler applications

Recently, EMMS-based models are being widely applied in simulations of high-throughput circulating fluidized beds (CFBs) with fine particles. Its use for low flux systems, such as CFB boiler (CFBB), still remains unexplored. In this work, it has been found that the original definition of cluster diameter in EMMS model is unsuitable for simulations of the CFB boiler with low solids flux. To remedy this, we propose a new model of cluster diameter. The EMMS-based drag model (EMMS/matrix model) with this revised cluster definition is validated through the computational fluid dynamics (CFD) simulation of a CFB boiler.     

Eulerian–Lagrangian simulation of distinct clustering phenomena and RTDs in riser and downer

Numerical simulation of fully developed hydrodynamics of a riser and a downer was carried out using an Eulerian–Lagrangian model, where the particles are modeled by the discrete element method (DEM) and the gas by the Navier–Stokes equations. Periodic flow domain with two side walls was adopted to simulate the fully developed dynamics in a 2D channel of 10 cm in width. All the simulations were carried out under the same superficial gas velocity and solids holdup in the domain, starting with a homogenous state for both gas and solids, and followed by the evolution of the dynamics to the heterogeneous state with distinct clustering in the riser and the downer. In the riser, particle clusters move slowly, tending to suspend along the wall or to flow downwards, which causes wide residence time distribution of the particles. In the downer, clusters still exist, but they have faster velocities than the discrete particles. Loosely collected particles in the clusters move in the same direction as the bulk flow, resulting in plug flow in the downer. The residence time distribution (RTD) of solids was computed by tracking the displacements of all particles in the flow direction. The results show a rather wide RTD for the solids in the riser but a sharp peak RTD in the downer, much in agreement with the experimental findings in the literature. The ensemble average of transient dynamics also shows reasonable profiles of solids volume fraction and solids velocity, and their dependence on particle density.     

Experimental study of the effect of friction phenomena on actual and calculated inventory in a small-scale CFB riser

Accurate information concerning riser inventory in a fluidized bed is required in some applications such as the calcium looping process,because it is related to the CO_2 capture efficiency of the system.In a circulating fluidized bed(CFB),the riser inventory is normally calculated from the riser pressure drop;however,the friction and the acceleration phenomena may have a significant influence on the total riser pressure drop.Therefore,deviation may occur in the calculation from the actual mass.For this reason the magnitude of the friction and the acceleration pressure drop in the entire riser is studied in small-scale risers.Two series of studies were performed:the first one in a scaled cold model riser of the 10 kW_(th)facility,and the second one in the 10kW_(th) fluidized bed riser under process conditions.The velocities were chosen to comply with the fluidization regimes suitable for the calcium looping process,namely,the turbulent and the fast.In cold-model experiments in a low-velocity turbulent fluidization regime,the actual weight(static pressure drop) of the particles is observed more than the weight calculated from a recorded pressure drop.This phenomenon is also repeated in pilot plant conditions.In the cold-model setup,the friction and acceleration pressure drop became apparent in the fast fluidization regime,and increased as the gas velocity rose.Within calcium looping conditions in the pilot plant operation,the static pressure drop was observed more than the recorded pressure drop.Therefore,as a conservative approach,the influence of friction pressure drop may be neglected while calculating the solid inventory of the riser.The concept of transit inventory is introduced as a fraction of total inventory,which lies in freefall zones of the CFB system.This fraction increases as gas velocity rises.     

The elimination of severe slugging-experiments and modeling

Severe slugging can occur in a pipeline-riser system operating at low liquid and gas rates. The flow of gas into the riser can be blocked by liquid accumulation at the base of the riser. This can cause formation of liquid slugs of a length equal to or longer than the height of the riser. A cyclic process results in which a period of no liquid production into the separator occurs, followed by a period of very high liquid production. This study is an experimental and theoretical investigation of two methods for eliminating this undesirable phenomenon, using choking and gas lift. Choking was found to effectively eliminate or reduce the severity of the slugging. However, the system pressure might increase to some extent. Gas lift can also eliminate severe slugging. While choking reduces the velocities in the riser, gas lift increases the velocities, approaching annular flow. It was found that a relatively large amount of gas was needed before gas injection would completely stabilize the flow through the riser. However, gas injection reduces the slug length and cycle time, causing a more continuous production and a lower system pressure. Theoretical models for the elimination of severe slugging by gas lift and choking have been developed. The models enable the prediction of the flow behavior in the riser. One model is capable of predicting the unstable flow conditions for severe slugging based on a static force balance. The second method is a simplified transient model based on the assumption of a quasi-equilibrium force balance. This model can be used to estimate the characteristics of the flow, such as slug length and cycle time. The models were tested against new severe slugging data acquired in this study. An excellent agreement between the experimental data and the theoretical models was found.     

STUDY ON THE OVERALL PRESSURE BALANCE OF A DOWNFLOW CIRCULATING FLUIDIZED BED SYSTEM

A pressure balance model for a circulating fluidized bed unit that incorporates a downer has been proposed. The model predictions were validated with the experimental data obtained from a special cold-model circulating fluidized bed. Comparison of the operation stability between a CFB downer and a CFB riser has been carried out. Only one critical gas velocity exists in the CFB-riser for a given riser solids flux, while there can be many critical gas velocities for the operation of a CFB downer. Therefore, it is possible to achieve high solids concentration in a CFB downer if appropriate operating conditions are used.     

STUDY ON THE OVERALL PRESSURE BALANCE OF A DOWNFLOW CIRCULATING FLUIDIZED BED SYSTEM

A pressure balance model for a circulating fluidized bed unit that incorporates a downer has been proposed. The model predictions were validated with the experimental data obtained from a special cold-model circulating fluidized bed. Comparison of the operation stability between a CFB downer and a CFB riser has been carried out. Only one critical gas velocity exists in the CFB-riser for a given riser solids flux, while there can be many critical gas velocities for the operation of a CFB downer. Therefore, it is possible to achieve high solids concentration in a CFB downer if appropriate operating conditions are used.     

Effect of baffles on performance of fluid catalytic cracking riser

Increasing demand of automobile fuel and a need to process heavier crude oil makes it imperative to find improvements to the design of existing fluid catalytic cracking (FCC) units. Several modifications to the design of the riser section of FCC units have been suggested in previous studies including: improved feed nozzle designs, multiple nozzle configurations, internal baffles, and novel two-stage-riser systems. In this study, we investigate the effects of baffles on the performance of FCC risers using computational fluid dynamics simulations. In this study, predictions from a basis model (without baffles) are compared with those from four different configurations including: (i) 5-cm baffles at 5-m spacing, (ii) 7.5-cm baffles at 5-m spacing, (iii) 10-cm baffles with 5-m spacing, (iv) 10-cm baffles at 2.5-m spacing, and (v) 10-cm baffles at 1-m spacing. The baffles force the catalyst away from walls toward the center of the riser, enhancing the radial dispersion of the catalyst and the heat transfer inside the riser. The use of longer baffles and smaller spacings further increases the dispersion, yielding more homogeneous radial profiles. The changes in the radial dispersion result in variations in the conversion, yields, and pressure drops. The baffles increase conversion of vacuum gas oil (VGO) and the yield of gasoline. However, the simulations showed that longer baffles and a larger number of baffles did not always give a higher yield or higher conversion. Among the simulated configurations, the 5-cm baffles at 5-m spacing gave the highest conversion of VGO, whereas the 10-cm baffles at 1-m spacing resulted in the highest yield of the gasoline. Thus, rational optimization of baffle configurations is required to achieve optimal performance.     

Numerical study on mitigating severe slugging in pipeline/riser system with wavy pipe

The gas/liquid two-phase flow in pipeline/wavy-pipe/riser systems was investigated numerically with CFD. A CFD model of the pipeline/wavy-pipe/riser system was obtained by adding a wavy pipe to the model of the pipeline/riser system verified by the experimental data previously. The effects of the geometrical parameters and location of the wavy pipe on its performance of slug mitigation and flow characteristics in pipeline/wavy-pipe/riser systems were examined through the CFD models. With the increase of the amplitude or length of the wavy pipe, the slug in the pipeline/riser system becomes shorter. The optimum location of the wavy pipe in the pipeline exists for a pipeline/riser system and a wavy pipe at given operating conditions. The CFD modelling provides a feasible and flexible way to investigate the effectiveness of the wavy pipes on mitigating severe slugging in pipeline/riser systems.     

Experimental study of the effect of friction phenomena on actual and calculated inventory in a small-scale CFB riser

Accurate information concerning riser inventory in a fluidized bed is required in some applications such as the calcium looping process, because it is related to the CO₂ capture efficiency of the system. In a circulating fluidized bed (CFB), the riser inventory is normally calculated from the riser pressure drop; however, the friction and the acceleration phenomena may have a significant influence on the total riser pressure drop. Therefore, deviation may occur in the calculation from the actual mass. For this reason the magnitude of the friction and the acceleration pressure drop in the entire riser is studied in small-scale risers. Two series of studies were performed: the first one in a scaled cold model riser of the 10 kW_th facility, and the second one in the 10 kW_th fluidized bed riser under process conditions. The velocities were chosen to comply with the fluidization regimes suitable for the calcium looping process, namely, the turbulent and the fast. In cold-model experiments in a low-velocity turbulent fluidization regime, the actual weight (static pressure drop) of the particles is observed more than the weight calculated from a recorded pressure drop. This phenomenon is also repeated in pilot plant conditions. In the cold-model setup, the friction and acceleration pressure drop became apparent in the fast fluidization regime, and increased as the gas velocity rose. Within calcium looping conditions in the pilot plant operation, the static pressure drop was observed more than the recorded pressure drop. Therefore, as a conservative approach, the influence of friction pressure drop may be neglected while calculating the solid inventory of the riser. The concept of transit inventory is introduced as a fraction of total inventory, which lies in freefall zones of the CFB system. This fraction increases as gas velocity rises.     

Effect of static bed height in the upper fluidized bed on flow behavior in the lower riser section of a coupled reactor

To study olefin reduction by using an auxiliary reactor for FCC naphtha upgrading, a large-scale cold model of a riser-bed coupled to an upper fluidized bed was established. The effect of static bed height in the upper fluidized bed on particle flow behavior in the lower riser was investigated experimentally. A restriction index of solids holdup was used to evaluate quantitatively the restrictive effect of the upper fluidized bed. Experimental results show that, under the restrictive effect of the upper fluidized bed, the riser could be divided into three regions in the longitudinal direction: accelerating, fully developed and restriction. The axial distribution of solids holdup in the riser is characterized by large solids holdup in the top and bottom sections and small solids holdup in the middle section. Overall solids holdup increased with increasing static bed height in the upper fluidized bed, while particle velocity decreased. Such restrictive effect of the upper fluidized bed could extend from the middle and top sections to the whole riser volume when riser outlet resistance is increased, which increases with increasing static bed height in the upper fluidized bed. The upper bed exerts the strongest restriction on the area close to the riser outlet.     

Effect of static bed height in the upper fluidized bed on flow behavior in the lower riser section of a coupled reactor

To study olefin reduction by using an auxiliary reactor for FCC naphtha upgrading, a large-scale cold model of a riser-bed coupled to an upper fluidized bed was established. The effect of static bed height in the upper fluidized bed on narticle flow behavior in the lower riser was investigated experimentally. A restriction index of solids holdup was used to evaluate quantitatively the restrictive effect of the upper fluidized bed. Experimental results show that, under the restrictive effect of the upper fluidized bed, the riser could be divided into three regions in the longitudinal direction： accelerating, fully developed and restriction. The axial distribution of solids holdup in the riser is characterized by large solids holdup in the top and bottom sections and small solids holdup in the middle section. Overall solids holdup increased with increasing static bed height in the upper fluidized bed, while particle velocity decreased. Such restrictive effect of the upper fluidized bed could extend from the middle and top sections to the whole riser volume when riser outlet resistance is increased, which increases with increasing static bed height in the upper fluidized bed. The upper bed exerts the strongest restriction on the area close to the riser outlet.     

Numerical simulation and experimental validation of gas–solid flow in the riser of a dense fluidized bed reactor

Gas–solid flow in the riser of a dense fluidized bed using Geldart B particles (sand), at high gas velocity (7.6–15.5 m/s) and with comparatively high solid flux (140–333.8 kg/m² s), was investigated experimentally and simulated by computational fluid dynamics (CFD), both two- and three-dimensional and using the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd and EMMS drag models. The results predicted by EMMS drag model showed the best agreement with experimental results. Calculated axial solids hold-up profiles, in particular, are well consistent with experimental data. The flow structure in the riser was well represented by the CFD results, which also indicated the cause of cluster formation. Complex hydrodynamical behaviors of particle cluster were observed. The relative motion between gas and solid phases and axial heterogeneity in the three subzones of the riser were also investigated, and were found to be consistent with predicted flow structure. The model could well depict the difference between the two exit configurations used, viz., semi-bend smooth exit and T-shaped abrupt exit. The numerical results indicate that the proposed EMMS method gives better agreement with the experimental results as compared with the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd models. As a result, the proposed drag force model can be used as an efficient approach for the dense gas–solid two-phase flow.     

Numerical study on hydrodynamics of gas–solids circulating fluidized bed with L-valve

L-valve is often used as a non-mechanical valve for the circulation of solids in gas–solids fluidized bed (GSFB) due to its advantages in simple construction and easy control. The information on solids circulation rate as well as the hydrodynamics performance of the CFB with L-valve is of great importance for its better control and design. This paper proposes a Eulerian-Eulerian approach based numerical model integrating the computational fluid dynamics (CFD) with turbulent model, the kinetic theory of granular flow (KTGF) and the drag model, thus the solids circulation rate and the local phase velocity as well as solids volume fraction can be predicted simultaneously. With this model, the hydrodynamics performance of the full loop GSCFB with a L-valve is analyzed in detail. It is found that the drag model affects the simulation significantly and the (energy minimization multiscale) EMMS method shows good performance in the full-loop simulation of GSCFB.     

Quantitative measurement of solids distribution in gas–solid riser flows using electrical impedance tomography and gamma densitometry tomography

An electrical impedance tomography (EIT) system has been developed to non-invasively measure particle distributions in the riser of a pilot-scale circulating fluidized bed (CFB). Although EIT systems have often been applied to yield qualitative information about gas–solid flows, the present EIT system yields quantitative information that is validated by comparison to a gamma densitometry tomography (GDT) system. EIT and GDT were applied to the CFB riser (14-cm inner diameter, 5.77-m height) containing fluid catalytic cracking particles in air. The flows examined were annular with a dilute core and had average and near-wall solids volume fractions up to 0.25 and 0.66, respectively. For all cases, the average and near-wall solids volume fractions from EIT and GDT agreed to within 0.03 and 0.07, respectively. This good agreement suggests that, where feasible, EIT can be used in place of GDT, which is advantageous since EIT systems are often safer, less expensive, and faster than GDT systems.