大颗粒气固流化床内两相流动的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气泡在床层底部的扰动及表面的抛洒作用．由此得出，气泡对流化床内颗粒及团聚物的横向运动有着重要的影响．

流化系统中的流动结构不仅呈现局部非均匀性，而且呈现整体非均匀性．局部非均匀性表现为稀

相和密相在同一点交替出现；整体非均匀性表现为系统内部不同空间位置可以出现稀相或密相两种完全不同的结构[1]．图5显示了流化床内颗粒体积分率瞬时分布，其中颗粒密度为906 kg／m3，h为床层轴向位置，H为固定床高度， 为床层径向位置，R为流化床的半径．图中(a)、(c)在初始流化床高度( ／H=1)的中心位置处( ／R：O)，(b)、(d)则在靠近壁面处( ／R：0．75)．由图5(a)、(b)曲线可知：在直型布风板流化床里，颗粒体积分率在靠近壁面处几乎都在0．25以上波动，在床层中心位置多个时间段内小于0．1，甚至有几段几乎为0，这些说明颗粒和颗粒微团主要集中在靠近床层壁面处，气泡则位于中心位置处．相对直型布风板流化床，颗粒在凹型布风板流化床床层底部有向中心运动的能力，但通过分析固定床高处颗粒体积分率瞬时分布却得知：和直型布风板流化床内颗粒体积分率分布一样，颗粒微团靠近壁面，气泡位于中心处，如图5(c)、(d)所示．

3．2 颗粒密度的影响

图6显示了在初始流化阶段，不同密度颗粒在凹型布风板流化床内的体积分率分布和流线图．由图可见，随着颗粒密度的增大，其体积分率分布及运动趋势都呈现出了很大的差异．在颗粒密度P ：906 kg／m3时，颗粒体积分率在流化床内主要分布在中心及壁面处．通过云图可知，出现这种现象的原因是由于流化床内存在两个大气泡。当颗粒密度增大到2 700 kg／m3时，大气泡也跟着增加到了3个；再当密度增加到3 600 kg／m3时，不但大气泡消失了，而且颗粒主要聚集在流化床中心位置．根据流线的切线方向即为颗粒速度方向的性质可知：床层上部的颗粒将从中心位置处下落，并与床层底部上升中的颗粒进行对冲．这一过程将对流化床带来诸多的不利影响，即增加固体颗粒的机械磨损及损坏床内零件等．

图7显示了凹型布风板流化床内颗粒径向平均速度( 。)分布．其中床内高h=0．2 m， 为床层径

向位置，R为流化床的半径．随着密度逐渐增大，颗粒速度在中心位置附近减小得并不大，而靠近壁面处则有很大的降低．之所以出现这一现象，是由于气泡与重力综合作用的结果．当操作气速一定时，颗粒密度增大，使得床层高度降低，颗粒体积分率，尤其是壁面处的体积分率增大，进而导致颗粒聚集倾向明显增强，团聚物增大．而团聚物受气体的影响主要发生在表面附近，所以当颗粒重力成倍地增长时，其在边壁处的速度便明显下降．这时大量气体将更加集中于床层中心，在维持气体通量沿截面平衡的同时，将以更快的速度携带颗粒运动．结果，颗粒在床层中心处的速度并没有很大变化．

4 结论

通过对比大颗粒在不同布风装置、曳力模型及密度情况下的数值模拟，可以得出以下结论：

a．相比于直型布风板流化床，凹型布风板流化床内气固混合将更加充分，主要表现为气泡产生快，颗粒横向运动能力强．b．颗粒密度增大时，凹型布风板流化床在初始流化阶段将呈现出不同的结构；充分流化后，在靠近边壁处的颗粒速度比中心位置处下降的快．C．采用Arastoopour、Gidaspow 和Syamlal-0’Brien 3种曳力模型模拟轴向空隙率分布规律，可以发现：在床层中心稀相区，Arastoopour模型差异较大；在靠近壁面处，Syamla1．0’Brien模型差异较大．另外，Arastoopour和Syamlal一0’Brien模型预测流化后的压降也比Gidaspow模型低；d．与床层膨胀比实验关联式相比，Arastoo—pour、Gidaspow和Syamla1．0’Brien 3种模型预测的值比实验关联式值要大一些；e．对3种曳力模型进行各种比较与研究后发现，Arastoopour、Gidaspow和Syamlal一0’Brien中Gidaspow模型相对适合于大颗粒气固流化床的数值模拟．

参考文献：

[1] 李静海，欧阳洁，高土秋，等．颗粒流体复杂系统的多

尺度模拟[M]．北京：科学出版社，2005．

[2] GELDART D．Types of gas fluidization fJ]．Powder

Techno1．1973(7)：285—297．

[3] 金涌，祝京旭，汪展文，等．流态化工程原理[M]．北

京：清华大学出版社，2001．