文章快速检索     高级检索
  浙江大学学报(工学版)  2018, Vol. 52 Issue (9): 1811-1820  DOI:10.3785/j.issn.1008-973X.2018.09.023
0

引用本文 [复制中英文]

董辉, 伍开松, 况雨春, 代茂林. 基于DEM-CFD水力旋流器的水合物浆体分离规律研究[J]. 浙江大学学报(工学版), 2018, 52(9): 1811-1820.
dx.doi.org/10.3785/j.issn.1008-973X.2018.09.023
[复制中文]
DONG Hui, WU Kai-song, KUANG Yu-chun, DAI Mao-lin. Study on separation law of hydrate slurry in hydrocyclone based on CFD-DEM[J]. Journal of Zhejiang University(Engineering Science), 2018, 52(9): 1811-1820.
dx.doi.org/10.3785/j.issn.1008-973X.2018.09.023
[复制英文]

基金项目

国家自然科学基金资助项目(L1322021);中国工程院自然基金资助项目(2013-04-ZQC-002)

作者简介

董辉(1991—),女,硕士,从事水合物浆体多相分离研究. orcid.org/ 0000-0002-6691-6463.
E-mail: 18280358469@163.com.

通信联系人

伍开松,男,教授. orcid.org/ 0000-0002-1092-1247.
E-mail: wks9998@163.com
.

文章历史

收稿日期:2017-06-23
基于DEM-CFD水力旋流器的水合物浆体分离规律研究
董辉, 伍开松, 况雨春, 代茂林     
西南石油大学 机电工程学院,四川 成都 610500
摘要: 采用离散元素法和计算流体力学耦合仿真的方法,研究不同入口速度下的水力旋流器分离效率. 以较佳入口速度的分离效率工况为例,研究不同入口区域的水合物颗粒的力学行为与特征. 结果表明:由入口上侧和入口内侧进入的水合物颗粒具有短路流特征. 由入口下侧和入口中心进入的水合物颗粒以螺旋轨迹线向下运动. 同时在离心浮力与离心惯性力的差值作用下,沿半径向中心运动. 进入内旋流场后,水合物颗粒则改变方向,以螺旋轨迹向上运动,直至从溢流管被排出. 这是典型的水合物颗粒分离运动轨迹. 由入口外侧进入的水合物颗粒中大多数也具有这一特征,只有少部分水合物颗粒紧贴水力旋流器的内壁,沿边界层或在边界层附近,以螺旋轨迹一直向下运动,最终从底流口排出. 这一部分水合物颗粒对总体分离效率影响最大.
关键词: 水力旋流器    水合物浆体    离散元法    计算流体力学    颗粒运动    
Study on separation law of hydrate slurry in hydrocyclone based on CFD-DEM
DONG Hui , WU Kai-song , KUANG Yu-chun , DAI Mao-lin     
College of Mechanical Engineering, Southwest Petroleum University, Chengdu 610500, China
Abstract: The separation efficiency of hydrocyclone at different inlet velocities was investigated by using the discrete element method and computational fluid dynamics coupling simulation method. Taking the optimum separation efficiency as an example, the mechanical behavior and characteristics of hydrate particles at different entrance regions were studied. Results show that the hydrate particles entering from upside and inside of the inlet have short circuit flow characteristics. The hydrate particles entering from underside and center of the inlet move downward with a spiral path, which at the same time move toward the center along the radius under the action of the buoyancy and the centrifugal force. The hydrate particles change direction into the swirling and move upward with the spiral path until discharged from the overflow pipe. This is the typical separation trajectory of hydrate particles. Most hydrate particles entering from outside of the inlet also has this feature. Only few of hydrate particles are close to the inner wall of hydrocyclone and move downward with a spiral path along the boundary layer, or near the boundary layer finally discharged from the underflow. This hydrate particle has the greatest influence on the overall separation efficiency.
Key words: hydrocyclone    hydrate slurry    discrete element method    computational fluid dynamics    particle motion    

在海洋非成岩水合物藏开发过程中,形成的水合物浆体是由水合物颗粒、砂颗粒和海水三相组成的. 在海底采用水力旋流器实现水合物浆体分离(即固-固-液分离)是一种比较经济的方式. 水力旋流器没有运动构件,具有结构简单、分离效率高、维护费用少和安装灵活等特点.

