搅拌管式反应器内流动特性的大涡模拟

采用实验和数值模拟的方法研究搅拌管式反应器内的混合过程,其中数值模拟采用大涡模拟的方法研究了反应器内流体的流动场,并就不同转速条件下流体的混合时间,将大涡模拟数值结果分别与标准k-ε模型的计算结果和实验测量值相比较,结果表明:管式搅拌反应器内的流动是非稳态的,具有不对称性。同时,大涡模拟方法可以预报漩涡,特别是桨叶背面的漩涡。与实验测量值相比,大涡模拟对混合时间的计算精度比标准k-ε模型计算精度高约22.8%,证明大涡模拟方法能够有效地模拟搅拌管式反应器内的流动特性。     

翼形桨搅拌槽内混合过程的数值模拟

采用FLUENT软件的多重参考系(MRF)及标准k-ε模型,将速度场与浓度场方程分开进行求解,对单层轴流式三叶CBY翼形桨搅拌槽内的混合过程进行了数值模拟,所得的混合时间的模拟结果与实验值相吻合。同时采用数值模拟的方法研究了不同的示踪剂加料点、监测点位置及操作条件对混合时间的影响规律;模拟结果表明,混合过程主要由搅拌槽内的流体流动所控制,混合时间与示踪剂加料点及监测点位置密切相关。上述的研究结果对于工业搅拌反应器的优化具有一定的参考意义。     

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

利用大涡模拟方法研究了涡轮桨搅拌槽内的流动特性,采用了三种亚格子模式:标准Smagorinsky-Lilly模式(SLM)、Smagorinsky-Lilly动力模式(DSLM)和亚格子动能动力模式(DKEM),并将模拟结果与标准k-ε模型及文献实验数据进行了详细的比较.结果表明:大涡模拟方法可获得搅拌槽内的瞬态流场;对桨叶区时均速度及湍流动能的预测与实验数据相吻合,比标准k-ε模型计算结果有明显改进,三种亚格子模型中DSLM和DKEM模拟结果更好.同时分析了大涡模拟中桨叶端部附近湍流动能估计偏差的原因,发现主要是由于对轴向湍流均方根速度的预测偏差造成的.大涡模拟方法为搅拌槽内非稳态、周期性的湍流流动和湍流特性的研究提供了强有力的工具.     

无挡板涡轮桨搅拌槽内湍流流动的分离涡模拟

采用分离涡模型对无挡板涡轮桨搅拌槽内的湍流流动进行了研究,重点分析了流场结构和速度分布,以检验该模型模拟搅拌槽内流体流动的有效性和正确性. 为了加快收敛,先采用标准k-e模型进行稳态流场计算,并以此结果为初始值进行分离涡模拟. 与现有文献大涡模拟及实验结果对比表明,分离涡模型能捕捉槽内流体的瞬时流动特征,获得的时均速度分布与大涡模拟及实验结果吻合较好,其中对切向速度分布的预测误差不超过7%,对径向速度分布的预测精度则低一些,局部误差接近12%. 分离涡模型适用于无挡板涡轮桨搅拌槽内湍流流动的模拟,能获得与大涡模拟相近的结果,且计算量更小(约为大涡模拟的1/3).     

多叶片组合式搅拌桨釜内流动特性和混合性能研究

开发可适用于较宽黏度范围的搅拌桨,强化釜内的流体流动和混合过程对于搅拌釜的节能增效具有重要的意义。实验与数值模拟相结合,在大涡模拟层面研究了多叶片组合式搅拌桨(MBC桨)从层流到湍流状态下,釜内的功率特性、流场分布、湍流特性和混合性能。结果表明:预测的功率曲线与实验结果一致;层流状态下釜内以切向流动为主,随着Reynolds数(Re)的增大,釜内轴向和径向流动逐渐增强,当Re达到486时,速度场分布与湍流状态下基本一致;在相同的能耗水平下,MBC桨对高黏度流体的混合性能优于商业Maxblend桨。桨叶的分散组合布置,强化了釜内的轴向和径向流动,使得MBC搅拌桨在从过渡流到湍流状态下均可实现较大的轴径向流动,湍动能分布较为均匀,混合过程显著加快。     

改进型INTER-MIG桨种分槽内混合过程的实验研究及大涡模拟

采用褪色实验法和数值模拟相结合的方法对双层改进型INTER-MIG桨种分槽内的混合过程进行研究,用高速相机记录实验中种分槽内的混合过程,结合大涡模型(LES)及动态Smagorinsky–Lilly模式亚格子模型求解湍流流动及示踪剂传递过程.结果表明,在近液面处加料时LES预测的混合过程与实验吻合,示踪剂呈螺旋状扩散,但预测的混合时间偏大;最佳加料点位于下层桨叶附近区域,其混合效率比在近液面处加料高17.16%;槽体中部区域监测到的混合时间最小,并分别向槽顶和槽底方向增大;改进型INTER-MIG桨种分槽槽底区域是混合困难区域.     

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

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

