Numerical simulation of counter-current flow field in the downcomer of a liquid–solid circulating fluidized bed

A comprehensive study on the hydrodynamics in the downcomer of a liquid–solid circulating fluidized bed (LSCFB) is crucial in the control and optimization of the extraction process using an ion exchange LSCFB. A computational fluid dynamics model is proposed in this study to simulate the counter-current two-phase flow in the downcomer of the LSCFB. The model is based on the Eulerian–Eulerian approach incorporating the kinetic theory of granular flow. The predicted results agree well with our earlier experimental data. Furthermore, it is shown that the bed expansion of the particles in the downcomer is directly affected by the superficial liquid velocity in downcomer and solids circulation rate. The model also predicts the residence time of solid particles in the downcomer using a pulse technique. It is demonstrated that the increase in the superficial liquid velocity decreases the solids dispersion in the downcomer of the LSCFB.     

Onset velocity of circulating fluidization and particle residence time distribution:A CFD-DEM study

Until now,the onset velocity of circulating fluidization in liquid-solid fluidized beds has been defined by the turning point of the time required to empty a bed of particles as a function of the superficial liquid velocity,and is reported to be only dependent on the liquid and particle properties.This study presents a new approach to calculate the onset velocity using CFD-DEM simulation of the particle residence time distribution(RTD).The onset velocity is identified from the intersection of the fitted lines of the particle mean residence time as a function of superficial liquid velocity.Our results are in reasonable agreement with experimental data.The simulation indicates that the onset velocity is influenced by the density and size of particles and weakly affected by riser height and diameter.A power-law function is proposed to correlate the mean particle residence time with the superficial liquid velocity.The collisional parameters have a minor effect on the mean residence time of particles and the onset velocity,but influence the particle RTD,showing some humps and trailing.The particle RTD is found to be related to the particle trajectories,which may indicate the complex flow structure and underlying mechanisms of the particle RTD.     

Onset velocity of circulating fluidization and particle residence time distribution: A CFD–DEM study

Until now, the onset velocity of circulating fluidization in liquid–solid fluidized beds has been defined by the turning point of the time required to empty a bed of particles as a function of the superficial liquid velocity, and is reported to be only dependent on the liquid and particle properties. This study presents a new approach to calculate the onset velocity using CFD–DEM simulation of the particle residence time distribution (RTD). The onset velocity is identified from the intersection of the fitted lines of the particle mean residence time as a function of superficial liquid velocity. Our results are in reasonable agreement with experimental data. The simulation indicates that the onset velocity is influenced by the density and size of particles and weakly affected by riser height and diameter. A power-law function is proposed to correlate the mean particle residence time with the superficial liquid velocity. The collisional parameters have a minor effect on the mean residence time of particles and the onset velocity, but influence the particle RTD, showing some humps and trailing. The particle RTD is found to be related to the particle trajectories, which may indicate the complex flow structure and underlying mechanisms of the particle RTD.     

A numerical investigation into the solids phase chromatography using a combined continuous and discrete approach

Solids phase chromatography for particle classification is based on different retention times of particles with different properties when they are elutriated through a confined geometry. This work aims at a fundamental understanding of such a technology by using the combined continuous and discrete method. A packed bed is employed as the model confined geometry. The numerical method is compared first with experimental observations, followed by a parametric analysis of the effects on the flow hydrodynamics and solids behaviour of various parameters including the number of injected particles, the superficial gas velocity, the contact stiffness and the diameter ratio of the packed column to the packed particles. The results show that the modelling captures some important features of the flow of an injected pulse of fine particles in a packed bed. An increase in the number of injected particles or the superficial gas velocity reduces the retention time, whereas the contact stiffness does not show much effect over the range of 5 × 10^2 to 5× 10^4 N/m. It is also found that the effect on the retention time of the diameter ratio of the packed column to the packed particles seems complex showing a non-monotonous dependence.     

Gas–liquid mass transfer in the gas–liquid–solid mini fluidized beds

