田梦琦，武慧铃，任文杰，周俊宇，冯思芳，毛熹*
四川农业大学水利水电学院，四川 雅安
收稿日期：2017年8月14日；录用日期：2017年8月29日；发布日期：2017年9月6日
采用标准k-ε湍流模型对竖缝式鱼道内部受加墩影响的流场进行二维数值模拟研究。研究的情况分为三种，即不加墩、加墩、改变墩位置，对这三种情况的流场进行数值模拟计算，并应用物理模型试验成果进行对比验证。结果表明：在竖缝式鱼道的主流路径处加墩，可以在一定程度上降低鱼池的流速。 This thesis is to make two dimensions numerical simulation on flow field of vertical slot fishway by standard k-ε turbulence model. The study is divided into three kinds of condition, that is, without plancon, with plancon, changing plancon’s position. The numerical simulation of flow field in these three cases is carried out, and the results of physical model test are compared and verified. The results indicate: adding the plancon to the mainstream path of vertical slot fishway can slow down the velocity in pool of fishway to a certain degree.
田梦琦，武慧铃，任文杰，周俊宇，冯思芳，毛熹*
四川农业大学水利水电学院，四川 雅安
收稿日期：2017年8月14日；录用日期：2017年8月29日；发布日期：2017年9月6日
采用标准k-ε湍流模型对竖缝式鱼道内部受加墩影响的流场进行二维数值模拟研究。研究的情况分为三种，即不加墩、加墩、改变墩位置，对这三种情况的流场进行数值模拟计算，并应用物理模型试验成果进行对比验证。结果表明：在竖缝式鱼道的主流路径处加墩，可以在一定程度上降低鱼池的流速。
关键词 :竖缝式鱼道，数值模拟，流场
Copyright © 2017 by authors and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
竖缝式鱼道是指水槽大部分被挡板拦截仅留一条竖缝过鱼的鱼道，其消能效果较充分且能适应水位较大幅度变化。目前，竖缝式鱼道应用广泛，如，北京上庄新闸鱼道 [
目前国内对竖缝式鱼道的研究日益增多，曹庆磊等分别采用标准k-ε湍流模型和RSM湍流模型对同侧竖缝式鱼道池室内的水流进行数值模拟研究 [
本文以竖缝式鱼道为研究对象，采用标准k-ε湍流模型进行数值模拟，研究加墩对水流流态的影响，并利用物理模型试验成果对数值模拟的结果进行验证。
数值模拟主要通过ANSYS软件进行，并简化为二维模型计算。
采用标准k-ε湍流模型进行数值模拟，该模型的主要控制方程由连续方程、动量方程和k-ε方程组成。
连续方程：
动量方程：
k-ε方程：
其中：
式(1)~(6)中：ρ为水体密度，视为常数，1000 kg/m3；p为时均压强，Pa；μ和μt分别表示水的分子粘性系数和紊动涡粘系数，m2/s；k为紊动能，m2/s2；ε为紊动耗散率，m2/s3；G为紊动能产生项。σk、C1ε、C2ε均为经验常数常数项，该模型中取
采用控制体积法对偏微分方程组进行离散，压力和速度的耦合采用SIMPEL法。计算模拟的进口条件采用速度进口，以试验资料的平均流速作为进口流速，出流采用自由出流，其它为固壁边界，初始流场静止。
网格划分采用以四边形为主的结构化网格(见图1)。时间步长取决于网格的尺寸和流速大小，本文中时间步长取0.002 s。
根据上述方法建立数学模型，对流量为0.263 L/s的相同工况下，三种不同结构型式的竖缝式鱼道进行了数值计算(见表1)。
鱼道共有4级池室，如图2所示，为便于和物理模型试验结果比较，都选取第3级鱼道池室的结果进行比较。
鱼道池室长2.4 m，宽2.0 m，高3.2 m。鱼道池室的具体结构尺寸，如图3所示，以下结构型式2和结构型式3的尺寸与1相同。
图1. 局部网格示意图
| 工况 | 结构型式 | |
|---|---|---|
| 流量0.263 L/s | 1 | 原始竖缝式鱼道，结构型式不做改变 |
| 2 | 在距离竖缝处0.3 m处加一个0.3 m × 0.3 m的墩 | |
| 3 | 在距离竖缝处0.6 m处加一个0.3 m × 0.3 m的墩 | |
表1. 数值模拟的三种结构型式
图2. 鱼道模型整体平面图
图3. 结构型式1的鱼道池室图(单位：米)
对Z = 0.072 m高度的流场进行数值模拟，可以从图4中看出，主流流经竖缝后形成射流，并在主流两侧形成两个比较明显的回流区，回流区的流速很小。鱼道池室内的最大流速为V1 = 1.05 m/s，平均流速为V1,av = 0.80 m/s，可以看出最大流速在竖缝处。
结构型式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的平面图如图8所示，墩与竖缝处的距离为0.6 m。
图4. 结构型式1的数值模拟结果
图5. 添加墩后的鱼道模型平面图
图6. 结构型式2的示意图(单位：米)
对Z = 0.072m高度的流场进行数值模拟，可以从图9中看出，主流被墩分为两条，在墩的分流消能作用下，流速得到一定程度的减小，最大流速为V3 = 0.86 m/s，最大流速同样出现在竖缝处，竖缝处平均流速为V3,av = 0.48 m/s。
图7. 结构型式2的数值模拟结果
图8. 结构型式3的示意图(单位：米)
结构型式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处)的能力较强。
鱼道模型主体为亚克力板，模型尺寸是在结构型式1、2、3的基础上进行相似原理计算得到的，同样将流量为0.263 L/s的工况进行相似原理计算，得到模型试验的工况。在相同工况下，分别进行结构型式1、2、3的试验，利用旋桨流速仪测量流速，将模型试验测得的数据通过相似原理计算还原为原始模型中的数据，与上节的数值计算结果进行对比。
图10分别为结构型式1、2、3在Z = 0.072 m时池室内的流速分布。
结构型式1的竖缝处平均流速为
将数值模拟结果与模型试验结果对比：在相同结构型式下，前者的竖缝处平均流速与后者比较接近，
图9. 结构型式3的数值模拟结果
图10. 结构型式1、2、3的流速示意图；(a) 结构型式1的流速示意图；(b) 结构型式2的流速示意图；(c) 结构型式3的流速示意图
对于结构型式1，数值模拟结果
应用标准k-ε湍流模型对鱼道内部受加墩影响的流场进行数值模拟计算，并利用模型试验的结果对数值模拟计算的结果进行了验证，得出如下结论：
1) 竖缝式鱼道池室内的最大流速一般出现在竖缝处。
2) 加墩对竖缝式鱼道内部流场有较大的影响。
3) 加墩可以减小竖缝式鱼道竖缝处和池室内的流速值。
4) 墩的位置会影响其消能、减小流速值的能力。
四川省教育厅自然科学项目(035Z1994)；大学生创新训练计划四川省级重点项目(04054693, 04070066)。
田梦琦,武慧铃,任文杰,周俊宇,冯思芳,毛 熹. 竖缝式鱼道内部流场受加墩影响的数值模拟研究 Numerical Simulation on Flow Field of Vertical Slot Fishway with Plancon[J]. 土木工程, 2017, 06(05): 459-466. http://dx.doi.org/10.12677/HJCE.2017.65055