锥形分布板射流流化床CFD模拟及参数分析

采用双流体模型结合颗粒动力学理论,对锥形分布板射流流化床内气固流动行为进行了三维的计算流体力学(CFD)模拟研究,系统分析了曳力模型、恢复系数和颗粒间摩擦力对射流流化床膨胀高度和气泡动力学行为的影响.结果表明,Syamlal-O'Brien和Gidaspow曳力模型低估了床内实际曳力,从而导致模拟膨胀高度低于实验值,而Modified Syamlal-O'Brien曳力模型能更好地预测床内实际曳力,计算结果与实验值吻合得较好;恢复系数对于射流流化床内气泡动力学行为有着重要影响,气泡的大小、上升速率和产生频率均随着恢复系数的增加而减小;颗粒间的摩擦力也是影响床层膨胀高度和床内气泡产生的重要因素.     

磷石膏颗粒湍动流化特性实验及模拟

针对磷石膏颗粒湍动流化体系曳力变化的问题,在实验的基础上,考虑非均匀结构对曳力的影响,引入修正因子φ修正Gidaspow曳力模型,对2D流化床进行了数值模拟研究。通过将Gidaspow模型在不同φ值下的模拟结果与实验结果进行对比,研究φ值的改变对模拟结果的影响规律及一定气速范围内磷石膏颗粒湍动流化体系曳力变化特性。结果表明,Gidaspow模型高估了实验体系曳力,对体系流化特性的预测效果较差;适当φ值的引入能明显提高Gidaspow模型对床层膨胀、压降及体系非均匀度的模拟精度。模拟结果反映出φ值越小,床层膨胀高度越低,床内颗粒浓度分布越不均匀,床层压降波动性越大。随着气速的升高(0.144~0.240 m/s),颗粒沿水平方向上的聚集程度加剧,φ值呈非线性减小(0.31~0.24)。流化体系的非均匀度随着气速增加而增大,颗粒浓度沿径向存在较大梯度,两侧边壁处附近出现环-核结构且流场分布对称性较差。     

浆态床内固含率轴向分布的数值模拟

采用DBS曳力模型计算气液相间作用,分别采用Gidaspow曳力模型、经Brucato修正的Gidaspow曳力模型和Schiller?Naumann曳力模型计算液固相间作用,忽略气固间的直接作用,对比了浆态床内不同颗粒粒径体系轴向固含率的模拟和实验结果. 结果表明,不同液固相间曳力模型对气含率的预测影响不大;在颗粒粒径较大(140 ?m)的体系中,较低表观气速下气液DBS与液固Schiller?Naumann曳力模型组合模拟的固含率随床高度增加而减小,与实验结果吻合,而其它曳力模型组合的模拟结果较差,轴向分布较均匀;在颗粒粒径较小(35 ?m)的体系中,几种曳力模型组合的模拟结果均与实验结果吻合较好,轴向分布较均匀.     

Geldart A类颗粒节涌床气固流动特性的实验及模拟

针对气固节涌床,在实验基础上,基于欧拉?欧拉双流体模型结合颗粒动力学理论,考虑Geldart A类颗粒聚团对气固间曳力的影响,采用修正后的Gidaspow曳力模型对气固节涌床进行数值模拟。结果表明,通过与实验结果及经验公式进行对比,修正的模型可准确合理地模拟流化床内节涌特性。表观气速0.09 m/s≤Ug≤0.39 m/s时,床层内部压力脉动标准偏差随表观气速增加而增加,流型由鼓泡转变为节涌直至节涌程度最大,床内气固流动主要受轴对称栓运动特性影响,床内压降、床层膨胀比、气栓平均上升速度、最大轴对称栓长度随表观气速增加而增加,最大轴对称栓产生位置随表观气速增加而降低;Ug>0.39 m/s后,床内压力脉动标准偏差随表观气速增加而降低,节涌程度降低至向湍动流态化流型转变,床内气固流动主要受壁面栓运动特性影响,增加表观气速,节涌床内压降变化幅度较小,气栓平均上升速度增加幅度加大,床层膨胀比及最大轴对称栓长度降低,最大轴对称栓产生的位置略有升高。     

稠密同轴气固射流的CFD模拟

同轴气固射流在能源领域具有广泛地应用,但大多数研究集中在颗粒浓度较低的工况。为了研究稠密同轴气固射流的流动特性,采用了稠密离散相模型(DDPM)耦合离散元模型(DEM)的方法对该体系进行计算流体动力学(CFD)模拟,该方法同时考虑了孔隙率对气固曳力的影响和颗粒间的碰撞作用。由于射流过程中,气体对颗粒的作用占主导,分别考虑了不同环形气体速度和气固曳力模型对气固流动的影响。模拟结果表明,该模型能合理地模拟在不同气速下稠密气固两相射流的颗粒弥散特性,与实验现象定性一致。在较高气速下,引入湍流模型对预测结果有显著影响,模拟得到的颗粒弥散程度较大。不同气固曳力模型对颗粒弥散的预测有明显影响,WenYu曳力模型下颗粒弥散程度较大,Gidaspow模型次之,SyamlalO’Brien模型给出的颗粒弥散程度较小。     

