用动网格模拟闸门开启过程非恒定水流特性
沙海飞1,2 周辉1,2吴时强1陈惠玲1
(1.南京水利科学研究院南京 210029;
2.水利部水科学与水工程重点实验室南京 210029)
【摘要】本文采用非结构动网格技术,采用网格变形与局部重构相结合的方法实现动网格,
网格的变形采用弹簧模型,对闸门开启过程非恒定水流特性进行了数值模拟,得到了不同开启
速度时的流量、压力分布以及不同流态的特征与发生、发展过程,结果令人满意。
【关键字】非结构动网格闸门非恒定水流数值模拟
0 前言
在水工建筑物中,闸门是用于控制水流最为常用的设施。在水坝、船闸等设计过程中,工程技术人员必须了解上下游水流的运动规律。由于闸门前后存在水位差,当水坝、船闸的闸门开启时,门后容易发生水流空化,时常造成空蚀破坏与闸门振动等危害[1]。文献[2]发现灌区水闸在闸门开启过程中,过闸水流为非恒定流,闸门开启的时段内闸下产生的冲坑深度,占一次放水过程产生冲坑深度的一半。可见在闸门开启过程非恒定水流是要引起足够重视的。而目前这方面的研究主要结合具体工程进行缩尺模型试验,闸门开启时的非恒定水流与复杂的漩涡运动增加了试验的难度,一般通过测量壁面压力来推断体型是否合理,对闸门后的内部流态了解很少。
随着计算机技术的飞速提高,计算流体力学得到了迅速发展,采用数值模拟闸门开启过程非恒定水流特性成为可能。文献[3]对假设闸门突然开启的情况进行了数值模拟,这和实际偏差较大。文献[1]对船闸弧门非恒定流态进行了二维计算,其计算网格是固定的,对弧门处网格进行了特殊处理,和实际也有一定的偏差。文献[4]用动网格技术对挑流冲刷过程进行了数值模拟,采用的是结构网格,用网格变形来实现动网格,这种方法很难用于闸门开启这种复杂的运动过程。随着非结构动网格技术的发展,以及在军事、仿生等领域的成功应用,使精确模拟闸门开启过程非恒定水流的特性成为可能。
本文采用非结构动网格技术,对平板闸门不同开启速度时非恒定水流进行了垂向二维数值模拟,探讨了闸门后流态发生、发展机理。
1 动网格技术
本文所处理的问题位移较大,所以采用网格变形与局部重构相结合的方法实现动网格。基本做法是,开设一个包围运动边界的窗口,窗口内的网格可以变形以适应边界的运动。以单元面积为判据,当窗口内即变形区出现严重扭曲的网格单元时,删除变形区内的网格单元,然后用阵面推进法重新生成该区网格。流动参数通过线性插值由旧网格映射到新网
中国水利学会第二首届青年科技论坛论文集
312
格。
控制网格变形的弹簧近似模型是将变形区网格看作一个弹簧网络系统,每一条边都认为是一根具有一定倔强系数的弹簧[5],网格点所受的合力为 ()1
i
N i ij j
i j f K x
x ==
-∑ (1)
其中ij K 为连结节点i 、j 的弹簧的倔强系数,i x 是节点i 的位置矢量,i N 是与节点i 相连的节点数。网格点的受力始终等于初始状态所受的合力,上式可表示节点i 的初始受力,则当运动边界上的点移动后,变形区内节点的新位置坐标可以通过迭代求解如下的线性系统得到。
111
222j j j
ik
i k j
sum sum ik
sumj i k
j K x f K a x f x f a K ≠≠??
-??????
??????-?
?????=?
?????????????????-??
?
?
∑∑
∑ (2) 0ik ik
k i k i a K ?=??若点不与点相连
若点与点相连,, (3)
该系统对角占优,Jacobi 迭代格式为: 111
i i
N N k k i
ij j i ij
j j x
K x f K
+==??=- ???
∑∑ (4)
当边界运动时,下个时刻的固定边界点和主动边界点的位置解已知,因此弹簧系统的边界
条件为Dirichlet 型,上式经过数次迭代就可达到满意的精度。
2 控制方程与数值求解方法
考虑是不可压缩水流流动问题,基本控制方程为Reynolds 方程,省略各变量的时间平均项标志“”,紊流模型采用ε-k 双方程模型,方程如下:
连续方程: 0u i x i
?=? (5)
动量方程: i
i j j u u u x Du p i f i Dt
x x i j γ???- ?
???
???=-+?? (6) k 方程: t k i i k Dk k G Dt x x ννεσ??
???? ? ???????
?
?=++-?? (7)
ε
方程: 212t k i i D C p C Dt x x k k εεενεεεενσ?????? ? ???????
??=++-?? (8) 式中,t 为时间;i u 和i x 分别为速度分量和坐标分量;ν为运动粘性系数;p 为修正压力;
i f 为质量力;2t u c k νε=为紊动粘性系数;方程中的经验常数0.09u c =, 1.0k σ=,
1.33εσ=,1 1.44C ε=,2 1.42C ε= k G 为平均速度梯度引起的紊动能产生项,
313 四、水力学 河流动力学
j i
i t k j j i u u u G x x x ν??????
? ???????
???+???=;对于雷诺应力项i j u u -的封闭问题,这里用k ε-双方程模型来处
理。
采用有限体积法对方程进行离散,对时间和空间均采用二阶精度格式,空间离散采用中心差分格式,对时间采用二阶精度的离散格式是很又必要的,因为研究的闸门开启过程水流随时间的演变是很重要的。压力校正采用SIMPLEC 算法,并采用显式校正步法对SIMPLEC 算法进行了显式校正。
上下游均是压力边界条件,给定上下游的总压,对于很长管段中闸门区域的局部模拟,可以采用一维计算为二维计算提供边界条件;壁面采用无滑移边界条件,采用壁面函数来处理壁面边界。
3 实例计算及分析
本文计算区域如图1所示,模拟闸门开启过程非恒定水流特性,上下游水深分别为30 m 和10 m ,闸门长6 m ,开启速度为v (假设闸门是匀速开启的),点A 和B 分别是闸后1.5 m 处上下壁面处的压力测点,计算中不考虑上下游水位的变化。计算时间步长为0.01s ,计
中国水利学会第二首届青年科技论坛论文集314 启时间变长后流量的波动明显加剧,这是由于开启速度的减缓,形成了较大尺度的漩涡,对流量形成了明显的影响。不同开启时间在达到全开时的流量也有明显的差异,开启时间短的过程在达到全开时的流量较小,之后流量的增加也较快,而开启时间长的过程在达到
3.2 流态特征分析
图4是开启时间为10s、20s和30s的不同时刻流态图,为了便于比较,3种开启时间列出的相对开度是一致的。不同开启时间在较小开度时,闸门后回流区随开度逐渐增大而逐渐形成、发展,在开度增加的过程中,主流厚度在增加,并不断的将原先处于主流上侧的低速水体加速汇入主流,从而能在启门后一段时间内维持门后漩涡的稳定,开启速度较慢时形成的漩涡较大较稳定。随着闸门的进一步开启,闸门后漩涡出现分裂、脱落,接着
e=时A点和B点压力为静水压,压力分别为和B点的压力随相对开度变化的情况。在0
10×9.8kPa和4×9.8kPa。在闸门开启的瞬间,压力陡然增大,开启速度越快压力增大值越
315 四、水力学河流动力学
e=附件,压力又开始高,随后在相对开度较小时压力有一个逐渐降低的过程,大概在0.5
回升。整个开启过程中压力都存在波动,在小开度时波动较小,这和小开度时流体较稳定时一致的;在大开度时,压力波动周期变小,波动范围变大,开启时间为30s时的情况更
图5 压力分布图
4 结论
本文采用非结构动网格对闸门开启过程进行了垂向二维数值模拟,采用网格变形与局部重构相结合的方法实现动网格,在闸门开启过程中网格质量较好,控制方程的离散采用有限体积法,对时间和空间的离散均采用二阶精度,得到的闸门开启过程中流量、压力等随时间变化的规律,以及在开启过程中闸后漩涡的生成、运动、合并、分裂、脱落等水流特性,在定性上可信的。该算例表明本文采用的方法对于边界运动较复杂的精细模拟问题有比较好的效果,为用数值模拟方法来研究闸门起闭时非恒定水流问题奠定了基础。但是,目前没有得到定量的数据来验证计算结果,对闸门开启过程非恒定水流的分析还是比较初步的。
参考文献
1 苑明顺,余子牛. 船闸输水阀门后底突扩廊道流态特性分析[J]. 水力发电学报,1997,(2) :43-51.
2 文恒,王永利. 闸门开启过程中非恒定流与闸下消能防冲的研究(一)[J]. 内蒙古农牧学院学报,
1998,19(2):94-98.
3 朱仁庆,杨松林,王志东. 闸门开启中水体流动的数值模拟. 华东船舶工业学院学报,1998,12(3):
18-21.
4 杨建明,吴建华. 动网格技术数值模拟挑流冲刷过程. 水动力学研究与进展A,2001,16(2):156-161.
5 Blom FJ, Considerations on the spring analogy[J], Journal of Aircraft,2000,32:647-668.
【作者简介】沙海飞,男,1979年9月, 硕士,2005年5月毕业于河海大学,现在南京水利科学研究院水工所,助理工程师。