由于水力旋流器内部流场的强烈旋转湍流特征[1-3],导致被分离颗粒的运动也具有一定的随机性[4]. 固相颗粒在水力旋流场中的分离运动显得特别复杂[5]. 颗粒运动的主要影响因素如下:固相颗粒的初始状态[6],包括入口速度和入口位置等;水力旋流场对固相颗粒的驱动效应[7];固相颗粒与器壁间的碰撞与摩擦[8];固相颗粒间的相互碰撞与摩擦[9-11]. 随着计算机技术的发展,许多学者利用计算流体力学(computational fluid dynamics, CFD)数值模拟的方法对水力旋流器进行了研究. 但是CFD模拟是建立在连续性假设的基础上,即认为研究对象是由相互连接没有间隙的大量微团构成,把固体颗粒作为流体,没有考虑颗粒对流场的作用,也没有考虑颗粒和颗粒、颗粒与壁面的碰撞,没有真实体现颗粒的特性,且离散相模型(DPM, discrete phase model)忽略粒子间的碰撞,粒子体积分数不能太大. 离散元素法(discrete element method, DEM)在处理颗粒碰撞及颗粒与流体之间相互作用问题上具有较大的优势,通过与CFD方法相耦合,可以较真实模拟出颗粒和流体的运动情况.

本文利用DEM-CFD耦合的方法,研究水力旋流器内水合物浆体的分离问题. 重点研究固相沙粒和固相水合物颗粒在水力旋流器内的运动规律以及不同入口区域固相颗粒的运动轨迹和力学行为特征. 研究固相颗粒运动与水力旋流器分离效率之间的关系,为水合物浆体的分离和水力旋流器的性能研究提供参考.

1 固液仿真模型

在DEM-CFD耦合仿真中,颗粒相按DEM方法进行求解,液相由CFD进行求解,并通过固液两相的受力交换实现耦合.

1.1 液相控制方程

水力旋流器中的液相可视为黏性不可压缩的流体,采用CFD求解[12],液相的质量守恒方程为

$\frac{{\partial {\rho _{\rm{f}}}}}{{\partial t}} + \nabla \cdot ({\rho _{\rm{f}}}\;{u_{\rm{f}}}) = 0,$ (1)

动量守恒方程为

$\begin{split}\displaystyle\frac{{\partial ({\rho _{\rm{f}}}{u_{\rm{f}}})}}{{\partial t}} &+ \nabla \cdot ({\rho _{\rm{f}}}\;{u_{\rm{f}}}\;{u_{\rm{f}}}) =\\ & - \displaystyle\frac{{\partial p}}{{\partial x}} + \nabla \cdot ({\mu _{\rm{f}}}\nabla {u_{\rm{f}}})+ {\rho _{\rm{f}}}\;g + F,\end{split}$ (2)

能量守恒方程为

$E = {E_{\rm{i}}} + \frac{1}{2}\left({u^2} + {v^2} + {w^2}\right) + {E_{\rm{f}}}.$ (3)

式中: ${\rho _{\rm{f}}}$ 为液相密度, ${u_{\rm{f}}}$ 为液相流动的速度, ${\mu _{\rm{f}}}$ 为液相的动力黏度, $p$ 为作用在液相微元体上的压力, $F$ 为作用在流体微元体上的体积力, $E$ 为单位质量液相的内能, ${E_{\rm{i}}}$ 为液相内能, ${E_{\rm{f}}}$ 为液相势能.

1.2 颗粒相控制方程

泥砂和水合物作为水合物浆体中的固体介质,在水合物浆体的流动过程中将泥砂和水合物离散化,均视为颗粒材料,即采用离散元素法(DEM)求解[13]。离散元素法主要受到接触模型和牛顿第二定律的影响,接触模型是用于单元间接触力的计算,牛顿第二定律则是用于求解单元的位移、速度和加速度。

离散元素法的颗粒模型是将颗粒与颗粒、颗粒与边界的接触采用振动运动方程进行模拟. 如图1所示为离散元素法中将颗粒接触模型表示成振动模型[14].图中,FtanFnor分别表示颗粒所受外力的切向分量和法向分量,ktkn分别表示接触模型中的切向、法向弹性系数,ctcn分别表示接触模型中的切向、法向阻尼系数.

图 1 颗粒接触模型表示成振动模型 Fig. 1 Scheme of particles contact model expressed as vibration model

将颗粒接触过程的振动运动进行法向和切向分解,颗粒接触过程的法向振动方程[14]

${m_{12}}\frac{{{{\rm{d}}^2}{x_{\rm{n}}}}}{{{\rm{d}}{t^2}}} + {c_{\rm{n}}}\frac{{{\rm{d}}{x_{\rm{n}}}}}{{{\rm{d}}t}} + {k_{\rm{n}}}{x_{\rm{n}}} = {F_{\rm{n}}},$ (4)

切向振动方程为