The gas–liquid–solid mini fluidized bed (GLSMFB) combines the advantages of fluidized bed and micro-reactor, and meets the requirements for safety and efficiency of green development of process industry. However, there are few studies on its flow performance and no studies on its mass and heat transfer performance. In this paper, the characteristics of gas–liquid mass transfer in a GLSMFB were studied in order to provide basic guidance for the study of GLSMFB reaction performance and application. Using CO₂ absorption by NaOH as the model process, the gas–liquid mass transfer performance of GLSMFB was investigated. The results show that the liquid volumetric mass transfer coefficient and the gas–liquid interfacial area both increase with the increase of the superficial gas velocity within the experimental parameter range under the same given superficial liquid velocity. At the same ratio of superficial gas to liquid velocity, the liquid volumetric mass transfer coefficient increases with the increase of the superficial liquid velocity. Fluidized solid particles strengthen the liquid mass transfer process, and the liquid volumetric mass transfer coefficient is about 13% higher than that of gas–liquid mini bubble column.     

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.     

Effects of various parameters on the attrition of bed material in a recirculating fluidized bed with a draft tube

We performed an experimental study to investigate the effects of various parameters on the attrition of bed material and its size distribution with increasing operation time in a recirculating fluidized bed (RCFB). The studied parameters included superficial velocity of fluidizing air, bed inventory, and spacing between the jet top and draft tube bottom (spacer height). The bed material was prepared from Indian Standard (IS) Grade I sand from sieves with a size range of 2.20–1.00 mm. Experiments were performed at ambient conditions, with the superficial air velocity ranging from 7.13–9.16 m/s, a bed inventory of 7–10 kg, spacing of 0.085 and 0.045 m between the jet top and draft tube bottom, and an operating time of 40 h. We investigated the influence of these parameters in terms of changes in the size distribution of particles, changes in the %-weight of particles of different size ranges, generation of particles with smaller diameters, the decrease of the downcomer bed height, variations in the coefficient of uniformity and coefficient of curvature, and material loss from entrainment of fines with increasing operation time. The mode of attrition was abrasion in all experiments. We found that with increasing operation time and other parameters (bed inventory, superficial air velocity, and spacer height) attrition of the bed material also increased. Generation and elutriation of fines were more pronounced at higher superficial air velocity, bed inventory, and spacer height.     

A numerical investigation into the solids phase chromatography using a combined continuous and discrete approach

Solids phase chromatography for particle classification is based on different retention times of particles with different properties when they are elutriated through a confined geometry.This work aims at a fundamental understanding of such a technology by using the combined continuous and discrete method.A packed bed is employed as the model confined geometry.The numerical method is compared first with experimental observations,followed by a parametric analysis of the effects on the flow hydrodynamics and solids behaviour of various parameters including the number of injected particles,the superficial gas velocity,the contact stiffness and the diameter ratio of the packed column to the packed particles.The results show that the modelling captures some important features of the flow of an injected pulse of fine particles in a packed bed. An increase in the number of injected particles or the superficial gas velocity reduces the retention time,whereas the contact stiffness does not show much effect over the range of 5×102 to5×104 N/m.It is also found that the effect on the retention time of the diameter ratio of the packed column to the packed particles seems complex showing a non-monotonous dependence.     

Local solids flow structure in a countercurrent liquid-upward and solids-downward fluidized bed (CCLSFB) system