基于EMMS介尺度模型的双分散鼓泡流化床的模拟

双分散气固鼓泡流化床中颗粒通常具有不同粒径或密度,导致产生颗粒偏析等现象,影响传递和反应行为。颗粒分离和混合与气泡运动密不可分,其中相间曳力起关键作用。最近Ahmad等提出了一种基于气泡结构的双分散介尺度曳力模型,能成功预测双分散鼓泡流化床的床层膨胀系数。本研究耦合该曳力模型与连续介质方法,模拟了两种不同的双分散鼓泡流化床,通过分析不同流化状态下的气泡运动、颗粒浓度比的轴向分布等参数,进一步检验模型的适用性。研究表明,当双分散颗粒处于完全流化状态时,耦合双分散介尺度曳力模型可合理预测不同颗粒的分离现象;而其处于过渡流化状态时,新曳力模型和传统模型均无法获得合理结果,此时调节固固曳力可改进模拟结果。     

基于EMMS介尺度模型的双分散鼓泡流化床的模拟

双分散气固鼓泡流化床中颗粒通常具有不同粒径或密度,导致产生颗粒偏析等现象,影响传递和反应行为。颗粒分离和混合与气泡运动密不可分,其中相间曳力起关键作用。最近Ahmad等提出了一种基于气泡结构的双分散介尺度曳力模型,能成功预测双分散鼓泡流化床的床层膨胀系数。本研究耦合该曳力模型与连续介质方法,模拟了两种不同的双分散鼓泡流化床,通过分析不同流化状态下的气泡运动、颗粒浓度比的轴向分布等参数,进一步检验模型的适用性。研究表明,当双分散颗粒处于完全流化状态时,耦合双分散介尺度曳力模型可合理预测不同颗粒的分离现象;而其处于过渡流化状态时,新曳力模型和传统模型均无法获得合理结果,此时调节固固曳力可改进模拟结果。     

基于EMMS模型的气固鼓泡床的模拟及气泡特性的分析

基于EMMS曳力模型,采用双流体的方法对气固鼓泡床内的气固流动特性进行模拟,建立基于图像处理气泡特性的分析方法,重点研究了不同表观气速下气泡在床层内分布特性,包括气泡平均当量直径、气泡速度和气泡球形度的轴向分布,以及气泡的生命周期。研究结果表明,小气泡多集中在床层底部和壁面区域,而大气泡多集中在床层中间区域。随着表观气速的增加,床层高度不断增加,气泡的球形度降低,气泡的大小、出现频率、上升速度以及生命周期均增加;然而,当表观气速增大到一定程度,继续增加气速对气泡的上升速度影响不大。     

基于EMMS模型的气固鼓泡床的模拟及气泡特性的分析

基于EMMS曳力模型,采用双流体的方法对气固鼓泡床内的气固流动特性进行模拟,建立基于图像处理气泡特性的分析方法,重点研究了不同表观气速下气泡在床层内分布特性,包括气泡平均当量直径、气泡速度和气泡球形度的轴向分布,以及气泡的生命周期。研究结果表明,小气泡多集中在床层底部和壁面区域,而大气泡多集中在床层中间区域。随着表观气速的增加,床层高度不断增加,气泡的球形度降低,气泡的大小、出现频率、上升速度以及生命周期均增加;然而,当表观气速增大到一定程度,继续增加气速对气泡的上升速度影响不大。     

基于气泡EMMS模型的鼓泡床两相流动CFD模拟

采用两步法将气泡EMMS模型拓展到计算网格层次,用于模拟鼓泡床两相流动,改进的EMMS曳力模型不仅与空隙率有关,且与网格内的速度直接关联.求解是在操作条件下利用宏观稳定性约束条件寻优得到乳化相内空隙率关系式,再用计算网格层次上的微观守恒方程结合乳化相空隙率关系式封闭求解其余结构参数及对应的曳力系数,最后借助用户自定义函数将气泡EMMS曳力嵌入双流体模型中,对含Geldart A类颗粒的鼓泡床内气固两相流动行为进行模拟与验证.结果表明,该模型能较好地模拟鼓泡床内的非均匀流动结构,准确捕捉到上稀下浓的轴向颗粒浓度分布特征,轴向颗粒浓度分布结果与实验值较接近,平均相对偏差为6.4%.床层径向上颗粒分布均呈明显的环-核结构.在不同床层高度(0.6,0.8和1.1 m)处的径向颗粒浓度分布与实验值基本吻合,仅在底部(0.4 m)略低,相对偏差小于14%.而Gidaspow曳力模型在4个床层高度与实验数据差距甚远,相对偏差可达91%.     

