节能和环保是我国的两个重大问题,循环流化床燃烧技术因其较高的环保性能,受到了广泛的关注。但由于炉内气固两相流动的复杂性,物料的多次循环造成了炉内各部位不同程度的磨损。因此,研究和预防炉内磨损已经成为一个迫切需要解决的问题。由于循环流化床内流动十分复杂以及实验器材的制约,传统实验方法在循环流化床研究中存在很多的局限性。随着计算机技术的迅猛发展,越来越多的研究采用数值模拟的方法。本文采用欧拉-拉格朗日坐标系下的DEM软球模型在二维空间对循环流化床内气固两相流动进行了数值模拟,其中颗粒相采用DEM(离散单元法),气相采用SIMPLER算法,气固两相之间采用基于牛顿第三定律的耦合模型。并用FORTRAN语言编制了相应的计算程序。采用编制的程序首先对单孔射流流化床的流动特性进行了数值模拟,假设不考虑传热、燃烧及颗粒转动,得到了不同时刻气固两相的流动情况,其中包括床高随入口空气流速的变化情况,颗粒的混合情况、空隙率的分布、气体相和固体相在空间中的速度分布、气固两相的流动特性等等。结果表明,随着入口空气流速的提高,流化床床层也会相应的升高。流化床中部的颗粒向上运动并向两边壁面方向扩散,在两边壁面处的颗粒向下运动,底层的颗粒受到喷动气流的卷吸作用而向中间汇聚。颗粒在空间的分布表现为上部为浓度较低的稀相区和下部为浓度较大的密相区。颗粒间相互碰撞、掺混,体现出循环流化床内气固两相流复杂而剧烈的运动。随后,根据数值模拟的结果,分析了循环流化床内颗粒浓度的分布规律以及颗粒碰撞壁面的位置、次数和速度,计算了壁面不同高度的磨损量。分析了冷态时循环流化床内的磨损特性,并提出了减少磨损的建议。结果表明,在循环流化床空间中,底部密相区的碰撞频率要高于稀相区,而稀相区中两边壁面的碰撞频率要高于中间区域。在密相区的颗粒对壁面的碰撞次数最多,磨损量也最大。本文的研究对今后运用数值模拟的方法进一步研究循环流化床内的传热及燃烧过程奠定了基础。