《武汉工程大学学报》 2015年01期
22-24
出版日期：2015-01-31
ISSN:1674-2869
CN:42-1779/TQ
圆柱段长度对水力旋流分离性能的影响
０ 引 言 水力旋流分离器广泛应用于化工、石油、选矿、环保、制药、食品、轻工等领域，其中固－液旋流器发展较早，技术也相对成熟［１］. 液－液旋流器在最近１０多年才用作于海上油气平台的油－水分离装置［２］. 利用油－水等两相互不相溶物质的密度差和水力旋流器内部强烈的回旋流产生的离心力，使油水发生分离. Ｔｈｅｗ［３］和Ｃａｌｄｅｎｔｙ［４］研究了液－液分离过程和液－固分离的不同之处，发现由于液体之间的密度差较小，液－液分离较固－液分离困难，因此对液－液旋流器的结构研究显得格外重要. 近年来，许多学者对分离器的结构进行了研究，通过优化入口形式［５］、锥角［２］以及溢流口［６］等方式提高油水旋流器的分离效率，但对圆柱段结构的研究相对较少. 本文将以圆柱段长度为出发点，通过数值模拟的方法，研究圆柱段长度对旋流分离器的性能影响. １ 物理模型与计算模型 １．１ 物理模型 图１为油－水旋流器的物理模型，基准模型采用的是Ｍａｒｔｉｎ Ｔｈｅｗ的液－液旋流器. 模型为对称双入口形式，因此具有良好的流场稳定性和对称性. 其基准尺寸为大小锥结合处直径Ｄ，大锥段锥角α为２０°，小锥段锥角β为１．５°. 圆柱段长度Ｌｃ的标准长度为２Ｄ，本文基于改变Ｌｃ长度进而分析油－水旋流分离的分离性能. 图１ Ｍａｒｔｉｎ Ｔｈｅｗ液－液旋流器模型图 Ｆｉｇ．１ Ｔｈｅ ｍｏｄｅｌ ｏｆ Ｍａｒｔｉｎ Ｔｈｅｗ’ｓ ｌｉｑｕｉｄ－ｌｉｑｕｉｄ ｈｙｄｒｏｃｙｃｌｏｎｅ １．２ 计算模型 由于旋流器的非稳态不可压缩流动特征，因此采用适于流线强烈弯曲和旋转三维流动的雷诺应力模型（ＲＳＭ）［７］. ■（?籽u′iu′j）＋■（?籽u′iu′j）＝Dｉｊ＋Pｉｊ＋Πｉｊ＋?着ｉｊ（１） 式（１）中Ｄｉｊ、Ｐｉｊ、Πｉｊ、εｉｊ分别为扩散项、产生项、压力应变项和耗散项，分别表示如下： Dｉｊ=■（■■）（２） Pij=-(u′iu′K■)+u′iu′K■（３） Πｉｊ=-C1?着aij+C2?着(aikakj-■akjaki?啄ij)+C3kSij+C4k(aikSjk+ajkSjk+ajkSjk-■akjaki?啄ij)+C5(aikWjk+ajkWik)（４） 压力应变项包括了雷诺应力的各向异性张量的二次方项，式（４）中： C1=3.4+1.8PKK/?着; C2=4.2; C3=■-1.3Π■■; C4=1.35; C5=0.4; aij=■-■?啄ij; Πa=aijaji; Sij=■(■+■); Wij=■(■-■) 多相流模型选择的是适用于强烈耦合的各向同性多相流模型Ｍｉｘｔｕｒｅ模型［８］. 对流项离散格式选择的是在复杂网格条件下具有良好收敛性的二阶迎风格式，压力插补格式采用的是ＰＲＥＳＴＯ格式，压力－速度耦合算法采用的Ｓｉｍｐｌｅ算法. １．３ 边界条件与数据处理 入口采用速度入口，入口含油率为５％，油的密度为８５０ ｋｇ／ｍ３，黏度为３．３２ ｃＰ；溢流口和底流口均为充分发展ｏｕｔｆｌｏｗ，其流量比分别为１５％和８５％；壁面采用无滑移条件. 分离效率：Ｅ＝１-?琢ｕ／?琢ｉｎ 其中；?琢ｕ为底流口含油溶度；?琢ｉｎ为入口含油浓度. 底流口压降：ΔＰｕ＝Ｐｉ - Ｐｕ 溢流口压降：ΔＰｏ＝Ｐｉ-Ｐｏ 其中：Ｐｉ为入口压力；Ｐｕ为底流口压力；Ｐｏ为溢流口压力. ２ 数值模拟计算 ２．１ 网格划分与独立性验算 网格采用ｇａｍｂｉｔ２．４．６进行划分，由于旋流器结构细长，且不同部位结构尺寸差距较大，因此采用分区划分. 以溢流口、入口段、圆柱段和大锥段、小锥段和底流管四个区域，采用六面体非结构网格，图２为网格划分整体图. 图２ 网格划分图 Ｆｉｇ．２ Ｇｅｎｅｒａｔｅｄ ｍｅｓｈ ｆｏｒ ｈｙｄｒｏｃｙｃｌｏｎｅ ｓｉｍｕｌａｔｉｏｎ 为保证计算结果准确性，对网格进行独立性验证，验证模型初始条件采用Ｂｅｌａｉｄｉ ａｎｄ Ｔｈｅｗ［９－１０］的实验数据和结果，分别在３种不同密度网格下，不同粒径油滴的分离效率与实验数据作对比. 图３为网格独立性验证图，由图3可知，在网格数量为５８万时，计算结果与实验数据较接近，而当网格数量进一步加大到９５万时，分离效率变化不大，且误差控制在１０％以内，计算结果有效. 图３ 网格独立性验证 Ｆｉｇ．３ Ｔｈｅ ｉｎｄｅｐｅｎｄｅｎｃｅ ｖｅｒｉｆｙ ｏｆ ｇｒｉｄ ２．２ 数值模拟结果 分别将圆柱段长度Ｌｃ为２Ｄ、１．５Ｄ、Ｄ、０．５Ｄ、０五种油水旋流器的模型导入ＦＬＵＥＮＴ中，以油滴粒径为３０ ?滋ｍ和入口速度分别为４ ｍ/ｓ、５ ｍ/ｓ、６ ｍ/ｓ进行计算，最终得到不同圆柱段长度油水旋流器的分离效率和压力降. 图４为不同圆柱段长度下的分离效率对比图，结果表明，分离效率随圆柱段缩短呈线性增长，在Ｌｃ＝０．５Ｄ时达到最大. 当Ｌｃ进一步减小到０，即无圆柱段时，分离效率基本无较大变化，并在入口速度为４ ｍ/ｓ时呈下降趋势. 图4 不同圆柱段长度分离效率对比图 Fig.4 De-oiling hydrocyclone separation efficiency of different length of cylindrical section 图５和图６分别为不同圆柱段长度下底流口压降和溢流口压降对比图，结果表明，随着圆柱段长度变化，溢流口和底流口的压力降变化趋势基本一致. 当圆柱段长度减小，压力降增大；当去除圆柱段时，压力降有所增加，但幅度较小. 图5 不同圆柱段长度底流口压降对比图 Fig.5 The underflow mouth pressure drop of different length of cylindrical section 图６ 不同圆柱段长度溢流口压降对比图 Ｆｉｇ．６ Ｔｈｅ ｏｕｔｆｌｏｗ ｍｏｕｔｈ ｐｒｅｓｓｕｒｅ ｄｒｏｐ ｏｆ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ３ 结果分析 ３．１ 速度分布 旋流器以Ｙ轴为中心轴，大锥段与圆柱段结合处为Ｙ＝０截面. 图７和图８分别为Ｙ＝－３０ ｍｍ和Ｙ＝－８０ ｍｍ（大锥段中间和小锥段上端）截面处的不同圆柱段长度的切向速度分布图. 由图７和图８可知，不同圆柱段长度旋流器的近壁面和内旋流区域的流体切向速度相差不大，而在外旋流区域，切向速度随着圆柱段的减少而增大，在主分离区域提供较大的离心力，增加油水旋流器的分离效率. 图９和图１０分别为Ｙ＝-３０ ｍｍ和Ｙ＝-８０ ｍｍ截面处的轴向速度分布图，由图9和图10可知，随圆柱段长度减小，外旋流区域的轴向速度变化不大，而内旋流区域轴向速度随圆柱段长度的减小而增大. 图７ Ｙ＝-３０ ｍｍ平面处切向速度对比图 Ｆｉｇ．７ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｔａｎｇｅｎｔｉａｌ ｖｅｌｏｃｉｔｙ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｃｅ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ｉｎ Ｙ＝－３０ ｍｍ 图８ Ｙ＝ -８０ ｍｍ平面处切向速度对比图 Ｆｉｇ．８ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｔａｎｇｅｎｔｉａｌ ｖｅｌｏｃｉｔｙ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ｉｎ Ｙ＝－８０ ｍｍ 图９ Ｙ＝－３０ ｍｍ平面处轴向速度对比图 Ｆｉｇ．９ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ａｘｉａｌ ｖｅｌｏｃｉｔｙ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ｉｎ Ｙ＝－３０ ｍｍ 图１０ Ｙ＝－３０ ｍｍ平面处轴向速度对比图 Ｆｉｇ．１０ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ａｘｉａｌ ｖｅｌｏｃｉｔｙ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ｉｎ Ｙ＝－８０ ｍｍ ３．２ 流场分布分析 图１１为不同圆柱段长度下的流场分布图，随圆柱段长度的减小，旋流器圆柱段和大锥段区域的循环流区域逐渐减小，且涡流核心向下迁移，在Ｌｃ＝０时，循环流消失. 循环流会影响油滴向内旋流迁移，并随循环流区域的增大，旋流器的分离效率会降低［１１］. 因此，随圆柱段长度的减小，其油水旋流器的分离效率增大. 图１１ 不同圆柱段长度下流场对比图 Ｆｉｇ．１１ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｆｌｏｗ ｆｉｅｌｄ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ３．３ 湍流分析 图１２为５种圆柱段长度下的湍流强度分布云图，由图12可知，因圆柱段和大锥段的中心区域循环流的改善，随圆柱段长度的减小，其湍流强度减小. 而在大小锥结合处，湍流强度随着Ｌｃ的减小而增大，且此区域的湍流强度达到最大. 图１２ 不同圆柱段长度湍流强度分布图 Ｆｉｇ．１２ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｔｕｒｂｕｌｅｎｃｅ ｉｎｔｅｎｓｉｔｙ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ 图１３为大小锥结合处的湍动能耗散率分布图，前文所知，随圆柱段减小，底流口和溢流口的压力降反而上升的结论和一般低雷诺数下的圆管流动规律相背离. 然而从图１３可看出随着圆柱段长度减小，其湍动能耗散率增大. 由湍流强度分布和湍动能耗散率的分布可解释随圆柱段长度的减小，其底流口和溢流口的压降增加的现象. 图１３ 大小锥结合处湍动能耗散分布 Ｆｉｇ．１３ Ｔｈｅ ｔｕｒｂｕｌｅｎｔ ｄｉｓｓｉｐａｔｉｏｎ ｒａｔｅ ｉｎ ｔｈｅ ｃｏｎｎｅｃｔｉｏｎ ｏｆ ｔｈｅ ｂｉｇ ａｎｄ ｓｍａｌｌ ｃｏｎｅ ｓｅｃｔｉｏｎｓ ｏｆ ｄｉｆｆｅｒｅｎｔ ｌｅｎｇｔｈ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｓｅｃｔｉｏｎ ４ 结 语 通过对不同长度的油水旋流器进行数值模拟研究，结果表明： ａ． 随着圆柱段长度的减小，其分离效率呈线性增长，但去除圆柱段后分离效率无变化，并呈减小趋势；底流口和溢流口的压降随着Ｌｃ的减小而增大. ｂ． 随着圆柱段长度减小，在其大锥段和小锥段处的切向速度和轴向速度逐渐增大，循环流区域逐渐减小，并在去除圆柱段后完全消失，从而提高了分离效率. ｃ． 圆柱段减小后，其大小锥结合处湍流强度增大，且其湍动能耗散率也随着Ｌｃ的较小而增大，解释了随着圆柱段减小其底流口和溢流口的压降增大的现象. 致 谢 本文研究工作得到了湖北省自然科学基金项目（２０１２ＦＦＢ０４７０７）、武汉工程大学科学研究基金项目（Ｋ２０１４１４）和武汉工程大学研究生教育创新基金项目（ＣＸ２０１３０８０）的资助，在此一并表示衷心感谢.