李 靜,李 丹,姜伯樂

(1.長江水利委員會長江科學(xué)院,湖北 武漢 430010;2.武漢大學(xué)水資源與水電工程國家重點實驗室,湖北 武漢 430072)

目前對帶閘墩溢流壩過壩水流的三維數(shù)值模擬比較少見,這主要是由于三維模擬的計算網(wǎng)格單元數(shù)急劇增加,閘墩的存在加大了網(wǎng)格劃分的難度并且使整個流場變得異常復(fù)雜。因此對過壩水流的研究工作通常較少考慮閘墩對流場的影響,而直接以二維流場為研究對象[1-3]。文獻[4]通過分別計算過壩水流和閘墩繞流2個二維模型來近似模擬三維帶閘墩的溢流壩流場,為解決三維過壩水流問題提供了一條途徑,但是將二維數(shù)值模型應(yīng)用于實際帶閘墩溢流壩過壩水流計算還有一定的差距。主要表現(xiàn)在閘墩的存在對水面線位置、壩面壓力、溢流流量等均產(chǎn)生顯著的影響。文獻[5]利用有限體積法對溢流壩三維湍流場進行數(shù)值模擬,但沒有考慮閘墩的影響。文獻[6-14]對溢流壩及湍流場進行了數(shù)值模擬研究。筆者以銀盤水電站帶閘墩溢流壩為例,以雷諾時均方程為基礎(chǔ),利用k-ε紊流模型模擬湍流流動,采用流體體積(VOF)法確定自由水面線的位置,采用有限體積法進行離散求解,并與物理模型試驗結(jié)果進行對比,比較準確地給出了過壩水流的水面線位置及壩面壓力分布,可為泄水工程的水力設(shè)計提供比較可靠的設(shè)計依據(jù)。

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

1.1 模型控制方程

VOF模型適用于跟蹤捕捉兩相或多相流中互不穿透流體間的界面變化。VOF模型通過引入體積分數(shù)、求解體積分數(shù)方程來捕捉界面[15]。

設(shè)某一控制單元內(nèi)第q相體積分數(shù)為αq(0≤αq≤1),則當 αq=0時,該控制單元內(nèi)無第q相流體;αq=1時,該控制單元內(nèi)充滿第q相流體;0＜αq＜1時,該控制單元內(nèi)存在相間界面。當控制單元中存在n相流體時,各控制單元內(nèi)各相體積分數(shù)之和為1,即

本文所研究問題屬于不可壓縮流體的湍流流動問題,因此模型中采用標準k-ε模型,其基本控制方程如下:

1.2 計算區(qū)域和網(wǎng)格劃分

溢流壩堰頂高程為195m,進口壩面頭部采用方程為x2/72+y2/4.0752=1的1/4橢圓曲線,后接WES曲線,曲線方程為y=0.038x1.85,其后由1∶1.8斜坡與半徑為69.80m的反弧段銜接。消力池底板高程為175m。閘孔寬度為15.5m,閘墩寬度為4.5m。本文研究工況:上游水位為215m,下游水位為199.66m,進口水深為40m,閘門開度為0.4,閘門開啟孔數(shù)為8個,下泄流量為10850m3/s。物理模型取2孔,按照重力相似準則設(shè)計,長度比尺為1∶50。為了與物理模型結(jié)果進行比較,二維和三維模型與物理模型尺度相同。根據(jù)WES曲線的特性,x方向坐標原點取在堰頂處。二維模型計算區(qū)域見圖1。

圖1 二維模型計算區(qū)域

由于閘墩附近的幾何形狀十分復(fù)雜,為了保證計算精度,對包含閘墩的三維區(qū)域進行網(wǎng)格劃分時對整個區(qū)域進行分區(qū),在各分區(qū)上采用不同的網(wǎng)格系統(tǒng)。三維模型計算區(qū)域和網(wǎng)格劃分分別見圖2和圖3。

圖2 三維模型計算區(qū)域

圖3 三維模型閘門局部網(wǎng)格劃分

1.3 邊界條件

