Fluidization behavior in a circulating slugging fluidized bed reactor. Part I: Residence time and residence time distribution of polyethylene solids

Square nosed slugging fluidization behavior in a circulating fluidized bed riser using a polyethylene powder with a very wide particle size distribution was studied. In square nosed slugging fluidization the extent of mixing of particles of different size depends on the riser diameter, gas velocity, hold up and solids flux in the riser. Depending on the operating conditions the particle residence time distribution of a riser in the slugging fluidization regime can vary from that of a plug flow reactor to that of a well-mixed system.Higher gas velocities cause shorter particle residence times because of a significant decrease in the hold-up of particles in the riser at higher gas velocities. A higher solids flux also shortens the average residence time. Both influences have been quantified for a given polyethylene-air system.Residence time and residence time distribution were determined for different particle size and the influence of gas velocity, solids flux, hold up and riser diameter was studied. When comparing data from segregation and residence time experiments it is clear that segregation data can predict the spread in residence time as a function of overall residence time, particle size and gas velocity. The differential velocity between small and large particles found in the segregation experiments can predict the spread in residence time as found in the residence time distribution experiments with a powder with a broad particle size distribution. Raining of particles through the slugs was studied as a function of plug length, gas velocity and pulse length. It was found that raining is not the determining mechanism for segregation of particles.     

CFD simulation of hydrodynamics of gas-liquid-solid fluidised bed reactor

A three dimensional transient model is developed to simulate the local hydrodynamics of a gas-liquid-solid three-phase fluidised bed reactor using the computational fluid dynamics (CFD) method. The CFD simulation predictions are compared with the experimental data of Kiared et al. [1999. Mean and turbulent particle velocity in the fully developed region of a three-phase fluidized bed. Chemical Engineering & Technology 22, 683-689] for solid phase hydrodynamics in terms of mean and turbulent velocities and with the results of Yu and Kim [1988. Bubble characteristics in the racial direction of three-phase fludised beds. A.I.Ch.E. Journal 34, 2069-2072; 2001. Bubble-wake model for radial velocity profiles of liquid and solid phases in three-phase fluidised beds. Industrial and Engineering Chemistry Research 40, 4463-4469] for the gas and liquid phase hydrodynamics in terms of phase velocities and holdup. The flow field predicted by CFD simulation shows a good agreement with the experimental data. From the validated CFD model, the computation of the solid mass balance and various energy flows in fluidised bed reactors are carried out. The influence of different interphase drag models for gas-liquid interaction on gas holdup are studied in this work.     

CFD simulation and experiments of dynamic parameters in gas–solid fluidized bed

Particle and bubble motion plays an important role in determining the hydrodynamic characteristics of a fluidized system. The dynamic parameters of a fluidized bed are reflection of the complex correlation between particle–particle and particle–bubble in a system. A two-dimensional Eulerian–Eulerian model integrating the kinetic theory of granular flow is used to simulate the bubble and linear low density polyethylene (LLDPE) particle dynamic behavior in a gas–solid fluidized bed. The simulated method is validated by pressure fluctuation experiment. The computed vertical turbulent energy spectrum of particles is applied to identify the particle motion intensity and the inhomogeneity of turbulent energy dissipation. The energy spectrum captures the Levy–Kolmogorov law in inertial range at high frequency. Furthermore, the flatness factors of wavelet decomposition coefficients of particle fluctuation velocity are for the first time introduced to analyze the intermittence caused by coherent structures in the flow field. The results show that the intermittence in dissipation range is much stronger than that in energy-containing and inertial range, and reinforces rapidly as the radial distance and the bed height increase. Moreover, the acoustic emission (AE) energy is found to be able to indicate the flow regimes. By combing granular temperature and AE energy, the relationship between the spatial distribution of granular temperature and the flow regimes is established. To get more detail of bubble motion behavior, the power spectrum of voidage fluctuation is analyzed. This work provides valuable insights into the dynamic characteristics and the flow field information of a gas–solid fluidized bed by CFD simulation.     

