颗粒物料临界速度和沉降速度的工程计算

碗磷设计与粉体工程

2009年第5期

SP&BMHRELATEDENGINEERING

・1・

颗粒物料临界速度和沉降速度的工程计算

丁德承

(南京三普造粒装备有限公司,江苏南京210041)

摘要:在流化床干燥器(或冷却器)、气流干燥器、流化床吸收器、流化床反应器及气力输送等的工程设计中,颗粒物料的临界速度、沉降速度是主要的工艺参数,计算公式多,但这些公式的局限性较大。介绍了临界速度、沉降速度和操作速度的常用计算公式和计算实例,认为用普拉诺夫斯基修正式进行工程计算不需要进行区间判断,适用于床层孔隙率s=0.4一1.O的任何场合。

关键词:流态化+颗粒物料;临界速度;沉降速度;工程计算;建议中图分类号:TQ028.015

文献标识码:B

文章编号:1009—1904(2009)05—0001—07

1概述

颗粒物料临界速度和沉降速度的工程计算,广泛用于固体物料(或液滴)的流化床干燥器、流化床冷却器、气流干燥器、流化床吸收器(湍球塔吸收器)及流化床反应器等设备的工艺计算。流化床干燥器和气流干燥器的流体力学范围仅限于气固两相,而流化床吸收器的流体力学范围牵涉到气、液、固三相,其工艺计算更为复杂。

颗粒物料临界速度是物料从固定床向流化床过渡的起始速度,即通常所谓的从床层孔隙率占由0.4向大于0.4过渡时的速度。大于l临界速度小于沉降速度就是流态化速度,此时床层孔隙率占=0.55—0.75。占=0.4时的临界速度和8=0.55—0.75时

超过70个计算颗粒物料临界速度的公式。但是,由于实验条件限制在一定范围内,这些公式应用范围的局限性较大。现把工程上常用的计算公式汇集如下。

1、M.李伐提出的公式

G:o.00923

式中:卜颗粒物料的质量流速,ks/(n12・s);卜粒子平均直径,nl;

d1.32盥≮掣竖

(1)

p。——粒子密度,ks/n13;以——气体密度,ks/m3;

r气体动力黏度,Pa・s即ks/(s・m)。

当如大于10时,式(1)要引入修正系数K,其值可按图1确定。如果粒子平均直径d=1/∑(xi/di)时,则式(1)对于多分散颗粒物料也具有足够的精确性。

2、O.M.托杰斯关系式

的流态化速度之比值称流态化数鼠。喷动床的床

层孔隙率占=0.75—0.95。

沉降速度是使用范围很广的一种速度,此时床层孔隙率占一1。颗粒物料(或液滴)在自由空间及有限空间内沉降速度的计算,对于气流干燥器操作速度的选定,流化床干燥器(或冷却器)、流化床反应器上部颗粒物料扬析及带出,气力输送装置中管道速度的确定,除沫器分离空间尺寸的确定是一个基本参数。极为重要。

\、

瞄%¨眈

、\

\、\

、、

2临界速度

在做完了沸腾干燥器物料衡算和热量平衡之后,动力学计算的第一步就是确定被干燥物料粒子的临界速度。

50

100

、、

500

10005000

Re/1

在超过半个世纪的时间里,各国学者曾提出了

图1

Re数与修正系数J|f的关系

硫磷设计与粉体工程

・2・

SP&BMHRELATEDENGlNEERING

2009年第5期

m/s2。

O.M.托杰斯关系式常在杂乱填平、固定床孔隙

率占=0.4时用。

r重力加速度,g=9.8l

3、图解法

Hk。=Rekp∥d

Rek。=Ar/(1

(2)(3)(4)

这是工程计算上比较简捷的计算方法,也是目前被广泛采用的。它由兀.B.李森科提出,并把如3与加之比称为李森科数(巧),即Ly=Re3/Ar,并由图2…求得。图中分别列出了8为O.4,0.5,0.6,0.7,0.8,0.9,1.O的曲线(实线);直线1为李森科数(如)最佳值;曲线2为速度最佳值;曲线(虚线)3,4间的范围为沸腾干燥器的操作区。

400+5.22徊r)

Ar=d3(P。一P。)g/(tJ2p。)

式中:u。。——颗粒物料l临界速度,m/s;

舶k。——临界雷诺数,1;t,——气体运动黏度,m2/s;Ar——阿基米德数,l;

图2李森科数母与阿基米德数Ar之关系

2009年第5期

丁德承.颗粒物料临界速度和沉降速度的工程计算・3・

MkD=7tyo蜀o。/g

(5)

4、通用关联式

近来n.r.罗曼科夫等对文献上已经发表的计算颗粒物料流态化临界速度的公式归纳整理后得出了通用关联式:

“。,/Ⅱ沉降=0.117

5一丁jIf毛粉(6)

式中,Ⅱ沉降为颗粒物料的沉降速度,m/s。

式(6)是对固定床孔隙率s=0.40一O.41而言。但实际上,固定床的孔隙率占值可能在0.35—0.45,那时候式中的0.1175值将为另一个值所代替。此值与占之关系,可由表1查出。表1中A、曰

值随Mkp/“沉降而异,当Ar>106时,取ukp/u沉降=A;

当Ar<106时,取“k。/u沉降=B。

表1系数A和口与固定层孔隙率的关系

5、普拉诺夫斯基通用公式

式中:如——粒子形状修正系数,1;

Re-鱼等掣

r床层孔隙率,1。

式(7)对于8=0.4~1.0范围都是适用的,此公式既能用于计算颗粒物料流态化临界速度U∽也能用于计算流态化速度u漉化和沉降速度u沉降。

6、计算举例

试计算粒子平均直径为2.5mm(最大5.0mm,

最小1.0mm)氯化钾沸腾干燥器的临界速度,主要

工艺条件如下:

进气温度t.800℃排气温度t2125

oC(物料温度)

物料密度P。

2000kg/m3

气体运动黏度t,

247×10“m2/s

气体密度P。=1.293x273/(273+125)

=0.887ks/n13

Ar:塑垒j业。一dSpmg

t,‘pg

v‘pg

2.53×10—9×2×103×9.81

2.472×10一xO.887

