燕尾挑坎溢洪道三維數(shù)值模擬研究　　夏鵬飛

三孔出流，此時(shí)剛過(guò)駝峰堰初入收縮段的水流受閘墩（常規(guī)半圓形式）和邊墩等因素影響，泄流在墩尾斜沖相遇，激浪高致使流態(tài)不佳，直觀的可看到水相體積分?jǐn)?shù)分布欠均勻，水面凹凸（）為收縮段中部ｘ＝８不夠平順；圖３ｂ０ｍ處水相體積分?jǐn)?shù)分布，水面中部有兩處突起，由此推斷收縮段水體此時(shí)已形成以（）為約收縮泄槽中心線為對(duì)稱(chēng)軸具有一定擾動(dòng)的沖擊波；圖３ｃ段末摻氣槽前ｘ＝１水流在收縮段末０３．７４ｍ處水的體積分?jǐn)?shù)，水面壅高、水的體積分?jǐn)?shù)分布略欠均勻，沖擊波的影響雖未盡去但是整體已趨于平順。綜合分析可知，三處圖形展示的流態(tài)特點(diǎn)鮮明，雖然模擬的自由水面與實(shí)際水面情況略有差異，但水流運(yùn)動(dòng)的總體形態(tài)和實(shí)際情況是相符的

。

３．２　收縮段水力特性

（）圖３為ｘ＝２因上游為ａ５．３１ｍ處收縮段水的體積分?jǐn)?shù)，

圖３　收縮段面水流形態(tài)

Ｆｉ．３Ｗａｔｅｒｆｏｒｍｏｆｔｈｅｎｅｃｋｅｄｄｏｗｎｓｅｃｔｉｏｎ　　　　　?　ｇ

３．３　摻氣空腔水力特性

（）圖４為縱斷面Ｙ＝－２即左邊壁）處，計(jì)算穩(wěn)定后ａ２．５ｍ（的摻氣空腔的流態(tài)，圖４中，兩條實(shí)線之間的為水體部分，下實(shí)線與溢洪道底之間的部分為摻氣空腔，從圖４中可看出水流再

次落入泄槽內(nèi)空腔前端亦無(wú)明顯回水。上實(shí)線與溢洪道之間為空氣，從圖４中亦可知，空氣的流速分布自水面向上方逐漸（）、（）遞減。由圖４圖４對(duì)比可見(jiàn)，模型試驗(yàn)與數(shù)值模擬的結(jié)ａｂ果吻合較好，均形成了空腔并且摻氣空腔穩(wěn)定完整，令人滿意

。

圖４　計(jì)算值與實(shí)驗(yàn)摻氣空腔對(duì)比

ＦｈｅｃｏｍｕｔｉｎａｎｄｅｘｅｒｉｍｅｎｔａｌｒｅｓｕｌｔｓｃｏｍａｒｉｓｏｎｏｎＮａｅｒｏｅｒｔｉ．４Ｔ　　　　　　　�。穑纾穑穑穑穑穑穑�ｇ　

燕尾挑坎溢洪道三維數(shù)值模擬研究　　夏鵬飛

（）模型試驗(yàn)中實(shí)測(cè)摻模擬摻氣空腔處的長(zhǎng)度為１９ｍ，ａ　　圖４

氣空腔長(zhǎng)度為２二者非常接近，吻合良好。模型試驗(yàn)的實(shí)際２ｍ，摻氣量比模擬的摻氣量大和摻氣空腔邊緣不容易確定等因素可（）還中展示了在摻氣空腔內(nèi)的能是產(chǎn)生誤差的主要原因。圖４ａ氣流速度矢量。圖４中顯示的摻氣空處的氣流漩渦清晰可辨，漩渦的方向均為逆時(shí)針?lè)较�，隨著遠(yuǎn)離摻氣槽位置的主漩渦后伴生的小漩渦有所減小，空氣漩渦的形成有助于增加氣流速度，這樣利于摻氣空腔的形成和向水流中混摻更多的空氣。

２７１

