《武汉工程大学学报》 2014年03期
48-53
出版日期：2014-03-31
ISSN:1674-2869
CN:42-1779/TQ
液液旋流分离器分离特性数值模拟
0引言水力旋流分离器（Hydrocyclone）是一种利用两相或多相之间的密度差，在高速旋转的流体中产生不同的离心力而分离的设备.水力旋流分离技术是一种较为成熟的分离技术，广泛应用于石油、化工、生物制药、煤矿、环保等行业［1］.液液旋流分离器是水力旋流分离器的一种，其分离介质为两相或多相不相溶的液体.液液旋流分离器研究起步较晚，并且通常液液两相的密度差较小，因此分离较固液或气液旋流分离器困难.由于液液旋流分离器复杂的结构和高速旋转的流体运动，其内部流场通常较为复杂，并且现有实验设备难以得到清晰的液液两相分离过程.本文将采用CFD（计算流体动力学）技术，对液液旋流分离器进行数值模拟分析，得到其流场分布和液液两相分离过程.1分析模型鉴于目前液液旋流分离器的应用情况和研究进展，选用实验数据较为权威的Martin Thew［2］型液液旋流分离器.图1为液液旋流分离器的模型图，基本由圆柱段、大锥段、小锥段、底流管、溢流口和入口组成，其入口为对称切向双入口，各部分的结构尺寸如表1所示，其基准尺寸Ds为20 mm，底流管长度Lu为235 mm.图1Martin Thew型液液旋流分离器模型Fig.1The model of liquidliquid Hydrocyclone by Martin Thew表1计算模拟尺寸表Table 1The size of modelDs/mmD/mmDo/mmDi/mmDu/mmLc/mmα/(°)β/(°)2040471040201.52数值模拟方法2.1网划划分及独立性验证采用FLUENT软件进行数值模拟，Gambit2.4.6软件进行建模和网格划分.综合国内外学者对旋流分离的数值模拟的研究，采用六面体非结构网格对模型进行划分，网格划分图如图2所示. 图2模型网格划分图Fig.2The grids of model为了保证数值模拟计算的准确性，将对网格进行独立性验证，将不同网格密度的分离器计算结果与Martin Thew［3］的实验数据进行对比，对比结果如图3所示.由图3可知，当网格数为58万时，其计算分离效率与实验数据相差为10%以内，而当网格进一步加密到95万时，计算结果与网格数为58万相差不大，为了提高计算效率，本文将采用网格数为58万的模型进行分析. 图3网格独立性验证Fig.3The independence verify of grids 注：Experiment;mesh density(950778);mesh density(580876);mesh density(313848)2.2计算模型由于旋流分离器的非稳态不可压缩流动特征，因此湍流模型采用适于流线强烈弯曲和旋转三维流动的雷诺应力模型(RSM)［4］.多相流模型选择的是适用于强烈耦合的各向同性多相流模型Mixture模型［5］.对流项离散格式选择的是在复杂网格条件下具有良好收敛性的二阶迎风格式，压力插补格式采用的是PRESTO［6］格式，压力速度耦合算法采用的Simple算法［7］.2.3介质物性及边界条件入口采用速度入口，速度为6 m/s，入口含油率为5%，油的密度为850 kg/m3,黏度为3.32 cP；溢流口和底流口均为充分发展出口边界条件，其流量比分别为15%和85%；壁面采用无滑移条件［8］.分离效率：E=1-αu/αin 其中；αu为底流口含油溶度；αin为入口含油浓度.第3期郑小涛，等：液液旋流分离器分离特性数值模拟武汉工程大学学报第36卷3计算结果分析 将网格导入FLUENT软件中，设计边界条件和计算控制模型，得到液液旋流分离器数值模拟计算结果.3.1流场分析图4为圆柱段和大锥段处的流场分布图，从图4可以看出：液液旋流分离器内部流场中向底流口流动的外旋流，向溢流口流动的内旋流分布明显；在圆柱段顶部及器壁处，由于流体的粘性作用存在短路流；在圆柱段，有大量的循环流存在，循环流的存在会影响油滴向中心的迁移，降低分离效率. 图4圆柱段及大锥段处流场分布图Fig.4The flow field distribution of cylindrical section and big cone section图5为大小锥结合处附近的流场分布图，从流体的流线可知，在小锥段上端依然存在明显向上的内旋流，并且由外旋流向内旋流迁移的流线非常明显，证明小锥段上端依然是油水的分离区域，并对油水的分离有一定作用.在大小锥结合处也存在明显的循环流，说明无论是圆柱段、大锥段还是小锥段，都存在着循环流. 图5大小锥结合处流场分布图Fig.5The flow field distribution of the juncture of small cone section with big small cone section3.2速度分布分析为分析液液旋流分离器内部流场分布，将在分离器各部分截取截面，从其各截面速度分布角度分析其流场及分离机理，截面截取示意图如图6所示. 图6截面及坐标系定位Fig.6The localization of coordinates and sectiion图7为各截面的切向速度分布图，切向速度是提供离心力的主要动力源.由图7可知，各截面切向速度均呈组合涡［9］特征，在外旋流为自由涡，内旋流为强制涡，并且在圆柱段和大锥段组合涡特征更为明显，与理论分析［10］相一致.从各截面切向速度数值比较来看，切向速度沿着中心轴向下依次减小，但是在内旋流区域下降梯度相等，这也是液液旋流分离器结构设计的精妙所在. 图7各截面切向速度分布图Fig.7The tangential velocity of different sections图8为各截面轴向速度分布，由于液液旋流分离有两个出口，轴向速度分布能够反映流体向各个出口的流动情况.由图8可知，圆柱段和大锥段的轴向速度分布呈3波峰形式，在外旋流区域存在波峰，此现象是由于圆柱段和大锥段的循环流导致的流体波动所致.从数值比较进行分析，在圆柱段和大锥段内旋流区域向上的轴向速度较大，所以此区域中心处流体向溢流口流动，而在小锥段上端轴向速度依然大于零，说明小锥段上端也为油水分离区域. 图8各截面轴向速度分布Fig.8The axial velocity of different sections3.3基于DPM模型油滴分离轨迹分析DPM模型是一种粒子随机轨道模型，即将油滴看作为固体颗粒，由于本文所分析的初始条件中含油率为5%，因此可以忽略油滴对流场的影响.图9为分离成功时油滴的运动轨迹.由图9可知油滴进入液液旋流分离器后，在圆柱段形成稳定的旋流轨迹，并向底流口运动，在外旋流区域随着流体做旋转运动，并且旋转半径不断减小，油滴颗粒到达小锥段上端时，进入沿溢流口流动的内旋流，并随着外旋流做旋转运动最后由溢流口排出. 图9分离成功油滴运动轨迹Fig.9The trail of the successfully separated oil droplets图10为分离失败时油滴的运动轨迹，与图9相比，其不同在于油滴在达到小锥段上端后并没有进入向溢流口流动的内旋流，而是随着流体向底流口排出. 图10分离失败油滴运动轨迹Fig.10The trail of the unsuccessfully separated oil droplets4结语a.液液旋流分离器内部流场分布特征明显，在圆柱段、大锥段和小锥段处存在大量的由循环流导致的漩涡.b.针对旋流分离器各截面的速度分布研究，发现其切向速度和轴向速度分布规律明显,并与理论分析相一致.切向速度分布中在圆柱段、大锥段和小锥段上端组合涡特征明显，轴向速度在圆柱段和大锥段由于循环流的存在呈现双Ｗ形式.c.采用DPM模型对油滴分离轨迹进行追踪，进一步验证小锥段上端也为分离区域，油滴是否更够分离出来，其关键在于油滴是否能够在分离区域内进入向上的内旋流.致谢本文研究工作得到湖北省自然科学基金项目（2012FFB04707）、武汉工程大学研究生教育创新基金项目（CX201226,CX2013080）的资助，在此一并表示衷心的感谢.