CFD–DEM study of effect of bed thickness for bubbling fluidized beds

The effect of bed thickness in rectangular fluidized beds is investigated through the CFD–DEM simulations of small-scale systems. Numerical results are compared for bubbling fluidized beds of various bed thicknesses with respect to particle packing, bed expansion, bubble behavior, solids velocities, and particle kinetic energy. Good two-dimensional (2D) flow behavior is observed in the bed having a thickness of up to 20 particle diameters. However, a strong three-dimensional (3D) flow behavior is observed in beds with a thickness of 40 particle diameters, indicating the transition from 2D flow to 3D flow within the range of 20–40 particle diameters. Comparison of velocity profiles near the walls and at the center of the bed shows significant impact of the front and back walls on the flow hydrodynamics of pseudo-2D fluidized beds. Hence, for quantitative comparison with experiments in pseudo-2D columns, the effect of walls has to be accounted for in numerical simulations.     

Time-series analysis of pressure fluctuations in gas–solid fluidized beds – A review

This work reviews methods for time-series analysis for characterization of the dynamics of gas–solid fluidized beds from in-bed pressure measurements for different fluidization regimes. The paper covers analysis in time domain, frequency domain, and in state space. It is a follow-up and an update of a similar review paper written a decade ago. We use the same pressure time-series as used by Johnsson et al. (2000). The paper updates the previous review and includes additional methods for time-series analysis, which have been proposed to investigate dynamics of gas–solid fluidized beds. Results and underlying assumptions of the methods are discussed.     

Characterization of gas–solid fluidized bed hydrodynamics by vibration signature analysis

Vibration measurement, as a non-intrusive technique, was used to characterize the hydrodynamics of fluidized beds. A series of experiments were performed in a lab-scale fluidized bed using two accelerometers for measuring the vibration of the bed and a pressure probe for measuring pressure fluctuations. The output signals were analyzed by statistical methods. The results show that the vibration technique can predict transition velocities at high velocities and indicate that analyzing the vibration signals can be an effective non-intrusive technique to characterize the hydrodynamics of fluidized beds. It was shown that transition from bubbling to turbulent velocity can be determined from the variation of standard deviation and kurtosis of vibration signals against superficial gas velocity of the bed. However, this point could be determined only from standard deviation of pressure fluctuations, and not from skewness or kurtosis of pressure fluctuations.     

A new drag model for TFM simulation of gas–solid bubbling fluidized beds with Geldart-B particles

In this work, a new drag model for TFM simulation in gas–solid bubbling fluidized beds was proposed, and a set of equations was derived to determine the meso-scale structural parameters to calculate the drag characteristics of Geldart-B particles under low gas velocities. In the new model, the meso-scale structure was characterized while accounting for the bubble and meso-scale structure effects on the drag coefficient. The Fluent software, incorporating the new drag model, was used to simulate the fluidization behavior. Experiments were performed in a Plexiglas cylindrical fluidized bed consisting of quartz sand as the solid phase and ambient air as the gas phase. Comparisons based on the solids hold-up inside the fluidized bed at different superficial gas velocities, were made between the 2D Cartesian simulations, and the experimental data, showing that the results of the new drag model reached much better agreement with experimental data than those of the Gidaspow drag model did.     

Dynamics of fine particles in liquid-solid fluidized beds

On the basis of the Local Equilibrium Model （LEM）, fine particles with large Richardson-Zaki exponent n show, under certain conditions during bed expansion and collapse, different dynamic behavior from particles with small n. For an expansion process there may be a concentration discontinuity propagating upward from the distributor, and, on the contrary, for a collapse process there may be a progressively broadening and upward-propagating continuous transition zone instead of discontinuity. The predictions of the bed height variation and the discontinuity trace have been validated experimentally.     

Two and three dimensional modeling of fluidized bed with multiple jets in a DEM–CFD framework

Fluidized beds with multiple jets have widespread industrial applications. The objective of this paper is to investigate the jet interactions and hydrodynamics of a fluidized bed with multiple jets. Discrete element modeling coupled with in-house CFD code GenIDLEST has been used to simulate a bed with nine jets. The results are compared with published experiments. Mono dispersed particles of size 550 μm are used with 1.4 times the minimum fluidization velocity of the particles. Both two and three dimensional computations have been performed. To the best of our knowledge, the results presented in this paper are the first full 3D simulations of a fluidized bed performed with multiple jets. Discrepancies between the experiment and simulations are discussed in the context of the dimensionality of the simulations. The 2D solid fraction profile compares well with the experiment close to the distributor plate. At higher heights, the 2D simulation over-predicts the solid fraction profiles near the walls. The 3D simulation on the other hand is better able to capture the solid fraction profile higher up in the bed compared to that near the distributor plate. Similarly, the normalized particle velocities and the particle fluxes compare well with the experiment closer to the distributor plate for the 2D simulation and the freeboard for the 3D simulation, respectively. A lower expanded bed height is predicted in the 2D simulation compared to the 3D simulation and the experiment. The results obtained from DEM computations show that a 2D simulation can be used to capture essential jetting trends near the distributor plate regions, whereas a full scale 3D simulation is needed to capture the bubbles near the freeboard regions. These serve as validations for the experiment and help us understand the complex jet interaction and solid circulation patterns in a multiple jet fluidized bed system.     

Influence of vertical heat exchanger tubes,their arrangement and the column diameter on the hydrodynamics in a gas–solid bubbling fluidized bed

The hydrodynamic behavior of a cold-flow gas–solid fluidized bed with an inner diameter of 22 cm is investigated by means of an ultra-fast X-ray tomographic setup. In the case of an exothermal reaction, heat exchanger tubes are required to remove the reaction heat out of the bubbling fluidized bed reactor. For the examined cold-flow model, the heat exchanger tubes are replaced by vertical internals that serve as placeholder. The influence of vertical internals on the bubble properties for different spatial configurations (square and circular arrangements) is investigated in addition to measurements without internals. Furthermore, the hydrodynamic results of the Ø 22 cm column are compared with an available data set which is based on measurements that were conducted in a column with an inner diameter of 14 cm. The objective of this paper is to provide measurement data for the scale-up process as well as for various computer models simulating a bubbling fluidized bed with vertical internals. It was found that the scale-up process from pilot plants to an industrial scale may be simplified if vertical internals are present, independently of the geometric arrangement.     

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.     

Particle shape effect on hydrodynamics and heat transfer in spouted bed: A CFD–DEM study

Spouted bed has drawn much attention due to its good heat and mass transfer efficiency in many chemical units. Investigating the flow patterns and heat and mass transfer inside a spouted bed can help optimize the spouting process. Therefore, in this study, the effects of particle shape on the hydrodynamics and heat transfer in a spouted bed are investigated. This is done by using a validated computational fluid dynamics–discrete element method (CFD–DEM) model, considering volume–equivalent spheres and oblate and prolate spheroids. The results are analysed in detail in terms of the flow pattern, microstructure, and heat transfer characteristics. The numerical results show that the prolate spheroids (Ar = 2.4) form the largest bubble from the beginning of the spouting process and rise the highest because the fluid drag forces can overcome the interlocking and particle–particle frictional forces. Compared with spherical particles, ellipsoidal spheroids have better mobility because of the stronger rotational kinetic energy resulting from the rough surfaces and nonuniform torques. In addition, the oblate spheroid system exhibits better heat transfer performance benefiting from the larger surface area, while prolate spheroids have poor heat transfer efficiency because of their orientation distribution. These findings can serve as a reference for optimizing the design and operation of complex spouted beds.     

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.     

Improvement,validation and application of CFD/DEM model to dense gas–solid flow in a fluidized bed

Dense gas–solid flow with solid volume fraction greater than 10% and at moderate Reynolds number is important in many industrial facilities such as fluidized beds. In this work, the Euler–Lagrange approach in combination with a deterministic collision model is applied to a laboratory-scale fluidized bed. The fluid–particle interaction is studied using a new procedure called the offset method, which results in several numbers of spatial displacements of the fluid grid. The proposed method is highly precise in determining porosity and momentum transfer, thus improving simulation accuracy. A validation study was carried out to assess the results using this in-house CFD/DEM code against 5-s operation of a Plexiglas spouted-fluidized bed, showing good qualitative correlation of solid distribution in the bed and acceptable quantitative agreement of pressure drops at different positions in the bed. In view of high computing cost, special emphasis is placed on effective program design, such as application of advanced detection algorithm for particle–particle/wall collisions, the multi-grid method and parallel calculation. In this context, the influence of increasing the processor number, up to 36, on calculation efficiency was investigated.     

Experimental study on reactor performance of gas–solids low-velocity fluidized beds

Reactor performance of bubbling fluidized bed (BFB) and turbulent fluidized bed (TFB) was carefully examined and systematically compared using catalytic ozone decomposition as a model reaction, based on a complete mapping of local flow structures and spatial distributions of ozone conversion and solids holdup. TFB clearly has a higher conversion and shows better reactor performance than BFB as a result of the vigorously turbulent flow and the relatively homogeneous gas–solids mixing in TEB. Besides, the intensive interaction between gas and solids in TFB leads to greater gas–solids contact efficiency of TFB over that of BFB. Due to gas bypassing and backmixing caused by bubbling behaviours and two-phase structure, BFB deviates significantly from a plug flow reactor and sometimes from a continuously stirred tank reactor. The flow structures essentially dictate the reactor performance in the low-velocity fluidized beds.     