不同曳力模型及颗粒碰撞恢复系数对短接触旋流反应器内气固流场的影响

为考察曳力模型和颗粒碰撞恢复系数对短接触旋流反应器内流动特性的影响,基于双流体模型, 结合颗粒动力学理论,对反应器内气固两相流场进行模拟研究。分别采用Gidaspow、Wen & Yu和Syamlal-O’Brien 3种曳力模型, 考察颗粒速度特性以及固含率径向分布。对比分析不同曳力模型的计算结果表明,Syamlal-O’Brien模型计算结果与实验结果误差较大,Wen & Yu模型在反应器边壁附近区域的计算结果误差较大,Gidaspow模型计算结果与实验结果最为吻合。此外,颗粒碰撞恢复系数较小时,所得计算值小于实验测量值,当恢复系数为0.95时颗粒扩散效果最好,计算结果与实验数据吻合度最高。     

Hydrodynamic modelling of dense gas-fluidised beds: comparison and validation of 3D discrete particle and continuum models

A critical comparison of a hard-sphere discrete particle model, a two-fluid model with kinetic theory closure equations and experiments performed in a pseudo-two-dimensional gas-fluidised bed is made. Bubble patterns, time-averaged particle distributions and bed expansion dynamics measured with a nonintrusive digital image analysis technique are compared to simulation results obtained at three different fluidisation velocities. For both CFD models, the simulated flow fields and granular temperature profiles are compared. The effects of grid refinement, particle-wall interaction, long-term particle contacts, particle rotation and gas-particle drag are studied. The mechanical energy balance for the suspended particles is introduced, and the energy household for both CFD models is compared. The most critical comparison between experiments and model results is given by analysis of the bed expansion dynamics. Though both models predict the right fluidisation regime and trends in bubble sizes and bed expansion, the predicted bed expansion dynamics differ significantly from the experimental results. Alternative gas-particle drag models result in significantly different bed dynamics, but the gap between model and experimental results cannot be closed. In comparison with the experimental results, the discrete particle model gives superior resemblance. The main difference between both CFD models is caused by the neglect of particle rotation in the kinetic theory closure equations embedded in the two-fluid model. Energy balance analysis demonstrates that over 80% of the total energy is dissipated by sliding friction. Introduction of an effective restitution coefficient that incorporates the additional dissipation due to frictional interactions significantly improves the agreement between both models.     

Experimental and numerical research for fluidization behaviors in a gas–solid acoustic fluidized bed

The effects of sound assistance on fluidization behaviors were systematically investigated in a gas–solid acoustic fluidized bed. A model modified from Syamlal–O'Brien drag model was established. The original solid momentum equation was developed and an acoustic model was also proposed. The radial particle volume fraction, axial root‐mean‐square of bed pressure drop, granular temperature, and particle velocity in gas–solid acoustic fluidized bed were simulated using computational fluid dynamics (CFD) code Fluent 6.2. The results showed that radial particle volume fraction increased using modified drag model compared with that using the original one. Radial particle volume fraction was revealed as a parabolic concentration profile. Axial particle volume fraction decreased with the increasing bed height. The granular temperature increased with increasing sound pressure level. It showed that simulation values using CFD code Fluent 6.2 were in agreement with the experimental data. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

CFD modeling and validation of the turbulent fluidized bed of FCC particles

An experimental and computational study is presented on the hydrodynamic characteristics of FCC particles in a turbulent fluidized bed. Based on the Eulerian/Eulerian model, a computational fluid dynamics (CFD) model incorporating a modified gas‐solid drag model has been presented, and the model parameters are examined by using a commercial CFD software package (FLUENT 6.2.16). Relative to other drag models, the modified one gives a reasonable hydrodynamic prediction in comparison with experimental data. The hydrodynamics show more sensitive to the coefficient of restitution than to the flow models and kinetics theories. Experimental and numerical results indicate that there exist two different coexisting regions in the turbulent fluidized bed: a bottom dense, bubbling region and a dilute, dispersed flow region. At low‐gas velocity, solid‐volume fractions show high near the wall region, and low in the center of the bed. Increasing gas velocity aggravates the turbulent disorder in the turbulent fluidized bed, resulting in an irregularity of the radial particle concentration profile. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

液固散式流态化CFD模拟中曳力模型的影响

