張?zhí)煜?，俞曉東，賀蔚*，張健，肖洋

(1. 河海大學(xué)水利水電學(xué)院，江蘇 南京 210098； 2. 河海大學(xué)水安全與水科學(xué)協(xié)同創(chuàng)新中心，江蘇 南京 210098)

明渠岸邊橫向取水是被廣泛應(yīng)用的取水方式之一，經(jīng)常被用于火電站、核電站和長(zhǎng)距離引水工程中.然而，取水口附近旋渦的產(chǎn)生是對(duì)運(yùn)行造成危害的主要原因之一，所以研究人員十分重視并有針對(duì)性的進(jìn)行了一些研究.主要從二維紊流、三維T型分叉結(jié)構(gòu)、三維層流進(jìn)行數(shù)值模擬研究[1-4].曹繼文等[5]在三維情況下對(duì)岸邊橫向取水口和附近區(qū)域進(jìn)行比較深入地研究，并應(yīng)用N-S方程對(duì)明渠橫向取水進(jìn)行了相應(yīng)的數(shù)值計(jì)算.同時(shí)楊帆[6]結(jié)合試驗(yàn)和數(shù)值模擬研究了取水角度對(duì)于明渠岸邊側(cè)向取水的影響.李濤等[7]則對(duì)不同形式明渠橫向分流口門(mén)的局部水頭損失進(jìn)行了研究.一些文獻(xiàn)也針對(duì)明渠水流進(jìn)行了深入研究，包括用VOF法表示明渠自由液面[8]，不同截面明渠水流的水力特性及流速分布規(guī)律[9]、進(jìn)水口體型優(yōu)化[10]等，這些對(duì)于明渠橫向取水旋渦的研究均起到一定幫助.

已有研究表明，明渠橫向取水旋渦會(huì)對(duì)取水工程帶來(lái)不利影響，為找到有效地削弱措施，需明確其影響因素，然而針對(duì)明渠岸邊橫向取水銜接處旋渦的研究還較少觸及，且相關(guān)水力學(xué)機(jī)理有待完善，所以對(duì)于明渠橫向取水旋渦影響因素的研究是必要的.文中擬以中國(guó)某橫向取水工程為例，結(jié)合數(shù)值模擬和模型試驗(yàn)，模擬取水過(guò)程中的水力學(xué)特性，分析旋渦產(chǎn)生機(jī)理及影響因素，并最終提出削弱措施.

1 研究對(duì)象及方法

1.1 研究區(qū)域

文中主要采用中國(guó)南方某橫向取水工程的相關(guān)工程參數(shù)，研究區(qū)域如圖1所示.在正常運(yùn)行工況下水位為-0.31 m，主河道平均流量為958.4 m3/s，平均流速約為0.52 m/s.橫向取水口位于主河道右岸，翼墻坡度為1∶2.5，取水口口門(mén)寬度為9.0 m，底部高程為-3.5 m，取水設(shè)計(jì)流量20 m3/s.在該工程中，水動(dòng)力條件對(duì)取水效率、口門(mén)泥沙淤積等都有重要影響.

圖1 研究區(qū)域示意圖

在進(jìn)行模型試驗(yàn)中發(fā)現(xiàn)橫向取水過(guò)程中在銜接處易產(chǎn)生旋渦，對(duì)該現(xiàn)象進(jìn)行數(shù)值模擬，分析其產(chǎn)生機(jī)理及影響因素，并進(jìn)行優(yōu)化設(shè)計(jì)和驗(yàn)證.

1.2 數(shù)學(xué)模型

1.2.1 模型建立

數(shù)學(xué)模型采用Gambit軟件建立并劃分網(wǎng)格，同時(shí)通過(guò)Fluent進(jìn)行計(jì)算.計(jì)算區(qū)域包括主河道上下游進(jìn)出口、固體邊壁、水面、取水口、主河道和分河道組成.其中主河道長(zhǎng)4.0 m，寬3.0 m，高0.2 m，分河道長(zhǎng)1.0 m，河底寬0.3 m，水面寬0.8 m，高0.2 m.上下游進(jìn)出口分別設(shè)在主河道上下游斷面，取水口設(shè)在分河道下游斷面.概化模型見(jiàn)圖2.

圖2 概化模型

