Statistical and frequency analysis of the pressure fluctuation in a fluidized bed of non-spherical particles

In this paper, the pressure fluctuation in a fluidized bed was measured and processed via standard deviation and power spectrum analysis to investigate the dynamic behavior of the transition from the bubbling to turbulent regime. Two types (Geldart B and D) of non-spherical particles, screened from real bed materials, and their mixture were used as the bed materials. The experiments were conducted in a semi-industrial testing apparatus. The experimental results indicated that the fluidization characteristics of the non-spherical Geldart D particles differed from that of the spherical particles at gas velocities beyond the transition velocity U_c. The standard deviation of the pressure fluctuation measured in the bed increased with the gas velocity, while that measured in the plenum remained constant. Compared to the coarse particles, the fine particles exerted a stronger influence on the dynamic behavior of the fluidized bed and promoted the fluidization regime transition from bubbling toward turbulent. The power spectrum of the pressure fluctuation was calculated using the auto-regressive (AR) model; the hydrodynamics of the fluidized bed were characterized by the major frequency of the power spectrum of the pressure fluctuation. By combining the standard deviation analysis, a new method was proposed to determine the transition velocity U_k via the analysis of the change in the major frequency. The first major frequency was observed to vary within the range of 1.5 to 3 Hz.     

Characteristics of pressure fluctuations and fine coal preparation in gas-vibro fluidized bed

To improve the separation efficiency of air dense medium fluidized beds for dry coal preparation, a gas-vibro fluidized bed has been proposed in which magnetic powder is used as the heavy medium. Pressure fluctuations in the gas-vibro fluidized bed were investigated using time- and frequency-domain analysis methods. The relationship between pressure fluctuations, bubble behavior, and separation efficiency was established. The low amplitude of the standard deviation, the power spectral density (PSD), the incoherent-output PSD, and the high amplitude of the coherent-output PSD, which corresponds to the bubble behavior in the bed, were improved for coal preparation. The coal ash content was reduced from 42.55% to 16.54% by using the gas-vibro fluidized bed.     

DEM simulation of polydisperse systems of particles in a fluidized bed

Numerical simulations based on three-dimensional discrete element model (DEM) are conducted for mono-disperse, binary and ternary systems of particles in a fluidized bed. Fluid drag force acting on each particle depending on its size and relative velocity is assigned. The drag coefficient corresponding to Ergun’s correlation is applied to the system of fluidized bed with particle size ratios of 1:1 for the mono-disperse system, 1:1.2, 1:1.4 and 1:2 for the binary system and 1:1.33:2 for the ternary system b...     

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.     

Bubbling behavior of cohesive particles in a two-dimensional fluidized bed with immersed tubes

Fluidization hydrodynamics are greatly influenced by inter-particle cohesive forces. This paper studies the fluidization of large cohesive particles in a two-dimensional fluidized bed with immersed tubes using “polymer coating” to introduce cohesive force, to gain better understanding of bubbling behavior when particles become cohesive and its effect on chemical processes. The results show that the cohesive force promotes bubble splitting in the tube bank region, thereby causing an increase in the number and a decline in the aspect ratio of the bubbles. As the cohesive force increases within a low level, the bubble number increases and the bubble diameter decreases, while the aspect ratio exhibits different trends at different fluidization gas velocities. The difference in the evolution of bubble size under various cohesive forces mainly takes place in the region without tubes. When the cohesive force is large enough to generate stable agglomerates on the side walls of the bed, the bubble number and the bed expansion sharply decrease. The tubes serve as a framework that promotes the agglomeration, thus accelerating defluidization. Finally, the bubble profile around tubes was studied and found to greatly depend both on the cohesive forces and the location of tubes.     

Optimization of multi-stage velocity gradients in a cylindrical fluidized bed flocculator

Flocculation time is conventionally believed to be proportional to the flocculation efficiency of a cylindrical fluidized bed flocculator. However, in a single-stage velocity gradient situation, the flocculation efficiency decreases when the optimal flocculation time is exceeded. A multi-stage velocity gradient was established in a cylindrical fluidized bed flocculator, based on the hydraulic classification theory. This multi-stage velocity gradient fluidized bed flocculator (MGF) created a more suitable environment for floc growth and protection, which was confirmed by the size distribution of flocs along the bed height. Correspondingly, the abatement efficiencies for Kaolin slurry and dyed wastewater treatment in the MGF were enhanced by 5–10%, and by 7–20%, respectively, compared with those in the single-stage velocity gradient fluidized bed flocculators (SGFs). The initial bed height distribution ratio along the velocity gradients was an important factor for MGF optimization.     

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 analysis of volatile liquid injection into a fluidized bed