CFD–DEM simulation of turbulence modulation in horizontal pneumatic conveying

A study is presented to evaluate the capabilities of the standard k–ε turbulence model and the k–ε turbulence model with added source terms in predicting the experimentally measured turbulence modulation due to the presence of particles in horizontal pneumatic conveying, in the context of a CFD–DEM Eulerian–Lagrangian simulation. Experiments were performed using a 6.5-m long, 0.075-m diameter horizontal pipe in conjunction with a laser Doppler anemometry (LDA) system. Spherical glass beads with two sizes, 1.5 and 2 mm, were used. Simulations were performed using the commercial discrete element method software EDEM, coupled with the computational fluid dynamics package FLUENT. Hybrid source terms were added to the conventional k–ε turbulence model to take into account the influence of the dispersed phase on the carrier phase turbulence intensity. The simulation results showed that the turbulence modulation depends strongly on the model parameter C_ε3. Both the standard k–ε turbulence model and the k–ε turbulence model with the hybrid source terms could predict the gas phase turbulence intensity trend only generally. A noticeable discrepancy in all cases between simulation and experimental results was observed, particularly for the regions close to the pipe wall. It was also observed that in some cases the addition of the source terms to the k–ε turbulence model did not improve the simulation results when compared with those of the standard k–ε turbulence model. Nonetheless, in the lower part of the pipe where particle loading was greater due to gravitational effects, the model with added source terms performed somewhat better.     

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.     

Modified correlation for minimum fluidization velocity of low-density particles in inverse liquid–solid fluidized beds

The minimum fluidization velocity (U_mf) is a key parameter for the scale-up of inverse liquid–solid fluidized beds. Theoretical predictions using common correlations were compared against experimental minimum fluidization velocity measurements of low density (28–638 kg/m³), 0.80–1.13 mm Styrofoam particles in a fluidized bed with a height of 4.5 m and 0.2 m diameter. The average absolute relative deviation for the predicted minimum fluidization velocity for particles below 300 kg/m³ was above 40% using the studied common correlations. A modified Wen and Yu correlation was thus proposed based on novel and past measurements with low-density and small-diameter particles, expanding the range for predicting U_mf. The new correlation predicted U_mf with deviations below 15% for ST028, ST122 and ST300. This modified correlation also improved U_mf predictions for comparable particles from a previous study, demonstrating its validity for a larger range of low-density particles.     

Effect of structure parameters on forces acting on baffles used in gas–solids fluidized beds

Effects of some important structural parameters, i.e. slat pitch, and layout position, on dynamic forces acting on the baffles were examined in the fluidized bed of FCC particles operating under different superficial gas velocities. The experimental baffles were made of multiple inclined slats. We found that the forces acting on the baffles decreased significantly with reducing pitch between the slats. For the baffles with a small slat pitch, the forces acting on the baffles increased slightly and then decreased with increasing superficial gas velocity, which is very different from the measured results of a single slat or tube immersed in fluidized beds. The different results are greatly related to the appearance of the “gas cushion” beneath the baffles, whose height increases with increasing superficial gas velocity. On the other hand, a region with stronger particle circulation induced by the inclined slat array was observed in the experiments. The slat near the wall and located below the region of downward-flowing particles was found to be subjected to the severest forces. Therefore, the slats located in similar locations of industrial baffles are suggested to be reinforced to increase their structural strength.     

A novel dual-material probe for in situ measurement of particle charge densities in gas–solid fluidized beds

Particle charge density is vitally important for monitoring electrostatic charges and understanding particle charging behavior in fluidized beds. In this paper, a dual-material probe was tested in a gas–solid fluidized bed for measuring the charge density of fluidized particles. The experiments were conducted in a two-dimensional fluidized bed with both single bubble injection and freely bubbling, at various particle charge densities and superficial gas velocities. Uniformly sized glass beads were used to eliminate complicating factors at this early stage of probe development. Peak currents, extracted from dynamic signals, were decoupled to determine charge densities of bed particles, which were found to be qualitatively and quantitatively consistent with charge densities directly measured by Faraday cup from the freely bubbling fluidized bed. The current signals were also decoupled to estimate bubble rise velocities, which were found to be in reasonable agreement with those obtained directly by analyzing video images.     

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.