模型內(nèi)包含2種材料，初始化設(shè)置時(shí)將模型上半部分設(shè)置為空氣，下半部分設(shè)置為液態(tài)水.上游流場(chǎng)進(jìn)口[11]分為2部分，上部分為空氣進(jìn)口，下部分為液體進(jìn)口.空氣進(jìn)口的邊界條件定為壓力進(jìn)口，液體進(jìn)口的邊界條件定為速度進(jìn)口.取水口的邊界條件定為速度進(jìn)口，并設(shè)置為負(fù)值.下游流場(chǎng)出口同樣分為2部分，空氣出口部分采用壓力出口，液體出口部分采用速度進(jìn)口，將速度設(shè)置為負(fù)值.關(guān)于底板、邊墻等均采用無(wú)滑移壁面邊界條件，而模型頂部設(shè)置為壓力進(jìn)口.

1.2.2 控制方程

文中主要研究的是不可壓縮流體，采用標(biāo)準(zhǔn)k-ε模型，該模型應(yīng)用較為廣泛，其連續(xù)方程、動(dòng)量方程、k方程和ε方程以及出現(xiàn)的經(jīng)驗(yàn)常數(shù)均可參考文獻(xiàn)[12].并運(yùn)用VOF法[13]追蹤自由液面，使用有限體積法進(jìn)行離散，并用SMIPLE算法對(duì)取水口及周?chē)拥肋M(jìn)行三維數(shù)值模擬.為研究取水口附近水力特性的變化情況，對(duì)試驗(yàn)特定平面內(nèi)的水體進(jìn)行速度監(jiān)測(cè).

時(shí)均不可壓縮流體的N-S方程為

?·u=0，

(1)

(2)

式中：u為速度；p為壓力；μ為動(dòng)力黏性系數(shù)；ρ為密度；g為重力加速度.

渦量是表述旋渦運(yùn)動(dòng)重要的物理量之一，渦量計(jì)算公式為

Ω=?×u，

(3)

式中：Ω為渦量矢量；u為速度矢量.

在VOF模型中，水氣兩相分別求解質(zhì)量方程、共享同一個(gè)動(dòng)量方程.當(dāng)出現(xiàn)水氣兩相共同存在于一個(gè)計(jì)算網(wǎng)格單元時(shí)，通過(guò)體積平均法來(lái)進(jìn)行計(jì)算

ρ=αaρa(bǔ)+αwρw，

(4)

式中：下標(biāo)a表示氣體；w表示液體.與此同時(shí)，氣體和液體的體積分?jǐn)?shù)相加等于1.

1.2.3 網(wǎng)格尺寸敏感性分析

模型主要采用混合網(wǎng)格，以四面體網(wǎng)格為主，其他網(wǎng)格為輔.針對(duì)不同數(shù)量的網(wǎng)格模型，對(duì)斷面b-b的平均水壓進(jìn)行對(duì)比分析，并繪制成圖3，觀察變化趨勢(shì).當(dāng)網(wǎng)格總數(shù)大于約350萬(wàn)時(shí)，斷面b-b上平均水壓變化差異較小(小于0.5%).綜合考慮后，最終選擇網(wǎng)格總數(shù)約為353.1萬(wàn)的數(shù)學(xué)模型進(jìn)行計(jì)算，物理模型如圖4所示.

圖3 不同網(wǎng)格劃分下斷面b-b平均水壓

圖4 河道取水口物理模型

1.3 物理模型試驗(yàn)

1.3.1 模型建立

物理模型以某橫向取水工程為實(shí)例進(jìn)行建造，按重力相似準(zhǔn)則設(shè)計(jì)，采用正態(tài)模型，模型幾何比尺為1∶30.模型模擬范圍包括取水口、上下游部分主河道等.物理模型寬度為90 m，上下游長(zhǎng)度各225 m.取水口為鋼筋混凝土結(jié)構(gòu)，糙率一般為0.014，比尺為1.763，有機(jī)玻璃的糙率為0.008符合阻力相似要求，所以采用有機(jī)玻璃制成.河道地形用水泥砂漿制作，采用專(zhuān)門(mén)的拉毛工藝處理，基本滿(mǎn)足阻力相似要求.選用多臺(tái)水泵進(jìn)行流量控制，試驗(yàn)前根據(jù)運(yùn)行工況對(duì)各泵流量進(jìn)行率定.

