CFD simulation of a gas–solid fluidized bed with two vertical jets

A computational fluid dynamics (CFD) model is used to investigate the hydrodynamics of a gas–solid fluidized bed with two vertical jets. Sand particles with a density of 2660 kg/m³ and a diameter of 5.0 × 10^?4 m are employed as the solid phase. Numerical computation is carried out in a 0.57 m × 1.00 m two-dimensional bed using a commercial CFD code, CFX 4.4, together with user-defined Fortran subroutines. The applicability of the CFD model is validated by predicting the bed pressure drop in a bubbling fluidized bed, and the jet detachment time and equivalent bubble diameter in a fluidized bed with a single jet. Subsequently, the model is used to explore the hydrodynamics of two vertical jets in a fluidized bed. The computational results reveal three flow patterns, isolated, merged and transitional jets, depending on the nozzle separation distance and jet gas velocity and influencing significantly the solid circulation pattern. The jet penetration depth is found to increase with increasing jet gas velocity, and can be predicted reasonably well by the correlations of Hong et al. (2003) for isolated jets and of Yang and Keairns (1979) for interacting jets.     

Orientation of cylindrical particles in gas–solid circulating fluidized bed

The orientation of cylindrical particles in a gas–solid circulating fluidized bed was investigated by establishing a three-dimensional Euler–Lagrange model on the basis of rigid kinetics, impact kinetics and gas–solid two-phase flow theory. The resulting simulation indicated that the model could well illustrate the orientation of cylindrical particles in a riser during fluidization. The influences of bed structure and operation parameters on orientation of cylindrical particles were then studied and compared with related experimental results. The simulation results showed that the majority of cylindrical particles move with small nutation angles in the riser, the orientation of cylindrical particles is affected more obviously by their positions than by their slenderness and local gas velocities. The simulation results well agree with experiments, thus validating the proposed model and computation.     

Numerical simulation of gas–solid flow with two fluid model in a spouted-fluid bed

The flow characteristics in a spouted-fluid bed differ from those in spouted or fluidized beds because of the injection of the spouting gas and the introduction of a fluidizing gas. The flow behavior of gas–solid phases was predicted using the Eulerian–Eulerian two-fluid model (TFM) approach with kinetic theory for granular flow to obtain the flow patterns in spouted-fluid beds. The gas flux and gas incident angle have a significant influence on the porosity and particle concentration in gas–solid spouted-fluid beds. The fluidizing gas flux affects the flow behavior of particles in the fountain. In the spouted-fluid bed, the solids volume fraction is low in the spout and high in the annulus. However, the solids volume fraction is reduced near the wall.     

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.     

Simulations of vertical jet penetration using a filtered two-fluid model in a gas–solid fluidized bed

The influence of a vertical jet located at the distributor in a cylindrical fluidized bed on the flow behavior of gas and particles was predicted using a filtered two-fluid model proposed by Sundaresan and coworkers. The distributions of volume fraction and the velocity of particles along the lateral direction were investigated for different jet velocities by analyzing the simulated results. The vertical jet penetration lengths at the different gas jet velocities have been obtained and compared with predictions derived from empirical correlations; the predicted air jet penetration length is discussed. Agreement between the numerical simulations and experimental results has been achieved.     

The characteristics of gas–solid flow and wall heat transfer in a fluidized bed reactor

Numerical study using computational fluid dynamics has been carried out to investigate the heat transfer characteristics of a laboratory fluidized bed reactor. The fluidized bed reactor of vTI (Johann Heinrich von Thünen-Institute)-Institute of Wood Technology and Wood Biology is modeled. For the simulation of multiphase flow and thermal fields, an Eulerian–Eulerian approach is applied. The flow and thermal characteristics of the reactor are fully investigated for the wide range of superficial gas velocities and two different particle diameters. In particular, the contributions of the gas bubble and emulsion phase flows on the wall heat transfer are scrutinized. From the predicted results, it is fully elucidated that particular near-wall bubble motions mainly govern the wall heat transfer.     

Apparent viscosity of mixed particle system in pulsed gas–solid separation fluidized bed