效果，解決了水舌集中與沖刷右岸的問(wèn)題。數(shù)值模擬結(jié)果采用水的體積分?jǐn)?shù)為０．中軸線３確定自由水面。水流流經(jīng)挑坎時(shí)，上的水流最先從挑坎缺口處射出，該股水流在重力作用下作拋體運(yùn)動(dòng)，高度逐漸降低，而兩側(cè)的水流由于圓弧底板的頂托，其）水面高度漸漸高于中軸線水面，這一點(diǎn)從圖５（中可以看出。ｂ當(dāng)中軸線附近的水流脫離底板，底部壓力突然降低為大氣壓，而在靠近邊墻一側(cè)底板上的壓力仍保持原來(lái)的動(dòng)水壓力，因此在邊墻與缺口之間會(huì)出現(xiàn)壓差，在該因素產(chǎn)生的壓差作用下，兩側(cè)圓弧底板上的水流會(huì)向中軸線運(yùn)動(dòng)，形成匯聚趨勢(shì)。通過(guò)觀察圖５可以看起挑之后呈一系列沿縱向不同程度拉伸的最最終兩翼水舌相匯，如圖５小正周期為２π的余弦函數(shù)圖像，

（）所示。數(shù)值模擬水舌形態(tài)與模型試驗(yàn)水舌形態(tài)吻合良好

。ａ

３．４　燕尾挑坎水舌特性

本工程下游河道較窄，屬于狹長(zhǎng)形河道，且右岸坡有一凸出平臺(tái)，因此主要考慮對(duì)挑坎進(jìn)行優(yōu)化使其形成窄長(zhǎng)型水舌。通過(guò)模型試驗(yàn)可以看出，燕尾型挑坎滿足了使水舌縱向擴(kuò)散的

圖５　水舌特性ａｅｒｏｅｒｔＦｉ．５Ｎ　�。穑穑穑穑�ｇ

　　數(shù)值模擬之于模型觀察的優(yōu)勢(shì)之一便是可以通過(guò)截取不同斷面提取下泄水流沿程水力特性的分布規(guī)律。本文中沿程截取了６個(gè)特征水舌斷面，對(duì)其水體形態(tài)進(jìn)行了分析。由文獻(xiàn)［］的研究成果可知水舌沿程變化可分為４個(gè)階段：初始段、過(guò)２

（）渡段、分裂段和破碎段。圖５中數(shù)值提取的水舌形態(tài)同樣反ｂ（）映了這一現(xiàn)象。初始段為圖５中的Ｘ＝１水體初ｂ９５ｍ斷面，入挑坎，水舌僅表層受空氣擾動(dòng)；過(guò)渡段為Ｘ＝２０５ｍ斷面至水舌外觀上是連續(xù)的緊密的，但表面有水顆粒１９ｍ斷面，Ｘ＝２

溢出且可發(fā)現(xiàn)空隙；分裂段為Ｘ＝２３４ｍ斷面至Ｘ＝２５５ｍ斷面，此階段水舌有明顯空隙存在；破碎段為Ｘ＝２水６５ｍ斷面，舌破碎成大水片，斷面輪廓模糊后續(xù)發(fā)展則水舌將完全破碎為小水片或水滴。

本文亦對(duì)燕尾型挑坎水舌挑距的數(shù)值計(jì)算和模型試驗(yàn)結(jié)果做了對(duì)比。采用該工程１∶５０比尺的整體水工模型試驗(yàn)數(shù)據(jù)作為依據(jù)。ＶＯＦ法是在假定水體不破碎和水汽間無(wú)混雜的

條件下進(jìn)行模擬的，因此對(duì)于實(shí)際情況下水舌表面有水顆粒溢出，空氣摻入的模擬有其局限性，圖６中水舌表面的模擬結(jié)果與實(shí)際情況略有偏差但與整體的運(yùn)動(dòng)形態(tài)和模型試驗(yàn)結(jié)果是一致的。聯(lián)系燕尾水舌的特殊形態(tài)，前后拉伸成較薄水簾撲入

圖６　水舌側(cè)視圖

