傅銘煥，張志昌

(1.浙江省水利水電勘測設(shè)計院，浙江杭州 310002； 2.西安理工大學(xué)，陜西西安 710048)

?

密排加糙床面消力池自由水躍躍長計算

傅銘煥1,2，張志昌2

(1.浙江省水利水電勘測設(shè)計院，浙江杭州 310002； 2.西安理工大學(xué)，陜西西安 710048)

水躍長度作為消力池設(shè)計的重要參數(shù)，對消力池安全穩(wěn)定以及經(jīng)濟(jì)合理的影響效果顯著。通過建立水躍區(qū)水體質(zhì)點的運動方程，研究密排加糙床面消力池水躍旋滾長度和水躍長度的變化規(guī)律，提出了水躍旋滾長度和水躍長度計算的理論方法，并根據(jù)已有文獻(xiàn)的試驗數(shù)據(jù)對所推公式涉及的物理參數(shù)進(jìn)行率定。計算研究結(jié)果表明，密排加糙床面消力池水躍旋滾長度和水躍長度均隨躍前斷面弗勞德數(shù)、躍前斷面水深和水躍共軛水深比的增大而增大，隨著床面當(dāng)量粗糙高度的增加而減小。經(jīng)驗證發(fā)現(xiàn)，水躍旋滾長度和水躍長度的計算值與實測值接近。

密排加糙床面； 運動方程； 水躍旋滾長度； 水躍長度

在泄水建筑物下游往往需要布設(shè)消力池來抬高下游水深，使下泄水流形成淹沒水躍，消耗高速水流能量。為了提高消力池消能效果，各種輔助消能工應(yīng)運而生，如消力墩、T形墩、齒墩等。近年來有人提出計算了粗糙壁面消力池，并對其水力特性進(jìn)行了試驗研究。研究表明，消力池底板上的連續(xù)凸體能顯著提升水躍區(qū)水體的剪切與紊動強(qiáng)度，有效降低消力池的共軛水深和水躍長度，從而獲得更好的消能效果[1- 4]。所謂的粗糙壁面消力池，就是將消力池底板做成連續(xù)的粗糙塊，現(xiàn)階段粗糙壁面消力池的布置形式主要有兩種，一種為波狀床面消力池，另一種為密排加糙床面消力池。文獻(xiàn)[5-6]通過邊界層理論和水躍區(qū)動量方程研究了波狀床面消力池壁面切應(yīng)力和水躍區(qū)邊界層沿程發(fā)展規(guī)律，給出了水躍共軛水深和水躍長度的理論計算方法。文獻(xiàn)[7]通過前人關(guān)于床面阻力的動量假設(shè)，提出了密排加糙床面水躍共軛水深計算的半理論公式。文獻(xiàn)[8]對文獻(xiàn)[1]和[3]作者提出的密排加糙床面水躍長度和水躍旋滾長度經(jīng)驗公式進(jìn)行了驗證，并提出了計算水躍長度和水躍旋滾長度新的經(jīng)驗公式。文獻(xiàn)[9]研究了均布糙體和孤立糙體的空化特性，結(jié)果表明，均布糙體的初生空化數(shù)遠(yuǎn)低于孤立糙體；文獻(xiàn)[10]在對明流反弧段加糙防蝕的研究中，認(rèn)為在反弧段進(jìn)行均勻加糙，對防止反弧末端的空蝕破壞十分有利。文獻(xiàn)[11]對波狀涵洞中流速的研究表明，波狀床面水流紊動加劇，雷諾切應(yīng)力增大，流速減小明顯，空蝕的可能性亦減小。由此可見，在溢流面上適當(dāng)?shù)木鶆蚣硬冢瑢Ψ乐箍瘴g破壞是有利的。但現(xiàn)階段對其水躍水力特性，尤其是水躍長度的研究仍停留在試驗層面。文獻(xiàn)[12]和[13]通過水躍區(qū)水體質(zhì)點的運動方程，分別從理論角度研究了平底二元水躍和漸擴(kuò)散水躍的水躍長度，獲得了很好的效果。作者亦試圖通過追蹤水躍區(qū)水體質(zhì)點的運動規(guī)律，探索密排加糙床面自由水躍長度理論計算的可行性，以供工程設(shè)計參考。

