Horizontal gas and gas/solid jet penetration in a gas–solid fluidized bed

In this paper, the real time, dynamic phenomena of the three-dimensional horizontal gas and gas/solid mixture jetting in a 0.3 m (12 in) bubbling gas–solid fluidized bed are reported. The instantaneous properties of the shape of the jets and volumetric solids holdup are qualified and quantified using the three-dimensional electrical capacitance volume tomography (ECVT) recently developed in the authors’ group. It is found that the horizontal gas jet is almost symmetric along the horizontal axis during its penetration. As the jet width expands, the total volume of the gas jet increases. A mechanistic model is also developed to account for the experimental results obtained in this study. Comparison of jet penetration length and width between the model prediction and ECVT experiment shows that both the maximum penetration length and the maximum width of the horizontal gas jet increase with the superficial gas velocity. When the horizontal gas jet coalesces with a bubble rising from the bottom distributor, it loses its symmetric shape and can easily penetrate into the bed. For the horizontal gas/solid mixture jet penetration in the bed, the tail of the jet at the nozzle shrinks and the jet loses its jet shape immediately when the jet reaches its maximum penetration length, which are different from the characteristics exhibited by the gas jet. The solids holdup in the core region of the gas/solid mixture jet is higher than that in the gas jet. The penetration length of the horizontal gas/solid mixture jet is also larger than that of the gas jet.     

The onset velocity of a liquid–solid circulating fluidized bed

The critical transition velocity, U_cr, previously defined by Liang et al. [W.-G. Liang, S.-L. Zhang, J.-X. Zhu, Y. Jin, Z.-Q. Yu and Z.-W. Wang, Flow characteristics of the liquid–solid circulating fluidized bed, Powder Technol., 90 (1997) 95–102.] to demarcate the liquid–solid conventional and circulating fluidization regimes, was found to vary with the total solids inventory and the solids feeding system. In this work, an onset velocity for circulating fluidization regime, U_cf, is proposed to give the lowest U_cr value and to provide a convenient demarcation velocity that is independent of system geometry. This liquid velocity is obtained by measuring the time required to empty all particles in a batch operated fluidized bed under different liquid velocities. This method can be used for a wide range of particles and involves less influence of the operating conditions such as the solids inventory and the solids feeding system. Compared to the critical transition velocity, this newly defined onset velocity is a more intrinsic parameter, only dependent on the liquid and particle properties. Based on the experimental results obtained in this work and other published results, the influence of particle properties and equipment setup on the onset velocity is also discussed.     

Investigation on gas–solid flow regimes in a novel multistage fluidized bed

Gas–solid flow regime in a novel multistage circulating fluidized bed is investigated in this study. Pressure fluctuations are first sampled from gas–solid flow systems and then are analyzed through frequency and time–frequency domain methods including power spectrum and Hilbert–Huang transform. According to the flow characteristics obtained from pressure fluctuations, it is found that the gas–solid motions in the multistage circulating fluidized bed exhibit two dominant motion peaks in low and high frequencies. Moreover, gas-cluster motions become intensive for the multistage circulating fluidized bed in comparison with the fast bed. Unlike the traditional methods, the fuzzy C-means clustering method is introduced to objectively identify flow regime in the multistage circulating fluidized bed on the basis of the flow characteristics extracted from bubbling, turbulent, fast, and multistage fluidized beds. The identification accuracy of fuzzy C-means clustering method is first verified. The identification results show that the flow regime in the multistage circulating fluidized bed is in the scope of fast flow regime under examined conditions. Moreover, the results indicate that the consistency of flow regime between two enlarged sections exists. In addition, the transition onset of fast flow regime in the multistage circulating fluidized bed is higher than that in the fast bed.     

Discrete particle simulation of solids motion in a gas–solid fluidized bed

