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.     

FlowVision在轮胎噪声仿真中的应用

Based on FlowVision simulation software,the fluid-structure interaction method was used to numerically simulate the near-field noise of a tire.The flow field visualization results showed the motion of the flow field around the rolling tire.Comparing the near-field noise simulation with the experimental sound pressure level results,it was found that the experimental and simulation results had obvious directivity.The error of peak frequency between the simulation results and the test results in the 90° direction was only 4～8 Hz,and the test,simulation and theoretical calculation results had good consistency.     

A CFD model for predicting the heat transfer in the industrial scale packed bed

Compared to the traditional lumped-parameter model,computational fluid dynamics (CFD) attracted more attentions due to facilitating more accurate reactor design and optimization methods when analyzing the heat transfer in the industrial packed bed.Here,a model was developed based on the CFD theory,in which the heterogeneous fluid flow was resolved by considering the oscillatory behavior of voidage and the effective fluid viscosity.The energy transports in packed bed were calculated by the convection and diffusion incorporated with gaseous dispersion in fluid and the contacting thermal conductivity of packed particles in solids.The heat transfer coefficient between fluid and wall was evaluated by considering the turbulence due to the packed particles adjacent to the wall.Thus,the heat transfer in packed bed can be predicted without using any adjustable semi-empirical effective thermal conductivity coefficient.The experimental results from the literature were employed to validate this model.     

Gas-liquid hydrodynamics in a vessel stirred by dual dislocated-blade Rushton impellers☆

Towards the objective of improving the gas dispersion performance, the dislocated-blade Rushton impeller was applied to the gas–liquid mixing in a baffled stirred vessel. The flow field, gas hold-up, dissolved oxygen, power consumption before and after gassing were studied using the computational fluid dynamics (CFD) technique. Dispersion of gas in the liquid was modelled using the Eulerian–Eulerian approach along with the dispersed k–εturbulent model. Rotation of the impeller was simulated with the multiple reference frame method. A modified drag coefficient which includes the effect of turbulence was used to account for the momentum exchange. The predictions were compared with their counterparts of the standard Rushton impeller and were validated with the experimental results. It is concluded that the dislocated-blade Rushton impeller is superior to the standard Rushton impeller in the gas–liquid mixing operation, and the findings obtained here lay the basis of its application in process industries.     

中等Reynolds数线性剪切泡状流的液相雷诺应力

Analytical investigation of liquid Reynolds stress in shear bubbly flow with intermediate Reynolds numbers is absent. In this paper, the velocity field around a sphere bubble in linear shear liquid is assumed to be the linear superposition of the velocity field of uniform liquid passing a sphere bubble and the linear shear velocity field. The formula of shear liquid Reynolds stress was derived by averaging the velocity field in the cell-ensemble averaging method, and the formula was corrected under conditions of intermediate Reynolds number. The formula was compared with that of Sato, and the predicted results of local liquid velocity of the fully developed upward bubbly flow in pipes were compared with the experimental data. The results show that the formula is valid and accurate in prediction.     

搅拌釜内固-液悬浮体系的CFD-DEM模拟(英文)

Computational fluid dynamics-discrete element method （CFD-DEM） coupled approach was employed to simulate the solid suspension behavior in a Rushton stirred tank with consideration of transitional and rotational motions of millions of particles with complex interactions with liquid and the rotating impeller. The simulations were satisfactorily validated with experimental data in literature in terms of measured particle velocities in the tank. Influences of operating conditions and physical properties of particles （i.e., particle diameter and density） on the two-phase flow field in the stirred tank involving particle distribution, particle velocity and vortex were studied. The wide distribution of particle angular velocity ranging from 0 to 105 r·min-1 is revealed. The Magnus force is comparable to the drag force during the particle movement in the tank. The strong particle rotation will generate extra shear force on the particles so that the particle morphology may be affected, especially in the bio-/polymer-product related processes. It can be concluded that the CFD-DEM coupled approach provides a theoretical way to under-stand the physics of particle movement in micro-to macro-scales in the solid suspension of a stirred tank.     

搅拌釡中层流流场的模拟

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.     

A Quantitative Structure Property Relationship for Prediction of Flash Point of Alkanes Using Molecular Connectivity Indices

Many structure-property/activity studies use graph theoretical indices, which are based on the topological properties of a molecule viewed as a graph. Since topological indices can be derived directly from the molecular structure without any experimental effort, they provide a simple and straightforward method for property prediction. In this work the flash point of alkanes was modeled by a set of molecular connectivity indices （Х）, modified molecular connectivity indices （ ^mХ^v ） and valance molecular connectivity indices （ ^mХ^v ）, with ^mХ^v calculated using the hydrogen perturbation. A stepwise Multiple Linear Regression （MLR） method was used to select the best indices. The predicted flash points are in good agreement with the experimental data, with the average absolute deviation 4.3 K.     

Interfacial Shear Stress of Stratified Flow in a Horizontal Pipe

