周 迎，鄧 歡，黃永艷

（1.山東省水利勘測設(shè)計(jì)院,山東 濟(jì)南 250013；2.陜西瀚泰水利水電勘測設(shè)計(jì)有限公司,陜西 西安 710065）

1 前言

隨著我國的工農(nóng)業(yè)的迅速發(fā)展,當(dāng)前很多灌區(qū)通過修建閘站合建樞紐進(jìn)行更新改造以滿足灌溉排水的需求。由于閘站合建樞紐設(shè)置在渠道的同一斷面上,水流通過樞紐時(shí),水閘和泵站各自單獨(dú)運(yùn)行,過水?dāng)嗝嫱豢s和突擴(kuò),上下游水流很容易出現(xiàn)偏流、回流和橫流等復(fù)雜流態(tài),甚至?xí)霈F(xiàn)惡性漩渦,以上這些不良流態(tài)會(huì)對渠道內(nèi)的水流產(chǎn)生不利影響,造成渠道沖刷和泥沙淤積等問題,對渠道的安全運(yùn)行有很嚴(yán)重的影響[1-2]。為解決流場中存在的問題,需要在閘下渠道中布設(shè)整流措施對流場進(jìn)行優(yōu)化,目前整流的工程措施主要有底坎、導(dǎo)流墩和壓水板等[3]。其中導(dǎo)流墩整流的主要原理是：通過合理的布置,水流通過導(dǎo)流墩時(shí),沿布設(shè)方向擴(kuò)散至整個(gè)流道,可以規(guī)范整個(gè)流道內(nèi)水量的分布、改善水流流態(tài)和避免渠道邊坡發(fā)生脫流[4],國內(nèi)很多學(xué)者通過物理模型試驗(yàn)或數(shù)值模擬對導(dǎo)流墩做了大量研究并取得了成果[5-6]。本文針對某閘站合建樞紐中存在的水流通過閘門后偏居于渠道內(nèi)的閘門側(cè),造成泵站側(cè)形成較大的漩渦等問題,擬選雙導(dǎo)流墩作為整流措施進(jìn)行分流,采用數(shù)值模擬的方法對雙導(dǎo)流墩的傾斜角進(jìn)行優(yōu)選,為同類工程提供參考。

2 計(jì)算模型及方案

2.1 模型建立

選取某閘站合建樞紐為研究對象,節(jié)制閘底板檻頂高程定為-0.50 m。啟閉機(jī)工作橋布置在閘首中部,寬4.5 m,長16 m,平臺(tái)高程13.0 m,啟閉機(jī)室樓梯間布置在節(jié)制閘東岸。泵站為雙向泵站,功能為排澇及從攔路港引水進(jìn)行水資源調(diào)度。泵型為豎井式貫流泵,單泵流量10 m3/s,設(shè)計(jì)流量40 m3/s。水泵安裝高程-1.15 m。泵站垂直水流方向總寬27.6 m,順?biāo)鞣较蜷L29.0 m。采用GAMBIT進(jìn)行建模,以閘門寬度為18 m為例,其模型見圖1。

圖1 閘站合建樞紐及渠道三維計(jì)算模型

2.2 網(wǎng)格剖分

在CFD進(jìn)行求解時(shí),對模型劃分的網(wǎng)格質(zhì)量和計(jì)算是否準(zhǔn)確密切相關(guān)。網(wǎng)格單元數(shù)較多時(shí),計(jì)算結(jié)果的準(zhǔn)確性能得到保證,但是計(jì)算量大,計(jì)算時(shí)間長；網(wǎng)格單元數(shù)較少時(shí),計(jì)算結(jié)果產(chǎn)生較大偏差。此外網(wǎng)格的類型包括結(jié)構(gòu)化網(wǎng)格和非結(jié)構(gòu)化網(wǎng)格,結(jié)構(gòu)化網(wǎng)格計(jì)算精度較高但劃分難度大,而非結(jié)構(gòu)化網(wǎng)格幾何適應(yīng)性較好,網(wǎng)格劃分簡便,但計(jì)算精度較結(jié)構(gòu)化網(wǎng)格有所降低。因此,在計(jì)算時(shí)應(yīng)合理的選擇網(wǎng)格的大小和類型?？紤]閘站合建樞紐構(gòu)造的復(fù)雜性較大,本文借助GAMBIT軟件對閘站合建樞紐模型進(jìn)行非結(jié)構(gòu)化網(wǎng)格劃分,在關(guān)鍵過流區(qū)域進(jìn)行網(wǎng)格加密,經(jīng)網(wǎng)格無關(guān)性分析確定網(wǎng)格單元數(shù)為300萬。 圖2為計(jì)算區(qū)域的網(wǎng)格剖分示意圖。