1 公式推導(dǎo)

1.1 水躍旋滾長度

圖1 密排加糙床面水躍Fig.1 Hydraulic jump on uniform rough height bed

1.1.1 理論分析 水流從急流到緩流的過渡過程中，上下游水流存在壓力差使得水體向上躍起形成回流，假設(shè)使水流向上躍起的加速度為f，對于密排加糙床面消力池，假定由于床面底部的凸起，額外作用于水流質(zhì)點豎直方向上的加速度為Mf，則密排加糙床面自由水躍水體質(zhì)點的受力情況如圖1所示。圖1中h1，h2分別為躍前斷面和躍后斷面水深；v1，v2分別為躍前斷面和躍后斷面平均流速；Lr為水躍旋滾長度；g為重力加速度；ks為密排加糙床面粗糙高度；M為由于床面的密排加糙額外在豎直方向上產(chǎn)生的加速度系數(shù)；x，y為坐標(biāo)軸，x軸方向同水體主流方向，y軸方向為豎直方向，坐標(biāo)原點為水躍躍前斷面水面處。則水躍區(qū)水體質(zhì)點在豎直方向y上沿時間的運動方程為：

d2y/dt2=(1+M)f-g

(1)

式中：t為水躍區(qū)水體質(zhì)點從躍前斷面開始運動的時間；y為水體質(zhì)點在t時間內(nèi)在豎直方向上運動的距離。

對式(1)求不定積分可得：

(2)

(3)

式中：c1，c2為積分常數(shù)。

假設(shè)在t=0時刻，水體質(zhì)點位于坐標(biāo)軸原點，此時質(zhì)點流速同主流方向，即水體質(zhì)點在豎直方向上的瞬時速度為零，可知：

(4)

將式(4)代入式(2)和(3)，可得c1=0，c2=0，故

(5)

根據(jù)總流方程，假定Q=v1bh1=vxb(h1+y)。其中vx為水體質(zhì)點流經(jīng)x斷面處的斷面平均流速；b為消力池寬度。假定水體質(zhì)點運動至x斷面時在x軸方向上的瞬時速度正比于該斷面的平均流速，可得

dx/dt=avx=av1/(1+y/h1)

(6)

式中：a為流速系數(shù)。

將式(5)代入式(6)，并對其變形可得

(7)

假定在t=0到t=tr時段內(nèi)，水體質(zhì)點從坐標(biāo)原點運動至水躍旋滾末端水面處，則對式(7)進(jìn)行積分，即

(8)

(9)

根據(jù)假定，t=tr時，水體質(zhì)點從原點運動至水躍旋滾末端水面，此時y=h2-h1。將其代入式(5)可得：

(10)

將式(10)代入式(9)，即：

(11)

(12)

當(dāng)密排加糙床面轉(zhuǎn)化成一般混凝土壁面時，Mf=0，此時式(12)變?yōu)橐话慊炷帘诿嫦Τ厮S旋滾長度，即：

(13)

文獻(xiàn)[14]分析認(rèn)為，在躍前斷面弗勞德數(shù)Fr1=1.82～7.67范圍內(nèi)，一般混凝土壁面水躍長度是水躍旋滾長度的1.973～1.614倍。取其均值，故一般混凝土壁面水躍長度可表示為：

(14)

文獻(xiàn)[15]提出的一般混凝土壁面水躍長度的計算式為：

(15)

聯(lián)立求解式(13)，(14)和(15)，可得：

(16)

將式(16)代入式(12)，得密排加糙床面水躍旋滾長度理論計算式為：

(17)