1.3.2 量測(cè)儀器

物理模型共有3個(gè)進(jìn)出水邊界.取水口利用水泵取水，并連接有閥門(mén)及LZS-P塑料轉(zhuǎn)子流量計(jì)，利用閥門(mén)來(lái)控制取水口流量.采用目測(cè)和水位測(cè)針進(jìn)行水位觀測(cè)，兩者相互驗(yàn)證，保證水位測(cè)量的準(zhǔn)確性.對(duì)于流速，采用多臺(tái)LGY-Ⅲ型多功能智能流速儀進(jìn)行量測(cè)，由計(jì)算機(jī)進(jìn)行數(shù)據(jù)高頻采集及記錄，監(jiān)測(cè)點(diǎn)位于取水口附近以及上下游河道均布置有多個(gè)點(diǎn)位，同時(shí)還會(huì)對(duì)一些特殊的點(diǎn)位使用手持流速儀進(jìn)行測(cè)量.對(duì)于一些流態(tài)方面的觀測(cè)，利用攝像機(jī)記錄流場(chǎng)圖像，對(duì)其進(jìn)行分析研究.

1.4 模擬工況設(shè)置

研究在設(shè)計(jì)運(yùn)行水位以及落潮平均流量下進(jìn)行，主要針對(duì)銜接段圓弧半徑對(duì)主河道與橫向取水口銜接段所產(chǎn)生的旋渦進(jìn)行研究，同時(shí)考慮主河道流量以及取水流量的影響.探究其相關(guān)水力學(xué)機(jī)理，提出了削弱措施.具體試驗(yàn)工況見(jiàn)表1.

表1 試驗(yàn)工況

2 模型驗(yàn)證

為了進(jìn)一步驗(yàn)證數(shù)學(xué)模型的可靠性，以模型試驗(yàn)中2個(gè)特定斷面作為研究對(duì)象，選取基本工況S1作為驗(yàn)證工況，在數(shù)學(xué)模型中選取與物理模型試驗(yàn)中相同位置的監(jiān)測(cè)點(diǎn)，并對(duì)水下4.0 cm(原型為水下1.2 m)進(jìn)行流速監(jiān)測(cè).將兩個(gè)斷面的流速分布與物理模型試驗(yàn)結(jié)果進(jìn)行對(duì)比驗(yàn)證，結(jié)果如圖5所示.

3 討 論

3.1 橫向取水引起旋渦水力特性分析

對(duì)基本工況S1進(jìn)行模型試驗(yàn)，工況S1銜接方式為直角銜接，示意圖及實(shí)物圖見(jiàn)圖6a,b.試驗(yàn)中河道來(lái)流較為平順，當(dāng)水流行至銜接段附近，會(huì)在銜接段形成渦帶，單個(gè)旋渦半徑較大，有一定的深度，渦帶長(zhǎng)度基本到達(dá)取水口中軸線，試驗(yàn)現(xiàn)象見(jiàn)圖6c.運(yùn)用數(shù)值模擬對(duì)該現(xiàn)象進(jìn)行模擬，分析橫向取水引起旋渦的水力特性問(wèn)題.

圖5 數(shù)學(xué)模型驗(yàn)證示意圖及斷面流速對(duì)比圖

圖6 工況S1物理模型及試驗(yàn)現(xiàn)象

根據(jù)圖7a得到：主河道與橫向取水口交界處的上游側(cè)存在高流速區(qū)域，由于取水口的影響，取水口內(nèi)的分河道存在局部低流速區(qū)，使得分河道內(nèi)流速分布復(fù)雜，流態(tài)較差，影響取水效率.根據(jù)圖7b得到：銜接段附近存在一條明顯的渦帶，渦帶長(zhǎng)度約占分河道寬度的54.6%，影響范圍較大.根據(jù)圖7c得到：分河道中存在有較大范圍的高湍動(dòng)能區(qū)，銜接段處有一明顯的高湍動(dòng)能帶，與渦帶位置基本相同.

圖7 工況S1系列云圖

