气升式环流反应器教程.ppt

气升式环流反应器内气液两相流动计算流体力学的模拟 姓名：承豪 专业：生物化工 学号：642081703002 目录 气升式环流反应器的几何结构 模拟结果及讨论 数值模拟过程 结论 气升式环流反应器区别于普通鼓泡床反应器的关键在于其内部有一个导流筒，气相从导流筒上升或从导流筒外的环隙上升，由此产生的导流筒内外的密度差驱动物料在反应器内进行循环，从而实现气液混合和必要的停留时间。 这个反应器具备了气升式反应器的基本特征—— 反应器外筒、导流筒以及气体分布器，反应器内充满常温液体水，气相从反应器底部的环形入口进入反应器，沿导流筒外的环隙空间(上升段)上升。反应器结构如图所示。 数值模拟过程 模拟过程参数的设定 欧拉一欧拉模型 欧拉模型方程组 欧拉一欧拉模型 目前使用较多的两流体模型是将每一种流体都看作是充满整个流场的连续介质，针对两相分别写出质量、动量和能量守恒方程，通过相界面间的相互作用(动量、能量和质量的交换)将两组方程耦合在一起。针对每一相可以用欧拉方法或拉格朗日方法描述，而采用基于两相均使用欧拉方法描述的便是欧拉一欧拉模型。该模型对每相求解一组运动方程和连续性方程，通过压力和相间交换系数进行耦合，耦合结果取决于所涉及的相的类型。当体系是气泡与液相(水)时，方程之间的耦合通过气泡和水之间的相互作用规律进行。 欧拉模型方程组 基于以上假设的欧拉模型方程组如下： 连续性方程： n为分散相的个数，本文研究的体系仅有一个分散相(气相)，所以，n=1；下标k分别代表气相(g)和液相(l)，本文不考虑气液之间的传质，所以 。 动量守恒方程： 模拟过程参数的设定 使用FLUENT6．0中的前处理软件Gambit用非结构化网格对反应器的物理模型进行了网格分割，基本尺寸为5 mm 时的网格划分图如下图。 算法：使用非耦合的求解器(segregated solver)进行非稳态计算，当模拟收敛且结果不再随时间变化时，说明已达到稳态，即通过非稳态计算得到稳态结果，对上述模型的模拟发现时间t=1180 S后，模拟结果随时间不再变化。 边界条件：气相进口为速度进口，反应器上部液面设为压力出口，反应器器壁及导流筒器壁均设为非滑移壁面。 初始条件：由于液相初始状态为静止状态，故其初始速度设为0，气相初始条件为环隙内的表观气速，其值分别为0.0191、0.0250、0.0315、0.0380、0.0440 m／s。 结果与讨论 气相表观速度对液相循环速度和气含率的影响 导流筒结构对气升式环流反应器内两相流动的影响 网格划分水平对模拟结果的影响 网格划分水平对模拟结果的影响 对反应器模型进行了由粗到细5种水平的网格划分，如下表所示。 由上图可知，网格水平4下的模拟结果与水平5下的结果很接近，即当网格划分达到水平4时，模拟结果不再随网格的细化而变化，说明这时候的模拟结果可看作网格独立的解。故后续模拟都是在网格水平4下进行的。 气相表观速度对液相循环速度和气含率的影响 气相作为液相循环的推动力来源，对液相循环速度和气含率会产生显著的影响，结果如图所示。 由图可知，环隙中液相表观速度开始随着气相入口速度的增大而增大，因为气量增大，导流筒内外的压力差亦增大，所以液相循环速度也增大。随着气量继续增大，液相循环速度增大趋缓，这是因为液相速度增大引起流动阻力增大，最终液相速度达到一个稳定值。 同时，导流筒内外(上升段)的气含率随着气相表观气速的增大而增大，液相循环速度的增大和上升段气含率增大，都会使气泡被液相夹带进入导流筒的量加大。根据导流筒外的气相表观速度UGR和导流筒外的液相表观速度ULR，分别由Hills和沈自求的关于上升段气含率εR与上升段表观气速和表观液速之间的关系，计算得到上升段气含率值。CFD 的模拟值和经验关系式的计算值均示于下图中。 由图可以看到，上升段气含率εR的CFD模拟值与沈自求关系式以及Hills关系式的计算值非常接近，也即操作表观气速 。 和CFD模拟得到的液相循环速度 和气含率很好地符合了这两种经验关系式，这验证了所使用的CFD模型具有很好的可靠性。 导流筒结构对气升式环流反应器内两相流动的影响 导流筒直径和导流筒位置高度对液相循环速度和气含率都会产生影响。 此图说明导流筒直径的增大并不能使液相循环量无限增大。 由图可以看出，导流简直径增大，上升段的气含率随着导流简直径的增大而增大。导流筒内的气含率很小，随着导流简直径的增大略有增大。 导流筒距离反应器底部位置升高，液相循环速增大。当然，导流筒顶端不能超过反应器液面，否则就无法形成有效的环流。 上升段气含率随着导流筒高度的增高而降低。导流筒内气含率接近于零，并无显著变化，说明该条件下，气泡被夹带进入导流筒的量很小。 (1)用两流体模型的模拟方法对气升式反