桨叶层数不同对搅拌槽内流动场的影响论文14471.doc

分类号 单位代码 密 级 学 号 学生毕业设计 题 目桨叶层数不同对搅拌槽内流动场的影响作 者院 (系)专 业指导教师答辩日期年 月 日摘 要 搅拌设备广泛应用于工业生产中，如石油、化工、制药、食品、采矿、冶金、废水处理等领域。搅拌操作是工业反应过程中的重要环节，在促进槽内物料流动，使得搅拌槽内物料均匀分布的同时，能够增大传热和传质系数，加速反应的进行。因此，研究搅拌槽的流动性能在工业生产中具有重要意义。 本文采用计算流体力学（cfd）对搅拌槽内的三维流动场进行了数值模拟，考察了六直叶圆盘桨（6dt）、六叶上斜叶桨（6pdtu）、六叶下斜叶桨（6pdtd）在桨叶层数不同时，搅拌槽内的流动场及搅拌功率。结果表明：相同类型的桨叶，在桨叶层数不同时，槽内的流场不同，双层桨时流体的速度分布比单层桨时更加均匀，功率消耗更大。（1）6dt桨在单层时，在整个搅拌槽内形成四个对称的漩涡，功率消耗比较大；6dt+6dt桨由于上下桨间的干涉作用，两层桨之间漩涡不明显，功率消耗更大。（2）6pdtu在单层时，其上扬作用使搅拌槽内形成四个大小不等的漩涡，功率消耗最小；6pdtu+6pdtu桨由于两层桨的相互作用，槽内流体混合均匀。（3）6pdtd桨在单层时，由于其下压作用使槽内流体在搅拌槽内仍然形成四个漩涡，且下部漩涡较小，上部漩涡较大；6pdtd+6pdtd桨由于两层桨相互干涉，在整个搅拌槽出现六个循环流动。研究结果对搅拌反应器的优化设计具有一定的参考价值。关 键 词：搅拌槽，桨叶层数，计算流体动力学（cfd），数值模拟inumerical analysis of the flow field of stirred tank with different impeller layersabstractmixing equipments are widely used in industrial production, such as petroleum, chemical industry, pharmaceutical industry, foodstuff, mining, metallurgy, waste-water treatment, etc. stirring operation is an important part of the industry reaction process, it can promote the flowing of materials and make materials distributing evenly in the stirred tank. at the same time, it can increase the coefficient of heat and mass transfer and accelerate the reaction.therefore, it has an important significance in industrial production to study the flow characteristics of stirred tank.three-dimensional flow fields in stirred tank has been simulated numerically by utlizing computational fluid dynamics (cfd), when the impeller layers are different, the flow characteristics and the mixing power demand of 6dt, 6pdtu and 6pdtd is researched. the results indicate that: with the same type of impeller, the flow field can be changed due to the different impeller layers. the liquid with dual impellers is mixed evenly than single impeller and it costs maximum power demand. (1) as for the 6dt with single impeller, four symmetrical vortexes are formed in the stirred tank and the power demand is larger. when the impeller is 6dt+6dt, because of the interference function,the vortexes between the upper and the lower impeller is not obvious, it costs maximum power demand.(2)as for the 6pdtu with single impeller, its upward function results in four vortexes with different size in the stirred tank; when the impeller is 6pdtu+6pdtu, because of the mutual interference function between the upper and the lower impeller, the liquid is mixed evenly and the power demand is the least in the stirred tank. (3)as for the 6pdtd with single impeller, its downward function also result in four vortexes in the stirred tank, and the range of the upper vortexes is smaller, while the range of the lower vortexes is larger; when the impeller is 6pdtd+6pdtd,there are six vortexes in the stirred tank because of the interference function between the upper and the lower impeller. the results are of importance to the design and optimization of the stirred reactor.key words: stirred tank, impeller layers, computational fluid dynamics(cfd), numerical simulation目 录1绪论11.1搅拌槽简介11.2 搅拌槽内流动特性研究现状11.2.1 宏观流动场11.2.2微观特性研究21.2.3 搅拌槽内流动特性的发展22计算流体动力学42.1 cfd技术原理介绍42.1.1 cfd简介42.2 cfd技术在搅拌槽中应用43研究目的、研究内容及研究方法73.1研究目的73.2研究的主要内容73.3研究方法83.3.1 湍流模型83.3.2 运动部件和静止部件之间相互作用的解决方法84搅拌槽内流场的数值模拟94.1建立几何模型94.2网格划分114.3数值求解步骤134.4数值模拟结果分析134.4.1 6dt桨时不同桨叶层数的数值模拟结果分析134.4.2 6pdtu桨时不同桨叶层数的数值模拟结果分析154.4.3 6pdtd桨时不同桨叶层数的数值模拟结果分析164.4.4 搅拌功率p和功率准数np的比较185总结与展望205.1总结20iii5.2展望20参考文献21致谢22v大学本科毕业论文1绪 论1.1 搅拌槽简介 搅拌设备在石油、化工、制药、食品、采矿、冶金、废水处理等领域应用广泛。在工业生产中，搅拌设备在许多场合是作为反应器来使用的，尤其是在三大合成材料（合成橡胶、合成纤维、合成塑料）的生产中，采用搅拌设备作为反应器的，约占反应器总数的85%以上，但更大量的搅拌设备是用于物料的混合、传热、传质以及制备乳液、悬浮液等1。 搅拌是通过桨叶的旋转向反应器内输入机械能，可以使两种或多种不同的物料互相分散，从而达到均匀混合；也可以加速传热和相间传质过程2。其结构主要包括搅拌槽、搅拌器以及挡板等。搅拌槽的作用主要是为物料反应提供合适的空间，其筒体基本上是圆筒。搅拌器又称搅拌桨或搅拌叶轮，是搅拌反应器的关键部件，其功能是提供过程所需要的能量和适宜的流动状态。在搅拌槽中设置挡板主要是为了消除流体的打旋现象，从而使搅拌槽中的流体流动更加均匀。1.2 搅拌槽内流动特性研究现状在搅拌设备中，由于旋转桨叶和静止挡板的相互作用，使搅拌设备的流场具有三维和高度不稳定的随机湍流特点，难以从理论方面来预测搅拌设备的优劣，也无法对搅拌设备的优化设计提供指导。因此对于搅拌槽的研究还处在半经验半理论阶段。近年来，许多国内外学者从实验研究和数值模拟两方面对搅拌设备进行了广泛的研究。在实验研究方面，早期的研究者由于测试手段的落后及差异，所得到的结果有较大的分歧。随着激光多普勒测速的发展和完善，对搅拌槽内流动特性的研究逐渐深入，实验结果也逐渐趋于一致，但实验方法存在着过程复杂、周期长、生产成本高等缺陷。而数值模拟方法可以详细获取搅拌设备内流体的流动状况，从更微观更本质的角度对搅拌设备进行研究，成为研究搅拌槽流动性能必不可少的工具。1.2.1 宏观流动场在搅拌槽内流体进行着复杂的三维流动，可以将其分为径向流、轴向流和混合流，根据搅拌桨所产生的宏观流动场的形态，将其分为径向流搅拌桨、轴向流搅拌桨和混合流搅拌桨。（1）对于径向流搅拌桨来说，流体在桨叶的作用下沿径向流动，到达槽壁后与槽壁碰撞，分成两部分流体分别在搅拌槽的底部和上部各产生一个循环区，如图1-1所示。对于所有的径向流搅拌桨来说，圆盘是产生径向流的主要原因，典型的径向流搅拌桨是rushton桨。（2）对于轴向流搅拌桨来说，流体在桨叶的作用下沿轴向运动，到达槽底与槽底碰撞后，方向改变沿槽壁向上流动，流经整个槽体后，最终沿轴向下流回叶轮区，形成简单的单循环上下流动，如图1-2所示。典型的轴向流搅拌桨如a310，cby等。 图1-1 径向流搅拌桨流型图 图1-2 轴向流搅拌桨流型图1.2.2 微观特性研究微观特性研究主要是指用湍流强度、湍流尺度、耗散率等物理量来描述湍流脉动特性。相比宏观流动场的研究，微观特性研究发展缓慢。近年来，随着cfd研究的广泛开展，微观特性研究也逐渐开展起来。1.2.3 搅拌槽内流动特性的发展从1982年harvey3第一次将计算流体力学（cfd）引入到搅拌槽内二维流动场的预测开始，对搅拌槽内流动特性的数值模拟研究迅速发展起来。随着工业技术的发展，对过程中流体混合过程、传热及传质的要求更高，传统的经验放大设计方法的可靠性受到了质疑。因此，对搅拌反应设备的优化设计需要建立更加可靠的放大准则。对搅拌槽内流动场特性的数值模拟研究，经历了从二维数值模拟到三维数值模拟，从单向流到多向流，边界条件设置从依靠实验数据的“黑箱”模型发展到动量源法，再到无需依靠实验数据的多重参考系法和滑移网格等技术。近年来非结构化网格技术极大的促进了cfd的发展，数值模拟以其优势成为一种独立有效的研究工具。2 计算流体动力学2.1 cfd技术原理介绍2.1.1 cfd简介计算流体力学cfd(computational fluid dynamics，cfd)是近代流体力学、数值流体力学和计算机科学结合的产物，是一门具有强大生命力的边缘科学，可以看做是在质量守恒方程、动量守恒方程和能量守恒方程控制下对流体的数值模拟。它以电子计算机为工具，应用各种离散化的计算方法，对流体力学的各种实际问题进行数值试验、计算机模拟和分析研究，以解决各种实际问题4。它从基本物理定理出发，在很大程度上替代了耗资巨大的流体动力学实验设备，在科学研究和工程技术中产生巨大影响。cfd是流体力学的一个分支，用于求解固定几何形状空间内流体的动量、热量和质量方程以及其他相关的方程，并通过计算机模拟获得某种流体在特定条件下的有关数据5。cfd计算相对于实验研究，具有成本低、速度快、资料完备、可以模拟真实及理想条件等优点，从而使cfd成为研究各种流体现象、流动过程和结构设计的有力工具。cfd以经典流体动力学和数值计算方法为基础，通过计算机数值模拟和图像显示，在时间和空间上定量描述流体流动，具有理论性和实践性双重特点。利用cfd对搅拌槽进行数值模拟，可以使槽内现象可视化，人们可以直观的了解槽内的流动情况，消除存在的问题，对搅拌槽进行优化设计，为新型高效搅拌设备的研制开辟一个新途径。2.2 cfd技术在搅拌槽中应用搅拌槽由于其内部流动的复杂性，搅拌混合尚未形成完善的理论体系，对搅拌槽等混合设备的放大设计，经验成分往往大于理论计算。随着其广泛应用，经验放大设计方法的可靠性受到前所未有的挑战，对搅拌槽内部流场有必要进行更深入的研究。近年来利用cfd方法研究搅拌槽内的流场发展很快，利用这种方法不仅可以节约大量的研究经费，而且还可以获得实验手段所不能得到的数据。夏建业等6模拟比较了不同桨叶组合搅拌槽内的气液氧传递系数以及在有氧消耗情况下搅拌反应器中氧的浓度分布，发现采用上层为轴流式的三宽叶搅拌桨，底层为径流式的六弯叶圆盘透平桨的组合形式氧供应能力最强。李良超等7采用在欧拉-欧拉双流体模型的基础上耦合气泡数密度(bnd)函数模型，对双层组合桨气液搅拌槽内气泡尺寸和局部气含率进行了cfd模拟，得到了双层桨搅拌槽内气液两相流场、气泡尺寸和气含率分布规律。高勇等8运用cfd方法，对中心龙卷流型搅拌槽内部流场进行数值模拟，分析比较了三种结构的导流板对槽内流体流动状况的影响，计算了搅拌功率准数并与已有研究成果进行了比较。搅拌通常在搅拌槽内进行，通过搅拌桨的旋转把机械能转化为流体的动能，带动槽内流体流动来完成传热传质等过程。如何准确地模拟旋转的搅拌桨与静止的挡板之间的相互运动是研究搅拌流场需要重点解决的问题。过去近三十年来，为了解决这一问题，提出几种不同的模拟方法，这个过程就是cfd技术不断、发展完善的过程。这些方法主要有：桨叶边界条件法，动量源法，内外迭代法，多重参考系法，滑移网格法等9。（1）桨叶边界条件法桨叶边界条件法（impeller boundary condition，ibc）又称“黑箱”模型法，是最早使用的一种搅拌数值模拟方法。该方法是在计算时不考虑桨叶区域内流动情况，而是以实验数据的平均值作为边界条件在桨叶表面进行赋值，以此来计算出整个搅拌槽内的流动情况。middleton等10采用ibc法首次对含有连续-竞争反应的搅拌槽内的流场进行了三维模拟，获得了与实验比较吻合的结果。ibc法是一种稳态算法，但是需要实验数据作为边界条件，摆脱不了对实验数据的依赖，模拟精度有限。（2）动量源法 动量源法（momentum source model，ms）是尝试摆脱实验数据的束缚而提出的第一种搅拌数值模拟方法，同ibc方法一样，ms法不考虑搅拌桨的具体结构，把桨叶对流体的作用看做流体动量的产生源，采用切向方向的附加“源”代替六直叶涡轮桨的作用，其边界条件的确立也不需要实验数据。morud11等采用ms法对有挡板的rushton桨搅拌槽内的气-液两相湍流流场进行了研究，模拟结果与实验值定性吻合。ms法无法研究桨叶区域内流体的流动情况，只能对搅拌槽其他区域内的流场进行研究，因而不能反映搅拌槽内的真实流动情况。（3）内外迭代法内外迭代法（inner-outer method，io）是brucato12 提出的一种无需依赖实验数据的搅拌数值模拟方法，它将计算域分为包含桨叶的内区和包含挡板的外区两部分。内区内的计算是采用搅拌桨速度旋转的坐标系来得到整个内区内的流动场，同时也得到了内区边界上的速度、湍流动能和耗散率。再以该边界上的值作为边界条件对外环进行计算，计算在静止坐标系下进行。通过在两个区域之间交替迭代计算，获得一个收敛结果。内外迭代法不需要实验数据，可以对搅拌槽的流动场进行整体模拟。但是这种方法需要试差迭代，收敛速度较慢。（4）多重参考系法多重参考系法（multiple reference frame，mrf）是luo等13 对io法进行了修正提出的，该方法将计算域划分为两个区域，桨叶所在的区域计算采用旋转坐标系，静止的挡板和其他区域计算使用静止参考系。多重参考系法划分的两个区域部分是没有重叠的，也不需要内外迭代过程，内外区域之间速度匹配直接通过在交界面上的坐标系转换来实现，使得计算变得更加简单。mrf是一种稳态算法，实现了对搅拌槽内整体流动场的模拟，而且不需要实验方法辅助，计算工作量小。适用于桨叶和挡板之间相互作用较小的体系。（5）滑移网格法滑移网格法（sliding mesh，sm）也是一种非稳态算法，适用于非定常流动问题的求解。该方法与mrf相比，对区域的划分是相同的，不同之处在于在两个区域的交界面处，滑移网格法有网格之间的相对滑移。micale等14 利用sm法准确地模拟了双层rushton桨搅拌槽内的流动类型。综上所述，旋转区域和静止区域的处理方法各有优缺点，就目前的计算机条件而言，多重参考系法（mrf）是一种既可行又可靠的处理旋转域和静止域的方法。3 研究目的、研究内容及研究方法3.1 研究目的由于搅拌槽的复杂性和多样性，在自动化选型和设计上，一直依靠经验知识人工完成，智能化水平不高，导致设计周期较长。采用cfd方法对实验过程进行数值模拟，在一定程度上有效的克服了实验研究的不足，可以使槽内的实际现象可视化。随着cfd预测技术的发展，利用数值模拟的方法获得搅拌槽内流场信息已成为现实。通过对搅拌槽内三维流场的数值模拟，深入了解搅拌槽内流动特性，通过分析不同桨叶层数对搅拌槽内流动场的影响，分析比较搅拌功率与搅拌准数，得出较为合理的桨叶层数，为后续的工业实际应用和搅拌设备的优化设计提供一定的参考价值。本文运用ansys-cfx对不同桨叶层数的搅拌槽流动场进行了数值模拟，研究了不同桨叶层数对流场的影响，对比分析了搅拌功率与搅拌准数，为提高搅拌槽效率以及优化搅拌设备的设计提供参考。3.2 研究的主要内容（1）搅拌槽及其流动性能研究现状简介；（2）介绍计算流体动力学(cfd)的原理；（3）研究计算流体动力学(cfd)在搅拌槽中的应用；（4）利用ansys workbench建立几何模型，并进行网格划分；（5）利用cfd对搅拌槽内的三维流场进行数值模拟：单层桨时，研究六直叶圆盘桨（6dt）、六叶上斜叶桨（6pdtu）、六叶下斜叶桨（6pdtd）时搅拌槽内的流场特性；双层桨时，研究6dt+6dt、6pdtu+ 6pdtu、6pdtd+6pdtd三种组合下搅拌槽内流体的流动特性，对比分析桨叶层数不同时对搅拌槽流动场的影响； 分析比较桨叶层数不同时的搅拌功率与搅拌准数。3.3 研究方法3.3.1 湍流模型分析流体流动最常见的模型是湍流模型，它是基于连续性方程、n-s方程和标准k-模型来模拟搅拌槽内的湍流流动。湍流模型可以更好的解决高应变率和流线弯曲程度较大的流体流动。而在壁面处，需要使用标准壁面函数法来确定固壁附近流体的流动。标准壁面函数法14的基本思想是：只在湍流核心区使用rng k-模型求解，对近壁区流体流动不进行求解，而是直接使用半经验公式将壁面处的物理量和湍流核心区内的求解变量联系起来，直接得到与壁面相邻控制体积的节点变量值。3.3.2 运动部件和静止部件之间相互作用的解决方法 运动部件和静止部件之间的相互作用采用多重参考系（mrf）法来处理，即将计算域分为内部旋转域和外部静止域两个区域，内部旋转域包括运动的桨叶和附近的流体，在旋转坐标系下求解连续方程和动量方程式；外部静止域包括静止的挡板和槽体，在静止坐标系下求解连续方程和动量方程式。静止区域和旋转区域分别各自求解方程，在两个域的交界面通过插值来实现质量、动量和能量交换16。在用此方法进行数值模拟时，将内部旋转区域内的流体设定为与搅拌桨相同转速进行旋转，将外部静止区域内的流体设定为静止的。搅拌轴处于静止区域，相对于区域内流体是运动的；搅拌桨处于运动区域，转速和周围的流体相同，相对于区域内流体是静止的。多重参考系法示意图如图3-1所示。 （a） 内部区域（采用旋转坐标系）（b） 外部区域（采用静止坐标系）es：两个域的界面图3-1 多重参考系法示意图4 搅拌槽内流场的数值模拟4.1 建立几何模型模拟对象为四周带有挡板的平底搅拌槽，其结构如图4-1所示，几何尺寸如表4-1所示。由于流动的对称性，模拟时选用1/2搅拌槽作为计算域，搅拌槽模型及桨叶模型如图4-2所示。用于模拟计算的工作介质为水，其密度为997 kg/m3，模拟时搅拌转速为300 r/min。表4-1 搅拌槽几何尺寸名称符号数值（mm）筒体直径t300筒体高度h450桨叶直径d100（=t/3）转轴直径d114叶轮圆盘直径d275圆盘厚度3叶片长度l25下层桨距槽底距离l1150两层桨间距l2150挡板数目nb4个挡板宽度wb25（=t/12） 图4-1 双层桨搅拌槽结构示意图 （a） 双层桨搅拌槽 （b） 6dt桨 （c） 6pdtu桨 （d） 6pdtd桨 图4-2 双层桨搅拌槽及桨叶模型图4.2 网格划分目前在复杂区域生成网格的方法总的来说可以分为两大类：结构化网格和非结构化网格。当计算区域的几何形状比较规则时采用结构化网格进行划分区域；当计算趋于复杂时，即使应用专门的网格技术也难以处理所求解的不规则区域，此时采用非结构化网格划分区域，使局部加密更加容易。由于模拟所选搅拌桨的结构不规则，从计算量、网格生成难易程度、收敛性三方面考虑，采用非结构四面体网格，对桨叶、挡板进行网格细化，以便能更好的捕捉其附近流体的流动特性。在所有壁面处采用inflation网格特性，保证了壁面处的速度梯度变化不至于太大。图4-3为搅拌槽及桨叶网格划分图。 (a) 双层桨搅拌槽网格 (b) 6dt桨网格 (c) 6pdtu桨网格 （d） 6pdtd桨网格图4-3 双层桨搅拌槽及桨叶网格4.3 数值求解步骤将划分好的网格导cfx中，进行数值求解。（1）设定边界条件由于应用mrf方法进行模拟，所以将动区域（fluid-move）内的流体设定以搅拌桨相同的转速进行转动，而静止区域（fluid-static）内的流体是静止的。将搅拌桨和搅拌轴定义为动边界，边界类型均为壁面边界（wall）。其中搅拌轴处于静止流体区域内，相对于区域内流体是运动的；搅拌桨处于运动流体区域，且和周围的流体以同样的转速进行转动，因此相对于区域内流体是静止的。将挡板的外表面与槽壁定义为静止壁面边界条件。（2）设定收敛残差对流动场进行模拟时，将各流动变量的收敛残差定义为。（3）设定时间步长应用隐式方法对非稳态问题进行模拟计算时，可以选取较大的时间步长，但是为了保证模拟结果的精确性，时间步长不应该太大，过大往往会导致计算的不收敛，为了计算的稳定性和精确性，采用固定的保守时间步长值，为0.5s。（4）迭代计算设定一定的迭代步数，进行迭代计算，直至质量、动量和湍流方程收敛。4.4 数值模拟结果分析4.4.1 6dt桨时不同桨叶层数的数值模拟结果分析6dt桨属于径向流搅拌桨，其功率消耗大，剪切应力强，又具有排出能力，因此适用于既要有强的剪切，又要有一定循环流量的场合。图4-4所示分别是6dt桨、6dt+6dt桨时搅拌槽内纵截面处速度矢量图。由图(a)可见：单层桨时，旋转桨叶把流体从轴向吸入而沿径向排出。由于挡板的作用，从桨叶排出的流体在遇到槽壁后沿槽壁一部分向上运动，一部分向下运动，在槽内形成四个循环流动。由图(b)可见：双层桨时，整个搅拌槽内的流型属于径向流，旋转桨叶在叶轮区产生径向射流，叶片端部的水平射流在流动过程中夹带周围流体流动。由于两桨的相互干涉，在两层桨之间形成的两个漩涡受到空间距离的制约，有连接的趋势，在上层桨和液面、下层桨和槽底之间分别形成四个大小基本相等的漩涡。 （a）6dt桨 （b）6dt+6dt桨 图4-4 6dt桨、6dt+6dt桨时搅拌槽内纵截面处速度矢量图图4-5所示分别是6dt桨、6dt+6dt桨时搅拌槽内纵截面处速度云图，由图（a）可见：单层桨时，搅拌槽上部速度较低，槽内速度不均匀。最大速度出现在桨叶端部附近，即桨叶端部附近流体的流动比较剧烈。由图（b）可见：双层桨时，流体在高速旋转时与周围速度较低的流体存在较大的速度梯度，槽内速度混合更均匀。最大速度出现在桨叶与流体接触的边缘，此处的动能最大，流体所受扰动最激烈。 （a）6dt桨 （b）6dt+6dt桨 图4-5 6dt桨、6dt+6dt桨时搅拌槽内纵截面处速度云图4.4.2 6pdtu桨时不同桨叶层数的数值模拟结果分析6pdtu桨属于轴向流搅拌桨，对流体具有上扬作用，使搅拌槽内的漩涡分界面在搅拌桨之上，流体产生高速向上的轴向流。将其应用于气液分散，可以将气体在整个搅拌槽内分散均匀。图4-6所示分别是6pdtu桨、6pdtu+6pdtu桨时搅拌槽内纵截面处速度矢量图，由图（a）可见：单层桨时，桨叶高速旋转将流体“甩出”，撞击槽壁后在槽内分成了上下两部分，由于上斜叶桨对流体具有上扬作用，因此在搅拌槽内纵截面形成两大两小四个漩涡。由图（b）可见：双层桨时，6pdtu桨对流体的上扬作用，使得流体产生高速向上的轴向流，将下层桨附近的流体送到上层桨周围，由于上层桨的继续上扬作用，使流体和上层桨甩出的流体混合后沿斜上方排出，到达槽壁与槽壁碰撞后沿槽壁向上、向下流动，在两层桨的相互作用下形成六个循环流动。 （a）6pdtu桨 （b）6pdtu+6pdtu桨图4-6 6pdtu桨、6pdtu+6pdtu桨时搅拌槽内纵截面处速度矢量图 图4-7所示分别是6pdtu桨、6pdtu+6pdtu桨时搅拌槽内纵截面处速度云图，由图（a）可见：单层桨时，由于桨叶对流体的挤压作用，因此槽内最大速度分布在桨叶和叶轮的上边缘。桨叶的上扬作用使流体在搅拌槽下部混合均匀。由图（b）可见：双层桨时，旋转桨叶在排出流体的同时，对流体有上扬作用，两层桨之间的轴向作用较强，有利于整个搅拌槽的混合。槽内速度较大区域集中分布在叶片和桨叶的边缘，双层桨间的综合速度要高于上桨与液面、下桨与槽底间的速度。 （a）6pdtu桨 （b）6pdtu+6pdtu桨图4-7 6pdtu桨、6pdtu+6pdtu桨时搅拌槽内纵截面处速度云图4.4.3 6pdtd桨时不同桨叶层数的数值模拟结果分析6pdtd桨对流体有下压作用，使搅拌槽内的漩涡分界面在搅拌桨之下，将其应用于固液反应中，能够将槽底的固体带动起来，避免固体沉积。图4-8所示分别是6pdtd桨、6pdtd+6pdtd桨时搅拌槽内纵截面处速度矢量图，由图(a)可见：下斜叶桨对流体具有下压作用，旋转桨叶使流体向斜下方流动，槽内流体在搅拌槽内仍然形成上下四个漩涡，但是下部漩涡较小，上部漩涡较大。由图(b)可见：双层桨时，由于6pdtd桨对流体的下压作用，使旋转桨叶甩出的流体大部分沿斜下方流动，到达槽壁后沿槽壁向上、向下流动，和下层桨叶甩出的流体混合后流动方向改变，沿径向流回槽壁，径向流动强烈。由于两层桨相互干涉，两层桨之间的涡环不明显，整个搅拌槽出现六个循环流动。 （a）6pdtd桨 （b）6pdtd+6pdtd桨 图4-8 6pdtd桨、6pdtd+6pdtd桨时搅拌槽内纵截面处速度矢量图图4-9所示分别是6pdtd桨、6pdtd+6pdtd桨时搅拌槽内纵截面处速度云图，由图(a)可见：单层桨时，由于桨叶的下压作用，在搅拌槽下部速度分布均匀,槽内最大速度在桨叶叶端附近。由图(b)可见：双层桨时，槽内速度较高区域分布在叶片和桨叶的边缘，双层桨间的综合速度高于其它区域的速度。 （a）6pdtd桨 （b）6pdtd+6pdtd桨图4-9 6pdtd桨、6pdtd+6pdtd桨时搅拌槽内纵截面处速度云图 4.4.4 搅拌功率p和功率准数np的比较搅拌功率是指搅拌槽在进行搅拌时，在单位时间里，输入槽内对物料做功并使之发生流动的能量，它不包括在轴封和传动装置中消耗的能量，搅拌功率与搅拌槽的结构，搅拌桨的形状、尺寸和转速，液体性质和内部附件（有无挡板和其它障碍物），搅拌桨在槽内位置以及重力加速度等有关17。通过模拟可以得到扭矩值，搅拌功率p与扭矩m的转换公式： （4-1）由公式（4-1）得出： (4-2)式中：p为搅拌功率（kw），m为搅拌扭矩（），n为搅拌转速（r/min）。功率准数np是雷诺数re的函数，其中 (4-3)式中：d为搅拌槽的直径，为流体的密度，为流体的黏度。对于一定几何结构的搅拌槽，np与re的函数关系可由实验测定，将这函数关系绘成曲线称为功率曲线。功率准数np与搅拌功率p之间的关系为： (4.4)式中：p为搅拌功率（w），n为搅拌转速（r/s），d为搅拌槽的直径（m），为流体的密度（kg/m3）。由公式（4-2）和公式（4-4）可以计算出各种情况下的搅拌功率和功率准数如表4-2所示。在计算时由于轴的扭矩与桨叶的扭矩相比太小，故在计算功率时只考虑桨叶的扭矩。表4-2 搅拌功率与功率准数 项目组合桨桨叶扭矩m（）功率p（w）功率准数np6dt6.4710-22.0326.696dt+6dt8.2610-22.6088.586pdtu1.4910-20.4681.546pdtu+6pdtu2.3010-20.7232.386pdtd1.5710-20.4921.626pdtd+6pdtd2.7010-20.8482.79由表可见：6dt桨功率准数为6.69，6dt+6dt功率准数较大为8.58，其功率准数之比约为0.78:1；6pdtu桨功率准数在单层桨情况下最小，为1.54，6pdtu+6pdtu桨功率准数在双层桨情况下也为最小，为2.38，但是比单层时功率准数大，其功率准数之比约为0.65:1；6pdtd桨在单层桨情况下功率准数居中，为1.62，在双层桨情况为2.79，其功率准数之比约0.58:1。由于搅拌层数的增加，导致在搅拌槽运转过程中所消耗的功率相对增加，上为了克服两层桨之间的干扰作用需要消耗的功率，因而总体来说，双层桨时功率准数比单层桨时功率准数大。因为6dt桨属于典型的径向流搅拌桨，它的剪切力和排出能力较强，故能量损耗也最大。5 总结与展望5.1 总结利用cfd方法对桨叶层数不同时搅拌槽内流体的流动特性进行了数值模拟，得出以下结论：（1）6dt桨在整个搅拌槽内形成四个漩涡，速度分布不均匀，功率消耗在单层桨中最大；6dt+6dt桨由于上下桨间的干涉作用，两层桨间的漩涡受到空间距离的制约，有连接的趋势，在上层桨和液面、下层桨和槽底之间形成的四个漩涡大小基本相等。槽内速度分布比较均匀，功率消耗比单层时大。（2）6pdtu桨的上扬作用使槽内流体在整个搅拌槽内形成两大两小四个漩涡，流体形成的流动场以轴向流为主。由于桨叶的上扬作用，使搅拌槽下部区域内速度分布均匀，在单层情况下功率准数最小；6pdtu+6pdtu桨由于其对流体的两次上扬作用可以使流体分散的更加均匀，搅拌槽内整体混合效果较理想，功率消耗比单层时略大。（3）6pdtd桨对流体的下压作用使槽内流体在搅拌槽内仍然形成上下两个漩涡，但是下部漩涡较小，上部漩涡较大；6pdtd+6pdtd桨由于两层桨相互干涉，两层桨之间的涡环不明显，整个搅拌槽出现六个循环流动，速度分布比单层时均匀。5.2 展望搅拌操作是化工过程中最普遍最重要的过程，因此对搅拌槽内流动特性的研究有重要的意义。前人对搅拌设备的研究为我们积累下了丰富的经验，我们才能够依靠这些经验进行有效的设计与放大。随着时代的发展，对搅拌设备的要求越来越高，通过对桨叶层数不同时搅拌槽内流动特性的研究，可以为搅拌槽的优化设计及应用提供一定的参考作用。本论文利用数值模拟的方法研究了在桨叶层数不同时，6dt桨、6pdtu桨以及6pdtd桨搅拌槽内流动场特性，取得了一定的成果，但由于条件和时间的限制，研究结果还存在很多不足。在以后的研究中，要继续扩展对搅拌槽流动特性的研究范围，从计算流体力学（cfd）角度研究其它因素的影响，以及从计算方法、几何模型等方面提高搅拌槽流动特性的预测精度，将是以后对搅拌槽数值模拟研究的主要发展方向。参考文献1 陈志平，章序文，林兴华，等.搅拌与混合设备设计选用手册m.北京：化学工业出版社，2004.2 王凯，虞军.搅拌设备m. 北京：化学工业出版社，2003.3 harvey p s，greaves m.turbulent flow in an agitated vessel.part i ：a predictive model. trans inst chem .eng，1982，60a: 195-1594.4 孙智一，吴晓蓉.计算流体力学数值模拟方法的探讨及应用j.水利科技与经济，2008，14（2）：126-128.5 黄永春，唐军，谢清若.计算流体力学在化学工程中的应用j.现代化工，2007，27（5）：65-74.6 夏建业，唐演，张嗣良.搅拌槽中不同搅拌系统气液氧传递的