哈尔滨工业大学循环流化床课件

第二章化床气固两相流流体动特性2.1化的义

：固体颗粒在气体的作用下，呈现出流体的特性，即流态化，其实质是气固两相流动。当颗粒处于流态化时，作用在固体颗粒上的重力与气流的拽力相平衡，颗粒处于一种悬浮状态，具有以下特点：1.任一高度的静压近似等于此高度上单位床截面内固体颗粒的重量。

)(1

)隙率）2.床层倾斜，床表面总是保持水平。3.床内固体颗粒可以像液体一样，从底部或侧面的孔中喷出流口）4.可以像液体一样从一个容器流入另一个容器，并保持水平。5.窄筛分颗粒流化后会非常均匀，宽筛分颗粒会分层流化。6.床内颗粒混合良好，整个床层的温度基本均匀。2.2定的降：流体通颗层阻固定床压降的研究是在渗流基础上发展起来的。研究表明，流体以层流的形式通过散料床层时，流体的空塔截面速度与床层压降△正比。当1<Re<20

;d

当Re>1000时

d

u

式中H为床层高度动力黏度pa.s为气体空截面速度m/s

为气体密度，kg/m。Ergun前人研究结果的基础上通过试验得出包含层流和湍流床层压降综合表达式：

(11.75()

式中为球形度，定义为与实际颗粒体积相等的球颗粒的表面积与其实际表面积之比为床层空隙率。该式第一项为黏性项，当流速较低时，它占主导作用；第二项为惯性项，当流速较高，流动为湍流时，该项起主要作用。2.3界化度临界流化是流化床操作的最低速度，是描述流化床的基本参数之一。确定临界流化速u的方法主要有理论计算和实验测定两种。临界流化速度u是当床层压降等于床层颗粒重量时所对应的流体速度，可由Ergun程导出AH)(

pg

)P

(

)(

)g150H

)

1umf)3

s

p

s

p150

(1

mf3

2

mf)

1.75

1mf

gmf

)(1pg

mf

gmf

p

mf

p两边同乘以

dp(1

)

得到150

)ugmf

11.75(

dd()p2g

u

d

)2r

C42.8571

)mf;Css

r12gsp；0.316kgm3r12gsp；0.316kgm3r

3()ggpg

，叫阿基米德数动力黏度，。计u公式中C、C系数值作

者

C

1

C

2

作

者

C

1

C

2Wen√Bourgeois(1968)Saxena(1977)

33.725.7

BabuThonglimp(1984)Grace(1982)√

19.927.2

对于宽筛分物料（0－10mm内采用下面公式计算临界流化速度

mf沸腾层温度下运动粘度，m

A2。yg

0.528r

y床层空隙率，若在体积为V粉体中，气体与颗粒所占的体积分别为V和。

VVsVVVsVbulk

sbulk

；

bulk

p细颗粒的空隙变化范围为0.26-0.48，通常为0.45左右。空隙大小与颗粒的形状、颗粒粒径分布、填充方式、流化条件有关，对于燃煤锅炉，一般

/3,kg/3m

－＝0.5为球形度，表示颗粒接近球体的程度，例如，球形颗粒的球形度为。一些典型物料的球形度物料

沙平均值

石灰石（压碎）煤（磨碎）

园砂

粉碎的砂石

0.750.650.73

例题1计算300砂粒的临界流化风速。已知砂粒密度床温825℃体密度体动力黏g

3，p4.49s/

p沸腾层温度下运动粘度，m/s。p沸腾层温度下运动粘度，m/s。Ar

gpg

0.00033(4.49)2

=103.7u

mf

(A)drg

0.3160.0003

[(27.2103.7)27.2]例题度－10mm为d，度kg3烟温为900℃已知烟气密度烟气动p力粘度×求临界流化速度。A0.528mfr

2y0.314

m

2

Ar

3)ggppg

0.3143)2

＝11537.70.0882mfr

0.088211537.7

148

m/s如果按窄筛分物料公式计算得到：u

d

(C2)12r1

g0.314