研究發(fā)現(xiàn)，在水流流動(dòng)過(guò)程中，壓力會(huì)逐漸下降補(bǔ)償水流運(yùn)動(dòng)產(chǎn)生的損失，橫向取水口與主河道交界處的上游側(cè)存在有一處低壓區(qū)，當(dāng)水流在流至交界處之前，處于加速減壓狀態(tài)，在交界處流速達(dá)到最大，壓力最小.但當(dāng)水流流過(guò)交界處時(shí)河道突然擴(kuò)散，流速突然降低并在分離點(diǎn)處流速降為0，分離點(diǎn)處無(wú)法再讓動(dòng)能轉(zhuǎn)換為壓能，主流將被迫脫離壁面，發(fā)生邊界層與邊界分離的現(xiàn)象，主流按照原來(lái)的流向繼續(xù)流動(dòng)，下游的水流則會(huì)流到主流所空出的部分，由此形成旋渦.

3.2 銜接段圓弧半徑對(duì)旋渦的影響

經(jīng)過(guò)上述數(shù)值模擬和模型試驗(yàn)研究發(fā)現(xiàn)，分離點(diǎn)的位置會(huì)對(duì)旋渦產(chǎn)生一定的影響.分離點(diǎn)越接近繞過(guò)物體的末端，旋渦就會(huì)越小，存在的長(zhǎng)度也會(huì)有所減短，對(duì)取水口的影響也會(huì)越小.改變銜接段的結(jié)構(gòu)會(huì)對(duì)分離點(diǎn)的位置產(chǎn)生一定影響，采用圓弧過(guò)渡，可以使水流擴(kuò)散趨勢(shì)有所緩解，減緩流速降低的速率，同時(shí)使分離點(diǎn)更靠近銜接段末端減緩發(fā)生邊界層分離現(xiàn)象，改善銜接段旋渦現(xiàn)象.由此，提出多種不同工況，并選取其中部分工況進(jìn)行模型試驗(yàn)，觀察不同銜接段圓弧半徑對(duì)于旋渦的影響.

工況S2,S3,S4,S5的區(qū)別主要是銜接段圓弧半徑長(zhǎng)度不同，從圖8a分析發(fā)現(xiàn)：隨著圓弧長(zhǎng)度不斷增加，銜接段的高流速區(qū)域不斷減小，分河道中局部低流速區(qū)也不斷減小，流速分布愈加均勻.從圖8b分析發(fā)現(xiàn)：由于銜接段圓弧不斷增加，使得水流到達(dá)銜接段時(shí)，擴(kuò)散趨勢(shì)有所減緩，旋渦會(huì)部分削弱，與S1工況相比，渦帶長(zhǎng)度分別減小了4.9%,11.2%,16.3%,21.7%.S5工況中渦帶影響范圍只占到分河道寬度的42.0%，渦帶的影響范圍得到明顯減小.從圖8c分析發(fā)現(xiàn)：原本分河道存在的高湍動(dòng)能區(qū)隨著銜接段圓弧半徑的增加在不停地減小，分河道中湍動(dòng)能最大值也有所減小，流態(tài)有所改善.

圖8 不同工況系列云圖

針對(duì)S2,S3工況進(jìn)行模型試驗(yàn)，示意圖和試驗(yàn)圖見(jiàn)圖9a,b.通過(guò)模型試驗(yàn)觀察到：隨著圓弧半徑的不斷增加，銜接段附近出現(xiàn)的旋渦在直徑和深度方面均有所減小，分離點(diǎn)的位置也向下游有所移動(dòng)，旋渦存在的時(shí)間也有所縮短，渦帶長(zhǎng)度也有所減小.其中，S3工況基本已由S1工況的小旋渦變?yōu)楸砻鏈u紋，深度有較明顯的減小，流態(tài)明顯好轉(zhuǎn)，對(duì)提高取水口工作效率有明顯幫助，試驗(yàn)現(xiàn)象見(jiàn)圖9c.

圖9 不同工況試驗(yàn)?zāi)Ｐ图霸囼?yàn)現(xiàn)象

3.3 主河道流速和取水流速對(duì)旋渦的影響

3.3.1 主河道流速的影響

S6,S7,S8工況主河道的流速如圖10所示.

圖10 不同工況云圖

