1. 引言
竖缝式鱼道是指水槽大部分被挡板拦截仅留一条竖缝过鱼的鱼道,其消能效果较充分且能适应水位较大幅度变化。目前,竖缝式鱼道应用广泛,如,北京上庄新闸鱼道 [1] ,老龙口水利枢纽工程鱼道 [2] ,西藏狮泉河鱼道等。虽然竖缝式鱼道应用广泛,但存在竖缝附近流速相对较大,不适合我国鱼类游泳能力的问题。
目前国内对竖缝式鱼道的研究日益增多,曹庆磊等分别采用标准k-ε湍流模型和RSM湍流模型对同侧竖缝式鱼道池室内的水流进行数值模拟研究 [3] ,郭维东等人对不同的相对竖缝宽度在同侧竖缝式鱼道水流流态产生的影响进行数值模拟研究 [4] ,诸韬等人对双侧竖缝式鱼道的水力特性进行三维数值模拟研究 [5] ,罗小凤等人对竖缝式鱼道导板的长度及导角大小在鱼道内流场产生的影响进行数值模拟研究 [6] ,毛熹等人通过模型实验提出了一种鱼道优化的新结构 [7] ,毛熹等人采用标准k-ε湍流模型对竖缝式鱼道的底孔、坡度在鱼道内流速、流态产生的影响进行数值模拟研究 [8] 。
本文以竖缝式鱼道为研究对象,采用标准k-ε湍流模型进行数值模拟,研究加墩对水流流态的影响,并利用物理模型试验成果对数值模拟的结果进行验证。
2. 数学模型
数值模拟主要通过ANSYS软件进行,并简化为二维模型计算。
2.1. 控制方程
采用标准k-ε湍流模型进行数值模拟,该模型的主要控制方程由连续方程、动量方程和k-ε方程组成。
连续方程:
(1)
动量方程:
(2)
k-ε方程:
(3)
(4)
其中:
(5)
(6)
式(1)~(6)中:ρ为水体密度,视为常数,1000 kg/m3;p为时均压强,Pa;μ和μt分别表示水的分子粘性系数和紊动涡粘系数,m2/s;k为紊动能,m2/s2;ε为紊动耗散率,m2/s3;G为紊动能产生项。σk、C1ε、C2ε均为经验常数常数项,该模型中取
,
,
,
。
2.2. 求解方法、边界条件和划分网格
采用控制体积法对偏微分方程组进行离散,压力和速度的耦合采用SIMPEL法。计算模拟的进口条件采用速度进口,以试验资料的平均流速作为进口流速,出流采用自由出流,其它为固壁边界,初始流场静止。
网格划分采用以四边形为主的结构化网格(见图1)。时间步长取决于网格的尺寸和流速大小,本文中时间步长取0.002 s。
3. 数值模拟
根据上述方法建立数学模型,对流量为0.263 L/s的相同工况下,三种不同结构型式的竖缝式鱼道进行了数值计算(见表1)。
鱼道共有4级池室,如图2所示,为便于和物理模型试验结果比较,都选取第3级鱼道池室的结果进行比较。
3.1. 结构型式1的数值计算
鱼道池室长2.4 m,宽2.0 m,高3.2 m。鱼道池室的具体结构尺寸,如图3所示,以下结构型式2和结构型式3的尺寸与1相同。
![](Images/Table_Tmp.jpg)
Table 1. Numerical simulation of three different kinds of structure
表1. 数值模拟的三种结构型式
![](//html.hanspub.org/file/4-2750399x21_hanspub.png)
Figure 2. The planar graph of fishway model
图2. 鱼道模型整体平面图
![](//html.hanspub.org/file/4-2750399x22_hanspub.png)
Figure 3. The pool diagram of structure 1 (unit: m)
图3. 结构型式1的鱼道池室图(单位:米)
对Z = 0.072 m高度的流场进行数值模拟,可以从图4中看出,主流流经竖缝后形成射流,并在主流两侧形成两个比较明显的回流区,回流区的流速很小。鱼道池室内的最大流速为V1 = 1.05 m/s,平均流速为V1,av = 0.80 m/s,可以看出最大流速在竖缝处。
3.2. 结构型式2的数值计算
结构型式2、3是在结构型式1的基础上添加墩,结构型式2和结构型式3的区别在于墩的位置不同。添加墩后的平面图如图5所示。
在结构型式2中,墩与竖缝处的距离为0.3 m,具体结构如图6所示。
对Z = 0.072 m高度的流场进行数值模拟,可以从图7中看出,在墩的分流消能作用下,竖缝处及鱼池内的流速减小,竖缝处的最大流速为V2 = 0.90 m/s,平均流速为V2,av = 0.46 m/s。主流也被墩分为两条,上面一条路径较长,下面一条路径较短。
3.3. 结构型式3的数值计算
结构型式3的平面图如图8所示,墩与竖缝处的距离为0.6 m。
![](//html.hanspub.org/file/4-2750399x23_hanspub.png)
Figure 4. The numerical simulation result of structure 1
图4. 结构型式1的数值模拟结果
![](//html.hanspub.org/file/4-2750399x24_hanspub.png)
Figure 5. The planar graph of fishway model with plancon
图5. 添加墩后的鱼道模型平面图
![](//html.hanspub.org/file/4-2750399x25_hanspub.png)
Figure 6. The diagram of structure 2 (unit: m)
图6. 结构型式2的示意图(单位:米)
对Z = 0.072m高度的流场进行数值模拟,可以从图9中看出,主流被墩分为两条,在墩的分流消能作用下,流速得到一定程度的减小,最大流速为V3 = 0.86 m/s,最大流速同样出现在竖缝处,竖缝处平均流速为V3,av = 0.48 m/s。
![](//html.hanspub.org/file/4-2750399x26_hanspub.png)
Figure 7. The numerical simulation result of structure 2
图7. 结构型式2的数值模拟结果
![](//html.hanspub.org/file/4-2750399x27_hanspub.png)
Figure 8. The diagram of structure 3 (unit: m)
图8. 结构型式3的示意图(单位:米)
3.4. 数值计算结果分析
结构型式1、2、3的最大流速均出现在竖缝处,最大流速值分别为:V1 = 1.05 m/s、V2 = 0.90 m/s、V3 = 0.86 m/s,平均流速值分别为:V1,av = 0.80 m/s、V2,av = 0.46 m/s、V3,av = 0.48 m/s。简单计算后可知:结构型式2比结构型式1竖缝处流速减小0.34 m/s,减小幅度约为42.5%;结构型式3比结构型式1竖缝处流速减小0.32 m/s,减小幅度约为40%。由此可以得出:结构型式2和结构型式3都可以较大幅度降低竖缝处流速值。结构型式2 (即将墩结构安放在距离竖缝0.3 m处),减小竖缝处流速值的能力比结构型式3 (即将墩结构安放在距离竖缝0.6 m处)的能力较强。
4. 模型试验
鱼道模型主体为亚克力板,模型尺寸是在结构型式1、2、3的基础上进行相似原理计算得到的,同样将流量为0.263 L/s的工况进行相似原理计算,得到模型试验的工况。在相同工况下,分别进行结构型式1、2、3的试验,利用旋桨流速仪测量流速,将模型试验测得的数据通过相似原理计算还原为原始模型中的数据,与上节的数值计算结果进行对比。
4.1. 模型实验结果与分析
图10分别为结构型式1、2、3在Z = 0.072 m时池室内的流速分布。
结构型式1的竖缝处平均流速为
,结构型式2的竖缝处平均流速为
,结构型式3的竖缝处平均流速为
,模型实验中的最大流速也是出现在竖缝处。
4.2. 数值模拟计算结果与模型试验结果对比分析
将数值模拟结果与模型试验结果对比:在相同结构型式下,前者的竖缝处平均流速与后者比较接近,
![](//html.hanspub.org/file/4-2750399x31_hanspub.png)
Figure 9. The numerical simulation result of structure 3
图9. 结构型式3的数值模拟结果
(a)
(b)
(c)
Figure 10. Velocity diagram of structure 1, 2, 3; (a) Velocity diagram of structure 1; (b) Velocity diagram of structure 2; (c) Velocity diagram of structure 3
图10. 结构型式1、2、3的流速示意图;(a) 结构型式1的流速示意图;(b) 结构型式2的流速示意图;(c) 结构型式3的流速示意图
对于结构型式1,数值模拟结果
,模型试验结果
,两者相差不大;对于结构型式2,数值模拟结果
,模型试验结果
,两者十分接近;对于结构型式3,数值模拟结果
,模型试验结果
,两者十分接近。数值模拟的流场与模型试验测得的流场在相同结构型式下也比较接近。试验结果较好地验证了数值模拟计算的结果。
5. 结论
应用标准k-ε湍流模型对鱼道内部受加墩影响的流场进行数值模拟计算,并利用模型试验的结果对数值模拟计算的结果进行了验证,得出如下结论:
1) 竖缝式鱼道池室内的最大流速一般出现在竖缝处。
2) 加墩对竖缝式鱼道内部流场有较大的影响。
3) 加墩可以减小竖缝式鱼道竖缝处和池室内的流速值。
4) 墩的位置会影响其消能、减小流速值的能力。
基金项目
四川省教育厅自然科学项目(035Z1994);大学生创新训练计划四川省级重点项目(04054693, 04070066)。