1、用动网格模拟闸门开启过程非恒定水流特性沙海飞 1，2 周辉 1，2 吴时强 1 陈惠玲 1（1.南京水利科学研究院 南京 210029；2.水利部水科学与水工程重点实验室 南京 210029）【摘 要】 本文采用非结构动网格技术，采用网格变形与局部重构相结合的方法实现动网格，网格的变形采用弹簧模型，对闸门开启过程非恒定水流特性进行了数值模拟，得到了不同开启速度时的流量、压力分布以及不同流态的特征与发生、发展过程，结果令人满意。【关键字】 非结构动网格 闸门 非恒定水流 数值模拟0 前言在水工建筑物中，闸门是用于控制水流最为常用的设施。在水坝、船闸等设计过程中，工程技术人员必须了解上下游水流的运
2、动规律。由于闸门前后存在水位差，当水坝、船闸的闸门开启时，门后容易发生水流空化，时常造成空蚀破坏与闸门振动等危害 1。文献2发现灌区水闸在闸门开启过程中，过闸水流为非恒定流，闸门开启的时段内闸下产生的冲坑深度，占一次放水过程产生冲坑深度的一半。可见在闸门开启过程非恒定水流是要引起足够重视的。而目前这方面的研究主要结合具体工程进行缩尺模型试验，闸门开启时的非恒定水流与复杂的漩涡运动增加了试验的难度，一般通过测量壁面压力来推断体型是否合理，对闸门后的内部流态了解很少。随着计算机技术的飞速提高，计算流体力学得到了迅速发展，采用数值模拟闸门开启过程非恒定水流特性成为可能。文献3对假设闸门突然开启的情况
3、进行了数值模拟，这和实际偏差较大。文献1对船闸弧门非恒定流态进行了二维计算，其计算网格是固定的，对弧门处网格进行了特殊处理，和实际也有一定的偏差。文献4用动网格技术对挑流冲刷过程进行了数值模拟，采用的是结构网格，用网格变形来实现动网格，这种方法很难用于闸门开启这种复杂的运动过程。随着非结构动网格技术的发展，以及在军事、仿生等领域的成功应用，使精确模拟闸门开启过程非恒定水流的特性成为可能。本文采用非结构动网格技术，对平板闸门不同开启速度时非恒定水流进行了垂向二维数值模拟，探讨了闸门后流态发生、发展机理。1 动网格技术本文所处理的问题位移较大，所以采用网格变形与局部重构相结合的方法实现动网格。基本
4、做法是，开设一个包围运动边界的窗口，窗口内的网格可以变形以适应边界的运动。以单元面积为判据，当窗口内即变形区出现严重扭曲的网格单元时，删除变形区内的网格单元，然后用阵面推进法重新生成该区网格。流动参数通过线性插值由旧网格映射到新网格。中国水利学会第二首届青年科技论坛论文集 312控制网格变形的弹簧近似模型是将变形区网格看作一个弹簧网络系统，每一条边都认为是一根具有一定倔强系数的弹簧 5，网格点所受的合力为(1)1iNijijfKx其中 为连结节点 、 的弹簧的倔强系数， 是节点 的位置矢量， 是与节点 相连ij ixiiNi的节点数。网格点的受力始终等于初始状态所受的合力，上式可表示节点 的初
5、始受力，则当运动边界上的点移动后，变形区内节点的新位置坐标可以通过迭代求解如下的线性系统得到。(2)1 11222jj jikj sumsuiksjjKxfafaK (3)0ikikiK若 点 不 与 点 相 连若 点 与 点 相 连， ，该系统对角占优，Jacobi 迭代格式为：(4)11i iNNkkijijj jxxfK当边界运动时，下个时刻的固定边界点和主动边界点的位置解已知，因此弹簧系统的边界条件为 Dirichlet 型，上式经过数次迭代就可达到满意的精度。2 控制方程与数值求解方法考虑是不可压缩水流流动问题，基本控制方程为 Reynolds 方程，省略各变量的时间平均项标志“ ”
6、，紊流模型采用 双方程模型，方程如下：k连续方程： (5)0uix动量方程： (6)ijjuxDpfit方程: (7)ktkiikGx方程： (8) 21tiiCpt式中， 为时间； 和 分别为速度分量和坐标分量； 为运动粘性系数； 为修正压tiup力； 为质量力； 为紊动粘性系数；方程中的经验常数 ，if2tck 0.9uc， ， ， 为平均速度梯度引起的紊动能产生项，1.0k.31.4C2.k313 四、水力学 河流动力学；对于雷诺应力项 的封闭问题，这里用 双方程模型来jitkjiuGx ijuk处理。采用有限体积法对方程进行离散，对时间和空间均采用二阶精度格式，空间离散采用中心差分格式
7、，对时间采用二阶精度的离散格式是很又必要的，因为研究的闸门开启过程水流随时间的演变是很重要的。压力校正采用 SIMPLEC 算法，并采用显式校正步法对SIMPLEC 算法进行了显式校正。上下游均是压力边界条件，给定上下游的总压，对于很长管段中闸门区域的局部模拟，可以采用一维计算为二维计算提供边界条件；壁面采用无滑移边界条件，采用壁面函数来处理壁面边界。3 实例计算及分析本文计算区域如图 1 所示，模拟闸门开启过程非恒定水流特性，上下游水深分别为30 m 和 10 m，闸门长 6 m，开启速度为 (假设闸门是匀速开启的)，点 A 和 B 分别是闸v后 1.5 m 处上下壁面处的压力测点，计算中不
8、考虑上下游水位的变化。计算时间步长为0.01s，计算区域的网格每个时间步长调整一次，经过单个时间步长的网格位移比较小，采用网格变形的方法处理，经过数步之后，网格扭曲就比较厉害，进行局部重构，图 2 是开启时间为 10s 的不同时刻的网格情况，网格质量基本较好。共模拟了开启时间为10s、20s、30s 和 60s 四种情况，也就是开启速度分别为 0.6m/s、0.3m/s、0.2m/s 和0.1m/s。 v301064.51B图 1 计算区域示意图图 2 不同时刻网格（开启总时间为 10s）3.1 流量变化规律图 3 是不同开启时间的开启过程流量曲线图，在闸门开启前，其流量都为零，随着闸门的逐渐
9、开启，流量都随之增大，在达到全开之后，流量都有一个平稳的增长过程。不同t=0.00s t=1.67s t=3.33s t=5.00s t=6.67s t=8.33s t=10.00s中国水利学会第二首届青年科技论坛论文集 314开启时间在开启过程中流量的变化有一定的差异，开启时间短的过程基本曾线性增加，开启时间变长后流量的波动明显加剧，这是由于开启速度的减缓，形成了较大尺度的漩涡，对流量形成了明显的影响。不同开启时间在达到全开时的流量也有明显的差异，开启时间短的过程在达到全开时的流量较小，之后流量的增加也较快，而开启时间长的过程在达到全开时的流量较大，之后的增速也较缓，慢慢接近稳态时的值。02
10、04060801001200 10 20 30 40 50 60 70 80 t(s)Q(m3/sm)图 3 开启过程流量曲线图3.2 流态特征分析图 4 是开启时间为 10s、20s 和 30s 的不同时刻流态图，为了便于比较，3 种开启时间列出的相对开度是一致的。不同开启时间在较小开度时，闸门后回流区随开度逐渐增大而逐渐形成、发展，在开度增加的过程中，主流厚度在增加，并不断的将原先处于主流上侧的低速水体加速汇入主流，从而能在启门后一段时间内维持门后漩涡的稳定，开启速度较慢时形成的漩涡较大较稳定。随着闸门的进一步开启，闸门后漩涡出现分裂、脱落，接着新的漩涡生成、运动、合并、分裂、脱落等引起水
11、流较剧烈变动，在门后不断有尺寸较小的次生涡产生，随着开度增大，门后回流区尺寸随开度增加而逐渐减小，漩涡产生分裂的周期也不断缩短。 XY20406080-430-21010t=3.0s,e152.,4.,2t=9.0s,e513.,6t=.s,e7Frame28Aug5shaife_035755-5.2.33t13456780-t=4.5s,e01t3.,6.,2t=13.5s,e049.,6t=21.s,eFrm8Aug05shaife_0- 7212XY3-4-t=.5s,e01.,2t=.5s,e04.,6t=.s,e07.,Frm8Aug25shaife_1046803.9.13257-
12、 523i3s(a) 开启时间 10s (b) 开启时间 20s (c) 开启时间 30s图 4 闸后流态分布图3.3 压力分布特征在闸门开启过程中闸门前后的压力变化较大，图 5 是不同开启时间情况下，闸后 A点和 B 点的压力随相对开度变化的情况。在 时 A 点和 B 点压力为静水压，压力分别e10s 20s 60s30s315 四、水力学 河流动力学为 109.8kPa 和 49.8kPa。在闸门开启的瞬间，压力陡然增大，开启速度越快压力增大值越高，随后在相对开度较小时压力有一个逐渐降低的过程，大概在 附件，压力0.5e又开始回升。整个开启过程中压力都存在波动，在小开度时波动较小，这和小开
13、度时流体较稳定时一致的；在大开度时，压力波动周期变小，波动范围变大，开启时间为 30s 时的情况更突出，在 处出现了负压，随后压力又迅速回升。0.5e-3.00.03.06.09.012.015.018.021.0-0.1 0.1 0.3 0.5 0.7 0.9 1.1闸 门 相 对 开 度 ep(9.8kPa)AB10s30s20s10s30s20s图 5 压力分布图4 结论本文采用非结构动网格对闸门开启过程进行了垂向二维数值模拟，采用网格变形与局部重构相结合的方法实现动网格，在闸门开启过程中网格质量较好，控制方程的离散采用有限体积法，对时间和空间的离散均采用二阶精度，得到的闸门开启过程中流
14、量、压力等随时间变化的规律，以及在开启过程中闸后漩涡的生成、运动、合并、分裂、脱落等水流特性，在定性上可信的。该算例表明本文采用的方法对于边界运动较复杂的精细模拟问题有比较好的效果，为用数值模拟方法来研究闸门起闭时非恒定水流问题奠定了基础。但是，目前没有得到定量的数据来验证计算结果，对闸门开启过程非恒定水流的分析还是比较初步的。参考文献1 苑明顺，余子牛. 船闸输水阀门后底突扩廊道流态特性分析J. 水力发电学报，1997，(2) ：43-512 文恒，王永利. 闸门开启过程中非恒定流与闸下消能防冲的研究(一)J. 内蒙古农牧学院学报，1998，19（2）：94-983 朱仁庆，杨松林，王志东. 闸门开启中水体流动的数值模拟. 华东船舶工业学院学报，1998，12（3）：18-214 杨建明，吴建华. 动网格技术数值模拟挑流冲刷过程. 水动力学研究与进展 A，2001，16（2）：156-1615 Blom FJ, Considerations on the spring analogyJ, Journal of Aircraft，2000，32：647-668【作者简介】 沙海飞，男，1979 年 9 月, 硕士，2005 年 5 月毕业于河海大学，现在南京水利科学研究中国水利学会第二首届青年科技论坛论文集 316院水工所，助理工程师。
