查日扣水電站豎井旋流泄洪洞水力學(xué)數(shù)值模擬研究

南 洪，賀 威，韓鵬輝，張麗花

(中國電建集團 西北勘測設(shè)計研究院有限公司，陜西 西安710065)

豎井旋流消能結(jié)構(gòu)是使泄洪洞高速水流在消能豎井內(nèi)旋轉(zhuǎn)，增大水流沿程損失，并且通過水流在井底水墊的摻和，達到水流消能的目的[1]，是泄水建筑物一種較新的布置形式[2－4]，近年來沙牌和溪古等水電工程采用了豎井旋流泄洪洞[5－9]。

豎井旋流泄洪洞的研究方法主要有水力學(xué)仿真計算和傳統(tǒng)的模型試驗。因為豎井內(nèi)水流復(fù)雜，常規(guī)的測量方法很難測量渦室及豎井內(nèi)的水力參數(shù)[10]。水力學(xué)仿真計算具有計算效率高、成本低、無比尺效應(yīng)等優(yōu)點。因此，采用數(shù)值仿真計算方法，研究豎井旋流泄洪洞的水力特性，并與水工模型試驗進行對比驗證，可為評判豎井旋流泄洪洞水力學(xué)數(shù)值仿真計算的適應(yīng)性提供參考。

1 數(shù)值計算原理

1.1 k－ε雙方程模型

本計算采用的紊流模型為k－ε模型，連續(xù)性方程為[10－15]:

動量方程為:

紊動動能k方程為:

紊動耗散率ε方程為:

式中:t為時間;ρ和μ分別為密度和分子黏性系數(shù);P為修正的壓力;ui和xi分別為速度和坐標分量;σε和σk分別為ε和k的紊流普朗特數(shù)，計算取σk=1.0，σε=1.3;G為平均速度梯度引起的紊動動能產(chǎn)生項;C1ε和C2ε為ε方程常數(shù)，計算取C1ε=1.44，C2ε=1.92;μt為紊流黏性系數(shù)，可由紊動耗散率 ε和紊動動能k求出:

式中:Cμ為常數(shù)，計算中 Cμ=0.09。

1.2 VOF模型

圖1 豎井旋流泄洪洞體型圖(單位:m)

采用VOF方法的k－ε紊流模型，分子黏性系數(shù)μ以及密度ρ用體積分數(shù)的平均值給出，μ和ρ是體積分數(shù)的函數(shù)，可用下式表示:

式中:αa和αw分別為單元中氣和水的體積分數(shù);ρa和ρw分別為氣和水的密度;μa和μw分別為氣和水的分子黏性系數(shù)。

水和氣的界面通過求解以下方程來完成:

2 數(shù)值計算模型

2.1 豎井旋流泄洪洞設(shè)計體型

查日扣水電站豎井旋流泄洪洞由引水道、渦室和豎井段、退水洞段、原導(dǎo)流洞組成，如圖1所示。引水道和原導(dǎo)流洞夾角為50°，引水道長約100 m，進口底板高程為3 950.00 m，檢修門孔口尺寸3 m×4 m，工作門孔口尺寸3 m×3.5 m，均為平板門，引水道正常運行為全有壓狀態(tài)。豎井段由通氣孔、渦室、漸變段、豎井、消力池組成。退水洞長20 m，上游斷面為3.8 m×7 m，下游斷面為3.8 m×4 m，中間有一長7.5 m，厚1.5 m的中墩。導(dǎo)流洞斷面尺寸為5.5 m×7 m，結(jié)合段長約343 m。校核洪水位為3 992.5 m時，最大設(shè)計泄量為200 m3/s。

2.2 計算模型

