D-T型双轴卧式搅拌装置半釜持液量时功率特性的数值模拟

1、D-T 型双轴卧式搅拌装置半釜持液量时功率特性的数值模拟冯惠生 1，杨腾 1，李文秀 2，张志刚 2，余国琮 1(1. 天津大学化工学院，天津 300072；2. 沈阳化工大学化学工程学院，沈阳 110142)-摘 要：应用动网格技术以及多相流模型中的 VOF 模型计算分析了一种带搅拌结构的卧式双轴再沸器的功率特*性CFD 模拟结果表明，总的功率准数与相对应的表观雷诺数的关系为 Np7.15( Re ) 0.775，表观雷诺数的指数为0.775，即功率准数的对数与表观雷诺数的对数呈线性关系，但并不成倒数关系；黏度、流速以及搅拌轴位置对功率的大小均有一定影响关键词：卧式搅拌设备；功率特性；动网格

2、；VOF 模型文章编号：0493-2137(2011)1-1-1003-06中图分类号：O351.2文献标志码：ANumerical Simulation of Power Performances in D T Shape HorizontalBiaxial Stirring Tank with Half Kettle Liquid HoldupFENG Hui-sheng1，YANG Teng1，LI Wen-xiu2，ZHANG Zhi-gang2，YU Guo-cong1(1. School of Chemical Engineering and Technology，Tianjin

3、University，Tianjin 300072，China；2. College of Chemical Engineering，Shenyang University of Chemical Technology，Shenyang 110142，China)Abstract：The power performances of a horizontal biaxial stirring device were computed and analyzed，using the-dynamic meshing model and VOF model. The CFD result shows t

4、hat the total power number versus Reynolds numberis Np7.15(Re ) 0.775That is to say，the logarithmic power number decreases linearly with the logarithmic appar-ent Reynolds number in the laminar regime on a log-log plot，but the index of the apparent Reynolds number is0.775 and therefore into a recipr

5、ocal relationship does not exist. It also shows that viscosities，flow rates and posi- tion of the shafts all have some impact on the power consumption.Keywords：horizontal stirring device；power performance；dynamic mesh；VOF model*带有搅拌结构的卧式再沸器是解决高黏度、强热敏性物系的分离纯化难题的重要途径1-3卧式双轴 设备特别适用于高黏、超高黏和带有粉体物系的混合，转动时

6、两轴上的搅拌构件之间以及它们与容器内 壁之间相互刮擦从而具有自清洁功能，更有利于强化 传质传热过程卧式双轴 T 型搅拌器4是一种常见 的结构简单且具有部分自清洁功能的高黏度搅拌设备；日本三菱重工公司生产了搅拌构件为椭圆盘状具 有自清洁功能的卧式双轴搅拌器5；瑞士 List 公司生 产的卧式双轴搅拌机，主搅拌轴上有许多被捏合杆连 在一起的盘片，捏合杆略有倾斜，使物料在进行径向混合的同时能受到一个轴向的输送力，清洁轴装有一排倾斜的捏合框，通过两轴上的元件相互啮合，从而使搅拌器具有自清洁功能6 高黏度流体的搅拌过程中需要很大的能耗，在搅拌设备的设计和放大中，搅拌功率是很重要的参数， 它直接影响搅拌槽

7、内流动状况和传质传热的强度近些年，研究者们主要通过实验的方法对各种卧 式搅拌装置的功率特性进行了研究冯连芳等4,7对卧式双轴 T 型自清洁反应器在非牛顿流体及牛顿流 体中搅拌功率特性进行了实验研究随着计算流体 力学的发展，通过 CFD 方法来预测搅拌功率能够获 得良好的结果Letellier 等8应用 CFD 方法获得了多级搅拌桨结构在 3 种不同放大比例下的功率曲线(功收稿日期：2010-05-05；修回日期：2010-10-12.作者简介：冯惠生（1962），男，博士研究生，副研究员， HYPERLINK mailto: 通讯 杨

8、腾， HYPERLINK mailto:yangteng20042880163 yangteng20042880163 .1004第 44 卷第 11 期天津大学学报率准数-雷诺数)和 Metzner 数，并通过试验进行了验证Kelly 等9建立的 CFD 模型准确预测了不同结构 的轴向流桨的功率准数冯惠生等10通过 CFD 方法预测了 D-T 型双轴卧式搅拌装置满釜时的功率特性 和传热特性，取得了较好结果可以看出 CFD 已经 开始广泛应用于搅拌过程的研究，但是对于卧式搅拌 设备特别是卧式双轴搅拌设备半釜持液量时的功率特性的研究还是很少笔者拟采用 CFD 的方法来研 究一种带推进叶片的双 T