=5665

(1)用图解法计算u。。

当s=0.4时查图2,得勿=7×10一。按式(5)计算临界速度。

//'kp=≈如gpm,pE

:沥i而了死■万可矿丐万百—万丽砑丽

=扔8.245×10。

=0.337m/s

(2)用0.M.托杰斯关系式计算配。。RekD=Ar/(1400+5.22√Ar)

=5

665/(1400+5.22/5665)

=3.16“k口=RekpV/d

=3.16×2.47×10—4/(2.5×10—3)

=0.312m/s

(3)利用普拉诺夫斯基通用公式计算Mb取k3=l,占=0.4。

Re:—V/3—67—+百ksA两r8派4"75—-一1

9.15

2————]曩丽■——一

、压万了死石i俨一19.15=(20.975—19.15)/0.588=3.104

ub=Rekpv/d

=3.104×2.47x10—4/(2.5×10—3)

=0.307m_/s

(4)结论

颗粒物料流态化临界速度‰。用上述三种方法的计算结果相差甚小,该差别是工程计算上能够接受的。

3沉降速度

关于球形颗粒在自由空间中沉降速度的确定,文献上已经有大量报道。这些公式的普遍缺点是适用的&数范围很小。所以,为求出沉降速度,必须要经过试算阶段,直到计算得到的沉降雷诺数与公式适用范围一致。现介绍较为简捷的区间判断法及

硫磷设计与粉体工程

・4・

SP&BMHRELATEDENGINEERING

2009年第5期

其应用实例。3.1区间判断法

用区间判断法计算颗粒物料沉降速度应首先算出A,数,继而按式(8)(0.M.托杰斯公式)确定沉降雷诺数廊沉降。

Re沉降=Ar/(18+0.61√五r)

(8)

根据船况降大小,颗粒物料在自由空间的沉降可分为以下三个区域,分别按式(9)、(10)、(11)计算沉降速度。

1、司托克斯区域

Re况降=0.000l一2.0,d=3~100斗m。

配沉降=g(P。一P。)d2/(1札)

(9)

2、阿莱因区域

Re沉降=2~500,d=100—1M沉降=【兰掣】V3(10)

500斗m。

3、牛顿区域

Re况降.=500—150000,d=1500—100000斗m。

u沉降:[3039(p.-pg)d11九(11)u沉降2.

P。

4、应用实例

已知某颗粒物料p。=2010kg/m3、d=0.001m,

试求在p。=0.68

ks/m34z=2.18×10。5

ks/(s・m)、

t,=247×10。m2/s气体中的沉降速度Ⅱ沉降。

则:

A,.:垡:!旦罢二旦s!曼

t,‘pg

一2二010×竺:墨!兰!xlO一9

—2.472×10—8×0.68=475.3

船:——jL=

18+0.61 ̄/Ar

酩沉降=【兰曼掣】V3

18+0.61弘75.3

=15.18(在阿莱因区域)

。惮.42旦2竺纠’3力

【225/40。Jr4×9.81

x20102×1×10-9]1n

一【225×2.18×10—5x0.68J

:泊:丽:7.75

m/s

此计算方法的最大缺点是没有考虑工程计算中的一些因素,例如,粒子形状因素,受阻沉降等。

要求出适合工程实际计算的通用公式,必须解决下列问题:首先在整个胁数范围内计算的精确度;其次,粒子形状因素及其运动受约束条件等影响应当在公式中得到充分反映;且用该式计算时又不复杂。

3.2通用公式(普拉诺夫斯基修正式)

普拉诺夫斯基等分析和整理了各个研究者的实验和计算数据的结果得到:

Re沉降=L{轰i_(12)

鼯/+翔厂一19.

-19.15367式中:k.——考虑到粒子相互作用引起的粒子运动

受约束条件的影响(也就是单位空间颗

粒物料体积分数X的影响),l;后:——考虑到由于管壁(器壁)影响引起的粒子运动受约束条件的系数,l;

后,——考虑粒子形状影响的系数,l。

系数k。即粒子群在有限空间内沉降时,颗粒物料体积分数X对沉降速度的影响:Re沉a-而等赫

.j}I=(1_X)4・75

因此,式(8)就可以写成:

(13)

18+0.61 ̄/Ar(1一X)‘o

式(12)可以写成(不考虑k。、后:、k3的影响):

m沉降=_、/3—67—+Arl(I忑-Xr)4—75-一19.15(14)

在受阻沉降条件下,须考虑气流中固体粒子质量浓度的影响,在气流干燥作业中,由于其质量浓度小于lkg/m3,所以可不必考虑。但是,在气力输送中,就必须考虑气流中颗粒物料体积分数对沉降速度的影响。如以P。=1

200kg/m3,p。=0.45kg/m3,

气体运动黏度u=80.4×10“m2/s为原始条件,计算直径d=1、3、5mm粒子在不同X时的沉降速度,可以看出X对沉降速度有很大的影响。

笔者用式(13)、式(14)分别计算d=1、3、5

him,X=0、0.Ol、0.05、0.10、0.20、0.30

m3/m3时的

船沉降(见表2),然后用式(2)算得它们的M沉降(见表

3)。根据笔者计算,在x等于0.Olm3/m3时,沉降速度较单一粒子在自由空间中约低3%,当x达到

0.05

m3/m3时,则要低12%一18%。

系数k:即管壁(或器壁)对粒子沉降速度的影

响,可以用下式表示:

k2=1/M2

式中,肘为与粒子直径和管径(设备直径)之比及与

2009年第5期

丁德承.颗粒物料临界速度和沉降速度的工程计算’

・5・

表2

Re数计算结果汇总

表3Ⅳ沉奠计算结果汇总

流动状态有关的修正乘积,1。除个别情况外,在工业实际中,颗粒直径均小于管径的10%,所以管壁对沉降速度的影响可以忽略,即M=1。

系数J},是粒子形状修正系数:||}3=1l—10饥

其中/为粒子形状因素,1。各种物料的形状系数见上海海运学院气力输送小组1977年发表的论文,现列于图3~4。

当球形粒子在无限空间沉降时,也即X=0,I|},=1,M=1,蠡3=l/=1,则式(12)就变成

Reg降=丛等簪堂