Experiments were conducted on a lab-scale fluidized bed to study the distribution of liquid ethanol injected into fluidized catalyst particles. Electrical capacitance measurements were used to study the liquid distribution inside the bed, and a new method was developed to determine the liquid content inside fluidized beds of fluid catalytic cracking particles. The results shed light on the complex liquid injection region and reveal the strong effect of superficial gas velocity on liquid distribution inside the fluidized bed, which is also affected by the imbibition of liquid inside particle pores. Particle internal porosity was found to play a major role when the changing mass of liquid in the bed was monitored. The results also showed that the duration of liquid injection affected liquid–solid contact inside the bed and that liquid–solid mixing was not homogeneous during the limited liquid injection time.     

Fluidization behavior of graphitized glassy particles in a fluidized carbon bed cooling process for investment casting

Fluidized Carbon Bed Cooling (FCBC) is an innovative investment casting process for directional solidification of superalloy components. It takes advantage of a fluidized bed with a base of small glassy carbon beads for cooling and other low-density particles that form an insulating layer by floating to the bed surface. This so-called “Dynamic Baffle” protects the fluidized bed from the direct heat input from the high-temperature heating zone and provides the basis for an improved bed microstructure. The prerequisites for a stable casting process are stable fluidization conditions where neither collapse of the bed nor particle blow out at excessive bubble formation occur.This work aimed to investigate the fluidization behavior of spherical carbon bed material in argon and air at temperatures between 20 to 350 °C. Systematic studies at reduced pressures using the FCBC prototype device were performed to understand the stable fluidization conditions at all stages of the investment casting process. The particle shape factor and size distribution characterization and the measurement of the powder’s minimum fluidization velocity and bed voidage show that this material can be fully utilized as a cooling and buoyancy medium during the FCBC process.     

Bubble characteristics in a bubbling fluidized bed with a rotating distributor

In this paper, the effect of a novel rotating distributor for fluidized beds on the bubble size is studied. The distributor is a perforated plate that rotates around the vertical axis of the column.The formation of the bubbles on the rotating distributor is theoretically analyzed. The pierced length of the bubbles ascending in the bed were measured using optical probes. The probability distribution of bubble diameter was inferred from these experimental measurements using the maximum entropy method. The radial profile of the bubble diameter is presented for the static and rotating configurations at different gas velocities. The frequency of bubble passage and the distribution of bubbles in the cross section of the bed are also reported. Results were finally shown for different heights above the distributor.A radial decrease in the bubble size when the distributor rotates is found. The bubble growth with the bed height is also lower in the rotating case.     

Effect of static bed height in the upper fluidized bed on flow behavior in the lower riser section of a coupled reactor

To study olefin reduction by using an auxiliary reactor for FCC naphtha upgrading, a large-scale cold model of a riser-bed coupled to an upper fluidized bed was established. The effect of static bed height in the upper fluidized bed on particle flow behavior in the lower riser was investigated experimentally. A restriction index of solids holdup was used to evaluate quantitatively the restrictive effect of the upper fluidized bed. Experimental results show that, under the restrictive effect of the upper fluidized bed, the riser could be divided into three regions in the longitudinal direction: accelerating, fully developed and restriction. The axial distribution of solids holdup in the riser is characterized by large solids holdup in the top and bottom sections and small solids holdup in the middle section. Overall solids holdup increased with increasing static bed height in the upper fluidized bed, while particle velocity decreased. Such restrictive effect of the upper fluidized bed could extend from the middle and top sections to the whole riser volume when riser outlet resistance is increased, which increases with increasing static bed height in the upper fluidized bed. The upper bed exerts the strongest restriction on the area close to the riser outlet.     

Effect of static bed height in the upper fluidized bed on flow behavior in the lower riser section of a coupled reactor