基于FLUENT软件的管式搅拌反应器流场的数值模拟

运用商业计算流体力学(CFD)软件FLUENT对一种新型管式搅拌反应器进行流场模拟,用GAMBIT建立流场实体模型,采用标准k-ε湍流模型以及多重参考系法(MRF)处理搅拌桨区.结果表明,计算所选模型能较准确地预测搅拌反应器的速度场、压力场及湍流动能分布,考察流量为1 m3/h时不同搅拌转速对反应器内部混合的影响,此时对应的最佳搅拌转速为50 rpm左右.模拟结果将为实验研究提供适当的操作参数,对搅拌反应器的优化和放大具有一定参考价值.     

煤浆混合槽混合性能的数值模拟

利用计算流体力学(CFD)方法研究煤浆混合槽的混合性能。采用Fluent 6.2计算软件,选用多重参考系法(MRF)、RNGk-ε湍流模型以及混合物多相流模型,对卧式煤浆混合槽的内部流动情况进行数值模拟,得到了湍流强度、湍动能、压力、速度分布以及出口各相分布等混合性能。模拟结果与实验结果进行了初步比较,吻合较好。研究结果表明,混合时在挡板和搅拌浆叶之间产生漩涡,有利于防止物料堆积,使混合更均匀。     

搅拌槽内多相流动数值模拟研究进展

综述了过程工业中广泛应用的搅拌槽反应器内多相流动数值模拟研究的进展.讨论多相湍流模型、相间作用力模型及搅拌桨处理方法等重要的数值模拟技术和方法,并对有关的计算模型进行了比较分析.针对搅拌槽内的各种多相体系,论述了不同研究者在桨区处理、相间动量传递描述和分散相的引入对体系湍流特性影响等方面的模拟方法并对结果作了比较.提出了需要进一步深入研究的课题.     

Large Eddy Simulations of Mixing Time in a Stirred Tank

Large eddy simulations （LES） of mixing process in a stirred tank of 0.476m diameter with a 3-narrow blade hydrofoil CBY impeller were reported. The turbulent flow field and mixing time were calculated using LES with Smagorinsky-Lilly subgrid scale model. The impeller rotation was modeled using the sliding mesh technique. Better agreement of power demand and mixing time was obtained between the experimental and the LES prediction than that by the traditional Reynolds-averaged Navier-Stokes （RANS） approach. The curve of tracer response predicted by LES was in good agreement with the experimental. The results show that LES is a reliable tool to investigate the unsteady and quasi-periodic behavior of the turbulent flow in stirred tanks.     

STUDY OF THE TURBULENT FLOW IN AN UNBAFFLED STIRRED TANK BY DETACHED EDDY SIMULATION

Detached eddy simulation (DES) of the liquid-phase turbulent flow in an unbaffled stirred tank agitated by a six-blade, 45°-pitched blade turbine was performed in this study. The tank wall is cylindrical with no baffle and the fluid flow problem was solved in a single reference frame (SRF) rotating with the impeller. For the purpose of comparison, computation based on large eddy simulation (LES) was also carried out. The commercial code Fluent was used for all simulations. Predictions of the phase-averaged turbulent flow quantities and power consumption were conducted. Results obtained by DES were compared with experimental laser Doppler velocimetry (LDV) data from the literature and with the predictions obtained by LES. It was found that numerical results of mean velocity and turbulent kinetic energy profiles as well as the power consumption are in good agreement with the LDV data. When performed on the same computational grid, which is under-resolved in the sense of LES, DES allows better accuracy than LES in that it works better in the boundary layers on the surface of the impeller and the stirred tank walls. It can be concluded that DES has the potential to predict accurately the turbulent flow in stirred tanks and can be used as an effective tool to study the hydrodynamics in stirred tanks.     

Stereo‐PIV experiments and large eddy simulations of flow fields in stirred tanks with Rushton and curved‐Blade turbines

The flow characteristics in pilot‐scale stirred tanks with Rushton and curved‐blade turbines were investigated by using stereoscopic particle image velocimetry (SPIV) experiments and large eddy simulation (LES) methods. The velocity and turbulent kinetic energy (TKE) in the impeller discharge regions were carefully resolved with a high resolution SPIV system, and the detailed phase‐resolved velocity and TKE profiles were used to validate the LES results. The effects of Reynolds number and blade shape on the flow characteristics were discussed. The LES results of velocity, TKE, and the evolution of trailing vortices were compared with the SPIV experimental data, and good agreement was obtained at various phase angles. The effects of subgrid scale model and hybrid grid with different mesh resolutions on the LES results were investigated. LES is a computationally affordable method for the accurate predictions of the complex flow fields in pilot‐scale stirred tanks is presented. © 2013 American Institute of Chemical Engineers AIChE J, 59: 3986–4003, 2013     