9、 型(D-T 型)卧式双轴搅拌 装置半釜持液量时的搅拌功率特性所以在该有限元 dA 上产生的扭矩为22dT = dFt (r cos) + z(2)而整个搅拌轴所承担的扭矩值为T = ( (r cos ) + z )dFt22(3)A式中 A 为搅拌轴(包括搅拌翅、推进叶片)上的表面面积1 模型建立1.1 几何模型拟采用单级搅拌桨进行计算，通过 SolidWorks构造结构模型，导入 Gambit 中，通过分离、合并等方 法，形成如图 1 所示的卧式双轴搅拌桨模型，实心搅 拌轴以外的区域为模拟区域搅拌轴直径为 44 mm，两搅拌轴的轴心距为 104 mm，搅拌翅顶端距轴心的距 离为 76 mm

10、 ，搅拌槽的直径是 160 mm ，槽 长 为210 mm翅数少的轴定义为清理轴，翅数多的轴定义为搅拌轴图 2 搅拌翅的受力分析Fig.2 Force analysis of stirring wing1.3 网格划分在模拟卧式双轴结构时，以 Tet/Hybrid 为网格单元，采用 Tet/Hybrid 网格化方法对模型进行了网格划分，模拟时的网格数为 3.251051.4 物理模型及边界条件搅拌轴的转动，导致模拟区域不断发生变化，且双轴的搅拌翅相互交叉，使得滑移网格法及多重参考系法不再适合 ，所以 选择 了动 网格 模型 (dynamic mesh)模拟区域内存在气液两相，属于多相流问题，

11、而多相流模型中 VOF 模型适合于分层的自由表面 流，所以选择 VOF 模型同时，该模拟过程属于非稳态过程模拟由于采用 VOF 模型，同时为高黏度流 体的 慢速 流动 ，选择 了非 耦 合 求 解器 (segregated solver)11-12通过 Profile 定义双轴的搅拌转速，操作条件在重力加速度为 9.81 m/s2 下进行计算，进出口边界条件类型为 velocity-inlet 所选用的物料模型气相为空气，液相为黏度改变而其他物性(与水相同)不变的虚构组分模拟过程中搅拌轴与清理轴的转速之比为 14图 3 所示为 t0.2 s 时(液相流速 0.01 m/s，清理轴转速为 19.

12、5 r/min)气液二相区域模拟结果，灰 色部分为液相区域，黑色部分为气相区域液相以相同的流速通过图 2 所示的进出口进入和离开模拟区域图 1 D-T 型卧式双轴搅拌装置的几何结构模型Fig.1 Geometry model of D-T horizontal biaxial stirring device1.2搅拌翅受力分析如图 2 所示，在推进叶片上任取一有限元 dA，其位置可用其中心点 M(r，z)表示作用在有限元上 的力可分为 2 部分，一部分为与中心轴平行的轴向力dFy，另一部分是经过 M 点与中心轴垂直的面上的作用力，该面上的作用力可分为在作用点 M 与中心轴连线方向上的力 dFr

