纳米颗粒C语言数值模拟_细颗粒鼓泡流化床内气固流动特性数值模拟研究

流态化技术在燃煤、化工等领域已广泛应用。随着循环流化床(CFB)锅炉机组的参数和容量不断提高，300MW及以上容量的CFB锅炉普遍设置了外置或内置式流化床换热器，用于布置中温过热器和再热器，其工作状况直接影响整个机组的安全经济运行。这些外置或内置式流化床实质为细颗粒鼓泡流化床，其流化特性与主燃烧室内粗颗粒流化床的流化特性存在较大差异。

周旭和孙献斌等通过试验发现，在外置式换热器中，气固流动存在流化滞区。卢啸风等从机理上对贴壁流化滞区的形成进行了分析，流化滞区的存在会造成外置流化床内换热不均，使外置流化床中部受热面超温爆管。同时，张文清等通过试验发现，受热面的布置方式也会对流化床内部的气固流动状态产生影响。当受热面管束垂直于灰流动方向布置时，外置流化床边壁处颗粒边界层厚度较小；当受热面管束平行于灰流动方向布置时，外置流化床边壁处颗粒边界层厚度较大。

目前已经出现了多种外置式或内置式流化床换热器的结构方案，有水平流动方式、溢流流动方式、自上而下流动方式等。不同结构的外置式或内置式换热器，其内部的颗粒流动方式不同。因此，为了确保细颗粒鼓泡流化床内的流化均匀性，减小流化滞区，保证受热面的安全运行，需要对细颗粒鼓泡流化床中的气固流动特性进行深入研究。本文以自行搭建的新型双流化床冷模试验系统中的副流化床为研究对象，考虑不同的回流方式和受热面布置方式，对细颗粒鼓泡流化床流动特性进行三维数值模拟研究。

1

模型介绍

  1 .1  控制方程

对于气体相g，连续性方程为

对于固体颗粒相s，连续性方程为

对于气体相g，动量方程为

式中：p为气体压力；为重力加速度；b为相间动量传递系数；为气体相应力张量，其表达式为

式中：mf为流体有效黏度，kg/(m·s)；为单位张量。

对于固体颗粒相s，动量方程为

其中：ps为固相压力，用来表示由于颗粒间的相互作用而产生的法向力；为固相应力张量，其表达式为

颗粒拟温度θ定义为θ=Cv2/3，其中，Cv为颗粒波动速度。颗粒拟温度方程为

式中：dp为颗粒粒径，m；e为颗粒间的碰撞恢复系数，本文取0.95；颗粒与墙壁间的碰撞恢复系数取0.9。

使用标准k-e 湍流模型对每一相的湍动能k</