臺(tái)階式溢洪道坡度對(duì)摻氣和起始點(diǎn)位置的影響研究

胡雋嘉

（上饒市水利電力勘測(cè)設(shè)計(jì)院,江西 上饒 334000）

1 引言

由于臺(tái)階式溢洪道的能量損失比光滑溢洪道高,因此臺(tái)階式溢洪道因能從水庫中釋放多余的水而受到極大的關(guān)注。由于能量損失增加,溢洪道底部的消力池長度比平滑溢洪道短。在階梯溢洪道的情況下,溢洪道底部（稱為消力池）的水躍長度也比光滑溢洪道短。由于消力池長度較短,工程師在設(shè)計(jì)時(shí)更偏向使用大量能量耗散的階梯式溢洪道。階梯式溢洪道邊界層厚度主要取決于粗糙度高度和流向位置。在特定位置之后,這種邊界層湍流到達(dá)自由表面；這被稱為表面起始點(diǎn)。如果這種湍流使表面張力喪失能力,空氣開始進(jìn)入水中,這被稱為充氣區(qū)的開始。在該點(diǎn)之后,充氣區(qū)厚度以分散到臺(tái)階的方式增加,直到總流量被充氣。

臺(tái)階式溢洪道摻氣是水流與臺(tái)階發(fā)生碰撞時(shí),水體發(fā)生破碎從而使水流摻氣。夾帶的空氣通過氣泡輸送,而截留的空氣在連貫水體的上方輸送。夾帶空氣和截留空氣之和稱為總空氣吸收量。由于上游氣蝕風(fēng)險(xiǎn),如果流量超過一定限度,需要在起始點(diǎn)上游曝氣,階梯式溢洪道中的起始點(diǎn)位置特別重要。

為此,本文利用Fluent 軟件對(duì)不同坡度的臺(tái)階式溢洪道模型進(jìn)行模擬。模擬設(shè)計(jì)對(duì)16 個(gè)不同的溢洪道進(jìn)行建模,以分析坡度、非摻氣水流區(qū)長度、起始點(diǎn)和流量之間的關(guān)系。為此,采用四種不同的流量值模擬不同坡度的臺(tái)階式溢洪道模型,研究起始點(diǎn)位置和非摻氣水流區(qū)長度隨流量和坡度的變化規(guī)律。

2 有限元計(jì)算分析

本文利用Fluent 軟件對(duì)四種臺(tái)階式溢洪道進(jìn)行數(shù)值模擬。階梯式溢洪道的坡度從平緩到陡峭不等。動(dòng)量和湍流動(dòng)能耗散率和壓力交錯(cuò)選項(xiàng)（PRESTO）使用二階逆風(fēng)和一階逆風(fēng)格式。臺(tái)階式溢洪道的數(shù)值模型由三個(gè)臺(tái)階式溢洪道組成,渠道坡度不同（12.5°、19°、29°和35°）。渠道寬度為0.4 m。整個(gè)階梯式溢洪道模型具有相同的階梯。

在模擬過程中應(yīng)用非常精細(xì)的網(wǎng)格以獲得良好的結(jié)果。網(wǎng)格元素的數(shù)量是287529,節(jié)點(diǎn)的數(shù)量是317868。網(wǎng)格劃分軟件采用AnsysICEM。非結(jié)構(gòu)化網(wǎng)格劃分采用高度平滑,最小尺寸為0.001 m。沿臺(tái)階應(yīng)用了十個(gè)充氣層,以在臺(tái)階附近獲得更好的結(jié)果。充氣層的應(yīng)用有助于研究臺(tái)階附近的流動(dòng)結(jié)果。

2.1 邊界條件

（1）入口邊界

階梯式溢洪道模型的進(jìn)口由一個(gè)速度進(jìn)口組成。階梯式溢洪道的流速根據(jù)流量而變化。階梯式溢洪道的流速取決于流量。流量越大,在入口處施加的速度越大。

（2）出口邊界

臺(tái)階式溢洪道的出口是有壓出口,條件是出口沒有回水。出口的回水流量值為0。

（3）墻壁

臺(tái)階式溢洪道的墻被視為具有無滑動(dòng)速度的靜止墻。所以這些被認(rèn)為是靜止的墻。

2.2 計(jì)算工況

本文采用流體體積（VOF）和k-ε realizable 模型模擬從緩坡到陡坡的四種不同坡度的臺(tái)階式溢洪道模型,以計(jì)算非摻氣水流區(qū)的起始點(diǎn)位置和長度,以及非摻氣水流區(qū)的起始點(diǎn)位置和長度隨臺(tái)階式溢洪道渠道坡度的變化規(guī)律。為此,以不同的流量率（0.625、0.75、0.875 和1）模擬坡度為12.5°、19°、29°和35°的四個(gè)階梯式溢洪道。

3 結(jié)果分析

3.1 河道坡度對(duì)起始點(diǎn)位置的影響

臺(tái)階式溢洪道的槽坡在溢洪道施工中起著重要的作用。臺(tái)階式溢洪道的坡度是計(jì)算起始點(diǎn)位置和非摻氣水流區(qū)長度的重要參數(shù)。隨著溢洪道坡度的增加,起始點(diǎn)位置向上移動(dòng),對(duì)于所有流量,隨著溢洪道坡度的增加,非摻氣水流區(qū)的大小減小。

