Axial segregation in bubbling gas‐fluidized beds with Gaussian and lognormal distributions of Geldart Group B particles

Bubbling, gas‐fluidized bed experiments involving Geldart Group B particles with continuous‐size distributions have been carried out. Sand of various widths of Gaussian or lognormal distributions were completely fluidized, then axial concentration profiles were obtained from frozen‐bed sectioning. Similar to previous works on binary systems, results show that mean particle diameter decreases with increasing bed height, and that wider Gaussian distributions show increased segregation extents. Surprisingly, however, lognormal distributions exhibit a nonmonotonic segregation trend with respect to distribution widths. In addition, the shape of the local‐size distribution is largely preserved with respect to that of the overall distribution. These findings on the nature of local‐size distribution provide experimental confirmation of previous results for granular and gas‐solid simulations. Lastly, an interesting observation is that although monodisperse Geldart Group D particles cannot be completely fluidized, their presence in lognormal distributions investigated still results in complete fluidization of all particles. © 2010 American Institute of Chemical Engineers AIChE J, 2010     

Transient multi gas—solid reactions in a bubbling fluidized bed

A “multimodel” for gas‐solid reactions in a reacting particle has been applied to a bubbling fluidized bed reactor. The particle is tracked and bed and particle variables are determined continuously. The conservation equations of mass and heat with auxiliary relations are solved in an accelerating particle, which may rise or fall. The effects of bulk pressure, velocity and temperature, and particle diameter are studied. Heat and mass transfer coefficients may fluctuate up to 75% and 148% respectively. Doubling the pressure changes h_c by 75% and k_c by ?45%. Increase in pellet diameter reduces both h_c and k_c.     

Full‐physics simulations of spray‐particle interaction in a bubbling fluidized bed

Numerical simulations of a gas‐particle‐droplet system were performed using an Euler‐Lagrange approach. Models accounting for (1) the interaction between droplets and particles, (2) evaporation from the droplet spray, as well as (3) evaporation of liquid from the surface of non‐porous particles were considered. The implemented models were verified for a packed bed, as well as other standard flow configurations. The developed models were then applied for the simulation of flow, as well as heat and mass transfer in a fluidized bed with droplet injection. The relative importance of droplet evaporation vs. evaporation from the particle surface was quantified. It was proved that spray evaporation competes with droplet deposition and evaporation from the particle surface. Moreover, we show that adopting a suitable surface coverage model is vital when attempting to make accurate predictions of the particle's liquid content. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2569–2587, 2017     

Link between bubbling and segregation patterns in gas‐fluidized beds with continuous size distributions

Experiments involving a bubbling, gas‐fluidized bed with Gaussian and lognormal particle‐size distributions (PSDs) of Geldart Group B particles have been carried out, with a focus on bubble measurements. Previous work in the same systems indicated the degree of axial species segregation varies non‐monotonically with respect to the width of lognormal distributions. Given the widely accepted view of bubbles as “mixing agents,” the initial expectation was that bubble characteristics would be similarly non‐monotonic. Surprisingly, results show that measured bubble parameters (frequency, velocity, and chord length) increase monotonically with increasing width for all PSDs investigated. Closer inspection reveals a bubble‐less bottom region for the segregated systems, despite the bed being fully fluidized. More specifically, results indicate that, the larger the bubble‐less layer is, the more segregated the system becomes. The direct comparison between bubbling and segregation patterns performed provides a more complete physical picture of the link between the two phenomena. © 2011 American Institute of Chemical Engineers AIChE J, 2011     

大涡模拟方法模拟鼓泡床内气固两相流动

采用大涡模拟（LES）方法模拟气相湍流,颗粒动理学方法考虑颗粒相碰撞产生的动量和能量传递和耗散,采用颗粒相大涡模拟方法（LESp）模拟颗粒脉动导致的能量耗散,同时考虑介观尺度对颗粒相压力的影响,建立了气体-颗粒LES-θ-LESp双流体模型,研究鼓泡流化床内气固两相流动的特性。数值模拟与文献实测颗粒速度和实测颗粒浓度结果具有相同的变化趋势。     

Effect of mixing state of binary particles on bubble behavior in 2D fluidized bed