CFD modelling of a liquid-solid fluidized bed

A multifluid Eulerian computational fluid dynamics (CFD) model with granular flow extension is used to simulate a liquid-solid fluidized bed. The numerical simulations are evaluated qualitatively by experimental data from the literature and quantitatively by comparison with new experimental data. The effects of mesh size, time step and convergence criteria are investigated. Varying the coefficient of restitution did not alter the results significantly. The Gidaspow drag relationship predicted a higher voidage than the Wen and Yu drag law. Two different liquid distributors (uniform and non-uniform) were simulated and compared, but a better representation of the geometry of the distributor plate did not greatly influence the results. Qualitatively, the simulations show trends similar to experimental trends reported by various authors. The predictions are also compared with new experimental results for 1.13 mm glass spheres at a wide variety of superficial liquid velocities (0.0085-0.110 m/s) and two different temperatures (12 and ) significantly affecting the liquid viscosity. The CFD model predictions are within 5% of the steady-state experimental data and show the correct trend with variation in viscosity.     

Simulation of bubbling fluidized bed of fine particles using CFD

Computational fluid dynamics (CFD) simulation for bubbling fluidized bed of fine particles was carried out. The reliability and accuracy of CFD simulation was investigated by comparison with experimental data. The experimental facility of the fluidized bed was 6 cm in diameter and 70 cm in height and an agitator of pitched-blade turbine type was installed to prevent severe agglomeration of fine particles. Phosphor particles were employed as the bed material. Particle size was 22 μm and particle density was 3,938 kg/m³. CFD simulation was carried by two-fluid module which was composed of viscosity input model and fan model. CFD simulation and experiment were carried out by changing the fluidizing gas velocity and agitation velocity. The results showed that CFD simulation results in this study showed good agreement with experimental data. From results of CFD simulation, it was observed that the agitation prevents agglomeration of fine particles in a fluidized bed.     

Flow field and residence time distribution simulation of a cross-flow gas-liquid wastewater treatment reactor using CFD

A three-dimensional Eulerian-Eulerian two-phase approach has been used for the simulation of a cross-flow gas-liquid wastewater treatment reactor. Two different turbulence models have been tested: the k-ε and Reynolds Stress Model (RSM) models. Bubble induced turbulence source terms have been added to these models. Numerical results have been validated using Laser Doppler Velocimetry (LDV) measurements. Simulations with both turbulence models successfully predicted the hydrodynamics of the reactor. Then particle tracking with a stochastic approach has been used to calculate residence time distributions (RTD) with the flow previously simulated. It has been shown that dispersion in the reactor is primarily due to turbulence. Results have been compared with experimental RTD for various liquid and gas flowrates both on a bench scale and full scale plant. The RSM model accurately predicted the dispersion whereas the standard k-ε model slightly underestimated the dispersion.     

Bed height and material density effects on fluidized bed hydrodynamics

Characterizing the hydrodynamics of a fluidized bed is of vital importance to understanding the behavior of this multiphase flow system. Minimum fluidization velocity and gas holdup are two of these key characteristics. Experimental studies addressing the effects of bed height and material density on the minimum fluidization velocity and gas holdup were carried out in this study using a 10.2 cm diameter cylindrical fluidized bed. Three different Geldart type-B particles were tested: glass beads, ground walnut shell, and ground corncob, with material densities of 2600, 1300, and 1000 kg/m³, respectively. The particle size range was selected to be the same for all three materials and corresponded to 500–600 μm. In this study, five different bed height-to-diameter ratios were investigated: H/D=0.5, 1, 1.5, 2, and 3. Minimum fluidization velocity was determined for each H/D ratio using pressure drop measurements. Local time-average gas holdup was determined using non-invasive X-ray computed tomography imaging. Results show that minimum fluidization velocity is not affected by the change in bed height. However, as the material density increased, the minimum fluidization velocity increased. Finally, local time-average gas holdup values revealed that bed hydrodynamics were similar for all bed heights, but differed when the material density was changed.     