圖2 計(jì)算區(qū)域網(wǎng)格示意圖

2.3 計(jì)算方案

在閘后渠道中設(shè)置導(dǎo)流墩來調(diào)整閘后的流態(tài)時(shí),考慮到本文為閘站分側(cè)布置,兩個(gè)導(dǎo)流墩布設(shè)方向見圖3。影響導(dǎo)流墩整流效果的因素很多,包括導(dǎo)流墩長度、 導(dǎo)流墩的夾角、相鄰導(dǎo)流墩的間距、導(dǎo)流墩與泵閘分隔墻的垂直水流流向距離、順?biāo)鞣较驅(qū)Я鞫斩疹^與泵閘隔墻的距離等因素。受文章篇幅所限,本文僅針對雙導(dǎo)流墩夾角的優(yōu)選進(jìn)行介紹,將與導(dǎo)流墻距離近的導(dǎo)流墩記為墩A,設(shè)墩A與水平方向的夾角為1,另外一根導(dǎo)流墩記為墩B,設(shè)墩B與水平方向的夾角為2,見圖3。由于墩A起到主要的分流作用,因此先假定墩B的傾斜夾角2不變,對墩A的傾斜角進(jìn)行研究,將墩B的傾斜角2固定為10°,將墩A和水平方向的傾斜夾角1設(shè)為10°, 15°, 20°和25°并對每個(gè)角度下的流場進(jìn)行數(shù)據(jù)計(jì)算并分析。當(dāng)?shù)玫蕉誂與水平方向的最優(yōu)夾角后使其固定不變,將墩B和水平方向的傾斜夾角 2設(shè)為-5°, 0°, 5°, 0°和13°并進(jìn)行流場數(shù)值模擬計(jì)算分析,得出墩B的傾斜角度的最優(yōu)取值范圍。

圖3 導(dǎo)流墩布設(shè)示意圖

3 流態(tài)評價(jià)指標(biāo)

對于閘站合建工程,水閘和泵站一般都是各自運(yùn)行,當(dāng)水閘運(yùn)行進(jìn)行引水排澇工作時(shí),由于水閘布置在渠岸的一側(cè),過水?dāng)嗝婷娣e減小,上游水流在泵站前逐漸收縮并向閘門側(cè)集中,水流被擠壓,原本順直的水流流向出現(xiàn)彎曲；水流出閘后,右側(cè)渠道突擴(kuò),在渠道的右側(cè)很容易形成回流區(qū),而且回流區(qū)的范圍較大,嚴(yán)重的最大寬度甚至能占到渠道寬度的一半,回流區(qū)在泵站側(cè),很容易造成泥沙沉積在泵站前,對泵站的運(yùn)行造成不利的影響。當(dāng)水流的流量、閘門的寬度發(fā)生變化時(shí),閘門對渠道水流的影響也會(huì)發(fā)生變化,引入動(dòng)能修正系數(shù)ak來研究相同水位條件下流量和閘門寬度變化對渠道水流的影響。

ak可反映水流流速在渠道內(nèi)的分布,當(dāng)水流流速在渠道內(nèi)一個(gè)斷面上的分布越不平均,在這個(gè)斷面上的值就越大,代表這個(gè)斷面的水流越紊亂[7]。當(dāng)水流流速在斷面上分布越平均,ak在這個(gè)斷面的值也就越小,趨近于1。水力學(xué)中對于正常穩(wěn)定的渠道,其ak的取值一般為1.05～1.10,因此當(dāng)水流流速計(jì)算的ak值變化不大,且接近于1.10左右時(shí),即可判斷水流流速在斷面上較穩(wěn)定,ak的計(jì)算大小取水文測驗(yàn)給出的經(jīng)驗(yàn)公式計(jì)算[8]：