基于无黏性双流体简化模型在商业化软件平台上,通过增加用户自定义子程序考察了Gidaspow、Syamlal-O’Brien、Di Felice、Gibilaro、Dallavalle和BVK曳力模型对液固散式流态化CFD模拟结果的影响行为,探讨了相应的影响机制。经与文献中不同颗粒Reynolds数的代表性实验数据对比后发现：BVK和Dallavalle曳力模型对床层膨胀高度和整体固含率的预测精度较高;BVK、Syamlal-O’Brien以及Dallavalle曳力模型给出的床内固含率径向分布较为准确;BVK曳力模型较为准确地再现了颗粒轴向速度的径向分布特征。BVK曳力模型的影响机制与液固散式流态化中颗粒动力学特性相符合,在所考察范围内其预测性能最优;Dallavalle曳力模型在其余5个传统模型中预测性能较优且形式简洁在程序中易于实现。     

三相流化床间歇反应器中固体轴向浓度分布

研究了在直径 4.2cm、高 130cm的垂直三相流化床中固体颗粒非粘附与粘附两种不同体系中固含率的分布特性 .基于混合长理论与实验数据 ,提出了将“沉降 -扩散”模型同时应用于描述非粘附与粘附体系 (进行气速的修正 )中固含率的轴向分布 ,模型预测与实验数据很好地吻合     

用改进的单元胞模型数值模拟液体穿过球形颗粒群的流动

The cell model developed since 1950s is a useful tool for exploring the behavior of particle assemblages,but it demands further careful development of the outer cell boundary conditions so that interaction in a particle swarm is better repressented.In this paper,the cell model and its development were reviewed,and the modifications of outer cell boundary conditions were suggested.Athe cell outer boundary,the restriction of uniform liquid flow was removed in our simulation conducted in the reference frame fixed with the particle.Zero shear stress condition was used to evaluate the outer boundary value of the stream function.Boundary vorticity was allowed to evolve to values compatible to existing stream function at the free shear outer boundary.The fore-aft symmetry of vorticity distribution at the outer boundary is thought critical to ensure the continuity of inflow and outflow between touching neighbor cells,and is also tested in the modified cell model.Numerical simulation in terms of stream function and vorticity based on the modified cell models was carried out to shed light on the interaction between liquid and particles.Lower predicted drag coefficient by the modified cell models was interpreted with the feature of flow structure.The drag coefficient from the simulation was also compared with correlations of drag coefficient reported in literature.It is found that the modified cell model with the uniformity of external flow relaxed and the fore-aft symmetry of boundary vorticity enforced was the most satisfactory of the overall performance of prediction.     

湍动流化床过渡段固含率分布特征的实验及数值模拟

在湍动流化床中,过渡段对于包括甲醇制烯烃在内的气固催化快反应有着重要的作用。采用PV6D反射型光纤探针对内径95 mm的湍动流化床内过渡段的固含率分布和脉动参数进行了测量,分别考察了表观气速和静床高的影响,并采用修正的基于颗粒动力学的三段曳力双流体模型进行模拟。实验表明,湍动流化床过渡段中固含率的轴向分布呈现S型和指数型两种类型,固含率轴向与径向分布都在过渡段内出现最大梯度,表明过渡段中固体浓度分布比稀相段和密相段更不均匀。表观气速和静床高的变化将导致S型和指数型分布的相互转变,并且对过渡段底部与壁面附近的固体高浓度区影响最为显著。局部固含率脉动概率密度分布表明,在静床高较小时,随着气速的增大,床层下部气含率最大值位置将从中心区移动至环隙区,呈现气含率的双峰型分布。本文提出的修正三段曳力模型考虑了颗粒团聚的影响,对过渡段中分布板影响区之外的固含率分布均能较好地模拟。     

Comparison and evaluation of interphase momentum exchange models for simulation of the solids volume fraction in tapered fluidised beds

An Eulerian computational fluid dynamics (CFD) model with granular flow extension was used to simulate a gas–solid fluidised bed in a tapered reactor. Various drag coefficient models were evaluated, which are used to calculate the drag force, describing the momentum transfer between the gas and solid phases. Comparison and evaluation between time-averaged solids volume fractions obtained from experiments and from simulations with several drag coefficient models were made. The predicted results obtained by the different drag models were verified using experimental data of Depypere et al. (2009). Initial results using a 2-phase Eulerian model showed poor agreement with experimental results. However, extending the Eulerian model to include 3 solid phases—with different mean particle diameter per phase in order to account for the particle size distribution of the fluidised solid material—yielded good agreement with experimental results. Furthermore, quantitative analyses showed that the modified Gidaspow drag model gave the best agreement between CFD simulations and experimental data.