Experimental and computational studies on flow behavior of gas–solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally. Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695-0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

Experimental and computational studies on flow behavior of gas-solid fluidized bed with disparately sized binary particles

This paper presents experimental and computational studies on the flow behavior of a gas-solid fluidized bed with disparately sized binary particle mixtures. The mixing/segregation behavior and segregation efficiency of the small and large particles are investigated experimentally.Particle composition and operating conditions that influence the fluidization behavior of mixing/segregation are examined. Based on the granular kinetics theory, a multi-fluid CFD model has been developed and verified against the experimental results. The simulation results are in reasonable agreement with experimental data. The results showed that the smaller particles are found near the bed surface while the larger particles tend to settle down to the bed bottom in turbulent fluidized bed. However, complete segregation of the binary particles does not occur in the gas velocity range of 0.695--0.904 m/s. Segregation efficiency increases with increasing gas velocity and mean residence time of the binary particles, but decreases with increasing the small particle concentration. The calculated results also show that the small particles move downward in the wall region and upward in the core. Due to the effect of large particles on the movement of small particles, the small particles present a more turbulent velocity profile in the dense phase than that in the dilute phase.     

Mixing and segregation of solid particles in a conical spouted bed: Effect of particle size and density

In this work, the mixing and segregation of binary mixtures of particles with different sizes and densities in a pseudo-2D spouted bed were studied experimentally. A binary mixture of solid particles including sand, gypsum, and polyurethane was used. To determine the particles mass fraction, and their mixing and segregation in the bed, an image-processing technique was developed and used. Important hydrodynamic parameters, such as the axial and radial segregation profiles of the solid particles, were measured. The effects of air velocity, particle size, and particle mass fraction were also evaluated. The flow regime in the spouted bed and the time required for reaching the equilibrium state of the solid particles were discussed. The results showed that the segregation of solid particles and the time to equilibrium both decreased when the air velocity increased to much larger than the minimum spouting velocity. The axial segregation increased with the diameter ratio of the particles. Upon completion of the test, coarse particles were concentrated mainly in the spout region, while fine particles were aggregated in the annulus region. Examination of the flow pattern in the spouted bed showed that the particles near the wall had longer flow paths, while those near the spout region had shorter flow paths.     

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.     

Discrete particle modeling of granular temperature distribution in a bubbling fluidized bed

The discrete hard sphere particle model (DPM) is applied in this work to study numerically the distributions of particle and bubble granular temperatures in a bubbling fluidized bed. The dimensions of the bed and other parameters are set to correspond to those of Müller et al. (2008). Various drag models and operational parameters are investigated to find their influence on particle and bubble granular temperatures. Various inlet superficial gas velocities are used in this work to obtain their effect on flow characteristics. It is found that the superficial gas velocity has the most important effect on granular temperatures including bubble granular temperature, particle translational granular temperature and particle rotational granular temperature. The drag force model affects more seriously the large scale variables such as the bubble granular temperature. Restitution coefficient influences all granular temperatures to some degree. Simulation results are compared with experimental results by Müller et al. (2008) showing reasonable agreement.     

Lagrangian–Eulerian simulation of slugging fluidized bed

This work studies gas–solid slugging fluidized beds with Type-D particles, using two-dimensional simulations based on discrete element model (DEM). DEM performance is quantitatively validated by two commonly accepted correlations for determining slugging behavior. The voidage profiles simulated with bed height corresponding to Baeyens and Geldart (1974) correlation for onset of slugging demonstrate a transitional flow pattern from free bubbling to slugging. The present calculated values for the maximum slugging bed height are in good agreement with the correlation from Matsen et al. (1969). Simulations show that fluidized beds with Type-D particles can operate in the round-nosed slugging regime and also shows that wall slugs and square-nosed slugs tend to be formed with increase in superficial gas velocity and in bed height, respectively.     

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.     

Eulerian–Lagrangian simulation of distinct clustering phenomena and RTDs in riser and downer