In this study, a thin 2D fluidized bed was used to investigate the effect of mixing state of the binary particles on bubble behavior through the analysis of images captured by a high-speed digital camera. Experimental results show that the mixing index increases gradually with increasing gas velocity and the binary particles are in different mixing states though they are in the steady fluidization state. The maximal bubble number is near the interface of the bed when the binary particles are in the segregation state, whereas the maximal bubble number is at the bottom when the binary particles are in the well mixing state. The small bubbles are position at the bottom and are adjacent to the bed wall, while the large bubbles are mainly located in the central regions of the bed. The average bubble diameter shows the different variation trends with the different mixing states of the binary particles. The correlations estimating bubble diameter according to the mixing state of the binary particles are developed, and the computing value agrees well with the experimental data.     

Filtered and heterogeneity‐based subgrid modifications for gas–solid drag and solid stresses in bubbling fluidized beds

Two different approaches to constitutive relations for filtered two‐fluid models (TFM) of gas–solid flows are deduced. The first model (Model A) is derived using systematically filtered results obtained from a highly resolved simulation of a bubbling fluidized bed. The second model (Model B) stems from the assumption of the formation of subgrid heterogeneities inside the suspension phase of fluidized beds. These approaches for the unresolved terms appearing in the filtered TFM are, then, substantiated by the corresponding filtered data. Furthermore, the presented models are verified in the case of the bubbling fluidized bed used to generate the fine grid data. The numerical results obtained on coarse grids demonstrate that the computed bed hydrodynamics is in fairly good agreement with the highly resolved simulation. The results further show that the contribution from the unresolved frictional stresses is required to correctly predict the bubble rise velocity using coarse grids. © 2013 American Institute of Chemical Engineers AIChE J, 60: 839–854, 2014     

Comparative analysis of subgrid drag modifications for dense gas‐particle flows in bubbling fluidized beds

Many subgrid drag modifications have been put forth to account for the effect of small unresolved scales on the resolved mesoscales in dense gas‐particle flows. These subgrid drag modifications significantly differ in terms of their dependencies on the void fraction and the particle slip velocity. We, therefore, compare the hydrodynamics of a three‐dimensional bubbling fluidized bed computed on a coarse grid using the drag correlations of the groups of (i) EMMS, (ii) Kuipers, (iii) Sundaresan, (iv) Simonin, and the homogenous drag law of (v) Wen and Yu with fine grid simulations for two different superficial gas velocities. Furthermore, we present an (vi) alternative approach, which distinguishes between resolved and unresolved particle clusters revealing a grid and slip velocity dependent heterogeneity index. Numerical results are analyzed with respect to the time‐averaged solids volume fraction and its standard deviation, gas and solid flow patterns, bubble size, number density, and rise velocities. © 2013 American Institute of Chemical Engineers AIChE J, 59: 4077–4099, 2013     

Effect of ring‐type internals on solids distribution in a dual circulating fluidized bed system—cold flow model study

The redistribution of solids in a counter‐current circulating fluidized bed (CFB) by effect of ring‐type internals was investigated in a downscaled cold‐flow model. The system consists of two interconnected CFB reactors, in which the primary reactor operates like a common riser while the secondary reactor operates in counter‐current. The unit works without circulation rate control devices and the inventory splits inherently between the two reactors by pressure balance and depending on the fluidization velocities. Previous studies showed an increment in the total pressure drop in the secondary reactor as result of the internals installation. With the purpose of obtaining comparable inventory in the secondary reactor with and without rings, a device for adjustment of total inventory was designed and installed. Effects of the aperture ratio, number of rings, fluidization velocity, and particles circulation rate were studied. The results obtained approach a guideline for the detailed design of similar configurations. © 2013 American Institute of Chemical Engineers AIChE J, 59: 3612–3623, 2013     

Borescopy in pressurized gas‐solid fluidized beds

A borescopic technique was used for finding the effect of pressure on the hydrodynamics of gas‐solid fluidized beds. The results showed that solids radial distribution may become more or less uniform with increasing pressure depending on the superficial gas velocity. Moreover, it is found that the solids volume fraction of the emulsion phase may decrease at relatively high pressures, only in the central region of the bed. Additionally, it is observed that with increasing pressure the bubble size generally decreased in the central regions and increased near the wall regions. This trend was more complicated at low excess gas velocities. The number of bubbles increased for the central regions and near the walls for all the performed experiments. However, this parameter showed a different trend at other radial positions. © 2018 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J, 64: 3303–3311, 2018     

Radial Liquid Dispersion and Bubble Distribution in Three‐Phase Circulating Fluidized Beds

