汪 厚 松， 刁 明 軍， 蔣 雷， 王 軍， 官 夏 菲， 任 昱

(四川大學(xué)水力學(xué)與山區(qū)河流開發(fā)保護國家重點實驗室，四川 成都 610065)

1 研究背景

在某些水利工程中，由于受地形和地質(zhì)條件的限制，泄水建筑物不可避免地布置成彎道形式，而泄水建筑物中的水流通常是高速水流，水流屬于急流狀態(tài)，因此彎道急流是水利工程中常見的水流現(xiàn)象。許多專家和學(xué)者[1-3]發(fā)現(xiàn)了彎道水流現(xiàn)象并對彎道水流的機理進行了深入研究。王韋、許唯臨[4]等對彎道水沙的運動進行了理論研究；王平義、蔡金德[5]等對彎曲河道的縱向垂線平均流速的平面分布進行了研究，推導(dǎo)出了變曲率動床彎道水流縱向垂線平均流速計算公式；刁明軍、楊海波[6]等對彎道水流的研究現(xiàn)狀進行了闡述，指出對彎道水流進行數(shù)值模擬研究將是一種趨勢；楊玲霞、張銀華[7]等通過水工模型試驗對改善彎道急流流態(tài)進行了研究；惠愛瑙、王飛虎[8]等采用水工模型試驗對大底坡急流彎道通氣減蝕所形成的摻氣空腔長度和水體摻氣濃度進行了研究；吳宇峰、伍超[9]等通過模型試驗對斜檻在急流彎道超高控制進行了研究，闡明了斜檻防治超高的原理，提出了斜檻的設(shè)計方法；許光祥、童思陳[10]等對單純彎道下的水面橫比降沿程分布進行了研究，得出了彎道橫比降沿程分布特性及其統(tǒng)一計算公式。橫向超高是彎道水流的一項重要指標(biāo)，改善彎道水流橫向超高具有重要的工程意義。在彎道前增設(shè)摻氣坎對彎道急流橫向超高的影響還缺乏研究，由于彎道的橫向超高主要因彎道縱向流速沿橫向分布不均勻所引起的，前置摻氣坎的設(shè)置必將影響到彎道縱向流速沿橫向分布，進而影響彎道橫向超高。因此，有必要研究前置摻氣坎對彎道急流橫向超高的影響。

肖鴻、周赤[11]等采用數(shù)值模擬和模型試驗對摻氣減蝕設(shè)施進行了研究，將水深和流速的模擬結(jié)果與模型試驗數(shù)據(jù)進行對比得出最大誤差不超過2%，證明了Eulerian模型和VOF模型模擬摻氣設(shè)施的可行性；羅永欽、刁明軍[12]等采用模型試驗和數(shù)值模擬對泄洪洞摻氣減蝕設(shè)施進行了研究，水面線的模擬結(jié)果最大誤差保持在3%以內(nèi)，數(shù)值模擬表明k-ε模型和VOF方法能夠準(zhǔn)確地對泄洪洞摻氣減蝕問題進行模擬與分析，特別是對自由水面捕捉和摻氣空腔的長度模擬精度較高；葉茂、伍平[13]和漆力健[14]等采用數(shù)值模擬的方法對摻氣坎后水流進行了研究，進一步證明了VOF方法在模擬自由水面和空腔長度的準(zhǔn)確性。基于數(shù)值模擬方法研究摻氣水流積累的成功經(jīng)驗，本研究以某工程的排洪隧洞的彎道為研究對象，采用標(biāo)準(zhǔn)k-ε模型與VOF模型對不同半徑下的彎道在有無有前置摻氣坎的情況下進行了數(shù)值模擬研究。

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

本文采用的標(biāo)準(zhǔn)k-ε紊流模型，其方程組[15]如下：

(1)

(2)

(3)

(4)

其中，ρ和μ分別為容積分?jǐn)?shù)平均的密度和動力粘性系數(shù),P考慮重力后的壓強項，μt為紊動粘性系數(shù)，可以由紊動能k和紊動能耗散率ε按下式求出：

(5)

式中Cμ為經(jīng)驗常數(shù);σk和σε分別為k和ε的紊流Prandtl數(shù)，C1ε和C2ε為ε方程常數(shù)。G為平均速度梯度引起的紊動能產(chǎn)生項，由下面的公式定義：

(6)

各常數(shù)取值見表1：

表1 紊流模型中常數(shù)的取值

采用VOF方法追蹤水面。該方法認(rèn)為計算域內(nèi)水氣交界面的網(wǎng)格空間由氣相、水相混合充填，密度ρ和動力粘性系數(shù)μ為氣相(水相)容積分?jǐn)?shù)的函數(shù)：

