管式搅拌反应器流动特性与混合特性的CFD数值模拟

【摘要】： 管式反应器广泛应用于化工冶金领域,根据我国一水硬铝石矿的管道化溶出中效率不高、管壁结疤现象严重等特点,本实验室自主开发了一种带搅拌的叠管式反应器,目前已通过冷态物理模拟的实验方法对这种新型搅拌反应器的流动和混合特性进行了初步研究,但有关的理论及设计方法仍不完善。近年来,利用计算流体力学(CFD)的方法研究搅拌设备内的流动和混合特性逐渐发展起来。CFD技术以其卓越的优势,对搅拌设备的研究、开发与设计提供新的研究方法。 本文以自主开发的单、叠管式搅拌反应器为研究对象,利用商业CFD软件FLUENT6.3和并行计算系统的图形工作站,在已有的研究基础上对其进行了CFD的系统研究,并与实验数值进行对比,验证了利用CFD对管式搅拌反应器进行流动和传质特性研究的可行性,本研究能为反应器的进一步设计和改善提供新的研究思路,也可以作为实验研究和进一步模拟的基础。 首先,采用雷诺时均方程(RANS)的标准k-ε湍流模型和多重参考系(MRF)法,对单、叠管式搅拌反应器进行数值模拟,根据传统(无搅拌)管式反应器与新型带有搅拌装置的管式反应器流场的比较,并通过对管式搅拌反应器速度场、湍流动能场的分析,带搅拌的管式反应器起到了增加管内湍流区,增强混合效果等作用。此外本文还讨论了不同转速对管内速度分布、压力分布、湍流动能分布的影响。 在用RANS法对单、叠管式反应器的流场研究的基础上,采用速度场和浓度场分开求解的方法,计算了该管式反应器的停留时间分布,其模拟计算结果与实验测量结果吻合较好。通过对管式搅拌反应器内的停留时间分布模拟计算的结果显示,反应器内接近活塞流,带搅拌装置通过增加扰动、避免死区等来改善反应器的性能,还讨论了进口流量及搅拌转速对平均停留时间、方差及RTD曲线的影响。 本文采用标准k-ε湍流模型和多重参考系(MRF)法,对单管式搅拌反应器的混合时间及功率准数Np进行了数值模拟,其模拟值与实验值吻合较好,并考察了不同搅拌转速以及不同监测点的位置对混合时间的影响规律。 鉴于搅拌设备的复杂性,RANS法模型对于桨叶区湍流预测偏低,影响管内流体的物质交换的预报,本文采用大涡模拟(LES)的方法,用滑移网格(SM)法来处理搅拌桨区,对单管式搅拌反应器的流场进行数值模拟,并用该方法计算了停留时间分布和混合时间。结果表明LES法比标准k-ε模型更能有效解决该管式搅拌反应器之间的物质交换问题,尤其是在搅拌转速较高时模拟得到的平均停留时间及混合时间更趋于实验值,其中平均停留时间的平均相对误差降低了5%,混合时间模拟值的平均相对误差降低了13%。 在单管式搅拌反应器内采用玻璃珠-水体系对固-液两相流场进行数值模拟,从CFD角度研究反应器内固-液悬浮状况下速度分布、浓度分布,并考察了不同搅拌转速的影响规律,同时,在固体体积2%和7.5%两种体系下,采用不同搅拌转速反应器内固相体积分率的的最大值来判断固体颗粒临界悬浮转速Njs,其值约为35rpm。