(15)

但是,在工程实际中,很少碰到这种理想状况,

也即多数是非球形颗粒群在有限空间内沉降。3.3计算举例

试求球形颗粒直径d=1.2×10一m、密度P。=

2640

kg/m3,在900℃烟道气中的自由沉降速度。

烟气密度P。=0.318kr,/m3,运动黏度l,=152.5X

10山n12/s,动力黏度/x=4.68×10。5kr,/(8・m)。

则:Ar:塑乌!出。一d3pmgv‘pg

u。pg

=幽裂警案警=605152

.。×2

10-12×.3180

1、区间判断法

鼋手裴冬璺

圆板片

圆柱体

七3=1.00

如=2.50

网球

表面不光滑的近似圜球

k3=3.00

椭圆体

匿|岛----5

|!|00

}蓬

方板片

图3颗粒物料形状系数b

他沉降2i石而;2i五百而亓

‰《掣卜1/3【.薪2

=92.45(在阿莱因区域)2d3】I/3

r4×9.812×2.642×1.23×10—311/3

一【225

x4.68×0.318x10—5

=浙了百砑=11.14m/s

2、0.M.托杰斯关系式

Jk沉降t,

92.45×152.5×10—6

Ⅱ沉降一d

1.2×10—3

=11.75Il∥8

3、普扭诺夫斯基修正式

m沉降=—V/-3—6ff瓦+A丽r-一19.15

:—v/3670.588

+60—51-19.15:103.68

胁沉降t,103.68

X152.5×10—6

Ⅱ沉降一d

1.2×10—3

=13.18m/8

硫磷设计与粉体工程

・6・

SP&BMHRELATEDENGINEERING

2009年第5期

坷三臼

b=1

00

圆球

佃●}-一●}—一

£:D=2

k3=1.66

圆柱体

訾陲譬

嗡三篚

口:b--O.9

b=I.76

kt=l86

棱形体正方形体

霉睡

h:a--025“=3.23

觑=1,】7

正方形板片

不规则球体

口:b=071

岛=I,06

“=2

27

不规则椭圆体不规则块体

图4实验测定的颗粒物料形状系数

4、结论

计算结果表明,用各种方法得到的颗粒物料沉降速度基本一致。如果计算球形粒子在自由空间沉降速度,可以利用式(9)、(10)、(11)直接求出。而托杰斯公式具有应用方便、区域广泛.数据可靠的优点。但是,从工程计算观点出发,托杰斯公式是不能满足需要的。因为在工程中碰到的物料都是不规则而需要引入形状系数,其次在工程中都是颗粒群在有限空间做自由沉降。普拉诺夫斯基修正式具有更大灵活性,既可以做单一粒子在无限空间自由沉降的计算,也可以进行不规则粒子群在有限空间内受阻沉降计算,适宜于工程设计计算。

4用通用公式计算实例

已知:球形粒子直径d=5toni、密度4000kg/m3、

流化空气温度215℃,气体密度0.723ks/m3,运动

黏度抄=3.39×10。5

m2/s。

应用普拉诺夫斯基通用公式计算颗粒物料的临界速度、吹失速度(沉降速度)和流态化速度(操作

速度)。

l、计算阿基米德数

Ar:丛生?出;一d3pmg

v‘pg

v。pg

4000×9.81×53×10—93.392×10—10×0.723=5.9×106

2、计算uk。

(1)用托杰斯关系式计算

尺e-一2而试iAF云石

5.9×105

400+5.22朽.9×106

=419

‰,2亍2——百面r一

尺e沉降”419×3.39×10—5=2.84m/s

(2)用普拉诺夫斯基通用式计算取占=0。4,缸=l。

Rekp=—、/3—67—+百k3A葡re丽4"7-5-一19.152—————百两河———一

:—、/7、压万i两i而畈百百蒋一19.15

6336.—1-19.15:468

0.588

u。一2丁2——百面厂

Rekpt,468×3.39×10—5

=3.17m/s

(3)两者误差

6=(3.17—2.84)/3.17=0.104=10.4%

3、计算H沉降

(1)用托杰斯关系式计算

m沉降2i五h丽r尹

5.9×106

18-I-0.61 ̄/5.9×106

139

2009年第5期

丁德承.颗粒物料临界速度和沉降速度的工程计算・7・

=业器=3

_

499

69

=’7.‘7lm/s

(2)用普拉诺夫斯基通用式计算

“沉降2丁2——1而面一

Re沉降p

3934×3.39×10—5

=26.67m/s

(2)用普拉诺夫斯基修正式计算

代e沉降2————1Il两函—一

√367+Ar—19.15

如流化=—v/3—67—+百Ar6丽4"7r5-一19.152—————弋焉丽———一

历石i曩丽一19.152———■F丙矿—~=2可4丽09.92=4

098.5

R%t降v

4098.5×3.39×10—5

:—,/521646.—3-19.15:1

0.588

√7j67+5.9×106×0.64・75—19.15

196

如沉降u

Ⅱ流化一d

1196×3.39×10—5

5×10—3

=8.11m/s

Ⅱ沉降2丁2——了i矿广

=27.79m/s

(3)两者误差

6=(8.11—7.71)/8.11=0.049

3=4.93%

5结束语

笔者通过长期工程设计,推介既可进行颗粒物料临界速度计算,又可进行操作速度和沉降速度计

(3)两者误差

6=(27.79—26.67)/27.79=0.0403=4.03%4、计算操作速度扯流化(设床层孔隙率占=0.6)(1)用托杰斯关系式计算

算的通用公式——普拉诺夫斯基修正式。用此公式

做工程计算不需要进行区间判断,适用于床层孑L隙占=0.4~1.0的任何场合。

参考文献:

[1]童景山,张克.流态化干燥技术[M].北京:中国建筑工业出

胎龇2

i五而

Ar苫

5.9×106×0.64・75

18+0.61撕.9

521279.3

106NO.64・75

版社,1985:61.

=——————————二======

18+0.61 ̄/521279.3

=l137

作者简介:丁德承(1936一),男,浙江定海人,化工工艺工程

师,从事磷复肥工艺设计和技术开发,我国第一套30kt/a、150kt/a料浆法磷铵装置设计的项目负责人。