氣體入口邊界為壓力邊界條件,入口邊界處壓強為大氣壓強。水的入口邊界為速度入口邊界條件。二維計算和三維計算中流速uin均為1.844m/s。入口紊動能k和耗散率ε的邊界條件分別由式(6)、式(7)得出:

式中:Hin為進口水深。

1.3.1 出口邊界

由于出口壓力分布已知,因此出口邊界采用壓力邊界條件,出口邊界上壓力分布為p=ρg(199.66-z)。過堰水流經(jīng)過消力池消能后流態(tài)仍然非常紊亂。水流過水斷面上流速分布需要經(jīng)過一段足夠長的距離才能逐漸調(diào)整至漸變流流速分布。因此,在數(shù)學(xué)模型建立以后,通過大量的重復(fù)測試,并與物理模型試驗對比,確定了合適的出口邊界位置。

1.3.2 壁面邊界條件

壁面邊界采用無滑移邊界,即在壁面上u=v=w=0。近壁區(qū)的黏性底層采用標準壁函數(shù)處理。

1.4 數(shù)值求解方法

本文采用基于非結(jié)構(gòu)網(wǎng)格的有限體積法對計算區(qū)域進行離散,并采用PISO算法進行求解。

2 物理模型

本文用于對比的試驗數(shù)據(jù)采用物理模型試驗獲得。物理模型采用重力相似原則,比尺為1∶50。物理試驗主要采集了水面線和壩面壓力數(shù)據(jù)。物理模型試驗流態(tài)見圖4。

圖4 物理模型試驗流態(tài)

3 計算結(jié)果比較與分析

3.1 水躍區(qū)流線

圖5 二維與三維模擬消力池內(nèi)流線分布

物理模型試驗結(jié)果表明,該控泄工況下消力池內(nèi)發(fā)生了淹沒水躍。圖5給出了消力池內(nèi)水躍的二維模擬與三維模擬結(jié)果。圖5(a)為二維模擬計算結(jié)果,圖5(b)～(d)分別為三維模擬第 1孔、閘墩和第2孔中軸線剖面上的流線分布。從圖5可以看出,三維計算中由于模擬的對稱性,第1孔與第2孔的流線分布基本相同,二維計算的流線分布中存在2個渦,而三維計算的流線分布中只存在1個渦。在三維模型中計算區(qū)域、網(wǎng)格劃分是完全對稱的。但從圖5(b)、5(d)中可以看出第 1孔與第 2孔中軸線剖面上的流線分布并不完全一致。這一現(xiàn)象看似與計算區(qū)域?qū)ΨQ性矛盾。實際上,在三維模擬中水流與空氣的摻混機理十分復(fù)雜,而且兩者的摻混過程屬于非恒定過程,因此三維模擬無法得出完全對稱的結(jié)果。

3.2 水面線位置

圖6給出了水面線的二維模擬、三維模擬以及與物理模型結(jié)果的對比。由圖6可以看出,二維模擬得出的水躍發(fā)生位置與物理模型結(jié)果有較大偏差,三維模擬得出的水面線與物理模型結(jié)果吻合較好,水躍發(fā)生位置與物理模型結(jié)果比較一致。在水躍發(fā)生點之前,二維與三維模擬得出的水面線與物理模型結(jié)果均十分接近。但是水躍發(fā)生位置2種模擬結(jié)果顯著不同。二維模擬的水躍發(fā)生位置比三維模擬提前許多,出現(xiàn)了較大的偏差。

圖6 二維模擬、三維模擬與物理模型水面線對比

3.3 壩面壓力分布

圖7 二維模擬、三維模擬與物理模型壩面壓力分布對比