The solids motion in a gas–solid fluidized bed was investigated via discrete particle simulation. The motion of individual particles in a uniform particle system and a binary particle system was monitored by the solution of the Newton's second law of motion. The force acting on each particle consists of the contact force between particles and the force exerted by the surrounding fluid. The contact force is modeled by using the analogy of spring, dash-pot and friction slider. The flow field of gas was predicted by the Navier–Stokes equation. The solids distribution is non-uniform in the bed, which is very diluted near the center but high near the wall. It was also found that there is a single solids circulation cell in the fluidized bed with ascending at the center and descending near the wall. This finding agrees with the experimental results obtained by Moslemian. The effects of the operating conditions, such as superficial gas velocity, particle size, and column size on the solids movement, were investigated. In the fluidized bed containing uniform particles better solids mixing was found in the larger bed containing smaller size particles and operated at higher superficial gas velocity. In the system containing binary particles, it was shown that under suitable conditions the particles in a fluidized bed could be made mixable or non-mixable depending on the ratios of particle sizes and densities. Better mixing of binary particles was found in the system containing particles with less different densities and closer sizes. These results were found to follow the mixing and segregation criteria obtained experimentally by Tanaka et al.     

Euler multiphase-CFD simulation on a bubble-driven gas–liquid–solid fluidized bed

An integrated flow model was developed to simulate the fluidization hydrodynamics in a new bubble-driven gas–liquid–solid fluidized bed using the computational fluid dynamic (CFD) method. The results showed that axial solids holdup is affected by grid size, bubble diameter, and the interphase drag models used in the simulation. Good agreements with experimental data could be obtained by adopting the following parameters: 5 mm grid, 1.2 mm bubble diameter, the Tomiyama gas–liquid model, the Schiller–Naumann liquid–solid model, and the Gidaspow gas–solid model. At full fluidization state, an internal circulation of particles flowing upward near the wall and downward in the centre is observed, which is in the opposite direction compared with the traditional core-annular flow structure in a gas–solid fluidized bed. The simulated results are very sensitive to bubble diameters. Using smaller bubble diameters would lead to excessive liquid bed expansions and more solid accumulated at the bottom due to a bigger gas–liquid drag force, while bigger bubble diameters would result in a higher solid bed height caused by a smaller gas–solid drag force. Considering the actual bubble distribution, population balance model (PBM) is employed to characterize the coalescence and break up of bubbles. The calculated bubble diameters grow up from 2–4 mm at the bottom to 5–10 mm at the upper section of the bed, which are comparable to those observed in experiments. The simulation results could provide valuable information for the design and optimization of this new type of fluidized system.     

Non-intrusive determination of bubble size in a gas–solid fluidized bed: An evaluation

Bubble sizes measured in a column of diameter 290 mm with FCC particles utilizing both an intrusive optical probe and non-intrusive pressure analysis are compared. The pressure signals were decoupled by differential pressure analysis and incoherence analysis. It is shown that pressure fluctuations induced by jetting/bubble formation can be effectively filtered out by differential pressure and incoherence analysis. The differential pressure signals measured across a vertical interval less than half the maximum bubble size unreasonably damps the power spectral density intensity, leading to underestimation of bubble size and overestimation of mean frequency. In the present work, the incoherence analysis tends to estimate greater bubble size than differential pressure analysis. Bubble chord lengths are overestimated by optical probe signals because small bubbles are not detected. Bubble sizes calculated by the equation of Horio and Nonaka (1987) agree reasonably well with that estimated by incoherence analysis at relatively high superficial gas velocities.     

Local solid mixing in gas–solid fluidized beds

Diffusivity of the solid particles in a 152-mm ID gas–solid fluidized bed was determined at different regimes of fluidization. The gas was air at room temperature and atmospheric pressure and the solids were 385 μm sand or 70 μm FCC particles. The experiments were done at superficial gas velocities from 0.5 to 2.8 m/s for sand and 0.44 to 0.9 m/s for FCC (in both bubbling and turbulent regimes). Movement of a tracer was monitored by radioactive particle tracking (RPT) technique. Once the time-position data became available, local axial and radial diffusivity of solids were calculated from these data. Calculated diffusivities are in the range of 3.3×10⁻³ to 5.6×10⁻² m²/s for axial and 2.6×10⁻⁴ to 1.5×10⁻³ m²/s for radial direction. The results show that the diffusivities, both axial and radial, increase with superficial gas velocity and are linearly correlated to the axial solid velocity gradient. Solid diffusivity in a bed of FCC was found to be higher than that of a bed of sand at the same excess superficial gas velocity.     

Development and verification of coarse-grain CFD-DEM for nonspherical particles in a gas–solid fluidized bed

