粗糙颗粒动理学及流化床内气固流动的数值模拟

基于气体分子动理学和颗粒动理学理论,考虑颗粒旋转流动对颗粒碰撞能量交换和耗散的影响,建立粗糙颗粒动理学。采用Chapman-Enskog颗粒速度分布函数,提出了颗粒相应力、热通量和颗粒碰撞能量耗散计算模型。采用欧拉-欧拉气固双相流模型,数值模拟鼓泡流化床内气体-颗粒两相流动特性。模拟结果得到了床内颗粒相速度和脉动速度分布,与Yuu等实验结果相吻合。分析不同的切向弹性恢复系数对颗粒相拟总温的变化规律,结果表明在低颗粒浓度时颗粒拟总温随切向弹性恢复系数而增加。     

颗粒旋转对鼓泡流化床内气固两相流动特性影响的数值模拟

运用考虑颗粒自旋转流动对颗粒碰撞能量交换和耗散影响的颗粒动理学方法，建立鼓泡流化床气固两相Euler-Euler双流体模型，数值模拟流化床内气体颗粒两相流动特性。分析表明，颗粒平动温度与旋转温度之比是法向和切向颗粒弹性恢复系数和摩擦系数的函数。与不考虑颗粒旋转效应计算结果相比，考虑颗粒旋转效应后床内较容易形成气泡，颗粒自旋转运动将导致床内非均匀结构更明显。并且床层平均空隙率和床层膨胀高度增加，床中心区域颗粒轴向速度提高，床内颗粒平动温度下降。考虑颗粒旋转效应后预测的颗粒轴向速度和颗粒脉动速度与文献实验结果基本吻合。考虑颗粒旋转效应后获得的气泡直径更接近于前人经验关联式。     

鼓泡床取热器流动特性的数值模拟

应用计算颗粒流体力学(CPFD)方法,采用点源注射的进气方式对气固鼓泡流化床取热器内流动特性进行数值模拟。考察不同气速下床层膨胀高度、轴径向时均固含率分布、颗粒轴向速度分布及床层颗粒内循环流率的分布。模拟结果与实验数据吻合较好,表明该模型可以用于描述鼓泡流化床取热器内气固两相的流动规律。     

基于DEM模拟的气固鼓泡床内流场间歇性及颗粒相干结构的分析

采用计算流体力学和离散元方法（CFD-DEM）对气固鼓泡床流动行为进行模拟研究,并基于模拟结果分析鼓泡床内气泡和颗粒微观运动特性。对颗粒速度的脉动能谱进行分析,发现鼓泡床流场中存在间歇性。通过对比鼓泡床不同轴、径向位置的颗粒脉动速度的平坦因子,发现鼓泡床内不同位置的流场间歇性不同,随着床层高度的增加,流场的间歇性减弱;在径向上,过渡区的流场间歇性明显大于边壁区和中心区。进一步采用连续小波分析方法揭示了相干结构（颗粒涡团）的分布以及演化过程,并分析了不同尺度下相干结构（颗粒涡团）的分布与鼓泡床内颗粒与气泡运动的关系。     

基于DEM模拟的气固鼓泡床内流场间歇性及颗粒相干结构的分析

采用计算流体力学和离散元方法(CFD-DEM)对气固鼓泡床流动行为进行模拟研究,并基于模拟结果分析鼓泡床内气泡和颗粒微观运动特性。对颗粒速度的脉动能谱进行分析,发现鼓泡床流场中存在间歇性。通过对比鼓泡床不同轴、径向位置的颗粒脉动速度的平坦因子,发现鼓泡床内不同位置的流场间歇性不同,随着床层高度的增加,流场的间歇性减弱;在径向上,过渡区的流场间歇性明显大于边壁区和中心区。进一步采用连续小波分析方法揭示了相干结构(颗粒涡团)的分布以及演化过程,并分析了不同尺度下相干结构(颗粒涡团)的分布与鼓泡床内颗粒与气泡运动的关系。     

三维单旋导向挡板鼓泡流化床内气固流动的CPFD模拟

单旋导向挡板是一种广泛应用于流化床反应器中的横向内构件,研究其对鼓泡流化床内气固流动的影响具有重要意义。采用计算颗粒流体力学(CPFD)方法耦合Igci曳力模型研究带有两层单旋导向挡板的鼓泡流化床内的三维气固流动,并与空筒鼓泡流化床作对比。首先将轴向和径向时均固含率分布的模拟结果和实验数据进行对比,验证了该计算方法的可靠性;随后考察了单旋导向挡板对鼓泡流化床内气泡平均当量直径沿床高的分布以及气体停留时间分布的影响,计算结果表明,单旋导向挡板对大气泡具有一定的破碎作用,在挡板影响区域内气泡平均当量直径变小。同时,有、无挡板流化床内的气体停留时间分布曲线与其数学特征基本一致,表明单旋导向挡板的存在对鼓泡流化床内的气体返混没有明显影响。     

