Two-fluid modeling of Geldart A particles in gas-fluidized beds

We have investigated the effect of cohesion and drag models on the bed hydrodynamics of Geldart A particles based on the two-fluid （TF） model. For a high gas velocity U0 = 0.03 m/s, we found a transition from the homogeneous fluidization to bubbling fluidization with an increase of the coefficient C1, which is used to account for the contribution of cohesion to the excess compressibility. Thus cohesion can play a role in the bed expansion of Geldart A particles. Apart from cohesion, we have also investigated the influence of the drag models. When using the Wen and Yu drag correlation with an exponent n = 4.65, we find an under-prediction of the bed expansion at low gas velocities （U0 = 0.009 m/s）. When using a larger exponent （n = 9.6）, as reported in experimental studies of gas-fluidization, a much better agreement with the experimental bed expansion is obtained. These findings suggest that at low gas velocity, a scale-down of the commonly used drag model is required. On the other hand, a scale-up of the commonly used drag model is necessary at high gas velocity （U0 = 0.2 and 0.06 m/s）. We therefore conclude that scaling the drag force represent only an ad hoc way of repairing the deficiencies of the TF model, and that a far more detailed study is required into the origin of the failure of the TF model for simulating fluidized beds of fine powders.     

Two-fluid modeling of Geldart A particles in gas–solid micro-fluidized beds

The fluidization behavior of Geldart A particles in a gas–solid micro-fluidized bed was investigated by Eulerian–Eulerian numerical simulation. The commonly used Gidaspow drag model was tested first. The simulation showed that the predicted minimum bubbling velocities were significantly lower than the experimental data even when an extremely fine grid size (of approximately one particle diameter) was used. The modified Gibilaro drag model was therefore tested next. The predicted minimum bubbling velocity and bed voidage were in reasonable agreement with the experimental data available in literature. The experimentally observed regime transition phenomena from bubbling to slugging were also reproduced successfully in the simulations. Parametric studies indicated that the solid-wall boundary conditions had a significant impact on the predicted gas and solid flow behavior.     

Development of a filtered drag model considering effect of the solid shear rate

This study investigates a 2D gas–solid fast fluidized bed of typical Geldart A particles using highly resolved simulations with two-fluid model. The results show that the solid shear rate has a considerable impact on the orientation of the meso-scale structures and hence on the filtered drag force. On the basis of the correlation for the filtered drag force established in the literature using the traditional markers (such as filtered solid volume fraction, filtered slip velocity and filter scale), a correction correlated with the solid shear rate in the direction of gravity is proposed for better prediction of the filtered drag force. The corrected model is shown to produce improved results in posterior tests of flows in different fluidization regimes including bubbling, turbulent and fast fluidization.     

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.     

鼓泡流化床中流动特性的多尺度数值模拟

鼓泡流化床因其较高的传热特性以及较好的相间接触已经被广泛应用于工业生产中,而对鼓泡流态化气固流动特性的充分认知是鼓泡流化床设计的关键.在鼓泡流化床中,气泡相和乳化相的同时存在使得床中呈现非均匀流动结构,而这种非均匀结构给鼓泡流化床的数值模拟造成了很大的误差.基于此,以气泡作为介尺度结构,建立了多尺度曳力消耗能量最小的稳定性条件,构建了适用于鼓泡流化床的多尺度气固相间曳力模型.结合双流体模型,对A类和B类颗粒的鼓泡流化床中气固流动特性进行了模拟研究,分析了气泡速度、气泡直径等参数的变化规律.研究表明,与传统的曳力模型相比,考虑气泡影响的多尺度气固相间曳力模型给出的曳力系数与颗粒浓度的关系是一条分布带,建立了控制体内曳力系数与局部结构参数之间的关系.通过模拟得到的颗粒浓度和速度与实验的比较可以发现,考虑气泡影响的多尺度曳力模型可以更好地再现实验结果.通过A类和B类颗粒的鼓泡床模拟研究发现,A类颗粒的鼓泡床模拟受多尺度曳力模型的影响更为显著.      

Two-fluid modeling of Geldart A particles in gas-solid micro-fluidized beds

The fluidizadon behavior of Geldart A particles in a gas-solid micro-fluidized bed was investigated by Eulerian-Eulerian numerical simulation.The commonly used Gidaspow drag model was tested first.The simulation showed that the predicted minimum bubbling velocities were significantly lower than the experimental data even when an extremely fine grid size(of approximately one particle diameter) was used.The modified Gibilaro drag model was therefore tested next.The predicted minimum bubbling velocity and bed voidage were in reasonable agreement with the experimental data available in literature.The experimentally observed regime transition phenomena from bubbling to slugging were also reproduced successfully in the simulations.Parametric studies indicated that the solid-wall boundary conditions had a significant impact on the predicted gas and solid flow behavior.     

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.     

