各向异性k-ε湍流模型在Rushton桨搅拌槽三维流场整体数值模拟中的应用

根据搅拌槽内的流动呈各向异性的特点,引入适用于强旋转流场的各向异性k-ε湍流模型,用改进的内外迭代法对有挡板的Rushton桨搅拌槽进行了整体数值模拟.利用文献中对搅拌槽内流场测定结果,给出了适用于Rushton 桨搅拌槽的各向异性湍流黏度系数值.模拟计算得到了搅拌槽内的流场分布和脉动速度分布,并同标准k-e湍流模型计算结果及文献数据进行比较.结果表明,各向异性k-ε湍流模型能成功反映Reynolds应力、湍流动能等湍流特征量,明显优于标准k-ε湍流模型.     

各向异性k-ε湍流模型在Rushton桨搅拌槽三维流场整体数值模拟中的应用

根据搅拌槽内的流动呈各向异性的特点 ,引入适用于强旋转流场的各向异性k -ε湍流模型 ,用改进的内外迭代法对有挡板的Rushton桨搅拌槽进行了整体数值模拟 .利用文献中对搅拌槽内流场测定结果 ,给出了适用于Rushton桨搅拌槽的各向异性湍流黏度系数值 .模拟计算得到了搅拌槽内的流场分布和脉动速度分布 ,并同标准k -ε湍流模型计算结果及文献数据进行比较 .结果表明 ,各向异性k -ε湍流模型能成功反映Reynolds应力、湍流动能等湍流特征量 ,明显优于标准k -ε湍流模型 .     

大涡模拟搅拌槽中的液相流动

采采用大涡模拟湍流模型对有档板的Rushton 桨搅拌槽进行了数值模拟研究。控制方程采用控制容积法进行离散，对流项用三阶QUICK格式，扩散项是二阶中心差分。压力 速度耦合方程在交错网格上采用SIMPLE算法进行求解。小尺度流动的模化采用动力学（dynamic）亚格子模型。搅拌桨与档板之间的相互作用采用改进的内外迭代法进行处理。计算结果和文献值吻合得很好。     

两种不同挡板的双层桨搅拌槽三维流场数值模拟

采用CFD方法分别对三角形挡板和长方体空心挡板的双层斜直叶桨搅拌槽在水中产生的三维流场进行数值模拟,选取多重参考系法和标准k-ε湍流模型进行模拟,对两者的轴向平面速度场和时均速度分布进行了详细的对比分析,研究结果对搅拌槽的设计和实际应用具有重要的参考价值。     

桨式搅拌器功率准数三种取得方法的对比

在有挡板条件下,对常用的桨式搅拌器(单层二叶平桨、二叶斜桨、四叶斜桨及双层四叶斜桨),桨槽径比为0.5—0.6,进行搅拌功率曲线的测绘。利用经验公式对功率准数进行了计算,通过关联值与实验值的对比发现,Nagata关联式在层流状态时关联值与实验值相差较小,在湍流时二者相差较大,而Kamei和Hiraoka关联式则在过渡流和湍流区与实验值比较吻合,在层流区的偏差比较大。利用计算流体力学模拟了搅拌器各种状态的功率准数值,模拟值与实验值对比发现,模拟值在不同的雷诺数时都与实验值吻合较好。     

搅拌反应器中计算流体力学数值模拟的影响因素研究进展

回顾了搅拌反应器中计算流体力学（CFD）数值模拟方法的研究历程,包括搅拌桨叶边界条件法、内外迭代法、滑移网格法、多重参考系法、快速照相法等方法;结合各种模拟方法的特点,对机械搅拌反应器数值模拟过程中的挡板的存在、尾涡的消除、复杂桨叶类型的模拟、网格精度的影响、湍流模型的选择、离散格式的差异等几个因素的影响作了分析。提出了搅拌反应器的CFD数值模拟仍有待深入研究。     

搅拌釜中湍流流场和细长颗粒取向的模拟

采用两流体模型和改进的内外迭代法模拟了搅拌釜中固相为细长颗粒的液固两相湍流流场,采用标准的κ-ε模型封闭,用SIMPLE算法求解压力.模拟结果与搅拌釜中固相为球形颗粒时的流场进行比较,速度场基本相同,仅釜中心处速度略有差别.对细长颗粒取向模拟时,与以往文献多以Jeffery方程为研究出发点不同,本文作者直接利用刚体颗粒的运动和转动方程,得到搅拌釜中细长颗粒取向的变化规律.     