The apparent viscosity reflects the resistance of the fluidized medium in the bed to the beneficiation particles, which directly affects the separation time and mismatch content. So, the falling-ball method was used to measure the apparent viscosity of a binary medium in a pulsed fluidized bed by varying the gas velocity, pulsation frequency, and fine particle content. The results show that with increasing gas velocity and fine particle content, the apparent viscosity of the bed gradually decreased, whereas it first decreased and then increased with pulsation frequency increasing and achieved a minimum value in the range of 4–6 Hz. Within limits, the adjustment of gas velocity and fine content can effectively reduce the apparent viscosity and improve the separation process. A model for predicting the apparent viscosity in a pulsation separation fluidized bed was established with good accuracy.     

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.     

Numerical simulations of gas–liquid–solid flows in a hydrocyclone separator

The flow behavior in hydrocyclones is quite complex. In this study, the computational fluid dynamics (CFD) method was used to simulate the flow fields inside a hydrocyclone in order to investigate its separation efficiency. In the computational fluid dynamics study of hydrocyclones, the air-core dimension is a key to predicting the mass split between the underflow and overflow. In turn, the mass split influences the prediction of the size classification curve. Three models, the model, the Reynolds stress model (RSM) without considering the air-core, and the Reynolds stress turbulence model with the volume of fluid (VOF) multiphase model for simulating the air-core, were compared in terms of their predictions of velocity, axial and tangential velocity distributions, and separation proportion. The RSM with air-core simulation model, since it reproduces some detailed features of the turbulence and multiphase, clearly predicted the experimental data more closely than did the other two models.     

Modeling of multiple noncatalytic gas–solid reactions in a moving bed of porous pellets based on finite volume method

A mathematical model is presented to simulate the multiple heterogeneous reactions with complex set of physicochemical and thermal phenomena in a moving bed of porous pellets. This model is based on both heat and mass transfer phenomena of gaseous species in a porous medium including chemical reactions at interfaces whose areas vary during the conversion. This model accounts for both the exothermic and endothermic reactions which can be equimolar or nonequimolar. Furthermore it considers simultaneously the reactions in the nonisothermal transient condition. A powerful technique based upon finite volume fully implicit approach has been implemented to solve the complicated governing equations numerically. The model has been validated by comparing with various experimental and analytical results in two cases: the single pellet scale as well as the counter current moving bed reactor.

Transportation characteristics of gas–solid two-phase flow in a long-distance pipeline

In this study, experiments on fly ash conveying were carried out with a home-made long-distance positive-pressure pneumatic conveying system equipped with a high performance electrical capacitance tomography system to observe the transient characteristics of gas–solid two-phase flow. The experimental results indicated that solids throughput increased with increasing solids–gas ratio when the conveying pipeline was not plugged. Moreover, the optimum operating state was determined for the 1000 m long conveying pipeline with a throttle plate of 26 orifices. At this state the solids throughput was about 12.97 t/h. Additionally, the transportation pattern of fly ash gradually changed from sparse–dense flow to partial and plug flows with increasing conveying distance because of the conveying pressure loss. These experimental results provide important reference data for the development of pneumatic conveying technology.     

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.     

Hydrodynamics in a new liquid–solid circulating conventional fluidized bed

A new type of liquid–solid fluidized bed, named circulating conventional fluidized bed (CCFB) which operates below particle terminal velocity was proposed and experimentally studied. The hydrodynamic behavior was systematically studied in a liquid–solid CCFB of 0.032 m I.D. and 4.5 m in height with five different types of particles. Liquid–solid fluidization with external particle circulation was experimentally realized below the particle terminal velocity. The axial distribution of local solids holdup was obtained and found to be fairly uniform in a wide range of liquid velocities and solids circulation rates. The average solids holdup is found to be significantly increased compared with conventional fluidization at similar conditions. The effect of particle properties and operating conditions on bed behavior was investigated as well. Results show that particles with higher terminal velocity have higher average solids holdup.     

Numerical investigations on gas–solid flow in circulating fluidized bed risers using a new cluster-based drag model

A cluster-based drag model is proposed for the gas–solid circulating fluidized bed (CFB) riser by including the cluster information collected from image processing and wavelet analysis into the calculation of system drag. The performance of the proposed drag model is compared with some commonly used drag models. A good agreement with the experimental data is achieved by the proposed cluster-based drag model. Error analysis of the proposed cluster-based drag model based on the local distributions of solids holdup and particle velocity is conducted. The clustering phenomenon in the low-density and high-density CFB risers and the effect of the cluster size on the simulation accuracy are also numerically studied by the proposed drag model.     

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.     

Effect of coal particle swarm properties on the fluidization characteristics and coal beneficiation in a dense-phase gas–solid fluidized bed

