Euler-Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method(DEM)using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model.The model takes into account detailed gas phase chemistry,zero-dimensional modeling of the pyrolysis and gasification of each individual particle,particle shrinkage,and heat and mass transfer between the gas phase and the particulate phase.We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.     

DEM simulation of particle percolation in a packed bed

The phenomenon of spontaneous particle percolation under gravity is investigated by means of the discrete element method. Percolation behaviors such as percolation velocity,residence time distribution and radial dispersion are examined under various conditions. It is shown that the vertical velocity of a percolating particle moving down through a packing of larger particles decreases with increasing the restitution coefficient between particles and diameter ratio of the percolating to packing particles. With the increase of the restitution coefficient,the residence time and radial dispersion of the percolating particles increase. The packing height affects the residence time and radial dispersion. But,the effect can be eliminated in the analysis of the residence time and radial dispersion when they are normalized by the average residence time and the product of the packing height and packing particle diameter,respectively.In addition,the percolation velocity is shown to be related to the vertical acceleration of the percolating particle when an extra constant vertical force is applied. Increasing the feeding rate of percolating particles decreases the dispersion coefficient.     

DEM simulation of particle percolation in a packed bed

The phenomenon of spontaneous particle percolation under gravity is investigated by means of the discrete element method. Percolation behaviors such as percolation velocity, residence time distribution and radial dispersion are examined under various conditions. It is shown that the vertical velocity of a percolating particle moving down through a packing of larger particles decreases with increasing the restitution coefficient between particles and diameter ratio of the percolating to packing particles. With the increase of the restitution coefficient, the residence time and radial dispersion of the percolating particles increase. The packing height affects the residence time and radial dispersion. But, the effect can be eliminated in the analysis of the residence time and radial dispersion when they are normalized by the average residence time and the product of the packing height and packing particle diameter, respectively. In addition, the percolation velocity is shown to be related to the vertical acceleration of the percolating particle when an extra constant vertical force is applied. Increasing the feeding rate of percolating particles decreases the dispersion coefficient.     

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.     

Particle convective heat transfer near the wall in a supercritical water fluidized bed by single particle model coupled with CFD-DEM

Supercritical water fluidized bed (SCWFB) is a promising reactor to gasify biomass or coal. Its optimization design is closely related to wall-to-bed heat transfer, where particle convective heat transfer plays an important role. This paper evaluates the particle convective heat transfer coefficient (h_pc) at the wall in SCWFB using the single particle model. The critical parameters in the single particle model which is difficult to get experimentally are obtained by the computational fluid dynamics-discrete element method (CFD-DEM). The contact statistics related to particle-to-wall heat transfer, such as contact number and contact distance, are also presented. The results show that particle residence time (τ), as the key parameter to evaluate h_pc, is found to decrease with rising velocity, while increase with larger thermal boundary layer thickness. τ follows a gamma function initially adopted in the gas–solid fluidized bed, making it possible to evaluate h_pc in SCWFB by a simplified single particle model. The theoretical predicted h_pc tends to increase with rising thermal gradient thickness at a lower velocity (1.5 U_mf), while first decreases and then increases at higher velocity (1.75 and 2 U_mf). h_pc occupies 30%–57% of the overall wall-to-bed heat transfer coefficient for a particle diameter of 0.25 mm. The results are helpful to predict the overall wall-to-bed heat transfer coefficient in SCWFB combined with a reasonable fluid convective heat transfer model from a theoretical perspective.     

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.     

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 GenlDLEST 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.     

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.     

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 devia- tion 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 Uo 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 flu- idized 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 Uk 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.     

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.     

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.     

Euler–Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor

We present an Euler–Lagrange method for the simulation of wood gasification in a bubbling fluidized bed. The gas phase is modeled as a continuum using the 2D Navier–Stokes equations and the solid phase is modeled by a Discrete Element Method (DEM) using a soft-sphere approach for the particle collision dynamic. Turbulence is included via a Large-Eddy approach using the Smagorinsky sub-grid model. The model takes into account detailed gas phase chemistry, zero-dimensional modeling of the pyrolysis and gasification of each individual particle, particle shrinkage, and heat and mass transfer between the gas phase and the particulate phase. We investigate the influence of wood feeding rate and compare exhaust gas compositions and temperature results obtained with the model against experimental data of a laboratory scale bubbling fluidized bed reactor.     

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.     

回转窑多粒径颗粒流动与偏析的DEM模拟与实验研究

