柳 楊,馬 飛,吳建華

(河海大學(xué)水利水電學(xué)院,江蘇 南京 210098)

窄縫坎是收縮式消能工的一種形式,它利用泄水建筑物末端兩側(cè)邊墻的收縮,使水流經(jīng)過窄縫坎后形成沿豎向和縱向充分拉伸的挑射水舌,增大水舌在空中的摻氣和消能,減輕對下游河床的沖刷。窄縫坎因其高消能率有廣泛的應(yīng)用,特別適用于高窄峽谷的水電工程。我國的龍羊峽、東江、安康、李家峽、天生橋一級、漫灣、水布埡、善泥坡等水電站都采用了窄縫坎的消能方式[1-2]。

明渠邊墻偏轉(zhuǎn)會(huì)引起沖擊波,窄縫坎因其兩側(cè)邊墻的收縮,改變了流線方向,試驗(yàn)中產(chǎn)生了沖擊波。早在1943年,Ippen[3]以力學(xué)原理為基礎(chǔ)建立了理想沖擊波理論,提出了沖擊波的基本關(guān)系式。Hager等[4-6]通過理論推導(dǎo)得到了沖擊波的簡化計(jì)算式。但是,前期的研究幾乎都依賴于靜水壓力假定,考慮動(dòng)水壓力影響的沖擊波研究較少。事實(shí)上,窄縫坎后動(dòng)水壓力較大,不再符合靜水壓力假定,劉亞坤[7]在Ippen關(guān)系式中的靜水壓力項(xiàng)引入了動(dòng)水壓力修正系數(shù),并通過物理模型試驗(yàn)對其進(jìn)行了率定。寧利中等[8]認(rèn)為窄縫坎流動(dòng)條件與明渠收縮段有明顯差別,窄縫坎收縮段長度對水流沖擊波波角產(chǎn)生一定影響,他們根據(jù)模型試驗(yàn)結(jié)果給出了計(jì)算波角的經(jīng)驗(yàn)公式。黃智敏等[9]分析認(rèn)為窄縫坎的波陣線為曲線,沖擊波的波角不是一個(gè)常數(shù),而是沿程增大的。可見,與明渠相比，窄縫坎偏轉(zhuǎn)角大、水流弗勞德數(shù)高,產(chǎn)生的沖擊波難以直接用傳統(tǒng)沖擊波理論進(jìn)行求解,有必要對窄縫坎沖擊波開展進(jìn)一步的研究。

在消力池的設(shè)計(jì)中,其結(jié)構(gòu)尺寸根據(jù)挑流入池水流的長度和寬度來確定。除了傳統(tǒng)挑流水舌運(yùn)動(dòng)軌跡的研究[10-12]和水舌橫向擴(kuò)散特性的研究[13-14],吳文平等[15]運(yùn)用自由拋射體理論研究了窄縫坎挑流水舌的縱向擴(kuò)散規(guī)律,得到了水舌入水的長度。此外,由于窄縫坎內(nèi)沖擊波交匯，在水舌頂部形成水冠,該部分水流在流動(dòng)過程中從水舌主體兩側(cè)散落,決定了水舌入水時(shí)的最大寬度。本文通過理論分析和物理模型試驗(yàn),研究窄縫坎沖擊波最大和最小波角;應(yīng)用沖擊波理論和自由拋射體運(yùn)動(dòng)原理,建立窄縫坎水舌最大入水寬度及其位置的估算公式,并與模型試驗(yàn)結(jié)果進(jìn)行比較。

1 試驗(yàn)裝置與測量參數(shù)