${m_{12}}\frac{{{{\rm{d}}^2}{x_{\rm{t}}}}}{{{\rm{d}}{t^2}}} + {c_{\rm{t}}}\frac{{{\rm{d}}{x_{\rm{t}}}}}{{{\rm{d}}t}} + {k_{\rm{t}}}{x_{\rm{t}}} = {F_{\rm{t}}},$ (5)

颗粒的滚动方程为

${I_{12}}\frac{{{{\rm{d}}^2}\theta }}{{{\rm{d}}{t^2}}} + \left({c_{\rm{t}}}\frac{{{\rm{d}}{x_{\rm{t}}}}}{{{\rm{d}}t}} + {k_{\rm{t}}}{x_{\rm{t}}}\right)s = M.$ (6)

式中: ${m_{12}}$ 为颗粒1、2的等效质量, ${I_{12}}$ 为颗粒1、2的等效转动惯量, $s$ 为旋转半径, ${x_{\rm{n}}}$ ${x_{\rm{t}}}$ 分别为颗粒的法向和切向相对位移, $\theta $ 为颗粒自身的旋转角度, ${F_{\rm{n}}}$ ${F_{\rm{t}}}$ 分别为颗粒所受外力的法向分量和切向分量, $M$ 为颗粒所受外力矩, ${k_{\rm{n}}}$ ${k_{\rm{t}}}$ 为接触模型中的法向、切向弹性系数, ${c_{\rm{n}}}$ ${c_{\rm{t}}}$ 为接触模型中的法向、切向阻尼系数.

本文假定颗粒受到重力、颗粒-颗粒接触力、颗粒-壁面接触力、颗粒-流体相互作用力等. 颗粒在t时刻的控制方程为

$m\;\frac{{{\partial ^2}u}}{{\partial {t^2}}} = mg + {F_{{\rm{p}} \text{-} {\rm{p}}}} + {F_{{\rm{p}} \text{-} {\rm{w}}}} + {F_{{\rm{p}} \text{-} {\rm{f}}}},$ (7)
$I\;\frac{{\partial \omega }}{{\partial t}} = T.$ (8)

式中: $m$ $u$ 分别为颗粒的质量和速度, ${F_{{\rm{p}} \text{-} {\rm{p}}}}$ 为颗粒间接触, ${F_{{\rm{p}} \text{-} {\rm{w}}}}$ 为颗粒与壁面产生的接触力, ${F_{{\rm{p}} \text{-} {\rm{f}}}}$ 为颗粒与流体产生的接触力, $I$ $\omega $ 分别为颗粒的转动惯量和角速度, $T$ 为颗粒在质心处受到的合外力矩.

1.3 CFD-DEM双向耦合

CFD与DEM双向耦合过程如下:在每个时间步长中,首先采用CFD求解计算域内的流体,直至迭代计算收敛,然后根据颗粒所在的流体单元内储存的流体的物理参数(如速度、压力),采用DEM计算出每个颗粒所受的力,并在相应时间步长内进行一次(或若干次)迭代计算. 颗粒离散元仿真计算结束后,再次进行流体CFD仿真计算. 将所有颗粒的动量汇作用在相应的流体网格单元上,用于计算流体与颗粒之间相互能量的转换,最终实现DEM-CFD双向耦合[14].

本文所采用的DEM-CFD耦合模块是基于CFD中的离散相模型(discrete phase model, DPM),将DEM计算完成的颗粒位置及物理参数(如速度、加速度)等信息赋予DPM模型中的每一个颗粒,即DEM中的颗粒与DPM模型中的颗粒建立一一对应的关系,颗粒受力的计算及动量汇的传递都在CFD中完成. 因此,DPM模型仅作为媒介的作用,能极大地提高仿真计算效率.

2 数值仿真模型与验证 2.1 几何模型及参数

建立水力旋流器流体流经的空间,即流体域,如图2所示. 图3为水力旋流器的结构示意图. 进料口为矩形. 水力旋流器的几何参数,如表1所示. 表中 $D$ 为圆柱直径, $a \text{、} b$ 为入口长、宽, ${d_{\rm{o}}}$ 为溢流管管径, ${d_{\rm{s}}}$ 为底流管管径, $H$ 为圆柱体高度, ${h_{\rm{o}}}$ 为溢流管插入深度, $\theta $ 为锥角.

图 2 水力旋流器的几何模型 Fig. 2 Geometry of hydrocyclone
图 3 水力旋流器结构示意图 Fig. 3 Structural schematic of hydrocyclone
表 1 水力旋流器几何参数 Table 1 Geometry parameters of hydrocyclone
2.2 初值与边界条件