(收稿日期:2008—05—08)

“漉化2—产2——五而r

尺8沉降t,

1137×3.39×10—5

'●’’,'’’’'',’','’l-’’’I,'’’’',,’’’,,’,,’’’,,,'’’'',’’,’'●'’,'’',,’’,’'’’’’,’,’,,,'’’●'',●,'’

正反浮选中低品位磷矿技术

磷矿石是磷化工工业的主要原料,特别对磷肥生产,关系到农业生产和粮食安全。我国磷矿石资源主要分布在西部和中部地区,而且中低品位磷矿多,高品位磷矿少。长期以来,我国磷化工和磷肥生产均采用高品位磷矿,特别是不科学、不合理地乱开滥采,造成高品位磷矿资源的浪费和损失。

我国多次强调发展经济必须立足可持续发展、爱护资源、节约资源。为此中低品位磷矿的利用提上议事日程,众多的科研院校、公司企业纷纷致力于中低品位磷矿的工艺研发和工业应用。

湖北洋丰中磷矿业有限公司与院校合作,对湖北放马山中低品位磷矿进行选矿技术的攻关。经过对各种工艺进行反复对比试验——选矿剂筛选和改良处理,采用正反浮选技术,使|t,(P205)为15.8%、埘(Mgo)为6.06%的原矿,经选矿

缓解我国高品位磷矿短缺

新工艺获得W(P205)为32.27%、”(Mgo)为1.42%的磷精矿,可以满足生产MAP、DAP的要求,回收率达到81.30%,通过中试考核,为工业生产提供了可行的科学依据。

该公司成立了1.2Mt/a磷矿石采选工程项目,总投资2.6亿元,其中固定资产投资2.3亿元,可以处理埘(P20,)为15.8%的原矿1.2Mt/a,年产30%P20;磷精矿500左右。

此项技术的成功应用,可使我国中低品位磷矿成为高品位磷精矿,顺利生产MAP、DAP和高浓度复混肥料。不仅为湖北省占有总量93%的中低品位磷矿石实现工业化的应用,而且也为全国利用中低品位磷矿提供技术依据和支撑,缓解我国短缺高品位磷矿资源的困境。

(冯元琦)

kt。

项目建成后,可实现年销售收入3.5亿元,就业人员1000名

颗粒物料临界速度和沉降速度的工程计算

作者:作者单位:刊名:英文刊名:年,卷(期):

丁德承, DING De-cheng

南京三普造粒装备有限公司,江苏,南京,210041

硫磷设计与粉体工程

SULPHUR PHOSPHORUS & BULK MATERIALS HANDLING RELATED ENGINEERING2009(5)

参考文献(1条)

1. 童景山;张克 流态化干燥技术 1985

本文读者也读过(10条)

1. 吴宁. 张琪. 曲占庆 固体颗粒在液体中沉降速度的计算方法评述[期刊论文]-石油钻采工艺2000,22(2)2. 非球形颗粒及颗粒群的自由沉降速度[期刊论文]-中国粉体技术2005,11(z1)3. 郑邦民. 夏军强 固体颗粒的群体沉降速度分析[期刊论文]-泥沙研究2004(6)

4. 曾庆川. 刘小兵. ZENG Qing-chuan. LIU Xiao-bing 紊流中颗粒沉降速度的测量方法[期刊论文]-四川工业学院学报2000,19(4)

5. 谢洪勇. 张大为. 张礼鸣 非球形颗粒及颗粒群的自由沉降速度[会议论文]-6. 孙远广 液固两相流固相浓度及沉降速度的测量研究[学位论文]2009

7. 付峥嵘. 李念平. 王汉青. FU Zheng-rong. LI Nian-ping. WANG Han-qing 风管中颗粒物沉降速度的解析法预测模型[期刊论文]-湖南大学学报(自然科学版)2008,35(2)

8. 梁政. 任连城. 吴世辉. LIANG-Zheng . REN Lian-cheng. WU Shi-hui 水力旋流器颗粒沉降速度与分离粒度研究[期刊论文]-西南石油大学学报2007,29(3)

9. H.Kalman. E.Rabinovich 颗粒-流体系统中的临界速度[会议论文]-2007

10. Irina V.Putilova. Viacheslav Y.Putilov 对粉体气力输送装置中管道腐蚀磨损的估算[会议论文]-2007

引证文献(1条)

1. 张尉杰. 丁德承 磷肥和磷复肥生产用冷却器的工艺计算[期刊论文]-硫磷设计与粉体工程 2011(4)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_llsjyftgc200905001.aspx

碗磷设计与粉体工程

2009年第5期

SP&BMHRELATEDENGINEERING

・1・

颗粒物料临界速度和沉降速度的工程计算

丁德承

(南京三普造粒装备有限公司,江苏南京210041)

摘要:在流化床干燥器(或冷却器)、气流干燥器、流化床吸收器、流化床反应器及气力输送等的工程设计中,颗粒物料的临界速度、沉降速度是主要的工艺参数,计算公式多,但这些公式的局限性较大。介绍了临界速度、沉降速度和操作速度的常用计算公式和计算实例,认为用普拉诺夫斯基修正式进行工程计算不需要进行区间判断,适用于床层孔隙率s=0.4一1.O的任何场合。

关键词:流态化+颗粒物料;临界速度;沉降速度;工程计算;建议中图分类号:TQ028.015

文献标识码:B

文章编号:1009—1904(2009)05—0001—07

1概述

颗粒物料临界速度和沉降速度的工程计算,广泛用于固体物料(或液滴)的流化床干燥器、流化床冷却器、气流干燥器、流化床吸收器(湍球塔吸收器)及流化床反应器等设备的工艺计算。流化床干燥器和气流干燥器的流体力学范围仅限于气固两相,而流化床吸收器的流体力学范围牵涉到气、液、固三相,其工艺计算更为复杂。

颗粒物料临界速度是物料从固定床向流化床过渡的起始速度,即通常所谓的从床层孔隙率占由0.4向大于0.4过渡时的速度。大于l临界速度小于沉降速度就是流态化速度,此时床层孔隙率占=0.55—0.75。占=0.4时的临界速度和8=0.55—0.75时