Numerical simulation of fully developed hydrodynamics of a riser and a downer was carried out using an Eulerian–Lagrangian model, where the particles are modeled by the discrete element method (DEM) and the gas by the Navier–Stokes equations. Periodic flow domain with two side walls was adopted to simulate the fully developed dynamics in a 2D channel of 10 cm in width. All the simulations were carried out under the same superficial gas velocity and solids holdup in the domain, starting with a homogenous state for both gas and solids, and followed by the evolution of the dynamics to the heterogeneous state with distinct clustering in the riser and the downer. In the riser, particle clusters move slowly, tending to suspend along the wall or to flow downwards, which causes wide residence time distribution of the particles. In the downer, clusters still exist, but they have faster velocities than the discrete particles. Loosely collected particles in the clusters move in the same direction as the bulk flow, resulting in plug flow in the downer. The residence time distribution (RTD) of solids was computed by tracking the displacements of all particles in the flow direction. The results show a rather wide RTD for the solids in the riser but a sharp peak RTD in the downer, much in agreement with the experimental findings in the literature. The ensemble average of transient dynamics also shows reasonable profiles of solids volume fraction and solids velocity, and their dependence on particle density.     

Effect of viscoelasticity on the rotation of a sphere in shear flow

When particles are dispersed in viscoelastic rather than Newtonian media, the hydrodynamics will be changed entailing differences in suspension rheology. The disturbance velocity profiles and stress distributions around the particle will depend on the viscoelastic material functions. Even in inertialess flows, changes in particle rotation and migration will occur. The problem of the rotation of a single spherical particle in simple shear flow in viscoelastic fluids was recently studied to understand the effects of changes in the rheological properties with both numerical simulations [D’Avino et al., J. Rheol. 52 (2008) 1331–1346] and experiments [Snijkers et al., J. Rheol. 53 (2009) 459–480]. In the simulations, different constitutive models were used to demonstrate the effects of different rheological behavior. In the experiments, fluids with different constitutive properties were chosen. In both studies a slowing down of the rotation speed of the particles was found, when compared to the Newtonian case, as elasticity increases. Surprisingly, the extent of the slowing down of the rotation rate did not depend strongly on the details of the fluid rheology, but primarily on the Weissenberg number defined as the ratio between the first normal stress difference and the shear stress.In the present work, a quantitative comparison between the experimental measurements and novel simulation results is made by considering more realistic constitutive equations as compared to the model fluids used in previous numerical simulations [D’Avino et al., J. Rheol. 52 (2008) 1331–1346]. A multimode Giesekus model with Newtonian solvent as constitutive equation is fitted to the experimentally obtained linear and nonlinear fluid properties and used to simulate the rotation of a torque-free sphere in a range of Weissenberg numbers similar to those in the experiments. A good agreement between the experimental and numerical results is obtained. The local torque and pressure distributions on the particle surface calculated by simulations are shown.     

Segregation characteristics of irregular binaries in gas solid fluidized beds——An ANN-approach

Binary mixtures of irregular materials of different particle sizes and/or particle densities are fluidized in a 15-cm diameter column with a perforated plate distributor. An attempt has been made in this work to determine the segregation characteristics of jetsam particles for both the homogeneous and heterogeneous binary mixtures in terms of segregation distance by correlating it to the various system parameters, viz. initial static bed height, height of a layer of particles above the bottom grid, superficial gas velocity and average particle size and/or particle densities of the mixture through the dimensional analysis. Correlation on the basis of Artificial Neural Network approach has also been developed with the above system parameters thereby authenticating the development of correlation by the former approach. The calculated values of the segregation distance obtained for both the homogeneous and heterogeneous binary mixtures from both the types of fluidized beds (i. e. under the static bed condition and the fluidized bed condition) have also been compared with each other.     

Transport modeling of sedimenting particles in a turbulent pipe flow using Euler–Lagrange large eddy simulation