为明晰回转窑内颗粒的运动行为及偏析机理,以绿豆、黄豆和黑豆为颗粒介质,依次对3种装填顺序下的颗粒流动过程进行离散元模拟与实验研究,以颗粒质量分数和平均粒度为判据,对颗粒偏析进行评价。结果表明,回转窑内颗粒流动区可分为自由滚落区、渗流呆滞区以及窑壁携带区,自由滚落区颗粒流速最大,而渗流呆滞区流速最小。窑内颗粒沿轴向输运过程发生径向偏析,形成夹层结构,小颗粒受渗流作用在渗流呆滞区中心形成内核,大粒径和中等粒径颗粒集中在自由滚落区和窑壁携带区。窑内颗粒力链分布不均匀,强力链分布于近窑壁区,弱力链分布于自由滚落区和渗流呆滞区,且渗流呆滞区力链细而密集。当窑头附近不同粒径颗粒存在轴向速度差时,颗粒在轴向发生掺混,并产生径向偏析。     

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.     

Circulation intensity and axial dispersion of non-cohesive solid particles in a V-blender via DEM simulation

In this study, discrete element method （DEM） was employed to simulate the movement of non-cohesive mono-dispersed particles in a V-blender along with particle-particle and particle-boundary interactions. To validate the model, DEM results were successfully compared to positron emission particle tracking （PEPT） data reported in literature. The validated model was then utilized to explore the effects of rotational speed and fill level on circulation intensity and axial dispersion coefficient of non-cohesive particles in the V-blender. The results showed that the circulation intensity increased with an increase in the rotational speed from 15 to 60 rpm. As the fill level increased from 20% to 46%, the circulation intensity decreased, reached its minimum value at a fill level of 34% for all rotational speeds, and did not change significantly at fill levels greater than 34%. The DEM results also revealed that the axial dispersion coefficient of particles in the V-blender was a linear function of the rotational speed. These trends were in good agreement with the experimentallv determined values reported bv previous researchers.     

Influence of cycle time distribution on coating uniformity of particles in a spray fluidized bed by using CFD-DEM simulations

Cycle Time Distribution (CTD) plays a critical role for determining uniformity of particle coating in spray fluidized beds. However, the CTD is influenced by both geometrical structure and operating conditions of fluidized bed. In this study, a spray fluidized bed of coating process is simulated by a comprehensive Computational Fluid Dynamics-Discrete Element Model (CFD-DEM). To achieve different behaviors of CTD, some modifications are designed on a pseudo-2D internally circulating fluidized bed, which traditionally composes of a high-velocity upward bed and low-velocity downward bed. These modifications include making the air distributor slope and/or laying a baffle in the downward bed. First, the CTD and evolution of particle size distribution under different bed structures are compared. The CTD directly influences the coating uniformity. By making the particles flowing along a parallel direction in the downward bed through the geometrical modifications, the CTD becomes narrower and the coating uniformity is significantly improved. Second, under the optimized bed structure, the influence of operating conditions on the coating uniformity is studied. Properly increasing the fluidization gas velocity and the fluidization gas temperature and reducing the liquid spray rate can improve the coating uniformity.     

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.     

Particle Measurement Sensor for in situ determination of phase structure of fluidized bed

Based on three-dimensional (3D) acceleration sensing, an intelligent particle spy capable of detecting, transferring, and storing data, is proposed under the name of Particle Measurement Sensor (PMS). A prototype 60-mm-dia PMS was tested to track its freefall in terms of velocity and displacement, and served as a particle spy in a fluidized bed delivering the in situ acceleration information it detects. With increasing superficial gas velocity in the fluidized bed, the acceleration felt by PMS was observed to increase. The variance of the signals, which reflect the fluctuation, increased at first, reaching a maximum at the gas velocity (Uc) which marks the transition from bubbling to turbulent fluidization. Through probability density distribution (PDD) analysis, the PDD peak can be divided into the emulsion phase peak and the bubble phase peak. The average acceleration of emulsion and bubble phase increased, while the variance of both phases reached a maximum at Oc, at the same time. However, the difference between the variances of two phases reached the maximum at Oc. Findings of this study indicate that PMS can record independent in situ information. Further, it can provide other in situ measurements when equipped with additional multi-functional sensors.     

Simulation-supported measurements in large circulating fluidized bed combustors

The concept of simulation-supported measurement is suggested for the elucidation of processes occurring in the combustion chambers of large-scale circulating fluidized bed combustors where the desired information cannot be obtained by direct measurements.The concept is illustrated with the example of secondary air injection where the way the air is released,the penetration depth and the evenness of air distribution over the cross-sectional area of the combustion chamber are of interest.The measured informat...