圖2 參數(shù)D的分布Fig.2 Distribution of parameter D

文獻(xiàn)[3]的作者進(jìn)行了5種體形的密排加糙床面消力池水躍特性試驗，粗糙高度ks分別為0.46，0.82，1.46，2.39和3.20 cm。根據(jù)文獻(xiàn)[3]的試驗工況，對密排加糙床面水躍區(qū)水體質(zhì)點相對加速度參數(shù)D重新分析，結(jié)果如圖2所示，圖中ds=2ks/(h1+h2)為消力池底板當(dāng)量粗糙高度。由圖2可知，相對加速度參數(shù)D是躍前斷面弗勞德數(shù)，水躍共軛水深比和當(dāng)量粗糙高度的函數(shù)，并隨著弗勞德數(shù)和當(dāng)量粗糙高度的增加而減小，隨著水躍共軛水深比的增大而增大。擬合各工況的試驗數(shù)據(jù)可得相對加速度參數(shù)D的計算式，即：

D=2.087 6(Fr1-1)-0.109 8[ds+2(η-1)/η]-1.109 8

(18)

將式(18)代入式(17)即可得密排加糙床面自由水躍旋滾長度半理論計算式，即：

(19)

式(19)的適用范圍為2.1

由式(19)可知，水躍旋滾長度隨著躍前斷面弗勞德數(shù)、躍前斷面水深和水躍共軛水深比的增大而增大，隨著床面當(dāng)量粗糙高度的增加而減小。

1.1.2 公式簡化 由于密排加糙床面水躍的躍后條件往往是未知的，故用式(19)來求解水躍旋滾長度存在一定困難。文獻(xiàn)[8]提出的密排加糙床面水躍共軛水深為：

(20)

將式(20)代入式(19)可得：

(21)

可見，用式(21)計算水躍旋滾長度較為復(fù)雜，根據(jù)式(21)中水躍旋滾長度的變化規(guī)律，將其轉(zhuǎn)化為：

Lr=7.62Fr1h1E/(Fr1-1)0.109 8

(22)

圖3 參數(shù)E的分布Fig.3 Distribution of parameter E

式中：E為床面粗糙高度、躍前斷面水深和躍前斷面弗勞德數(shù)的函數(shù)。

根據(jù)文獻(xiàn)[3]的試驗工況，分析得E隨床面粗糙高度、躍前斷面水深和躍前斷面弗勞德數(shù)的變化規(guī)律如圖3所示。由圖3可得

(23)

將式(23)代入式(22)，可得水躍旋滾長度的簡化公式，即：

(24)

式(24)的適用范圍為2.10

1.2 水躍長度

1.2.1 理論分析 對于密排加糙床面，文獻(xiàn)[8]分析發(fā)現(xiàn)，當(dāng)弗勞德數(shù)Fr1=2.4～10.5時，水躍長度為水躍旋滾長度的1.30～1.46倍，取其均值1.38，則水躍長度可表示為：

Lj=1.38Lr+erf1(Lr)

(25)

式中：erf1(Lr)為密排加糙床面水躍長度的偏離誤差。

將式(25)代入式(17)，可得：

(26)

圖4 參數(shù)A的分布Fig.4 Distribution of parameter A

式中：A為考慮水躍區(qū)相對加速度和水躍長度偏離誤差的參數(shù)。式(26)即為密排加糙床面水躍長度的理論公式。

文獻(xiàn)[1]研究了3種體形的密排加糙床面消力池水躍特性，粗糙高度ks分別為0.43，0.64和1.13 cm。根據(jù)文獻(xiàn)[1]的試驗工況，計算參數(shù)A如圖4所示。由圖4可知，參數(shù)A的分布規(guī)律與相對加速度參數(shù)D相似，仍是躍前斷面弗勞德數(shù)，水躍共軛水深比和當(dāng)量粗糙高度的函數(shù)。對各工況的試驗數(shù)據(jù)進(jìn)行擬合，可得

(27)