超过70个计算颗粒物料临界速度的公式。但是,由于实验条件限制在一定范围内,这些公式应用范围的局限性较大。现把工程上常用的计算公式汇集如下。

1、M.李伐提出的公式

G:o.00923

式中:卜颗粒物料的质量流速,ks/(n12・s);卜粒子平均直径,nl;

d1.32盥≮掣竖

(1)

p。——粒子密度,ks/n13;以——气体密度,ks/m3;

r气体动力黏度,Pa・s即ks/(s・m)。

当如大于10时,式(1)要引入修正系数K,其值可按图1确定。如果粒子平均直径d=1/∑(xi/di)时,则式(1)对于多分散颗粒物料也具有足够的精确性。

2、O.M.托杰斯关系式

的流态化速度之比值称流态化数鼠。喷动床的床

层孔隙率占=0.75—0.95。

沉降速度是使用范围很广的一种速度,此时床层孔隙率占一1。颗粒物料(或液滴)在自由空间及有限空间内沉降速度的计算,对于气流干燥器操作速度的选定,流化床干燥器(或冷却器)、流化床反应器上部颗粒物料扬析及带出,气力输送装置中管道速度的确定,除沫器分离空间尺寸的确定是一个基本参数。极为重要。

\、

瞄%¨眈

、\

\、\

、、

2临界速度

在做完了沸腾干燥器物料衡算和热量平衡之后,动力学计算的第一步就是确定被干燥物料粒子的临界速度。

50

100

、、

500

10005000

Re/1

在超过半个世纪的时间里,各国学者曾提出了

图1

Re数与修正系数J|f的关系

硫磷设计与粉体工程

・2・

SP&BMHRELATEDENGlNEERING

2009年第5期

m/s2。

O.M.托杰斯关系式常在杂乱填平、固定床孔隙

率占=0.4时用。

r重力加速度,g=9.8l

3、图解法

Hk。=Rekp∥d

Rek。=Ar/(1

(2)(3)(4)

这是工程计算上比较简捷的计算方法,也是目前被广泛采用的。它由兀.B.李森科提出,并把如3与加之比称为李森科数(巧),即Ly=Re3/Ar,并由图2…求得。图中分别列出了8为O.4,0.5,0.6,0.7,0.8,0.9,1.O的曲线(实线);直线1为李森科数(如)最佳值;曲线2为速度最佳值;曲线(虚线)3,4间的范围为沸腾干燥器的操作区。

400+5.22徊r)

Ar=d3(P。一P。)g/(tJ2p。)

式中:u。。——颗粒物料l临界速度,m/s;

舶k。——临界雷诺数,1;t,——气体运动黏度,m2/s;Ar——阿基米德数,l;

图2李森科数母与阿基米德数Ar之关系

2009年第5期

丁德承.颗粒物料临界速度和沉降速度的工程计算・3・

MkD=7tyo蜀o。/g

(5)

4、通用关联式

近来n.r.罗曼科夫等对文献上已经发表的计算颗粒物料流态化临界速度的公式归纳整理后得出了通用关联式:

“。,/Ⅱ沉降=0.117

5一丁jIf毛粉(6)

式中,Ⅱ沉降为颗粒物料的沉降速度,m/s。

式(6)是对固定床孔隙率s=0.40一O.41而言。但实际上,固定床的孔隙率占值可能在0.35—0.45,那时候式中的0.1175值将为另一个值所代替。此值与占之关系,可由表1查出。表1中A、曰

值随Mkp/“沉降而异,当Ar>106时,取ukp/u沉降=A;

当Ar<106时,取“k。/u沉降=B。

表1系数A和口与固定层孔隙率的关系

5、普拉诺夫斯基通用公式

式中:如——粒子形状修正系数,1;

Re-鱼等掣

r床层孔隙率,1。

式(7)对于8=0.4~1.0范围都是适用的,此公式既能用于计算颗粒物料流态化临界速度U∽也能用于计算流态化速度u漉化和沉降速度u沉降。

6、计算举例

试计算粒子平均直径为2.5mm(最大5.0mm,

最小1.0mm)氯化钾沸腾干燥器的临界速度,主要

工艺条件如下:

进气温度t.800℃排气温度t2125

oC(物料温度)

物料密度P。

2000kg/m3

气体运动黏度t,

247×10“m2/s

气体密度P。=1.293x273/(273+125)

=0.887ks/n13

Ar:塑垒j业。一dSpmg

t,‘pg

v‘pg

2.53×10—9×2×103×9.81

2.472×10一xO.887

=5665

(1)用图解法计算u。。

当s=0.4时查图2,得勿=7×10一。按式(5)计算临界速度。

//'kp=≈如gpm,pE

:沥i而了死■万可矿丐万百—万丽砑丽

=扔8.245×10。

=0.337m/s

(2)用0.M.托杰斯关系式计算配。。RekD=Ar/(1400+5.22√Ar)

=5

665/(1400+5.22/5665)

=3.16“k口=RekpV/d

=3.16×2.47×10—4/(2.5×10—3)

=0.312m/s

(3)利用普拉诺夫斯基通用公式计算Mb取k3=l,占=0.4。

Re:—V/3—67—+百ksA两r8派4"75—-一1

9.15

2————]曩丽■——一

、压万了死石i俨一19.15=(20.975—19.15)/0.588=3.104

ub=Rekpv/d

=3.104×2.47x10—4/(2.5×10—3)

=0.307m_/s

(4)结论

颗粒物料流态化临界速度‰。用上述三种方法的计算结果相差甚小,该差别是工程计算上能够接受的。

3沉降速度

关于球形颗粒在自由空间中沉降速度的确定,文献上已经有大量报道。这些公式的普遍缺点是适用的&数范围很小。所以,为求出沉降速度,必须要经过试算阶段,直到计算得到的沉降雷诺数与公式适用范围一致。现介绍较为简捷的区间判断法及

硫磷设计与粉体工程

・4・

SP&BMHRELATEDENGINEERING

2009年第5期

其应用实例。3.1区间判断法

用区间判断法计算颗粒物料沉降速度应首先算出A,数,继而按式(8)(0.M.托杰斯公式)确定沉降雷诺数廊沉降。

Re沉降=Ar/(18+0.61√五r)

