您好,欢迎访问三七文档
当前位置:首页 > 临时分类 > 旋流器内流特性分析的思路和方法
旋转条件下旋流器内流特性分析的思路和方法对于静态旋流分离器而言,混合液由切向入口进入旋流器内部,形成两个相同方向的涡流区,及内部在近轴区的强制涡和外部近壁区的自由涡。外部涡流区向下流动,从底流出口流出,内部涡流区从反流点开始向上流动,从溢流出口流出。同时,在入口强旋流区出现循环流动,流体在这个区域停留时间较长,轴向速度很小,当流入小锥部分时,循环流基本消失。而在旋转条件下,旋流器内流体做复杂的旋转湍流流动,管内流体微团在一个有限的区域内发生强烈的旋转、剪切与变形,加之壁面旋转的复合效应,使得旋流器内部流场变得十分复杂。1基本方程为合理的简化数值计算模型,对旋转条件下旋流器内流体做如下假设:(1)流体为不可压缩牛顿流体;(2)管内流体与外界无热量交换,且其温度不变;(3)管内流体的加速度远远大于重力加速度,故不考虑流体的重力。流体流动的质量守恒方程即连续性方程为:div0tu(1-1)式中u是速度矢量,div表示散度。动量守恒方程表述为:divdivgradiipStxuuuu(1-2)式中,ix表示x,y,z三维坐标,iS表示uS,vS,wS是动量方程在上述三个方向的广义源项。依照上述假设条件,系统不存在热交换,无需满足能量守恒定律。同时,在仅考虑不可压流动状态下,忽略重力作用的影响,湍流瞬时控制方程可以简化为如下形式:div0u=(1-3)1divdivgradiptxuuuu(1-4)2输运方程目前常用的湍流模型包括涡粘模型与雷诺应力模型两大类,由于两方程的kε模型及其改进形式并不适用于强旋流动的数值模拟,故选取了雷诺应力方程模型(RSM)。在均匀不可压湍流中,无量纲雷诺应力的输运方程表示为:1Re2ReijijijijijijjiijkjikikjikkjkkkkPDjjiijikkDuuuuuuuuupupuuuuuDtxxxxuuuupxxxx(2-1)其中瞬态项与对流项表示成为物质导数的形式,ijD、ijP、ij与ij分别为扩散项,剪应力产生项,压力应变关联项与耗散项。上式中,ijP只包含二阶关联项,无需进行处理;而ijD、ij与ij均包含未知的关联项,必须要建立合理的模型,才能使整个计算模型封闭。扩散项ijD可以分解为湍动扩散项,TijijkjikikjkDuuupupux与粘性扩散项,1ReijLijkuuDx,湍动扩散项一般采用下式模化,,ijtTijkkkuuDxx(2-2)耗散项表示分子粘性对雷诺应力产生的耗散。建立耗散项的计算公式时,认为耗散过程主要发生在小尺度涡区,小尺度涡团结构接近各向同性,可以忽略各向异性耗散,即认为是局部各向同性,根据上述假设,建立耗散项计算公式:23ijij(2-3)压力应变关联项中,二阶应力(SSG)模型被认为是更适合预测高应变率及流线弯曲程度较大的流动,其表达式如下:()()*112*33451()()32()()()3sijrijijijikjkklklijijijijikjkjkikklklijikjkjkikCCPbCbbbbCCbbkSCkbSbSbSCkbb(2-4)式中:123ijijijuubk为雷诺应力各向异性张量;13.4C,*11.80C,24.2C,30.8C,*31.30C,41.25C,50.40C为模型系数;2kkPP;12iiuku为湍动能;12jiijijuuSxx为应变率张量;jiijijxuxu21为旋转率张量。3离散方程由应变量在节点之间的分布假设及推导离散方程的方法不同,离散化方法一般可分为有限差分法、有限元法和有限体积法等。其中,有限差分方法不适合于求解旋转变径圆管这类边界条件复杂的流动,而有限元法因求解速度远低于有限差分法与有限体积法,一般不用于求解本文这一类网格数量较大且流动规律复杂的湍流问题。因此,选择有限体积法作为旋转条件下旋流器数值计算的离散化方法。图3-1给出了直角坐标系下旋流器计算空间内某一节点P的控制体积,E,W,N,S,F,B分别代表了P点在三个坐标轴方向的邻点,e,w,n,s,f,b代表了控制体积P的六个界面。图3-1P点控制体积根据Guass离散定律有:divVAdVdAn(3-1)yzxPWENSBFebfswn其中V是控制体体积,A是控制体积的表面积,为向量形式的物理量,n为表面外法线方向的单位矢量。展开后divewnsfbewnsfbVAAAAAAdVdAdAdAdAdAdA(3-2)对上述瞬态湍流通用方程在控制体积P内积分,其结果的通用形式如下:PewnsfbVIIIIIISdV(3-3)其中,iI的每一项都由两部分组成:一个是对流项,表示为cI;另一个是扩散项,表示为dI;方程右端表示源项,通常采用简化方法处理。一、对流项的离散输运方程中对流项的离散是要将对流项转化成为控制体积界面上的值,根据湍流的控制方程。以界面e为例,对流项的离散形式为:eceeAIudAuA(3-4)其余界面上的对流项也可以用同样的方法计算得到。此处定义物理量F,F表示通过界面上单位面积的对流质量通量(convectivemassflux),简称对流质量通量,其表示式为:Fu(3-5)F在界面e上的值为:eeFu(3-6)这样,六个界面的对流项可以表示为:ceeeeIFA(3-7)c(3-8)cssssIFA(3-9)cnnnnIFA(3-10)cffffIFA(3-11)cbbbbIFA(3-12)二、扩散项的离散采用与对流项相同的方法对扩散项进行离散,同样由控制方程出发,界面e上的扩散项离散形式为:edeAeIdAAxx(3-13)采用同样的方法求出其他界面上的扩散项:dwwIAx(3-14)dssIAx(3-15)dnnIAx(3-16)dfbIAx(3-17)dbbIAx(3-18)三、源项的离散把稳态输运方程中所有未列入扩散项和对流项的其余项都归纳到源项中。一般情况下,源项不为常数,而是所求未知量的函数。因此,源项的处理十分重要,甚至是数值求解成败的关键所在。在有限体积法离散化过程中,应用较广泛的一种处理方法是将源项局部线性化,假定在未知量微小的变动范围内,源项可以表示成该未知量的线性函数。故在控制体积P内,源项可以表示为:CPPSSS(3-19)其中,CS为常数部分,PS是S随变量变化的曲线在P处的斜率。为了保证代数方程迭代求解过程收敛,PS应为负值。4计算方法SIMPLE(Semi-ImplicitMethodforPressure-LinkedEquations)算法是求解压力耦合方程的半隐式方法,多年来一直被广泛使用。旋转条件下旋流器为不可压缩稳态流动问题,计算网格以结构网格为主。这种情况下,SIMPLE算法与PISO及SIMPLEC等算法相比,有着计算稳定性强、收敛性好等优点。因此,选择SIMPLE算法对离散后的控制方程组进行求解计算。为解决在普通网格上离散控制方程时给计算带来的严重问题,SIMPLE算法一般采用交错网格来实现。5边界条件一、入口边界条件已知管内流体流量,故选择速度入口边界条件,假设流体由变径圆管两侧的切向入口垂直流入管内,根据流量计算速度值为:2inQuA=(5-1)入口的湍流流动参数采用湍流强度与水力直径进行定义:8/1)(Re16.0HDI(5-2)ininHCAD4(5-3)以上各式中:u为垂直于入口边界的流体速度,m/s;Q为流量,m3/h;inA为单个入口截面积,mm2;HD为水力直径,mm;inC为湿周,mm。二、出口边界条件出口边界按流动充分发展,选取为自由出流出口,通过设定两个出口的分流比来分配流体流量,即:0z(Φ,,,,uvwkε)(5-4)三、壁面边界条件旋流器变径圆管是按一定转速绕自身轴线旋转的。因此,壁面按照无滑移旋转壁面条件处理,即壁面上流体的速度与壁面具有相同的速度;压力取为第二类边界条件,即0pn。6小结针对旋转条件下旋流器内的特点,建立了其内部流体的数值计算模型,主要内容包括:(1)通过合理的假设,在直角坐标系下建立了旋转条件下旋流器的基本控制方程和输运方程。(2)基于有限体积法,对所建立的湍流控制方程进行离散化。采用SIMPLE算法对离散后的湍流控制方程进行计算。(3)确定合适的边界条件,构成了定解问题。
本文标题:旋流器内流特性分析的思路和方法
链接地址:https://www.777doc.com/doc-2358273 .html