化工学报 ›› 2018, Vol. 69 ›› Issue (8): 3443-3451.DOI: 10.11949/j.issn.0438-1157.20171666

• 流体力学与传递现象 • 上一篇    下一篇

快速流化床颗粒团絮特征的MP-PIC数值模拟

孙子文1, 陈岱琳2, 钟文琪1,3, Aibing Yu3,4   

  1. 1 东南大学能源热转化及其过程测控教育部重点实验室, 江苏 南京 210096;
    2 福建省发展和改革委员会工业发展处, 福建 福州 350003;
    3 东南大学-蒙纳士大学联合研究院颗粒系统模拟与仿真中心, 江苏 苏州 215123;
    4 澳大利亚蒙纳士大学化学工程系, 克莱顿 Vic 3800
  • 收稿日期:2017-12-19 修回日期:2018-03-14 出版日期:2018-08-05 发布日期:2018-08-05
  • 通讯作者: 钟文琪
  • 基金资助:

    国家重点研发计划项目(2016YFB0600802);国家自然科学基金重大项目(51390492)。

MP-PIC simulation of particle clusters in fast fluidized bed risers

SUN Ziwen1, CHEN Dailin2, ZHONG Wenqi1,3, YU Ai-Bing3,4   

  1. 1 Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, Southeast University, Nanjing 210096, Jiangsu, China;
    2 Department of Industrial Development, Fujian Provincial Development and Reform Commission, Fuzhou 350003, Fujian, China;
    3 Centre for Simulation and Modelling of Particulate Systems, Southeast University-Monash University Joint Research Institute, Suzhou 215123, Jiangsu, China;
    4 Department of Chemical Engineering, Monash University, Clayton, Vic 3800, Australia
  • Received:2017-12-19 Revised:2018-03-14 Online:2018-08-05 Published:2018-08-05
  • Supported by:

    supported by the National Key Research and Development Program of China (2016YFB0600802) and the National Natural Science Foundation of China (51390492).

摘要:

为了研究快速流化床颗粒的团絮特征,建立了基于多相质点网格法的快速流化床气固多相流三维数理模型,气相场采用大涡湍流模型,通过求解牛顿运动方程得到颗粒相运动信息,气固间相互作用力采用Gidaspow曳力模型,固体间作用力通过计算颗粒应力梯度得到。基于该模型,对三维快速流化床上升管(H=3 m、d=0.1 m)气固流动开展了数值模拟,并与实验进行了校正,研究了在气速工况Ug=5.28 m·s-1下的颗粒(ρp=2650 kg·m-3dp=250 mm)团絮性质,实现了对上升管内颗粒团絮的基本类型(条形团絮、马鞍形团絮、U形团絮)的成功预测,并揭示了不同类型团絮在上升管内形成、发展、聚并直至破碎的演化规律。结果表明,上升管径向颗粒团絮的平均颗粒浓度分布呈现中间低两边高的环核结构,颗粒团絮速度的分布与其相反;随着轴向高度的增加,颗粒团絮的颗粒浓度逐渐降低而速度逐渐增加,但达到一定高度后变化减缓。

关键词: 欧拉-朗格朗日, 多相质点网格法, 数值模拟, 快速流化床, 颗粒团絮, 多相流, 介尺度

Abstract:

A three-dimensional Eulerian-Lagrangian model, based on MP-PIC method, was developed to simulate particle clusters in fast fluidized bed risers. In the model, gas phase was simulated by large eddy simulation (LES), real particles were displaced by numerical-particles and particles was described by Newton movement equations, drag between gas and particles was calculated by Gidaspow model, and normal stress of particles was calculated and interpolated into Eulerian grids to describe particle contraction. The model was used to simulate gas solid movement in a 3D fast fluidized bed riser with dimension of height H=3 m and diameter d=0.1 m and adjusted according to experimental data. Cluster characteristics of particles (ρp=2650 kg·m-3 and dp=250 mm) was studied at operational condition of Ug=5.28 m·s-1gas velocity. The results indicate that the MP-PIC simulation can recognize different shapes of clusters, such as stripe-shaped, saddle-shaped and U-shaped, as well as disclose cluster development including formation, expansion, coalescence and breakup. Mean particle concentration distribution of clusters in the riser cross section shows a core-annular structure with low middle and high edge concentrations, which is opposite to mean cluster velocity distribution. The mean particle concentration of clusters decrease while the mean cluster velocity increase with increase of riser height, but the change of both concentration and velocity slows down above certain riser height.

Key words: Eulerian-Lagrangian, MP-PIC, numerical simulation, fast fluidized bed, particle clusters, multiphase flow, mesoscale

中图分类号: