脉动流化床内颗粒分离行为的数值模拟

第39卷,总第225期2021年1月,第1期
《节能技术》
ENERGY CONSERVATION TECHNOLOGY Vol.39,Sum.No.225
Jan.2021,No.1
脉动流化床内颗粒分离行为的数值模拟
张晓光
(中国特种设备检测研究院,北京　100029)
摘　要:脉动流化床可以实现多组分颗粒的分离而受到广泛关注。

基于此,这里采用双流体模型,对脉动流化床多组分混合颗粒分离行为开展了数值模拟研究。

通过与实验结果比较,模拟可以较好地预测颗粒的分离行为。

比较了脉动流和连续流对分离行为的影响,结果表明,脉动流可以极大地促进颗粒的分离效率。

同时评估了脉动频率和颗粒粒径的影响,研究发现,脉动频率的增大可以促进颗粒的分离效果,然而过高的频率会阻碍颗粒的分离。

关键词:脉动气流;流化床;分离;数值模拟;双流体模型
中图分类号:TK229 文献标识码:A 文章编号:1002-6339(2021)01-0050-04
Simulation of Segregation Behaviors of Binary Mixture in
a Fluidized Bed with Pulsed Flow
ZHANG Xiao -guang
(China Special Equipment Inspection and Research Institute,Beijing 100029,China)
Abstract :Pulsation -assisted fluidized beds can achieve the separation of multi -component particles,
which has attracted more and more concerns.Based on this,the two -fluid model is employed to numeri⁃cally investigate the segregation of binary mixture in pulsed fluidized bed.By a comparison with the ex⁃
perimental data,the model can give a good prediction.A comparison of pulsed flow and continuous flow is also conducted.The results reveal that the pulsed flow can greatly improve the segregation efficiency of particles.Meanwhile,the effects of pulsation frequency and particle size are evaluated.It is found that the increase of pulsation frequency can enhance the segregation degree.Whereas an excessive frequency
will hinder the particle segregation.
Key words :pulsed flow;fluidized bed;segregation;numerical simulation;two -fluid model 收稿日期　2020-02-19 修订稿日期　2020-07-15作者简介院张晓光(1986~),男,硕士研究生,高级工程师,主要
研究方向为锅炉设计与数值模拟。

流化床通过对不同物性颗粒的流化,可以实现
混合颗粒的分离,同时在流化分级过程中可以完成
干燥等过程,得到了研究者的广泛关注[1-3]。

袁竹
林等[4]采用离散颗粒方法模拟了流化床中不同密度颗粒分布。

研究发现,颗粒分离与表观气速和固
体流率有关,当风速较大和固体流率增大时,分离程度逐渐减小。

Lu 等[5]通过实验和理论研究了气体鼓泡流化床中不同粒径二元混合物流化行为,研究表明,颗粒大小,小颗粒质量分数和气速对二元混合体系分离有相当大的影响。

尹炜迪等[6]对选煤流化床内气固流动和颗粒分层进行了数值模拟。

研究
表明,静压梯度力占浮升作用的80%以上是主要作用力。

脉动气流的引入可以有效地调控流化床中的颗
·
05·
粒流化行为,当处理不同种类或不同性质的颗粒物料时,会涉及到双组分或多组分物料颗粒混合分离
等问题。

李占勇等[7]实验探究了不同脉宽比下脉动气流辅助流化对颗粒混合程度的影响。

结果表明,添加辅助脉动气流有助于提高颗粒混合程度,颗粒密度差别对双组分颗粒混合的影响大于颗粒直径差别对其的影响。

大多数脉动气流的引入是为了改善流化质量,提高颗粒的混合。

然而脉动气流同时也可以为颗粒分离提供一种有效的方法。

Saidi 等[8]通过实验的手段研究了气流脉动对流化床混合颗粒分离的影响。

实验表明,通过改变固定流速下脉冲气流的频率,混合二元粒子的分离效率显著提高。

基于双流体模型,结合多组分颗粒动理学理论,对脉动流化床双组分颗粒分离行为开展数值模拟,探讨了脉动气流对分离行为的影响,分析了脉动频率和颗粒粒径对分离效率的作用规律。

1　数学模型
为了研究双组分颗粒分离行为,双流体模型被采用。

假设颗粒是球形颗粒,连续性方程和动量守恒方程被表达如下
∂∂t (εg ρg )+∂
∂x i
(εg ρg u g )=0(1)∂
∂t
(εm ρm )+∂∂x i (εm ρm u mi )=0
(2)
∂∂t (εg ρg
u
gi )+
∂
∂x j (εg ρg u gj u gi
)
=-εg
∂P g
∂x j +εg ρg g i +∂τgij ∂x j
-∑n s
m =1
β
gm
(u gi -u mi )
(3)
∂∂t (εm ρm u mi )+∂∂x j
(εm ρm u mj u mi )
=-εm
∂P g ∂x j +εm ρm g i +∂τmij
∂x j +∑n
s
l =1,l ≠m
βlm (u li -u mi )+βgm (u gi -u mi )
(4)
为了描述颗粒的脉动能量的输运,颗粒拟温度被引入,颗粒拟温度输运方程表达如下
32εm ρm ∂Θm ∂t +U mj ∂Θm
∂x j
=∂∂x j κm ∂Θm ∂x j
+τmij
∂U mi ∂x j +Πm -J m (5)
为了实现模型的封闭,多组分颗粒动理学理论被采用[9]。

气固相间作用力和固固相间动量交换分
别采用Gidaspow 曳力模型
[10]
和Syamlal -O'B rien 模
型
[11]。

2　模拟对象
基于Saidi 等[8]实验,脉动流化床的三维模拟被
实施,反应器高度为0.5m,床径为0.11m,底部为速度入口,压力出口位于反应器顶部,压力为1atm。

对于壁面,无滑移边界条件被采用。

颗粒的密度为
2650kg /m 3,粒径为0.92mm 和0.33mm,属于D 类和B 类颗粒。

具体的操作条件和参数见表1。

表1　主要模拟参数
名称
数值反应器高度/m 0.5反应器管径/m 0.11颗粒密度/kg ·m -32650颗粒直径/mm 0.33/0.92初始颗粒高度/m 0.06初始颗粒份额1∶1进口速度/m ·s -10.16出口压力/Pa 1.013×105
脉动频率/Hz
1-9
模拟基于Ansys -Fluent 商业软件,通过自定义函数对边界条件进行修改以描述脉动流,这里采用
正弦波的脉动形式,其进口速度表达式如下
v =v 0+v 0sin(2πft )
(6)
计算时间为40s,时间步长为2×10-4s,通过
网格无关性验证,考虑计算精度和成本,14000个网格被采用。

3　结果与讨论
为了验证模型的可行性,图1给出了速度为
0.16m /s 时模拟得到的分离效率与实验值的比较。

图1　模拟结果与实验的比较
由图可见,模型能够对实验结果给出一个合理的预测。

尽管有一些差异,可能是由于在模拟中布风板的影响没有被考虑,同时在实验中分离层厚度的测量有一定的误差,总的来说,模型具有一定的可行性。

·
15·
图2给出了在连续流中,颗粒在不同时刻的颗粒浓度瞬时云图。

从图中可以看到,气泡的运动导
致颗粒床层的变化。

随着时间的推移,颗粒逐渐被分离,粗颗粒逐渐沉积在床层底部,细颗粒集中在床层上部。

图2　不同时刻下连续流的流化床颗粒浓度分布
图3进一步给出了在脉动流中,颗粒在不同时刻的颗粒浓度瞬时云图。

相比之下,连续流引起的分离在后期逐渐趋于稳定,随时间变化不是很显著。

图3　不同时刻下脉动流的流化床颗粒浓度分布(f =1Hz)
图4给出了气流停止进入,床层稳定后的颗粒浓度分布。

相较于连续流,脉动流的粗颗粒沉积层更厚,细颗粒和粗颗粒交界面颗粒浓度变化更加显著,而连续流,粗细颗粒分界面较为模糊,颗粒床层上部颗粒浓度依然较大,这说明脉动流有益于颗粒的分离。

图4　不同床型下颗粒浓度分布
图5给出了脉动流化床大小颗粒份额沿高度的变化。

从图中可以看到,大颗粒集中在床层底部
0.01m 以下,而随着高度进一步增加进入到了大小颗粒混合区,此时大颗粒浓度迅速减小,而小颗粒浓度剧增,在床层顶部,主要以小颗粒为主,尽管此时大颗粒仍有一些,总的来说颗粒分离效果较好。

图5　大小颗粒份额沿高度的变化(f =1Hz)
脉动频率是脉动流化床控制的重要因素。

图6给出了脉动频率对流化床分离效率的影响。

可以看到,随着脉动频率的增加,分离效率的变化并非线性的,而是呈现先增大后减小的趋势。

尽管脉动流会促进颗粒的分离,然而过高的脉动频率会加大布风板上方气体的扰动,进而使得多组分颗粒混合程度被促进,导致分离效率的降低。

图6　脉动频率对分离效率的影响
图7给出了小颗粒粒径的变化对双组份颗粒分离效率的影响。

可以看到随着小颗粒粒径的增大,分离效率逐渐减小,这主要是由于颗粒粒径差异变小,进而使两者的终端速度差异变小,颗粒的流动行为趋于相似而使分离效果变差。

图7　颗粒粒径对分离效率的影响
·
25·
4　结论
基于双流体模型,对脉动床双组分颗粒分离行为进行模拟研究,模拟结果与实验值吻合较好,比较了连续流和脉动流对流化床颗粒分离行为,同时分析了不同脉动频率对双组份颗粒分离的影响。

主要研究结论如下:
(1)相较于连续流,脉动流可以极大的提高颗粒的分离效率。

(2)颗粒的分离效率与脉动频率并非线性的,过高的脉动频率由于会增加床内的气流扰动造成分离效率下降。

(3)通过研究小颗粒粒径对分离效率的影响发现,颗粒的增大会降低颗粒的分离效果。

参考文献
[1]Wang G,Bai X,Wu C,Li W,Liu K,Kiani A,Recent Advances in the Beneficiation of Ultrafine Coal Particles[J].Fu⁃el Process Technol,2018(178):104-125.
[2]汪新智,王天宇,李春霆,等.鼓泡流化床内颗粒分离行为模拟研究[J].节能技术,2012(3):203-206. [3]FanX,Zhou C,Zhao Y,Duan C,Liu Q,Flow Pattern Transition and Coal Beneficiation in Gas Solid Fluidized Bed with Novel Secondary Distributor[J].Adv.Powder Technol,2018(29):1255-1264.
[4]袁竹林,徐益谦.用数值模拟对循环流化床颗粒分离特性的研究[J].工程热物理学报,2001(s1):185-188. [5]Lu H,He Y,Gidaspow D,Yang L,Size Segregation of Binary Mixture of Solids in Bubbling Fluidized Beds[J].Powder Technol,2003(134):86-97.
[6]尹炜迪,王庆功,吕俊复,等.选煤流化床内气固流动和颗粒分层的数值模拟[J].中国矿业大学学报,2019 (48):430-436.
[7]李占勇,潘波,高新源,等.脉动气流辅助流化下双组分颗粒的混合特性研究[J].农业机械学报,2015(46): 247-253.
[8]Saidi M,Basirat Tabrizi H,Chaichi S,et al.Pulsating Flow Effect on the Segregation of Binary Particles in a Gas-sol⁃id Fluidized Bed[J].Powder Technol,2014(264):570-576.
[9]Srivastava A,Sundaresan S.Analysis of a Frictional-kinetic Model for Gas-particle Flow[J].Powder Technol,2003 (129):72-85.
[10]Gidaspow D.Multiphase Flow and Fluidization:Con⁃tinuum and Kinetic Theory Descriptions[M].Academic Press, 1994.
[11]Syamlal M.The Particle-Particle Drag Term in a Multiparticle Model of Fluidization[R].Springfield:National Technical Information Service,1987.
(上接第43页)
参考文献
[1]王志轩,潘荔,刘志强,等.中国煤电清洁发展现状及展望[J].电力科技与环保,2018.34(2):1-8. [2]龙辉,黄晶晶.“十三五”燃煤发电设计技术发展方向分析[J].发电技术,2018.39(1):13-17.
[3]裴哲义,王新雷,董存,等.东北供热机组对新能源消纳的影响分析及热电解耦措施[J].电网技术,2017,41 (6):1786-1792.
[4]谢天,孙永春,付青山,等.200MW机组热电解耦方案研究[J].节能技术,2018,36(6):88-91.
[5]韩立,郭涛.350MW供热机组低压缸零出力经济运行研究[J].节能技术,2019,37(1):61-63,85. [6]李树明,刘青松,朱小东,等.350MW超临界热电联产机组灵活性改造分析[J].发电技术,2018,39(5):65-70.
[7]唐树芳,唐郭安,刘帅.300MW机组新型凝抽背供热改造技术分析[J].工程技术研究,2018,31(15):88-89.
[8]沈士一.汽轮机原理[M].北京:中国电力出版社.1996,164-175.
[9]A.T.Koxtkvt.汽轮机末级动叶栅在非设计工况下的振动[J].郑叔琛译.国外透平,1884,(24):58-61. [10]曹祖庆.汽轮机变工况特性[M].北京:水利水电出版社,1991,156-158.
[11]李晓勇,张明智,赵志芳.影响空冷机组背压及叶片的因素及参数选取方法[J].电力建设,2005,26(8):1-3. [12]王仲博.小容积流量下汽轮机末级长叶片可靠性的试验研究[J].中国电机工程学报,1987(4):45-51. [13]刘帅,高新勇,孙士恩,等.一种用于切除低压缸进汽的热电联产系统:201820288317.3[P].2018-03-01. [14]刘帅.丹东金山热电有限公司切除低压缸进汽改造可行性研究报告[R].杭州:华电电力科学研究院有限公司,2018,44-49.
[15]刘帅.丹东金山热电有限公司切除低压缸进汽工程总结报告[R].杭州:华电电力科学研究院有限公司, 2019,3-9.
·35·。