豎井旋流泄洪洞計算模型如圖2所示，泄洪洞進水口段、引水道段、旋流豎井段、退水洞和導(dǎo)流洞均按0.5 m剖分單元網(wǎng)格，上、下游水體部分適當(dāng)變疏，按1 m剖分，共有節(jié)點368 433個，單元328 233個。上游水體四周設(shè)為進水邊界;下游水體四周設(shè)為出水邊界;通氣孔頂設(shè)為通風(fēng)邊界，模型其它外邊界均為流道邊界。設(shè)計洪水位為3 992.00 m，相應(yīng)下游水位為3 900.62 m;校核洪水位為3 992.5 m，相應(yīng)下游水位為3 901.34 m。計算采用有限體積法隱式迭代求解方法，壓力場和速度場的耦合采用PISO法。

圖2 計算模型圖

3 計算結(jié)果及分析

3.1 流量結(jié)果

計算泄洪洞設(shè)計及校核水位泄量分別為206 m3/s、208 m3/s，水工模型試驗確定設(shè)計及校核水位泄量分別為210 m3/s、211 m3/s，二者差別較小，均滿足最大泄量200 m3/s流量要求。

3.2 流態(tài)結(jié)果

泄洪洞校核工況流態(tài)計算結(jié)果如圖3所示，設(shè)計工況流量僅比校核工況少2 m3/s。由流態(tài)計算結(jié)果可以看出，水流在進入豎井前的引水道內(nèi)全部為水;進入豎井后，在豎井邊壁附近形成水氣混合體，從上到下水流含氣量逐漸變大，豎井中心為空氣;水流在豎井出水口段為含氣量在45%～55%之間的水氣混和體;導(dǎo)流洞頂部為空氣，底部為含氣量在45%～55%之間的水氣混和體。

圖3 泄洪洞流態(tài)計算結(jié)果云圖

水工模型試驗結(jié)果表明，校核工況時泄洪洞引水洞段為滿流;觀察水流在豎井內(nèi)的流態(tài)，也表明水流摻氣量從上到下逐漸增大，但其含氣量值很難測量;導(dǎo)流洞內(nèi)頂部為空氣，底部為水氣混和體，在最大水深時，其摻氣水面線較低。

3.3 流速結(jié)果

校核工況流速計算結(jié)果如圖4所示，由計算結(jié)果可以看出，水流在進入豎井前，流速在13 m/s～21 m/s范圍之間，平均流速約19 m/s;水流進入豎井后，在豎井進水口部位以水平向流速為主，平均流速約13 m/s，在豎井中部以豎直向流速為主，最大值約20 m/s;水流在豎井出水口部位的流速在13 m/s～20 m/s范圍之間，平均流速約16 m/s;水流在導(dǎo)流洞內(nèi)的流速從上游到下游逐漸減小，在導(dǎo)流洞末端流速約13 m/s;通氣孔風(fēng)速約45 m/s。

水工模型試驗結(jié)果表明，校核水位時引水道段的平均流速在11.9 m/s～20.0 m/s之間;退水洞孔口至泄0+322.562 m段，斷面最大流速為15.4 m/s～18.6 m/s;泄0+322.562 m～泄洪洞末端(泄0+472.562 m)，斷面最大流速為13.4 m/s～14.2 m/s，流速沿程減小;通氣孔風(fēng)速試驗很難測量;說明水工模型試驗流速結(jié)果與數(shù)值計算差別較小。

圖4 泄洪洞流速計算結(jié)果矢量圖

3.4 壓強結(jié)果

校核工況壓強計算結(jié)果如圖5所示，由計算結(jié)果可以看出，水流在泄洪洞工作閘門以前以及豎井消力池底部壓強較大，其它部位壓強均較小。在工作門以前，壓強范圍為3 kPa～300 kPa之間，隨水流流向逐漸減小;豎井消力池底部壓強最大值約377 kPa，在豎井和出水口交接處的頂部有最大值約－30.8 kPa的壓強。

水工模型試驗結(jié)果表明，泄洪洞內(nèi)水流壓強分布規(guī)律與數(shù)值計算結(jié)果基本一致，豎井消力池底部壓強最大值約415 kPa;在豎井和出水口交接處的頂部有最大值約－37.2 kPa的壓強，說明水工模型試驗壓強與數(shù)值計算差別較小。