The local solids holdup and local particle velocity in a Countercurrent Liquid-upward and Solids-downward Fluidized Bed (CCLSFB) were investigated in details using optical fiber probes with two different models in a Plexiglas column of 1.5 m in height and 7.0 cm in inner diameter. A new flow regime map including fluidized bed, transition, and flooding regimes was established. The axial solids holdup distribution is almost uniform at low liquid velocity and/or solids flowrate and becomes less uniform with higher solids holdup at the top of the column after the operating liquid velocity is reaching the flooding velocity. The radial solids holdup profile is also nearly flat with a slightly lower solids holdup in the near-wall region at low liquid velocity and solids flowrate but becomes nonuniform as the operating liquid velocity approaches the flooding velocity. Two equations were also proposed to correlate radial local solids holdups. The descending particle velocity in CCLSFB increases with the decrease of the liquid velocity and the increase of the solids flowrate. A generally uniform particle velocity distribution was found in the axial direction, as well as in the radial direction except for a small decrease near the wall. These results on the local solids flow structure would provide basic information and theoretical supports for the design and industrial application of CCLSFB.     

Hydrodynamics of a bubble-induced inverse fluidized bed reactor with a nanobubble tray

This paper reports on the hydrodynamics of a bubble-induced inverse fluidized bed reactor, using a nanobubble tray gas distributor, where solid particles are fluidized only by an upward gas flow. Increasing the gas velocity, the fixed layer of particles initially packed at the top of the liquid starts to move downwards, due to the rise of bubbles in this system, and then gradually expands downwards until fully suspended. The axial local pressure drops and standard deviation were examined to delineate the flow regime comprehensively under different superficial gas velocities. Four flow regimes (fixed bed regime, initial fluidization regime, expanded regime, and post-homogeneous regime) were observed and three transitional gas velocities (the initial fluidization velocity, minimum fluidization velocity, and homogeneous fluidization velocity) were identified to demarcate the flow regime. Three correlations were developed for the three transitional velocities. As the fine bubbles generated from the nanobubble tray gas distributor are well distributed in the entire column, the bed expansion process of the particles is relatively steady.     

Anti-wear beam effects on gas–solid hydrodynamics in a circulating fluidized bed

Anti-wear beams installed on water walls of circulating fluidized bed (CFB) boilers are one of the most effective ways to protect against water-wall erosion. Beam effects from, for example, beam size and superficial gas velocity were investigated on gas–solid hydrodynamics in a CFB test rig using CFD simulations and experimental methods. The downward flow of the wall layer solids is observed to be disrupted by the beam but is then restored some distance further downstream. When falling solids from the wall layer hit the anti-wear beam, the velocity of the falling solids decreases rapidly. A fraction of the solids accumulates on the beam. Below the beams, the falling solids have reduced velocities but upward-moving solids were observed on the wall. The effect of the beam increases with width and superficial gas velocity. Wear occurs mainly above the beam and its variation with width is different above to below the beam. There is an optimum width that, when combined with beam height, results in less erosion.     

Experimental analysis of volatile liquid injection into a fluidized bed

Experiments were conducted on a lab-scale fluidized bed to study the distribution of liquid ethanol injected into fluidized catalyst particles. Electrical capacitance measurements were used to study the liquid distribution inside the bed, and a new method was developed to determine the liquid content inside fluidized beds of fluid catalytic cracking particles. The results shed light on the complex liquid injection region and reveal the strong effect of superficial gas velocity on liquid distribution inside the fluidized bed, which is also affected by the imbibition of liquid inside particle pores. Particle internal porosity was found to play a major role when the changing mass of liquid in the bed was monitored. The results also showed that the duration of liquid injection affected liquid–solid contact inside the bed and that liquid–solid mixing was not homogeneous during the limited liquid injection time.     

FLOW DYNAMICS OF GAS-SOLID FLUIDIZED BEDS WITH EVAPORATIVE LIQUID INJECTION

1. Introduction Injections of evaporative liquids into fluidized solid parti- cles are routinely practiced in industrial processes involv- ing gas-solid fluidization systems such as fluid catalytic cracking, polymerization, and plastic coating (Fan et al., 2001). In the FCC riser system, heavy oil is injected into the system to evaporate rapidly by contact with the hot catalyst particles. Simultaneously, thermal and catalytic cracking reactions take place. During a polymerization process, a …     

Flow structure and bubble dynamics in supercritical water fluidized bed and gas fluidized bed: A comparative study

