大颗粒气固流化床内两相流动的CFD模拟

大颗粒气固流化床内两相流动的CFD模拟

摘要：采用欧拉双流体模型和颗粒动力学方法，数值模拟了大颗粒流化床在不同密度、布风装置及曳力模型情况下的气固两相流动，考察了大颗粒流化床流化和流动特点，颗粒体积分率分布，床层压力瞬时变化，床层碰撞比，以及颗粒速度径向和空隙率轴向分布规律．研究结果表明，与直型布风板流化床比较，凹型布风板流化床内的气泡产生快，颗粒横向运动能力强；随着颗粒密度的增大，其在凹型布风板流化床边壁处的速度比中心位置处减小的快；比较3种曳力模型，发现其模拟的轴向空隙率分布和床层压力存在较大差异，且与床层膨胀比实验关联式相比，3种模型预测的值比实验关联式要大一些．通过研究，3个曳力模型中Gidaspow模型相对适用于大颗粒气固流化床的数值模拟．

关键词：流化床；欧拉双流体模型；并行计算；大颗粒

近年来，随着流态化技术的发展，大颗粒流化床在煤粉流态化燃烧和水泥熟料流态化煅烧等领域的应用也越来越广泛．由于流化床内两相流动情况复杂，使得人们对气固两相间的作用、固相应力本构方程的建立、两相湍流的认识以及多种因素的相对控制和协调的理解等变得很困难[】]．实际上大多数流化床反应器都是根据经验设计的，大颗粒流化床的设计更是如此．文献[2]在研究颗粒的粒度及颗粒的表观密度等对流化特性影响后，将颗粒分成了A(30～100 tma)、B(100～600 tLm)、C(一般情况下粒度小于20 tLm)、D(600 Fm以上)4类_3]．依据此分类，粒度在600肿以上的颗粒称为过粗颗粒．然而由于颗粒的

表观密度与气体密度之差不同，本文所用颗粒直径为855 可能为B类(鼓泡颗粒)，也有可能为D类(喷动用颗粒)．其中，D类颗粒流化时极易产生大气泡或节涌，使实验难以操作，然而数值模拟可以克服这一困难，而且D类颗粒粒度在1．5 rain以下时，是完全

可以流化的[3]．文献[4]用粒径为3 mm的颗粒进行了模拟与实验，研究了气体进口速度和温度对床内含湿量、颗粒温度等的影响，得出模拟与实验的结果大体是一致的．文献[5]研究了表观气速、床内有无管道及布风方式对大颗粒流动的影响．模拟和试验的结果都表明，布风方式对颗粒体积分率及速度径向分布有着很大的影响，而且不论有无管道，某些布风方式都有助于气固形成环核流动结构．文献[6]通过改变颗粒粒径(从o．25 mm到1 mm)、密度、进口气速等参数后进行了模拟，结果表明：颗粒的粒径和进口气速对颗粒滑移速度的影响较大；合适的进口气速对减少能

耗起着很重要的作用．本文借助CFD软件FLUENT对大颗粒气固流化床进行了模拟计算．对比并分析了不同密度颗粒、曳力模型及布风装置对流化床流动特性的影响．有些曳力模型采用皿F(用户自定义函数)实现．通过这些研究，从数值计算的角度揭示出了一些大颗粒的流化及流动特性．

1 控制方程及曳力系数模型

1．1 流体控制方程

由于气固间没有质量交换，且升力、附加质量力等对流化床的影响很小，故气固两相流动所遵循的连续方程和动量方程可以简化成如下形

式：

动量方程

1．2 曳力系数模型

颗粒在流场中受到的作用力包括曳力、重力、浮力和其他作用力(如Basset力、Magnus力和Saff．man力等)．若忽略其他力的作用，则可认为气固间作用主要为曳力作[1]．Syamla1．0BriencAmstoopour~。]和GidaspowC。]等人先后对气固曳力作了大量的研究，并给出了反应相间作用强弱程度的曳力系数，如下所示：

2 流化床几何结构和模拟参数设定

图1(a)、(b)分别给出了两种不同布风板流化床的几何结构．两个流化床的床高和床宽分别为

1．2 m 和0．186 m(其中，凹型布风板的圆心距离床底0．116 m)．由于计算需要在一组离散的网格点上进行，所以模拟采用了G 舳IT

2．2．30(Fluent 公司，美国)进行网格生成．与表1等人相比，网格

不但比其他的要小得多，而且数目也多．为了提高数值计算的速度，本次模拟使用了基于Linux操作系统的联想并行机进行求解计算．求解器采用Fluent6．3．26(Fluent公司，美国)，并选取一￡湍流模型及速度一压力耦合的“SIMPLE”算法．初始时固定床床层的高度日为0．372 m，固相体积分率为0．622，最大时达到0．65．流化床床内无任何构件，气体的密度为1．225 kg／m3，黏度为1．785×10～，颗粒直径为855 m．流化床的边界条件为：两壁面均假定成无滑移，下部为气体匀速进口，速度大小为1 m／s；上部为压力出口，表压为0．依据颗粒动力论，将固相压力及黏度分别选用Lun。et．al和Syam．1a1．0’Brien模型，并把颗粒的弹性恢复系数假定为0．90(1代表完全弹性碰撞，0则表示完全非弹性碰撞)．

¨

3 结果与讨论

3．1 布风方式的影响

图2(见下页)描述了颗粒在直型布风板流化床内不同时刻的体积分率分布．可以看出，流化床内部可以分为底部床层和上部自由空间．0 S 时颗粒静止地堆积在床底，当匀速气体从进口流进时，床内逐渐形成了典型的中心区颗粒体积分率小，边壁区颗粒体积分率大的环一核流动结构[1引，图2中0．3 s～0．6 S时间内显示了其部分变化过程．从图2中0．5 s～1．117 S可看出，气泡主要在床中心与壁面之间的区域内产生与上升．随着气泡的上升其形状和大小都在发生变化．在

1．117 S气泡发生破裂，颗粒被抛人自由空间．

图3显示了颗粒在不同时刻速度矢量场．其中t表示时间，颗粒密度为2 000 kg／m3，最右端t=

2．5 S为床层表面情况，其他的均为床层底部．通过图形可以看出，随着时间变化，直型布风板流化床中颗粒及团聚物逐渐有了横向运动．除了t=0．2 S主要是由于重力与床壁作用外，其余是由气泡引起的．t=0．5 S气泡在床层底部的产生；t=2．5 S气泡在床层底部的扰动及表面的抛洒作用．由此得出，气泡对流化床内颗粒及团聚物的横

向运动有着重要的影响．