Bubble behavior in corrugated‐wall bubbling fluidized beds—Experiments and CFD simulations

A new concept to harness bubble dynamics in bubbling fluidization of Geldart D particles was proposed. Various geometrical declinations of a cold‐prototype corrugated‐wall bubbling fluidized bed were compared at different flow rates (U_g) to conventional flat‐wall fluidized bed using high‐speed digital image analysis. Hydrodynamic studies were carried out to appraise the effect of triangular‐shaped wall corrugation on incipient fluidization, bubble coalescence (size and frequency), bubble rise velocity, and pressure drop. Bubble size and rise velocity in corrugated‐wall beds were appreciably lower, at given U_g/U_mb, than in flat‐wall beds with equal flow cross‐sectional areas and initial bed heights. The decrease (increase) in size (frequency) of bubbles during their rise was sustained by their periodic breakups while protruding through the necks between corrugated plates. Euler‐Euler transient full three‐dimensional computational fluid dynamic simulations helped shape an understanding of the impact of corrugation geometry on lowering the minimum bubbling fluidization and improving gas distribution. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Experimental and numerical study on a bubbling fluidized bed with wet particles

As liquid bridge between particles acts an important role in the particle system, it is of considerable significance to analyze the flow hydrodynamics of wet particles in fluidized beds, which will improve the reactor design and process optimization. Thus, experimental and numerical investigations on wet particles in a bubbling fluidized bed are conducted in current work. On experimental side, particle image velocimetry (PIV) technology is employed with a designed bubbling fluidized bed. The silicone oil is used in this work because it is nonvolatile and transparent. On numerical side, a modified discrete element method (DEM) numerical method is developed by compositing an additional liquid‐bridge module into the traditional soft‐sphere interaction model. Most of the physical parameters are chosen to correspond to the experimental settings. Good agreements of particle velocity are found between the DEM simulation and PIV measurement. The performance of different liquid contents and superficial gas velocities are examined. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1970–1985, 2016     

Verification of fluidized bed electrical capacitance tomography measurements with a fibre optic probe

Previous studies aimed at determining the spatial accuracy of electrical capacitance tomography (ECT) have employed phantoms placed within the ECT measurement space. No previous studies have compared ECT with a second independent measurement technique in an operating fluidized bed. In the present work, radial voidage profiles have been measured with ECT in the 0.14-m I.D. riser of a circulating fluidized bed (CFB) and in a bubbling fluidized bed with a 0.19-m I.D. The dynamic and time-averaged radial voidage profiles have been compared with measurements taken with a fibre optic probe in the same riser and in a slightly narrower (0.15-m I.D.) bubbling fluidized bed. In spite of the intrusiveness of the latter technique, the time-averaged radial profiles in the CFB riser fall within 10% of each other when the CFB is operated at high-flux conditions that lead to a very dense wall region. Iterative reconstruction of the ECT images is not needed in this case. Similar agreement is found between the two techniques in the bubbling fluidized bed, but off-line iterative image reconstruction is clearly necessary in this fluidization regime. These results suggest that ECT, which is often described as a tomographic imaging technique with low spatial resolution, can in fact provide semi-quantitative time-averaged images of the cross-section of fluidized beds of diameter comparable to or less than that used here.     

Particle‐resolved direct numerical simulation of gas–solid dynamics in experimental fluidized beds

Particle‐resolved direct numerical simulations (PR‐DNS) of a simplified experimental shallow fluidized bed and a laboratory bubbling fluidized bed are performed by using immersed boundary method coupled with a soft‐sphere model. Detailed information on gas flow and individual particles’ motion are obtained and analyzed to study the gas–solid dynamics. For the shallow bed, the successful predictions of particle coherent oscillation and bed expansion and contraction indicate all scales of motion in the flow are well captured by the PD‐DNS. For the bubbling bed, the PR‐DNS predicted time averaged particle velocities show a better agreement with experimental measurements than those of the computational fluid dynamics coupled with discrete element models (CFD‐DEM), which further validates the predictive capability of the developed PR‐DNS. Analysis of the PR‐DNS drag force shows that the prevailing CFD‐DEM drag correlations underestimate the particle drag force in fluidized beds. The particle mobility effect on drag correlation needs further investigation. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1917–1932, 2016     

Effective particle diameters for simulating fluidization of non‐spherical particles: CFD‐DEM models vs. MRI measurements

Computational fluid dynamics—discrete element method (CFD‐DEM) simulations were conducted and compared with magnetic resonance imaging (MRI) measurements (Boyce, Rice, and Ozel et al., Phys Rev Fluids. 2016;1(7):074201) of gas and particle motion in a three‐dimensional cylindrical bubbling fluidized bed. Experimental particles had a kidney‐bean‐like shape, while particles were simulated as being spherical; to account for non‐sphericity, “effective” diameters were introduced to calculate drag and void fraction, such that the void fraction at minimum fluidization (ε_mf) and the minimum fluidization velocity (U_mf) in the simulations matched experimental values. With the use of effective diameters, similar bubbling patterns were seen in experiments and simulations, and the simulation predictions matched measurements of average gas and particle velocity in bubbling and emulsion regions low in the bed. Simulations which did not employ effective diameters were found to produce vastly different bubbling patterns when different drag laws were used. Both MRI results and CFD‐DEM simulations agreed with classic analytical theory for gas flow and bubble motion in bubbling fluidized beds. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2555–2568, 2017     

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     

Fictitious particle method: A numerical model for flows including dense solids with large size difference

Large solids coexist with small solids in a number of dense gas‐solid flow applications such as fluidized beds and pneumatic conveyers. A new numerical model that is based on the discrete element method–computational fluid dynamics mesoscopic model and extended by introducing an idea appearing in volume penalization method is presented. In computational cells including large and small solids, the amount of momentum exchange between the fluid and the solids is estimated by assuming that a large solid consist of small, dense fictitious particles. We describe the proposed model in detail and show the optimal model parameters found through a number of parameter‐dependency studies. Validation study is performed for the motion of a large sphere in a bubbling fluidized bed and good agreements are confirmed for floating and sinking motions of the sphere between the present model and the experiment. © 2014 American Institute of Chemical Engineers AIChE J, 60: 1606–1620, 2014     

Using Electrical Capacitance Tomography to Investigate Gas Solid Contacting

An ECT system has been shown to be useful in observing differences in particle behaviour in a bubbling fluidised bed. An image analysis technique is further described that utilises the construction of solids concentration profiles at three key values, x = 0.2, x = 0.5 and x = 0.8, which have been identified as important solids concentrations in a new Bubble Structure Model. Significant differences in the bubble structure are shown to result when the concentration of fine particles is slightly increased in a bubbling fluidised bed. Changes of this type would seriously alter the gas solid contact efficiency in the fluidised bed, which would significantly influence selectivity of in‐bed catalytic reactions and gas residence time.     

Evaporative liquid jets in gas–liquid–solid flow system

Some aspects of the fundamental characteristics of evaporative liquid jets in gas–liquid–solid flows are studied and some pertinent literature is reviewed. Specifically, two conditions for the solids concentration in the flow are considered, including the dilute phase condition as in pneumatic convey and the dense phase condition as in bubbling or turbulent fluidized beds. Comparisons of the fundamental behavior are made of the gas–solid flow with dispersed non-evaporative as well as with evaporative liquids.For dilute phase conditions, experiments and analyses are conducted to examine the individual phase motion and boundaries of the evaporative region and the jet. Effects of the solids loading and heat capacity, system temperature, gas flow velocity and liquid injection angle on the jet behavior in gas and gas–solid flows are discussed. For dense phase conditions, experiments are conducted to examine the minimum fluidization velocity and solids distribution across the bed under various gases and liquid flow velocities. The electric capacitance tomography is developed for the first time for three-phase real time imaging of the dense gas–solid flow with evaporative liquid jets. The images reflect significantly varied bubbling phenomenon compared to those in gas–solid fluidized beds without evaporative liquid jets.     

Heat Transfer from Immersed Vertical Cylinders in Gas‐Liquid and Gas‐Liquid‐Solid Fluidized Beds