To study olefin reduction by using an auxiliary reactor for FCC naphtha upgrading, a large-scale cold model of a riser-bed coupled to an upper fluidized bed was established. The effect of static bed height in the upper fluidized bed on narticle flow behavior in the lower riser was investigated experimentally. A restriction index of solids holdup was used to evaluate quantitatively the restrictive effect of the upper fluidized bed. Experimental results show that, under the restrictive effect of the upper fluidized bed, the riser could be divided into three regions in the longitudinal direction： accelerating, fully developed and restriction. The axial distribution of solids holdup in the riser is characterized by large solids holdup in the top and bottom sections and small solids holdup in the middle section. Overall solids holdup increased with increasing static bed height in the upper fluidized bed, while particle velocity decreased. Such restrictive effect of the upper fluidized bed could extend from the middle and top sections to the whole riser volume when riser outlet resistance is increased, which increases with increasing static bed height in the upper fluidized bed. The upper bed exerts the strongest restriction on the area close to the riser outlet.     

Application of CFD-DEM to the study of solid exchange in a dual-leg fluidized bed

The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtain the dominant frequency （F） of the simulated differential particle number （APLR） between the two half beds. Effects of fluidization velocity （u） and bed material inventory （H） on the solid exchange behavior were studied using the CFD-DEM model. Not only snapshots of the simulated particle flow patterns using the OpenGL code but also the dominant frequency of APLR was similar to the experimental results. The simulation results show that higher fluidization velocity assists the exchange of more particles between the two half beds, but the dispersion of clusters on the bed surface into single particles decreases the cluster exchange frequency. A greater bed material inventory results in more intense cluster exchange. The cluster exchange frequency decreases with an increase of the bed material inventory.     

Comparison of the standard Euler-Euler and hybrid Euler-Lagrange approaches for modeling particle transport in a pilot-scale circulating fluidized bed

Particle transport phenomena in small-scale circulating fiuidized beds （CFB） can be simulated using the Euler-Euler, discrete element method, and Euler-Lagrange approaches. In this work, a hybrid Euler-Lagrange model known as the dense discrete phase model （DDPM）, which has common roots with the multiphase particle-in-cell model, was applied in simulating particle transport within a mid-sized experimental CFB facility. Implementation of the DDPM into the commercial ANSYS Fluent CFD package is relatively young in comparison with the granular Eulerian model. For that reason, validation of the DDPM approach against experimental data is still required and is addressed in this paper. Additional difficulties encountered in modeling fluidization processes are connected with long calculation times. To reduce times, the complete boiler models are simplified to include just the combustion chamber. Such simplifications introduce errors in the predicted solid distribution in the boiler. To investigate the conse- quences of model reduction, simulations were made using the simplified and complete pilot geometries and compared with experimental data. All simulations were performed using the ANSYSFLUENT 14.0 package. A set of user defined functions were used in the hybrid DDPM and Euler-Euler approaches to recirculate solid particles.     

Development and assessment of algorithms for DEM-LES simulations of fluidized bed

The use of high-fidelity Discrete Element Method (DEM) coupled with Computational Fluid Dynamics (CFD) for particle-scale simulations demands extensive simulation times and restricts application to small particulate systems. DEM-CFD simulations require good performance and satisfactory scalability on high-performance computing platforms. A reliable parallel computing strategy must be developed to calculate the collision forces, since collisions can occur between particles that are not on the same processor, or even across processors whose domains are disjoint. The present paper describes a parallelization technique and a numerical verification study based on a number of tests that allow for the assessment of the numerical performance of DEM used in conjunction with Large-Eddy Simulation (LES) to model dense flows in fluidized beds. The fluid phase is computed through solving the volume-averaged four-way coupling Navier-Stokes equations, in which the Smagorinsky sub-grid scale tensor model is used. Furthermore, the performance of Sub-Grid Scale (SGS) turbulence models applied to Fluidized Bed Reactor (FBR) configurations has been assessed and compared. The developed numerical solver represents an interesting combination of techniques that work well for the present purpose of studying particle formation in fluidized beds.     

Three-dimensional dynamic characterization of square-nosed slugging phenomena in a fluidized bed