式中: vm為某一斷面水流的最大速度,v為這一斷面水流的平均速度,ε為斷面的流速不均勻系數(shù)。

4 計(jì)算結(jié)果分析

4.1 墩A的傾角 1優(yōu)選

將墩B傾斜角 2固定為10°,對墩A的傾斜角 1為10°, 15°, 20°和25°的模型分別進(jìn)行數(shù)值計(jì)算,其流速分布見圖4～圖7。從圖4 中可以看出,出閘主流通過墩A和墩B后被分為兩股水流,從泵閘隔墻和墩A之間流道中經(jīng)過的水流對泵站側(cè)形成的回流區(qū)進(jìn)行沖擊,使得回流區(qū)影響范圍縮短；通過左側(cè)翼墻和墩B之間渠道的水流對渠道左側(cè)岸坡形成的脫壁區(qū)進(jìn)行沖擊,使得水流脫壁現(xiàn)象幾乎消失。由于兩根導(dǎo)流墩將水流分成兩股,故水流主要集中在渠道的兩側(cè),因此在閘后80 m～120 m處的流道中間形成了一段小流速區(qū)域,對流場產(chǎn)生不利影響。

將圖5 與圖4 進(jìn)行對比,可知隨著墩A的傾斜角1的增加,泵站側(cè)形成的回流區(qū)明顯變小,回流區(qū)影響范圍從閘后80 m縮短至閘后70 m處,而且渠道左側(cè)岸坡處出現(xiàn)脫壁區(qū)的寬度也有了一定程度的減小,由于過閘水流被分成三部分,通過墩A和墩B中的水流對渠道中間的靜水區(qū)進(jìn)行了沖擊,使得存在于渠道中間的死水區(qū)范圍變小。

圖4 墩A與水平方向的夾角為10°的流速分布圖

由圖6可以看出,當(dāng)流速繼續(xù)增大時(shí),整個(gè)流場中水流的流速分布情況完全發(fā)生了紊亂,水流集中于一側(cè)的現(xiàn)象更加集中,墩A的傾斜角1的繼續(xù)增大并沒有減小泵站側(cè)形成的回流區(qū),反而使得泵站側(cè)回流區(qū)影響長度增大至閘后80 m處,而且造成左側(cè)渠道岸坡水流脫壁范圍的延長。由于墩A的傾斜角 1太大,造成墩后形成大范圍的死水區(qū),其寬度已經(jīng)占到渠道寬度的一半多,水流幾乎全被擠壓至左側(cè)渠道,水流流場直到閘后180 m處才趨于平順,與圖5 對比,水流經(jīng)過墩A和墩B后流速分布不均勻的狀況更加嚴(yán)重。

圖5 墩A與水平方向夾角為15°的流速分布圖

圖6 墩A與水平方向的夾角為20°的流速分布圖

由圖7 可以看出,隨著墩A角度的繼續(xù)增加,泵站側(cè)形成的回流區(qū)寬度逐漸減小,但是對回流區(qū)的長度變化沒有太大影響,渠道左側(cè)岸坡處的水流脫壁現(xiàn)象得到很好的改善；墩A傾斜角的繼續(xù)增加使得通過墩A和墩B中間的水流方向發(fā)生突變,導(dǎo)致墩后的水流流速分布發(fā)生紊亂,形成折沖水流,會(huì)對渠道岸邊形成沖刷,嚴(yán)重不利于渠道的運(yùn)行安全。

圖7 墩A夾角為25°的流速分布圖