(8)

根据船况降大小,颗粒物料在自由空间的沉降可分为以下三个区域,分别按式(9)、(10)、(11)计算沉降速度。

1、司托克斯区域

Re况降=0.000l一2.0,d=3~100斗m。

配沉降=g(P。一P。)d2/(1札)

(9)

2、阿莱因区域

Re沉降=2~500,d=100—1M沉降=【兰掣】V3(10)

500斗m。

3、牛顿区域

Re况降.=500—150000,d=1500—100000斗m。

u沉降:[3039(p.-pg)d11九(11)u沉降2.

P。

4、应用实例

已知某颗粒物料p。=2010kg/m3、d=0.001m,

试求在p。=0.68

ks/m34z=2.18×10。5

ks/(s・m)、

t,=247×10。m2/s气体中的沉降速度Ⅱ沉降。

则:

A,.:垡:!旦罢二旦s!曼

t,‘pg

一2二010×竺:墨!兰!xlO一9

—2.472×10—8×0.68=475.3

船:——jL=

18+0.61 ̄/Ar

酩沉降=【兰曼掣】V3

18+0.61弘75.3

=15.18(在阿莱因区域)

。惮.42旦2竺纠’3力

【225/40。Jr4×9.81

x20102×1×10-9]1n

一【225×2.18×10—5x0.68J

:泊:丽:7.75

m/s

此计算方法的最大缺点是没有考虑工程计算中的一些因素,例如,粒子形状因素,受阻沉降等。

要求出适合工程实际计算的通用公式,必须解决下列问题:首先在整个胁数范围内计算的精确度;其次,粒子形状因素及其运动受约束条件等影响应当在公式中得到充分反映;且用该式计算时又不复杂。

3.2通用公式(普拉诺夫斯基修正式)

普拉诺夫斯基等分析和整理了各个研究者的实验和计算数据的结果得到:

Re沉降=L{轰i_(12)

鼯/+翔厂一19.

-19.15367式中:k.——考虑到粒子相互作用引起的粒子运动

受约束条件的影响(也就是单位空间颗

粒物料体积分数X的影响),l;后:——考虑到由于管壁(器壁)影响引起的粒子运动受约束条件的系数,l;

后,——考虑粒子形状影响的系数,l。

系数k。即粒子群在有限空间内沉降时,颗粒物料体积分数X对沉降速度的影响:Re沉a-而等赫

.j}I=(1_X)4・75

因此,式(8)就可以写成:

(13)

18+0.61 ̄/Ar(1一X)‘o

式(12)可以写成(不考虑k。、后:、k3的影响):

m沉降=_、/3—67—+Arl(I忑-Xr)4—75-一19.15(14)

在受阻沉降条件下,须考虑气流中固体粒子质量浓度的影响,在气流干燥作业中,由于其质量浓度小于lkg/m3,所以可不必考虑。但是,在气力输送中,就必须考虑气流中颗粒物料体积分数对沉降速度的影响。如以P。=1

200kg/m3,p。=0.45kg/m3,

气体运动黏度u=80.4×10“m2/s为原始条件,计算直径d=1、3、5mm粒子在不同X时的沉降速度,可以看出X对沉降速度有很大的影响。

笔者用式(13)、式(14)分别计算d=1、3、5

him,X=0、0.Ol、0.05、0.10、0.20、0.30

m3/m3时的

船沉降(见表2),然后用式(2)算得它们的M沉降(见表

3)。根据笔者计算,在x等于0.Olm3/m3时,沉降速度较单一粒子在自由空间中约低3%,当x达到

0.05

m3/m3时,则要低12%一18%。

系数k:即管壁(或器壁)对粒子沉降速度的影

响,可以用下式表示:

k2=1/M2

式中,肘为与粒子直径和管径(设备直径)之比及与

2009年第5期

丁德承.颗粒物料临界速度和沉降速度的工程计算’

・5・

表2

Re数计算结果汇总

表3Ⅳ沉奠计算结果汇总

流动状态有关的修正乘积,1。除个别情况外,在工业实际中,颗粒直径均小于管径的10%,所以管壁对沉降速度的影响可以忽略,即M=1。

系数J},是粒子形状修正系数:||}3=1l—10饥

其中/为粒子形状因素,1。各种物料的形状系数见上海海运学院气力输送小组1977年发表的论文,现列于图3~4。

当球形粒子在无限空间沉降时,也即X=0,I|},=1,M=1,蠡3=l/=1,则式(12)就变成

Reg降=丛等簪堂

(15)

但是,在工程实际中,很少碰到这种理想状况,

也即多数是非球形颗粒群在有限空间内沉降。3.3计算举例

试求球形颗粒直径d=1.2×10一m、密度P。=

2640

kg/m3,在900℃烟道气中的自由沉降速度。

烟气密度P。=0.318kr,/m3,运动黏度l,=152.5X

10山n12/s,动力黏度/x=4.68×10。5kr,/(8・m)。

则:Ar:塑乌!出。一d3pmgv‘pg

u。pg

=幽裂警案警=605152

.。×2

10-12×.3180

1、区间判断法

鼋手裴冬璺

圆板片

圆柱体

七3=1.00

如=2.50

网球

表面不光滑的近似圜球

k3=3.00

椭圆体

匿|岛----5

|!|00

}蓬

方板片

图3颗粒物料形状系数b

他沉降2i石而;2i五百而亓

‰《掣卜1/3【.薪2

=92.45(在阿莱因区域)2d3】I/3

r4×9.812×2.642×1.23×10—311/3

一【225

x4.68×0.318x10—5

=浙了百砑=11.14m/s

2、0.M.托杰斯关系式

Jk沉降t,

92.45×152.5×10—6

Ⅱ沉降一d

1.2×10—3

=11.75Il∥8

3、普扭诺夫斯基修正式

m沉降=—V/-3—6ff瓦+A丽r-一19.15

:—v/3670.588

+60—51-19.15:103.68

胁沉降t,103.68

X152.5×10—6

Ⅱ沉降一d

1.2×10—3

=13.18m/8

硫磷设计与粉体工程

・6・

SP&BMHRELATEDENGINEERING

2009年第5期

坷三臼

b=1