“波纹”内壁搅拌槽内的湍流流动特性

设计了一种“波纹”内壁搅拌槽,通过在内壁上安装一定数量的半圆管,以增强搅拌效果,并采用计算流体动力学方法对其内部流动特性进行了分析.通过与文献中的实验及模拟结果的比较,验证了所建数值模型及模拟方法的可靠性,随后对不同半圆管尺寸(半径)、数量和安装位置时搅拌槽内的流场结构和功率消耗进行了研究.结果表明,与标准搅拌槽相比,“波纹”内壁结构尽管不会明显增强流体的湍流程度,但能提高流体速度分布的均匀程度,而且功率消耗没有明显变化;当在搅拌槽内壁上均匀安装4根半径为搅拌槽直径1/20的半圆管时[LL1] ,能显著减少挡板后方及槽顶部流体的低速区.该结果为“波纹”内壁搅拌槽在过程工业中的应用奠定了基础.     

轴流桨搅拌槽三维流场数值模拟

利用k -ε湍流模型预测了搅拌槽在不同操作条件下宏观速度场 ,模型成功预测了搅拌槽内速度分布 ,计算结果与实验结果吻合较好 .模型预测结果表明 ,搅拌槽内宏观流动场受搅拌桨槽径比影响较大 .对单层搅拌桨 -槽体系 ,挡板前后宏观流动场差别很大 ,在挡板以前区域 ,轴向流动较强 ,在整个r -z断面上形成一个整体循环 ;而在挡板后面区域 ,流体在桨叶安装位置高度附近转向轴心流动 ,槽体上半部区域形成二次循环区域 ,且二次循环区域内流体以向下流动为主 .     

涡轮桨搅拌槽内单循环流动特性的大涡模拟

利用大涡模拟方法研究了涡轮桨搅拌槽内的单循环流动特性,采用Smagorinsky-Lilly动力亚格子模式,与文献实验及模拟数据进行了详细的比较. 结果表明,叶片后方的双尾涡偏向槽底运动,上尾涡在30o处已开始衰减. 800000个非均匀分布的计算网格和30个桨叶旋转周期的样本数据统计可获得准确的大涡模拟数据. 时均速度、均方根速度和湍流动能的大涡模拟值与实验数据一致,而k-e模型的模拟值与实验不符. 桨叶区呈现较强的各向异性,这是导致k-e模型预测不准确的主要原因. 对于搅拌槽内的复杂流动,大涡模拟方法是一个非常有效的工具.     

用大涡模拟结合改进的内外迭代法数值模拟搅拌槽中的液相流动

1 INTRODUCTION Stirred tank reactors are widely encountered in the chemical, pharmaceutical, and hydrometallurgical proc- esses. The fluid motion in stirred tanks is three-dimensional, complex, and covers a wide range of spatial and temporal scales. In the area surrounding the impeller, the flow is highly turbulent and swirling. The numerical simulation of such reactor systems is helpful in quantifying the effects of the impeller type, geometry, and the operational conditions in order to …     

用各向异性代数应力模型数值模拟搅拌槽中的三维全流场

In accordance to the anisotropic feature of turbulent flow, an anisotropic algebraic stress model is adopted to predict the turbulent flow field and turbulent characteristics generated by a Rushton disc turbine with the improved inner-outer iterative procedure. The predicted turbulent flow is compared with experimental data and the simulation by the standard κ-ε turbulence model. The anisotropic algebraic stress model is found to give better prediction than the standard κ-ε turbulence model. The predicted turbulent flow field is in accordance to experimental data and the trend of the turbulence intensity can be effectively reflected in the simulation. The distribution of turbulent shear rate in the stirred tanks was simulated with the established numerical procedure.     

搅拌釡中层流流场的模拟

Stirred tanks are used extensively in process industry and one of the most commonly used impellers in stirred tanks is the R.ushton disk turbine. Surprisingly few data are available regarding flow and mixing in stirred-tank reactors with Rushton turbine in the laminar regime, in particular the laminar flow in baffled tanks.In this paper, the laminar flow field in a baffled tank stirred by a standard R.ushton turbine is simulated with the improved inner-outer iterative method. The non-inertial coordinate system is used for the impeller region, which is in turn used as the boundary conditions for iteration. It is found that the simulation results are in good agreement with previous experiments. In addition, the flow number and impeller power number calculated from the simulated flow field are in satisfactory agreement with experimental data. This numerical method allows prediction of flow structure requiring no experimental data as the boundary conditions and has the potential of being used to scale-up and design of related process equipment.     