ａｔｅｒａｌｖｉｅｗｏｆｔｈｅｎａｅｆｏｒｍＦｉ．６Ｌ　　　　　�。穑穑�

２８１

水中，故水舌挑距考慮其最近及最遠(yuǎn)挑距更為合理。通過(guò)圖６可得數(shù)值模擬水舌挑距。表１比較了水舌挑距，數(shù)值計(jì)算水舌挑距與模型試驗(yàn)實(shí)測(cè)挑距基本吻合。近、遠(yuǎn)距的計(jì)算值由于ＶＯＦ法本身的局限性及網(wǎng)格不夠密等原因很難精確確定水舌入水點(diǎn)位置，誤差不可避免。

表１　計(jì)算值與試驗(yàn)值水舌挑距對(duì)比Ｔａｂ．１Ｔｈｅｃｏｍｕｔｉｎａｎｄｅｘｅｒｉｍｅｎｔａｌｒｅｓｕｌｔｓ　　　�。穑纾稹。悖铮恚幔颍椋螅铮睿瑁铮颍椋�ｏｎｔａｌｌｅｎｔｈｏｆｗａｔｅｒｔｏｎｕｅ　　　　�。穑纾缥恢锰羯渌�舌近距／ｍ　挑射水舌遠(yuǎn)距／ｍ　

數(shù)值計(jì)算挑距

１９．７�。罚罚�０　

模型實(shí)測(cè)挑距

１９．０　７４．０　

相對(duì)誤差／％

３．７４．０

燕尾挑坎溢洪道三維數(shù)值模擬研究　　夏鵬飛

參考文獻(xiàn)：

［］成都：四川大學(xué)，１１２．Ｄ］．２０　錢(qián)小燕．燕尾型挑坎水舌特性研究［［］張文周．空中水舌運(yùn)動(dòng)特性研究［水力發(fā)電學(xué)報(bào)，２］Ｊ．　劉宣烈，

（）：１８８，２４．４６－５９

［］］劉　鈞．空中摻氣水舌運(yùn)動(dòng)軌跡及射距［天津大學(xué)學(xué)３Ｊ．　劉宣烈，

（）：報(bào)，８９，２０．２３－３１９

［］許唯臨，雷　軍．高水頭岸邊泄洪洞水力特性的數(shù)值模擬４　鄧　軍，

（）：［］水利學(xué)報(bào)，０５，３６１０１２０９－１２１２．Ｊ．２０　　

［］］楊永泉，王玉蓉．挑流消能水氣二相流數(shù)值模擬［水利５Ｊ．　刁明軍，

（）：學(xué)報(bào)，２０３，９２．７７－８０

［］刁明軍，岳書(shū)波，等．溢洪道出口斜切挑坎挑射水舌三維６　薛宏程，

（）：］數(shù)值模擬［水利學(xué)報(bào)，１３，４４６７０３－７０９．Ｊ．２０

［］陳永燦，李永紅．三維ＶＯ７Ｆ模型及其在溢洪道水流計(jì)算　李　玲，

（）：］中的應(yīng)用［水力發(fā)電學(xué)報(bào)，０７，２６２８３－８７．Ｊ．２０

［］］吳時(shí)強(qiáng)，陳慧玲，等．多孔溢洪道泄流三維數(shù)值模擬［８Ｊ．　沙海飛，

（）：水利水電技術(shù)，２０５，３６１０４２－４６．０

［］許唯臨，鄧　軍．窄縫消能工水力特性的數(shù)值模擬與試驗(yàn)９　陳華勇，

（）：］研究［水利學(xué)報(bào)，１２，４３４４４５－４５１．Ｊ．２０

［］］劉少斌，楊曉池．１０ｌｕｅｎｔ軟件在泄水工程中的應(yīng)用［Ｊ．Ｆ　巨　江，

（）：水力發(fā)電學(xué)報(bào)，０９，２８２１１０－１１４．２０

［］］汪　德．不同墩型溢流壩水動(dòng)力特性研究［河海大學(xué)１１Ｊ．　王志東，

（）：學(xué)報(bào)，２０４，３２４３９１－３９４．０