The heat transfer coefficient, h, was measured using a cylindrical heater vertically immersed in liquid‐solid and gas‐liquid‐solid fluidized beds. The gas used was air and the liquids used were water and 0.7 and 1.5 wt‐% carboxymethylcellulose (CMC) aqueous solutions. The fluidized particles were sieved glass beads with 0.25, 0.5, 1.1, 2.6, and 5.2 mm average diameters. We tried to obtain unified dimensionless correlations for the cylinder surface‐to‐liquid heat transfer coefficients in the liquid‐solid and gas‐liquid‐solid fluidized beds. In the first approach, the heat transfer coefficients were successfully correlated in a unified formula in terms of a modified j_H‐factor and the modified liquid Reynolds number considering the effect of spatial expansion for the fluidized bed within an error of 36.1 %. In the second approach, the heat transfer coefficients were also correlated in a unified formula in terms of the dimensionless quantities, Nu/Pr^1/3, and the specific power group including energy dissipation rate per unit mass of liquid, E^1/3D^4/3/ν_l, within a smaller error of 24.7 %. It is also confirmed that a good analogy exists between the surface‐to‐liquid heat transfer and mass transfer on the immersed cylinder in the liquid‐solid and gas‐liquid‐solid fluidization systems.     

Bubble Splitting in a Pseudo‐2D Gas‐Solid Fluidized Bed for Geldart B‐Type Particles

Bubble splitting in 2D gas‐solid freely bubbling fluidized beds is experimentally investigated using digital image analysis. The quantitative results can be applied for the development of a new breakage model for bubbly fluidized beds, especially discrete bubble models. The variation of splitting frequency with bubble diameter, new resulting bubble volumes, positions, and also the assumptions of mass and momentum conservation for bubbles after breakage are studied in detail. Small bubbles are found to be more stable than large ones and nearly all mother bubbles split into two almost equally sized daughter bubbles. The momentum of gas bubbles in the vertical direction remains approximately constant after breakage, whereas that of bubbles in the horizontal direction changes with no clear trend. The effect of fluidizing gas velocity in breakage frequency is also examined.     

Staugrenze bei der Begasung mit Lochböden

In gas‐liquid contacting with perforated plates a certain gas velocity exists, above which the weeping rate will decrease noticeably. This loading point was identified as the transition between the bubbling and jetting regime, based on a physical model derived in previous studies. The same criterion is valid for gas‐solid fluidized beds.     

Validation study on spatially averaged two‐fluid model for gas‐solid flows: II. Application to risers and fluidized beds

In our prior study (Schneiderbauer, AIChE J. 2017;63(8):3544–3562), a spatially averaged two‐fluid model (SA‐TFM) was presented, where closure models for the unresolved terms were derived. These closures require constitutive relations for the turbulent kinetic energies of the gas and solids phase as well as for the subfilter variance of the solids volume fraction. We had ascertained that the filtered model do yield nearly the same time‐averaged macroscale flow behavior in bubbling fluidized beds as the underlying kinetic‐theory‐based two‐fluid model, thus verifying the SA‐TFM model approach. In the present study, a set of 3D computational simulations for validation of the SA‐TFM against the experimental data on riser flow and bubbling fluidized beds is performed. Finally, the SA‐TFM predictions are in fairly good agreement with experimental data in the case of Geldart A and B particles even though using very coarse grids. © 2018 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J, 64: 1606–1617, 2018     

Assessment of mesoscale solid stress in coarse‐grid TFM simulation of Geldart A particles in all fluidization regimes

This study focused on assessing the effect of mesoscale solid stress in the coarse grid two‐fluid model (TFM) simulation of gas–solid fluidized beds of Geldart Group A particles over a broad range of flow regimes, including bubbling, turbulent, fast, and pneumatic transport fluidization regimes. Particularly, the impact of mesoscale solid pressure, mesoscale solid viscosity, and mesoscale solid stress anisotropy were investigated by comparing six different coarse‐grid TFM settings. Compared with the available experimental data, it is found that both the kinetic theory‐based TFM with only drag correction and the filtered TFM can predict the flow behavior in all fluidization regimes. Mesoscale solid pressure and viscosity have the opposite impact on flow hydrodynamics; they compete and offset each other, which confirms the assumption employed in many previous studies that the mesoscale solid stress could be neglected in coarse‐grid TFM simulation. Published 2018. This article is a U.S. Government work and is in the public domain in the USA. AIChE J, 64: 3565–3581, 2018     

Heat transfer within dynamically structured bubbling fluidized beds subject to vibration: A two-fluid modeling study

