Mass and heat transfer behavior of oscillating helical coils in relation to heterogeneous reactor design

Rates of mass and heat transfer at vibrating helical coils were studied by the electrochemical technique with the object of using helical coils as heat exchanger/reactor for conducting liquid–solid diffusion controlled reactions. Variables studied were frequency and amplitude of vibration, tube diameter, and superimposed axial flow velocity. The data for vibrating coil (batch reactor) were correlated for 59 < < 4965; Sc = 2314 by the equation: The data were found to be consistent with the analogy model. For vibrated helical coils with superimposed axial flow, the data were correlated by the equation: Importance of the present results in the design and operation of heterogeneous reactors used to conduct diffusion controlled exothermic reactions involving heat sensitive materials was pointed out. Also the importance of the present results in the design and operation of shell and helical tube heat exchangers of improved performance was highlighted. © 2017 American Institute of Chemical Engineers AIChE J, 63: 3141–3149, 2017     

Modeling of cavern formation in yield stress fluids in stirred tanks

Prediction of cavern formation in yield stress fluids in stirred tanks is of great importance for optimization. A new torus model is developed and then validated by experimental data and computational fluid dynamics simulation. Unlike existing mathematical models, the new torus model assumes that the circular center of the torus should not be outside the impeller swept region as the Reynolds number (Re) increases. Hence the cavern boundary is shaped like an apple torus rather than a horn torus. The new model also considers the cavern‐vessel interactions. At relatively high Re, the new model predicts cavern shape and size better than other models. It correctly captures the cavern outline at various Re, which verified the assumption about torus center. The new model is then used to identify the influence of rheological parameters on cavern formation, and further extended to the cavern prediction of the dual‐impeller system. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3057–3070, 2014     

一种计算搅拌槽混合时间的新方法

基于对混合时间定义的思考,提出了一种新的定义方法,在湍流流场数值计算的基础上通过求解示踪剂的浓度输运方程,研究了单层涡轮桨搅拌槽内的混合过程。结果表明：搅拌转速和搅拌桨安装位置都影响混合时间的大小,而进料位置对混合时间的影响不大。对于不同的搅拌转速而言,随搅拌转速的增大,相同体积分数对应的混合时间逐渐减小。当搅拌桨安装在槽中间位置时所对应的混合时间最小。利用适宜的尺寸和安装位置的导流筒可有效降低混合时间。     

搅拌槽中酯化反应热失控的CFD模拟

热失控是化工过程中常见的安全风险之一。在间歇釜式反应器中,桨叶的机械转动可以增强流体的循环流动、湍流强度、混合程度以及传热,进而有效防范热失控。防控效果与反应器结构和搅拌桨型密切相关。针对丙酸异丙酯酯化反应,采用计算流体力学模拟研究了桨型(Rushton桨、30o PBT桨及60o PBT桨)、转动方向和挡板对釜式反应器内温度演化的影响,从流动结构方面分析了原因。基于散度的失控判据比较了三种搅拌桨抑制热失控的能力,抑制能力为Rushton桨＞30° PBTD桨＞60° PBTD桨。本研究可为搅拌反应器热失控的优化设计提供一定的理论依据。     

Transitional flow in a Rushton turbine stirred tank

The way in which the single phase flow of Newtonian liquids in the vicinity of the impeller in a Rushton turbine stirred tank goes through a laminar‐turbulent transition has been studied in detail experimentally (with Particle Image Velocimetry) as well as computationally. For Reynolds numbers equal to or higher than 6000, the average velocities and velocity fluctuation levels scale well with the impeller tip speed, that is, show Reynolds independent behavior. Surprising flow structures were measured—and confirmed through independent experimental repetitions—at Reynolds numbers around 1300. Upon reducing the Reynolds number from values in the fully turbulent regime, the trailing vortex system behind the impeller blades weakens with the upper vortex weakening much stronger than the lower vortex. Simulations with a variety of methods (direct numerical simulations, transitional turbulence modeling) and software implementations (ANSYS‐Fluent commercial software, lattice‐Boltzmann in‐house software) have only partial success in representing the experimentally observed laminar‐turbulent transition. © 2017 American Institute of Chemical Engineers AIChE J, 63: 3610–3623, 2017     

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

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

