第33卷第5期Vol.33No.5水利水电科技进展
AdvancesinScienceandTechnologyofWaterResources
2013年9月Sep.2013
DOI:10.3880/j.issn.10067647.2013.05.011
兴隆水利枢纽工程鱼道水力学数值模拟
汪红波1,2,王从锋1,刘德富1,赵摇萍1,向经文1,童摇迪3
(1.三峡大学水利与环境学院,湖北宜昌摇443002;2.广西电力工业勘察设计研究院,广西南宁摇530023;
3.长江水利委员会长江勘测规划设计研究院,湖北武汉摇430010)
摘要:应用RNGk鄄着湍流模型和有限差分法对兴隆水利枢纽工程鱼道的设计模型和另外2种对比模拟模型进行了数值模拟计算,分析了过鱼池的流态、紊动能和紊动能耗散率等指标,并对比分析了不同孔口布置的淹没孔口式横隔板鱼道的水力特性。结果表明:兴隆水利枢纽鱼道具有良好的消能效果,相邻隔板孔口异侧布置比同侧布置消能效果好;可以通过改变淹没孔口式鱼道隔板的孔口位置、数目及尺寸进行体型优化。
关键词:鱼道;横隔板;流态;消能;数值模拟;兴隆水利枢纽工程
中图分类号:S956郾3;TV131.4摇摇摇文献标志码:A摇摇摇文章编号:10067647(2013)05004705
NumericalsimulationofhydrauliccharacteristicsoffishwayinXinglongHydro鄄JunctionProject//WANGHongbo1,2,WANGCongfeng1,LIUDefu1,ZHAOPing1,XIANGJingwen1,TONGDi3(1.CollegeofHydraulicandEnvironmentalEngineering,ChinaThreeGorgesUniversity,Yichang443002,China;2.GuangxiElectricPowerIndustryInvestigationDesignandResearchInstitute,Nanning530023,China;3.ChangjiangInstituteofSurvey,Planning,DesignandResearch,ChangjiangWaterResourcesCommission,Wuhan430010,China)
Abstract:ThedesignmodelandothertworeferencemodelsofthefishwayinXinglongHydro鄄JunctionProjectwerenumericallysimulatedbytheRNGk鄄着turbulencemodelandthefinitedifferencemethod.Theflowpatternandparameters,suchasturbulentkineticenergyandenergydissipationrate,wereanalyzed.Thehydrauliccharacteristicsofthesubmergedorificetransversediaphragmplatefishwaywithdifferentorificelayoutswerecomparedandanalyzed.TheresultsshowthatenergydissipationeffectofthefishwayinXinglongHydro鄄JunctionProjectisgood.Theenergydissipationeffectwiththeoppositearrangementofadjacentbaffles爷orificeisbetterthanthatwiththeipsilateralarrangement.Thesubmergedorificetransversediaphragmplatefishwaycanbeoptimizedbychangingtheposition,amountandsizeoftheorifices.Ithasastrongplasticityandgoodapplicationprospect.
Keywords:fishway;baffle;flowpattern;energydissipation;numericalsimulation;XinglongHydro鄄JunctionProject
河流的连通性遭到破坏,生境片段化,阻隔了鱼类的洄游通道,使洄游及半洄游性鱼类的生存、繁衍受到
摇摇水利水电工程的建设改变了河流的连续性,使
如倾斜式、三角堰等,通过堰流水垫消能,它是针对拥有较强游泳能力的鱼类设计的,适用于水位变幅很小,有足够布置场地的条件;淹没孔口式鱼道在横隔板上开孔,通过孔口水流的对流扩散消能,可以开表孔和底孔,适用于喜欢在中、底层洄游的中、大型鱼类;竖缝式鱼道水流通过隔板的竖缝下泄,类似射流,可利用池内水流的对冲作用消能,一般适用于能适应较复杂流态的大、中型鱼类,且常用于施工期及天然障碍处过鱼,其结构简单,能适应较大变幅的水位,过鱼效率较高,加拿大、美国等众多国家较多采用竖缝式3者组合的鱼道,能够较好地发挥这3类鱼道了这种形式[3鄄4];组合式鱼道是堰流式、淹没孔口式和
类资源的重要工程设施[2]。
能创造良好的水流条件的过鱼通道,是一种保护鱼
严重威胁[1]。鱼道是一种通过特殊结构形式消能,
鱼道按其隔板形式可以分为横隔板式鱼道、丹
尼尔式鱼道、涵洞式鱼道和仿自然式鱼道4种。横隔板式鱼道是一种国内应用最广泛的鱼道,它又可以分为堰流式鱼道、淹没孔口式鱼道、竖缝式鱼道和组合式鱼道4种。堰流式鱼道在水槽内间隔一定距离布置溢流堰式的挡墙,堰顶可以采取不同的形式,
基金项目:水利部公益性行业科研专项(201201030)
作者简介:汪红波(1987—),女,湖北红安人,硕士研究生,主要从事鱼道及其他过鱼设施研究。E鄄mail:[email protected]通信作者:王从锋(1974—),男,山东曹县人,教授,博士,主要从事鱼道与生态混凝土材料研究。E鄄mail:[email protected]
的优点,具有较好的水流条件,过鱼效果良好。
鱼道内水流条件的好坏是鱼道设计成功与否的关键,加拿大阿尔伯特大学的Puertas等[5]对竖缝式鱼道的流场结构进行了试验研究,发现当鱼道各级Ead等[6]和Yagci[7]分别对堰流式和淹没孔口式鱼道的水力特性进行了研究,发现当流量和堰上水头较小时,鱼道内水流为暴跌流,水流经过堰顶跌向池底,在水池中部垂直面上形成逆时针漩涡;当流量和堰上水头较大时,鱼道内水流为连续流,堰上水流大部分在表面流动,水池中部水流在垂直面上形成顺时针漩涡。我国对鱼道水力特性的研究起步较晚,徐体兵等[8鄄9]对竖缝式鱼道的水力特性进行了研究,认为隔板导角对竖缝式鱼道的流态有较明显的影响,[10]的基础上,对兴隆水利枢纽工程淹没孔口式鱼道导角越大主流衰减越快,弯曲程度越大。本文在文献水池的长宽比为10颐8时可以获得较稳定的流态。
维设计模型如图1(a)所示,由3孔a型隔板(图2
设计模型过鱼池隔板间距为2郾6m,休息池隔板间距10级过鱼池布置1级休息池。鱼道的设计流速为0郾4~0郾8m/s,池室设计水深为0郾6~2郾5m。根据已调查的汉江洄游性鱼类、半洄游性鱼类现状,结合兴隆水利枢纽工程特点,主要过鱼对象的过坝时段为每年的5—8月,主要过鱼对象为洄游性鱼类鳗鲡、长颌鲚的亲体和成体,半洄游性鱼类草鱼、青鱼、鲢、鳙、铜鱼、鳤、鳡等的亲体和成体。取A、B2种模型(图1(b)(c))作为设计模型数值模拟计算的参照模型,池室尺寸与设计模型一致,仅隔板形式不同。A模型由3孔a型隔板相邻隔板孔口同侧布置而成,B模型由2孔式隔板(图2(c))相邻隔板孔口异侧布置而成。图2为隔板尺寸详图,3孔式a型隔板和3孔式b型隔板底部小孔位置相同,大孔位30cm伊30cm;2孔式隔板上下异侧布置2个大孔,孔口尺寸均为80cm伊80cm。
置相反,大孔尺寸为100cm伊100cm,小孔尺寸为为5郾2m。过鱼池底坡为2%,休息池底坡为0,间隔
(a))和3孔b型隔板(图2(b))相邻交错布置而成。
的水力特性进行初步研究,深入分析3种不同孔口布置的淹没孔口式鱼道的水力学特性。
1摇物理模型
位于汉江中下游河段湖北省潜江、天门市境内的兴隆水利枢纽工程的开发任务主要是灌溉和航运,兼顾发电。兴隆水利枢纽工程鱼道是一种典型的淹没孔口式鱼道,宽2m,隔板厚0郾2m、高2郾5m,开有3孔,上部布置1个大孔,下部布置2个小孔,其局部三
2摇数学模型
2.1摇控制方程
本文研究的问题属于不可压缩流体的湍流流动,数值模拟选用RNGk鄄着湍流模型,采用VOF模
图1摇3
种模型的局部三维示意图
图2摇隔板尺寸(单位:cm)
型追踪自由液面,控制方程如下:
a.RNGk鄄着湍流模型方程:k方程
鄣(籽k)鄣(籽kui)
+=鄣t鄣xi
化网格,总网格单元数为3821502个,网格单元平均长、宽、高均为6cm。对隔板处网格进行加密后,网格单元x、y、z尺寸最大比例为y颐z=1,x颐y=x颐z=1郾8,A模型和B模型均采用与兴隆水利枢纽工
(1)
程鱼道相同的网格尺寸和网格数量。3种模型的局部网格如图3所示。2.3摇边界及初始条件
(2)
模型的进口采用流量边界,取1m3/s,进口水位
鄣鄣k+G+籽着琢k滋effk
鄣xixj
()
着方程
着为湍动能耗散率;滋eff为扩散系数;Gk为湍动能k1郾39;t为时间;xi、xj为坐标张量。
b.VOF方程:1F+F鄣t
式中:籽为流体密度;k为湍流动能;ui为速度张量;的产生项;C1着、C2着均为经验常数;琢k、琢着为常数
鄣(鄣
FAxu)+R(FAyv)+xy
DIF
SOR
2
鄣鄣着+C1着着G-C籽着琢滋2着
xi着eff鄣xjkkk
(
籽着)鄣(籽着ui)
+=
鄣xi鄣t
)
为4郾26m(模型内控制净水深为2郾40m,以模型底部水平面为基准面,高程为0m,进口底坡最高点高程为1郾86m),紊动能为2郾69伊10-4m2/s2,紊动能耗出口水位为2郾7m;出口流速为0郾208m/s,边墙和底
散率为7郾33伊10-6m2/s3;出口采用流速出口边界,部设置为固壁边界,粗糙度为0郾8mm。3种模型采用相同的边界条件和初始条件。
[
FDIF=
FAxu鄣(
FAzw)+孜鄣zx
F
x
3摇计算结果分析
每种模型均取位于孔口高度中间位置的2个切片进行分析对比,以顺水流方向为x轴正向,鱼道宽度向左延伸方向为y轴正向,垂直向上为z轴正向,如图1所示。面向水流左手侧为左侧,右手侧为右侧。主要分析鱼道内的流态、紊动能和紊动能耗散率。3.1摇流态分析
选取设计模型的一个过鱼池进行流态分析,图4为设计模型过鱼池的上、下层孔口切片的流速和流线分布(图中流线箭头只表示方向,不表示大小,u1为合速度)。此过鱼池上游隔板为3孔a型隔板,下
(3)]=F+F
F+R鄣鄣F+
UAR[鄣鄣x(UA鄣)()鄣y鄣x鄣z
鄣F+孜UAF(4)UA])鄣z(x鄣z
F
y
F
z
F
x
式中:F为流体体积分数;Ax、Ay、Az分别为对应x、y、z方向的微元面积;u、v、w分别为对应x、y、z方向的流体速度;孜为坐标系数,采用直角坐标时,孜=0,采体积分数密度源项;R为系数,采用直角坐标时为常数;UF为扩散系数;CP为常数,是紊流施密特数2.2摇计算区域及网格划分的倒数;滋是动能消散系数。
用圆柱坐标时,孜=1;FDIF为体积分数扩散项;FSOR为
C滋
UF=p
籽
(5)
游隔板为3孔b型隔板。由图4可知,在大孔层,水流由上游隔板大孔流入,由下游隔板大孔流出,主流呈“L冶形分布,忽略孔口边界处对鱼类影响不大的流速,最大流速为1郾2m/s,产生在池室主流中间段,并0郾75m/s,产生在左侧边墙处。大孔层流速大于在此过鱼池的左侧形成回流区,回流区最大流速为
所有数学模型计算区域均选取鱼道直段的一部分,全长88郾4m,含30级过鱼池和2级休息池。用有限差分法对模型进行离散求解,能大量减少网格数量,用VOF方法追踪自由液面。采用六面体结构
0郾8m/s的区域面积占整个池室面积的23郾0%。在小涡,回流区最大流速为0郾6m/s,右侧小孔水流受上方
孔层,上游隔板左侧小孔水流向左偏转形成回流漩
图3摇3种模型局部网格
大孔水流影响产生分流现象,向左偏转并分为2股,分别流向下游隔板的2小孔,小孔的流速云图呈带状分布,最大流速产生在孔口处。上游隔板左、右2小孔的最大流速分别为1郾35m/s、1郾20m/s;下游隔板20郾8m/s的区域面积占整个池室面积的11郾5%。小孔的最大流速均为1郾35m/s。小孔层流大于
另一侧为顺时针漩涡区,流速不大于0郾30m/s。小孔层分成4个流速区,从左至右依次为漩涡区、直流区、漩涡区、直流区。漩涡区流速不大于0郾15m/s,左侧的小孔水流受大孔水流影响存在分流现象,分成2股分别流向下游隔板2小孔,最大流速产生在孔口处,上游隔板2小孔的最大流速均为0郾90m/s,下游隔板
图4摇设计模型流速和流线
图5为A模型上、下层孔口切片的流速和流线分布。由图5可知,大孔层可以分成2个流速区,大孔侧为直流区,流速大于0郾80m/s的区域面积为整个池室面积的1/2,流速在0郾75~1郾20m/s范围内。
图5摇A模型流速和流线
20郾小孔的最大流速均为80m图/s0郾75m/s。小孔层流速大于6的区域面积占整个池室面积的为B模型上、下层孔口切片的流速和流线11郾5%。
分布。由图6可知,上孔层主流呈“L冶形分布,忽略1郾孔口边界处对鱼类影响不大的流速20m/s,产生在池室右侧中间主流处,0郾30m/,最大流速为成回流区,回流区最大流速为s。并在左侧形上孔层流速大于0郾80m/s的区域面积占整个池室面积的忽略孔口边界处对鱼类影响不大的流速31郾0%。下孔层主流亦呈“L冶形分布,与上孔层类似大于05m0郾/s,回流区最大流速为0郾30m/s。,下孔层流速最大流速为
,1郾15郾4%。
80m/s的区域面积占整个池室面积的
图6摇B模型流速和流线
由以上分析可知,设计模型和A模型相比,就大流速区(流速大于0郾80m/s)的面积而言,相邻隔板孔口异侧布置(设计模型)比同侧布置(A模型)小,A模型的大流速区面积为整个池室的1/2,不利于鱼类的上溯;从最大流速方面看,大孔层流速相差不大,小孔层设计模型的最大流速大于A模型的最大流速。B模型上孔层大流速区面积较比A模型小,比设计模型稍大;下孔层大流速区面积与设计模
型相差不大,但设计模型下孔层上游隔板左、右2小
2小孔的最大流速均为1郾35m/s,不利于鱼类上溯。3.2摇消能效果对比因此,B模型流态较好。
孔的最大流速分别为1郾35m/s、1郾20m/s,下游隔板布置大孔比布置小孔流态好;堰流比淹没出流流态稳定。参考文献:
[1]刘志雄,周赤,黄明海.鱼道应用现状和研究进展[J].
各种水域中广泛存在不同程度的紊流流态(如各种漩涡、小尺度紊动、回流等),紊流剪切应力会对鱼类造成各种各样的伤害,如击坏眼睛,扭伤身体,使鱼类迷失方向,导致平衡能力下降,不易躲避其他肉食动物的捕食等[10],因此,紊动能及紊动能
ZHOUChi,HUANGMinghai.Situationanddevelopment(inChinese))
RiverScientificResearchInstitute,2010,27(4):28鄄31.
长江科学院院报,2010,27(4):28鄄31.(LIUZhixiong,
offishwayresearchandapplication[J].JournalofYangtze
耗散率是鱼道水力特性研究的重要指标。紊动能是反映水流因紊动而引起的能量损失的重要指标,紊动能耗散率是紊动能转化成分子热运动动能的速咬[2]KAMULAR,BARTHELJ.Effectsofmodificationsontheresearch,2000,5(1):67鄄79.
hydraulicsofDenilfishways[J].Borealenvironment
率,紊动能不变,紊动能耗散率越大,消能效果越好。由表1可知,兴隆水利枢纽工程鱼道设计模型与A模型相比0郾倍多094,0但紊动能耗散率比异侧布置小~,0郾上层孔口同侧布置1150m是异侧布置(A模型()紊动能为
2/s2,;设计模型下孔层2)种指的9
标相差不大,说明相邻隔板孔口异侧布置比同侧布
置消能效果好。A模型与B模型相比,上孔层紊动能大小相近,但紊动能耗散率相差较大,B模型的紊动能耗散率约是A模型的2倍,B模型下层孔口水流受上层孔口水流影响,紊动能比上层孔口大,紊动能耗散率也大。
表1摇3种模型紊动能及紊动能耗散率最大值
模型切片位置0郾0105~0郾0125(m2·s-3/设计模型上层孔口(m紊动能2·s-2/
)
紊动能耗散率0郾070~0郾090)下层孔口0郾0080~0郾01400郾040~0郾080A模型上层孔口0郾0940~0郾11500郾050~0郾060下层孔口0郾0500~0郾0700
0郾040~0郾100B模型
上层孔口0郾09500郾080~0郾141下层孔口
0郾1150~0郾1910
0郾111~0郾191
4摇结摇论
a.兴隆水利枢纽工程鱼道结构设计较优,选择了较好的孔口尺寸和孔口位置,使紊动能较小,紊动能耗散率比0郾b.A模A型模型大流态较,消能效果良好顺直。
利于鱼类上溯8m/s的区域面积约占整个过鱼池面积的,但大孔层流速。小孔层大流速区面积也比兴隆水利1/大2,于
不枢纽工程鱼道小孔层的大c.从流速方面看,B。
模型下孔层优于设计模型;从消能效果上看d.相邻隔板孔口异侧布置比同侧布置消能效
,B模型上孔层优于设计模型。
果好,大流速区(>0郾80m/s)面积小,可以减小鱼类上溯的流速障碍;上、下层孔口水流相互影响,下层
[3]RAJARATNAMN,KATOPODISC.HydraulicsofDenil[4]王兴勇(9):1219鄄1233.
fishways[J].JournalofHydraulicEngineering,1984,110电科学,研郭军究.院国内外鱼道研究与建设学报,2005,3(3):222鄄228.[J].中国水利水
(WANGconstructionXingyong,GUOoffishwaysJun.atBriefhomereviewandabroad[J].onresearchJournalandof[5]PUERTASResearch,2005,3(3):222鄄228.ChinaInstituteofWaterResourcesJ,PENAL,TEIJEIRO(inChinese))
andHydropower
ofverticalslotT.fishwaysExperimental[J].
[6]EADCanadianapproachtoMetallurgicalthehydraulicsQuarterly,2004,130(1):10鄄23.poolSandA,weirKATOPODISfishwaysC.[JFlow].Journalregimesofandstructurein
[7]EngineeringYAGCIO.HydraulicandScience,aspects2004.of3(5):pool鄄weir379鄄390.Environmentalecologicallyfriendlywaterstructure[J].fishwaysEcologicalas[8]徐体兵Engineering,水利学报,孙双科2010,,2009,40.竖缝式鱼道水流结构的数值模拟36(1):36鄄46.
(11):1386鄄1391.(XUTibing,SUN[J].
vertical
Shuangke.slot
Numericalfishwaysimulation[J].
ofJournal
theflowof
structureEngineering,2009,40(11):1386鄄1391.(inChinese))
Hydraulic
in
[9]罗小凤江科学院院报,李嘉.竖缝式鱼道结构及水力特性研究,2010,27(10):50鄄54.(LUOXiaofeng,LI[J].长
Jia.vertical鄄slotStudyonstructureandhydrauliccharacteristicsofScientificChinese))
ResearchfishwayInstitute,2010,27[J].Journalof(10Yangtze):50鄄54.River
(in[10]汪红波值模拟研究,王从锋[J].,刘德富水电能源科学,等.横隔板式鱼道水力特性数
,2012,30(5):65鄄68.on(WANGtransversenumericalHongbo,WANGdiaphragmsimulationCongfeng,LIUplateoffishwayhydraulic[J].characteristicsDefu.etal.StudyWaterResources
ofandPower,2012,30(5):65鄄68.(收稿日期:2012(in11Chinese))
20摇编辑:熊水斌)
第33卷第5期Vol.33No.5水利水电科技进展
AdvancesinScienceandTechnologyofWaterResources
2013年9月Sep.2013
DOI:10.3880/j.issn.10067647.2013.05.011
兴隆水利枢纽工程鱼道水力学数值模拟
汪红波1,2,王从锋1,刘德富1,赵摇萍1,向经文1,童摇迪3
(1.三峡大学水利与环境学院,湖北宜昌摇443002;2.广西电力工业勘察设计研究院,广西南宁摇530023;
3.长江水利委员会长江勘测规划设计研究院,湖北武汉摇430010)
摘要:应用RNGk鄄着湍流模型和有限差分法对兴隆水利枢纽工程鱼道的设计模型和另外2种对比模拟模型进行了数值模拟计算,分析了过鱼池的流态、紊动能和紊动能耗散率等指标,并对比分析了不同孔口布置的淹没孔口式横隔板鱼道的水力特性。结果表明:兴隆水利枢纽鱼道具有良好的消能效果,相邻隔板孔口异侧布置比同侧布置消能效果好;可以通过改变淹没孔口式鱼道隔板的孔口位置、数目及尺寸进行体型优化。
关键词:鱼道;横隔板;流态;消能;数值模拟;兴隆水利枢纽工程
中图分类号:S956郾3;TV131.4摇摇摇文献标志码:A摇摇摇文章编号:10067647(2013)05004705
NumericalsimulationofhydrauliccharacteristicsoffishwayinXinglongHydro鄄JunctionProject//WANGHongbo1,2,WANGCongfeng1,LIUDefu1,ZHAOPing1,XIANGJingwen1,TONGDi3(1.CollegeofHydraulicandEnvironmentalEngineering,ChinaThreeGorgesUniversity,Yichang443002,China;2.GuangxiElectricPowerIndustryInvestigationDesignandResearchInstitute,Nanning530023,China;3.ChangjiangInstituteofSurvey,Planning,DesignandResearch,ChangjiangWaterResourcesCommission,Wuhan430010,China)
Abstract:ThedesignmodelandothertworeferencemodelsofthefishwayinXinglongHydro鄄JunctionProjectwerenumericallysimulatedbytheRNGk鄄着turbulencemodelandthefinitedifferencemethod.Theflowpatternandparameters,suchasturbulentkineticenergyandenergydissipationrate,wereanalyzed.Thehydrauliccharacteristicsofthesubmergedorificetransversediaphragmplatefishwaywithdifferentorificelayoutswerecomparedandanalyzed.TheresultsshowthatenergydissipationeffectofthefishwayinXinglongHydro鄄JunctionProjectisgood.Theenergydissipationeffectwiththeoppositearrangementofadjacentbaffles爷orificeisbetterthanthatwiththeipsilateralarrangement.Thesubmergedorificetransversediaphragmplatefishwaycanbeoptimizedbychangingtheposition,amountandsizeoftheorifices.Ithasastrongplasticityandgoodapplicationprospect.
Keywords:fishway;baffle;flowpattern;energydissipation;numericalsimulation;XinglongHydro鄄JunctionProject
河流的连通性遭到破坏,生境片段化,阻隔了鱼类的洄游通道,使洄游及半洄游性鱼类的生存、繁衍受到
摇摇水利水电工程的建设改变了河流的连续性,使
如倾斜式、三角堰等,通过堰流水垫消能,它是针对拥有较强游泳能力的鱼类设计的,适用于水位变幅很小,有足够布置场地的条件;淹没孔口式鱼道在横隔板上开孔,通过孔口水流的对流扩散消能,可以开表孔和底孔,适用于喜欢在中、底层洄游的中、大型鱼类;竖缝式鱼道水流通过隔板的竖缝下泄,类似射流,可利用池内水流的对冲作用消能,一般适用于能适应较复杂流态的大、中型鱼类,且常用于施工期及天然障碍处过鱼,其结构简单,能适应较大变幅的水位,过鱼效率较高,加拿大、美国等众多国家较多采用竖缝式3者组合的鱼道,能够较好地发挥这3类鱼道了这种形式[3鄄4];组合式鱼道是堰流式、淹没孔口式和
类资源的重要工程设施[2]。
能创造良好的水流条件的过鱼通道,是一种保护鱼
严重威胁[1]。鱼道是一种通过特殊结构形式消能,
鱼道按其隔板形式可以分为横隔板式鱼道、丹
尼尔式鱼道、涵洞式鱼道和仿自然式鱼道4种。横隔板式鱼道是一种国内应用最广泛的鱼道,它又可以分为堰流式鱼道、淹没孔口式鱼道、竖缝式鱼道和组合式鱼道4种。堰流式鱼道在水槽内间隔一定距离布置溢流堰式的挡墙,堰顶可以采取不同的形式,
基金项目:水利部公益性行业科研专项(201201030)
作者简介:汪红波(1987—),女,湖北红安人,硕士研究生,主要从事鱼道及其他过鱼设施研究。E鄄mail:[email protected]通信作者:王从锋(1974—),男,山东曹县人,教授,博士,主要从事鱼道与生态混凝土材料研究。E鄄mail:[email protected]
的优点,具有较好的水流条件,过鱼效果良好。
鱼道内水流条件的好坏是鱼道设计成功与否的关键,加拿大阿尔伯特大学的Puertas等[5]对竖缝式鱼道的流场结构进行了试验研究,发现当鱼道各级Ead等[6]和Yagci[7]分别对堰流式和淹没孔口式鱼道的水力特性进行了研究,发现当流量和堰上水头较小时,鱼道内水流为暴跌流,水流经过堰顶跌向池底,在水池中部垂直面上形成逆时针漩涡;当流量和堰上水头较大时,鱼道内水流为连续流,堰上水流大部分在表面流动,水池中部水流在垂直面上形成顺时针漩涡。我国对鱼道水力特性的研究起步较晚,徐体兵等[8鄄9]对竖缝式鱼道的水力特性进行了研究,认为隔板导角对竖缝式鱼道的流态有较明显的影响,[10]的基础上,对兴隆水利枢纽工程淹没孔口式鱼道导角越大主流衰减越快,弯曲程度越大。本文在文献水池的长宽比为10颐8时可以获得较稳定的流态。
维设计模型如图1(a)所示,由3孔a型隔板(图2
设计模型过鱼池隔板间距为2郾6m,休息池隔板间距10级过鱼池布置1级休息池。鱼道的设计流速为0郾4~0郾8m/s,池室设计水深为0郾6~2郾5m。根据已调查的汉江洄游性鱼类、半洄游性鱼类现状,结合兴隆水利枢纽工程特点,主要过鱼对象的过坝时段为每年的5—8月,主要过鱼对象为洄游性鱼类鳗鲡、长颌鲚的亲体和成体,半洄游性鱼类草鱼、青鱼、鲢、鳙、铜鱼、鳤、鳡等的亲体和成体。取A、B2种模型(图1(b)(c))作为设计模型数值模拟计算的参照模型,池室尺寸与设计模型一致,仅隔板形式不同。A模型由3孔a型隔板相邻隔板孔口同侧布置而成,B模型由2孔式隔板(图2(c))相邻隔板孔口异侧布置而成。图2为隔板尺寸详图,3孔式a型隔板和3孔式b型隔板底部小孔位置相同,大孔位30cm伊30cm;2孔式隔板上下异侧布置2个大孔,孔口尺寸均为80cm伊80cm。
置相反,大孔尺寸为100cm伊100cm,小孔尺寸为为5郾2m。过鱼池底坡为2%,休息池底坡为0,间隔
(a))和3孔b型隔板(图2(b))相邻交错布置而成。
的水力特性进行初步研究,深入分析3种不同孔口布置的淹没孔口式鱼道的水力学特性。
1摇物理模型
位于汉江中下游河段湖北省潜江、天门市境内的兴隆水利枢纽工程的开发任务主要是灌溉和航运,兼顾发电。兴隆水利枢纽工程鱼道是一种典型的淹没孔口式鱼道,宽2m,隔板厚0郾2m、高2郾5m,开有3孔,上部布置1个大孔,下部布置2个小孔,其局部三
2摇数学模型
2.1摇控制方程
本文研究的问题属于不可压缩流体的湍流流动,数值模拟选用RNGk鄄着湍流模型,采用VOF模
图1摇3
种模型的局部三维示意图
图2摇隔板尺寸(单位:cm)
型追踪自由液面,控制方程如下:
a.RNGk鄄着湍流模型方程:k方程
鄣(籽k)鄣(籽kui)
+=鄣t鄣xi
化网格,总网格单元数为3821502个,网格单元平均长、宽、高均为6cm。对隔板处网格进行加密后,网格单元x、y、z尺寸最大比例为y颐z=1,x颐y=x颐z=1郾8,A模型和B模型均采用与兴隆水利枢纽工
(1)
程鱼道相同的网格尺寸和网格数量。3种模型的局部网格如图3所示。2.3摇边界及初始条件
(2)
模型的进口采用流量边界,取1m3/s,进口水位
鄣鄣k+G+籽着琢k滋effk
鄣xixj
()
着方程
着为湍动能耗散率;滋eff为扩散系数;Gk为湍动能k1郾39;t为时间;xi、xj为坐标张量。
b.VOF方程:1F+F鄣t
式中:籽为流体密度;k为湍流动能;ui为速度张量;的产生项;C1着、C2着均为经验常数;琢k、琢着为常数
鄣(鄣
FAxu)+R(FAyv)+xy
DIF
SOR
2
鄣鄣着+C1着着G-C籽着琢滋2着
xi着eff鄣xjkkk
(
籽着)鄣(籽着ui)
+=
鄣xi鄣t
)
为4郾26m(模型内控制净水深为2郾40m,以模型底部水平面为基准面,高程为0m,进口底坡最高点高程为1郾86m),紊动能为2郾69伊10-4m2/s2,紊动能耗出口水位为2郾7m;出口流速为0郾208m/s,边墙和底
散率为7郾33伊10-6m2/s3;出口采用流速出口边界,部设置为固壁边界,粗糙度为0郾8mm。3种模型采用相同的边界条件和初始条件。
[
FDIF=
FAxu鄣(
FAzw)+孜鄣zx
F
x
3摇计算结果分析
每种模型均取位于孔口高度中间位置的2个切片进行分析对比,以顺水流方向为x轴正向,鱼道宽度向左延伸方向为y轴正向,垂直向上为z轴正向,如图1所示。面向水流左手侧为左侧,右手侧为右侧。主要分析鱼道内的流态、紊动能和紊动能耗散率。3.1摇流态分析
选取设计模型的一个过鱼池进行流态分析,图4为设计模型过鱼池的上、下层孔口切片的流速和流线分布(图中流线箭头只表示方向,不表示大小,u1为合速度)。此过鱼池上游隔板为3孔a型隔板,下
(3)]=F+F
F+R鄣鄣F+
UAR[鄣鄣x(UA鄣)()鄣y鄣x鄣z
鄣F+孜UAF(4)UA])鄣z(x鄣z
F
y
F
z
F
x
式中:F为流体体积分数;Ax、Ay、Az分别为对应x、y、z方向的微元面积;u、v、w分别为对应x、y、z方向的流体速度;孜为坐标系数,采用直角坐标时,孜=0,采体积分数密度源项;R为系数,采用直角坐标时为常数;UF为扩散系数;CP为常数,是紊流施密特数2.2摇计算区域及网格划分的倒数;滋是动能消散系数。
用圆柱坐标时,孜=1;FDIF为体积分数扩散项;FSOR为
C滋
UF=p
籽
(5)
游隔板为3孔b型隔板。由图4可知,在大孔层,水流由上游隔板大孔流入,由下游隔板大孔流出,主流呈“L冶形分布,忽略孔口边界处对鱼类影响不大的流速,最大流速为1郾2m/s,产生在池室主流中间段,并0郾75m/s,产生在左侧边墙处。大孔层流速大于在此过鱼池的左侧形成回流区,回流区最大流速为
所有数学模型计算区域均选取鱼道直段的一部分,全长88郾4m,含30级过鱼池和2级休息池。用有限差分法对模型进行离散求解,能大量减少网格数量,用VOF方法追踪自由液面。采用六面体结构
0郾8m/s的区域面积占整个池室面积的23郾0%。在小涡,回流区最大流速为0郾6m/s,右侧小孔水流受上方
孔层,上游隔板左侧小孔水流向左偏转形成回流漩
图3摇3种模型局部网格
大孔水流影响产生分流现象,向左偏转并分为2股,分别流向下游隔板的2小孔,小孔的流速云图呈带状分布,最大流速产生在孔口处。上游隔板左、右2小孔的最大流速分别为1郾35m/s、1郾20m/s;下游隔板20郾8m/s的区域面积占整个池室面积的11郾5%。小孔的最大流速均为1郾35m/s。小孔层流大于
另一侧为顺时针漩涡区,流速不大于0郾30m/s。小孔层分成4个流速区,从左至右依次为漩涡区、直流区、漩涡区、直流区。漩涡区流速不大于0郾15m/s,左侧的小孔水流受大孔水流影响存在分流现象,分成2股分别流向下游隔板2小孔,最大流速产生在孔口处,上游隔板2小孔的最大流速均为0郾90m/s,下游隔板
图4摇设计模型流速和流线
图5为A模型上、下层孔口切片的流速和流线分布。由图5可知,大孔层可以分成2个流速区,大孔侧为直流区,流速大于0郾80m/s的区域面积为整个池室面积的1/2,流速在0郾75~1郾20m/s范围内。
图5摇A模型流速和流线
20郾小孔的最大流速均为80m图/s0郾75m/s。小孔层流速大于6的区域面积占整个池室面积的为B模型上、下层孔口切片的流速和流线11郾5%。
分布。由图6可知,上孔层主流呈“L冶形分布,忽略1郾孔口边界处对鱼类影响不大的流速20m/s,产生在池室右侧中间主流处,0郾30m/,最大流速为成回流区,回流区最大流速为s。并在左侧形上孔层流速大于0郾80m/s的区域面积占整个池室面积的忽略孔口边界处对鱼类影响不大的流速31郾0%。下孔层主流亦呈“L冶形分布,与上孔层类似大于05m0郾/s,回流区最大流速为0郾30m/s。,下孔层流速最大流速为
,1郾15郾4%。
80m/s的区域面积占整个池室面积的
图6摇B模型流速和流线
由以上分析可知,设计模型和A模型相比,就大流速区(流速大于0郾80m/s)的面积而言,相邻隔板孔口异侧布置(设计模型)比同侧布置(A模型)小,A模型的大流速区面积为整个池室的1/2,不利于鱼类的上溯;从最大流速方面看,大孔层流速相差不大,小孔层设计模型的最大流速大于A模型的最大流速。B模型上孔层大流速区面积较比A模型小,比设计模型稍大;下孔层大流速区面积与设计模
型相差不大,但设计模型下孔层上游隔板左、右2小
2小孔的最大流速均为1郾35m/s,不利于鱼类上溯。3.2摇消能效果对比因此,B模型流态较好。
孔的最大流速分别为1郾35m/s、1郾20m/s,下游隔板布置大孔比布置小孔流态好;堰流比淹没出流流态稳定。参考文献:
[1]刘志雄,周赤,黄明海.鱼道应用现状和研究进展[J].
各种水域中广泛存在不同程度的紊流流态(如各种漩涡、小尺度紊动、回流等),紊流剪切应力会对鱼类造成各种各样的伤害,如击坏眼睛,扭伤身体,使鱼类迷失方向,导致平衡能力下降,不易躲避其他肉食动物的捕食等[10],因此,紊动能及紊动能
ZHOUChi,HUANGMinghai.Situationanddevelopment(inChinese))
RiverScientificResearchInstitute,2010,27(4):28鄄31.
长江科学院院报,2010,27(4):28鄄31.(LIUZhixiong,
offishwayresearchandapplication[J].JournalofYangtze
耗散率是鱼道水力特性研究的重要指标。紊动能是反映水流因紊动而引起的能量损失的重要指标,紊动能耗散率是紊动能转化成分子热运动动能的速咬[2]KAMULAR,BARTHELJ.Effectsofmodificationsontheresearch,2000,5(1):67鄄79.
hydraulicsofDenilfishways[J].Borealenvironment
率,紊动能不变,紊动能耗散率越大,消能效果越好。由表1可知,兴隆水利枢纽工程鱼道设计模型与A模型相比0郾倍多094,0但紊动能耗散率比异侧布置小~,0郾上层孔口同侧布置1150m是异侧布置(A模型()紊动能为
2/s2,;设计模型下孔层2)种指的9
标相差不大,说明相邻隔板孔口异侧布置比同侧布
置消能效果好。A模型与B模型相比,上孔层紊动能大小相近,但紊动能耗散率相差较大,B模型的紊动能耗散率约是A模型的2倍,B模型下层孔口水流受上层孔口水流影响,紊动能比上层孔口大,紊动能耗散率也大。
表1摇3种模型紊动能及紊动能耗散率最大值
模型切片位置0郾0105~0郾0125(m2·s-3/设计模型上层孔口(m紊动能2·s-2/
)
紊动能耗散率0郾070~0郾090)下层孔口0郾0080~0郾01400郾040~0郾080A模型上层孔口0郾0940~0郾11500郾050~0郾060下层孔口0郾0500~0郾0700
0郾040~0郾100B模型
上层孔口0郾09500郾080~0郾141下层孔口
0郾1150~0郾1910
0郾111~0郾191
4摇结摇论
a.兴隆水利枢纽工程鱼道结构设计较优,选择了较好的孔口尺寸和孔口位置,使紊动能较小,紊动能耗散率比0郾b.A模A型模型大流态较,消能效果良好顺直。
利于鱼类上溯8m/s的区域面积约占整个过鱼池面积的,但大孔层流速。小孔层大流速区面积也比兴隆水利1/大2,于
不枢纽工程鱼道小孔层的大c.从流速方面看,B。
模型下孔层优于设计模型;从消能效果上看d.相邻隔板孔口异侧布置比同侧布置消能效
,B模型上孔层优于设计模型。
果好,大流速区(>0郾80m/s)面积小,可以减小鱼类上溯的流速障碍;上、下层孔口水流相互影响,下层
[3]RAJARATNAMN,KATOPODISC.HydraulicsofDenil[4]王兴勇(9):1219鄄1233.
fishways[J].JournalofHydraulicEngineering,1984,110电科学,研郭军究.院国内外鱼道研究与建设学报,2005,3(3):222鄄228.[J].中国水利水
(WANGconstructionXingyong,GUOoffishwaysJun.atBriefhomereviewandabroad[J].onresearchJournalandof[5]PUERTASResearch,2005,3(3):222鄄228.ChinaInstituteofWaterResourcesJ,PENAL,TEIJEIRO(inChinese))
andHydropower
ofverticalslotT.fishwaysExperimental[J].
[6]EADCanadianapproachtoMetallurgicalthehydraulicsQuarterly,2004,130(1):10鄄23.poolSandA,weirKATOPODISfishwaysC.[JFlow].Journalregimesofandstructurein
[7]EngineeringYAGCIO.HydraulicandScience,aspects2004.of3(5):pool鄄weir379鄄390.Environmentalecologicallyfriendlywaterstructure[J].fishwaysEcologicalas[8]徐体兵Engineering,水利学报,孙双科2010,,2009,40.竖缝式鱼道水流结构的数值模拟36(1):36鄄46.
(11):1386鄄1391.(XUTibing,SUN[J].
vertical
Shuangke.slot
Numericalfishwaysimulation[J].
ofJournal
theflowof
structureEngineering,2009,40(11):1386鄄1391.(inChinese))
Hydraulic
in
[9]罗小凤江科学院院报,李嘉.竖缝式鱼道结构及水力特性研究,2010,27(10):50鄄54.(LUOXiaofeng,LI[J].长
Jia.vertical鄄slotStudyonstructureandhydrauliccharacteristicsofScientificChinese))
ResearchfishwayInstitute,2010,27[J].Journalof(10Yangtze):50鄄54.River
(in[10]汪红波值模拟研究,王从锋[J].,刘德富水电能源科学,等.横隔板式鱼道水力特性数
,2012,30(5):65鄄68.on(WANGtransversenumericalHongbo,WANGdiaphragmsimulationCongfeng,LIUplateoffishwayhydraulic[J].characteristicsDefu.etal.StudyWaterResources
ofandPower,2012,30(5):65鄄68.(收稿日期:2012(in11Chinese))
20摇编辑:熊水斌)