圖5 泄洪洞壓強計算結(jié)果云圖

3.5 流線結(jié)果

校核工況流線計算結(jié)果如圖6所示，由計算結(jié)果可以看出，由進水口進入豎井的水體，其上表部在豎井內(nèi)旋轉(zhuǎn)約兩圈，底部在豎井內(nèi)旋轉(zhuǎn)約一圈，與水工模型試驗成果基本一致。

圖6 泄洪洞上游水體流線計算結(jié)果

4 結(jié)論

通過對查日扣水電站豎井旋流泄洪洞進行水力學(xué)三維數(shù)值計算，以及與水工模型試驗成果的對比驗證，結(jié)論如下:(1)水力學(xué)數(shù)值計算在流量、流速、壓強等方面與水工模型試驗成果誤差較小，計算精度比較高;(2)在計算水流摻氣量、通風(fēng)流速、水流流線等模型試驗很難測量的方面，數(shù)值仿真計算更有優(yōu)勢，能簡單直觀的給出結(jié)果;(3)數(shù)值計算在計算水流的脈動壓力、河道的沖淤變形等方面還有很大的不足，要參考模型試驗成果。

[1]安盛勛，王君利，陳念水，等.水平旋流消能泄洪洞設(shè)計與研究[M].北京:中國水利水電出版社，2008.

[2]董興林，高季章，魯慎吾，等.導(dǎo)流洞改為旋渦式豎井溢洪道綜合研究[J].水力發(fā)電，1995，22(3):32-37.

[3]郭琰，倪漢根.旋流式豎井溢洪道的過流能力及蝸室與錐形漸縮段的水力特性研究[J].水動力學(xué)研究與進展，1995，10(1):97-105.

[4]郭 琰，倪漢根.旋流式豎井溢洪道豎井的水流特性研究[J].水動力學(xué)研究與進展，1995，10(2):146-154.

[5]李貴信，周 恒.公伯峽水電站導(dǎo)流洞改建為旋流泄洪洞方案研究[J].水力發(fā)電，2002，29(8):41-44.

[6]劉善均，許唯臨，王 韋.高壩豎井泄洪洞、龍?zhí)ь^放空洞與導(dǎo)流洞“三洞合一”優(yōu)化布置研究[J].四川大學(xué)學(xué)報，2003，35(2):10-14.

[7]邵敬東.漩流式豎井泄洪洞在沙牌工程中的應(yīng)用[J].水電站設(shè)計，2003，19(4):61-68.

[8]方衛(wèi)紅，李學(xué)強.溪古水電站豎井旋流溢洪洞設(shè)計與分析[J].西北水電，2012(4):41-45.

[9]林崇勇.豎井旋流泄洪洞在甲巖水電站的應(yīng)用[J].云南水力發(fā)電，2013(5):26-28.

[10]張曉東，劉之平，高季章，等.豎井旋流泄洪洞數(shù)值模擬[J].水利學(xué)報，2003，34(8):58-63.

[11]楊朝暉，吳守榮，余 挺，等.豎井旋流泄洪洞三維數(shù)值模擬研究[J].四川大學(xué)學(xué)報，2007，39(2):41-45.

[12]錢小燕，鄧 軍，張葉林，等.豎井旋流泄洪洞渦室整流挑坎試驗與數(shù)值模擬[J].人民長江，2011，42(6):20-22.

[13]付 波，牛爭鳴，李國棟，等.豎井進流水平旋轉(zhuǎn)內(nèi)消能泄洪洞水力特性的數(shù)值模擬[J].水動力學(xué)研究與進展，2009(3):164-171.

[14]陳華勇，鄧 軍，胡 靜，等.旋流式豎井泄洪洞水力要素的數(shù)值模擬與試驗研究[J].水力發(fā)電，2008(3):79-82.

[15]牛爭鳴，程慶迎，譚立新.豎井進流水平旋轉(zhuǎn)內(nèi)消能泄洪洞流場數(shù)值模擬[J].西安理工大學(xué)學(xué)報，2005，21(2):113-117.