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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Effect of Gaussian size distribution width on minimum fluidization velocity in tapered gas–solid fluidized beds

This paper investigated the effect of Gaussian distribution width, average particle diameter, particle loading, and the tapered angle on minimum fluidization velocity (U_mf) by conducting extensive experiments in tapered fluidized beds. Three powders with Gaussian size distribution and different distribution widths were used in the experiments. An increase in U_mf with increasing the average particle diameter, particle loading, and the tapered angle was observed. There was also a nonmonotonic behavior of U_mf as the Gaussian distribution width increased. An empirical correlation including dimensionless groups for predicting U_mf in tapered beds was developed in which the effect of distribution width was considered. The proposed correlation predictions were in good agreement with the experimental data, with a maximum deviation of 16.5% and average and standard deviations of, respectively, 6.4% and 7.4%. The proposed correlation was also compared with three earlier models, and their accuracy was discussed.     

Determination of hydrodynamic behavior of gas–solid fluidized beds using statistical analysis of acoustic emissions

In the processes involving the movement of solid particles, acoustic emissions are caused by particle friction, collision and fluid turbulence. Particle behavior can therefore be monitored and characterized by assessing the acoustic emission signals. Herein, extensive measurements were carried out by microphone at different superficial gas velocities with different particle sizes. Acoustic emission signals were processed using statistical analysis from which the minimum fluidization velocity was determined from the variation of standard deviation, skewness and kurtosis of acoustic emission signals against superficial gas velocity. Initial minimum fluidization velocity, corresponding to onset of fluidization of finer particles in the solids mixture, at which isolated bubbles occur, was also detected by this method. It was shown that the acoustic emission measurement is highly feasible as a practical method for monitoring the hydrodynamics of gas–solid fluidized beds.     

X-ray measurements on the influence of optical probes on gas–solid fluidized beds

1. Download : Download full-size image

     

Detailed parametric design methodology for hydrodynamics of liquid–solid circulating fluidized bed using design of experiments

A design-of-experiments methodology is used to develop a statistical model for the prediction of the hydrodynamics of a liquid–solid circulating fluidized bed. To illustrate the multilevel factorial design approach, a step by step methodology is taken to study the effects of the interactions among the independent factors considered on the performance variables. A multilevel full factorial design with three levels of the two factors and five levels of the third factor has been studied. Various statistical models such as the linear, two-factor interaction, quadratic, and cubic models are tested. The model has been developed to predict responses, viz., average solids holdup and solids circulation rate. The validity of the developed regression model is verified using the analysis of variance. Furthermore, the model developed was compared with an experimental dataset to assess its adequacy and reliability. This detailed statistical design methodology for non-linear systems considered here provides a very important tool for design and optimization in a cost-effective approach.