Computational fluid dynamics coupled with discrete element method (CFD-DEM) has been widely used to understand the complicated fundamentals inside gas–solid fluidized beds. To realize large-scale simulations, CFD-DEM integrated with coarse-grain model (CG CFD-DEM) provides a feasible solution, and has led to a recent upsurge of interest. However, when dealing with large-scale simulations involving irregular-shaped particles such as biomass particles featuring elongated shapes, current CG models cannot function as normal because they are all developed for spherical particles. To address this issue, a CG CFD-DEM for nonspherical particles is proposed in this study, and the morphology of particles is characterized by the super-ellipsoid model. The effectiveness and accuracy of CG CFD-DEM for nonspherical particles are comprehensively evaluated by comparing the hydrodynamic behaviors with the results predicted by traditional CFD-DEM in a gas–solid fluidized bed. It is demonstrated that the proposed model can accurately model gas–solid flow containing nonspherical particles, merely the particle dynamics are somewhat lost due to the scaleup of particle size. Finally, the calculation efficiency of CG CFD-DEM is assessed, and the results show that CG CFD-DEM can largely reduce computational costs mainly by improving the calculation efficiency of DEM. In general, the proposed CG CFD-DEM for nonspherical particles strikes a good balance between efficiency and accuracy, and has shown its prospect as a high-efficiency alternative to traditional CFD-DEM for engineering applications involving nonspherical particles.     

An assessment of the ability of computational fluid dynamic models to predict reactive gas–solid flows in a fluidized bed

Fine grid, two dimensional simulations of reactive gas–solid flows occurring in a fluidized bed reactor were carried out using the Eulerian multi-fluid kinetic theory of granular flow (KTGF) approach in the commercial flow solver, ANSYS FLUENT 12.1. The fuel reactor of a pilot scale Chemical Looping Combustion rig, operated in the bubbling fluidization regime at the Vienna University of Technology, was simulated. Grid dependence studies were carried out as well as sensitivity studies to the fuel inlet condition and the inclusion of gas phase turbulence. Simulations could not accurately reproduce the experimental trend for the case when highly reactive nickel oxide was used as the oxygen carrier material, but in general satisfactory quantitative agreement was observed. The failure to correctly capture the experimental trend was primarily attributed to the fine length-scales at the feed gas inlets not being adequately resolved even at the finest grid investigated. The trend quickly worsened when coarser grids were used, indicating that the generality of the model is lost when grid dependence effects are present. A number of possible dimensional effects were also discussed. Subsequently, the model was used to successfully capture another experimental trend obtained with a much less reactive ilmenite oxygen carrier material. The model captured this trend correctly because the reaction was now limited by the reaction rate and not by species transfer to the large scale gas-emulsion interfaces. Results were therefore not as sensitive to the correct hydrodynamic modelling of the interface, especially near the gas inlets, and the model retained its generality over a wide range of operating conditions.     

Numerical simulation of local and global mixing/segregation characteristics in a gas–solid fluidized bed

Researches on solids mixing and segregation are of great significance for the operation and design of fluidized bed reactors. In this paper, the local and global mixing and segregation characteristics of binary mixtures were investigated in a gas–solid fluidized bed by computational fluid dynamics-discrete element method (CFD-DEM) coupled approach. A methodology based on solids mixing entropy was developed to quantitatively calculate the mixing degree and time of the bed. The mixing curves of global mixing entropy were acquired, and the distribution maps of local mixing entropy and mixing time were also obtained. By comparing different operating conditions, the effects of superficial gas velocity, particle density ratio and size ratio on mixing/segregation behavior were discussed. Results showed that for the partial mixing state, the fluidized bed can be divided into three parts along the bed height: complete segregation area, transition area and stable mixing area. These areas showed different mixing/segregation processes. Increasing gas velocity promoted the local and global mixing of binary mixtures. The increase in particle density ratio and size ratio enlarged the complete segregation area, reduced the mixing degree and increased the mixing time in the stable mixing area.     

Flow regimes in a gas–liquid–solid three-phase moving bed