圖7給出了壩面壓力分布的二維模擬、三維模擬與物理模型結(jié)果的對比。壩面壓力采集范圍為樁號-11.70～37.81m。從圖7中可以看出,由于閘墩的影響,進口壓力二維模擬與三維模擬結(jié)果有較大差別,三維模擬結(jié)果與物理模型結(jié)果吻合得更好,但在水躍區(qū)二維模擬與三維模擬得出的壓力相差不大。壩面壓力的數(shù)學(xué)模擬結(jié)果與物理模型結(jié)果相比還存在一定的差距,主要原因是數(shù)學(xué)模型是在物理模型的基礎(chǔ)上對進口條件進行簡化后建立的。進口條件的簡化對計算結(jié)果造成的影響在二維模擬中表現(xiàn)得非常突出。由于沒有閘墩、進口側(cè)收縮的影響,堰頂上游段水流流速比物理模型試驗所測得的流速大,因此造成二維模擬結(jié)果中進口段壩面壓力偏小。

4 結(jié) 語

以銀盤水電站帶閘墩溢流壩為例,用紊流模型模擬了控泄工況下二維和三維數(shù)值模擬水面線位置、壩面壓力分布以及流線分布情況,并與物理模型結(jié)果進行了對比。與二維模擬結(jié)果相比,三維模擬結(jié)果與物理模型試驗結(jié)果吻合程度較好,特別是對進口壓力與水躍發(fā)生點的模擬。但由于數(shù)學(xué)模型進出口邊界條件與物理模型存在差異,數(shù)學(xué)模擬結(jié)果與物理模型試驗結(jié)果存在一定的差異。相比二維模擬中進出口邊界條件簡化造成的影響,三維模擬受到的影響較少,因此推薦采用三維模型進行模擬。

:

[1]蘇遠波,馬震岳,張宏戰(zhàn),等.太平哨水電站溢流壩閘墩裂縫成因數(shù)值分析[J].水利水電科技進展,2009,29(S1):103-106.

[2]陳群,戴光清,劉浩吾.階梯式溢流壩研究綜述[J].水利水電科技進展,2002,22(1):47-50.

[3]魏文禮,劉玉玲,金忠青.溢流壩面明流流動水力特性數(shù)值分析[J].河海大學(xué)學(xué)報:自然科學(xué)版,2000,28(3):107-110.

[4]譚仁忠,王世夏.有閘墩的溢流壩過壩水流數(shù)學(xué)模型[J].水利學(xué)報,1988,19(8):40-44.

[5]馬福喜,王金瑞.三維紊流數(shù)值研究[J].水動力學(xué)研究與進展:A輯,1995,10(2):115-124.

[6]王志東,汪德.含閘墩溢流壩三維過壩水流數(shù)值模擬[J].水科學(xué)進展,2004,15(6):735-738.

[7]王志東,汪德.不同墩型溢流壩水動力特性研究[J].河海大學(xué)學(xué)報:自然科學(xué)版,2004,32(4):391-394.

[8]陳群,戴光清,劉浩吾.帶有曲線自由水面的階梯溢流壩面流場的數(shù)值模擬[J].水利學(xué)報,2002,33(9):20-26.

[9]陳永明,吉慶豐.階梯溢流壩水流數(shù)值模擬及消能分析[J].灌溉排水學(xué)報,2006,25(2):68-71.

[10]劉曉川,楊紀偉,李書芳,等.階梯溢流壩自摻氣水流的數(shù)值模擬[J].水利科技與經(jīng)濟,2008,14(3):190-192.

[11]程香菊,羅麟,趙文謙,等.階梯溢流壩自由表面摻氣特性數(shù)值模擬[J].水動力學(xué)研究與進展:A輯,2004,19(2):152-157.

[12]FERZIGER J H,PERIC M.Computational methods for fluid dynamics[M].Berlin:Springer,1996.

[13]YAKHOT V,ORSZAG S A.Renomalization group analysis of turbulence[J].Journal of Scientific Computing,1986,1(1):3-5.

[14]LAUNDER B E,SPALDING D B.The numerical computation of turbulent flows[J].Computer Methods in AppliedMechanics and Engineering,1974,3:269-289.

[15]LAUNDER B E,SPALDING D B.Lectures in mathematical models of turbulence[M].London:Academic Press,1972.

[16]HIRT CW,NICHOLS BD.Volume of fluid(VOF)method for the dynamics of free boundaries[J].J Comput Phys,1981,39:201-225.