[(27.2211537.7)

温度对不同粒径颗粒的临界流化速度的影响粒径

）

27空气中

(cm/s)

800℃烟气中

(cm/s)4080100

0.05

ofpiofpi20040060080010002000300040005000

3.96

1.295.16350300

-

800C250

-

27

o

Csmc

200150

m

100500010002000300040005000d(p不同温度下临界流化速度随粒径的变化关系2.4均粒度d国内流化床床料平均粒径计算方法采用量平均径d,

0.600.62p

s

i

i

sx为两筛之间物料的重量份额。

d

i

i

i

，

dii

分别为相邻两筛的筛

imfbftfimfbftf孔直径。国外采用调和平均（sauter或称比表面积体积平均径。

ii

，

di

dii2在工业中，粒度分布有时采表示，表示有50重量的样品低于这个粒度。例如，循环流化床锅炉煤入炉粒度要求<1.25mm。2.5界化度

mf

实测确定临界流化速度u

mf

的最好方法是通过实验测定。降低流速使床层自流化床缓慢地复原至固定床，同时记下相应的气体流速u和床层压eq\o\ac(△,降)P，在双对数坐标纸上标绘得到如图－3曲线。通过固定床数据区和流化床数据区的点各自画（撇开中间区数据条直线的交点即是临界流化点，对应的横坐标的速度即是临界流化速度。图u为起始流态化速度，此时床层中有部分颗粒进入流化状态。图中为完全流态化速度，此时床层中所有颗粒进入流化状态。对于粒径分布较窄的床层，bf、mfutf三者非常接近，有时很难区分。2.6粒分1973年Geldart将固体颗粒分为四大类，即A、B、C、四类。颗粒的属于哪一类主要取决于颗粒的尺寸和密度，同时也取决于流化介质的性质。如图所示。

C类粒子细，一般小于

kg/子间作用力与重力相近，难以流化，如果要流化该类粒子，则需特殊技术，否则常会造成沟流。A类粒：一般20－100

（kg/m类粒子能很好流化，但表观速度在超出临界流化速度后，气泡不立即出现，床层会发生明显膨胀，进一步增加气速，膨胀的床层会出现气泡，存在最小鼓泡速度,u。很多循

环流化床分离器分离下的颗粒属于A类子。B类粒：粒径通常在100－500

m（一说为650

kg们能够流化得很好气泡在大于临界流化速u后立即出现大部分CFB锅炉采用这类粒子。D类粒：所有颗粒中最粗的，通常达到1mm更大。虽然它们也会鼓泡，但固体颗粒的混合相对较差，更容易产生喷射流。需要较高的气速去流化。下表给出GraceJ.1984给出的不同类型的颗粒特性比较。

mbmbA类粒子鼓泡速度u可由和Deldart(1980)提出的公式计算：

的最小

0.716Fd(

0.347)

gF是颗粒小于

的质量份额。2.7式态和泡态和郭慕孙首先采用弗鲁德准则数来区分这两种流化形式。散式流态化

Frmf

umfd

鼓泡或聚式流态化：

Fr

pu2dg

p也有建议采用4个无因次群来表征流态化形式散式流态化：

(Fr)(mfmf

)(

HmfD

)100聚式流态化：

(Fr

)

u

bbpbbp2.8泡模简最著名的鼓泡床模型是将床层处理为气泡相和乳相。许多流化床模型都假定多于临界流化的气体量都以气泡的形式通过床层，即

VA

mf经过大量的研究后发现：(1)的实测值通常小于u-u的计算值，在气栓流（节涌）条件下，与实测值比较接近，但在布风板上方附近和高速条件下，的实测值与-差别特别大。(2)对上式还不能做出经验或半经验的修正。(3)存在无法实测的两相之间气体流动分量。对于气泡相，气体能够通过气泡边界；对乳化相，气体从颗粒间隙流入气泡。目前尚不能定量描述这一气体的流动分布。(4)对上式可以谨慎地用于近似估算。有许多作者将平均气泡直径表示为高度的函数，其中最好的Mori与Wen和Darton提出的模。在与模型中，首先计算最大气泡直径

(u)]

气泡沿床高的分布为d)de

该公式主要应用于0.5≤，60≤450，u-u≤，当量床直径D≤130cm的操作范围。其主要特点是考虑了床径对气泡大小的影响。Darton出了一个气泡沿床高几何合并的模型，得出：)A)g

tD浮DtD浮D式中，A/N是每个小孔对应的布风板面积。在某一床高处，气泡的速度可用下式估算()(z)u

Staub提出了计算床层平均空隙率的计算式如下：

u

当u接近u时计算误差较大，当速度较高，接近过渡到湍流方式时，其计算结果令人满意。2.9粒端度当一个颗粒在无限静止的介质中，在重力作用下作自由落体运动时，颗粒在重力作用下加速，降落一段时间后，当流体对颗粒的阻力＋浮力＝颗粒的重量时，颗粒即以等速降落，这个速度称为颗粒的终端速度或自由沉降速度。FF

浮FD

4

uG2

tF浮

6

d

G

6

d

F+F=G6

dg(

dt8

tttt4(ut3C曳力系数C是Re的函数，可分为个区t24层流区Re，=tRet18.5过渡区Re=2~500，=tRet湍流区Re=500~150000，≈0.44t

t

代计算公式，得到球形颗粒终端速度的分区解析式t

tut

gd2pg18

Re<2tAr)0.6667.5dp

g

Rett

()ppgg

Re=500~150000t以上公式是基于球形颗粒，对于非球形颗粒，要对u加以修正，乘以修t

t正系数k后即为非球形颗粒的终端速度。0.67

时当Re<2时，修正系数

[]t

t

当Re修正系数：tt

gdv3(5.31

]dv:等效体积直径，就是与颗粒有相同体积的球形颗粒直径。当Re在－1000之间，k采用插值求得。t

slugmfkckslugmfkck节涌(尺寸D于0.6倍的径D(Geldart,1986)。最小节涌速度u=D为床直径，对于矩形床截面，采用当量直径来代替，D

＝4面积周长。当量当鼓泡流化床的气体流速继续增加时，气泡相互作用加剧，气泡的合并和分裂更加频繁，压力波动的幅度也增大。当气速超过某一气速，压力波动开始减小，大的气泡开始消失；超过另一气速u压力波动的幅度值要比这时床压小得多，同时波动频率非常高，通常已经看不到气泡。操作速度在－范围之内的流化床即为湍流流化床。eq\o\ac(△,P)u

k

u

c目前尚无通用的从鼓泡床向湍流床转变速度的计算公式根据小直径床的实验结果，具有一定的适用范围（J.Grace,1982u3.0

u

适用条件：

d0.7/。

Horio和morishita(1988)提出不同的公式：Ar

rrsslugrrsslugAr

（A4

）

r

（A

）

湍流流化床的运行风速会高于细粒子的终端速度，而低于粗粒子的终端速度。在湍流流态化状态下，气固接触良好，反应器接近理想的返混式反