將式(27)代入式(26)，可得水躍長度的半理論計算式為：

(28)

式(28)的適用范圍為2.4

可見，密排加糙床面水躍旋滾長度隨著躍前斷面弗勞德數(shù)、躍前斷面水深和水躍共軛水深比的增大而增大，隨著床面當(dāng)量粗糙高度的增加而減小。

1.2.2 公式簡化 將式(20)代入式(28)，并根據(jù)式(28)表明的函數(shù)關(guān)系，對其重新整理可得：

(29)

圖5 參數(shù)F的分布Fig.5 Distribution of parameter F

式中：F為床面粗糙高度、躍前斷面水深和躍前斷面弗勞德數(shù)的函數(shù)。

根據(jù)文獻(xiàn)[3]的工況，對F隨床面粗糙高度、躍前斷面水深和躍前斷面弗勞德數(shù)的的變化規(guī)律進(jìn)行分析，結(jié)果如圖5所示，可得：

(30)

將式(30)代入式(29)，則水躍長度計算式可簡化為

(31)

式(31)的適用范圍為2.4

2 公式分析與驗證

圖6給出了文獻(xiàn)[3]實測的水躍旋滾長度試驗值與式(19)和式(24)計算的水躍旋滾長度計算值，圖中弗勞德數(shù)范圍為2.10～9.89。從圖6可以看出，式(19)和式(24)計算的水躍旋滾長度均與實測值接近。在277組試驗工況中，式(19)計算的水躍旋滾長度平均誤差為6.60%，誤差超過10%的共有57組，其中誤差15%～20%的有12組，20%～30%的有9組，超過30%的有3組。式(24)計算的水躍旋滾長度平均誤差為7.4%，誤差超過10%的共有66組，其中誤差15%～20%的有12組，20%～30%的有8組，超過30%的有3組。

圖6 旋滾長度的比較Fig.6 Comparison of vortex roll lengths

文獻(xiàn)[3]亦給出過考慮壁面粗糙高度的水躍旋滾長度經(jīng)驗計算式，即：

Lr=ch1(Fr1-1)

(32)

式中：c=6.525exp(-0.6ks/h1)。

式(32)的水躍旋滾長度平均誤差為15%，誤差大于10%的有179組，其中10%～20%的有110組，20%～30%的有45組，30%～40%的有15組，40%～50%的有7組，誤差大于50%的有2組[8]。

文獻(xiàn)[8]對文獻(xiàn)[3]的試驗成果重新進(jìn)行優(yōu)化分析，給出的水躍旋滾長度經(jīng)驗計算式為：

(33)

式(33)的平均誤差4.25%，超過10%的54組，其中有20組超過了15%[8]。由此可知，式(19)和式(24)計算的水躍旋滾長度精度遠(yuǎn)高于式(32)的計算值，而與式(33)的計算精度較為接近。

圖7給出了文獻(xiàn)[1]實測的水躍長度試驗值與式(28)和式(31)計算的水躍長度計算值，圖中弗勞德數(shù)范圍為2.4～10.5。由圖7可知，式(28)和(31)的計算值均與實測值接近，式(28)計算的平均誤差為7.46%。在92組試驗工況中，誤差超過10%的共有22組，其中10%～15%的有12組，超過20%的有4組，最大誤差為28.1%。式(31)計算的平均誤差為7.75%，誤差超過10%的共有30組，其中24組數(shù)據(jù)誤差均小于15%，最大誤差為23.7%。

圖7 水躍長度的比較Fig.7 Comparison of hydraulic jump lengths

可見，本文提出的水躍旋滾長度與水躍長度的半理論公式和簡化公式精度均較高，能很好反映水躍旋滾長度和水躍長度的實際情況，表明用水躍區(qū)水體質(zhì)點的運動方程研究水躍旋滾長度和水躍長度是可行的。分析認(rèn)為因水躍區(qū)強(qiáng)烈紊動，使得水躍旋滾末端前后波動明顯，施測困難造成公式計算值的局部偏離。