為驗(yàn)證墩A在不同角度下特征斷面的動(dòng)能變化,分別選取閘后80 m、130 m、180 m處的斷面動(dòng)能修正系數(shù)ak,并將其繪制成如圖8 所示的折線圖。從3 個(gè)不同斷面的動(dòng)能修正系數(shù)趨勢圖中可以看出,在導(dǎo)流墩的其它影響參數(shù)保持不變的前提下,隨著墩A傾斜角度的增加,ak呈現(xiàn)先減小后增大的趨勢,這是由于隨著墩A角度的增加,導(dǎo)流墩的分流效果發(fā)揮明顯,但是當(dāng)墩A傾斜角過于大時(shí),水流通過墩A后,流向發(fā)生大的突變,以至在渠道中間形成折沖水流,導(dǎo)致ak增加。結(jié)合流速分布圖,當(dāng)墩A傾斜角度為15°時(shí),導(dǎo)流墩的整流效果最優(yōu),即減小了泵站側(cè)形成的回流區(qū)面積,也減小了左側(cè)渠道邊坡的水流脫壁現(xiàn)象,同時(shí)調(diào)整了水流流速均勻性的分布。

圖8 不同角度三個(gè)斷面動(dòng)能修正系數(shù)折線圖

4.2 墩B的傾角 2優(yōu)選

由4.1 研究可知,墩A的最優(yōu)傾角 1為15°,因此固定墩A的傾角1為15°,將墩B和水平方向的傾斜夾角2設(shè)為-5°, 0°, 5°, 0°和13°并對每個(gè)角度下的流場進(jìn)行數(shù)值計(jì)算,其流速分布如圖9～圖13。

由上圖9可以看出墩B向左傾斜雖然減小了渠道左側(cè)岸坡處水流脫壁區(qū),但是對閘下水流流態(tài)惡化,水流通過兩根導(dǎo)流墩的中間區(qū)域后,在墩A后的渠道中形成了大范圍的低流速區(qū)域,水流主要被擠壓至左側(cè)渠道中,同時(shí)存在于中間渠道的低速區(qū)造成水流的流向發(fā)生偏折,嚴(yán)重影響了渠道中水流流速分布的均勻性,而且如果水流的含沙率高,死水區(qū)的存在會(huì)造成泥沙沉積于泵站側(cè)的渠道中,影響泵站工作時(shí)的效率。

由圖10 與圖9 進(jìn)行對比分析,我們可以看出,通過兩根導(dǎo)流墩中間的水流流出后對墩A后形成的小流速區(qū)進(jìn)行沖擊,使得低流速區(qū)的范圍減少,由閘后200 m縮短至閘后120 m左右,泵站側(cè)形成的回流區(qū)也有了明顯的減小,圖10中的水流流態(tài)分布與圖9 相比更加均勻。從圖9 與圖10 兩個(gè)流速分布圖的對比得出,布置墩B時(shí)應(yīng)向泵站側(cè)傾斜,才能達(dá)到良好的整流效果,否則不僅不會(huì)達(dá)到預(yù)期的效果,還會(huì)造成渠道內(nèi)水流流態(tài)的紊亂。

圖9 墩B傾斜角度為-5°的流速分布圖

圖10 墩B傾斜角度為0°的流速分布圖

由圖11 可以看出,墩B角度為5°時(shí),流場內(nèi)流速的分布比較均勻,雖然在通過墩A后形成了一個(gè)小面積的小流速區(qū)域,但是相比角度為-5°和0°的狀況,小流速區(qū)域已經(jīng)大大減小,而且通過小流速區(qū)域的水流流向平順,未發(fā)生大角度的偏折。對比以上三種不同角度的流速分布圖,發(fā)現(xiàn)隨著墩B傾斜角度的增加,導(dǎo)流墩后渠道的流態(tài)得到了一定程度的改善,但是泵站側(cè)形成的回流區(qū)沒有明顯變化,由此可以得出,對泵站側(cè)形成的回流區(qū)起到主要影響的是墩A,墩B的作用主要是改善渠道水流的流態(tài)。

將圖12 與圖11 對比發(fā)現(xiàn),渠道中流態(tài)的變化并不明顯,但隨著墩B的傾斜角度增加至10°時(shí),通過墩B與翼墻之間的水流對墩A后形成的小流速區(qū)進(jìn)行的沖擊,使得墩A后的小流速區(qū)明顯減小。由于墩B傾斜角度的增加,在墩B后的小流速區(qū)域范圍增加,通過小流速區(qū)域的水流流向發(fā)生大的偏折。對比角度為5°和10°的兩個(gè)流速分布圖,水流流態(tài)并未發(fā)生大的改變,只是在墩B的角度增加時(shí),墩B后小流速區(qū)范圍會(huì)增大。