13、 及与该连线垂直方向的力而扭矩的产生为 dFt 作用的结果M 点距中心轴的距离为dFt，l = (r cos ) + z22(1)2011 年 11 月冯惠生等：D-T 型双轴卧式搅拌装置半釜持液量时功率特性的数值模拟1005倍，所以功率消耗主要在清理轴上在搅拌过程中，功率的消耗来自2 方面：一方面，轴(包含搅拌翅、推进叶片)受到压力作用，如果要使它按一定转速转动，就必须提供一定力矩，此部 分功率简称压力功率；另一方面，由于物料自身的黏 度使得轴在转动时需要克服黏性对流体流动产生的 阻力，此部分功率简称黏度功率每个轴的压力功率和黏度功率之和称为此轴的总功率如图 5 和图 6 所示，随着黏度的增

14、加，压力的阻 碍 所产生的功率 消 耗 呈 线性 增长 ；在 黏 度 小 于7 Pas 时，黏性的阻碍所产生的功率消耗与黏度呈 线性关系，而在黏度大于 7 Pas 后，黏性的阻碍所 产生的功率消耗有一定程度的下降；但由于压力部 分产生的功率消耗较多，使得总功率消耗仍呈上升 趋势图 3 t0.2 s 时气液二相区域的模拟结果Fig.3 Simulation results of gas-liquid phase region whent= 0.2 s2结果与讨论通过模拟计算可以求得各种条件下的卧式双轴结构每个轴上的扭矩绝对值，然后分别由下式求算搅 拌功率P = TN(4)式中：T 为力矩在转动轴

15、上的分量，Nm；N 为转速，rad/s 对于卧式双轴搅拌釜来说，由于其双轴搅拌构件是不完全对称的，2 个轴消耗的功率是不相等的2个轴的转速不同，分别对清理轴、搅拌轴的功率以及 两轴的总功率的特性进行研究，又因为两轴的速度成 一定的比例(41)，且转动方向相反(搅拌轴为逆时 针方向转动，而清理轴为顺时针方向转动)，以清理轴为主动轴，则可以以清理轴的转速对总功率进行研 究其中 PtolP1P2，P1 和 P2 分别是清洗轴和搅拌 轴的功率2.1 黏度对功率的影响如图 4 所示，随着黏度的增加，双轴所消耗的总功率及两轴各自的功率都增加，在黏度小于 7 Pas 时，功率消耗与黏度呈线性关系，在黏度大于

16、 7 Pas 时，功率消耗的增长速度有所减缓从模拟的数据结 果可以看出，由于清理轴的转速是搅拌轴转速的 4图 5 清理轴功率消耗与黏度关系Fig.5 Power consumptions of clean-up shaft vsviscosities图 6 搅拌轴功率消耗与黏度关系Fig.6 Power consumptions of main shaft vs viscosities2.2 流体流速对功率的影响如图 7 所示，随着流速的增大，双轴所需的总功率及两轴各自的功率都有减小，且与流速呈近似线性关系如图 8 和图 9 所示，随流速增大，压力作用部 分消耗的功率有所减小，但在流速大于 0

17、.05 m/s 时， 开始趋于平稳；而随流速增大，黏性作用部分消耗的 功率一直在减小，且在流速大于 0.05 m/s 时，有加速减小的趋势图 4 双轴功率消耗与黏度关系Fig.4 Power consumptions of two shafts vs viscosities1006第 44 卷 第 11 期天津大学学报翅的体积未发生变化，但转动出液相的搅拌翅的部分已经附着液相，阻碍转动的进行，从而使功率消耗增 加当清理轴、搅拌轴转动到 和 值分别约等于150和 38时，出现一定波动，主要与清理轴、搅拌轴 (包括搅拌翅、推进叶片部分)的结构以及相对位置有 关由于在总功率的消耗中，清理轴的功率消耗

18、占了 主要部分，总功率的变化趋势与清理轴的功率消耗变 化趋势相似随着黏度的增加，双轴总功率消耗以及 每个轴的功率消耗随着位置的变化有所增加，主要是 由于随黏度的增大，在搅拌过程中，搅拌翅所划过的 位置，液相补入的速度减缓，使得液相整体偏离理想 状态（液面完全水平，在水平面下，液相与清理轴和 搅拌轴完全接触)的程度增大图 7 双轴功率消耗与液相流速关系Fig.7 Power consumptions of two shafts vs flow speeds of liquid图 8 清理轴功率消耗与流速关系Fig.8 Power consumptions of clean-up shaft vs