双层涡轮桨搅拌反应器内混合时间的大涡模拟

采用计算流体力学(CFD)方法对直径为0.476m双层涡轮桨搅拌反应器内的流动及混合进行了数值模拟,并实验测试了混合过程。利用大涡模拟(LES)及Smagorinsky-Lilly亚格子模型求解湍流流动与示踪剂传递过程,桨叶区域采用滑移网格技术。研究结果表明,大涡模拟得到的示踪剂响应曲线和混合时间与实验结果吻合良好,其预测精度明显优于基于雷诺平均(Reynolds-averaged Navier-Stokes,RANS)的标准k-ε模型的模拟结果。大涡模拟是研究搅拌反应器内非稳态及周期性湍流流动的有效方法。     

Theoretical prediction of gas-liquid mass transfer coefficient, specific area and hold-up in sparged stirred tanks

A prediction method for calculating the volumetric mass transfer coefficient, k_La, in gas-liquid sparged stirred tanks is proposed. A theoretical equation based on Hibie's penetration theory and the isotropic turbulence theory of Kolmogoroff is used for k_L determination. The values of the interfacial area have been calculated from a hold-up theoretical equation and the mean size of the gas bubble. Both Ostwald-De Waele and Casson models are used to describe the rheological properties of the fluid. The model predicts the mass transfer coefficient and the interfacial area values in stirred tank reactors, analysing the influence of different variables. The values of the volumetric mass transfer coefficient can be calculated for different geometries of the reactor, different physicochemical properties of the liquid and under different operational conditions. The capability of prediction has been examined using experimental data available in the literature for Newtonian and non-Newtonian fluids, for very different vessel sizes, different numbers and types of stirrers and a wide range of operational conditions, with very good results.     

基于EMMS模型的搅拌釜内气液两相流数值模拟

采用欧拉-欧拉模型对搅拌釜内气液两相流进行了三维CFD模拟,重点研究了采用不同曳力模型时CFD模拟对搅拌桨附近排出流区两相流动的预测能力。模拟结果表明CFD能准确地预测排出流区的液相速度分布,但采用传统的Schiller-Naumann曳力一定程度上低估了排出流区的气液相间曳力,导致在完全扩散区CFD预测的分布器和桨叶下方区域气含率偏小,而基于气液非均匀结构和能量最小多尺度（EMMS）方法得到的DBS-Global曳力模型能更准确地描述完全扩散区气液搅拌釜内流动情况。与传统曳力模型相比,采用DBS-Global曳力模型能显著提高对气含率的预测。     

Modeling of the 3D hydrodynamics of 2-blade impellers in stirred tanks filled with a highly viscous fluid

This paper deals with mixing in stirred tank reactors by paddle agitators and two-blade impellers with different blade widths. Computational Fluid Dynamics (CFD) is used to solve the 3D hydrodynamics and to obtain at every point the flow patterns, the stress components and the viscous dissipation function. From the latter information, the power consumption is calculated. These results are compared with available experimental data: good agreement is observed. For example, for a paddle agitator, the flow is essentially plane far away from the horizontal walls. That means that, in order to create an axial circulation in the volume of the tank, the impeller height has to be decreased. This work documents the ability of CFD to model changes in flow patterns for three dimensional flows in stirred tanks.     

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

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

基于螺旋片强化的套管换热器性能