The liquid dispersion and bubble distribution in the radial direction have been investigated in the riser of a three‐phase circulating fluidized bed whose diameter is 0.102m and 3.5m in height. Effects of gas and liquid velocities and solid circulation rate have been determined. It has been found that the radial distribution of bubbles is related closely to the liquid dispersion in the radial direction. The size and rising velocity of bubbles tend to increase as the radial position approaches to the center of the riser. The bubble size increases with increasing U_G, but it decreases with increasing U_L or G_S in all radial positions. The radial dispersion coefficient of the liquid phase increases with increasing U_G or G_S, however, it tends to decrease with increasing U_L. The value of Dr has been well correlated in terms of dimensionless groups based on the isotropic turbulence model.     

Bubble‐separation dynamics in a planar cyclone: Experiments and CFD simulations

A planar cyclone is designed for visualizing bubbles in the cross‐section of a degassing hydrocyclone. The pressure distribution is studied through a series of experiments and Reynolds stress model simulations. The velocity distribution of the planar cyclone mostly exhibits the quasi‐forced vortex zone and boundary layer zone. The bubble dynamics are simulated using both Euler‐Euler and Euler‐Lagrange approaches, and the output is compared with the imaging results. The Euler‐Euler simulation provides more accurate predictions of the bubble trajectory. The histograms of residence time and traveling distance given by the Euler‐Lagrange approach exhibit a reasonably regular pattern. With higher values of the inlet Reynolds number, stronger forces acting on the bubbles lead to a decreased but more uniform residence time. © 2018 American Institute of Chemical Engineers AIChE J, 64: 2689–2701, 2018     

Experimental study of hydrodynamics and thermal behavior of a pseudo‐2D spout‐fluidized bed with liquid injection

A novel nonintrusive technique is presented to investigate hydrodynamic and thermal behavior of gas–solid spout‐fluidized beds with liquid injection, by simultaneously capturing visual and infrared images. Experiments were performed in a pseudo‐2D bed with draft plates filled with glass or γ‐alumina particles to investigate the effect of liquid injection and particle properties on the flow characteristics. For the glass particles under dry and wet conditions, time‐averaged particle velocities show similar quasi‐steady‐state behavior. However, under wet conditions, lower particle velocities were observed in both spout and annulus as compared with the dry system. Whereas, γ‐alumina particles do not show considerable variation in the particle velocities under dry and wet conditions and fluidize well at higher liquid injection rates. Additionally, for the glass particles, the particle temperature significantly decreases as compared to the γ‐alumina particles. © 2015 American Institute of Chemical Engineers AIChE J, 61: 1146–1159, 2015     

Bubble dynamics in a 3‐D gas–solid fluidized bed using ultrafast electron beam X‐ray tomography and two‐fluid model

Bubble characteristics in a three‐dimension gas‐fluidized bed (FB) have been measured using noninvasive ultrafast electron beam X‐ray tomography. The measurements are compared with predictions by a two‐fluid model (TFM) based on kinetic theory of granular flow. The effect of bed material (glass, alumina, and low linear density polyethylene (LLDPE), d_p ～1 mm), inlet gas velocity, and initial particle bed height on the bubble behavior is investigated in a cylindrical column of 0.1‐m diameter. The bubble rise velocity is determined by cross correlation of images from dual horizontal planes. The bubble characteristics depend highly upon the particle collisional properties. The bubble sizes obtained from experiments and simulations show good agreement. The LLDPE particles show high gas hold‐up and higher bubble rise velocity than predicted on basis of literature correlations. The bed expansion is relatively high for LLDPE particles. The X‐ray tomography and TFM results provide in‐depth understanding of bubble behavior in FBs containing different granular material types. © 2014 American Institute of Chemical Engineers AIChE J, 60: 1632–1644, 2014     

On the accuracy of Landweber and Tikhonov reconstruction techniques in gas‐solid fluidized bed applications

As electrical capacitance tomography technique needs a sophisticated reconstruction, the accuracy of two of the most widely used reconstruction techniques (Landweber and Tikhonov) for gas‐fluidized bed applications were assessed. For this purpose, the results of two‐fluid model simulations were used as an input of reconstruction. After finding the optimum reconstruction parameters for the studied system, it is found that both techniques were able to obtain the radial profile and overall value of average volume fraction very well. Conversely, both methods were incapable to determine bubble sizes accurately especially small bubble sizes, unless the Landweber technique with inverted Maxwell concentration model is applied. The probability distribution of the reconstructed results was also smoother in transition between the emulsion and bubble phases compared to the reality. Finally, no significant differences in noise immunity of these two techniques were observed. © 2015 American Institute of Chemical Engineers AIChE J, 61: 4102–4113, 2015     

Modeling and control of solar‐grade silicon production in a fluidized bed reactor