3 結(jié) 語

通過水躍區(qū)水體質(zhì)點的運動方程，理論研究密排加糙床面水躍旋滾長度和水躍長度的變化規(guī)律。研究表明，水躍區(qū)相對加速度參數(shù)D和A是躍前斷面弗勞德數(shù)，水躍共軛水深比和當(dāng)量粗糙高度的函數(shù)；水躍旋滾長度和水躍長度隨著躍前斷面弗勞德數(shù)、躍前斷面水深和水躍共軛水深比的增大而增大，隨著床面當(dāng)量粗糙高度的增加而減小。水躍旋滾長度半理論計算式(19)計算的水躍旋滾長度在弗勞德數(shù)2.1～9.89范圍內(nèi)，平均誤差為6.60%；簡化計算式(24)計算的平均誤差為7.4%。水躍長度半理論計算式(28)計算的水躍長度在弗勞德數(shù)2.4～10.5范圍內(nèi)，平均誤差為7.46%；簡化計算式(31)計算的平均誤差為7.75%。

[1]HUGHES W C,ERNEST FLACK J.Hydraulic jump properties over a rough bed[J].Journal of Hydraulic Engineering,1984,110(12)：1755- 1771.

[2]EAD S A,RAJARATNAM N.Hydraulic jumps on corrugated beds[J].Journal of Hydraulic Engineering,2002,128(7)：656- 663.

[3]CAROLLO F G,FERRO V,PAMPALONE V.Hydraulic jumps on rough beds[J].Journal of Hydraulic Engineering,2007,133(9)：989- 999.

[4]ABBASPOUR A,DALIR A H,FARSADIZADEH D,et al.Effect of sinusoidal corrugated bed on hydraulic jump characteristics[J].Journal of Hydro-environmental Research,2009(3)：109- 117.

[5]張志昌,傅銘煥,李若冰.波狀床面消力池共軛水深和水躍長度的計算[J].水力發(fā)電學(xué)報,2014,33(5)：120- 127.(ZHANG Zhi-chang,FU Ming-huan,LI Ruo-bing.Calculation of conjugate depths and length of jump on corrugated beds of stilling basin[J].Journal of Hydroelectric Engineering,2014,33(5)：120- 127.(in Chinese))

[6]張志昌,傅銘煥,趙瑩,等.波浪形底板消力池自由水躍特性的探討[J].應(yīng)用力學(xué)學(xué)報,2013,30(6)：870- 875.(ZHANG Zhi-chang,FU Ming-huan,ZHAO Ying,et al.Study on properties of free hydraulic jump on corrugated beds of stilling basin[J].Chinese Journal of Applied Mechanics,2013,30(6)：870- 875.(in Chinese))

[7]張志昌，傅銘煥，李若冰,等.加糙消力池共軛水深和水躍長度的試驗分析[J].水利水運工程學(xué)報，2013(6)：61- 66.(ZHANG Zhi-chang,FU Ming-huan,LI Ruo-bing,et al.Test analysis of conjugate depth and hydraulic jump length of roughened stilling basin[J].Hydro-Science and Engineering，2013(6)：61- 66.(in Chinese))

[8]張志昌,傅銘煥,趙瑩,等.人工粗糙壁面的水躍特性研究[J].應(yīng)用力學(xué)學(xué)報,2014，31(2)：270- 275.(ZHANG Zhi-chang,FU Ming-huan,ZHAO Ying,et al.Characteristics of hydraulic jump on artificial roughened beds[J].Chinese Journal of Applied Mechanics，2014，31(2)：270- 275.(in Chinese))

[9]ARNDT R A,IPPEN A T.Rough surface effects on cavitation inception[J].Journal Basic Engneering，1968，90(2)：249- 261.