基于RNG k-ε模型对螺旋片强化的套管换热器的传热进行模拟,通过模拟结果与文献中的实验结果进行对比来验证模拟的可行性;分析了Reynolds数为2362～16860范围内的螺旋升角α变化对Nusselt数和摩擦阻力系数f的影响;并考察了等泵功率下的综合传热性能PEC值的变化规律。结果表明：Nu和f的平均误差分别为7.1%和1.3%,证明所采用的研究方法是可行的;α在15°～75°范围内,Nu和f均随着α的减小而增大,特别地,当α小于35°时,f随α的减小剧烈增大;在等泵功率下,PEC值为0.84～1.93;α在15°～45°时,α为35°具有较好的综合传热性能,α为55°、65°和75°时,虽然其PEC值比35°时略高,但其Nu与35°时相比要小很多,实际应用中考虑到传热速率的问题,选择35°的螺旋升角较为合适,此时,PEC值为1.26～1.62。另外,为减小f,提出倾斜螺旋片强化的方法;螺旋升角α为35°、螺旋片倾斜角β为10°时,与普通螺旋片相比,Nu基本一致,甚至略大,而f减小了12.5%～14.5%,此时,PEC值为1.38～1.71;场协同理论也很好地验证了这一结果。     

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     

一种新型自吸式搅拌装置及其设计计算方法

提出一种新型的能将釜底液体连续喷洒在容器内壁形成液膜以提高传热速率的自吸式搅拌装置。从流体动力学的基本方程出发,推导出搅拌装置的设计计算方法。用水为介质在直径200 mm,高300 mm的外表面电热套加热的圆柱形敞口容器中分别在自吸式搅拌装置和普通螺旋桨叶轮搅拌下进行了对比蒸发实验。结果表明,采用自吸式搅拌装置可减少蒸发时间47%以上。     

湿法提钒浸出段搅拌反应器结构的优化

清洁提钒工艺中,熟料的湿法浸出是重要的操作单元.浸出搅拌反应器的合理设计与优化,可缩短浸矿时间以及提高浸矿效率.本文通过改变搅拌桨桨叶间的层间距、搅拌桨的安装层数以及安装导流筒等方法,对攀钢集团公司的浸出搅拌反应器进行优化和改进;并结合计算流体动力学(CFD)Fluent商业软件,分别模拟了原浸出搅拌反应器和改进后浸出搅拌反应器的宏观流场结构.结果表明:在层间距C₂为1100mm的原双层搅拌桨浸出搅拌反应器内,流体轴向速度较小,“死区”现象较严重;与原反应器相比,调整双层搅拌桨桨叶之间的层间距C₂为1800mm以及安装3层搅拌桨或导流筒,都可加强反应器内流体的轴向流动,减小“死区”范围,进而改善流场结构的均匀分布,有助于强化流体混合.     

矩形截面螺旋管内气液两相流型转换数值模拟

采用CLSVOF(coupled level set and volume of fluid method)方法,以空气和水为介质对矩形截面螺旋管内气液两相流动进行数值模拟,气相折算速度U_G为0.1～2.5 m·s^-1,液相折算速度U_L为0.09～4.5 m·s^-1.研究螺旋直径、螺旋升角对流型转换边界的影响,并绘制了不同螺旋直径、不同螺旋升角下的流型图.数值结果表明,与传统VOF方法相比,CLSVOF可以得到更精确的相界面;随着螺旋升角的增加,塞状流向泡状流的转换边界向U_L减小的方向进行,但是幅度很小,塞状流向弹状流的转换边界向U_L减小的方向进行;随着螺旋直径的增加,塞状流向泡状流的转换边界向U_L减小的方向进行,塞状流向弹状流的转换边界向U_L减小、U_G增大的方向进行;与Murai流型图相比,流型转换边界有所差异.     

Validation of CFD simulations of a stirred tank using particle image velocimetry data

Methods for validating CFD simulations based on the Reynolds Average Navier-Stokes equation (RANS) against Particle Image Velocimetry (PIV) measurements are investigated and applied to one of the most common problems in the chemical process industry — the prediction of flow field in a stirred vessel. A total of 1024 sequential instantaneous 2D velocity fields along the central axial plane of a stirred vessel with a P-4 axial impeller are obtained through PIV measurement. From the PIV data, the mean velocity, turbulent kinetic energy, Reynolds stresses and dissipation rate fields are extracted. By introducing several tools to quantify the similarities and differences between two-dimensional fields, CFD predictions of the flow field are validated against PIV data. Furthermore, using PIV and LDV data, the effect of boundary conditions on CFD simulation results is examined. The effect of different Reynolds stress closures on the flow prediction is also studied.     

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

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