A volume-filtered Euler–Lagrange large eddy simulation methodology is used to predict the physics of turbulent liquid–solid slurry flow through a horizontal periodic pipe. A dynamic Smagorinsky model based on Lagrangian averaging is employed to account for the sub-filter scale effects in the liquid phase. A fully conservative immersed boundary method is used to account for the pipe geometry on a uniform cartesian grid. The liquid and solid phases are coupled through volume fraction and momentum exchange terms. Particle–particle and particle–wall collisions are modeled using a soft-sphere approach. Three simulations are performed by varying the superficial liquid velocity to be consistent with the experimental data by Dahl et al. (2003). Depending on the liquid flow rate, a particle bed can form and develop different patterns, which are discussed in light of regime diagrams proposed in the literature. The fluctuation in the height of the liquid-bed interface is characterized to understand the space and time evolution of these patterns. Statistics of engineering interest such as mean velocity, mean concentration, and mean streamwise pressure gradient driving the flow are extracted from the numerical simulations and presented. Sand hold-up calculated from the simulation results suggest that this computational strategy is capable of predicting critical deposition velocity.     

Segregation characteristics of irregular binaries in gas solid fluidized beds—An ANN-approach

Binary mixtures of irregular materials of different particle sizes and/or particle densities are fluidized in a 15-cm diameter column with a perforated plate distributor. An attempt has been made in this work to determine the segregation characteristics of jetsam particles for both the homogeneous and heterogeneous binary mixtures in terms of segregation distance by correlating it to the various system parameters, viz. initial static bed height, height of a layer of particles above the bottom grid, superficial gas velocity and average particle size and/or particle densities of the mixture through the dimensional analysis. Correlation on the basis of Artificial Neural Network approach has also been developed with the above system parameters thereby authenticating the development of correlation by the former approach. The calculated values of the segregation distance obtained for both the homogeneous and heterogeneous binary mixtures from both the types of ftuidized beds （i.e. under the static bed condition and the ftuidized bed condition） have also been compared with each other.     

Direct numerical simulation of particle behaviour in the wall region of turbulent flows in horizontal channels

The motion of small particles in the wall region of turbulent channel flows has been investigated using direct numerical simulation. It is assumed that the particle concentration is low enough to allow the use of one-way coupling in the calculations, i.e. the fluid moves the particles but there is no feedback from the particles on the fluid motion. The velocity of the fluid is calculated by using a pseudospectral, direct solution of the Navier-Stokes equations. The calculations indicate that particles tend to segregate into the low-speed regions of the fluid motion near the wall. The segregation tendency depends on the time constant of the particle made non-dimensional with the wall shear velocity and kinematic viscosity. For very small and very large time constants, the particles are distributed more uniformly. For intermediate time constants (of the order 3), the segregation into the low-speed fluid regions is the highest. The finding that segregation occurs for a range of particle time constants is supported by experimental results. The findings regarding the more uniform distributions, however, still remain to be verified against experimental data which is not yet available. For horizontal channel flows, it is also found that particles are resuspended by ejections (of portions of the low-speed streaks) from the wall and are, therefore, primarily associated with low-speed fluid. The smaller particles are flung further upwards and, as they fall back towards the wall, they tend to be accelerated close to the fluid velocity. The larger particles have greater inertia and, consequently, accelerate to lower velocities giving higher relative velocities. This velocity difference, as a function of wall-normal distance, follows the same trend as in experiments but is always somewhat smaller in the calculations. This appears to be due to the Reynolds number for the numerical simulation being smaller than that in the experiment. It is concluded that the average particle velocity depends not only on the wall variables for scaling, but also on outer variables associated with the mean fluid velocity and fluid depth in the channel. This is because fluid depth in combination with the wall shear velocity determines how much time a particle, of a given size and density, spends in the outer flow and, hence, how close it gets to the local fluid velocity.     

Onset velocity of circulating fluidization and particle residence time distribution:A CFD-DEM study

Until now,the onset velocity of circulating fluidization in liquid-solid fluidized beds has been defined by the turning point of the time required to empty a bed of particles as a function of the superficial liquid velocity,and is reported to be only dependent on the liquid and particle properties.This study presents a new approach to calculate the onset velocity using CFD-DEM simulation of the particle residence time distribution(RTD).The onset velocity is identified from the intersection of the fitted lines of the particle mean residence time as a function of superficial liquid velocity.Our results are in reasonable agreement with experimental data.The simulation indicates that the onset velocity is influenced by the density and size of particles and weakly affected by riser height and diameter.A power-law function is proposed to correlate the mean particle residence time with the superficial liquid velocity.The collisional parameters have a minor effect on the mean residence time of particles and the onset velocity,but influence the particle RTD,showing some humps and trailing.The particle RTD is found to be related to the particle trajectories,which may indicate the complex flow structure and underlying mechanisms of the particle RTD.     