試驗(yàn)在河海大學(xué)高速水流實(shí)驗(yàn)室進(jìn)行。試驗(yàn)?zāi)Ｐ拖到y(tǒng)包括水泵、平水塔、鋼板水箱、窄縫坎、尾水渠及回流系統(tǒng)等(圖1)。模型試驗(yàn)以某工程泄洪中孔為原型,模型比尺1∶40。該中孔模型采用短有壓進(jìn)口,有壓段長64.5 cm,有壓段出口頂板設(shè)置1∶10壓坡,出口斷面尺寸為15.0 cm×18.3 cm(寬×高),明渠段長48.3 cm,寬15.0 cm,邊墻高38.0 cm,明渠末端為水平底板窄縫坎,底板距下游尾水渠渠底高147.0 cm,兩側(cè)邊墻寬度由B=15.0 cm收縮至b=8.0 cm。尾水渠中設(shè)尾坎以模擬水墊塘,尾坎高67.0 cm,距模型中孔出口723.0 cm。窄縫坎結(jié)構(gòu)參數(shù)中,收縮角θ是沖擊波波角的主要影響因素,本文研究了收縮角θ分別為17.30°、12.0°和8.0°的情況,相應(yīng)的模型編號(hào)為M1、M2和M3。試驗(yàn)中工作水頭為80.0～172.0 cm,來流弗勞德數(shù)Fr=2.17～3.85。

圖1 模型試驗(yàn)系統(tǒng)示意圖

試驗(yàn)涉及流量、沖擊波波角、挑流水舌最大寬度以及最大寬度所在位置等參數(shù)，其中,流量采用矩形薄壁堰測量,水舌最大寬度及其相應(yīng)的位置用鋼尺測量,沖擊波波角由試驗(yàn)測得的相關(guān)參數(shù)計(jì)算求得。試驗(yàn)中挑射水舌有一定的擺動(dòng),相應(yīng)地，沖擊波波角、挑射水舌最大入水寬度及其位置的測量誤差約為±5%。

2 流態(tài)觀察

圖2為M1模型窄縫坎水舌形態(tài)的試驗(yàn)照片,可以看到,水流進(jìn)入窄縫坎收縮段后,由于窄縫坎寬度束窄,產(chǎn)生沖擊波,水深沿程迅速增大;窄縫坎挑流水舌在縱向和豎向上充分拉伸,水舌與空氣接觸面積大,摻氣充分,整個(gè)水舌呈白色;受沖擊波影響,水舌頂部形成摻氣更為劇烈的水冠,其顏色顯得更白些。圖3是不同來流條件下M1模型水流俯視照片,可以看出,在沖擊波交匯點(diǎn)前,窄縫坎內(nèi)靠近邊墻的水面高于中線水面,即水流表面是從中線到兩側(cè)邊墻逐漸上升的曲面;窄縫坎沖擊波波陣線不是水平的,而是沿流向傾斜向上的;與明渠偏轉(zhuǎn)引起的沖擊波不同,隨著窄縫坎收縮段內(nèi)水面逐漸升高,沖擊波波陣線不是一條直線,而是一組直線;窄縫坎內(nèi)靠近兩側(cè)邊墻的水流形成的波陣線最高,對應(yīng)波角最小,這兩股水流交匯后相互穿越并帶著中線上的水流從主體水舌兩側(cè)跌落,做近似自由拋物體運(yùn)動(dòng),由于主體水舌厚度較薄,這兩股水流入水時(shí)所確定的寬度(圖2中點(diǎn)A處的橫向?qū)挾?即為水舌入水時(shí)的最大寬度。試驗(yàn)表明,在所有來流條件下,每組模型都呈現(xiàn)類似的水流流態(tài)。

圖2 M1模型窄縫坎水舌形態(tài)(Fr=3.8)

圖3 M1模型水流頂部流態(tài)(沿圖2中箭頭方向拍攝)

3 水舌入水寬度的計(jì)算

根據(jù)試驗(yàn)現(xiàn)象,繪制了窄縫坎及流態(tài)示意圖,如圖4所示。圖中,窄縫進(jìn)口寬度為B,出口斷面寬度為b,收縮段長度為L,收縮角為θ。當(dāng)水流以速度v1流經(jīng)收縮段時(shí),產(chǎn)生沖擊波,波陣線在水平面的投影和來流方向的夾角為波角,將最大波角記為βM,最小波角記為βm。O′點(diǎn)為沖擊波交匯點(diǎn),從O′點(diǎn)到收縮段起始斷面的水平距離定義為L0,O′點(diǎn)水深為h0。