鼓泡流化床内高灰煤气化过程的颗粒尺度特性研究

基于多相流体质点网格方法(MP-PIC)对高灰煤在三维鼓泡流化床气化过程进行了数值模拟研究。在欧拉-拉格朗日框架下将气相和固相分别视作连续介质和离散相处理。首先,将模拟得到的出口处气体组分结果与实验数据进行对比,实验数据与模拟结果具有良好的一致性。其次,研究了煤颗粒在气化炉内的温度、传热系数、速度和停留时间,从颗粒尺度揭示了鼓泡流化床气化炉内的颗粒分布特性和气固流动特征。结果表明:在气化炉入口附近煤颗粒与床层温差最大,传热系数最大;由于流化床内强非线性的气固流动,床中煤温度和传热系数的空间分布不均匀;煤颗粒和床料的瞬时速度具有稳定的波动幅度,其中垂直方向速度波动最明显,且煤颗粒的瞬时速度比床料的瞬时速度略大;由于颗粒间的剧烈碰撞,延长了煤颗粒停留时间。此外,对鼓泡流化床中煤气化过程颗粒尺度的研究,有助于深入了解固体颗粒的流动行为以及气固相相互作用特性,对鼓泡流化床反应器的设计优化具有重要意义。     

鼓泡流化床离散模拟中的一个局部空隙率模型

气固流化床离散颗粒模拟中通常采用面积加权平均法计算局部空隙率,不能较好地反映流化床中显著的非均匀结构特性.为了考虑非均匀结构对局部空隙率的影响,文中提出一个适用鼓泡流化床的局部空隙率模型.将流场划分为稀区(气泡区)和密区(乳相区);非均匀结构对密区局部空隙率的影响通过引入局部空隙率下限和非均匀影响系数描述;将提升管流动...     

静电效应对有无埋管气固鼓泡床内气泡特性的影响分析

采用基于双流体模型(TFM)耦合静电模型的方法,研究颗粒的静电对有无埋管气固鼓泡床内气固流动特性和气泡特性的影响。首先在无静电场存在的条件下,利用双流体模型对自由鼓泡床和埋管鼓泡床内的流动情况进行模拟并与实验结果进行对比;进一步耦合静电模型,考察静电对自由鼓泡床和埋管鼓泡床内床层的整体性质和气泡特性的影响。研究结果表明,在无静电场条件下采用双流体模型能较好地预测自由鼓泡床和埋管鼓泡床内的气固流动状况以及气泡的平均直径和气泡的上升速度。埋管的存在使鼓泡床内气固流动发生强烈扰动,并使气泡的平均直径和气泡的上升速度均呈振荡分布。静电的存在对自由鼓泡床和埋管鼓泡床内床层的平均固含率影响不大,但对气泡分布规律影响较大,使得自由鼓泡床内气泡数目减少,而埋管鼓泡床下部区域的气泡分布比较集中,上部有大气泡出现。     

基于DEM模拟气固鼓泡床中颗粒碰撞参数对流场间歇性的影响

采用计算流体力学和离散元(CFD-DEM)方法研究鼓泡床内的气固流动状态, 考察了颗粒弹性系数和恢复系数对流场间歇性的影响, 并利用小波变化分析方法分析了颗粒弹性系数和恢复系数对流场相干结构的影响。研究结果表明:颗粒弹性系数和恢复系数对颗粒速度脉动能、床层平均高度、平坦因子以及流场间歇性有一定影响。随着颗粒弹性系数取值的变大, 高频区能量和平坦因子先降低后增加, 床层高度先增加后降低, 流场间歇性先减弱后增强;颗粒恢复系数取值越大, 高频区能量和平坦因子越低, 床层高度越大, 流场间歇性越弱。     

Practical validation of the two-fluid model applied to dense gas-solid flows in fluidized beds

This paper discusses the simulation of bubbling gas-solid flows by using the Eulerian two-fluid approach. Predictions of particle motion, bed expansion, bubble size and bubble velocity in bubbling beds containing Geldart B particles are compared with experimental results and correlations found in the literature. In addition, gas mixing in a bed of Geldart A particles is investigated.An in-house code has been developed based on the finite-volume method and the time-splitting approach using a staggered grid arrangement. The velocities in both phases are obtained by solving the 2D Reynolds-averaged Navier/Stokes equations using a partial elimination algorithm (PEA) and a coupled solver. The k-ε turbulence model is used to describe the turbulent quantities in the continuous phase.In general, the model predictions are in good agreement with experimental data found in the literature. Most important observations are: the level of the restitution coefficient was found to be crucial in order to obtain successful results from 2D axisymmetric simulations of a system containing Geldart B particles. Bubble size and bubble rise velocities are not as sensitive to the restitution coefficient. The turbulence model is of outmost importance concerning gas mixing in a fluidized bed of Geldart A particles.From these numerical analyzes an optimized granular flow two-fluid model can be designed for the purpose of simulating reactive systems in fluidized bed reactors.     