Onset velocity of circulating fluidization and particle residence time distribution: A CFD–DEM study

Until now, the onset velocity of circulating fluidization in liquid–solid fluidized beds has been defined by the turning point of the time required to empty a bed of particles as a function of the superficial liquid velocity, and is reported to be only dependent on the liquid and particle properties. This study presents a new approach to calculate the onset velocity using CFD–DEM simulation of the particle residence time distribution (RTD). The onset velocity is identified from the intersection of the fitted lines of the particle mean residence time as a function of superficial liquid velocity. Our results are in reasonable agreement with experimental data. The simulation indicates that the onset velocity is influenced by the density and size of particles and weakly affected by riser height and diameter. A power-law function is proposed to correlate the mean particle residence time with the superficial liquid velocity. The collisional parameters have a minor effect on the mean residence time of particles and the onset velocity, but influence the particle RTD, showing some humps and trailing. The particle RTD is found to be related to the particle trajectories, which may indicate the complex flow structure and underlying mechanisms of the particle RTD.     

Three-dimensional multi-phase simulation of the mixing and segregation of binary particle mixtures in a two-jet spout fluidized bed

This study presents a three-dimensional numerical study of the mixing and segregation of binary particle mixtures in a two-jet spout fluidized bed based on an Eulerian–Eulerian three-fluid model. Initially, the particle mixtures were premixed and packed in a rectangular fluidized bed. As the calculation began, the gas stream was injected into the bed from the distributor and jet nozzles. The model was validated by comparing the simulated jet penetration depths with corresponding experimental data. The main features of the complex gas–solid flow behaviors and the mechanism of mixing and segregation of the binary mixtures were analyzed. Moreover, further simulations were carried out to evaluate the effects of operating conditions on the mixing and segregation of binary particle mixtures. The results illustrate that mixing can be enhanced by increasing the jet velocity or enlarging the difference of initial proportions of binary particle mixtures.     

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.     

Effect of magnetic field orientation on fluidized beds of magnetic particles: Theory and experiment

Geldart-A fluidized beds of fine particles experience a jamming transition between a fluid-like state and a solid-like state at a certain superficial gas velocity, that depends on the relative strength of interparticle attractive forces with respect to particle weight. Interparticle forces provide the bed with a certain tensile strength in the jammed state. In the work presented here we analyze the behavior of a fluidized bed of magnetic particles subjected to an externally applied magnetic field, which contributes to enhance interparticle forces. The importance of the magnetic contribution to interparticle forces is measured by the changes in the tensile strength and the superficial gas velocity at the jamming transition. The link of the field orientation with the microstructure of the bed is discussed.     

Effect of magnetic field orientation on fluidized beds of magnetic particles： Theory and experiment

Geldart-A fluidized beds of fine particles experience a jamming transition between a fluid-like state and a solid-like state at a certain superficial gas velocity, that depends on the relative strength of interparticle attractive forces with respect to particle weight, lnterparticle forces provide the bed with a certain tensile strength in the jammed state. In the work presented here we analyze the behavior of a fluidized bed of magnetic particles subjected to an externally applied magnetic field, which contributes to enhance interparticle forces. The importance of the magnetic contribution to interparticle forces is measured by the changes in the tensile strength and the superficial gas velocity at the jamming transition. The link of the field orientation with the microstructure of the bed is discussed,     

Prediction of segregation characterization based on granular velocity and concentration in rotating drum

In a binary granular system composed of two types of particles with different granule sizes and the same density, particle sorting occurs easily during the flow process. The segregation pattern structure is mainly affected by the granular velocity and granular concentration in the flow layer. This paper reports on the experimental velocity and concentration measurement results for spherical particles in a quasi-two-dimensional rotating drum. The relationship between the granular velocity along the depth direction of the flow layer and granular concentration was established to characterize structures with different degrees of segregation. The corresponding relationships between the granular velocity and concentration and the segregation pattern were further analyzed to improve the theoretical models of segregation (convection–diffusion model and continuous flow model) and provide a reference for granular segregation control in the production process.