CFD+技术用于搅拌生化反应器液固两相混合的研讨.pdfVIP

J020 CFD 技术用于搅拌生化反应器液固两相混合的研究 洪厚胜，张庆文，万红贵，欧阳平凯 南京工业大学制药与生命科学学院，江苏 南京 210009 摘 要 本文应用 CFD 技术对机械搅拌生化反应器中液固两相混合进行了数值模拟。应用旋转参 考坐标系处理搅拌桨与挡板之间的相对运动，建立液固两相k−ε 双流体模型，采用有限体积法求 解流体流动控制方程组。数值计算结果与实验较为吻合，偏差在 10%以下。本文的研究结果对生 化反应器的优化设计及放大和混合的基础研究具有参考价值。 关键词： 计算流体力学；湍流；数学模型；分散度 中图分类号： 文献标识码： 1 引 言 在医药、食品、环境等领域存在着大量的液固非均相生化反应，其中固体物质常作为催 化剂、脱水剂、缩合剂等存在于反应液中[1] 。固体在液体中的分散好坏直接影响着反应的效 率。目前已有多种形式的反应器用于生化反应的流体混合，但机械搅拌式反应器由于其搅拌 器结构多样、操作范围宽，在生物技术实验室及工业领域始终占着主导地位[2] 。然而这并不 能说明机械搅拌式反应器就是一种混和性能优良的生化反应器，事实上，随着反应器规模的 增大，相间混合不均匀问题逐渐显现且愈演愈烈[3] 。要解决这一问题，必须从反应器设计入 手，但传统的宏观关联设计方法已显得无济于事。随着计算机技术、数理方法和流体力学的 进展，近年来，出现了一门新的交叉学科，即计算流体力学（Computational Fluid Dynamics， 简称 CFD ），CFD 技术为解决上述问题提供了新思路和新方法。本文通过建立液固两相湍流 双流体模型，对机械搅拌生化反应器中流体混合行为进行三维 CFD 数值计算。在模拟出局 部分散度并与文献实验结果对比的基础上，给出了详细的流场信息，便于对反应器中液固两 相混合行为与其结构和操作参数之间进行偶联分析，不仅对反应器优化设计具有经济意义， 而且对放大和混合的基础研究具有理论价值[4,5] 。 2 搅拌桨的 CFD处理 对于搅拌生化反应器中搅拌桨与挡板之间相对移动的动界面处理，在过去相当长的时间 [6] 内，一般采用将搅拌桨作为黑箱（black box ）处理的方法 ，即由实验测得搅拌桨周围虚构 表面的速度场作边界条件或将桨叶对流体的作用看作流体动量的产生源[7] 。这种方法的主要 缺陷是缺少有效的、稳定的边界条件，同时不能考虑到搅拌桨本身几何形状对流体的作用， 从而达不到优化搅拌桨和搅拌槽及操作参数的目的。因而必需要有一种新方法来加以解决。 近年来滑移网格思想得到广泛应用，整个流动区域被分成两部分，应用两不同的块网格 作相对滑移运动来模拟搅拌桨与挡板之间的相对运动[6] 。这是一种应用方便、行之有效的方 法，它不需要实验手段辅助。本文采用类似的处理方法——旋转参考坐标系的方法，流动区 - 1 - 域同样被划分成两部分，代表搅拌桨部分流体流动在旋转的参考坐标系下计算，其它部分在 固定的参考坐标系下计算，两者之间的交界面进行一定的速度修正。 收稿日期： ；修订日期： 作者简介：洪厚胜（1965 -），男，副教授，博士，通讯联系人，hhs@njut.edu.cn 。 - 2 - 3 数学模型 两相流计算模型主要有：（1）连续介质模型；（2 ）连续介质-轨道模型；（3 ）随机轨道 模型[8] [9] ；（4 ）双流体模型（颗粒相似流体模型） 。近年来，在两相流的数值模拟中，双流 体模型应用越来越广泛。本文采用双流体