水力旋流器入口水合物浆体中颗粒的体积分数为0.6%,其中水合物颗粒体积分数占总固体颗粒体积分数的75%,砂颗粒体积分数占总颗粒体积分数的25%. 基于DEM-CFD耦合仿真的方法对水力旋流器进行仿真分析,设置入口为速度边界,流体与颗粒的入口速度均为7 m/s;溢流口和底流口均为压力出口,出口为压力条件;壁面采用无滑移壁面边界;采用RSM湍流模型. 颗粒与颗粒、颗粒与壁面的碰撞采用Hertz-Mindlin接触模型. 将水合物颗粒和砂颗粒视为球形颗粒处理,颗粒直径为2 mm,液相海水的密度为1 025 kg/m3,黏度为0.001 72 Pa·s,颗粒的参数如表2所示. 表中 $\rho $ 为颗粒密度, $\gamma $ 为泊松比, ${f_{\rm{s}}}$ 为颗粒静摩擦系数, $f$ 为颗粒滚动摩擦系数, ${f_{\rm{r}}}$ 为材料恢复系数.

表 2 固相颗粒的材料物性参数 Table 2 Material properties of solid particles
2.3 网格独立性验证

水力旋流器计算域采用六面体结构网格,如图4所示. 为验证仿真模型的网格独立性[16],采用上述相同的物理参数和边界条件对4种不同网格数量的模型进行数值模拟,得到计算结果如表3所示. 表中, ${Q_{{\rm{\beta NGH}}}}$ 为溢流口水合物质量流量, ${Q_{{\rm{\gamma Sand}}}}$ 为底流口砂质量流量,pmax为最大压力, ${v_{\max }}$ 为最大速度. 从表3中可以看出,网格数量对仿真结果有一定的影响,溢流口水合物质量流量和底流口砂质量流量随网格数量增加到一定程度后计算结果趋向稳定. 网格数量越多仿真计算所需要的时间越长,因此综合考虑仿真时间和仿真结果,选用网格数量为194 819个的模型可满足网格独立性.

图 4 水力旋流器计算网格 Fig. 4 Grid of hydrocyclone current simulation
表 3 不同网格数量的模型计算结果比较 Table 3 Comparison of calculation results of models with different grid numbers
2.4 CFD-DEM耦合可信性

为验证耦合仿真的可信性,采用文献[4]中的水力旋流器的实验数据与DEM-CFD耦合仿真计算结果进行对比,如图5所示为不同截面径向位置R的切向速度 ${v_{\rm{t}}}$ 和轴向速度 ${v_{\rm{a}}}$ . 图中为距顶部为60 mm和120 mm截面上速度比较图. 图中切向速度仿真数据与实验数据较接近,仅在壁面处和中心区域有些不同,这是因为实验测量误差等原因造成的. 图中轴向速度仿真数据和实验数据吻合较高,表明所采用的耦合计算方法是可信的.

图 5 不同截面的速度分布仿真结果与实验结果比较 Fig. 5 Comparison of simulation results and experimental results of velocity distribution at different sections
3 分离效率与颗粒运动讨论 3.1 分离效率

分离效率是衡量水力旋流器分离过程进行完善程度的一个技术指标,取值范围在0~1.0,计算出的分离效率值越大,表明其分离效率越高,分离更加彻底. 本文利用汉考克效率计算式[15]来对水力旋流器的分离效率进行计算.

分离效率表达式为

$\eta = \frac{{(\alpha - \gamma )(\beta - \alpha )}}{{\alpha (\beta - \gamma )(1 - \alpha )}},$ (9)
$\alpha = \frac{{{Q_{{\rm{\alpha NGH}}}}/{\rho _{{\rm{NGH}}}}}}{{{Q_{{\rm{\alpha NGH}}}}/{\rho _{{\rm{NGH}}}} + {Q_{{\rm{\alpha sand}}}}/{\rho _{{\rm{sand}}}}}},$
$\beta = \frac{{{Q_{{\rm{\beta NGH}}}}/{\rho _{{\rm{NGH}}}}}}{{{Q_{{\rm{\beta NGH}}}}/{\rho _{{\rm{NGH}}}} + {Q_{{\rm{\beta sand}}}}/{\rho _{{\rm{sand}}}}}},$
$\gamma = \frac{{{Q_{{\rm{\gamma NGH}}}}/{\rho _{{\rm{NGH}}}}}}{{{Q_{{\rm{\gamma NGH}}}}/{\rho _{{\rm{NGH}}}} + {Q_{{\rm{\gamma sand}}}}/{\rho _{{\rm{sand}}}}}},$