圖11 墩B傾斜角度為5°的流速分布圖

圖12 墩B傾斜角度為10°的流速分布圖

從圖13 中可以看出,渠道中的流態(tài)由于墩B傾斜過大發(fā)生了流態(tài)紊亂現(xiàn)象,水流通過墩B時(shí),墩B不僅沒有起到導(dǎo)流的作用還阻礙了水流的正常流動(dòng),以至于在墩B后形成了大范圍的低流速區(qū)域,同時(shí)低流速區(qū)域中出現(xiàn)了一定程度的回流現(xiàn)象。水流被分流成兩股水流,分別存在于左側(cè)渠道和右側(cè)渠道。隨著墩B的角度增大至13°,泵站側(cè)回流區(qū)的影響范圍不僅沒有減小反而增加了,由以上分析可以看出,墩B的傾斜角度不能太大,否則將造成水流主要分布于渠道的兩側(cè),而中間渠道的水流流速很小的不良流場分布。

圖13 墩B傾斜角度為13°的流速分布圖

同樣采用動(dòng)能修正系數(shù)的值對上面用流速分布圖得到的結(jié)論進(jìn)行驗(yàn)證,分別選取閘后80 m、130 m和180 m處的斷面對墩B不同傾斜角 2下的流場計(jì)算值,將三個(gè)不同斷面的ak值繪成如圖14所示的折線圖。從圖中可以看出,在墩A和墩B其它影響流態(tài)的因素保持不變的前提下,隨著墩B向泵站側(cè)傾斜角度的增加,ak呈現(xiàn)先減小后增大的趨勢,這是由于隨著墩B角度的增加,通過墩B的水流對墩A后形成的小流速區(qū)域進(jìn)行了沖擊,使得墩A后的小流速區(qū)域范圍減小,但是當(dāng)墩B角度過于傾斜時(shí),不僅達(dá)不到導(dǎo)流的作用,還會(huì)阻隔水流的正常流動(dòng),因此水流通過墩B后,水流的流向發(fā)生一定的彎折,而且會(huì)在墩B后的渠道中形成大范圍的小流速區(qū)域,導(dǎo)致ak增加,影響渠道中的水流流態(tài),不利于渠道的運(yùn)行安全。對墩B進(jìn)行布置時(shí),墩B應(yīng)向泵站側(cè)傾斜,否則不僅達(dá)不到整流的效果,還會(huì)造成渠道內(nèi)水流的紊亂。墩B傾斜角度可設(shè)置在0°和10°之間,通過對墩B角度為0°、5°、10°進(jìn)行比較,發(fā)現(xiàn)角度為5°時(shí)導(dǎo)流墩的整流效果較好,即減小了泵站側(cè)形成的回流區(qū)的范圍,也調(diào)整了水流流速均勻性的分布,降低了動(dòng)能修正系數(shù)ak的值。

圖14 墩B不同角度下三個(gè)斷面的ak折線圖

5 結(jié)論

基于N-S方程和湍流模型對布置有閘站合建樞紐的渠道建立了三維數(shù)學(xué)模型,通過CFD數(shù)值模擬軟件對整流措施雙導(dǎo)流墩在不同傾斜角下流場進(jìn)行模擬,通過分析速度場分布和渠道特征斷面上的動(dòng)能變化優(yōu)選雙導(dǎo)流墩的傾斜角。研究表明通過調(diào)整雙導(dǎo)流墩的傾斜角,可以達(dá)到改善流態(tài)、消減回流區(qū)的目的。其中,與導(dǎo)流墻較近的導(dǎo)流墩的傾斜角取值范圍應(yīng)選擇在10°～20°之間,15°時(shí)達(dá)到最優(yōu)；與導(dǎo)流墻較遠(yuǎn)的導(dǎo)流墩的傾斜角取值范圍應(yīng)選擇在5°～10°之間,5°時(shí)達(dá)到最優(yōu),否則會(huì)造成渠道中水流紊亂。