CFD simulation of gas-solid bubbling fluidized bed containing FCC particles

The hydrodynamics of a bubbling gas-solid fluidized bed of 57.4 μm FCC particles was simulated by using a state-of-the-art two-fluid model integrating the kinetic theory of granular flow for particulate phase stresses. The overestimation of the bed expansion was resolved by using a suitable scale factor in the drag model as suggested by McKeen and Pugsley (T.R. McKeen, T.S. Pugsley, Powder Technol., 129, 139 (2003)). This study showed that the method was appropriate in simulation of a gas-solid fluidized bed of Geldart A particles at high gas velocities (0.3 to 0.61 m/s). The reduction of computational time especially for simulation of large-scale systems was achieved. The time-averaged local voidage was compared with the experimental data and the trend of varying several parameters on the hydrodynamic of the bed was investigated. The simulation results showed both qualitative and quantitative agreement with the literature.     

Numerical simulation of gas-particle flow with a second-order moment method in bubbling fluidized beds

Flow behavior of gas and particles is performed by means of gas-solid two-fluid model with the second-order moment model of particles in the bubbling fluidized bed. The distributions of velocity and moments of particles are predicted in the bubbling fluidized beds. Predictions are compared with experimental data measured by Jung et al. (2005) in a bubbling fluidized bed and Patil et al. (2005) experiments in a bubbling fluidized bed with a jet. The simulated second-order moment in the vertical direction is on average 1.5-2.3 times larger than that in the lateral direction in the bubbling fluidized bed (Jung et al., 2005). For a bubbling fluidized bed with a jet, the ratio of normal second-order moment in the vertical direction to in the lateral direction is in the range of 0.5-2.5 (Patil et al., 2005). The bubblelike Reynolds normal stresses per unit bulk density used by Gidaspow et al. (2002) are computed from the simulated hydrodynamic velocities. The simulated bubblelike Reynolds normal stresses in the vertical direction is on average 4.5-6.0 times larger than that in the lateral direction in the bubbling fluidized bed (Jung et al., 2005). The predictions are in agreement with experimental second-order moments measured by Jung et al. (2005) and porosity measured by Patil et al. (2005).     

A bubbling fluidization model using kinetic theory of rough spheres

Collisional motion of inelastic rough spheres is analyzed on the basis of the kinetic theory for flow of dense, slightly inelastic, slightly rough sphere with the consideration of gas–solid interactions. The fluctuation kinetic energy of particles is introduced to characterize the random motion of particles as a measure of the translational and rotational velocities fluctuations. The kinetic energy transport equation is proposed with the consideration of the redistribution of particle kinetic energy between the rotational and translational modes and kinetic energy dissipation by collisions. The solid pressure and viscosity are obtained in terms of the particle roughness and restitution coefficient. The partition of the random‐motion kinetic energy of inelastic rough particles between rotational and translational modes is shown to be strongly affected by the particle restitution coefficient and roughness. Hydrodynamics of gas–solid bubbling fluidized beds are numerically simulated on the basis of the kinetic theory for flow of rough spheres. Computed profiles of particles are in agreement with the experimental measurements in a bubbling fluidized bed. The effect of roughness on the distribution of energy dissipation, kinetic energy, and viscosity of particles is analyzed. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Size segregation in gas-solid fluidized beds with continuous size distributions

Discrete-particle simulations of a gas-solid fluidized bed are used to investigate the species segregation (de-mixing) behavior of systems with continuous particle size distributions. Both Gaussian and lognormal distributions are investigated over a range of distribution widths, restitution and friction coefficients, and gas velocities. The results indicate that: (i) the average particle diameter decreases as the height within the bed increases, (ii) the level of segregation increases with an increase in the width of the particle size distribution, and (iii) segregation is attenuated as bubbling becomes more vigorous. Furthermore, the shape of the local size distribution (i.e., Gaussian or lognormal) is found to mimic that of the overall size distribution in most regions of the fluidized bed.     

Hydrodynamics of Geldart group A particles in gas-solid fluidized beds

Geldart group A particles were fluidized in a 10 cm i.d.×1.8 m high Plexiglas-made bed with ambient air to determine the hydrodynamic properties in a gas-solid fluidized bed. The effects of static bed heights, position of pressure measuring points, differential and absolute pressure fluctuations on the hydrodynamic behavior of a Geldart group A particles in a gas-solid fluidized bed were investigated. The particles used in this study were 80 micrometer FCC powders and 60 micrometer glass beads. The variance of pressure fluctuations was used to find the minimum bubbling velocity. The obtained minimum bubbling velocity was compared with the other methods available in the literature. This method was found to be much easier and had better data reproducibility than the classical visual method or sedimentation method. The variance of pressure fluctuations increased due to the increase of superficial gas velocity and static bed height. The obtained minimum bubbling velocity and pressure fluctuations were found to depend on the measuring position along the axial direction. The effect of measuring position was discussed. Cross-correlation of two pressure signals was used to find the delay time, then the bubble rising velocity.