式中: $\eta $ 为汉考克综合分离效率, $\alpha $ 为入口中水合物的体积分数, $\beta $ 为溢流口水合物的体积分数, $\gamma $ 为底流口水合物的体积分数, ${Q_{{\rm{\alpha NGH}}}}$ ${Q_{{\rm{\alpha sand}}}}$ 分别表示入口中水合物与砂的质量流量, ${Q_{{\rm{\beta NGH}}}}$ ${Q_{{\rm{\beta sand}}}}$ 分别表示溢流口中水合物与砂的质量流量, ${Q_{{\rm{\gamma NGH}}}}$ ${Q_{{\rm{\gamma sand}}}}$ 分别表示底流口中水合物与砂的质量流量.

在仿真结果中,提取溢流口和底流口处水合物质量流量分别为0.018 92 kg/s2、0.000 84 kg/s2,砂的质量流量分别为0 kg/s2和0.019 66 kg/s2. 由分离公式(式(9))计算得到分离效率为95.8%. 改变入口速度,分析入口速度对水合物浆体分离效率的影响,不同速度溢流口和底流口水合物和砂质量流量如表4所示. 表中 ${Q_{{\rm{\beta NGH}}}}$ 为溢流口水合物质量流量, ${Q_{{\rm{\gamma NGH}}}}$ 为底流口水合物质量流量. 由分离公式计算得到不同入口速度 $v_{\rm{inlet}}$ 的分离效率 $\eta $ ,如图6所示.

表 4 不同速度下溢流口和底流口水合物和砂的质量流量 Table 4 Mass flow of NGH and sediment of overflow and underflow under different velocities

图6中可以看出,随着入口速度的增加,水力旋流器分离效率呈现先增加,然后略微有所下降的规律. 入口速度越低,水力旋流器内部的切向速度场的均值越低. 由于水合物颗粒密度小于液体密度,在离心浮力和离心惯性力的共同作用下,水合物颗粒通常会沿径向向水力旋流器的中心轴线方向运动. 入口速度增加,水合物颗粒的离心浮力与离心惯性力的差值增加,使得水合物颗粒沿径向向中心轴线运动的速度也会增加,最终到达内旋流场的可能性也相应增加,因此水合物颗粒分离的效率也会增高. 如图6所示的AB段. 当入口速度增加到某一特定值时,水力旋流器的分离效率增加达到峰值. 当进一步提高入口速度时,由于水力旋流场(包括内旋流场和外旋流场)的湍流强度过大[17-19],使得水合物颗粒径向运动的干扰增加,同时水合物颗粒的分散程度也会增加,反而会使水合物颗粒分离效率有所下降[20-21],如图6所示的BC段.

图 6 不同入口速度水力旋流器分离效率曲线 Fig. 6 Curve of hydrocyclone separation efficiency under different velocity

图7为不同入口速度下水合物浆体处于稳定流动状态时,水合物浆体中水合物颗粒群砂颗粒群分布图. 图中黑色颗粒表示水合物颗粒,灰色颗粒表示砂颗粒. 水合物颗粒大多数分布在水力旋流器中心轴线附近,处于内旋流,有向上运动从溢流管排出的趋势. 沙颗粒大多数分布在水力旋流器的内壁面,处于外旋流,有向下运动从底流口排出的趋势. 为了定量地描述在水力旋流器内部的水合物颗粒和沙颗粒的分布状态,从坐标原点位置开始沿z轴负方向每隔30 mm取截面,统计各截面的颗粒数,如图8所示. 图中NNGH为水合物颗粒数,Nsand为砂颗粒数. 当入口速度较低只有4 m/s时,由于水合物颗粒的离心浮力与离心惯性力的差值较小,其径向速度也较小,以致于有些水合物颗粒径向运移还没有进入内旋流,就被外旋流场中的液体带入了底流口,导致水合物颗粒沿轴线方向分布比较分散,因此入口速度较低时分离效率较低,如图8(a)中方块实线所示. 如图7(b)所示,随着入口速度增加到7 m/s时,可以看出水合物颗粒在水力旋流器中部的中心区域聚集的越来越多,并且大部分已经进入内旋流,如图8(a)圆形虚线的轴向距离为200~350 mm水合物颗粒数,因此水力旋流器的水合物颗粒分离效率越来越高. 当入口速度达到8 m/s时,分离效率达到峰值. 继续增加入口速度时,分离效率反而会有所降低. 可以看出,水力旋流器内砂颗粒分布情况基本一致,由于砂的离心惯性力远远大于离心浮力,使得砂颗粒会沿壁面螺旋向下运动,最终从底流口排出,如图8(b)所示.

图 7 不同入口速度水合物浆体的颗粒分布 Fig. 7 Flow distribution of NGH slurry under different velocities
图 8 不同速度水力旋流器内颗粒分布曲线图 Fig. 8 Curve of particle distribution in hydrocyclone under different velocities
3.2 水合物颗粒运动分析