Multifluid Eulerian modeling of dense gas-solids fluidized bed hydrodynamics: Influence of the dissipation parameters

Computational fluid dynamic (CFD) models must be thoroughly validated before they can be used with confidence for designing fluidized bed reactors. In this study, validation data were collected from a fluidized bed of (Geldart's group B) alumina particles operated at different gas velocities involving two fluidization hydrodynamic regimes (bubbling and slugging). The bed expansion, height of bed fluctuations and frequency of fluctuations were measured from videos of the fluidized bed. The Eulerian-Eulerian two fluid model MFIX was used to simulate the experiments. Two different models for the particle stresses—Schaeffer [Syamlal, M., Rogers, W., O’Brien, T.J., 1993. MFIX documentation: theory guide. Technical Report DOE/METC-94/1004 (DE9400087), Morgantown Energy Technology Centre, Morgantown, West Virginia (can be downloaded from Multiphase Flow with Interphase eXchanges (MFIX) website 〈http://www.mfix.org〉); Schaeffer, D.G., 1987. Instability in the evolution equations describing incompressible granular flow. Journal of Differential Equations 66, 61-74.] and Princeton [Srivastava, A., Sundaresan, S., 2003. Analysis of a frictional-kinetic model for gas-particle flow. Powder Technology 129(1-3), 72-85.] models—and different values of the restitution coefficient and internal angle of friction were evaluated. 3-D simulations are required for getting quantitative and qualitative agreement with experimental data. The results from the Princeton model are in better agreement with data than that from the Schaeffer model. Both free slip and Johnson-Jackson boundary conditions give nearly identical results. An increase in coefficient of restitution (e) from 0.8 to 1 leads to larger bed expansions and lower heights of fluctuations in the bubbling regime, whereas it leads to unchanged bed expansion and to a massive reduction in the height of fluctuations in the slugging regime. The angle of internal friction (φ) in the range 10-40^° does not affect the bed expansion, but its reduction significantly reduces the height of fluctuations.     

CFD simulation of gas and solids mixing in FCC strippers

The gas and solid mixing in fluid catalytic cracking strippers with and without internals were investigated using computational fluid dynamics simulations. The Eulerian–Eulerian two‐fluid model coupled with the modified Gidaspow drag model was used to simulate the gas‐solid flow behavior. The grid independency study and the comparison of 2D and 3D simulations were carried out first. The residence time distribution model and axial dispersion model were utilized to obtain the parameters indicating the back‐mixing degree, such as mean residence time, dimensionless variance and Peclet number of gas and solids. Moreover, the influence of bubble size and gas/solid flow distribution on the mass transfer between the bubble and emulsion phase were also analyzed. The results show that the baffles in the V‐baffle stripper can efficiently enhance the gas and solids mixing, reduce the back‐mixing degree of gas and solids, strengthen the mass transfer between the bubble and emulsion phase, and hence improve the stripping efficiency. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Axial dispersion of gas in a circulating fluidized bed

An axial dispersion of gas in a circulating fluidized bed was investigated in a fluidized bed of 4.0 cm I.D. and 279 cm in height. The axial dispersion coefficient of gas was determined by the stimulus-response method of trace gas of CO₂. The employed particles were 0.069 mm and 0.147 mm silica-sand. The results showed that axial dispersion coefficients were increased with gas velocity and solid circulation rates as well as suspension density. The experimentally determined axial dispersion coefficients in this study were in the range of 1.0-3.5 m²/s.     

Hydrodynamics of a superheated steam vacuum fluidized bed

Hydrodynamics of a superheated steam vacuum fluidized bed was experimentally studied. In these experiments, eight different types of large particles (1970–7430 μm) were used. In all cases, a behavior similar to that found in an air fluidized bed was observed. The minimum fluidization velocity was found to be increasing with decreasing operating pressure. In the case of employing superheated steam, the minimum fluidization conditions are established at a lower velocity than using air as the fluidizing medium. These tendencies are attributed to the variation of the mean free path of molecules. On the other hand, the experiments showed that the bed voidage in the minimum fluidization conditions is almost insensitive to the variation of the operating pressure. Several equations were developed to predict the minimum fluidization velocity. The values provided by these equations were compared with the experimental data as well as with the predictions of the correlations presented in the technical literature.     

Hydrodynamic study of a toroidal fluidized bed reactor

Hydrodynamic behavior of a newly developed toroidal fluidized bed reactor is studied in this work. The reactor has a gas distributor consisting of angled blades in an annular ring at the reactor bottom. The driving force for particles to move over the distributing blades comes from the velocity head of gas jets accelerated upon entering the blade spacing. Relevant hydrodynamic behaviors are measured with various inert materials in a pilot scale 400-mm toroidal fluidized bed reactor. The observed hydrodynamic behavior is found to be essentially predictable at ambient temperature by conventional hydrodynamic models. Fine particle tracking on the reactor wall is clearly observed through oxidation of zinc dross at a bed temperature of around 1120°C, and is simulated on the basis of a simplified mathematical model. Hydrodynamic issues, such as particle flying trajectory and retention time in the reactor, are discussed based on the developed model.     

Hydrodynamic similarity in circulating fluidized bed risers

Hydrodynamic similarity in the fully developed zone of co-current upward gas-solid two-phase flow systems under different operating conditions was investigated by measuring the axial profiles of pressure gradient, radial profiles of solid concentration and particle velocity in two circulating fluidized bed (CFB) risers of 15.1 and 10.5 m high, with FCC and sand particles, respectively. The experimental data obtained from this work and in the literature show that when the scaling parameter, G_s/(ρ_pU_g), is modified as , a detailed hydrodynamic similitude of the gas-solid flow in the fully developed zone of the risers under different operating conditions can be achieved. Furthermore, the experimental results from different gas-solid flow systems also show that as long as remains constant, there is the same solid concentration in the fully developed zone of different CFB risers with different particles. With the same , the local solid concentrations, the descending particle velocities, the cluster frequencies and the solid concentrations inside clusters in the fully developed zone of the risers all display the same axial and radial distribution, respectively. In other words, the empirical similarity parameter, , appears to have incorporated the effects of operating parameters (G_s and U_g), so that, the gas-solid flow in the fully developed zone of CFB risers under those different operating conditions but having the same shows similar micro- and macro-hydrodynamic characteristics. The study shows that the empirical similarity parameter, , is also independent of the upward gas-solid flow systems.     

泰勒反应器中流体流动及停留时间分布研究

以水为介质对泰勒反应器中的流动状况和停留时间及其分布(RTD)进行了研究,并应用计算流体力学(CFD)技术对反应器进行了流场模拟和RTD计算。结果表明,在实验范围内,泰勒反应器中停留时间分布受内筒转速、轴向流动速率等因素影响,基于流体力学计算结果与实验结果基本相当。     

Phase hold-up and critical fluidization velocity in a three-phase inverse fluidized bed

We studied the hydrodynamic characteristics of a three-phase inverse fluidized bed made of a transparent acrylic column of 0.115 m inner diameter and 2 m heights. Air, water and polyethylene particles were used as the gas, liquid and solid phase, respectively. We used both hydrophobic low density polyethylene (LDPE) and hydrophilic LDPE as solid phase, and distilled water as liquid phase, and filtered air as gas phase. The LDPE was chemically treated by chlorosulfonic acid to change the surface property from hydrophobic to hydrophilic. We tried to solely investigate the effect of the surface hydrophilicity of polymeric particles on the phase holdup and the critical fluidization velocity of three-phase inverse fluidization. Thus, we measured the static pressure and eventually observed critical fluidization velocity. Critical fluidization velocity became smaller in case of using MDPE hydrophobic particles than LDPE hydrophilic particles. This was thought to be due to the retardation of rising bubbles near hydrophobic particles and, subsequently, the increase of gas hold-up.     

Interpretation of the hydrodynamic behaviour in a conical fluidized bed dryer

Wet batches of placebo pharmaceutical granule were dried at inlet superficial gas velocities of 0.64 and 1.3 m/s in a Glatt GPCG-1 fluidized bed. Using pressure fluctuation analysis, the hydrodynamic behaviour indicates a transition from a multiple bubbling regime to a coalescence dominated regime as drying proceeds. The transitional fluidization behaviour is linked to the physical mechanisms associated with the constant and falling rate periods of drying porous materials. Excess surface moisture present during the constant rate period increases interparticle forces through liquid bridging. These liquid bridges stabilize the bed structure which limits bubble formation in the bed. Once the falling rate period is reached, the liquid bridges cannot be maintained and bubble coalescence increases. The resulting bubbling bed hydrodynamics can be explained using the simple two-phase model proposed by Toomey and Johnstone [1952. Gas fluidization of solid particles. Chemical Engineering Progress 48, 220-226] using the full support velocity and bed voidage characteristics of the granule at varying moisture contents.     

CFD simulation of barium carbonate precipitation in a fluidized bed reactor

CFD techniques are used to study the precipitation of barium carbonate in a solid–liquid fluidized bed reactor. Experimental analysis of the hydrodynamic behaviour for a neutralization reaction in the fluidized bed column, followed by CFD simulations is carried out using different reaction models. The Eddy Dissipation model, the Eddy Dissipation model-MTS and the Eddy Dissipation Concept micro-mixing models are tested in order to simulate the acid–base instantaneous reaction.     

Flow structure of the solids in gas-solid fluidized beds

Existence of clusters in dense fluidized beds was investigated by analyzing the time-position data of a tracer obtained in several radioactive particle tracking experiments. It was found that in the case of sand particles, more gas passes through the bed as bubbles with increasing the superficial gas velocity and in the case of FCC powder, flow of the gas through the bed as bubbles does not increase in the turbulent fluidization regime. Cluster diameters were estimated from their velocities and found that descending clusters are generally larger than ascending ones and the size of both increases with increasing the superficial gas velocity. Bubble velocities evaluated in this work are in good agreement with the correlations in the bubbling regime of the fluidization available in the literature.     

流化床密相区颗粒扩散系数的CFD数值预测

应用离散颗粒模型直观获得颗粒运动情况,并从单个颗粒和气泡作用的角度分析颗粒运动和混合,证实气泡在床层中上升、在床层表面爆破以及气泡上升引起的乳化相下沉运动对颗粒混合起关键作用。应用基于颗粒动理学的双流体模型系统地对床宽分别为0.2、0.4、0.8 m的二维流化床在鼓泡区和湍动区的气固两相流动行为进行数值模拟。受离散颗粒模型启发,在双流体模型计算结果基础上,引入理想示踪粒子技术计算床内平均颗粒扩散系数。计算结果表明,颗粒横向扩散系数(D_x)总体上随流化风速增大而增大,但受床体尺寸影响较大;颗粒轴向扩散系数随流化风速增大而增大,受床体尺寸影响较弱。文献报道的密相区颗粒横向扩散系数分布在10^-4～10^-1 m²·s^-1数量级。本文提出的计算方法在数量级上与文献实验结果吻合,表明在大尺寸流化床且高流化风速下,颗粒横向扩散系数远大于小尺寸鼓泡流化床,为不同研究者实验结果的分歧提供了理论依据,也为预测大型流化床内颗粒扩散速率提供了放大策略。