旋流分离器液固分离数值模拟研究

1、旋流分离器液固分离数值模拟研究杨博文郑小涛/rj瀚浩文王浩宇黄苏武汉工程大学机电工程学院摘要：旋流分离器被发明至今己有一百多年，是一种古老而现代的多相流分离技术。被 广泛用于各个行业。水力旋流器的内部进行着复杂的两相流运动，是液相以及固 体颗粒相分离的一个过程，而固体颗粒的在旋流器内部的运动，很大程度上由 液相的运动所决定。本文通过采用cfd软件fluent对旋流分离器进行模拟研究, 液体的连续相采用rng k- e模型，i占i体离散相采用dpm模型，通过数值分析可 以得到旋流器内部流场的情况，探究进口速度以及颗粒大小对分离效率的影响。关键词：旋流分离器;分离效率;粒径大小;入口速度;模拟;旋

2、流分离设备lq被广泛应用于石油化工、采矿等行业2-3,自从e. bretney 于1891年申请了首个水利旋流分离器的专利后，旋流分离技术的研究开始了飞 速的发展。分离效率是评价一个旋流分离器好坏的一个很重要的参数，影响分离 效率的因素有许多，主要分为三类:操作参数、结构参数和物性参数。本文利用 cfd4,5软件fluent6, 7探究混合物进口速度以及分离颗粒的粒径对旋流分 离器效率的影响。1旋流分离器的工作原理1. 1基本原理图1旋流分离器内部流动情况下载原图离心沉降作用是旋流分离设备的主要工作原理，在一定的压强作用之下，由入 口进入的混合流体（常见的有固-液混合物或油水混合物等）沿着切向

3、方向进入 到旋流分离器的内腔z后，由于高质相和低质相存在着密度差，所以受到不同 的的向心浮力、离心力以及流体阻力，所以在旋流分离器内壁的限制下，部分流 体一边向下一边向中心做着旋转运动，这些便形成了外旋流，由于旋流分离器 内部中心低压区的吸引以及分离器圆锥段的限制，内部的流体变为向上做着旋 转运动，这些便形成了内旋流，如图1所示。随着流动半径的减小，压力能逐渐转变为流体的动能和其他的能量损失。高质相 在径向主要受到:离心惯性力，向心浮力和流体阻力8-10,这3个不同的力促 使不同密度不同粒径的微粒在旋流分离器的径向上有序的分布，并最终随着外 旋流从底部被排出，或随内旋流从溢流口排出。1. 2

4、rng k- £湍流模型yakhot和orszag等人利用重整化群理论发展并改进而成的rng k- e模型，在 k-£模型的基础上，增加了一个附加生成项，从而在一定程度上弥补了标准 k-e模型上的一些不足。rng k- s模型考虑了有效雷诺数对湍流输运的影响，通过对旋转流动的研究， 使rngk- £模型适合处理低雷诺数、旋转度人的流动以及高应变率的流体。所以, 本文中将采用该模型进行数值计算。2数值模拟2.1计算模型本次计算采用的主要几何参数如下:d=275mm, s二50mm, b=32mm, de二68mm, hl=200mm, h2=530mm,入口段内径

5、50mm,溢流管长度 160mmo将分离器分为四个区域:入口段、圆柱段、圆锥段、溢流管。旋流器的结构图如 图2所示。图2旋流分离器模型下载原图 2.2网格的划分网格的生成为数值模拟前处理的一个过程。模型的网格划分是进行数值计算非常 关键的一步，网格质量的好坏，不仅仅会影响到计算结果的准确度，而且会让 计算结果岀现较大的错误。针对旋流分离器内部的流场特点，此次计算将采用ansys自带的icem cfd对旋 流器网格进行划分。经过大量模型计算，比较后证明：计算网格划分得太稀疏， 不能准确反应流场的特征，误差大；网格数目达到一定程度时，己经可以较好地 反映岀真实的情况，继续增大网格数目，不再能加强计

6、算的精度，却增加了计 算消耗的资源，降低了计算的速度。木章模型将整个计算域划分为53万左右的 网格数目，包含156万个面，54万多个节点，可以得到很好的计算结果。网格 划分如图3所示图3分离器网格划分下载原图2.3计算的材料以及边界条件材料参数如表1,所说的边界条件是在求解域的边界上所求解的变量或其一阶 导数随地点及时间变化的规律。(1) 入口：速度进口 (velocity inlet),大小为 20. 00m/s;(2) 出口：自由出流(outflow)，流量为1;(3) 颗粒捕集口：自由出流(outflow),流量为0;(4) 本次模拟的离散相为沙粒，密度为1500kg/mo表1材料参数下