在水力旋流器入口截面上,水合物颗粒所处的位置区域不同,在水力旋流器内部流场中的运动轨迹不同,直接影响着水合物颗粒的分离效率. 为此,本文以较佳入口速度的分离效率工况为例,讨论入口速度为8 m/s的水合物颗粒在水力旋流器内部流场中的运动轨迹.

将入口截面分为5个典型位置区域,如图2所示,第I区域为入口上侧,第Ⅱ区域为入口外侧,第Ⅲ区域为入口下侧,第Ⅳ区域为入口内侧,第Ⅴ区域为入口中心区域. 为了考虑水合物颗粒运动的随机性影响,在耦合仿真结果中,从每个区域里分别提取5个水合物颗粒的运动轨迹数据并用不同曲线绘制三维空间轨迹线图,如图9所示. 从图9中可以下列基本规律.

图 9 入口截面不同区域进入的水合物颗粒运动轨迹 Fig. 9 Trajectories of NGH particles from different areas of inlet cross-section

1)如图9(a)、(d)所示,入口上侧I区和入口内侧Ⅳ区的水合物颗粒运动轨迹的总体规律类似. 从这2个区域进入的水合物颗粒,由于水合物颗粒的离心浮力大于离心惯性力,使得水合物颗粒会沿径向向中心运动. 位于入口上侧I区的水合物颗粒进入旋流场后,由于靠近水力旋流器的顶面,受到流体向下驱动的力很小,水合物颗粒在靠近顶面附近,以盘状螺旋轨迹形式向溢流管的外壁靠拢,再沿外壁向下运动,形成了短路流,进入溢流管排出. 位于入口内侧Ⅳ区的水合物颗粒进入旋流场后,由于靠近中心轴线距离最近,最容易在离心浮力作用下向溢流管外壁靠拢,形成短路流. 如果这2个区域的水合物颗粒浓度越高,水力旋流器的分离效率也会越高.

2)如图9(c)、(e)所示,入口下侧的Ⅲ区和入口中心的Ⅴ区的水合物颗粒运动轨迹的总体规律也类似,属于水合物颗粒的正常分离轨迹. 从这2个区域进入的水合物颗粒,由于进入了水力旋流器的外旋流场,在外旋流场的带动下,水合物颗粒是以空间螺旋轨迹向下运动. 水合物颗粒同时在向水力旋流器的中心运动,当其突破轴向零速面,到达水力旋流器的内旋流场后,水合物颗粒开始边螺旋向上运动,最终进入溢流管排出.

3)如图9(b)所示,入口外侧的Ⅱ区的5个水合物颗粒,其中有4个属于正常分离轨迹,只有一个颗粒被从底流口排出. 整体5个水合物颗粒在外旋流带动下,一直以空间螺旋方式向下运动,有4个水合物颗粒进入内旋流,开始向上部溢流管运动,最终被分离. 只有一个颗粒被外旋流带入了底流口. 由此说明从入口外侧进入的水合物颗粒分离难度最大,高于所有其他入口区域.

3.3 入口外侧水合物颗粒力学行为分析

入口外侧Ⅱ区的水合物颗粒有4个属于典型水合物颗粒分离轨迹,只有一个颗粒从底流口排出. 为了研究其运动变化规律和力学行为,选择从溢流口排出的水合物颗粒和从底流口排出的水合物颗粒作为代表. 根据计算结果,绘制其位置L、速度v和受力F随时间t变化曲线图,如图10~12所示.

图 10 入口截面外侧进入典型水合物颗粒位置变化 Fig. 10 Variations of typical NGH particle position from outside areas of inlet cross-section
图 11 入口截面外侧进入的典型水合物颗粒速度变化 Fig. 11 Variations of typical NGH particle velocity from outside areas of inlet cross-section
图 12 入口截面外侧进入典型水合物颗粒受力变化 Fig. 12 Variations of typical NGH particle force from outside areas of inlet cross-section