圖4 窄縫坎及流態(tài)示意圖

為方便計(jì)算,在圖4中建立坐標(biāo)系xOy。坐標(biāo)原點(diǎn)O位于收縮段起始處,距離底板高度為h1,x軸為水流方向,y軸豎直向上。窄縫坎底板距下游水面的高度為z,取z=80.0 cm。

如圖4所示,在沖擊波作用下,窄縫坎收縮段兩側(cè)邊墻處的兩股水流在水舌頂部夾角為2βm。中線上的頂部水流由于劇烈的沖擊作用,以最大的出射角α并跟隨邊墻附近的兩股水流從主體水舌兩側(cè)跌落,落點(diǎn)為A,則A點(diǎn)處水舌的橫向?qū)挾燃礊檎p坎水舌最大入水寬度。記此時(shí)落點(diǎn)A與沖擊波交匯點(diǎn)O′的水平距離為L1,相應(yīng)的水流落入下游時(shí)的水舌寬度,記為BM,即挑流水舌入水時(shí)的最大寬度,定義窄縫坎出口距水舌最大入水寬度處的水平距離為LM。則由幾何關(guān)系可得

(1)

根據(jù)式(1)即可求得水舌最大入水寬度BM。在計(jì)算時(shí),需要先求解βm和L1。

3.1 波角βM和βm的計(jì)算

窄縫坎邊墻偏轉(zhuǎn)角較大,沖擊波后不再符合靜水壓力假定,直接用Ippen理論研究窄縫挑坎中的沖擊波誤差較大,需要進(jìn)行修正。Hager 等[5]認(rèn)為,窄縫坎收縮段內(nèi)形成的沖擊波波角可以表示為來流弗勞德數(shù)Fr和收縮角θ的函數(shù):

β=f(1/Fr,θ)

(2)

通過分析本文窄縫坎系列模型試驗(yàn)數(shù)據(jù),在不同來流弗勞德數(shù)條件下分別得到了βM和βm與Fr和θ的關(guān)系,如圖5和圖6所示。

圖5 βm-0.65θ與1/Fr的關(guān)系

圖6 βM-1.24θ與1/Fr的關(guān)系

根據(jù)式(2),并利用模型試驗(yàn)資料可得

(3)

(4)

式(3)(4)的適用范圍為:0.139≤θ≤0.302,2.17≤Fr≤3.85。

3.2 L1的近似計(jì)算

忽略空氣阻力,把挑射水流的運(yùn)動(dòng)看作質(zhì)點(diǎn)的自由拋射體運(yùn)動(dòng),L1可近似按自由拋射體的運(yùn)動(dòng)軌跡來計(jì)算。根據(jù)自由拋射體運(yùn)動(dòng)原理可以求得

(5)

式中:v、α分別為沖擊波交匯點(diǎn)O′處水流的流速及出射角;h0為交匯點(diǎn)O′處水深。確定了v、α及h0后,可得到L1的值。

3.2.1 v的求解

Ippen沖擊波理論給出了沖擊波波后水流流速:

(6)

由式(6)可知,當(dāng)來流流速及收縮角一定時(shí),波后流速的大小取決于波角。經(jīng)試算,近似以水沖擊波波角βM的值代入式(6),計(jì)算沖擊波波后水流流速,即

(7)

3.2.2 h0的求解

為了求解h0,列來流斷面和O′斷面的連續(xù)性方程,并以O(shè)′點(diǎn)的流速v代表O′斷面的平均流速。由于O′點(diǎn)距離窄縫坎出口很近,該斷面處水舌橫向收縮和底部向下擴(kuò)散很小,因此,將O′斷面近似認(rèn)為是寬B0、高h(yuǎn)0的矩形斷面,則

v1h1B=vh0B0

(8)

