基于cfd的空气柱固液分离水力旋流器内清水流场三维数值模拟研究.docxVIP

基于cfd的空气柱固液分离水力旋流器内清水流场三维数值模拟研究 0 水力旋流器在固-液分离应用中的应用 水流补偿器是利用各相密度差制备的液体非均匀相混合物的分离装置。尽管早在1891年水力旋流器已经问世,然而其有关的研究与应用所涉及的大都是分级方面,应用于固-液分离方面的则比较少。与重力分离设备相比,水力旋流器结构紧凑、体积小、质量轻、易于设计与安装、维修费用低、易于调节与控制。鉴于这些优点,水力旋流器在固-液分离领域的应用正越来越受到关注。 严格地讲,由于水力旋流器内进行的是固-液分离现象,水相和颗粒相是相互影响的。然而,由于水力旋流器内部存在着内旋流、外旋流、循环流和短路流以及空气柱,流动情况十分复杂。因此,从研究方法的角度,通常先研究清水流场,在此基础之上再进一步研究固-液分离现象。鉴于旋流器内流场的复杂性,对水力旋流器流场进行精确的解析求解十分困难。而单纯通过实验的方法,周期长,费用高,因此本文作为研究固-液分离现象的基础,采用数值模拟技术,利用CFD商业软件包FLUENT软件,对固-液分离用水力旋流器内的清水流场进行了数值模拟研究。 国内的相关文献中对水力旋流器内清水流场的模拟大都忽略了空气柱的影响,或是利用“刚盖法”,在建模的时候直接将旋流器中心部分挖去一圆柱形来取代空气柱的存在,然而这对于实际流场的深入研究是不准确的。本文首次使用“mixture”模型对空气柱的最终形成状态进行了模拟,其效果良好。 1 结构参数的确定 根据生产能力Q=60m3/h,对应的名义直径D选择了250mm。应用里特马关系式进行验算,即: 式(1)中,Q—生产能力(进料量,L/S),Q为60m3/h,即16.67 L/S;△P—压强降,取△P=105N·m-2;ρi—悬浮液的密度,由于针对的是泥沙和水的分离,所以,当泥沙密度为3000kg/m3,泥沙体积浓度最高为10%的情况下,ρi=0.91000+0.13000=1200 kg/m3;K—系数,常取K=1.3。将各值代入式(1)中,计算得:D=246mm,取标准直径D=250mm。 本文的模拟结构参数选择如图一所示,其中主要几何参数如下: 为更符合实际工况,底流口延长一平行直管段,长度为50 mm。旋流器锥角α=15°。 2 数学模型 2.1 水力旋流器的结构简化模型 由于水力旋流器内的流动十分复杂,为便于对其内部流动细节进行细致研究,在建立清水流场数学模型之前,对实际使用的水力旋流器作了一些简化,即: (1)不考虑由于流体与壁面的摩擦而产生的热效应; (3)假设流场恒温,不考虑能量传递。 2.2 不可压缩流体的动力学方程 实际工程应用的水力旋流器一般工作在湍流区,因此描述其行为的模拟必然会涉及到湍流方程。湍流的基本方程组包括:连续性方程、Navier-Stokes方程、雷诺应力方程、湍动能方程和湍动能散率方程。对于不可压缩流体: 连续性方程 动量方程 式(3)中:为未知雷诺应力(i,j=1,2,3)。 由于方程式(3)中出现了新的脉动速度的二阶关联项,即雷诺应力项,致使它与连续性方程(2)不能封闭,应选择适当的湍流模型来对项进行模拟,以使方程组封闭。 2.3 可降压器模型的求解 工程中常用的湍流模型有标准k-ε模型、RNGk-ε、雷诺应力模型(RSM)。前两者都没有突破涡粘性假设下的各向同性的框架,而旋流器内是基于各向异性的强三维旋流,所以本文采用了雷诺应力模型(RSM)。 在该模型中,确定雷诺应力各分量的输运方程如下: Suser是用户自定义的源项 经过模型化以后,湍动能方程和湍动能耗散率方程就变成如下形式: 式(6)中:Cμ=0.09,σk=0.82,σε=1.0,Cε1=1.44 Cε2=1.92; Sk,Sε是用户自定义的源项。 由此便构成了雷诺应力模型(RSM)的基本封闭方程组。 2.4 适当降低网格数量 水力旋流器内的实际流动是在不规则的区域中进行的,要实现任意区域内的数值模拟,较精确地拟合边界形状和边界条件,则要使用贴体坐标(body-fitted coordinates)法。在划分网格时,要求网格能够适应计算区域中物理量的变化情形,即在物理量变化比较剧烈的地方适当地加密网格数量以保证计算的精度,在物理量变化平缓的地方适当地减少网格数量以减少计算量。考虑边界条件的离散化,希望网格线与物理区域的边界线正交,以利于准确地计算边界上的物理量,同时应尽量使得网格分布与计算域的几何形状保持一致。 在靠近外壁面和中心内筒的区域,网格取得较密,远离壁面的区域网格取得较为稀疏。计算网格的总数目由沿高度方向上的网格密度控制。同时为更好地捕捉实际流动中的边界特征,在横截面上划分了边界层网格。最终,兼顾考虑到计算精度和计算机运算速度,将整个三维模型划分了160622个网格单元,如图二、三、四