00

圆球

佃●}-一●}—一

£:D=2

k3=1.66

圆柱体

訾陲譬

嗡三篚

口:b--O.9

b=I.76

kt=l86

棱形体正方形体

霉睡

h:a--025“=3.23

觑=1,】7

正方形板片

不规则球体

口:b=071

岛=I,06

“=2

27

不规则椭圆体不规则块体

图4实验测定的颗粒物料形状系数

4、结论

计算结果表明,用各种方法得到的颗粒物料沉降速度基本一致。如果计算球形粒子在自由空间沉降速度,可以利用式(9)、(10)、(11)直接求出。而托杰斯公式具有应用方便、区域广泛.数据可靠的优点。但是,从工程计算观点出发,托杰斯公式是不能满足需要的。因为在工程中碰到的物料都是不规则而需要引入形状系数,其次在工程中都是颗粒群在有限空间做自由沉降。普拉诺夫斯基修正式具有更大灵活性,既可以做单一粒子在无限空间自由沉降的计算,也可以进行不规则粒子群在有限空间内受阻沉降计算,适宜于工程设计计算。

4用通用公式计算实例

已知:球形粒子直径d=5toni、密度4000kg/m3、

流化空气温度215℃,气体密度0.723ks/m3,运动

黏度抄=3.39×10。5

m2/s。

应用普拉诺夫斯基通用公式计算颗粒物料的临界速度、吹失速度(沉降速度)和流态化速度(操作

速度)。

l、计算阿基米德数

Ar:丛生?出;一d3pmg

v‘pg

v。pg

4000×9.81×53×10—93.392×10—10×0.723=5.9×106

2、计算uk。

(1)用托杰斯关系式计算

尺e-一2而试iAF云石

5.9×105

400+5.22朽.9×106

=419

‰,2亍2——百面r一

尺e沉降”419×3.39×10—5=2.84m/s

(2)用普拉诺夫斯基通用式计算取占=0。4,缸=l。

Rekp=—、/3—67—+百k3A葡re丽4"7-5-一19.152—————百两河———一

:—、/7、压万i两i而畈百百蒋一19.15

6336.—1-19.15:468

0.588

u。一2丁2——百面厂

Rekpt,468×3.39×10—5

=3.17m/s

(3)两者误差

6=(3.17—2.84)/3.17=0.104=10.4%

3、计算H沉降

(1)用托杰斯关系式计算

m沉降2i五h丽r尹

5.9×106

18-I-0.61 ̄/5.9×106

139

2009年第5期

丁德承.颗粒物料临界速度和沉降速度的工程计算・7・

=业器=3

_

499

69

=’7.‘7lm/s

(2)用普拉诺夫斯基通用式计算

“沉降2丁2——1而面一

Re沉降p

3934×3.39×10—5

=26.67m/s

(2)用普拉诺夫斯基修正式计算

代e沉降2————1Il两函—一

√367+Ar—19.15

如流化=—v/3—67—+百Ar6丽4"7r5-一19.152—————弋焉丽———一

历石i曩丽一19.152———■F丙矿—~=2可4丽09.92=4

098.5

R%t降v

4098.5×3.39×10—5

:—,/521646.—3-19.15:1

0.588

√7j67+5.9×106×0.64・75—19.15

196

如沉降u

Ⅱ流化一d

1196×3.39×10—5

5×10—3

=8.11m/s

Ⅱ沉降2丁2——了i矿广

=27.79m/s

(3)两者误差

6=(8.11—7.71)/8.11=0.049

3=4.93%

5结束语

笔者通过长期工程设计,推介既可进行颗粒物料临界速度计算,又可进行操作速度和沉降速度计

(3)两者误差

6=(27.79—26.67)/27.79=0.0403=4.03%4、计算操作速度扯流化(设床层孔隙率占=0.6)(1)用托杰斯关系式计算

算的通用公式——普拉诺夫斯基修正式。用此公式

做工程计算不需要进行区间判断,适用于床层孑L隙占=0.4~1.0的任何场合。

参考文献:

[1]童景山,张克.流态化干燥技术[M].北京:中国建筑工业出

胎龇2

i五而

Ar苫

5.9×106×0.64・75

18+0.61撕.9

521279.3

106NO.64・75

版社,1985:61.

=——————————二======

18+0.61 ̄/521279.3

=l137

作者简介:丁德承(1936一),男,浙江定海人,化工工艺工程

师,从事磷复肥工艺设计和技术开发,我国第一套30kt/a、150kt/a料浆法磷铵装置设计的项目负责人。

(收稿日期:2008—05—08)

“漉化2—产2——五而r

尺8沉降t,

1137×3.39×10—5

'●’’,'’’’'',’','’l-’’’I,'’’’',,’’’,,’,,’’’,,,'’’'',’’,’'●'’,'’',,’’,’'’’’’,’,’,,,'’’●'',●,'’

正反浮选中低品位磷矿技术

磷矿石是磷化工工业的主要原料,特别对磷肥生产,关系到农业生产和粮食安全。我国磷矿石资源主要分布在西部和中部地区,而且中低品位磷矿多,高品位磷矿少。长期以来,我国磷化工和磷肥生产均采用高品位磷矿,特别是不科学、不合理地乱开滥采,造成高品位磷矿资源的浪费和损失。

我国多次强调发展经济必须立足可持续发展、爱护资源、节约资源。为此中低品位磷矿的利用提上议事日程,众多的科研院校、公司企业纷纷致力于中低品位磷矿的工艺研发和工业应用。

湖北洋丰中磷矿业有限公司与院校合作,对湖北放马山中低品位磷矿进行选矿技术的攻关。经过对各种工艺进行反复对比试验——选矿剂筛选和改良处理,采用正反浮选技术,使|t,(P205)为15.8%、埘(Mgo)为6.06%的原矿,经选矿

缓解我国高品位磷矿短缺

新工艺获得W(P205)为32.27%、”(Mgo)为1.42%的磷精矿,可以满足生产MAP、DAP的要求,回收率达到81.30%,通过中试考核,为工业生产提供了可行的科学依据。

该公司成立了1.2Mt/a磷矿石采选工程项目,总投资2.6亿元,其中固定资产投资2.3亿元,可以处理埘(P20,)为15.8%的原矿1.2Mt/a,年产30%P20;磷精矿500左右。