19、 flow speeds of liquid图 10 双轴转动角度示意Fig.10 Diagram of rotation angle for two shafts图 9 搅拌轴功率消耗与流速关系Fig.9Power consumptions of main shaft vs flow speeds of liquid图 11 双轴总功率消耗随 值变化Total power consumptions of two shafts vs ro-tation angle Fig.112.3不同角度位置功率的变化情况如图10 所示实体部分为初始位置，随着搅拌的进行，搅拌轴与清理轴都各自偏离初始位置，偏

20、离初始位置的值分别用图 10 中所显示的角度 和表示，由于搅拌的周期性及轴的对称性，本文研究搅拌轴转动45以及清理轴转动 180周期内的功率变化情况 如图 11图 13 所示，随着转动角度的变化，双 轴的位置发生变化，所需总功率消耗及各轴功率在周期范围内呈规律性变化当转动角度 很小时，也就是刚开始搅拌时，功率较小，主要原因是：在初始阶段，搅拌翅只有一半浸没在液相中，另一半并未附着 液体，而随着搅拌的进行，虽然浸没在液相中的搅拌图 12 清理轴功率消耗随 值变化Fig.12 Power consumptions of clean-up shaft vsrotation angle 2011 年

21、11 月冯惠生等：D-T 型双轴卧式搅拌装置半釜持液量时功率特性的数值模拟1007置满釜时7, 13结果不同产生该结果的可能原因为：搅拌是在半釜持液量时进行，搅拌轴自身所带的推进 叶片具有一定倾斜角；同时流体存在轴向流动，即以一定速度进入和流出搅拌槽，并在推动叶片作用下产 生更复杂的轴向速度分布图 13 搅拌轴功率消耗随 值变化Fig.13 Power consumptions of main shaft vs rotation angle 功率准数 Np 与表观雷诺数 Re 的关系*2.4各个轴及其总的功率准数和表观雷诺数定义为d，Re * =2N PN =1(5)1图 15 搅拌轴功率准数