图10(a)(b)所示,左边为水合物颗粒位置曲线,右边为不同时刻水合物颗位置示意图,即颗粒运动轨迹. 在水力旋流器柱锥交界z=–180 mm处,图中z方向位置变化曲线上点A的时间分别为t=3.072 s、2.911 s,此时,对应径向位置变化曲线于点B. 如图10(a)(b)两图方块实线所示,点B处半径值均为73.8 mm,而且点B之前的阶段,颗粒的径向位置基本保持恒定,说明水合物颗粒都在靠近壁面的边界层中运动. 由于边界层附近液体的流速较低,当水合物颗粒因惯性力进入流体边界层后,大部分颗粒就会被束缚在边界层中随流体慢速流动. 从总体情况来看:从入口外侧Ⅱ区进入的水合物颗粒到达柱锥交界处以前,都是贴着水力旋流器内壁螺旋向下运动. 在运动过程中,水合物颗粒与壁面之间存在着随机碰撞力,如图12所示. 水合物颗粒的总体速度值也都在逐渐减小,如图11中黑色实线所示,说明水合物颗粒的动能在减小. 经过水力旋流器柱锥交界处时,由于受到几何形状突变的影响,使得入口外侧进入的五个水合物颗粒中的4个颗粒发生了突变,跳出了水力旋流器内壁的边界层流. 在离心浮力的作用下,水合物颗粒快速沿径向向中心运动. 直到突破轴向零速面(如图10(a)中的点C图11(a)中点C1),水合物颗粒进入内旋流场,然后以螺旋向上的运动轨迹,向溢流口方向运动. 在经过溢流管入口处D时,由于流体的局部湍流和旋涡作用,水合物颗粒的受力发生剧烈波动(如图12(a)中的点D2附近),速度的各个分量也发生剧烈波动(如图11(a)中的点C1附近),但是总体平均速度值呈现增加,如图11(a)所示.

入口外侧进入的5个水合物颗粒中只有一个颗粒受柱锥交界面突变的影响不大,仍然紧贴着水力旋流器锥面段的内壁,继续螺旋向下运动,直到图10(b)中的点C,进入底流排出流阶段. 在水合物颗粒垂直向上浮力和垂直向下的砂水裹夹力的双重作用,如图12(b)中的方块实线的z方向受力曲线,水合物颗粒z方向速度有微小的正负交替,如图11(b)中的三角形实线的z向速度曲线,最终结果是水合物颗粒z方向位置在微小波动中,向底流口方向移动并被排出.

4 结 论

(1)相比CFD分析方法,CFD-DEM耦合多相流分析方法的计算结果更接近于实验结果,可以获取多相流中的固相颗粒的空间轨迹、速度和受力的全部数据.

(2)用汉考克效率公式作为评价标准,根据入口、溢流口和底流口的水合物和砂的质量流量得到入口速度对水力旋流器分离效率的影响. 入口速度为8 m/s时水合物的分离效率较好,原因是入口速度越高,水合物颗粒离心浮力越大,分离效率越高;但是过高的入口速度会导致湍流过于激烈,反而会降低水合物颗粒的分离效率. 因此,必须选择合适的入口速度,才能达到最佳分离效果.

(3)当水合物颗粒从入口上侧和入口内侧进入时,水合物颗粒容易形成短路流而直接从溢流管排出. 这2个区域水合物颗粒浓度越高,分离效率也会越高. 当从入口下侧和入口中心进入时,以空间螺旋轨迹线向下运动. 同时沿径向向水力旋流器中心运动. 到达水力旋流器的内旋流场后,水合物颗粒则改变方向,螺旋向上运动,最终进入溢流管被排出. 这一类轨迹属于典型的水合物颗粒分离运动轨迹. 当水合物颗粒从入口外侧进入时,其中大多数水合物颗粒运动属于典型分离运动. 但是有少部分水合物颗粒,紧贴水力旋流器的内壁沿边界层或在边界层附近,螺旋向下运动,直至从底流口排出.