此项技术的成功应用,可使我国中低品位磷矿成为高品位磷精矿,顺利生产MAP、DAP和高浓度复混肥料。不仅为湖北省占有总量93%的中低品位磷矿石实现工业化的应用,而且也为全国利用中低品位磷矿提供技术依据和支撑,缓解我国短缺高品位磷矿资源的困境。

(冯元琦)

kt。

项目建成后,可实现年销售收入3.5亿元,就业人员1000名

颗粒物料临界速度和沉降速度的工程计算

作者:作者单位:刊名:英文刊名:年,卷(期):

丁德承, DING De-cheng

南京三普造粒装备有限公司,江苏,南京,210041

硫磷设计与粉体工程

SULPHUR PHOSPHORUS & BULK MATERIALS HANDLING RELATED ENGINEERING2009(5)

参考文献(1条)

1. 童景山;张克 流态化干燥技术 1985

本文读者也读过(10条)

1. 吴宁. 张琪. 曲占庆 固体颗粒在液体中沉降速度的计算方法评述[期刊论文]-石油钻采工艺2000,22(2)2. 非球形颗粒及颗粒群的自由沉降速度[期刊论文]-中国粉体技术2005,11(z1)3. 郑邦民. 夏军强 固体颗粒的群体沉降速度分析[期刊论文]-泥沙研究2004(6)

4. 曾庆川. 刘小兵. ZENG Qing-chuan. LIU Xiao-bing 紊流中颗粒沉降速度的测量方法[期刊论文]-四川工业学院学报2000,19(4)

5. 谢洪勇. 张大为. 张礼鸣 非球形颗粒及颗粒群的自由沉降速度[会议论文]-6. 孙远广 液固两相流固相浓度及沉降速度的测量研究[学位论文]2009

7. 付峥嵘. 李念平. 王汉青. FU Zheng-rong. LI Nian-ping. WANG Han-qing 风管中颗粒物沉降速度的解析法预测模型[期刊论文]-湖南大学学报(自然科学版)2008,35(2)

8. 梁政. 任连城. 吴世辉. LIANG-Zheng . REN Lian-cheng. WU Shi-hui 水力旋流器颗粒沉降速度与分离粒度研究[期刊论文]-西南石油大学学报2007,29(3)

9. H.Kalman. E.Rabinovich 颗粒-流体系统中的临界速度[会议论文]-2007

10. Irina V.Putilova. Viacheslav Y.Putilov 对粉体气力输送装置中管道腐蚀磨损的估算[会议论文]-2007

引证文献(1条)

1. 张尉杰. 丁德承 磷肥和磷复肥生产用冷却器的工艺计算[期刊论文]-硫磷设计与粉体工程 2011(4)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_llsjyftgc200905001.aspx


相关文章

  • 卧螺离心机生产能力的FCPS理论_梁毅
  • 第 51 卷第 1 期2014 年 2 月化 工 设 备 与 管 道PROCESS EQUIPMENT & PIPINGVol. 51 No. 1Feb. 2014 卧螺离心机生产能力的FCPS理论 梁毅1,孙恒2,王东琪2,何雨谦 ...查看


  • 粉体工程与设备复习题
  • 粉体工程习题 一.选择题(以下各小题均有4或3个备选答案,请圈出唯一正确的答案) 1.RRB粒度分布方程中的n是 . A.功指数 B.旋涡指数 C.均匀性指数 D.时间指数 2.粒度分析中常采用RR坐标来绘制粒度分布曲线.该坐标的横坐标为颗 ...查看


  • 生物工程设备习题集
  • 生物工程设备习题集 一. 单项选择题: (每题1分) 1. 目前啤酒厂的糖化锅中利用_____D_____进行搅拌. A.圆盘平直叶涡轮搅拌器 B.螺旋浆式搅拌器 C.醪液内二氧化碳的密度梯度 D. 二折叶旋桨式搅拌器 2. 3. 4. 空 ...查看


  • 化工原理第四版思考题答案
  • 3-2 球形颗粒在流体中从静止开始沉降,经历哪两个阶段?何谓固体颗粒在流体中的沉降速度?沉降速度受那些因素的影响? 答:1.加速阶段和匀速阶段:2.颗粒手里平衡时,匀速阶段中颗粒相对于流体的运动速度t3.影响因素由沉降公式确定dp.p. ...查看


  • 陶瓷机械习题集和答案
  • 陶瓷工业机械与设备习题集 一.选择题 1. 当研磨体在脱离角处具有最大能量时,脱离角应是: A A.54度44分:B.45度:C.60度45分:D.90度 2. 球磨机临界转速的定义: C A. 最外层研磨体获得最大粉碎能量时的转速:B. ...查看


  • [生物分离工程]课程教学大纲
  • <生物分离工程>课程教学大纲 课程编号:02062 英文名称:Bioseparations Engineering 一.课程说明 1. 课程类别 专业课程 2. 适应专业及课程性质 制药工程专业 必修 3. 课程目的 本课程主要 ...查看


  • 循环流化床锅炉原理
  • 循环流化床锅炉基本原理 肖 平 西安热工研究院有限公司 循环流化床锅炉技术部 2008-9-21 第一章.概述 第一节 什么叫循环流化床 l 国外循环流化床锅炉 l 1979年,芬兰奥斯龙(Ahlsltrom)公司,世界首台20t/h CF ...查看


  • 关于流化床技术
  • 关于流化床技术 1 什么是制药机械?与其它机械的区别 1.1 制药机械是完成和辅助完成制药工艺的生产机械与设备即将原料及辅料按生产工艺要求,制成可以直接用于临床医疗或作为药品制成品的机械与设备.它包括原料药机械及设备.制剂机械.药用粉碎机械 ...查看


  • 食品工程原理
  • 食品工程原理复习 1.单元操作与三传理论的概念及关系. 2.粘度的概念及牛顿内摩擦(粘性)定律.牛顿黏性定律的数学表达式是dudy ,服从此定律的流体称为牛顿 流体. 3.理想流体的概念及意义. 4.热力体系:指某一由周围边界所限定 ...查看


热门内容