双曲面搅拌机流场的数值模拟研究

双曲面搅拌机流场的数值模拟研究彭珍珍1,　赵恒文1,　郭聪聪1,　汪文生2,　曾德全2

(1.河海大学水利水电工程学院,江苏南京210098;2.南京贝特环保通用设备制造

有限公司,江苏南京210098)

摘　要:　应用计算流体力学(CF D)软件、利用多重参考系法(MRF)对新型的污水搅拌设备———双曲面搅拌机的流场进行数值模拟,考察了搅拌槽内水体的流态特性,并提出了优化建议。结果表明,未加挡板时,切向剪切是搅拌槽内水体的主要流动形式;在搅拌槽周围加上挡板后,可使槽内的水流形态重新分布,使轴向流速显著增加,立体循环运动更明显,更符合污水处理的搅拌形。

关键词:　计算流体力学(CF D);　双曲面搅拌机;　数值模拟

中图分类号:X703 文献标识码:C 文章编号:1000-4602(2009)19-0091-04

Nu m er i ca l S im ul a ti on of Flow F i eld i n Hyperbolo i d Sti rrer PENG Zhen2zhen1,　ZHAO Heng2wen1,　G UO Cong2cong1,　WANG W en2sheng2,

　ZENG De2quan2

(1.College of W ater Conservancy and Hyd ropo w er Engineering,Hohai U niversity,N anjing210098, China;2.N anjing B eite Environm en tal P rotection GE M anufacture Co.L td.,N anjing210098,China) Abstract:　The computati onal fluid dyna m ics(CF D)s oft w are and multi2reference fra me(MRF) method were app lied t o nu merically si m ulate the fl ow field in a ne w waste water stirring equi pment—hy2 perbol oid stirrer.The fluid fl ow characteristic in the stirred tank was investigated,and the op ti m izati on suggesti on was put f or ward.The results show that the tangential shear is the main fl ow pattern of the waste water in the stirred tank without the baffle.After putting the baffle ar ound the stirred tank,the fl ow patterns in the stirred tank can be redistributed,the axial vel ocity can be increased significantly,and the vertical circulati on of the waste water can be more obvi ous,which is more suitable for the require ment of the stirring pattern and vel ocity in the waste water treat m ent.

Key words:　co mputati onal fluid dyna m ics(CF D);　hy perbol oid stirrer;　nu merical si m ulati on

搅拌混合设备作为一种在污水处理中常用的工艺设备,可通过创建水流起到强化搅拌混合作用,防止活性污泥沉淀,被广泛用于混合池、厌氧池、水解生化池等工艺中。但由于搅拌设备在工作时会产生复杂的流体运动,目前搅拌设备的设计和应用对实测和工程师经验的依赖性比较强,这样就存在设计周期长、耗费大、精确度低等问题。将计算流体力学(CF D)数值模拟技术应用到搅拌设备的设计分析中,通过数值模拟软件求解描述过程,可以实现过程设计、优化以及放大。目前国内对CF D用于搅拌过程的研究很多,但是对双曲面搅拌机这种新型的污水处理设备的数值模拟研究还未见报道。因此,笔者采用CF D方法对叶轮直径为300mm的双曲面搅拌机的搅拌流场进行数值模拟,并提出了优化意见,以提升其搅拌混合效果。

1　数值模拟

111　计算方程的选取

旋转的搅拌桨作用于流体产生了复杂的三维湍

第25卷　第19期2009年10月

中国给水排水

CH I N A WATER&WASTE WATER

Vol.25No.19

Oct.2009

流流动,因此选择合适的湍流模型成为模拟的关键。搅拌器三维流场数值模拟采用了三维雷诺平均N-S方程,并用标准k-ε湍流模型使方程组闭合。已有研究表明,标准k-ε模型稳定、简单,可在较大的工程范围内得到足够的精度,在定常条件下搅拌流场的不可压流动可用以下方程组描述:

①　连续方程

　9

9t+

9(ρu)

9x+

9(ρv)