INTERPARTICLE FORCES IN HIGH TEMPERATURE FLUIDIZATION OF GELDART A PARTICLES

Previous reports and current studies show that fluidization of some Geldart A particles is enhanced by in-creasing bed temperature.Both the averaged local particle concentration and the particle concentration in the dense phase decrease with increasing bed temperature,at constant superficial gas velocities.However,conventional models fail to predict these changes,because the role of interparticle forces is usually neglected at different bed temperatures.Here.the temperature increases,the interparticle attractive forces decrease while the interparticle repulsive forces increase.Consequently.fluidization behaviors of some Geldart A particles seem to increasingly shift from typical Geldart A towards B with increasing temperature.     

Hydrodynamics of a bubble-induced inverse fluidized bed reactor with a nanobubble tray

This paper reports on the hydrodynamics of a bubble-induced inverse fluidized bed reactor, using a nanobubble tray gas distributor, where solid particles are fluidized only by an upward gas flow. Increasing the gas velocity, the fixed layer of particles initially packed at the top of the liquid starts to move downwards, due to the rise of bubbles in this system, and then gradually expands downwards until fully suspended. The axial local pressure drops and standard deviation were examined to delineate the flow regime comprehensively under different superficial gas velocities. Four flow regimes (fixed bed regime, initial fluidization regime, expanded regime, and post-homogeneous regime) were observed and three transitional gas velocities (the initial fluidization velocity, minimum fluidization velocity, and homogeneous fluidization velocity) were identified to demarcate the flow regime. Three correlations were developed for the three transitional velocities. As the fine bubbles generated from the nanobubble tray gas distributor are well distributed in the entire column, the bed expansion process of the particles is relatively steady.     

Flow structure and bubble dynamics in supercritical water fluidized bed and gas fluidized bed: A comparative study

Supercritical water (SCW) fluidized bed is a new reactor concept for hydrogen production from biomass or coal gasification. In this paper, a comparative study on flow structure and bubble dynamics in a supercritical water fluidized bed and a gas fluidized bed was carried out using the discrete element method (DEM). The results show that supercritical water condition reduces the incipient fluidization velocity, changes regime transitions, i.e. a homogeneous fluidization was observed when the superficial velocity is in the range of the minimum fluidization velocity and minimum bubbling velocity even the solids behave as Geldart B powders in the gas fluidized bed. Bubbling fluidization in the supercritical water fluidized bed was formed after superficial velocity exceeds the minimum bubbling velocity, as in the gas fluidized bed. Bubble is one of the most important features in fluidized bed, which is also the emphasis in this paper. Bubble growth was effectively suppressed in the supercritical water fluidized bed, which resulted in a more uniform flow structure. By analyzing a large number of bubbles, bubble dynamic characteristics such as diameter distribution, frequency, rising path and so on, were obtained. It is found that bubble dynamic characteristics in the supercritical water fluidized bed differ a lot from that in the gas fluidized bed, and there is a better fluidization quality induced by the bubble dynamics in the supercritical water fluidized bed.     

Vibration time series analysis of bubbling and turbulent fluidization

A non-intrusive vibration monitoring technique was used to study the hydrodynamics of a gas–solid fluidized bed. Experiments were carried out in a 15 cm diameter fluidized bed using 226, 470 and 700 μm sand particles at various gas velocities, covering both bubbling and turbulent regimes. Auto correlation function, mutual information function, Hurst exponent analysis and power spectral density function were used to analyze the fluidized bed hydrodynamics near the transition point from bubbling to turbulent fluidization regimes. The first pass of the autocorrelation function from one half and the time delay at which it becomes zero, and also the first minimum of the mutual information, occur at a higher time delay in comparison to stochastic systems, and the values of time delays were maximum at the bubbling to turbulent transition gas velocity. The maximum value of Hurst exponent of macro structure occurred at the onset of regime transition from bubbling to turbulent. Further increase in gas velocity after that regime transition velocity causes a decrease in the Hurst exponent of macro structure because of breakage of large bubbles to small ones. The results showed these methods are capable of detecting the regime transition from bubbling to turbulent fluidization conditions using vibration signals.     

LBM-DEM四向耦合喷动床内介尺度模拟与分析

研究喷动床内颗粒的流动特性对于喷动床的设计和优化具有重要意义。基于气固两相流流动的LBM-DEM四向耦合模型,对单孔射流喷动床中颗粒的流动进行数值模拟。其中,气相采用修正的格子玻尔兹曼方法,颗粒相采用离散单元法,流固之间受力的双向耦合基于牛顿第三定律,颗粒与颗粒及颗粒与壁面的受力双向耦合采用软球模型。模拟得到了流化过程、颗粒与气体的速度分布、床层膨胀高度变化以及床宽对流化过程的影响。结果表明,喷动床内存在强烈的内循环,床宽增加导致颗粒运动能力减弱,射流速度增加使颗粒运动更加剧烈,床层膨胀高度增加。     

