朱偉國(guó)
(江西久源建設(shè)工程有限公司,南昌 330038)
溢洪道臺(tái)階水平面末端增設(shè)尾坎的坎式階梯溢洪道在水利水電工程中得到廣泛應(yīng)用,但相關(guān)研究仍較少,現(xiàn)有的研究主要集中在摻氣特性、水氣二相流結(jié)構(gòu)、消能效果、摩擦因子等方面。對(duì)于這種新興溢洪道設(shè)置后的水流流態(tài)、壓力分布、三維流場(chǎng)等均未見(jiàn)報(bào)道;尾坎的位置也只能局限在臺(tái)階水平面中部與末端之間,施工及養(yǎng)護(hù)維修均較為不便。文章依托具體工程,應(yīng)用三維數(shù)值模擬技術(shù)對(duì)水電站溢洪道尾坎位置影響溢洪道內(nèi)水流流態(tài)、壓強(qiáng)、流場(chǎng)、紊動(dòng)能耗散率等的程度展開(kāi)比較分析,為該水電站溢洪道消能率的提升及階梯抗空蝕安全性的提高提供借鑒參考。
某大型水電站溢洪道泄槽寬10.0m,工作水頭最大為140m,單寬流量最大達(dá)到75m3/(s·m),溢洪道進(jìn)口處設(shè)置WES實(shí)用堰。下游河道較為狹窄,且巨大滑坡體存在于溢洪道出口河段左岸,在挑流水舌的影響下穩(wěn)定性很難保證。在溢洪道出口無(wú)法展開(kāi)大規(guī)模調(diào)整的情況下,階梯式溢洪道設(shè)置主要面臨溢洪道單寬流量大、工作水頭高等問(wèn)題。結(jié)合該工程實(shí)際情況及類似工程設(shè)計(jì)經(jīng)驗(yàn),在前7級(jí)階梯水深尚未達(dá)到均勻流狀態(tài),且溢洪道單寬流量超76m2/s時(shí),立面負(fù)壓取-24.5kPa,再加上前幾級(jí)階梯底部存在較低摻氣濃度,空蝕破壞十分容易發(fā)生。為此,必須在恰當(dāng)位置增設(shè)摻氣減蝕裝置,以提升消能率及階梯安全性。
應(yīng)用湍流模型中k-ε模型展開(kāi)坎式階梯溢洪道復(fù)雜水流流態(tài)模擬分析效果較好,而建立在歐拉網(wǎng)格下的Volume of Fluid界面追蹤方法能較為準(zhǔn)確地追蹤自由面,為此,文章采用k-ε模型與Volume of Fluid界面追蹤方法相結(jié)合的三維數(shù)值模擬技術(shù),對(duì)不同體型階梯溢洪道水力特性展開(kāi)模擬分析[1]。
文章采用控制體積法展開(kāi)三維數(shù)值模擬方程求解,水流壓力和流速的耦合均采用PIOS算法。為確保來(lái)流平穩(wěn),在WES實(shí)用堰入口1.5m處設(shè)置庫(kù)區(qū)入口,并以此為速度進(jìn)口。上邊界采用與大氣相通的壓力進(jìn)口;在與最后一級(jí)臺(tái)階相距1.5m處設(shè)置下游出口,以避免水躍的不利擾動(dòng);邊墻和底板均采用無(wú)滑移固壁邊界條件。
結(jié)合該水電站坎式階梯溢洪道設(shè)計(jì)思路,擬采取的階梯溢洪道結(jié)構(gòu)見(jiàn)圖1。其階梯溢洪道主要包括上游庫(kù)區(qū)、WES實(shí)用堰、階梯段、出水渠等部分。其中WES實(shí)用堰長(zhǎng)Lcrest和寬W分別為1.01m和0.52m;上游進(jìn)口圓角半徑r為0.08m;溢洪道共設(shè)置10級(jí)臺(tái)階,且各級(jí)臺(tái)階尺寸完全一致,長(zhǎng)度l均為0.2m,寬h全部為0.1m;臺(tái)階水平面尾坎高w為0.031m,厚度lw為0.015m;尾坎和臺(tái)階水平面末端凸角間距用d表示。
為展開(kāi)尾坎具體設(shè)置位置的研究,文章擬定出4種階梯溢洪道體型,分別表示為體型1~4,其尾坎與臺(tái)階水平面末端凸角之間的距離依次為d=3l/20、d=6l/20、d=9l/20、d=12l/20。
圖1 坎式階梯溢洪道結(jié)構(gòu)
坎式階梯溢洪道中水流流態(tài)根據(jù)來(lái)流條件及臺(tái)階幾何參數(shù)的不同可以劃分成跌落水流、過(guò)渡水流和滑移水流[2]。文章主要研究最為常見(jiàn)的滑移水流流態(tài),并假定其滿布于各級(jí)臺(tái)階,無(wú)空腔。通過(guò)對(duì)不同體型階梯溢洪道中軸面水面線分布情況(圖2)的分析,溢洪道臺(tái)階水流流態(tài)穩(wěn)定,且水面線基本平行于虛擬底部;在尾坎頂托的作用下,體型1水面線比其余體型水面線略高;但整體而言,尾坎移動(dòng)對(duì)水面線的影響并不明顯。
通過(guò)對(duì)不同體型溢洪道臺(tái)階上流線分布情況的分析,一定底部上方的流線呈平行態(tài)勢(shì);而底部以下區(qū)域尾坎上下游均存在持續(xù)旋轉(zhuǎn)的橫軸旋渦,此類旋渦是階梯溢洪道水流能量削減的主要形式。臺(tái)階上方主流和渦流持續(xù)發(fā)生能量交換,在兩種形式的摻雜作用下,水路動(dòng)能得到有效削減。隨著尾坎持續(xù)向上游的移動(dòng),其上游臺(tái)階立面間所預(yù)留的空間越來(lái)越小,渦流受到持續(xù)擠壓后尺度減小。尤其在體型4下,所預(yù)留的空間最為狹小,渦流只能形成于尾坎上方,底部幾乎不存在渦流,整個(gè)過(guò)程中均難以形成大尺度旋渦。
圖2 不同體型階梯溢洪道中軸面水面線分布
通過(guò)對(duì)臺(tái)階水平面凸角流速分布情況的比較分析看出,對(duì)于不同體型的坎式階梯溢洪道而言,越是臨近虛擬底部區(qū)域,臺(tái)階水平面末端凸角上流速分布差別越大;且隨著尾坎向上游的持續(xù)移動(dòng),流速也隨之增大;但自由面周?chē)魉俜植驾^為接近。造成以上流速分布趨勢(shì)的原因在于,尾坎與臺(tái)階凸角距離較近時(shí),尾坎所形成的虛擬底部高出臺(tái)階凸角所形成的虛擬底部,故尾坎底部下方流速小;隨著尾坎的移動(dòng),其所形成的虛擬底部持續(xù)下降,滑移主流也逐漸靠近臺(tái)階凸角,底部流速持續(xù)增大;當(dāng)尾坎虛擬底部低出臺(tái)階凸角底部,滑移主流直接接觸凸角,尾坎所在位置對(duì)流速的影響已微乎其微。以上情況見(jiàn)圖3a。
在體型不變情況下,不同臺(tái)階凸角流速分布規(guī)律基本一致,見(jiàn)圖3b,圖中y為從臺(tái)階凸角開(kāi)始起算的與虛擬底部垂直的水深,ymax為相應(yīng)水深最大值。沿著水深的變化,流速呈先增后減的趨勢(shì),且最大和最小流速依次出現(xiàn)在y/ymax=0.78及y/ymax=0處。
(a)不同體型 (b)不同臺(tái)階
將三維數(shù)值模擬的流速結(jié)果和幾何模型試驗(yàn)結(jié)果展開(kāi)比較,階梯溢洪道內(nèi)水流流速數(shù)值模擬結(jié)果和試驗(yàn)結(jié)果吻合度高,不同體型及不同臺(tái)階下平均相對(duì)誤差最大值僅為12%和8%,說(shuō)明文章分析結(jié)果穩(wěn)定可靠[3]。
根據(jù)分析,可將臺(tái)階水平面壓強(qiáng)分布劃分成尾坎上游和尾坎下游兩個(gè)區(qū)域。①上游區(qū)流線分布顯示出該區(qū)域內(nèi)橫軸旋渦較大,且在旋渦、尾坎、上游臺(tái)階立面的綜合作用下使區(qū)域內(nèi)凹角壓強(qiáng)增大;隨著尾坎逐漸向上游的移動(dòng),溢洪道體型1~3水流壓強(qiáng)持續(xù)減小,而體型4的壓強(qiáng)呈直線形分布,原因在于該體型下的尾坎和上游臺(tái)階立面間空間較小,不存在明顯旋渦,水流也幾乎呈靜止?fàn)顟B(tài)。②下游區(qū)域臺(tái)階水平面受到滑移主流直接沖擊的位置處壓強(qiáng)最大;體型1和體型2下尾坎形成的虛擬底部比臺(tái)階凸角形成的虛擬底部高,故相應(yīng)體型下滑移主流無(wú)法沖擊到臺(tái)階水平面,臨近主流區(qū)的水流壓強(qiáng)自然較大;而體型3和體型4種尾坎所形成的虛擬底部比臺(tái)階凸角虛擬底部低,尾坎移動(dòng)并不影響主流沖擊臺(tái)階水平面的位置,故相應(yīng)體型下壓強(qiáng)最值及沿程分布十分接近。綜上所述,體型1和體型2的臺(tái)階立面壓強(qiáng)分布規(guī)律相似,即隨著水深的增大,壓強(qiáng)基本呈降低趨勢(shì),負(fù)壓區(qū)基本不存在;體型3和體型4臺(tái)階立面壓強(qiáng)分布規(guī)律也較為相似,即隨著水深的增大,壓強(qiáng)呈先降后增態(tài)勢(shì),且在y/h=0.9的區(qū)域內(nèi)達(dá)到峰值。
臺(tái)階水平面末端尾坎起到一定的凸體作用,此處水流很容易脫離邊壁,降低該區(qū)域內(nèi)局部壓強(qiáng),甚至引發(fā)空化危險(xiǎn)。為此文章以壁面壓強(qiáng)較低的下游立面為對(duì)象展開(kāi)其壓強(qiáng)分布趨勢(shì)規(guī)律的探析。根據(jù)尾坎下游立面壓強(qiáng)分布趨勢(shì)圖,體型1和體型2尾坎下游立面壓強(qiáng)在凸體范圍內(nèi)變化劇烈,且存在極小值;而體型3和體型4下的壓強(qiáng)則隨著水深的增大而呈持續(xù)減小趨勢(shì)。造成這種差異的原因在于,體型1和2下尾坎受到滑移主流的直接沖擊,故尾坎下游立面周?chē)鲊?yán)重脫離邊壁,相應(yīng)范圍內(nèi)的壓強(qiáng)驟降,但并未表現(xiàn)出負(fù)壓;隨著尾坎持續(xù)向上游滑動(dòng),臺(tái)階水平面又進(jìn)而遭到滑移主流的直接沖擊,部分水流從下游流走,未流走的水流則回流至尾坎凹角處加入旋渦運(yùn)動(dòng),且水流始終未表現(xiàn)出脫離邊壁的情況。由此,溢洪道體型3和體型4下的水流壓強(qiáng)明顯比體型1和體型2大。
文章以紊動(dòng)能耗散率指標(biāo)體現(xiàn)并衡量紊動(dòng)能向分子熱運(yùn)動(dòng)動(dòng)能轉(zhuǎn)化的速率。該指標(biāo)取值越大,意味著消能效果越好[4]。通過(guò)對(duì)不同體型坎式階梯溢洪道中軸面紊動(dòng)能耗散率分布云圖的分析看出,對(duì)于相同體型而言,上游紊動(dòng)能耗散率較小,越往下游越大;體型1和體型2中相同臺(tái)階最大紊動(dòng)能耗散率均位于尾坎頂;而體型3和體型4最大紊動(dòng)能耗散率則出現(xiàn)在臺(tái)階凸角上游。造成這種現(xiàn)象的原因在于臺(tái)階凸角周?chē)部彩艿交浦髁鞯闹苯記_擊,而引起紊動(dòng)能的較大耗散。
該坎式階梯溢洪道中軸面不同臺(tái)階及尾坎最大紊動(dòng)能耗散率模擬結(jié)果見(jiàn)表1,由表可知,臺(tái)階面和尾坎的最大紊動(dòng)能耗散率均隨水流方向呈增大趨勢(shì);體型1和體型2因尾坎與臺(tái)階凸角距離較近,尾坎紊動(dòng)能耗散率最大值明顯高于臺(tái)階面;而對(duì)于體型3和體型4尾坎與臺(tái)階凸角距離較遠(yuǎn)的情況,臺(tái)階面紊動(dòng)能耗散率最大值明顯較高[5-6]。
表1 最大紊動(dòng)能耗散率
溢洪道消能效果及水流能量損失程度主要通過(guò)消能率指標(biāo)予以體現(xiàn),消能率是溢洪道上下游斷面水頭差和上游斷面總水頭之比。根據(jù)計(jì)算,體型1~4階梯溢洪道消能率依次為48.8%、47.1%、44.0%和44.1%。從結(jié)果看出,體型1和體型2溢洪道的消能率較為接近,而體型3和體型4溢洪道的消能率無(wú)明顯差別。結(jié)合以上對(duì)最大紊動(dòng)能耗散率的分析,在尾坎和臺(tái)階凸角距離較小時(shí),發(fā)揮主要消能作用的結(jié)構(gòu)是尾坎,反之則主要由臺(tái)階面發(fā)揮消能作用。整體而言,雖然尾坎和臺(tái)階凸角距離較小時(shí)消能率會(huì)增大,但尾坎移動(dòng)對(duì)消能率的作用微乎其微。
文章所得出的坎式溢洪道內(nèi)水力特性三維數(shù)值模擬結(jié)果和幾何模型試驗(yàn)結(jié)果吻合度較高,表明應(yīng)用湍流模型中的k-ε模型展開(kāi)坎式階梯溢洪道復(fù)雜水流流態(tài)模擬分析效果較好。文章研究還表明,尾坎設(shè)置位置對(duì)階梯溢洪道水面線影響較小,且隨著尾坎位置的上移,其上游側(cè)旋渦尺度不斷減小,旋渦中心也持續(xù)移動(dòng)至上游,但尾坎下游側(cè)旋渦的尺寸所受影響并不大。對(duì)于不同體型的階梯溢洪道結(jié)構(gòu)而言,隨著尾坎設(shè)置位置向上游的移動(dòng),同一臺(tái)階凸角流速呈增大趨勢(shì);對(duì)于同一體型而言,下游臺(tái)階凸角流速明顯比上游快,但越往底部這種變化趨勢(shì)越弱化。尾坎設(shè)置位置越向上游,臺(tái)階水平面最大壓強(qiáng)、最大紊動(dòng)能耗散率等均表現(xiàn)出不同的變化趨勢(shì)及規(guī)律性,均可作為該水電站溢洪道工程尾坎設(shè)計(jì)與設(shè)置的可靠依據(jù)。