矩形明渠跌水壓力分布特性數(shù)值模擬研究

喬 丹

（楊凌職業(yè)技術(shù)學(xué)院，陜西 咸陽 712100）

0 引言

我國丘陵塬坡地區(qū)較多，在修建灌溉渠系工程時(shí)，跌水常作為上下游水流銜接的落差建筑物而被應(yīng)用其中。在明渠中設(shè)置跌坎后，水流流經(jīng)跌坎時(shí)形成自由跌水現(xiàn)象，跌坎上游附近的水流受到重力作用的影響而產(chǎn)生明顯的加速流動，從緩流轉(zhuǎn)變?yōu)榧绷?，水面降低，跌坎處流線急劇彎曲，流態(tài)較為復(fù)雜?；诿髑乃髁鲬B(tài)和水力特性，采用能量方程、動量方程等方法對跌水測流問題進(jìn)行了理論分析研究。試驗(yàn)證實(shí)，跌坎上游流速均勻斷面到跌坎處斷面之間，沿程各橫斷面壓力分布不再符合靜水壓力分布規(guī)律。因此在明渠跌水流量理論計(jì)算方法的推求中，其關(guān)鍵問題是如何解決跌坎處壓力分布非靜壓分布規(guī)律的難題，單純依靠試驗(yàn)難以準(zhǔn)確測定壓力分布特性，可借助紊流數(shù)值模擬來獲得跌坎處壓力分布特性，為明渠跌水測流水力計(jì)算提供基礎(chǔ)。

本文采用k-ε雙方程紊流模型，結(jié)合VOF方法對矩形明渠自由跌水的水流流態(tài)進(jìn)行了三維數(shù)值模擬，得到了矩形明渠自由跌水的渠道橫斷面壓力分布特性，并將計(jì)算結(jié)果同模型試驗(yàn)值進(jìn)行了比較，為矩形明渠跌水測流理論計(jì)算奠定基礎(chǔ)。

1 數(shù)值模擬計(jì)算模型及工況

1.1 計(jì)算模型介紹

為了與試驗(yàn)數(shù)據(jù)進(jìn)行對比，對計(jì)算結(jié)果進(jìn)行驗(yàn)證，計(jì)算模型與試驗(yàn)?zāi)Ｐ突疽恢?。選擇的模型依然是矩形渠道，在跌坎處分為兩個(gè)部分，上游段為矩形明渠，選擇的渠道的長度為40 m，明渠的寬度為0.3 m，渠道高度為0.1～1.0 m；跌落水舌不受限制，所以在跌坎處兩側(cè)各外擴(kuò)0.15 m，跌坎下游的渠道的寬度為0.6 m，跌坎高度為0.6 m，下游段長度取為5 m。計(jì)算模型如圖1所示（以渠道高度0.4 m為例）。

圖1 矩形明渠自由跌水三維計(jì)算模型

1.2 模擬工況

選取十個(gè)不同的坡度，每個(gè)坡度上對應(yīng)五個(gè)流量進(jìn)行計(jì)算?？偣步⒘?0個(gè)矩形明渠跌水模型進(jìn)行模擬計(jì)算。具體的模擬計(jì)算工況見表1。

表1 計(jì)算工況表

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

在計(jì)算結(jié)束后，用Teclopt對計(jì)算結(jié)果進(jìn)行處理。根據(jù)計(jì)算結(jié)果，對矩形明渠自由跌水過程中的壓力分布特性進(jìn)行分析。

2.1 渠道橫斷面壓力分布特性

本文計(jì)算了矩形明渠自由跌水的渠道縱、橫斷面壓力分布，并同模型試驗(yàn)進(jìn)行了比較。

計(jì)算所得渠道沿程不同橫斷面的壓力分布等值線圖如圖2所示（S=0.005，Q=7 0L/s）。

圖2 跌坎上游不同橫斷面壓力分布等值線圖（S=0.005，Q=70 L/s）