The gas–liquid–solid three-phase moving beds could supply a potential solution for multiphase reactions with catalyst easily deactivated, and the flow regimes in it were studied by optical method and pressure drop measurement. Results showed that taking the trickle flow as the initial flow regime, the flow channels were more obvious as the particle velocity increased. When the initial flow regimes were pulse flow and bubble flow respectively, the pulse-to-trickle and bubble-to-pulse flow transitions mainly occurred at moderate-to-high particle velocities (0.01–0.04 m s⁻¹ under conditions used in this work). Moreover, the flow regime map in the three-phase moving bed was constructed and shown that the region of trickle flow increased and the region of bubble flow decreased. Finally, the application of three-phase moving beds was discussed, and it could be suitable for those reactions, which had to operate in the pulse flow, bubble flow, and transition zone.     

Hydrodynamics of gas–solid flow in the circulating fluidized bed reactor for dry flue gas desulfurization

Process design and scale-up require a fundamental understanding of the hydrodynamics of gas–solid flow in the circulating fluidized bed flue gas desulfurization (CFB-FGD) reactor although the CFB system has been widely used in flue gas desulfurization and flue gas cleaning processes. The hydrodynamics in the CFB-FGD reactor model was investigated by pressure measurements and specially designed sampling probe based on three dimensionless groups for practicable similarity of industrial CFB-FGD process. The results show that the pressure drop in the venturi section is predominant as high as 60% of the total pressure drop and the total pressure drop significantly increases with the increasing external solid circulating rates at the same superficial gas velocity. Moreover, the measurements of radial solid mass fluxes show that the flow pattern in the CFB-FGD reactor is a typical core–annulus flow and this flow structure prevails until the top of the reactor. Reflux ratios are used to quantitatively evaluate the internal solid reflux in the reactor and the values in the low section of the reactor are much higher than those in the upper section.     

Effect of the injection-nozzle geometry on the interaction between a gas–liquid jet and a gas–solid fluidized bed

In industrial fluid cokers, bitumen is first mixed with steam in a premixer, and then fed to the atomization nozzle. The objective of this work was to evaluate the impact of both the premixer and the nozzle geometrical configuration on the quality of the liquid–solid contact resulting from injections of liquid into a gas–solid fluidized bed. To assess the quality of the liquid–solid contact a method based on electric conductance measurements of the bed material previously developed by the authors [9] was used. Liquid atomization efficiency in open air, spray geometry, and spray stability were also characterized to evaluate their effects on the nozzle spraying performance within the fluidized bed. This study indicated that spray stability is highly beneficial to the liquid–solid contact efficiency. In particular, fluid constrictions such as the series of converging and diverging sections within the nozzle have a stabilizing effect on the spray. Future optimization of the existing liquid-injection systems should consider alternative gas–liquid premixers and nozzle geometries to enhance the jet stability.     

Dynamic characteristics of the volatile cloudy zone in a gas–solid fluidized bed with hydrocarbon liquid spray

In a gas–solid fluidized bed with continuous hydrocarbon liquid spray, a volatile “cloudy zone” could be formed, defined as a dynamically steady liquid-affected zone, including droplets, wet particles, and the gas which passes through the zone. A new flow pattern with the dynamic coexistence of cloudy zone and non-cloudy zone (gas–solid zone), is accordingly established. The temperature, particle concentration, and particle velocity fields are measured in real-time via infrared thermography and particle imaging velocimetry, respectively. Results show that the area and range of central position of the cloudy zone illustrate a heavier fluctuant trend with the increasing velocity of liquid spray, and the main frequency of area fluctuation is close to that of the bubble rising. Moreover, the particle concentration and particle velocity in the cloudy zone are lower than those in the non-cloudy zone, breaking the conventional symmetrical distributions of hydrodynamic parameters of particles in a gas–solid fluidized bed.     

Supersonic attrition nozzles in gas–solid fluidized beds

Gas–solid fluidized beds are used in both catalytic and non-catalytic processes, and some of the industrial applications are fluid catalytic cracking, polyethylene production, drying and classification, coating, and granulation. In some applications, the size distribution of the bed particles must be controlled in order to maintain good fluidization, and attrition nozzles can be used for this purpose. Supersonic attrition nozzles are more efficient than subsonic nozzles, and, in this study, different geometries of the Laval nozzle, a convergent–divergent (C–D) nozzle, have been investigated. The geometry of this type of nozzles gives supersonic velocities under the right operating conditions.     

Computational fluid dynamics for dense gas–solid fluidized beds: a multi-scale modeling strategy