Slugging represents one of the major regimes in fluidization, which occurs in small diameter beds with large bed height-to-diameter ratio or in large diameter beds with internals that resemble multiple small diameter fluidized beds. Slug types include round-nosed slug, wall slug and square-nosed slug. Studies of the slugs have been mainly focused on round-nosed or wall slugs known as half slug, typically occurring in Geldart group A particle fluidization. The square-nosed slug typically occurring for Geldart group D particles appears to be regarded as simple in its structure. The Electrical Capacitance Volume Tomography (ECVT) imaging of the square-nosed slugging phenomena conducted in this study reveals otherwise. That is the structure of the square-nosed slug is, in fact, complex, particularly with respect to its dynamic variation in fluidization. More broadly, this study examines experimentally the hydrodynamic characteristics of the square-nosed fluidization regime. Specifically, simultaneous measurements from multiple ECVT sensors provide non-invasive, continuous, 3-dimensional imaging of the entire flow region of the slugging bed and hence enabling the dynamic characterization of the evolution of the slugs. The analysis of the 3D images reconstructed for real-time gas–solid volume fraction profile of the slugging fluidized bed indicates that there are three different zones, namely, the bottom fluidization zone, the gas slug zone, and the solid slug zone, co-existing in the bed. The three zones present different hydrodynamic characteristics during the slug evolution. It is found that varying the gas velocity of the slugging bed mainly varies the maximum length of the gas slug zone, while it only has a minor effect on the lengths of the bottom fluidization zone and solid slug zone. It also has an insignificant effect on the solid volume fraction of the three zones.     

Rotation characteristic and granular temperature analysis in a bubbling fluidized bed of binary particles

In this work, a discrete particle model (DPM) was applied to investigate the dynamic characteristics in a gas–solid bubbling fluidized bed of binary solid particles. The solid phase was simulated by the hard-sphere discrete particle model. The large eddy simulation (LES) method was used to simulate the gas phase. To improve the accuracy of the simulation, an improved sub-grid scale (SGS) model in the LES method was also applied. The mutative Smagorinsky constant case was compared with the previously published experimental data. The simulation by the mutative Smagorinsky constant model exhibited better agreement with the experimental data than that by the common invariant Smagorinsky constant model. Various restitution coefficients and different compositions of binary solids were investigated to determine their influences on the rotation characteristics and granular temperatures of the particles. The particle translational and rotational characteristic distributions were related to certain simulation parameters.     

DEM investigation on conveying of non-spherical particles in a screw conveyor

Screw conveyors are extensively used in modern industry such as metallurgy, architecture and pharmaceutical due to their high-efficiency in the transportation of granular materials. And substantial efforts have been devoted to the study of the screw conveyors. Numerical method is an effective way to study screw conveyor. However, previous studies have mainly focused in the regime of spherical particles while the in-depth investigations for non-spherical particles that should be the most encountered in practical applications are still limited. In view of the above situations, discrete element method (DEM), which has been widely accepted in simulating the discrete systems, is utilized to investigate the conveying process of non-spherical particles in a horizontal screw conveyor, with particles being modeled by super-ellipsoids. In addition, a wear model called SIEM (Shear Impact Energy Model) is incorporated into DEM to predict the wear of screw conveyor. The DEM simulation results demonstrate that the particle shape is influential for the flow behaviors of particles and the wear of conveyor. The conveying performance evaluated quantitatively of both mass flow rate and power consumption is subsequently obtained to investigate the effect of sphericity of particle with different operation parameters. Moreover, particle collision frequency and collision energy consumption are acquired to investigate the possible particle breakage between particles and screw blade. The comparisons between particle–particle collision and particle–wall collision reveal that particles with large shape index have more possibility to be damaged in particle–wall impingement.     

Hydrodynamics of gas-solid fluidization of a homogeneous ternary mixture in a conical bed： Prediction of bed expansion and bed fluctuation ratios

Hydrodynamic characteristics of fluidization in a conical or tapered bed differ from those in a columnar bed because the superficial velocity in the bed varies in the axial direction. Fixed and fluidized regions could coexist and sharp variations in pressure drop could occur, thereby giving rise to a noticeable pressure drop-flow rate hysteresis loop under incipient fluidization conditions. To explore these unique properties, several experiments were carried out using homogeneous, well-mixed, ternary mixtures with three dif- ferent particle sizes at varying composition in gas-solid conical fluidized beds with varying cone angles. The hydrodynamic characteristics determined include the minimum fluidization velocity, bed fluctuation, and bed expansion ratios. The dependence of these quantities on average particle diameter, mass fraction of the fines in the mixture, initial static bed height, and cone angle is discussed. Based on dimensional analysis and factorial design, correlations are developed using the system parameters, i.e. geometry of the bed （cone angle）, particle diameter, initial static bed height, density of the solid, and superficial velocity of the fluidizing medium. Experimental values of minimum fluidization velocity, bed fluctuation, and bed expansion ratios were found to agree well with the developed correlations.