從圖中可以看出，在距離跌坎較遠(yuǎn)處，渠道內(nèi)的水流未受到跌坎的影響，各個(gè)橫斷面壓力等值線互相平行。在跌坎附近由于受到重力和慣性力的作用，壓力分布發(fā)生變化，橫斷面壓力等值線向兩側(cè)傾斜，同一水深處靠邊墻兩側(cè)的壓力小于中間位置處的壓力。跌坎處距離渠底一定距離處出現(xiàn)壓力最大值。

2.2 沿程不同斷面豎向壓力分布特性

從圖2可知，明渠跌水時(shí)，各橫斷面豎向壓力分布發(fā)生了明顯變化，因此，圖3～圖5給出了更加直觀的各個(gè)斷面壓力豎向分布變化規(guī)律。由于試驗(yàn)中測壓孔布置有限，無法對豎向壓力分布進(jìn)行詳細(xì)研究，因此通過數(shù)值計(jì)算對此進(jìn)行分析。

由圖4和圖5中可知，壓力數(shù)值計(jì)算值與實(shí)測值吻合良好。

由圖3～5可知，跌坎邊緣（x=0）處的壓力分布近似于拋物線型，在距離渠底某一高度處出現(xiàn)最大值，該斷面的平均壓力值小于靜水壓力分布的平均壓力。距離跌坎越遠(yuǎn)，水流所受跌坎的影響越小，在距離跌坎一定距離時(shí)，壓力的垂向分布近似直線分布。在負(fù)坡和平坡中由于不存在均勻流，因此其壓力分布并非靜水壓力分布，只是呈線性分布，其壓力值小于靜水壓力；而正坡渠道中，距跌坎一定距離后產(chǎn)生了均勻流，因此其壓力分布為靜水壓力分布。從圖中可以看出，在距跌坎5 cm的范圍內(nèi)，水壓力分布由跌坎處的拋物線型分布逐漸向直線型分布過渡。

在同一水深位置處，距離跌坎越遠(yuǎn)，壓力越大，并且越向上游，壓力增大的幅度越來越小。在負(fù)坡渠道中，由于渠道中的水深從跌坎向上游沿程一直增加，距渠底同一高度的位置其與水面的距離沿程是不斷增加的，因此，距渠底同一高度位置處的壓強(qiáng)沿程一直增加，如圖3所示。在平坡渠道中，距離跌坎約5 m之后，由于渠道中的水深變化較小，因此渠道中各橫斷面壓力分布變化較小，如圖4所示。對于正坡渠道而言，當(dāng)S=0.0240時(shí)，距跌坎0.22 m之后，橫斷面壓力分布已無明顯變化，斷面x=-0.42 m的壓力分布與斷面x=-8.00 m的壓力分布幾乎完全相同（見圖5）。

圖3 S=-0.0112、Q=80.00 L/s時(shí) 不同斷面壓力豎向分布

圖4 S=0、Q=66.12 L/s時(shí) 不同斷面壓力豎向分布

圖5 S=0.0240、Q=58.91 L/s時(shí)不同斷面壓力豎向分布

3 結(jié)語

（1）通過對不同底坡的矩形明渠在不同流量下的自由跌水水流特性的三維數(shù)值模擬，可以得出其跌坎上游渠道沿程各個(gè)橫斷面壓力分布等值線圖以及不同橫斷面壓力豎向分布特性，所得數(shù)據(jù)同模型試驗(yàn)結(jié)果吻合良好，表明本文所采用的紊流數(shù)學(xué)模型是可靠的，可獲得較為滿意的結(jié)果。

（2）壓力分布結(jié)果顯示，上游渠道各橫斷面壓力等值線互相平行，跌坎附近壓力分布發(fā)生變化，同一水深處兩側(cè)的壓力小于中間位置處的壓力。跌坎邊緣處的壓力分布近似于拋物線型，該斷面的平均壓力值小于靜水壓力分布的平均壓力；在距離跌坎一定距離后，壓力的垂向分布近似直線分布。

（3）本文借助紊流數(shù)值模擬獲得矩形明渠自由跌水壓力分布特性，這為后續(xù)的矩形明渠流量理論計(jì)算研究分析提供依據(jù)。