Experimented data are presented for the void fraction aud the shear stresses of stratified gas-liquid flow in a pipe, A new technique was used to measure the interface shear strew. The interfacial shear stress was determined by using two methods: a momentum balance of gas and an extrapolation of the Reynolds shear stress prone at the gas-liquld interface. A new formula ,relatiog to the interfacial friction factor with the void fraction and superficiol gas Rcynold number, was dewloped to predict the interface shear stress. The predicted values are in good agreement with experimental data.     

HEAT TRANSFER STUDIES OF HIGHLY VISCOUS NON-NEWTONIAN FLUIDS IN VERTICAL TUBES BY FINITE ELEMENT METHOD

A numerical method capable is developed for handling steady laminar flow and heat trans-fer of a highly viscous power-law fluid whose density,viscosity,specific heat and thermalconductivity,vary with temperature.The governing equations are found to be continuity,monmentumand energy expressions.Important effects such as varying viscosity,natural convection and viscousdissipation are incorporated in the theoretical model.These equations are being attracted by employing a decoupled finite element method.Galerkin'sprinciple is used to handle the momentum and continuity equations.Consistent(SU/PG)andnon-consistent(SU)streamline upwind methods are employed for the energy equation.Comparisonof calculated results and experimental data shows good agreement.Similar results are obtained withSU and SU/PG methods.Velocity and temperature profiles which provide insights into the processare also given.     

最大叶片式桨流体动力学性能的数值模拟

运用CFD方法研究了最大叶片式桨在层流区域内的流体动力学性能,模拟体系为高粘度牛顿流体和非牛顿流体,主要考察桨叶的功耗特性、Metzner常数、剪切性能和排液性能. 结果表明,功耗计算值与文献实验值基本一致;桨叶的Metzner常数ks=10与流变行为指数n无关;搅拌形成双循环流型,釜中部桨叶所在区剪切速率大、排液量大,产生的漩涡流也大,导致剪切效率低于0.5;随雷诺数增加(Ren=1.4, 5.0, 7.6),全釜平均剪切速率(3.40, 9.91, 15.05 s-1)和全釜平均排液量(0.0014, 0.0033, 0.0052 m3/s)逐渐增加,尤其是桨叶下端两翼区平均剪切速率(4.36, 11.48, 16.35 s-1)和平均排液量(0.0026, 0.0064, 0.0095 m3/s)增加相对较大. 说明Ren增加,搅拌混合作用加强,剪切速率大产生的界面积大,排液量大使高低剪切区内流体快速循环,有利于流体高效混合.     

Numerical simulation of three-dimensional flow field of spiral ribbon-screw impeller

应用计算流体力学数值模拟方法,研究了螺带-螺杆组合式搅拌器对高黏度非牛顿流体的搅拌流场特性。结合实例,考察了宏观流场、速度场、剪切速率和表观黏度的分布,分析了流变指数对搅拌特性的影响,并计算了不同工况下的搅拌功率。结果表明：螺带-螺杆组合搅拌器针对高黏度流体,具有良好的全局搅拌特性,在本文的模拟条件下,组合桨的Metzner常数k_s约为13.8,而流变指数变化对k_s的影响不大,为螺带-螺杆搅拌桨在高黏度非牛顿流体搅拌中的应用提供有益指导和参考。     

锚式搅拌槽中高粘牛顿流体的流速分布及其粘性耗散

本文用照相法测定了锚式搅拌槽中高粘度流体的流型和速度分布.利用最小二乘法获得速度分布拟合方程以及槽中剪切率γ和耗散能Φ分布方程.比较了槽壁区域剪切率和全槽平均剪切率、Metzner常数k_S和搅拌器特性参数C_3,指出k_S是表征槽中高剪切区域剪切率大小的特征量,C_3表征了全槽平均剪切率的大小.对于均匀剪切场,两者是一致的.可以用单位体积剪切率分布系数φ_V来描述槽中剪切率分布的均匀程度.本文还研究了Re数对速度分布流型的影响,发现在低Re数下,流型是对称的.高Re数下,流型则不对称.     

螺带型叶轮搅拌槽内流场分析与计算

采用粒子图像测速与数值计算方法研究了螺带桨流场,选用11种几何结构桨型及5种流变性质的流体(0.3≤n≤1),考察了层流流动(Re<100)状态下,桨叶结构、转速及流体假塑性对流型、速度分布及剪切速率分布的影响. 计算表明,核心区无量纲旋转角速度(a)随桨叶直径及螺距增大而减小,其值约在0.48~0.89;无量纲切向、轴向及径向速度分别约为0.44~0.47, 0.06~0.09及0.04,且前两者随螺距的增大分别呈线性减小及增大. a随流体粘度降低及假塑性增强而减小,非牛顿流体的流速小于牛顿流体;剪切速率随桨叶直径增大及流体假塑性增强而增大.     

Using electrical resistance tomography and computational fluid dynamics modeling to study the formation of cavern in the mixing of pseudoplastic fluids possessing yield stress