Supercritical water (SCW) fluidized bed is a new reactor concept for hydrogen production from biomass or coal gasification. In this paper, a comparative study on flow structure and bubble dynamics in a supercritical water fluidized bed and a gas fluidized bed was carried out using the discrete element method (DEM). The results show that supercritical water condition reduces the incipient fluidization velocity, changes regime transitions, i.e. a homogeneous fluidization was observed when the superficial velocity is in the range of the minimum fluidization velocity and minimum bubbling velocity even the solids behave as Geldart B powders in the gas fluidized bed. Bubbling fluidization in the supercritical water fluidized bed was formed after superficial velocity exceeds the minimum bubbling velocity, as in the gas fluidized bed. Bubble is one of the most important features in fluidized bed, which is also the emphasis in this paper. Bubble growth was effectively suppressed in the supercritical water fluidized bed, which resulted in a more uniform flow structure. By analyzing a large number of bubbles, bubble dynamic characteristics such as diameter distribution, frequency, rising path and so on, were obtained. It is found that bubble dynamic characteristics in the supercritical water fluidized bed differ a lot from that in the gas fluidized bed, and there is a better fluidization quality induced by the bubble dynamics in the supercritical water fluidized bed.     

Experimental study on filtration performance of the moving bed granular filter with axial flow

The filtration performance of the moving bed granular filter with axial flow (MBGF-AF) is investigated through a large cold experiment. The effect of different operation parameters on the filtration performance (collection efficiency, pressure drop) of the axial-flow moving bed filter is investigated in combination with the dust deposition effect and the mechanism of trapping dust by the capturing particles. The results show that the collection efficiency of MBGF-AF is enhanced by decreasing the superficial gas velocity, increasing the inlet dust concentration properly, or decreasing the moving velocity of the capturing particles. A model covering the above operation parameters is established to calculate the collection efficiency of the moving bed granular filter. It is used in a wide range of operating parameters for the MBGFs.     

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.     

Transport modeling of sedimenting particles in a turbulent pipe flow using Euler–Lagrange large eddy simulation

A volume-filtered Euler–Lagrange large eddy simulation methodology is used to predict the physics of turbulent liquid–solid slurry flow through a horizontal periodic pipe. A dynamic Smagorinsky model based on Lagrangian averaging is employed to account for the sub-filter scale effects in the liquid phase. A fully conservative immersed boundary method is used to account for the pipe geometry on a uniform cartesian grid. The liquid and solid phases are coupled through volume fraction and momentum exchange terms. Particle–particle and particle–wall collisions are modeled using a soft-sphere approach. Three simulations are performed by varying the superficial liquid velocity to be consistent with the experimental data by Dahl et al. (2003). Depending on the liquid flow rate, a particle bed can form and develop different patterns, which are discussed in light of regime diagrams proposed in the literature. The fluctuation in the height of the liquid-bed interface is characterized to understand the space and time evolution of these patterns. Statistics of engineering interest such as mean velocity, mean concentration, and mean streamwise pressure gradient driving the flow are extracted from the numerical simulations and presented. Sand hold-up calculated from the simulation results suggest that this computational strategy is capable of predicting critical deposition velocity.     

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.     

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.     

Lagrangian–Eulerian simulation of slugging fluidized bed

This work studies gas–solid slugging fluidized beds with Type-D particles, using two-dimensional simulations based on discrete element model (DEM). DEM performance is quantitatively validated by two commonly accepted correlations for determining slugging behavior. The voidage profiles simulated with bed height corresponding to Baeyens and Geldart (1974) correlation for onset of slugging demonstrate a transitional flow pattern from free bubbling to slugging. The present calculated values for the maximum slugging bed height are in good agreement with the correlation from Matsen et al. (1969). Simulations show that fluidized beds with Type-D particles can operate in the round-nosed slugging regime and also shows that wall slugs and square-nosed slugs tend to be formed with increase in superficial gas velocity and in bed height, respectively.