Dense gas–particle flows are encountered in a variety of industrially important processes for large scale production of fuels, fertilizers and base chemicals. The scale-up of these processes is often problematic, which can be related to the intrinsic complexities of these flows which are unfortunately not yet fully understood despite significant efforts made in both academic and industrial research laboratories. In dense gas–particle flows both (effective) fluid–particle and (dissipative) particle–particle interactions need to be accounted for because these phenomena, to a large extent, govern the prevailing flow phenomena, i.e. the formation and evolution of heterogeneous structures. These structures have significant impact on the quality of the gas–solid contact and as a direct consequence thereof strongly affect the performance of the process.Due to the inherent complexity of dense gas-particles flows, we have adopted a multi-scale modeling approach in which both fluid–particle and particle–particle interactions can be properly accounted for. The idea is essentially that fundamental models, taking into account the relevant details of fluid–particle (lattice Boltzmann model (LBM)) and particle–particle (discrete particle model (DPM)) interactions, are used to develop closure laws to feed continuum models which can be used to compute the flow structures on a much larger (industrial) scale. Our multi-scale approach (see Fig. 1) involves the LBM, the DPM, the continuum model based on the kinetic theory of granular flow, and the discrete bubble model. In this paper we give an overview of the multi-scale modeling strategy, accompanied by illustrative computational results for bubble formation. In addition, areas which need substantial further attention will be highlighted.     

A fundamental CFD study of the gas–solid flow field in fluidized bed polymerization reactors

A Eulerian–Eulerian model incorporating the kinetic theory of granular flow was applied to describe the gas–solid two-phase flow in fluidized bed polymerization reactors. The model parameters were examined, and the model was validated by comparing the simulation result with the classical calculated data. The effects of distributor shape, solid particle size, operational gas velocity and feed manner on the flow behavior in the reactor were also investigated numerically. The results show that with the increase of solid particle diameter, the bubble numbers decrease and the bubble size increases, resulting in a smaller bed expansion ratio. Bed expansion ratio increases with increasing the gas inlet velocity. Moreover, the final fluidized qualities are almost the same for the plane distributor case and the triangle distributor case. There exists a tempestuous wiggle from side to side in the bed at the continuous feed manner, which could not be obtained at a batch feed manner.     

Experimental research on steady-state operation characteristics of gas–solid flow in a 15.5 m dual circulating fluidized bed system

Gas-solid hydrodynamic steady-state operation is the operating basis in a chemical looping dual-reactor system. This study reported the experimental results on the steady-state operation characteristics of gas–solid flow in a 15.5 m high dual circulating fluidized bed(CFB) cold test system. The effects of superficial gas velocity, static bed material height and solid returning modes on the steady-state operation characteristics between the two CFBs were investigated. Results suggest that the solid distributions in the dual CFB test system was mainly determined by the superficial gas velocity and larger solid inventory may help to improve the solid distributions. Besides, cross-returning mode coupled with self-returning is good for steady-state running in the dual-reactor test system.     

Total and gas convective heat transfer from a vertical tube to a mixed particle gas—solid fluidized bed

Experimental data were obtained for the average gas convective and total heat transfer coefficients for a vertical tube immersed in an air-fluidized bed of narrowly as well as widely distributed particle size mixtures. The gas convective heat transfer coefficient was determined by measuring the rate of mass loss from a vertical naphthalene tube 0.0262 m in diameter and 0.1012 m in length and using a heat and mass transfer analogy. These data were obtained at a bed temperature of about 330 K and superficial velocity of 0.1 to 1.1 m/s. The total heat transfer coefficients were measured under identical conditions using an electrically heated vertical tube. The total heat transfer coefficient decreased with an increase in particle diameter from 0.237 to 1.35 mm. The addition of fines was found to increase the total heat transfer coefficient. The gas convective heat transfer coefficient increased with increase in particle size and fluidizing velocity. The dependence of the gas convective heat transfer coefficient on gas velocity was more pronounced for large particles. The addition of fines resulted in decrease in gas convective coefficient. The relative contribution of the gas convective component of heat transfer coefficient was found to increase with increase in particle diameter. Its dependency on fluidizing velocity was found to be more complex. The experimental data were compared with the existing heat transfer models and correlations.