22、与表观雷诺数关系Fig.15 Power number of main shaft vs apparent Reynolds number1p1 N 3d 4 L1Pd，Re * =2N N=2(6)22p2 N d L3 42P + PN =(7)12N + Np2 )3 d 4 L (12d 2 N1 + N2 2Re* =(8)式中：N1 和 N2 分别为清理轴和搅拌轴的转速；P1 和P2 分别为清理轴和搅拌轴的功率；d 为搅拌桨直径；L 为釜长；为物料密度；为物料的黏度考察清理轴、搅拌轴的功率准数以及总的功率准数与对应的雷诺数的关系，并进行关联从图 14图*16 中可以看出，在 1 R

23、e 50 时，功率准数的对数-与雷诺数的对数呈线性关系图 16 双轴总功率准数与表观雷诺数关系Fig.16 Total power number of two shafts vs apparentReynolds number3结论-Np12.19(Re1*) 0.775(9)(10) (11)(1) 随着黏度的增加，双轴所需的总功率消耗及-Np26.92(Re2*) 0.773Np7.15(Re*) 0.775两轴各自的功率消耗都增加，在黏度小于 7 Pas 时，功率消耗与黏度呈线性关系，在黏度大于 7 Pas 时 功率消耗的增长速度有所减缓(2) 随着流速的增大，双轴所需的总功率及两轴 各

24、自的功率消耗都有减小 ，且与流速呈近似线性关系 (3) 随着转动角度的变化，轴所需功率消耗在周期范围内呈规律性变化(4) 在 1Re*50 时，清理轴、搅拌轴的功率准-数以及总的功率准数的对数与相对应的表观雷诺数但可以看出在此处功率准数与雷诺数并不成倒数关系，这一点与立式搅拌装置以及一些卧式搅拌装-的对数呈线性关系，关联式为：Np1=2.19( Re1*) 0.775，Np2 6.92(Re2*) 0.773 ，Np 7.15(Re*) 0.775 ，但在 此处功率准数与雷诺数并不成倒数关系-图 14 清理轴功率准数与表观雷诺数关系Fig.14 Power number of clean up

25、 shaft vs apparentReynolds number符号说明T力矩，Nm；1008天津大学学报第 44 卷 第 11 期F力，N；l 距离，m； N转速，rad/s； P功率，W； d搅拌桨的直径，m； L釜长，m； 物料密度，kg/m3；物料的黏度，Pas； Re*表观雷诺数，Re*d 2N/； Np功率准数，NpP/N 3d 4L6焦海亮，包雨云，黄雄斌，等. 高黏度流体混合研究进展J. 化工进展，2007，26(11)：1574-1582.Jiao Hailiang ， Bao Yuyun ， Huang Xiongbin ， et al. Recent research

26、progress in mixing of high viscous fluidJ. Chemical Industry and Engineering Progress ，2007，26(11)：1574-1582(in Chinese).赵 轶，冯连芳，顾雪萍，等. 卧式双轴 T 型搅拌器 在牛顿流体中的功率消耗J. 化学工程，2000，28 (4)：32-35.Zhao Yi，Feng Lianfang，Gu Xueping，et al. Power consumption of horizontal two-shaft t-shape self-cleaning agitator in

27、 newtonian fluid J. Chemical Engineering，2000，28(4)：32-35(in Chinese).Letellier B，Xuere C. Scale-up in laminar and transient regimes of a multi-stage stirrer，a CFD approachJ.Chem Eng Sci，2002，57(21)：4617-4632.Kelly W，Gigas B. Using CFD to predict the behaviour of power law uid near axial-ow impeller

28、s operating inthe transitional ow regimeJ. Chem Eng Sci，2003，58(10)：2141-2152.7参考文献：1冯惠生. 高凝固点热敏物系精馏节能新技术J. 化工进展，2006，25(2)：6-9.Feng Huisheng. Energy-saving distillation technology on the heat-sensitive complicated multi-components mixture with high melting-point and high-viscidity J. ChemicalIndust

29、ry and Engineering Progress，2006，25(2)：6-9(in Chinese).冯惠生，李小平，姜浩锡，等. 邻对苯二酚产品分离 纯化过程的热敏性问题J. 化学工业与工程，2003，20(1)：45-50.Feng Huisheng ， LI Xiaoping ， Zhang Haoxi ， et al. Problem of heat-sensitive materials during separation and purification process for producing pyrocatechol and hy-droquinoneJ. Chemi

30、cal Industry and Engineering ，2003，20(1)：45-50(in Chinese).冯惠生，魏月娥，陈毅烈，等. 真空精馏提纯顺酐的 研究J. 精细石油化工，2007，24(2)：50-52.Feng Huisheng，Wei Yuee，Chen Yilie，et al. Purifica- tion of maleic anhydride by vacuum distillation J. Spe- ciality Petrochemicals ， 2007 ， 24(2) ： 50-52(in Chi-nese).冯连芳，赵 铁，顾雪萍，等. 卧式双轴 T

31、 型搅拌 器在非牛顿流体中的搅拌功率特性J. 化学反应工程 与工艺，2000，16(4)：331-336.Feng Lianfang，Zhao Yi，Gu Xueping，et al. Power consumption of horizontal two-shaft t-shape self-cleaningagitators in the non-Newtonian fluid J. Chemical Reac- tion Engineering and Technology，2000，16(4)：331-336(in Chinese).冯连芳，曹松峰，顾雪萍，等. 高黏搅拌聚合反应装 置J合成橡胶工业，2001，24(5)：257-261.89210 冯惠生，刘 荣，杨 腾. 应用 CFD 对 D-T 型双轴卧式搅拌装置功率特性和传热特性的模拟J. 化工进 展，2009，28(增 2)：318-322.Feng Huisheng，Liu Rong，Yang Teng.