多相反应器内流动和传递的介观与宏观尺度数值模拟.pdfVIP

多相反应器内流动和传递的介观与宏观 尺度数值模拟 杨超∗ 毛在砂 中国科学院过程工程研究所，绿色过程与工程重点实验室，北京 100190 摘要：多相搅拌槽等反应器是化工/冶金常用设备，本文介绍了多相反应器内流动和 传递的介观颗粒尺度和宏观设备尺度的数学模型和数值计算方法，给出对颗粒（液 滴、气泡）尺度、颗粒群尺度和宏观反应器设备尺度的多相流场、相含率分布、宏 观和微观混合状况、浓度场等数值模拟结果。利用多尺度数学模型化方法和数值模 拟，可以较准确地预测多相反应器内的流动、传递和反应特性，实现反应器的科学 设计放大和优化调控，并指导工业多相反应器和结晶器的优化改造。 关键词：多相反应器；多相流动和传递；模拟计算；多尺度 1. 引言 资源利用率低、能耗高、污染重已成为制约我国经济发展的瓶颈。在创建新的 过程工艺和改造传统产业时，不仅需要从工艺的角度注重原子经济性和环境友好， 更需要从工程共性的角度开发高效、节能和安全的反应分离设备。以化工领域最常 用的搅拌槽为例，搅拌桨和其它内构件产生的旋转剪切流场和多相传递过程的时空 结构复杂、影响因素多，具有典型的非均匀性、非平衡、非线性的多尺度特征，导 致工业设计和放大困难，大量工业反应设备的效率有待提升，研究反应器内多相流 体运动、混合、质量传递、能量传输和化学反应的复杂耦合过程，已成为近年的热 门课题。微观尺度（微混合、传质、剪切率、微粒碰撞等）、介观尺度（进料分布、 传热、颗粒和桨叶碰撞等）和宏观尺度（宏观混合、颗粒悬浮和分布等）的流体力 学和传递特性，直接决定反应器内的浓度（过饱和度）分布和反应、成核、生长、 团聚、破碎的速率，以及最终的产物分布、结构和纯度。 除了对多相体系进行系统、深入的实验研究和理论分析以外，越来越多的研究 人员也开展了基于数学模型的数值模拟。数值模拟可以部分地替代实验，而非线性 的多相流动、传递、反应过程耦合的复杂体系的理论分析往往十分困难，日益发展 的计算机和数值计算却能有效地求解涉及的偏微分方程模型，以认识多相体系的化 学工程性质。颗粒（包括液滴、气泡和固体颗粒）尺度的流动、传递和反应是在反 应器内的宏观多相过程形成的局部环境中进行，颗粒介观尺度上的性能在全反应器 的积分则是反应器的总体操作特性，多相化学反应器等过程设备的数学模型应该在 设备尺度和颗粒尺度上耦合进行，因此，颗粒尺度（包括单颗粒和多颗粒）的流动 感谢国家自然科学基金 (Nos. 、973 计划 (2004CB217604, 2007CB613507)和国家科技支撑计划(2008BAF33D03) 的经费资助。 *通讯联系人。Email: chaoyang@home.ipe.ac.cn; Tel: 010 和传递机理及规律是建立多尺度的多相反应器模型与数值模拟的基础。 在实验测量方面，PIV （Particle Image Velocimetry ）[1]等先进测量方法已被广 泛用于反应器内单相和低相含率的两相流场和相含率的测定，但对于高相含率或三 相体系的测量还很困难。在数学模型和数值模拟方面，利用Fluent 等商业软件平台、 采用基于雷诺平均（RANS ）的k–ε模型的单相或低相含率两相（气液或液固）的模 拟报道占近年研究的主体[2,3]。商业软件对单相体系的预测较好，但处理两相或三 相复杂体系还有待改进。大涡模拟[4]可获得优于RANS 模型的较准确瞬时湍流流场， 计算速度明显快于直接模拟法，是目前较有前途的反应器模拟方法。与反应器型式、 安装位置、加料、操作条件等密切相关的宏观混合和微观混合的模型和模拟一直是 多相研究和工业设计的核心问题之一[5]，但CFD （计算流体力学）耦合微观及宏观 混合模型、多相传递和反应模型的多尺度反应器模型和数值计算还有较大困难。研 究者已发展出了多种模拟颗粒尺度运动的数值方法，可以归为两大类，即拉格朗日 型和欧拉型，如标志网格法、移动网格法、前沿追踪法、边界积分法、离散粒子法、 连续对流法、流体体积法（VOF ）、水平集方法(level set method)等。本文将简要介 绍我们对反应器内多相流动和传质的介观颗粒尺度和宏观设备尺度的数值模拟方面 的研究进展。 2. 数学模型和数值计算方法 对于反应器内不可压缩流体的多相