[10]聶孟喜.明流反弧段加糙防蝕的初步試驗研究[J].水力發(fā)電學(xué)報，1987，6(4)：79- 88.(NIE Meng-xi.Preliminary experimental study of roughening bucket for cavitation protection in open flow[J].Journal of Hydroelectric Engineering，1987，6(4)：79- 88.(in Chinese))

[11]EAD S A,RAJARATNAM N.Plane turbulent wall jets on rough boundaries with limited tailwater[J].Journal of Engineering Mechanics，2004，130(10)：1245- 1250.

[12]詹明安，周建國，黃文政.關(guān)于水躍躍長的計算[J].人民長江，1990(3)：18- 21.(ZHAN Ming-an,ZHOU Jian-guo,HUANG Wen-zheng.Calculation of length of hydraulic jump[J].Yangtze River，1990(3)：18- 21.(in Chinese))

[13]吳宇峰，伍超，劉小兵.漸擴(kuò)散水躍躍長的研究[J].水科學(xué)進(jìn)展，2007，18(2)：210- 215.(WU Yu-feng,WU Chao,LIU Xiao-bing.Length of hydraulic jump in gradually broadened channel[J].Advances in Water Science，2007，18(2)：210- 215.(in Chinese))

[14]張志昌，趙瑩，傅銘煥.矩形平底明渠水躍長度公式的分析與應(yīng)用[J].西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版)，2014，42(11)：190- 198.(ZHANG Zhi-chang,ZHAO Ying,FU Ming-huan.Analysis and application of rectangular flat-bottomed open channel hydraulic jump length formulas[J].Journal of Northwest A & F University(Nat Sci Ed)，2014，42(11)：190- 198.(in Chinese))

[15]劉沛清.矩形明渠水躍段速度分布的理論分析[J].水利學(xué)報,1993(9)：48- 54.(LIU Pei-qing.Theoretical analysis on rectangular open channel velocity distribution in hydraulic jump region[J].Journal of Hydraulic Engineering,1993(9)：48- 54.(in Chinese)

Length calculation of free hydraulic jump on rough bed of stilling basin

FU Ming-huan1,2，ZHANG Zhi-chang2

(1.ZhejiangDesignInstituteofWaterConservancyandHydroelectricPower,Hangzhou310002,China; 2.Xi’anUniversityofTechnology,Xi’an710048,China)

As the important parameters in designing the stilling basin,calculation of the length of the free hydraulic jump plays a major role in the safety,stability and economization for the stilling basin.The research on the variation laws of the vortex roll length of the free hydraulic jump and the hydraulic jump length occurred on the uniform rough height bed of the stilling basin was carried out according to establishing the motion equations for water body particle in the hydraulic jump zone,and developing a theoretical method for calculation of the vortex roll length of the free hydraulic jump and the hydraulic jump length,and the physical parameters are calibrated by the given expressions based on the testing data from the existing literatures.Research results show that both of the hydraulic jump vortex roll length and the hydraulic jump length on the uniform rough height bed of the stilling basin increase with the increase of the upstream Froude number,upstream supercritical water depth and the conjugate depth ratio,and at the same time decrease with the increase of the equivalent roughness of the riverbed.From the verification results it is found that the calculated values are well approximate to the measured values of the hydraulic jump vortex roll length and the hydraulic jump lengths.

uniform rough height beds; motion equations; hydraulic jump vortex roll length; length of hydraulic jump

10.16198/j.cnki.1009-640X.2016.06.010

傅銘煥，張志昌.密排加糙床面消力池自由水躍躍長計算[J].水利水運工程學(xué)報,2016(6):69-75.(FU Ming-huan,ZHANG Zhi-chang.Length calculation of free hydraulic jump on rough bed of stilling basin[J].Hydro-Science and Engineering,2016(6):69-75.)

2015-12-28

傅銘煥(1989—)，男，浙江杭州人，碩士，主要從事水工水力學(xué)計算工作。E-mail:fuminghuan2007@163.com

TV135.2+1

A

1009-640X(2016)06-0069-07