涡轮桨搅拌槽内湍流特性的V3V实验及大涡模拟

分别用体三维速度测量技术（volumetric three-component velocimetry measurements,V3V）和大涡模拟（large eddy simulation,LES）方法对涡轮桨搅拌槽内流场进行研究,发现在完全湍流状态下,涡轮桨搅拌槽内流场的量纲1相平均速度及湍动能分布同Reynolds数无关。用V3V方法实现了Rushton桨叶附近三维流场的重构;探讨尾涡的三维结构及运动规律;分析了叶片后方30°截面轴向、径向和环向速度沿径向分布规律。用V3V实验结果对比了2D-PIV（particle image velocimetry）数据中的尾涡涡对位置和涡量,涡对位置吻合度较好,但2D-PIV中涡量较V3V小37.5%;通过大涡模拟得到完整的尾涡结构,发现在叶片上边缘后侧存在一个和尾涡形成方式相同但不成对出现的涡结构;将大涡模拟结果和2D-PIV及V3V实验结果对比发现,大涡模拟在速度分布及尾涡运动轨迹方面均同实验结果吻合较好,表明大涡模拟能较好地预测涡轮桨搅拌槽内流场。     

PIV experiments and large eddy simulations of single-loop flow fields in Rushton turbine stirred tanks

The single-loop flow fields in Rushton turbine stirred tanks with clearance C=0.15T (T is tank diameter) were investigated by using particle image velocimetry (PIV) experiments and large eddy simulation (LES) methods. The velocity and turbulent kinetic energy (TKE) were carefully measured and resolved with high resolution camera. The regions with high TKE are affected by the movement of the trailing vortices generated behind the impeller blades. The effects of both geometrical configuration and Reynolds number were discussed. It is found that the Reynolds number has little effect on the mean flow for the configuration of impeller diameter D=T/3, C=0.15T. However, the single-loop flow pattern is changed into a double-loop one if D is increased from T/3 to T/2. The LES results were compared with the PIV experiments and the laser Doppler anemometry (LDA) data in the literature. The effect of the grid was validated, and the levels of local anisotropy of turbulence near the impeller discharge regions were investigated. Both the phase-averaged and phase-resolved LES results are in good agreement with the PIV experimental data, and are better than the predictions of the k–ε model. The agreement shows that the LES method can be used to simulate the complex flow fields in stirred tanks.     

搅拌釜三维流场的雷诺应力模型数值模拟

引言 近年来,雷诺平均Navier-Stokes方程及其封闭方程组成的湍流模型在应用于搅拌釜复杂湍流现象的模拟上取得了成功,主要是采用标准k-ε[1]、RNG k-ε[2]等模型.     

多层翼形桨搅拌槽内混合过程的CFD模拟

The mixing process in a stirred tank of 0.476m diameter with single, dual and triple 3-narrow blade hydrofoil CBY impellers was numerically simulated by using computational fluid dynamics （CFD） package FLU-ENT6.1. The multi-reference frame （MRF） and standard k-ε turbulent model were used in the simulation. The shaft power and the mixing time predicted by CFD were in good agreement with the experiment. The effects of tracer feeding and detecting positions on mixing time were investigated. The results are of importance to the optimum design of industrial stirred tank/reactors.     

CFD simulation of liquid-phase mixing in solid-liquid stirred reactor

A comprehensive CFD model was developed to gain an insight into solid suspension and its implications on the liquid-phase mixing process in a solid-liquid stirred reactor. The turbulent solid-liquid flow in a stirred reactor was simulated using a two-fluid model with the standard k-ε turbulence model with mixture properties. The multiple reference frames (MRFs) approach was used to simulate impeller rotation in a fully baffled reactor. The computational model with necessary sub-models was mapped on to a commercial solver FLUENT 6.2 (of Fluent Inc., USA). The predicted solid concentration distribution was compared with the experimental data of Yamazaki et al. [1986. Concentration profiles of solids suspended in a stirred tank. Powder Technology 48, 205-216]. The computational model was then further extended to simulate and understand the implications of the suspension quality on liquid-phase mixing process. The computational model and the predicted results discussed here will be useful for understanding the liquid-phase mixing process in stirred slurry reactors in various stages of solid suspension.