9y+

9(ρw)

9z=0

(1)

式中　u、v、w———相对速度分量

②　动量方程

　9u i

9t+

9

9x i(u i u j)=F i-

1

ρ

9ρ

9x i+v

92u j

9x i9x j

(2)

③　湍动能方程

　9u i k

9x i-

9

9x i[(v+

v t

σ

k

)

9k

9x i]=P r-

ε(3)

式中　u

i

———沿i方向的速度分量

　f

i

———沿i方向的质量力

112　搅拌设备的CF D分析方法

从数值模拟的角度来看,模拟搅拌槽的一大难题是如何处理好运动的桨叶和静止的挡板及槽体之间的相互作用。为此,许多研究者提出了各自不同的解决方法,目前应用最广泛的是多重参考系法(MRF)和滑移网格法(SG)[1]。在同样的计算条件下,两种方法的计算结果基本一致,但MRF法的计算量要小一个数量级,而SG法对计算所需的硬件设备要求很高。因此笔者选用MRF法进行模拟。MRF法将计算域分成两个不重叠的区域,即搅拌桨区域和桨外区域,搅拌桨区域里的流体运动采用旋转坐标系,并且以搅拌桨的相同转速进行旋转,桨外区域的流体则采用静止坐标系,两个不同区域速度的匹配直接通过交界面上的插值转换来实现。

113　模型和网格划分

双曲面搅拌机的搅拌槽为圆柱形平底,直径为680mm,槽内液位高为

400mm,以水为工作介质。搅拌桨形式为双曲面叶轮(见图1),该曲面是曲线xy=b沿y轴旋转而构成的曲面体,在渐开双弧面上均匀布置了八条导流叶片,叶轮的最大直径为300mm。该形式的叶轮结构简单、便于维护保养,且其独特的设计能够有效地防止水中污染物的缠绕;另外,该叶轮的水力损失较小,搅拌所需要的功率小,达到了节能的目的[2]。

图1　双曲面搅拌机的叶轮

Fig.1　I m peller of hyperbol oid stirrer

利用Ga mbit软件对两个模型水域进行网格划分,采用非结构化网格,通过对不同数量网格的试算和比较,得出了网格无关解,最终确定了网格方案,此时网格单元的总数约为1185764个,并对叶轮表面进行了网格加密处理。

114　边界条件的设定

双曲面搅拌机的水流表面为自由面,CF D对自由面边界类型的处理一般有两种方法:刚体近似模拟和VOF模型模拟[3]。从理论上讲,VOF方法更接近于实际流动,但该方法所需时间较长,而刚体近似方法所需时间要少很多,且该方法也是目前应用数值模拟处理自由表面的主要措施。所以笔者对自由边界条件的处理采用刚体近似方法。

2　结果与分析

211　宏观流动特性

双曲面搅拌器内的流动区域可分为5部分:包含叶轮在内的叶轮区、叶轮下方的排出区、槽底与槽壁交接处的角区、槽中心轮毂下方的锥形体区域和占槽中大部分体积的主循环区域。未加挡板时,双曲面搅拌器槽内流体在叶轮的作用下形成径向射流,水流运动主要表现为沿着叶轮的径向和切向运动,形成整个水体的轴向循环运动。

212　时均速度分布规律

根据实际应用惯例,槽内液体流速≥0.3m/s 可作为搅拌混合效果的重要判定依据。以搅拌转速为230r/m in为例,在槽中间的最大纵截面上,分别在不同的流区内取4个等高线进行流场比较,这4个等高线的位置分别为:紧靠叶轮正下方Z=30 mm处、叶轮附近Z=75mm处、叶轮上方紧靠叶轮Z=140mm处和叶轮上方的主循环区中Z=300 mm处。通过这些位置的速度分布情况来初步判定其搅拌效果。搅拌槽中心轴向截面内的速度分布规律见图2(R为搅拌槽半径)。

由图2(a)可知,四个不同高度的绝对速度随半

第25卷　第19期 中国给水排水