計算斜切挑坎沖坑深度的方法探討

曾 貝 佳

(重慶市水利電力建筑勘測設(shè)計研究院，重慶 渝北 401120)

?

計算斜切挑坎沖坑深度的方法探討

曾 貝 佳

(重慶市水利電力建筑勘測設(shè)計研究院，重慶 渝北 401120)

通過對斜切挑坎的數(shù)值模擬計算，將常規(guī)水力學(xué)參數(shù)計算結(jié)果與試驗結(jié)果進行對比分析后發(fā)現(xiàn)，數(shù)值計算在研究斜切坎的水深、壓強等方面可以得到與試驗研究相一致的結(jié)果，因此可以認(rèn)為，其挑流水舌入水長度的計算值也同樣與實驗結(jié)果具有良好的吻合度，從而避免了物理模型試驗由于測量水平的限制造成的較大的實驗誤差，進而可以根據(jù)該入水長度計算出水舌入水單寬流量并運用各家估算公式，得出沖坑平衡水墊深度。

斜切挑坎；數(shù)學(xué)模型；網(wǎng)格劃分；結(jié)果比較

0 引 言

斜切挑坎作為一種相對較為新型的挑坎形式，在工程界得到了廣泛的應(yīng)用。斜切挑坎是將傳統(tǒng)等寬挑坎沿某一角度斜切，從而使原來單一挑角的挑流鼻坎演變?yōu)橐幌盗杏胁煌羯浣堑亩嗵艚翘艨病Ｓ捎诟叩涂蔡艚堑牟町悾瑢?dǎo)致擴散坎內(nèi)水深明顯不同，出坎水深呈現(xiàn)出高坎大，低坎小的分布規(guī)律；另外，斜切挑坎不僅是一個簡單的挑流結(jié)構(gòu)，更重要的是一種散流結(jié)構(gòu)，盡管立面上斜切坎后的水流仍以一族不同的拋物線拋向下游河床，而平面上水舌卻發(fā)生了明顯的擴散，水舌有效入水長度大幅增加，入水單寬流量得到了極大程度地減小，從而明顯降低對下游河床的沖刷[1]。

國內(nèi)許多學(xué)者對斜切挑坎進行了物理模型試驗[2,3,4,5]以及數(shù)值模擬[6,7]的工作，得到了許多較好的成果。但是，對于斜切坎下游沖坑的計算，許多學(xué)者都提出了不同的估算公式，由于每個公式都包含有挑流水舌入水單寬流量，亦即需要事先得到水舌入水長度，而一般物理模型都是根據(jù)試驗測量得出，由于試驗手段限制的緣故，難以準(zhǔn)確獲得水舌入水長度，也就造成了相對較大的誤差。本文根據(jù)數(shù)值模擬計算得到水舌入水長度，并據(jù)此準(zhǔn)確算出水舌入水單寬流量，進而根據(jù)各家計算公式計算得出沖坑平衡水墊深度，從而使之與物理模型試驗進行對比，找出斜切挑坎下游沖坑平衡水墊深度的有效計算公式。

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

本文采用標(biāo)準(zhǔn)k-ε雙方程紊流模型，該紊流模型對雷諾應(yīng)力各分量采用各向同性假設(shè)。目前雙方程模型在實際工程和科研中使用較多，標(biāo)準(zhǔn)的k-ε模型比零方程和單方程模型有了很大的改進，在工程實際中得到了較為廣泛的應(yīng)用和驗證，它考慮了紊動速度比尺和紊動長度比尺的輸運。對于大多數(shù)水流問題，標(biāo)準(zhǔn)雙方程模型能得到較為滿意的結(jié)果。對于不可壓非定常流，標(biāo)準(zhǔn)k-ε紊流模型的連續(xù)方程、動量方程和k、ε方程分別為：

(1)

動量方程：

(2)

k方程：

(3)

ε方程：

(4)

式中 ρ和μ分別為體積分?jǐn)?shù)平均的密度和分子粘性系數(shù)。P為修正壓力；μt為紊流粘性系數(shù)，它可由紊動能k和紊動耗散率ε求出：

(5)

其中，Cμ為經(jīng)驗常數(shù)，Cμ取0.09。

σk和σε分別是k和ε的紊流普朗特數(shù)，σk=1.0，σε=1.3。C1ε和C2ε為ε方程常數(shù)，C1ε=1.44， C2ε=1.92。G為由平均流速梯度引起的紊動能產(chǎn)生項，它可以由下式定義：

(6)