ρ=αwρw+(1-αw)ρa

(7)

μ=αwμw+(1-αw)μa

(8)

式中αw為水的容積分?jǐn)?shù)，ρw、ρa分別為水和空氣的密度，μw、μa分別為水和空氣的動力粘性系數(shù)。通過對水的容積分?jǐn)?shù)αw的迭代求解，ρ和μ可以通過公式(7)、公式(8)求出。

本次模擬的區(qū)域從上游庫區(qū)到彎道結(jié)尾直線段末端。溢流堰頂處樁號為(渠)0+000，摻氣坎處樁號為(渠)0+320，彎道起始處樁號為(渠)0+370，明渠收縮段末端以下寬度均為5 m，如圖1所示。摻氣坎上游底坡為0.11，下游底坡為0.08，如圖2所示。模型區(qū)域的整體網(wǎng)格數(shù)量為35萬。由于彎道處為本次研究的重點，對彎道區(qū)域進行了加密處理，如圖1所示。

圖1 計算域三維模型

圖2 摻氣坎尺寸圖

3 計算工況

本次數(shù)值模擬計算工況如表2所示：

4 結(jié) 果

由于數(shù)據(jù)提取和結(jié)果分析的需要，對于轉(zhuǎn)彎半徑為R=150 m，轉(zhuǎn)彎角θ=47°的彎道，從彎道起始處開始，彎道全程等間距設(shè)置21個監(jiān)測斷面(編號為1到21)，各監(jiān)測斷面之間間距為6.15 m；對于轉(zhuǎn)彎半徑為R=500 m，轉(zhuǎn)彎角θ=26°的彎道，從彎道起始處開始，彎道全程等間距設(shè)置37個監(jiān)測斷面(編號為1到37)，監(jiān)測斷面間距為6.33 m。由于彎道區(qū)域的流速較高，屬于高速水流，彎道區(qū)域的流速均在20 m/s以上，存在強烈的摻氣現(xiàn)象，故本文選取水氣界面的體積分?jǐn)?shù)為0.2的等值線作為自由水面，并對每個監(jiān)測斷面的左岸水深hl和右岸水深hr進行了提取，得到各工況下彎道段的橫向超高Δh=hl-hr沿程分布，如圖3所示。

表2 計算工況

圖3 橫向超高沿程分布