This paper analyzes the influence of different coal mass fraction in an air dense medium fluidized bed (ADMFB). The effect of the low density particles layer on heavy sedimentation increased with increasing material layer thickness. The thickness of the low density particles layer also affected the final settling time of the high density particles. Increasing the thickness of the low density particles layer by Δh provoked an increase in the settling of high density particles that was related to their diameter (Δh/D). The pressure gradient across the bed was lower than that observed for the control experiment, which had only the dense material, owing to a decrease in the pressure gradient in Zones 1 and 5 (at the top and bottom of the bed, respectively). Introducing different coal sizes resulted in different fluidization environments, particle accumulation layers, and changes to the surrounding zone. However, the influence of the coal particles on the local bed characteristics was related to its concentration. The feeding mass fraction of 6–13 mm size and 13–25 mm size coal should be limited to10% and 13%, respectively. The ranges of possible deviation were found to be 0.08–0.15 and 0.07–0.10 for the respective samples.     

Particle size distribution in CPFD modeling of gas–solid flows in a CFB riser

A computational particle fluid dynamics (CPFD) numerical method to model gas–solid flows in a circulating fluidized bed (CFB) riser was used to assess the effects of particle size distribution (PSD) on solids distribution and flow. We investigated a binary PSD and a polydisperse PSD case. Our simulations were compared with measured solids concentrations and velocity profiles from experiments, as well as with a published Eulerian-Eulerian simulation. Overall flow patterns were similar for both simulation cases, as confirmed by experimental measurements. However, our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser, as seen in other numerical studies. Above this bottom region, distributions of particle volume fraction and particle vertical velocity were consistent with our experiments, and the simulated average particle diameter decreased as a power function with riser height. Interactions between particles and walls also were successfully modeled, with accurate predictions for the lateral profiles of particle vertical velocity. It was easy to implement PSD into the CPFD numerical model, and it required fewer computational resources compared with other models, especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Dynamics of annular gas–liquid two-phase swirling jets

The dynamics of annular gas–liquid two-phase swirling jets have been examined by means of direct numerical simulation and proper orthogonal decomposition. An Eulerian approach with mixed-fluid treatment, combined with an adapted volume of fluid and a continuum surface force model, was used to describe the two-phase flow system. The unsteady, compressible, three-dimensional Navier–Stokes equations have been solved by using highly accurate numerical methods. Two computational cases have been performed to examine the effects of liquid-to-gas density ratio on the flow development. It was found that the higher density ratio case is more vortical with larger spatial distribution of the liquid, in agreement with linear theories. Proper orthogonal decomposition analysis revealed that more modes are of importance at the higher density ratio, indicating a more unstable flow field. In the lower density ratio case, both a central and a geometrical recirculation zone are captured while only one central recirculation zone is evident at the higher density ratio. The results also indicate the formation of a precessing vortex core at the high density ratio, indicating that the precessing vortex core development is dependent on the liquid-to-gas density ratio of the two-phase flow, apart from the swirl number alone.     

Numerical study of gas–solid two-phase flow and erosion in a cavity with a slope

Gate valve is mainly used to turn on or turn off the pipeline in pneumatic conveying. When the gate valve is fully open, the particles are easy to collide with the cavity rear wall and enter into the cavity, resulting in particles’ accumulation in the cavity. The particles in cavity will accumulate between the cavity bottom and the flashboard bottom wall and prevent the gate from turning off normally. Meanwhile, the particles’ collision with cavity rear wall will cause serious erosion. Both the particles’ accumulation and erosion will cause the poor sealing of the gate valve, further resulting in the leakage of the pipeline system. To reduce the particles’ accumulation in cavity and erosion on cavity when the gate valve is fully open, we simplify the gate valve into a cavity structure and study it. We find that adding a slope upstream the cavity can effectively reduce the particles’ accumulation in the cavity and the erosion on the cavity rear wall. In this work, Eulerian–Lagrangian method in commercial code (FLUENT) was used to study the gas–solid two-phase flow and erosion characteristics of a cavity with a slope. The particle distribution shows that the particles with Stokes number St = 1.3 and St = 13 cannot enter the cavity due to the slope, but the particles with St = 0.13 enter the cavity following the gas. For St = 13, the particles collide with the wall many times in the ideal cavity. Erosion results show that the slope can transfer the erosion on cavity rear wall to the slope and reduce the maximum erosion rate of the wall near the cavity to some degrees.