Electrical resistance tomography (ERT) provides a non-intrusive technique to examine, in three dimensions, the homogeneity and flow pattern inside the mixing tank. In this study, a 4-plane 16-sensor ring ERT system was employed to study the shape and the size of cavern generated around a radial-flow Scaba 6SRGT impeller in the mixing of xanthan gum solution, which is a pseudoplastic fluid possessing yield stress. The size of cavern measured using ERT was in good agreement with that calculated using Elson's model (cylindrical model). The 3D flow field generated by the impeller in the agitation of xanthan gum was also simulated using the commercial computational fluid dynamics (CFD) package (Fluent). The CFD model provided useful information regarding the impeller pumping capacity, flow pattern, and the formation of cavern around the impeller. CFD results showed good agreement with the experimental data and theory.     

Effect of impeller type and scale-up on spatial distribution of shear rate in a stirred tank

The main spatial distribution features of shear rate in a stirred tank operated with five different radial and axial flow impellers were presented with particle image velocimetry (PIV) experiments. Not only the average shear rate in the whole tank but also the local value in the vicinity of impeller increases linearly with impeller speed. Furthermore, the shear coefficient (K_s,imp) at the impeller outlet is linearly related to the impeller flow number (N_q) and decreases with the increase of N_q in general at the constant power consumption per unit volume (P_v). During scale-up based on the constant P_v and geometric similarity, CFD results show that the volume-averaged shear rate (γ_avg) for RDT decreases faster than that of other impellers with the impeller tip velocity (U_tip). The novel multi-blade combined (MBC) impeller with the increased height-to-diameter ratio of the stirred tank is able to more effectively improve the distribution uniformity of shear rate at the same P_v after scale-up. These studies provide a data basis for selecting the impeller types and improving the shear rate environment in the large-scale stirred tank.     

Analysis of shear-induced coagulation in an emulsion polymerisation reactor using computational fluid dynamics

Shear-dependent coagulation is a costly problem for the latex manufacturing industry, due to product degradation and reactor downtime. In this study, a method for calculating the shear-dependent coagulation rate in emulsion polymerisation is developed. The method combines simple models for coagulation (only binary collisions being considered) with the effects of rheology on the flow field, using computational fluid dynamics (CFD) to solve the detailed flow field in the reaction vessel. By using the local shear rates (LSR), the method developed provides a more detailed and system-specific assessment compared with using an average shear rate (ASR) for calculating the coagulation rate. The difference in the predictions between the ASR and the proposed LSR method was investigated. It was found that the ASR and LSR methods predict different coagulation rates, especially for more sophisticated coagulation models where the coagulation rate is not linearly dependent on the shear rate. The LSR method was also used to study the effect of the rheology of the latex, of the impeller speed and of the reactor design on the coagulation rate. It was found that the LSR method is useful for providing both visual and numerical means to identify regions with elevated coagulation rates in the modelled reaction vessel. The treatment provides estimates of the amounts of coagulum formed on the vessel walls and on the impeller.     

最大叶片式桨在假塑性流体中的搅拌流场模拟

为研究最大叶片式桨在高黏假塑性流体中的搅拌流动行为,以黄原胶溶液为研究体系,采用计算流体力学方法重点研究了釜内流体的功耗特性、速率分布、剪切速率、表观黏度分布和总体流动状况。结果表明:最大叶片式桨具有与大多数径流桨相似的"双循环"流型结构,且预测的功耗特性与实验数据一致性良好。最大叶片式桨适用于高黏假塑性流体的混合,而对于高黏牛顿流体的混合则效果不佳。釜内的剪切速率分布较宽泛,且受转速影响较大。转速可作为该桨改善黄原胶体系混合效率的重要参数之一。     

Cavern Boundary Determination with Pseudoplastic Fluid Based on the Apparent Viscosity Method

The cavern characteristics of xanthan gum solutions stirred by impellers of a perturbed six-bent-bladed turbine (6PBT) were investigated numerically using computational fluid dynamics (CFD) with laminar flow model. The apparent viscosity method was proposed to determine the boundary of the cavern. The cavern sizes predicted by apparent viscosity were in good agreement with the calculated results, and this method was not influenced by speed and impeller configuration. However, the predicted results of the cavern by the traditional speed method usually displayed a great deviation, especially at high speed. Therefore, it is feasible to determine the cavern boundary with the apparent viscosity, i.e., 0.25 times the yield viscosity of a pseudoplastic fluid, as the unified standard.     

内外单螺带式搅拌器的Metzner常数

对１１种不同向何尺寸的内外单累带搅拌器（ＩＯ－ＨＲ），在６６中不同流动行为指数（ｎ＝０．３７－１）的幂律流体中的剪切特性进行了实验研究。发现内外单螺带 的Ｍｅｔｚｎｅｒ常数，Ｋｓ与流动行为指数ｎ、桨的几何尺寸均相关。提出改进的数据处理方法来求算Ｋｓ，并详细分析了桨的几何尺寸和流体剪切稀化特性对Ｋｓ的影响。最后的对包括经典的Ｍｅｔｚｎｅｒ方法在内的不同求解Ｋｓ的方法做了分析比较。