7、载原表3模拟结果3. 1单向流结果3. 1. 1单向流速度云图由于本次采用的单侧入口，导致了旋流器内部的漩涡核心并不在轴心。从图4 中可得到，圆柱段的速度对称性差一些，而圆锥段的速度对称性较好，等速度 线基本上是圆环。从全局来看，分离器内部的液体旋流是比较对称的，由速度分 布的云图颜色可以看出，上部分的速度大小的对称性略差，入口一侧的速度较 大;而越趋向下部分，速度大小分布的对称性越好。图5为纵截面上的速度矢量图，从图中可以看岀，水在进入旋流分离器之后， 开始向下流动，形成外旋流，到达分离器底部之后再向上流动，形成内旋流。在 主流上伴有许多局部二次涡流：(1)进入旋流器后，部分液体沿壁面向上运

8、动， 在顶部又反向沿着溢流口段的外壁向下流动，形成环形空间内的纵向涡流，这 是最主要的二次涡流，这种涡流会把一部分颗粒向上带到顶部形成“上灰环” 导致无法分离；（2）溢流管段下端附近有较大的速度，紊流严重。由于径向速度 较大，小部分流体从进口进入后直接进入溢流管，形成短路流；（3）在圆锥段 也有一些不规则的小的局部二次涡流。图4速度分布图下载原图图5纵截面速度矢量图下载原图 3.1.2单向流速度分布图6纵截面切向速度图下载原图图6为纵截面不同高度处的切向速度图，切向速度分布的轴对称性很好，其分 布曲线呈“驼峰”形分布，就是每侧都存在一个最大切向速度点。屮心区域旋流 运动产生的离心力有利于将颗粒

9、甩向外部;而外部由于旋流强度低，对颗粒的携 带作用减弱（尤其是一些大粒径的颗粒），这样利于颗粒在壁面被捕集。下图7为纵截面不同高度处的轴向速度分布图，轴向速度的分布对称性较好。可 将液体分为下行外旋流和上行内旋流两区域。上、下旋流的分界点与旋风分离器 的形状有关，在圆柱段部分，此分界面大致也呈圆柱状，其半径一般要稍大于 溢流管的半径。下载原图图7纵截面轴向速度分布图3. 2固液分离模拟结果 旋流分离器的单相液相流场已经通过计算得到。但旋流分离器实际进行的是液、 固两相的分离，因而必须引入相应的两相流动模型，对其中的颗粒相研究。本次 模拟采用dpm模型进行模拟。设置跟踪的步长为100000步，有

10、足够的吋间让所有可以离开分离器的颗粒被分 离。颗粒在所设置的时间或步长内离开了分离器，便不再跟踪;若在所设定的时 间或步长内没有离开分离器，便认为分离器不能对该颗粒进行分离，也不再跟 踪。这样，分离器的总分离效率为：通过整理计算之后的结果可以得到表2各种粒径的分离效率以及图8各种粒径分 离效率柱状图。表2各种粒径的分离效率下载原表图8各种粒径分离效率柱状图下载原图通过对比可以得到：仃)当入口流速增大时，由于离心力场增强，提高了对沙粒的分离效率。随着入 口速度的加大，旋流器的分离效率也跟着增加，尤其是小粒径和中等粒径颗粒 分离效率的增幅更大，大颗粒的增幅较小。所以，入口速度的提高有利于颗粒的 捕

11、集，分离器分离效率也越高。然而，入口流量不能太大，随着流量的增加，压 力的损失也相应增加。另外，过高的流速会加大颗粒的破碎以及分离器内壁面的 磨损。(2) 相同的流量下，颗粒粒径越大，其分离效率越大，但增加的幅度越来越小， 最后趋向于稳定。从理论上说，旋流分离器都有一个固定的临界粒径，小于临界 粒径的微粒是无法被捕获的，但实际上，颗粒进入了分离器z后，因为存在颗 粒间的互相碰撞，颗粒的聚集夹带及静电和分子引力等因素，颗粒的运动有较 大的随机性，使得一部分小于临界粒径的颗粒也可以被捕集，同样的，一部分 大于临界粒径的大颗粒也可能会逃逸。参考文献ljchine b, concha f. flow

12、patterns in conical and cylindrical hydro-cyclonesjchemical engineering journal, 2000, 80:267-273.2 辛周，龚俊，黄建龙，等旋流分离技术在黄河泥沙分离中的应用j.兰州: 兰州理工大学学报，2004, 30 (4) : 67-69.3 王福军计算流体动力学分析:cfd软件原理与应用m.北京:清华大学出版 社，2004, 30 (4) : 67-69.4 cullivan j c, willaams r a. new understanding of a hydrocyclone flow field

13、 and separation mechanism from computer fluid dynamics j. minerals engineering, 2004, 17:650-660.5 徐继润罗茜水力旋流器流场理论m北京:科学出版社，1998.6 韩占忠，王敬，兰小平.fluent流体工程仿真计算实例与应用m 北京:北 京理工大学出版社，2004.7 a w b, j n t. mineral processing technologym.amsterdam:elsevier science&technology books, 2006:212-223.8 chu k w, kuang sb, yu a b, et al. peirticle scale modelling of the multi phase flow in a dense medium cyclonc:effcct of fluctuation of solids flowratej. minerals engi