式中:v1和h1分別為來流斷面流動(dòng)的平均流速和水深;B0為O′斷面處水流寬度,并且近似B0≈b,則由式(8)可得

(9)

3.2.3 α的求解

在xOy坐標(biāo)系中(圖4),求解沖擊波交匯點(diǎn)處的水流出射角α,圖中虛線為窄縫坎中線截面上的水面線,可用拋物線函數(shù)x2=2py近似表達(dá)[16]。在出射角α的計(jì)算中,近似以水深h0代替O′點(diǎn)到出口底板的高度,O點(diǎn)坐標(biāo)(0,0),O′點(diǎn)坐標(biāo)為(L0,h0-h1),則p=L02/(2(h0-h1),進(jìn)而得到拋物線方程為x2=L02y/(h0-h1)。對方程在O′點(diǎn)求導(dǎo)，得該點(diǎn)的斜率tanα=2(h0-h1)/L0,因此,出射角α可以近似表示為

(10)

將式(7)(9)(10)代入式(5),可以計(jì)算得出L1的值。最后,將式(3)和(5)代入式(1),可得到窄縫坎挑流水舌最大入水寬度BM及最大寬度位置LM。

4 試驗(yàn)驗(yàn)證

采用河海大學(xué)高速水流實(shí)驗(yàn)室窄縫坎模型試驗(yàn)系統(tǒng)來驗(yàn)證上述計(jì)算過程，其中窄縫坎來流流量已知,h1可近似等于有壓段出口高度18.3 cm,根據(jù)流量-水位關(guān)系和結(jié)構(gòu)參數(shù),可得到來流流速v1及水流弗勞德數(shù)Fr。利用本文的計(jì)算方法求得水舌最大入水寬度BM和窄縫坎出口距水舌最大入水寬度處的水平距離LM，在模型試驗(yàn)中測量每種工況下的BM以及LM，將BM和LM的計(jì)算值及實(shí)測值點(diǎn)繪在圖7和圖8中，可以看出,兩者吻合較好，表明本文計(jì)算方法理論基礎(chǔ)可靠，計(jì)算過程簡便，具有較高的精度。

圖7 BM實(shí)測值與計(jì)算值的比較

圖8 LM實(shí)測值與計(jì)算值的比較

5 結(jié) 語

窄縫坎是一種常用的水力消能單元,試驗(yàn)表明,通過窄縫坎的流動(dòng)受沖擊波的作用,在其收縮段表面產(chǎn)生多條傾斜向上的波陣線,對應(yīng)著不同的波角,根據(jù)本文試驗(yàn)結(jié)果得到,波角隨來流弗勞德數(shù)的減小及窄縫坎收縮角的增大而增大。水舌頂部的沖擊波交匯后分成左右兩股水流向下游跌落,決定了水舌的最大入水寬度。試驗(yàn)研究表明,隨著收縮角或來流弗勞德數(shù)的增大,窄縫坎挑射水舌最大入水寬度增大且其位置向下游偏移。本文通過理論分析和試驗(yàn)研究,得到了窄縫坎挑射水舌最大入水寬度及其所在位置的估算公式。經(jīng)與模型試驗(yàn)結(jié)果比較,兩者吻合情況較好，表明本文計(jì)算方法具有較高的精度。

參考文獻(xiàn)：

[1] WU Jianhua,MA Fei,YAO Li.Hydraulic characteristics of slit-type energy dissipaters [J].Journal of Hydrodynamics,2012,24(6):883-887.

[2] 黃秋君,吳建華.收縮式消能工的研究現(xiàn)狀及進(jìn)展[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,2008,36(2):219-223.(HUANG Qiujun,WU Jianhua.Research on contraction energy dissipators [J].Journal of Hohai University:Natural Sciences,2008,36(2):219-223.(in Chinese)

[3] IPPEN A T.Gas-wave analogies in open channel flow [C]//Proceedings of the Second Hydraulics Conference.Iowa:University of Iowa,1943:248-265.

[4] HAGER W H.Supercritical flow in channel junctions [J].Journal of Hydraulic Engineering,1989,115(5):595-616.

[5] HAGER W H,SCHWALT M,JIMENEZ O,et al.Supercritical flow near an abrupt wall deflection [J].Journal of Hydraulic Research,1994,32(1):103-118.

[6] NI Hangen,LIU Yakun.A generalized explicit solution of the sequent depth ratio for the hydraulic jump [J].Journal of Hydrodynamics,2005,17(5):596-600.

[7] 劉亞坤.沖擊波與收縮式消能工若干問題的研究[D].大連:大連理工大學(xué),2006.

[8] 寧利中,戴振霖.窄縫挑坎收縮段急流沖擊波及控制水深的計(jì)算[J].西安理工大學(xué)學(xué)報(bào),1992,8(3):197-240,215.(NING Lizhong,DAI zhenlin.Calculations of shock waves and controlling water depth in region of contraction on slit-type flip bucket [J].Journal of Xi’an University of Technology,1992,8(3):197-240.(in Chinese)

[9] 黃智敏,翁情達(dá).窄縫挑坎收縮段急流沖擊波特性的探討[J].水利水電技術(shù),1989(8):11-15.(HUANG Zhimin,WENG Qingda.Discussions on characteristics of shock waves on slit-type flip bucket [J].Water Resources and Hydropower Engineering,1989(8):11-15.(in Chinese)

[10] 劉宣烈,張文周.空氣阻力對挑流水舌的影響[J].天津大學(xué)學(xué)報(bào),1982(2):67-77.(LIU Xuanlie,ZHANG Wenzhou.Effect of air resistance to a trajectory jet [J].Journal of Tianjin University,1982(2):67-77.(in Chinese)

[11] 劉宣烈,劉鈞,姚仲達(dá),等.空中摻氣水舌運(yùn)動(dòng)軌跡及射距[J].天津大學(xué)學(xué)報(bào),1989(2):23-30.(LIU Xuanlie,LIU Jun,YAO Zhongda,et al.The moving orbit and the horizontal length of aerated jet flow in open air [J].Journal of Tianjin University,1989(2):23-30.(in Chinese)

[12] 寧利中.挑流水舌挑距及其影響因素概述[J].水資源與水工程學(xué)報(bào),2004,15(3):35-39.(NING Lizhong.Discussion on the horizontal length of jet flow and its factors of influence [J].Journal of Water Resources & Water Engineering,2004,15(3):35-39.(in Chinese)

[13] 劉宣烈,劉鈞.三元空中水舌摻氣擴(kuò)散的試驗(yàn)研究[J].水利學(xué)報(bào),1989,17(11):10-17.(LIU Xuanlie,LIU Jun.Experimental study on the diffusion and aeration of three-dimensional jet [J].Journal of Hydraulic Engineering,1989,17(11):10-17.(in Chinese)

[14] 傅佩芬.溢流拱壩挑流水舌入水寬度計(jì)算[J].四川水力發(fā)電,2000,19(增刊1):70-71.(FU Peifen.The calculation of the jet width of overflow arch dams [J].Sichuan Hydropower,2000,19(Sup1):70-71.(in Chinese)

[15] 吳文平,張彥法.窄縫挑流水舌的運(yùn)動(dòng)擴(kuò)散規(guī)律及應(yīng)用[J].水力發(fā)電學(xué)報(bào),1989(4):71-76.(WU Wenping,ZHANG Yanfa.The theory and appliance of jet diffusion for slit-type flip bucket [J].Journal of Hydroelectric Engineering,1989(4):71-76.(in Chinese)

[16] 張葉林,鄧軍,劉軍,等.丹達(dá)河水電站窄縫消能研究[J].人民黃河,2011,33(6):126-131.(ZHANG Yelin,DENG Jun,LIU Jun,et al.Research on slit-type energy dissipaters of Dandahe Hydropower Station [J].Yellow River,2011,33(6):126-131.(in Chinese)