引入水氣分層兩相流的VOF[7]模型后，在控制體內(nèi)對第q相流體的容積分?jǐn)?shù)規(guī)定為：αq=0表示控制體內(nèi)無q相流體；αq=1表示控制體內(nèi)充滿q相流體；0<αq<1表示控制體內(nèi)部分充滿q相流體，對所有流體的容積分?jǐn)?shù)總和為1，即：∑αq=1。在水氣分層兩相流中，αω為水的體積分?jǐn)?shù)，αa為氣的體積分?jǐn)?shù)，αω+αa=1，ρ和μ就是體積分?jǐn)?shù)的函數(shù)，而不是一個常數(shù)。它們可由下式表示：

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

(7)

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

(8)

式中 ρω和分別是水和氣的密度；μω和μa分別是水和氣的分子粘性系數(shù)，通過對水的體積分?jǐn)?shù)αω的迭代求解，ρ和μ都可以由式(7)、(8)求出。

2 計算體型及網(wǎng)格劃分

數(shù)學(xué)模型模擬范圍從挑坎起始端上游20 m位置開始，包括20 m長具有一定坡度的順直泄槽，斜切挑坎，以及挑流水舌在空中可能運動至的所有空間，即要求建立的計算域須要將所有空中水體質(zhì)點包含其內(nèi)。一般來講網(wǎng)格越細計算結(jié)果精度就越高，但過細的網(wǎng)格就意味著占用的內(nèi)存資源就越大，其所需計算時長也越長，由于本組計算模型計算區(qū)域相對較小，故而將網(wǎng)格劃分得比較細以保證足夠的計算精度，網(wǎng)格尺度的變化范圍0.08 m～0.2 m。

采用控制容積法[9]對偏微分方程組進行離散，壓力-速度耦合采用PISO算法。由于計算的下游出流位置流動已基本發(fā)展成穩(wěn)定狀態(tài)，可給定在出流邊界上的法向梯度為零，即：?v/?x= 0，?k/?x= 0，?ε/?x= 0；在固壁上給定法向的速度為零和無滑移條件，近壁的黏性底層采用壁函數(shù)法[9]處理。整個計算區(qū)域如圖1所示，水流入口采用速度入口邊界條件，根據(jù)模型試驗測得進口位置斷面水深h=3.55 m以及法向速度vn=34 m/s，出口及所有的氣體邊界均采用壓力邊界條件，壓力為大氣壓值。

圖1 數(shù)學(xué)模型計算域

3 計算結(jié)果與試驗結(jié)果的比較

3.1 計算結(jié)果

圖2是斜切挑坎水流流態(tài)圖，可以看到，數(shù)值計算流態(tài)與模型試驗流態(tài)非常吻合，出挑水舌空中形態(tài)良好，沿程連續(xù)平順，無不良形態(tài)出現(xiàn)。由于斜切的原因，水舌有向左側(cè)偏轉(zhuǎn)的趨勢，但是偏轉(zhuǎn)角度相對較小。同時，由于左側(cè)挑坎擴散的原因，水舌受到壓力梯度的影響，部分水體向左側(cè)偏轉(zhuǎn)，這有利于出挑水體的橫向擴散。通過后處理軟件得到水舌入水長度為39.66 m，而泄槽順直段寬度為10 m，這充分體現(xiàn)了斜切挑坎的對于水體的擴散作用，即由于高低坎的作用使得水流入水長度大幅增加。

圖2 流態(tài)圖

圖3為泄槽沿程水深計算值與試驗值的對比圖，可以看到，二者吻合較好。水流在前部順直渠道時沿程水深幾乎沒有變化，在進入反弧挑坎以后沿程水深逐漸減小，在出挑位置達到最小，水流由于斜切坎的存在在橫向上得到擴散，導(dǎo)致水體在立向上的水深逐漸減小。同時，沿程水深在進入挑坎前保持一致，進入挑坎以后由于左側(cè)墻逐漸向外側(cè)擴散以及左側(cè)挑角相對于右側(cè)為小的緣故，其左側(cè)水深低于右側(cè)水深。

圖3 泄槽沿程水深試驗值與計算值對比圖

圖4為泄槽中軸線沿程壓強的計算值與試驗值對比圖，可以看到，二者吻合良好。水流在前部順直渠道時沿程壓強幾乎沒有變化，在進入反弧挑坎以后呈現(xiàn)先增大后減小的趨勢，壓強在反弧段最低位置達到最大。水流進入反弧挑坎以后，受到挑坎施加的反向力作用，壓強增大，水深減小，直至挑坎最低位置壓強變到最大，其后由于挑坎橫向擴散的緣故其壓強又沿程逐漸減小。