IMPROVEMENT OF FLUIDIZABILITY OF FINE POWDERS ——A COMPUTER STUDY

Simulations of the gas fluidization of a cohesive powder were performed using the Stokesian Dynamics method and an agglomeration-deagglomeration model to investigate methods of improving the fluidizability of fine powders. Three techniques (a) high gas velocity (b) vibration-assisted fluidization and (c) tapered fluidizer were used in the simulations which provided detailed information on the bed microscopy such as the motion of 1 O0 particles in a fluidizing vessel along with the formation and destruction of cohesive bonds dudng collisions. While all three techniques were found to effectively improve the fluidizability of a strongly cohesive powder, we suggest a combination of high velocity fluidization assisted by extemal vibration of the fluidized bed to minimize entrainment of particles.     

Behavior of mixed ZnO and SiO2 nano-particles in magnetic field assisted fluidization

The fluidization behavior of ZnO nano-particles in magnetic fluidized bed （MFB） by adding coarse magnetic particles was investigated, followed by the co-fluidization of mixtures of ZnO and SiO2 nano-particles. For such co-fluidization, bed expansion was found to change smoothly with gas velocity through a range of stable operation. By measuring the bed expansion ratio and pressure drop, a stability diagram for the mixture in MFB was obtained. Within this stable operation range, with increasing gas velocity the pressure drop hardly changes as the bed expands, up to an expansion ratio of more than 4.     

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.     

A state-of-the-art review on transition between the fast fluidization and pneumatic transport of Geldart group B particles

The fluidization state in the circulating fluidized bed (CFB) boiler is crucial to its stable and safe operation. However, up to now, the research field has not reached unanimity on whether the fluidization regime that the upper furnace of the boiler operates in is the fast fluidization or pneumatic transport. To this end, this paper reviewed relevant research on the transition between the fast fluidization and pneumatic transport of Geldart group B particles, including the flow characteristics of the fast fluidization, the transition condition between the fast fluidization and pneumatic transport, the determination methods of the transport velocity u_tr and saturation carrying capacity G_s1 and the influencing factors on these two parameters. Previous research findings can provide certain guidelines for the design and optimization of the CFB boiler, and result in plenty of prediction correlations for u_tr and G_s1. Nonetheless, owing to insufficient data available on Geldart group B particles, especially the ones obtained under high temperature or pressure conditions and in large-scale CFB apparatuses, the existing correlations are not well suited for the prediction of u_tr and G_s1 of Geldart group B particles. Thus, further efforts are urgently demanded on the fast fluidization transition of Geldart group B particles.     

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.     

基于CFD-DEM的高炉风口回旋区形状和传热特性数值模拟研究

采用计算流体力学-离散单元法(CFD-DEM)耦合方法,对高炉风口回旋区进行了数值模拟研究。首先通过与实验结果对比,验证了CFD-DEM模型的正确性;然后考察了不同气速对风口回旋区形状和传热特性及颗粒接触的影响。数值模拟结果表明:风口回旋区的大小和形状均受气速影响较大,在较大进气速度下,颗粒受到的曳力大于颗粒间的摩擦阻力并破坏颗粒间的桥力,形成较大尺寸的回旋区;且颗粒之间接触力较小,形成较大的空隙结构,更有利于热气体向周围扩散以强化传热。目前考察的三种气速结果表明:当气速为11m/s时,热量向下方传递速度最快;当气速为13m/s时,热量向上方传递速度最快;而当气速为15m/s时,热量向右方传递速度最快;此外,气速越大流态化越明显,颗粒间接触越少,接触力也越小。     

Numerical simulation and experimental validation of gas–solid flow in the riser of a dense fluidized bed reactor

Gas–solid flow in the riser of a dense fluidized bed using Geldart B particles (sand), at high gas velocity (7.6–15.5 m/s) and with comparatively high solid flux (140–333.8 kg/m² s), was investigated experimentally and simulated by computational fluid dynamics (CFD), both two- and three-dimensional and using the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd and EMMS drag models. The results predicted by EMMS drag model showed the best agreement with experimental results. Calculated axial solids hold-up profiles, in particular, are well consistent with experimental data. The flow structure in the riser was well represented by the CFD results, which also indicated the cause of cluster formation. Complex hydrodynamical behaviors of particle cluster were observed. The relative motion between gas and solid phases and axial heterogeneity in the three subzones of the riser were also investigated, and were found to be consistent with predicted flow structure. The model could well depict the difference between the two exit configurations used, viz., semi-bend smooth exit and T-shaped abrupt exit. The numerical results indicate that the proposed EMMS method gives better agreement with the experimental results as compared with the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd models. As a result, the proposed drag force model can be used as an efficient approach for the dense gas–solid two-phase flow.