A multiscale model predicts silicon production yield and powder loss in a fluidized bed reactor for solar silicon production. The reaction module calculates the silicon vapor deposition and powder scavenging rates. A computational fluid dynamics model predicts temperature and bed density. A population balance model calculates the particle‐mass distribution functions on silicon yield. The model results are validated against industrial data. Furthermore, we conduct a sensitivity analysis to investigate the effect of gas flow rate and inlet silane concentration. Finally, a control strategy is proposed to maintain the process at the desired operating point. © 2014 American Institute of Chemical Engineers AIChE J, 60: 1740–1751, 2014     

Realization and control of multiple temperature zones in liquid‐containing gas–solid fluidized bed reactor

Fluidized bed reactors (FBRs) have been developed to establish multiple temperature zones for various industrial processes. To overcome the common weakness, this work proposed to spray liquid into bottom and upper zones, respectively, to realize multiple temperature zones FBR (MTZFBR). Temperature, pressure, and acoustic emission techniques were applied to fully characterize liquid interaction and hydrodynamics. Compared with the bottom liquid‐spraying approach, the upper liquid‐spraying approach showed higher temperature difference (ΔT) and better fluidization stability, thus was selected for further control studies. Effects of liquid flow rate, static bed height, and inlet gas temperature on MTZFBR were studied systematically. The results showed that increasing liquid evaporation behavior or decreasing liquid bridge behavior enhance ΔT and fluidization stability and vice versa. G–L–S fluidization pattern depended mostly on the liquid behaviors and fluidization stability, and thus the stabilized MTZFBR could be regarded as a coexisted mode of two distinctive G–L–S fluidization patterns. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1454–1466, 2016     

Dependence of particle fluctuation velocity on gas flow, and particle diameter in gas fluidized beds for monodispersed spheres in the Geldart B and A fluidization regimes

In this paper we present new experimental data on the steady-state, mean squared, fluctuation velocity, or granular temperature, of Geldart B polymer, glass, nickel, and stainless steel monodispersed spheres averaged over the wall of a gas fluidized bed, as a function of gas flow and sphere diameter. The granular temperature is obtained by Acoustic Shot Noise technology—namely power spectral analysis of the steady state vibrational energy of the wall excited by random sphere impact, and calibrated by hammer excitation over the wall. The new data extends to polymer and metallic spheres the experimental discovery of a 1996 paper of Cody et al. that the fluctuation velocity of Geldart B glass spheres when scaled to the gas superficial velocity, U_s, is inversely proportional to sphere diameter, directly proportional to a fundamental length scale, D_o^B, and is a universal function of U = (U_s / U_mf). We also demonstrate that the new data is consistent with the diameter dependence of the fluctuation velocity that can be derived from both the 1997 paper of Menon and Durian, who measured random sphere motion near the wall through the spectroscopy of scattered laser light, and the 1992 paper of Rahman and Campbell, who measured the average granular pressure of random sphere impact on a porous steel membrane. While the inverse scaling of the fluctuation velocity with sphere diameter, and the existence of a fundamental length scale for gas fluidization, D_o^B, had not been a feature of any published fundamental model, or computer simulation, of the steady state granular temperature of spheres in gas fluidized beds, we show that it is a feature of two recent dense kinetic fluidization models published in 1999, by Buyevich and Kapbasov, and Koch and Sangani. Both theories implicitly define a fundamental length scale for the fluctuation velocity, D^? = (μ_f² / ρ_p²g)^1 / 3, where ρ_p is the sphere density, μ_f is the gas viscosity, and g is the laboratory gravitational field. The new data for polymer, glass, nickel and stainless steel spheres presented in this paper, defines D_o^B = (56 ± 2)D^?. We use the Anderson-Jackson stability model to show that the length scale D_o^B, also defines a stability length scale, such that for D < D_o^B(D > D_o^B), the uniform dense phase of the fluidized bed is stable (unstable), against one dimensional, first order fluctuations in sphere concentration. The length scale, D_o^B is thus the theoretical equivalent to the empirical scaling length introduced by Geldart, D_B/A, to distinguish spheres (D > D_B/A) that bubble at fluidization, from spheres (D < D_B/A) that fluidize before bubbling. Finally, we present new experimental data, on the remarkable changes in the granular temperature, bed expansion, and bed collapse time, between Geldart B and Geldart A monodispersed glass spheres, and compare that data to granular temperature, and bed expansion, for Geldart A rough, non-spherical, log-normal dispersed diameter catalytic particles.