搅拌槽内液-固两相体系的数值研究(Ⅰ)流场的数值模拟

Three-dimensional solid-liquid flow is mathematically formulated by means of the “two-fluid” approach and the two-phase k-ε-Ap turbulence model. The turbulent fluctuation correlations appearing in the Reynolds time averaged governing equations are fully incorporated. The solid-liquid flow field and solid concentration distribution in baffled stirred tanks with a standard Rushton impeller are numerically simulated using an improved “inner-outer” iterative procedure. The flow pattern is identified via the velocity vector plots and a recirculation loop with higher solid concentration is observed in the central vicinity beneath the impeller. Comparison of the simulation with experimental data on the mean velocities and the turbulence quantities of the solid phase is made and quite reasonable agreement is obtained except for the impeller swept volume. The counterpart of liquid phase is presented as well. The predicted solid concentration distribution for three experimental cases with the average solid concentration up to 20% is also found to agree reasonably with the experimental results published in the literature.     

CFD Prediction of the Critical Agitation Speed for Complete Dispersion in Liquid‐Liquid Stirred Reactors

A computational fluid dynamics (CFD) model is adopted to simulate the turbulent immiscible liquid‐liquid flow in a stirred vessel based on a two‐fluid model with a k‐ϵ‐A_P turbulence model. An improved inner‐outer iterative procedure is adopted to deal with the impeller rotation in a fully baffled stirred tank. Different drag formulations are examined, and the effect of the droplet size on both the dispersed phase holdup distribution and the velocity field is analyzed. Two different numerical criteria are tested for determining the critical impeller speed for complete dispersion. The simulated critical impeller speeds are generally in good agreement with the correlations in the literature when the fixed droplet size is properly selected. This demonstrates that the modeling approach and the numerical criteria proposed in this work are promising for predicting the dispersion characteristics in liquid‐liquid stirred tanks.     

Numerical Simulation of Gas‐Liquid Flow in a Stirred Tank with Swirl Modification

Although the standard k‐? model is most frequently used for turbulence modeling, it often leads to poor results for strongly swirling flows involved in stirred tanks and other processing devices. In this work, a swirling number, R_S, is introduced to modify the standard k‐? model. A Eulerian‐Eulerian model is employed to describe the gas‐liquid, two‐phase flow in a baffled stirred tank with a Rushton impeller. The momentum and the continuity equations are discretized using the finite difference method and solved by the SIMPLE algorithm. The inner‐outer iterative algorithm is used to account for the interaction between the rotating impeller and the static baffles. The predictions, both with and without R_S corrections, are compared with the literature data, which illustrates that the swirling modification could improve the numerical simulation of gas‐liquid turbulent flow in stirred tanks.     

Numerical simulation of fluid dynamics in the stirred tank by the SSG Reynolds Stress Model

The Speziale, Sarkar and Gatski Reynolds Stress Model (SSG RSM) is utilized to simulate the fluid dynamics in a full baffled stirred tank with a Rushton turbine impeller. Four levels of grid resolutions are chosen to determine an optimised number of grids for further simulations. CFD model data in terms of the flow field, trailing vortex, and the power number are compared with published experimental results. The comparison shows that the global fluid dynamics throughout the stirred tank and the local characteristics of trailing vortices near the blade tips can be captured by the SSG RSM. The predicted mean velocity components in axial, radial and tangential direction are also in good agreement with experiment data. The power number predicted is quite close to the designed value, which demonstrates that this model can accurately calculate the power number in the stirred tank. Therefore, the simulation by using a combination of SSG RSM and MRF impeller rotational model can accurately model turbulent fluid flow in the stirred tank, and it offers an alternative method for design and optimisation of stirred tanks.     

Numerical simulation of a pitched‐blade turbine stirred tank with mirror fluid method

Mirror fluid method [Yang and Mao, Phys. Rev. 2005; E 71:036704] combined with local grid refinement is proposed to deal with the numerical simulation of turbulent flow in a pitched‐blade turbine stirred tank. By such a novel method, the domain occupied by the impeller is assigned suitable flow parameters explicitly by the mirror relation, so that the correct shear and normal forces on the fluid side of an interface segment is eventually guaranteed. Satisfactory agreement between our predictions and the reported experimental data is achieved both in single‐phase baffled or unbaffled stirred tanks and solid–liquid two‐phase systems. © 2012 Canadian Society for Chemical Engineering