圖4 泄槽底板中軸線沿程壓強分布

3.2 沖坑平衡水墊深度

從挑坎挑射出來的水體攜帶的巨大能量的耗散大體上分為三部分：即空中耗散消能、水墊強剪切與旋滾消能以及沖刷基巖消能。挑射水流在空中充分摻氣，實現(xiàn)了水流的擴散，事實上，其消耗能量的比例相對較小，它的最主要功能是為了通過自由擴散增大下游水面有效入水面積，減小入水單寬流量，為進入下游水體消能“做準(zhǔn)備”；跌入下游水墊后，入射水舌沿程擴散，斷面流速逐漸減小，主流與下游水墊之間產(chǎn)生強剪切作用并在水舌邊緣產(chǎn)生極為強烈的水躍旋滾以實現(xiàn)充分的能量耗散；然而，射流穿過水墊之后，仍然匯聚了相當(dāng)大的能量，有研究表明，其剩余能量相當(dāng)于總水頭的50%左右。若射流主要沖擊位置是可以被沖動的，那么河床將被逐漸掏刷以增加消能效率，直到?jīng)_刷坑內(nèi)部的水墊深度達到足以消除全部剩余能量位置，此時，沖刷坑里的最大水深被稱之為沖刷平衡水墊深度。

眾多的原型觀測以及理論研究表明，沖坑形成與發(fā)展的最大影響因素是挑射水體的單寬流量q，上下游水頭差Z以及基巖綜合抗沖系數(shù)Kr。

根據(jù)前文所述，斜切坎沿程水深以及壓強均與物理模型試驗具有較好的吻合度，故而可以認(rèn)為其出挑水舌形態(tài)以及入水長度也同樣具有較高的精度。根據(jù)前文得到的水舌入水長度為39.66 m，此時泄流總量Q=1 205.64 m3/s，從而計算得出水舌入水單寬流量q=30.40 m3/(m·s)。表1列舉了各家公式對于上下游水頭差Z=74 m時的下游河床沖刷平衡水墊深度計算結(jié)果。

表1 沖坑平衡水墊深度計算結(jié)果

王世夏公式中Kr=0.7～1.8，而余常昭公式中Kr=1.05～1.90，實驗過程中先對比了部分實驗結(jié)果與公式計算結(jié)果，以期得到一個固定的Kr值，從而實現(xiàn)減少實驗工作量的目的。通過對比分析，得出在本文設(shè)計挑坎情況下王世夏公式Kr=1.35時二者吻合較好，余常昭公式Kr=1.40，而其余公式則具有相對較大的誤差。物理模型試驗測得H0=18.34，可以看出，此時用王世夏公式以及余常昭公式均具有較好的吻合度，其余公式則出現(xiàn)相對較大的誤差。

4 結(jié) 語

通過對斜切挑坎的數(shù)值模擬計算，將常規(guī)水力學(xué)參數(shù)計算結(jié)果與試驗結(jié)果進行對比分析后發(fā)現(xiàn)，數(shù)值計算在研究斜切挑坎的水深、壓強等方面可以得到與試驗研究相一致的結(jié)果，因此可以認(rèn)為，其挑流水舌入水長度的計算值也同樣與實驗結(jié)果具有良好的吻合度，從而避免了物理模型試驗由于測量水平的限制造成的較大的實驗誤差，進而可以根據(jù)該入水長度計算出水舌入水單寬流量并運用各家估算公式，得出沖坑平衡水墊深度。通過對比分析，王世夏公式以及余常昭公式對于斜切挑坎下游沖坑水墊平衡深度有著較高的計算精度。

同時，根據(jù)本文所述方法通過數(shù)值模擬計算得出斜切挑坎入水長度并進一步算出入水單寬流量，以此數(shù)據(jù)帶入經(jīng)驗公式得到了與試驗數(shù)據(jù)吻合較好的下游沖坑深度，這從反向驗證了數(shù)值模擬入水長度結(jié)果的可靠性，并且表明利用數(shù)值模擬的入水長度進行下游沖坑深度計算的方式是可行的。

[1] 楊首龍.單圓弧斜切雙擴散坎消能工的水力設(shè)計[J].人民長江.2004,35(9):20-21.

[2] 張守磊,陳和春,張坤.斜切挑流鼻坎水舌挑距水力計算研究[J].中國水運.2008,8(12):160-161.

[3] 張東明,王媛.舌形坎與短邊墻坎泄流水舌空中擴散規(guī)律研究[J].中國農(nóng)村水利水電.2007,(3):106-109.

[4] 杜青，王曉燕，劉希成. 斜向挑坎在溢洪道設(shè)計中的應(yīng)用[J].山東水利.2004,(10):34-35.

