胡雋嘉
(上饒市水利電力勘測(cè)設(shè)計(jì)院,江西 上饒 334000)
由于臺(tái)階式溢洪道的能量損失比光滑溢洪道高,因此臺(tái)階式溢洪道因能從水庫(kù)中釋放多余的水而受到極大的關(guān)注。由于能量損失增加,溢洪道底部的消力池長(zhǎng)度比平滑溢洪道短。在階梯溢洪道的情況下,溢洪道底部(稱為消力池)的水躍長(zhǎng)度也比光滑溢洪道短。由于消力池長(zhǎng)度較短,工程師在設(shè)計(jì)時(shí)更偏向使用大量能量耗散的階梯式溢洪道。階梯式溢洪道邊界層厚度主要取決于粗糙度高度和流向位置。在特定位置之后,這種邊界層湍流到達(dá)自由表面;這被稱為表面起始點(diǎn)。如果這種湍流使表面張力喪失能力,空氣開始進(jìn)入水中,這被稱為充氣區(qū)的開始。在該點(diǎn)之后,充氣區(qū)厚度以分散到臺(tái)階的方式增加,直到總流量被充氣。
臺(tái)階式溢洪道摻氣是水流與臺(tái)階發(fā)生碰撞時(shí),水體發(fā)生破碎從而使水流摻氣。夾帶的空氣通過(guò)氣泡輸送,而截留的空氣在連貫水體的上方輸送。夾帶空氣和截留空氣之和稱為總空氣吸收量。由于上游氣蝕風(fēng)險(xiǎn),如果流量超過(guò)一定限度,需要在起始點(diǎn)上游曝氣,階梯式溢洪道中的起始點(diǎn)位置特別重要。
為此,本文利用Fluent 軟件對(duì)不同坡度的臺(tái)階式溢洪道模型進(jìn)行模擬。模擬設(shè)計(jì)對(duì)16 個(gè)不同的溢洪道進(jìn)行建模,以分析坡度、非摻氣水流區(qū)長(zhǎng)度、起始點(diǎn)和流量之間的關(guān)系。為此,采用四種不同的流量值模擬不同坡度的臺(tái)階式溢洪道模型,研究起始點(diǎn)位置和非摻氣水流區(qū)長(zhǎng)度隨流量和坡度的變化規(guī)律。
本文利用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è)階梯式溢洪道模型具有相同的階梯。
在模擬過(guò)程中應(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é)果。
(1)入口邊界
階梯式溢洪道模型的進(jìn)口由一個(gè)速度進(jìn)口組成。階梯式溢洪道的流速根據(jù)流量而變化。階梯式溢洪道的流速取決于流量。流量越大,在入口處施加的速度越大。
(2)出口邊界
臺(tái)階式溢洪道的出口是有壓出口,條件是出口沒(méi)有回水。出口的回水流量值為0。
(3)墻壁
臺(tái)階式溢洪道的墻被視為具有無(wú)滑動(dòng)速度的靜止墻。所以這些被認(rèn)為是靜止的墻。
本文采用流體體積(VOF)和k-ε realizable 模型模擬從緩坡到陡坡的四種不同坡度的臺(tái)階式溢洪道模型,以計(jì)算非摻氣水流區(qū)的起始點(diǎn)位置和長(zhǎng)度,以及非摻氣水流區(qū)的起始點(diǎn)位置和長(zhǎng)度隨臺(tái)階式溢洪道渠道坡度的變化規(guī)律。為此,以不同的流量率(0.625、0.75、0.875 和1)模擬坡度為12.5°、19°、29°和35°的四個(gè)階梯式溢洪道。
臺(tái)階式溢洪道的槽坡在溢洪道施工中起著重要的作用。臺(tái)階式溢洪道的坡度是計(jì)算起始點(diǎn)位置和非摻氣水流區(qū)長(zhǎng)度的重要參數(shù)。隨著溢洪道坡度的增加,起始點(diǎn)位置向上移動(dòng),對(duì)于所有流量,隨著溢洪道坡度的增加,非摻氣水流區(qū)的大小減小。
所有模型的非摻氣水流區(qū)長(zhǎng)度隨渠道坡度、表面粗糙度和弗勞德表面粗糙度的變化情況。臺(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)系,見(jiàn)圖1。由圖1 可知,隨著弗勞德數(shù)Fr的增加,Li增大。對(duì)于Li的最小值,存在弗勞德數(shù)(Fr)最小值的度量。
圖1 不同坡度臺(tái)階式溢洪道流量q 與起始點(diǎn)長(zhǎng)度的關(guān)系
Li/ks是從溢洪道頂部到起始點(diǎn)的標(biāo)準(zhǔn)化距離。它是Li與表面粗糙度ks的比值。歸一化Li隨弗勞德數(shù)的變化見(jiàn)圖2。歸一化的Li隨著表面粗糙度ks的增加而減小。因此歸一化Li與弗勞德數(shù)直接相關(guān),Li/ks隨著弗勞德數(shù)的增加而增加。Li/ks的最小值顯示在最小弗勞德數(shù)處。
圖2 不同坡度臺(tái)階式溢洪道模型歸一化Li 隨弗勞德數(shù)的變化
臨界深度(dc)是特定流量值下能量最小的水流深度,對(duì)于每個(gè)流量值,都有一個(gè)唯一的臨界深度值,可通過(guò)以下公式確定。這意味著隨著流量的增加,臨界深度值增加。隨著臺(tái)階高度的增加,相同流量下的相對(duì)水深減小。相對(duì)流動(dòng)深度與Li的關(guān)系見(jiàn)圖3。隨著相對(duì)水深值的增加,所有臺(tái)階式溢洪道模型的Li均增加。相對(duì)水深,表明在較低的河道坡度值下,相對(duì)水深增加,非摻氣水流區(qū)的跨度也增加。
圖3 不同坡度臺(tái)階式溢洪道起始點(diǎn)長(zhǎng)度隨相對(duì)臨界深度(dc/h)的變化
不同坡度臺(tái)階式溢洪道的相對(duì)臨界深度與標(biāo)準(zhǔn)化起始點(diǎn)(Li/ks)的關(guān)系見(jià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)系
(1)采用VOF 和k-εrealizable 模型模擬從緩坡到陡坡的四種不同坡度的臺(tái)階式溢洪道模型,以計(jì)算非摻氣水流區(qū)的起始點(diǎn)位置和跨度。
(2)臺(tái)階式溢洪道邊坡流速的增加可能會(huì)引起撇流現(xiàn)象。臺(tái)階式溢洪道中的撇渣流態(tài)取決于流量、臺(tái)階高度和臺(tái)階長(zhǎng)度。對(duì)于確定撇渣流態(tài)的所有排放,排放值應(yīng)大于臨界值。
(3)隨著流量的增加,摻氣位置(起始點(diǎn))下移,無(wú)摻氣區(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ū)的長(zhǎng)度Li與臨界水深有直接關(guān)系,隨著臨界水深的增加,臺(tái)階式溢洪道所有坡面渠道的Li均增加。
(7)相對(duì)流動(dòng)深度是臨界深度與臺(tái)階高度的商。相對(duì)流動(dòng)深度表明,在較低的渠道坡度值下,相對(duì)流動(dòng)深度增加,非摻氣流動(dòng)區(qū)的長(zhǎng)度也增加。