Bubbling fluidized beds are often used to achieve a uniform particle temperature distribution in industrial processes involving gas and particles. However, the chaotic bubble dynamics pose significant challenges in scale-up. Recent work (Guo et al., 2021, PNAS 118, e2108647118) has shown that using vibration can structure the bubbling pattern to a highly predictable manner with the characteristic bubble properties independent of system width, opening opportunities to address key issues associated with conventional bubbling fluidized beds. Herein, using two-fluid modeling simulations, we studied heat transfer characteristics within the dynamically structured bubbling fluidized bed and compared to unstructured bubbling fluidized beds and packed beds. Simulations show that the structured bubbling fluidized bed can achieve the most uniform particle temperature distribution because it can achieve the best particle mixing while maintaining a global heat transfer coefficient similar to that of a freely bubbling fluidized bed.     

Development of novel two-interconnected fluidized bed system

Two-interconnected fluidized bed systems are widely used in various processes such as Fisher-Tropsch, hot gas desulfurization, CO₂ capture-regeneration with dry sorbent, chemical-looping combustion, sorption-enhanced steam methane reforming, chemical-looping hydrogen generation system, and so on. However, conventional two-interconnected fluidized bed systems are very complex, large, and difficult to operate because most of these systems require a riser and/or pneumatic transport line for solid conveying and loopseals or seal-pots for gas sealing, recirculation of solids to the riser, and maintaining of pressure balance. To solve these problems, a novel two-interconnected fluidized bed system has been developed. This system has two bubbling beds, solid injection nozzles, solid conveying lines, and downcomers. In this study, the effects of operating variables on solid circulation rate and gas leakage between two beds have been investigated in a cold mode two-interconnected fluidized bed system. The solid circulation rate increased as the hole diameter on the injection nozzle, the diameter of the injection nozzle, the solid height above the holes, and the number of holes on the injection nozzle increased. The gas leakage between the beds was negligible. Moreover, long-term operation of continuous solid circulation up to 60 hours was performed to check the feasibility of stable operation. The pressure drop profiles in the system loop were maintained steadily and solid circulation was smooth and stable.     

Granular Temperature in Bubbling Fluidized Beds

In the CFD simulation of industrial‐scale bubbling fluidized beds, coarse grids are necessary due to limited computation resources, which creates a problem about how to characterize the effect of meso‐scale bubbling structures on the particulate phase stresses and inter‐phase interaction force. In this paper, the two‐phase theory of bubbling fluidization is applied in the analysis of the fluctuation characteristics of solid volume fractions and then it was used to establish a simple model for granular temperature in bubbling fluidized beds, which is based on the local pseudo‐thermal energy balance. It was shown that the granular temperature obtained experimentally is well predicted by the present model.     

Characteristics of a Pressurized Gas‐Solid Magnetically Fluidized Bed

The hydrodynamic, heat and mass transfer characteristics of a pressurized co‐current gas‐solid magnetically fluidized bed (MFB) were systematically investigated considering major influence factors, such as magnetic field strength, superficial gas velocity, and operating pressure. It was shown that this pressurized gas‐solid MFB has the advantages of a wider operation range of the superficial gas velocity under bubble‐free particulate fluidization, a larger bed voidage with smaller pressure drop across the bed, and larger heat transfer efficiency, compared with a conventional fluidized bed. Moreover, the minimum bubbling velocity, gas‐solid mass, and heat transfer coefficients were correlated at high accuracy within the investigated range of operating conditions.     

Effect of bed size on hydrodynamics in 3‐D gas–solid fluidized beds

It is well known that hydrodynamics observed in large scale gas–solid fluidized beds are different from those observed in smaller scale beds. In this article, an efficient two‐fluid model based on kinetic theory of granular flow is applied, with the goal to highlight and investigate hydrodynamics differences between three‐dimensional fluidized beds of diameter 0.10, 0.15, 0.30, 0.60, and 1.0 m, focusing on the bubble and solids flow characteristics in the bubbling regime. Results for the 0.30 m diameter bed are compared with experimental results from the literature. The bubble size evolution closely follows a correlation proposed by Werther for small beds, and a correlation proposed by Darton for sufficiently large beds. The bubble size increases as the bed diameter is increased from 0.10 to 0.30 m, and remains approximately constant for bed diameters from 0.30 to 1.0 m. Concurrently, an increase in bubble rise velocity is observed, with a much high bubble rise velocity in the largest bed of diameter 1.0 m due to gulf stream circulations. The dynamics in shallow and deep beds is predicted to be different, with marked differences in bubble size and solids circulation patterns. © 2015 American Institute of Chemical Engineers AIChE J, 61: 1492–1506, 2015