参考文献
[1]
李慧. 固-液分离用水力旋流器的三维数值模拟研究[D]. 郑州: 郑州大学, 2006.
LI Hui. The 3D numerical simulation study on solid-liquid hydrocyclone [D]. Zhengzhou: Zhengzhou University, 2006. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y947345
[2]
DOBY M J, NOWAKOWSKI A F, YIU I, et al. Understanding air core in hydrocyclone by studying pressure distribution as a function of viscosity[J]. International Journal of Mineral Processing, 2008, 86(1–4): 18-25.
[3]
王剑刚.三维旋转湍流场激光测速研究[D]. 上海: 华东理工大学, 2016: 37–39.
WANG Jian-gang. Three-dimensional measurement of turbulent hydrocyclone flow field with optical methods [D]. Shanghai: East China University of Science and Technology, 2016: 37–39. http://cdmd.cnki.com.cn/Article/CDMD-10251-1017035693.htm
[4]
HSIEH K T. Phenomenological model of the hydrocyclone [D]. Salt Lake City: The University of Utah, 1988: 228–234. http://www.researchgate.net/publication/222265370_Phenomenological_model_of_the_hydrocyclone_Model_development_and_verification_for_single-phase_flow
[5]
BHASKAR K U, MURTHY Y R, RAJU M R, et al. CFD simulation and experimental validation studies on hydrocyclone[J]. Minerals Engineering, 2007, 20(1): 60-71. DOI:10.1016/j.mineng.2006.04.012
[6]
王志斌. 水力旋流器分离过程非线性随机特性研究[D]. 四川: 四川大学, 2006: 68–71.
WANG Zhi-bin. Research on nonlinear stochastic characteristics of separation processes in hydrocyclones [D]. Sichuan: Sichuan University, 2006: 68–71. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y1002212
[7]
徐继润, 罗茜, 邓常烈. 水力旋流器的径向速度[J]. 有色金属(选矿部分), 1985(5): 12-17.
XU Ji-run, LUO Qian, DENG Chang-lie. Radial velocity of hydrocyclone[J]. Nonferrous Metals (Mineral Processing Section), 1985(5): 12-17.
[8]
MOUSAVIAN S M, NAJAFI A F. Numerical simulations of gas-liquid-solid flows in a hydrocyclone separator[J]. Archive of Applied Mechanics, 2009, 79(5): 395-409. DOI:10.1007/s00419-008-0237-2
[9]
QIU L C, WU C Y. A hybrid DEM/CFD approach for solid-liquid flows[J]. Journal of Hydrodynamics, 2014, 26(1): 19-25. DOI:10.1016/S1001-6058(14)60003-2
[10]
CHU K, CHEN J, YU A. Applicability of a coarse-grained CFD-DEM model on dense medium cyclone[J]. Minerals Engineering, 2016, 90: 43-54. DOI:10.1016/j.mineng.2016.01.020
[11]
CHU K W, WANG B, XU D L, et al. CFD–DEM simulation of the gas-solid flow in a cyclone separator[J]. Chemical Engineering Science, 2011, 66(5): 834-847. DOI:10.1016/j.ces.2010.11.026
[12]
VERSTEEG H K, MALALASEKERA W. An introduction to computational fluid dynamics: the finite volume method [M]. New York: Wiley, 1995: 9–16. http://www.mendeley.com/catalog/introduction-computational-fluid-dynamics-3/
[13]
CUNDALL P A, HART R D. Numerical modeling of discontinue, engineering computations, 1992, (9)2: 101–113. In: Comprehensive Rock Engineering, 1993, 2: 231–243. https://www.researchgate.net/publication/242055475_Numerical_modeling_of_discontinua
[14]
胡国明. 颗粒系统的离散元素法分析仿真[M]. 武汉: 武汉理工大学出版社, 2010: 8–11. http://www.bookask.com/book/1108320.html
[15]
庞学诗, 水力旋流器理论与应用[M]. 长沙: 中南大学出版社, 2005: 52–57 http://www.bookask.com/book/1305910.html
[16]
冯静安, 唐小琦, 王卫兵, 等. 基于网格无关性与时间独立性的数值模拟可靠性的验证方法[J]. 新疆: 石河子大学学报, 2017, 35(1): 52-56.
FENG Jing-an, TANG Xiao-qi, WANG Wei-bing, et al. Reliability verification method of numerical simulation based on grid independence and time independence[J]. Xinjiang: Journal of Shihezi University, 2017, 35(1): 52-56.
[17]
黄波, 陈晶晶. 重介质旋流器颗粒运动特性分析及数值模拟[J]. 煤炭工程, 2015, 47(5): 115-117.
HUANG Bo, CHEN Jing-jing. Characteristic analysis and numerical simulation of particle motion in dense medium cyclone[J]. Coal Engineering, 2015, 47(5): 115-117.
[18]
喻黎明, 邹小艳, 谭弘, 等. 基于CFD-DEM耦合的水力旋流器水沙运动三维数值模拟[J]. 农业机械学报, 2016, 47(1): 126-132.
YU Li-ming, ZOU Xiao-yan, TAN Hong, et al. 3D numerical simulation of water and sediment flow in hydrocyclone based on coupled CFD-DEM[J]. Transactions of The Chinese Society of Agricultural Machinery, 2016, 47(1): 126-132.
[19]
胡坤. ANSYS CFD疑难问题实例详解 [M]. 北京: 人民邮电出版社, 2017. http://www.bookask.com/book/1662082.html
[20]
黄思. 流体机械的数值仿真研究及应用[M]. 广州: 华南理工大学出版社, 2015. http://www.bookask.com/book/2134570.html
[21]
陈云富. 旋转流场中颗粒运动与分布的特性研究[D]. 江苏: 东南大学, 2016.
CHEN Yun-fu. Investigation on the particle motion and distribution characteristics in rotational flow field [D]. Jiangsu: Southeast University, 2016. http://cdmd.cnki.com.cn/Article/CDMD-10286-1016328323.htm