[5] 莫海春，陳和春，羅倫，鄧偉.斜切挑坎水舌入水寬度隨切角變化規(guī)律研究[J].2013,31(11):124-126.

[6] 薛洪程，刁明軍，岳書波，徐蘭蘭.溢洪道出口斜切型挑坎挑流水舌三維數(shù)值模擬[J].2013,44(6):703-709.

[7] 陳大宏,李煒.自由表面流動數(shù)值模擬方法的探討[J].水動力學(xué)研究與進展: SerA,2001, 16(2): 216 -224.

[8] 葉茂，黃武林，伍超，陳開杰.導(dǎo)流洞挑流水舌的模型試驗及數(shù)值模擬研究[J].2011,30(4):98-102.

[9] Patankar S V. Numerical heat transfer and fluid flow [M].New York:Hemisphere Corporation and Mc Graw Hill Book Company,1980.

(責(zé)任編輯：卓政昌)

四川攀枝花1-8月太陽能發(fā)電增長227.8%

攀枝花市的清潔能源主要指水力發(fā)電、風(fēng)力發(fā)電和太陽能發(fā)電。1-8月，攀枝花市清潔能源發(fā)電117.91億千瓦小時，同比增長14.5%，占全市發(fā)電量的比重89%，比去年同期提升8.4個百分點。其中，水力發(fā)電量115.76億千瓦小時，同比增長12.5%；風(fēng)力發(fā)電1.81億千瓦小時，增長100%；太陽能發(fā)電0.34億千瓦小時，增長227.8%。

國電大渡河流域沙坪二級水電智慧工程鑒定成果達到國際領(lǐng)先水平

8月28日，國電大渡河公司《沙坪智慧工程管理模型研究與管控平臺開發(fā)應(yīng)用》項目科技鑒定會議在北京召開。鑒定會議認(rèn)為，該項目研究成果具備三項開拓性創(chuàng)新科技，整體達到了國際領(lǐng)先水平。一是提出了含有風(fēng)險(偏差)評價體系、事件預(yù)警體系、決策(糾偏)體系及考核評價體系的智慧工程管控理念，形成了覆蓋水電工程建設(shè)從宏觀整體到過程細節(jié)全方位、全過程的智慧化管理模型。二是基于集成了數(shù)據(jù)采集層(智能終端)、數(shù)據(jù)服務(wù)層(云數(shù)據(jù)中心)、智慧單元層(生產(chǎn)監(jiān)控)、智慧決策層(項目管控)的智慧工程一體化管控平臺總體架構(gòu)，利用統(tǒng)一通訊技術(shù)、BIM、云計算、物聯(lián)網(wǎng)、大數(shù)據(jù)等先進技術(shù)實現(xiàn)了混凝土生產(chǎn)監(jiān)控系統(tǒng)、3D數(shù)字廠房系統(tǒng)、施工安全監(jiān)測系統(tǒng)等多個智慧單元的應(yīng)用。三是提出了適用于水電流域梯級開發(fā)建設(shè)的智慧工程技術(shù)標(biāo)準(zhǔn)體系框架，編制了《智慧工程管控平臺技術(shù)導(dǎo)則》和《智慧工程三維信息模型技術(shù)標(biāo)準(zhǔn)》，推動智慧工程整體解決方案的應(yīng)用。

四川：5億補償火電 委屈了水電

9月19日，四川省發(fā)改委公示關(guān)于利用2015年水電站超發(fā)電量收益清算資金對公用燃煤火電機組欠發(fā)電量進行電價補償情況。據(jù)統(tǒng)計，2015年四川電網(wǎng)統(tǒng)調(diào)統(tǒng)分水電站超發(fā)電量約100億千瓦時，按超發(fā)電量的50%、每千瓦時0.10元進行清算，形成收益清算資金近5億元。 近年來，四川水電連年棄水，消納形勢越發(fā)嚴(yán)峻，加之電價持續(xù)下調(diào)，水電利潤收縮明顯，經(jīng)營壓力陡增。作為緩解四川部分水火矛盾的一項舉措，此次5億元補償火電機組，對當(dāng)?shù)厮娖髽I(yè)來說已經(jīng)有些力不從心。

2016-09-11

V211.76；O242.1；C32

B

1001-2184(2016)05-0104-04

曾貝佳(1990-)，女，苗族，重慶酉陽人，畢業(yè)于四川大學(xué)水工結(jié)構(gòu)專業(yè)，助理工程師，從事水庫、大壩、堤防、飲水管道設(shè)計工作.