［］１２ｉｒｔＣ　Ｗ，ＮｉｃｈｏｌｓＢＤ．Ｖｏｌｕｍｅｏｆｆｌｕｉｄ（ＶＯＦ）ｍｅｔｈｏｄｆｏｒｔｈｅ�。取　　　　　　�

［］，ｄｎａｍｉｃｓｏｆｆｒｅｅｂｏｕｎｄａｒＪ．ＪｏｕｒｎａｌｏｆＣｏｍｕｔｅＰｈｓｉｃｓ１９８１　　　　　�。�ｙｐｙ（）：３９２０１－２２５．

［］１３ｈａｒｌｅｓＳｅｚｉａｌｅＧ，ＴｈｏｍａｓＧａｔｓｋｉＢ，ＮｅａｓｓａｎＦｉｔｚｍａｕｒｉｃｅ．Ａｎ　Ｃ　　　　�。�

ａｎａｌｓｉｓｏｆＲＮＧｂａｓｅｔｕｒｂｕｌｅｎｔｍｏｄｅｌｓｆｏｒｈｏｍｏｅｎｅｏｕｓｓｈｅａｒ　　　　　　　　ｙｇ］，（）：ｆｌｏｗ［Ｊ．Ｐｈｓ．Ｆｌｕｉｄｓ．１９９１，３９２２８－２２８１．７　　ｙ

４　結(jié)　語(yǔ)

對(duì)含有多孔多閘墩、應(yīng)用標(biāo)準(zhǔn)的ｋ～Ｆ法，ＶＯε紊流模型、收縮段、摻氣槽及燕尾型挑坎的復(fù)雜陡坡溢洪道的水力特性進(jìn)行了三維數(shù)值模擬研究，數(shù)值模擬所得結(jié)果得到了物理模型試驗(yàn)的驗(yàn)證，計(jì)算值與試驗(yàn)測(cè)量值吻合較好。研究表明：Ｆ①ＶＯ模型對(duì)復(fù)雜體型的溢洪道自由水面能夠進(jìn)行比較理想的追蹤，本文中對(duì)流速紊動(dòng)強(qiáng)烈的收縮段水面特性亦能進(jìn)行定性分析，模擬的摻氣空腔內(nèi)的氣流漩渦清晰可辨，故該種方法可有效指導(dǎo)復(fù)雜體型溢洪道水面線研究和摻氣槽內(nèi)設(shè)置摻氣設(shè)施。②計(jì)算結(jié)果表明，本文的數(shù)值模擬方法能夠較好的模擬泄洪洞出口燕尾型挑坎挑射水舌形狀并獲得水舌挑距等重要水力參數(shù)，可為泄水建筑體型優(yōu)化設(shè)計(jì)提供依據(jù)，精度亦能滿足設(shè)計(jì)要求。

□

�崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳崳�

（上接第１１９頁(yè)）

［］王二平，董志慧，等．矩形斷面明渠流速分布的研究及應(yīng)３　孫東坡，

（）：］，用［０４，２４－１５１．１４Ｊ．水動(dòng)力學(xué)研究與進(jìn)展（Ａ輯）２０［］金　輝，張艷艷，等．矩形明渠流速分布特征及其在流量４　王二平，

（）：］量測(cè)中的應(yīng)用［０８，４８．２５－２Ｊ．灌溉排水學(xué)報(bào)，２０

［］韓金旭，彭世彰，等．Ｕ形渠道流速分布特性分析與試驗(yàn)５　楊士紅，

［］（）：Ｊ．農(nóng)業(yè)機(jī)械學(xué)報(bào)，２１２，４３１１９２－９６．０

［］韓金旭，彭世彰，等．Ｕ形渠道水流流速垂向分布規(guī)律及６　楊士紅，

（）：］模擬［１２，３０３３０９－３１４．Ｊ．排灌機(jī)械工程學(xué)報(bào)，２０

［］］呂德生，史菊蘭．Ｕ形渠道斷面測(cè)流方法［７Ｊ．中國(guó)農(nóng)村　呂宏興，