泵站出水流道基本流态的分析

泵站出水流道根本流态的分析泵站出水流道根本流态的分析泵站出水流道根本流态的分析 泵站出水流道根本流态的分析摘要:为了改进泵站出水流道的水力设计方法,采用三维紊流数值模拟的方法,模拟了虹吸式、直管式和斜式三种型式出水流道内的流动形态;发现出水流道平面方向上的扩散情况较好,而立面方向在出口断面附近那么不同程度地存在着旋涡,该旋涡对流道出口断面的流速分布有很明显的影响;提出在流道设计时应最大限度地利用流道宽度方向的扩散,以免出口断面的有效面积过多地被旋涡挤占。关键词:泵站出水流道流态近十几年来,我国对水泵装置作了大量的研究工作,特别是对低扬程轴流泵水力模型和进水流道优化水力设计的研究已获得很多进展,有许多成果已经在泵站工程中得到成功的应用。由于种种原因,人们对进水流道内的流态比较注意、比较理解,而对出水流道内的流态那么缺乏较为深化的理解。在过去相当长的一段时期内,只做过一些关于出水流道水力损失方面的试验研究。出水流道是水泵装置的一个重要组成部分,对水泵装置的性能有非常明显的影响。出水流道的水力设计至今仍建立在传统的一维流动理论的根底上,这种理论与出水流道实际的三维流动情况出入很大。近些年来,人们对出水流道在水泵装置,尤其是在低扬程水泵装置中的作用,已经有了愈来愈清楚的认识,提出了重视研究出水流道水力设计理论和方法的要求。1997年9月1日公布施行的国家标准?泵站设计标准?所规定:“出水流道布置对泵站的装置效率影响很大,因此流道的型线变化应比较均匀〞[1]。流道的外特性是由其内特性决定的,对流道内特性的认识应是更为本质的认识。本文采用三维紊流数值模拟的方法,对虹吸式、直管式和斜式出水流道内的根本流态进展了初步的分析计算,力图提醒这三种形式出水流道内的三维流动形态,为认识和解决各类有关出水流道的水力学问题奠定必要的根底。1出水流道流动模拟的数学模型泵站出水流道三维流动模拟采用了雷诺平均N-S方程,并以标准κ-ε紊流模型使方程组闭合。选用这种模型的原因,是因为试验证明,标准κ-ε紊流模型对三维流动是非常适用的[2]。1.1控制方程在定常条件下,泵站出水流场的不可压流动可用以下一组方程描绘:1.连续方程(1)2.动量方程(2)3.紊动能方程(3)4.紊动能耗散率方程(4)(4a)式中,vt为涡粘性系数,可采用下式计算:vt=Cμ(κ2/ε)(4b)κ-ε模型中的有关常数为:这里,κ为vonKarman常数[3-5]。(1)~(4)式也可统一地用以下椭圆型守恒方程式表示:(5)式中,Φ表示具有守恒型的通用变量,ΓΦ为扩散系数,RΦ为源项。对应于Φ的特定含义,ΓΦ与RΦ相应地具有特定的形式。1.流厨口计算流场的进口设置在水泵后导叶出口断面,这里无疑是充分开展的流动。专门设计的试验说明,在设计条件下,导叶出口的环量很小[6],这里认为导叶出口的环量为零。因此,流厨口的边界条件仅提进口流速垂直于流厨口断面。2.流场出口计算流场的出口设置在距流道出口有一定间隔的出水池内,这里的边界条件近似按静水压力分布给出[7],即(6)3.固体壁面出水流道边壁、出水池底部等处均为固体壁面,其边界条件按固壁定律处理[8]。固壁边界条件的处理中对所有固体壁面的节点应用了无滑移条件,而对紧靠固体壁面节点的紊流特性,那么应用了所谓对数式固壁函数处理之。固壁处的摩擦速度可表为:u*=κuw/[ln(zw/dw)](7)式中,zw为靠近固壁的单元的中心至固壁的间隔,dw为固壁的绝对糙度,uw为该单元平行于固壁的速度分量。κ和ε可分别表为:(8)ε=u*3/κzw(9)4.自由外表出水池的外表为自由水面,假设忽略水面风引起的切应力及与大气层的热交换,那么自由面对速度和紊动能均可视为对称平面处理,而紊动能耗散率为(10)式中,zs为靠近自由外表单元的中心至自由外表的间隔,κs为该单元的紊动能量,CBE=0.07[9]。图1虹吸式出水流道的网格剖分图2直管式出水流道的网格剖分图315°斜式出水流道的网格剖分图4错列式网格示意图(X-Y平面)图5三维网格的方向通用微分方程(5)式的离散形式可表为:ΦP=(aEΦE+aWΦW+aNΦN+aSΦS+aHΦH+aLΦL+S)/(aE+aW+aN+aS+aH+aL+aP)(11)上式中的aE、aW、aN、aS、aH、aL分别为控制体各相邻网格点的系数,脚标中的大写字母E、W、N、S、H、L分别表示与控制体的东、西、北、南、上、下各个外表相邻的控制体的网格点(见图5),下同。这些系数的大小与网格的体积有关,它们表示控制体各相应交界面上对流和扩散作用产生的对网格点上有关因变量的影响。影响系数的计算分因变量为标量(如p、κ、ε等等)和因变量为速度矢量两种情况:1.因变量为标量时的系数计算(12a)(12b)(12c)aS=max(0,ds-12ms+max(0,-ms)(12d)(12e)(12f)在以上各式中,m表示通过控制体外表的对流强度,可为正,也可为负,决定于水流的方向(按迎风格式);d表示通过控制体外表的Φ的扩散率,恒为正;脚标中的小写字母e、w、n、s、h、l那么分别表示控制体的东、西、南、北、上、下等外表(见图4),下同。1)控制体外表对流强度的计算me=ρueAe(13a)mw=ρuwAw(13b)mn=ρunAn(13c)ms=ρusAs(13d)mh=ρuhAh(13e)ml=ρulAl(13f)2)控制体外表扩散率的计算de=ΓeAe/|PE|(14a)dw=ΓwAw/|WP|(14b)dn=ΓnAn/|PN|(14c)ds=ΓsAs/|SP|(14d)dh=ΓhAh/|PH|(14e)dl=ΓlAl/|LP|(14f)式中,A表示控制体外表的面积,PE那么表示网格点P至