從圖10a發(fā)現(xiàn)：受主河道流速不斷增加的影響，銜接段處高流速區(qū)域有所增加，但是由于取水流量不變，導(dǎo)致分河道處的低流速區(qū)域也在不斷增加，流態(tài)也越加復(fù)雜.從圖10b發(fā)現(xiàn)：隨著主河道流速的增加，銜接段處的渦帶長(zhǎng)度在不斷的增加，同時(shí)旋渦的強(qiáng)度也有所增加.從圖10c發(fā)現(xiàn)：隨著流速的變大，高湍動(dòng)能區(qū)面積也在不斷增加，側(cè)面驗(yàn)證了流態(tài)也愈加復(fù)雜.主河道流速的增加對(duì)橫向取水工程的流態(tài)較為不利.

隨著主河道流速的增加，渦帶長(zhǎng)度分別增加約3.80%,9.10%,17.83%，旋渦的強(qiáng)度也在不斷增加，銜接段處的流態(tài)趨向復(fù)雜.

3.3.2 取水流速的影響

S9,S10,S11工況的差距則主要在取水流量上，取水流量越大，流速越大，如圖11所示.從圖11a發(fā)現(xiàn)：隨著取水流速的增加，銜接段處高流速區(qū)域有所增加，分河道內(nèi)低流速區(qū)域在不斷減小，但在一定程度上會(huì)影響到主河道流速以及流態(tài).從圖11b發(fā)現(xiàn)：隨著取水流速的增加，銜接段處的渦帶方向會(huì)更偏向于取水口方向，雖然渦帶的長(zhǎng)度和強(qiáng)度都有所增加，但影響范圍占分河道寬度的百分比不斷減小.從圖11c發(fā)現(xiàn)：隨著流速變大，高湍動(dòng)能區(qū)面積不斷增加，同時(shí)高湍動(dòng)能區(qū)更集中在分河道中的上游側(cè)，銜接段處的高湍動(dòng)能帶也有所增加，側(cè)面驗(yàn)證了渦帶的增強(qiáng).

取水流速的增加可能會(huì)在一定程度上對(duì)主河道流態(tài)產(chǎn)生影響，同時(shí)隨著取水流速的增加，銜接段處的旋渦會(huì)更偏向于取水口方向，旋渦的強(qiáng)度和長(zhǎng)度有所增加，渦帶長(zhǎng)度分別增加約5.5%,26.5%,34.1%.但渦帶占分河道寬度的百分比有所減小，分別減小約0.6%,3.6%和6.5%.

圖11 不同工況云圖

4 結(jié) 論

文中對(duì)明渠岸邊橫向取水過(guò)程的旋渦進(jìn)行了研究，分析旋渦的產(chǎn)生機(jī)理、影響因素以及削弱措施，并通過(guò)Fluent對(duì)不同工況進(jìn)行三維數(shù)值模擬，結(jié)果表明：

1) 主河道與橫向取水口銜接段附近旋渦的產(chǎn)生原因是：水流在流動(dòng)到交界處時(shí)河道突然擴(kuò)散，流速突然降低并在分離點(diǎn)處流速降為0，分離點(diǎn)處無(wú)法再讓動(dòng)能轉(zhuǎn)換為壓能，主流將被迫脫離壁面，下游的水流則會(huì)流到主流所空出的部分，由此形成旋渦.

2) 改變銜接段的結(jié)構(gòu)銜接方式會(huì)對(duì)分離點(diǎn)的位置產(chǎn)生影響，采用圓弧過(guò)渡，可以使水流擴(kuò)散趨勢(shì)有所緩解，減緩流速降低的速率，同時(shí)使分離點(diǎn)更靠近銜接段末端減緩發(fā)生邊界層分離現(xiàn)象，改善銜接段旋渦現(xiàn)象.

3) 基本工況S1在銜接段附近存在明顯的渦帶，渦帶長(zhǎng)度約占分河道寬度的54.6%，影響范圍較大.隨著銜接段圓弧不斷增加，水流擴(kuò)散趨勢(shì)有所減緩，與S1工況相比，渦帶長(zhǎng)度分別減小了4.9%,11.2%,16.3%,21.7%，其中S5工況中的影響范圍得到明顯減小，影響范圍只占到分河道寬度的42.0%.

4) 主河道流速和取水流速在一定程度上也會(huì)對(duì)銜接段處的旋渦產(chǎn)生部分影響.隨著主河道流速和取水流速的增加，渦帶長(zhǎng)度明顯增加，旋渦的強(qiáng)度也在不斷增加，銜接段處的流態(tài)更為復(fù)雜，對(duì)橫向取水工程的流態(tài)較為不利.