((9)

5 研 究

由圖3可知，各彎道橫向超高沿程呈波狀分布，對于R=150 m的彎道，1～6號監(jiān)測斷面區(qū)域橫向超高沿程逐漸增加，6～13號監(jiān)測斷面橫向超高沿程超高逐漸減小，13～19號監(jiān)測斷面橫向超高沿程逐漸增加，19～21號監(jiān)測斷面橫向超高逐漸減小，波長特性約為73.8 m；對于R=500 m彎道而言，同樣表現(xiàn)出類似于R=150 m時的沿程變化規(guī)律，相應(yīng)的波長特性約為101.28 m。增設(shè)摻氣坎和未設(shè)摻氣坎的橫向超高沿程波狀分布同步，同一彎道不同流量下的波長特性相同，波幅特性隨流量的減小而減小，波長特性主要受邊壁條件控制。彎道內(nèi)的急流由于彎道邊壁偏轉(zhuǎn)而產(chǎn)生擾動，所形成的沖擊波經(jīng)過彎道壁面反射向下游傳播，導(dǎo)致橫向超高沿程呈波狀分布[16]。擾動因邊壁的變化而產(chǎn)生，并經(jīng)邊壁反射向下游傳播，因此，波長特性主要受邊壁條件控制。增設(shè)摻氣坎對彎道橫向超高沿程波狀分布規(guī)律基本無影響。

圖4 各工況下橫向超高平均降低值

在彎道前增設(shè)置摻氣坎后，各個監(jiān)測斷面處的橫向超高比未增設(shè)摻氣坎時對應(yīng)工況下的橫向超高有所降低，橫向超高降低值在橫向超高沿程的第一個單調(diào)區(qū)間內(nèi)表現(xiàn)不顯著，基本為0，隨后向下游逐漸表現(xiàn)顯著，在R=150 m的彎道中橫向超高最大降低值為0.507 m(1 000年一遇)、0.529 m(100年一遇)和0.275 m(5年一遇)，在R=500 m的彎道中橫向超高最大降低值為0.436 m (1 000年一遇)、0.405 m (100年一遇)和0.376 (5年一遇)，如圖3所示。彎道斷面的水面形態(tài)取決于彎道縱向流速沿橫向分布，由于水流經(jīng)過摻氣坎到跌落至底板的過程中，增加了水體與空氣以及水體與水體之間的剪切摩擦，改變了水體原有的狀態(tài)，使進入彎道的水體流場進行了重新分布，彎道縱向流速沿橫向分布相對于未增設(shè)前置摻氣坎時的流場發(fā)生了變化，強迫摻入的大量氣體沿程外溢，使水流的紊動增強，改變了彎道區(qū)域原有流場，降低了流速分布的不均勻性，從而降低了彎道區(qū)域的橫向超高。

由于監(jiān)測斷面取值存在離散性，因此本研究采用橫向超高的平均降低值來反映前置摻氣坎對彎道急流橫向超高的整體改進程度。由圖4可知，對R=150 m的彎道而言，摻氣坎對彎道橫向超高的整體改進程度在100年一遇工況時最大為0.26 m，在1 000年一遇工況時次之為0.24 m，在5年一遇工況時最小為0.14 m；對R=500 m的彎道而言，摻氣坎對彎道橫向超高的整體改進程度在5年一遇工況時最大為0.12 m，在1 000年一遇工況時次之為0.11 m，在100年一遇工況時最小為0.09 m。每種工況下，摻氣坎對R=150 m的彎道的改進程度比對R=500 m的彎道改進的程度要大。在同一摻氣坎作用下，彎道半徑越小，彎道橫向超高整體改進效果越顯著；摻氣坎對同一彎道的橫向超高改進效果還受流量的影響。

6 結(jié) 論

本研究運用數(shù)值模擬的方法，對不同半徑的彎道急流在有無前置摻氣坎條件下的橫向超高進行了研究，得出了以下結(jié)論:

(1)在彎道前增設(shè)摻氣坎可降低彎道急流的橫向超高，并且對彎道急流橫向超高的沿程波狀分布規(guī)律基本無影響，此措施可為類似的工程提供改善彎道急流橫向超高的經(jīng)驗；

(2)前置摻氣坎對彎道急流橫向超高的整體改進程度與彎道半徑有關(guān)，半徑越小，改進程度越顯著，摻氣坎對同一彎道橫向超高的改進程度還受流量的影響。

[1] 羅索夫斯基著，尹學(xué)良譯.彎道水流的研究[J].泥沙研究，1958，3(1)：83-95.

[2] THOMSONJ.On the origin and winding of rivers in alluvialplains[M].Proc.Royal Society of London，1876，25(3):83-88.

[3] 劉煥芳.彎道自由水面形態(tài)的研究[J].水利學(xué)報，1990(4)：46-50.

[4] 王韋，許唯臨，等.彎道水沙運動理論及應(yīng)用[M].成都：成都科技大學(xué)出版社，1994(7)：27-50.

[5] 王平義，蔡金德，等.彎曲河道縱向垂線平均流速平面分布的研究[J].水動力學(xué)研究與進展，1994(6)：267-275.

[6] 刁明軍，楊海波，等.彎道水力學(xué)研究現(xiàn)狀與進展[J].西南民族大學(xué)學(xué)報(自然科學(xué)版)，2007，33(3)：597-599.

[7] 楊玲霞，張銀華，等.改善彎道急流流態(tài)的實驗研究[J].中國農(nóng)村水利水電，2006(7)：77-79.

[8] 惠愛瑙，王飛虎，等.大底坡急流彎道通氣減蝕水力特性研究[J].西北水資源與水工程，1997，8(2)：55-59.

[9] 吳宇峰，伍超，等.斜檻在急流彎道控制超高的設(shè)計研究[J].水力發(fā)電學(xué)報，2007，26(3)：78-81.

[10] 許光祥，童思陳，等.彎道水面橫比降沿程分布特性研究[J].水力發(fā)電學(xué)報，2009，28(4):115-118.

[11] 肖鴻，周赤，等.摻氣減蝕設(shè)施數(shù)值模擬研究探討[J].長江科學(xué)院院報，2014，31(10)：115-119.

[12] 羅永欽，刁明軍，等.高壩明流泄洪洞摻氣減蝕三維數(shù)值模擬分析[J].水科學(xué)進展，2012，23(1)：111-115.

[13] 葉茂，伍平，等.泄洪洞摻氣水流的數(shù)值模擬研究[J].水力發(fā)電學(xué)報，2014，33(4)：106-108.

[14] 漆力健，廖華勝，等.低弗勞德數(shù)摻氣坎空腔回水問題研究[J].水利學(xué)報，2007，38(7)：819-825.

[15] 王福軍.計算流體動力學(xué)分析——CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社，2004：120-123.

[16] 吳持恭.水力學(xué)：下冊(第四版)[M].北京：高等教育出版社，2008：311-324.