所有模型的非摻氣水流區(qū)長度隨渠道坡度、表面粗糙度和弗勞德表面粗糙度的變化情況。臺(tái)階式溢洪道的表面粗糙度取決于臺(tái)階高度和渠道坡度。臺(tái)階高度和渠道坡度越大,臺(tái)階式溢洪道的表面粗糙度越大。弗勞德數(shù)計(jì)算公式如下：

式中：ks為表面粗糙度；q為單位流量；g為重力加速度。

從溢洪道頂部到起始點(diǎn)的間距（Li）隨流量（q）的變化。隨著單位流量（q）的增加,Li增加,并隨渠道坡度而變化。而在相同的流量值下,Li隨著溝道斜率的增大而減小。這表明,在恒定流量值下,較高的渠道坡度具有較少的非摻氣水流區(qū)。

Li和弗勞德數(shù)（Fr）之間的關(guān)系,見圖1。由圖1 可知,隨著弗勞德數(shù)Fr的增加,Li增大。對(duì)于Li的最小值,存在弗勞德數(shù)（Fr）最小值的度量。

圖1 不同坡度臺(tái)階式溢洪道流量q 與起始點(diǎn)長度的關(guān)系

Li/ks是從溢洪道頂部到起始點(diǎn)的標(biāo)準(zhǔn)化距離。它是Li與表面粗糙度ks的比值。歸一化Li隨弗勞德數(shù)的變化見圖2。歸一化的Li隨著表面粗糙度ks的增加而減小。因此歸一化Li與弗勞德數(shù)直接相關(guān),Li/ks隨著弗勞德數(shù)的增加而增加。Li/ks的最小值顯示在最小弗勞德數(shù)處。

圖2 不同坡度臺(tái)階式溢洪道模型歸一化Li 隨弗勞德數(shù)的變化

3.2 初始長度隨臨界深度的變化

臨界深度（dc）是特定流量值下能量最小的水流深度,對(duì)于每個(gè)流量值,都有一個(gè)唯一的臨界深度值,可通過以下公式確定。這意味著隨著流量的增加,臨界深度值增加。隨著臺(tái)階高度的增加,相同流量下的相對(duì)水深減小。相對(duì)流動(dòng)深度與Li的關(guān)系見圖3。隨著相對(duì)水深值的增加,所有臺(tái)階式溢洪道模型的Li均增加。相對(duì)水深,表明在較低的河道坡度值下,相對(duì)水深增加,非摻氣水流區(qū)的跨度也增加。

圖3 不同坡度臺(tái)階式溢洪道起始點(diǎn)長度隨相對(duì)臨界深度（dc/h）的變化

不同坡度臺(tái)階式溢洪道的相對(duì)臨界深度與標(biāo)準(zhǔn)化起始點(diǎn)（Li/ks）的關(guān)系見圖4。在所有四個(gè)臺(tái)階式溢洪道模型中,歸一化Li隨著相對(duì)水深的增加而增加。歸一化Li的最大值出現(xiàn)在最大相對(duì)臨界深度,最小斜率為12.5。

圖4 不同坡度臺(tái)階式溢洪道的相對(duì)臨界深度與標(biāo)準(zhǔn)化起始點(diǎn)（Li/ks）的關(guān)系

4 結(jié)論

（1）采用VOF 和k-εrealizable 模型模擬從緩坡到陡坡的四種不同坡度的臺(tái)階式溢洪道模型,以計(jì)算非摻氣水流區(qū)的起始點(diǎn)位置和跨度。

（2）臺(tái)階式溢洪道邊坡流速的增加可能會(huì)引起撇流現(xiàn)象。臺(tái)階式溢洪道中的撇渣流態(tài)取決于流量、臺(tái)階高度和臺(tái)階長度。對(duì)于確定撇渣流態(tài)的所有排放,排放值應(yīng)大于臨界值。

（3）隨著流量的增加,摻氣位置（起始點(diǎn)）下移,無摻氣區(qū)Li的跨度隨著流量的增加而增大。

（4）隨著坡度由緩變陡,摻氣開始靠近壩頂,隨著臺(tái)階式溢洪道坡度的增大,摻氣開始點(diǎn)向壩頂移動(dòng)。

（5）歸一化的Li隨著表面粗糙度ks的增加而減小。因此歸一化Li與弗勞德數(shù)有直接關(guān)系,Li/ks隨著弗勞德數(shù)的增加而增加。Li/ks的最小值顯示在最小弗勞德數(shù)處。

（6）臨界深度（dc）是流量增加導(dǎo)致深度增加的水流深度。非摻氣水流區(qū)的長度Li與臨界水深有直接關(guān)系,隨著臨界水深的增加,臺(tái)階式溢洪道所有坡面渠道的Li均增加。

（7）相對(duì)流動(dòng)深度是臨界深度與臺(tái)階高度的商。相對(duì)流動(dòng)深度表明,在較低的渠道坡度值下,相對(duì)流動(dòng)深度增加,非摻氣流動(dòng)區(qū)的長度也增加。