陳 明，梁應(yīng)辰，宣國祥，陳明棟

(1.重慶交通大學(xué) 水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室，重慶 400074；2.中華人民共和國交通運(yùn)輸部，北京 100736； 3.南京水利科學(xué)研究院，江蘇 南京 210029)

?

船閘短廊道格柵消能室消能特性研究

陳 明1，梁應(yīng)辰2，宣國祥3，陳明棟1

(1.重慶交通大學(xué) 水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室，重慶 400074；2.中華人民共和國交通運(yùn)輸部，北京 100736； 3.南京水利科學(xué)研究院，江蘇 南京 210029)

針對帶格柵消能室的船閘短廊道頭部輸水系統(tǒng)，建立了整體輸水系統(tǒng)三維模型；通過RNGk～ε雙方程紊流模型，廊道工作閥門的開啟過程利用動網(wǎng)格技術(shù)，并結(jié)合自由水面處理技術(shù)的VOF法，對船閘正常運(yùn)行工況的灌水全過程進(jìn)行了三維流動動態(tài)仿真研究。以試驗(yàn)觀測得到的水位、流量和流速過程線，驗(yàn)證了數(shù)學(xué)模型和數(shù)值方法。根據(jù)消能率和閘室自由水面紊動能總量，分析了格柵消能室的消能效果；依據(jù)紊動能、紊動耗散率及流速分布揭示了消能機(jī)理。結(jié)果表明：格柵消能室可消耗90%以上的水流能量，且主要的能量消耗發(fā)生在消能室內(nèi)大量水流能量的交換過程中。

航道工程；船閘；輸水系統(tǒng)；消能；數(shù)值模擬

0 引 言

帶格柵消能室的短廊道輸水型式系一種典型的船閘集中輸水系統(tǒng)，應(yīng)用較為廣泛，我國的連江、高砂、丹竹、芒稻、興隆等船閘均采用了該類輸水型式。目前，對于集中輸水系統(tǒng)消能工消能特性的研究仍以物理模型試驗(yàn)為主。早在20世紀(jì)70年代，對于集中輸水系統(tǒng)消能工的設(shè)計(jì)和消能效果，荷蘭Delft水工實(shí)驗(yàn)室就進(jìn)行了相關(guān)的試驗(yàn)研究[1]。在我國，周華興[2]試驗(yàn)探討了船閘的垂直豎井輸水系統(tǒng)封閉式消能工的設(shè)計(jì)及其消能效果；連恒鐸[3]通過模型試驗(yàn)對一組不同消能設(shè)施的消能效果進(jìn)行了比較，結(jié)果表明采用倒口消能設(shè)施不但出流流態(tài)最平穩(wěn)而且獲得了基本消除險(xiǎn)象紊動的效果；鄒炳生等[4]根據(jù)短廊道輸水消能機(jī)理，結(jié)合工程經(jīng)驗(yàn)及輸水模型試驗(yàn)，探討了短廊道輸水雙層格柵消能室及其消能工設(shè)計(jì)中的問題；吳英卓等[5]針對興隆船閘采用的短廊道格柵消能室的輸水系統(tǒng)，通過物模試驗(yàn)分析了消能室出口處局部消力設(shè)施的消能效果。

采用模型試驗(yàn)方法，不僅耗時(shí)長，投入經(jīng)費(fèi)多，而且結(jié)果受試驗(yàn)條件及觀測技術(shù)影響，存在一定的尺度效應(yīng)和誤差。更重要的是難于得到輸水系統(tǒng)內(nèi)較為復(fù)雜水流結(jié)構(gòu)及與消能有關(guān)的水力要素(如紊動能，紊動耗散率等)。基于此，筆者針對某船閘帶格柵消能室的短廊道集中輸水系統(tǒng)，采用數(shù)值模擬的方法對船閘正常運(yùn)行工況的灌水全過程進(jìn)行三維水流動態(tài)仿真，研究了格柵消能室的消能效果并揭示了消能機(jī)理。

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

1.1 基本方程

船閘灌水過程的不可壓非恒定流是通過數(shù)值求解RNGk～ε紊流模型[6]而得，該模型的不可壓縮非定常流控制方程張量形式如下：

連續(xù)方程

(1)

動量方程

(2)

紊動動能k輸運(yùn)方程

(3)

紊動動能耗散ε輸運(yùn)方程

(4)

式(1)～式(4)采用控制體積法進(jìn)行離散，壓力、速度的耦合求解利用SIMPLEC算法。對無滑移的固壁邊界利用壁面函數(shù)法。閥門井及閘室自由液面采用水汽兩相流的VOF模型[7]進(jìn)行處理。

1.2 網(wǎng)格剖分及邊界條件

計(jì)算模型選取某船閘工程帶格柵消能室的短廊道頭部輸水系統(tǒng)(圖1)，閘室有效尺寸為140 m × 14 m × 2.5 m(長 × 寬 × 檻上最小水深)，上游最高和最低通航水位分別為45.00，40.50 m，下游最高和最低通航水位分別為43.37，36.00 m。因此其正常工作水頭H=9 m。格柵消能室正面布置12個格柵孔，頂面布置兩排格柵孔，共24個。

圖1 船閘輸水系統(tǒng)Fig.1 Shiplock filling and emptying system

該模型計(jì)算區(qū)域較大且復(fù)雜，包含了上引航道在內(nèi)共199 m的范圍，因此采用分塊四面體和六面體混合網(wǎng)格剖分計(jì)算區(qū)域，同時(shí)對復(fù)雜和重點(diǎn)研究部位(閥門廊道段和格柵式消能室及其附近區(qū)域)進(jìn)行了網(wǎng)格加密處理，剖分網(wǎng)格單元總數(shù)約920 000個，節(jié)點(diǎn)總數(shù)約360 000個。計(jì)算區(qū)域及工作閥門廊道段網(wǎng)格剖分如圖2。

圖2 計(jì)算區(qū)域及網(wǎng)格剖分Fig.2 Grid of computational area

鑒于上游引航道斷面平均流速較小，其產(chǎn)生的動水壓力值較小，所以進(jìn)口斷面可采用沿水深分布的靜水壓力進(jìn)口邊界；兩側(cè)閥門井頂部自由液面直接與大氣相通，采用空氣壓力進(jìn)口；閘室自由液面同樣與大氣相通，采用空氣壓力出口。設(shè)定的計(jì)算邊界條件如圖2所示。針對閥門開啟過程，為保證網(wǎng)格質(zhì)量，采用動網(wǎng)格技術(shù)調(diào)整和更新閥門運(yùn)動后的網(wǎng)格[8]。

根據(jù)閘室底部高程(32.50 m)和下游最低通航水位(36.00 m)可知，在T=0 s時(shí)，閘室初始水深為3.5 m。此外，閥門開啟方式與開啟速度對輸水流量及閘室水流條件影響較大。進(jìn)行數(shù)值計(jì)算時(shí)，閥門按7 min勻速開啟方式，速度為6.4 mm/s。

1.3 數(shù)學(xué)模型驗(yàn)證結(jié)果

采用船閘整體水工物理模型(比尺為λL=36)實(shí)測的灌水過程水力特性(包括水位、流量及流速過程線)對數(shù)學(xué)模型進(jìn)行驗(yàn)證。

1.3.1 水位、流量驗(yàn)證

圖3 水位、流量實(shí)測與計(jì)算值比較Fig.3 Comprison of calculated and experimental H and Q

圖3對比了灌水過程閘室中心點(diǎn)的水位及輸水流量計(jì)算值與實(shí)測值。由圖3可知，兩者灌水時(shí)間相差甚微，計(jì)算為535 s，實(shí)測為540 s，誤差僅為0.93%。水位變化過程線吻合較好，流量變化曲線與實(shí)測曲線基本一致。最大流量計(jì)算值為58.06 m3/s，實(shí)測值為60.58 m3/s，誤差為4.16%。此外，在420 s附近時(shí)段(即閥門開度約為0.79)，計(jì)算流量值相對偏大，原因在于當(dāng)工作閥門開度大于等于0.8后，原型流量系數(shù)大于模型流量系數(shù)所致[9]。

1.3.2 流速驗(yàn)證

在測定船閘灌水過程閘室的非恒定三維流速時(shí)，選取距閘室起始斷面57 m的橫斷面上不同水深處的任意3點(diǎn)，具體位置見表1。

表1 測點(diǎn)位置

閘室非恒定瞬態(tài)流速的采集采用三維點(diǎn)式聲學(xué)多普勒ADV流速儀進(jìn)行測量。根據(jù)主流方向(x方向)的流速計(jì)算值與實(shí)測值間的過程線對比(圖4)知，在整個灌水過程中，數(shù)值計(jì)算的非恒定流速總體變化趨勢與實(shí)測值基本吻合，較好地捕捉了閘室水流特性。

圖4 主流方向(x方向)流速過程線比較Fig.4 Comprison of calculated and experimental x-velocity of measure points

2 消能效果

2.1 消能率

以左右廊道出水口(1-1和2-2)為格柵消能室的入流進(jìn)口，以正面和頂面格柵孔(3-3和4-4)為格柵消能室的出流口(圖5)，通過水流能量可建立格柵消能室的消能率計(jì)算公式。

圖5 格柵消能室的消能率計(jì)算斷面位置Fig.5 Cross sections used for computation of energy dissipation rate

在格柵消能室段，水流能量主要以動能為主，因此消能率η可表示為：

(5)

式中：E1-1～E4-4分別為左右廊道出水口水流動能以及正面和頂面格柵孔出流剩余動能，均按式(6)進(jìn)行計(jì)算。

由文獻(xiàn)[10]可知，某一斷面上的水流動能E可按式(6)計(jì)算：

(6)

式中：Ac為某一研究斷面；v為水流流速大??；S為積分曲面。

根據(jù)式(5)可求出在整個灌水過程中關(guān)于格柵消能室的消能率，見圖6?？梢娤苈手稻?0%以上，說明該類集中輸水系統(tǒng)90%以上的水流能量主要倚靠格柵消能室予以消耗。

圖6 格柵消能室的消能率時(shí)程曲線Fig.6 The time-history curve of energy dissipation rate of room

2.2 閘室自由水面紊動情況

閘室水面的紊動或涌動情況是判別消能設(shè)施的消能效果的另一個重要方面。筆者以有無設(shè)置格柵消能室時(shí)的閘室自由表面紊動能進(jìn)行對比研究格柵消能室的消能效果。

若以閘室自由水面上的紊動能總量來分析整體紊動情況或水面涌動情況，可通過式(7)計(jì)算某一面上的紊動能總量：

(7)

式中：Af為所研究的閘室自由水面；k為閘室自由水面上任意處的紊動能。

對集中輸水系統(tǒng)，閘室自由水面紊動較為強(qiáng)烈的部位主要位于閘室的上半部分，筆者重點(diǎn)研究靠近格柵消能室附近的一段閘室區(qū)域(14.8 m≤x≤65 m)。

圖7為在有無格柵消能室情況下灌水期內(nèi)前340 s(包括出現(xiàn)最大流量的時(shí)刻在內(nèi))閘室自由表面紊動能總量。

圖7 有無格柵消能閘室自由表面紊動能總量的比較Fig.7 Comparison of total turbulent kinetic energy on free surface in lock chamber with and without energy dissipation room

從圖7中可看出，設(shè)有格柵消能室后的水面紊動明顯小于未設(shè)置格柵消能室的情形，且最大差值達(dá)55.83 m4/s2，說明格柵消能室的設(shè)置與否對該類輸水系統(tǒng)的設(shè)計(jì)和正常運(yùn)營起關(guān)鍵作用。同時(shí)還應(yīng)注意，盡管格柵消能室的消能效果相對于無消能設(shè)施的情況較為顯著，但由于消能室內(nèi)未布設(shè)任何型式的輔助消能工，在模型試驗(yàn)時(shí)對閘室流態(tài)的觀察發(fā)現(xiàn)，閘室水面的紊動對船舶的停泊條件依然造成一定影響。紊動能總量的最大值9.05 m4/s2亦可解釋實(shí)驗(yàn)觀測的現(xiàn)象，說明采用的格柵消能室尚存在較大的改進(jìn)余地。因此，為了達(dá)到更加優(yōu)越的消能效果，需對格柵消能的消能工布置型式或體型作進(jìn)一步的研究。

3 消能機(jī)理

揭示格柵消能室的消能機(jī)理，既可充分解釋其消能效果，又對其水力設(shè)計(jì)具有重要指導(dǎo)意義。實(shí)際上，格柵消能室消能過程即為水流能量的傳遞和轉(zhuǎn)換過程，而紊動能k和紊動耗散率ε是反映能量傳遞和轉(zhuǎn)換特性的重要參數(shù)。若紊動能和紊動耗散率的值越大，則能量傳遞和轉(zhuǎn)換速率也越快。

限于篇幅，圖8僅給出了在最大閘室斷面平均流速時(shí)水平剖面y=-5 m上的紊動能和紊動耗散率分布。從圖8中可知，紊動能和紊動耗散率主要分布在格柵消能室所處的位置(14.8 m≤x≤20.2 m)，由此也說明水流能量的傳遞和交換主要集中在格柵消能室內(nèi)發(fā)生。

圖8 閘室特征水平剖面紊動能、紊動耗散率分布Fig.8 Turbulent kinetic energy distribution and Turbulent dissipation rate distribution on horizontal section in lock chamber

格柵消能室通過水流能量的傳遞和交換進(jìn)行消能是左右兩側(cè)廊道內(nèi)的兩股高速水流進(jìn)入消能室后發(fā)生強(qiáng)烈碰撞和對沖，大部分水體參與并形成兩個較為強(qiáng)烈的反向漩渦，在左側(cè)形成順時(shí)針漩渦，在右側(cè)形成逆時(shí)針漩渦，如流速矢量圖9，與此同時(shí)，水流與消能室壁面發(fā)生劇烈摩擦。在以上多種運(yùn)動方式的助推下發(fā)生水流剪切、摩擦、摻混等現(xiàn)象。

圖9 閘室特征水平剖面流速矢量(T=209 s)

由于黏性作用，由時(shí)均動能轉(zhuǎn)化來的脈動動能變?yōu)闊崮芎纳ⅲ鲹p失了能量，從而引起水流動能向熱能的交換，以取得能量消剎的效果。事實(shí)上，格柵消能室的消能過程可分為三步：

1)在高速水流由廊道流入消能室時(shí)，由于斷面的突擴(kuò)，形成了第1次水流能量擴(kuò)散過程；

2)在消能室內(nèi)發(fā)生大量水流能量交換過程；

3)當(dāng)水流流經(jīng)格柵孔時(shí)，由于斷面的突縮和突擴(kuò)，致使水流能量得到進(jìn)一步消耗。

4 結(jié) 論

船閘帶有格柵消能室的短廊道頭部輸水系統(tǒng)，其輸水過程產(chǎn)生的水流屬非定常流，流態(tài)相當(dāng)復(fù)雜。筆者采用三維數(shù)值模擬研究手段，得到流速分布、紊動能和紊動耗散率分布等詳細(xì)的水力特性，研究得到以下結(jié)論：

1)數(shù)值模擬結(jié)果與物理模型實(shí)測數(shù)據(jù)對比結(jié)果表明：兩者吻合較好。驗(yàn)證了數(shù)值計(jì)算結(jié)果的合理性和可靠性，說明利用文中提出的數(shù)值計(jì)算方法進(jìn)行船閘輸水全過程三維流動的動態(tài)仿真研究是可行的。

2)通過對格柵消能室的消能效果分析得知，格柵消能室的消能率較高，可消耗90%以上的水流能量，但由于消能室內(nèi)未布設(shè)任何型式的輔助消能工，閘室內(nèi)仍存在一定程度的水面紊動。因此采用何種消能工布置型式以穩(wěn)定閘室自由水面需作進(jìn)一步的研究。

3)在對消能機(jī)理的研究中，提出格柵消能室的消能量主要包括高速水流由廊道流入消能室時(shí)消耗的能量、在消能室內(nèi)發(fā)生大量水流能量交換過程消耗的能量，以及水流流經(jīng)格柵孔時(shí)消耗的能量3部分。而且，主要的水流能量消剎是發(fā)生在消能室內(nèi)進(jìn)行大量水流能量的交換過程中。

[1] Delft Hydraulics Laboratory.Low-Head Navigation Lock Door Filling and Emptying Systems Developed Hydraulic Investigation [R].Delft:Delft Hydraulics Laboratory,1973.

[2] 周華興.船閘豎井短廊道輸水系統(tǒng)的合理布置[J].水道港口,1980 (4):38-43. Zhou Huaxing.Reasonable layout type of shiplock filling and emptying system with short culverts and vertical shaft [J].Journal of Waterway and Harbor,1980(4):38-43.

[3] 連恒鐸.倒口出流消能技術(shù)在船閘集中輸水系統(tǒng)設(shè)計(jì)中的應(yīng)用[J].水運(yùn)工程,1996(4):15-18. Lian Hengduo.Stilling technique of inverte doutlet flow in design of concentrat ed filling & emptying system of shiolock [J].Prot & Waterway Engineering,1996(4):15-18.

[4] 鄒炳生,王召兵.雙層格柵消能室在中水頭短廊道輸水船閘中的應(yīng)用[J].水運(yùn)工程,2002(8):54-57,62. Zou Bingsheng,Wang Zhaobing.Application of double-deck grid stilling chamber in middle water head short-culvert conveyance shiplock [J].Prot & Waterway Engineering,2002(8):54-57,62.

[5] 吳英卓,江耀祖,王智娟,等.興隆船閘關(guān)鍵水力學(xué)問題研究[J].水利水運(yùn)工程學(xué)報(bào),2012(4):6-12. Wu Yingzhuo,Jiang Yaozu,Wang Zhijuan,et al.Research of key hydraulic problems of Xinglong Shiplock[J].Hydro-Science and Engineering,2012(4):6-12.

[6] Yakhot V,Orzag S A.Renormalizationn group analysis of turbulence:basic theory [J].Journal of Scientific Computation,1986,1:3-11.

[7] Hirt C W,Nichols B D.Volume of fluid (VOF) method for the dynamics of free boundary [J].Journal of Computational Physics,1981,39:201-225.

[8] 楊忠超,楊斌,陳明棟,等.高水頭船閘閥門開啟過程水力特性仿真研究[J].重慶交通大學(xué)學(xué)報(bào):自然科學(xué)版,2009,28(4):758-762. Yang Zhongchao,Yang Bin,Chen Mingdong,et al.Simulation study on hydraulics characteristic of high head lock during valve opening process based on dynamic mesh technology [J].Journal of Chongqing Jiaotong University:Natural Science,2009,28(4):758- 762.

[9] 須清華.船閘輸水系統(tǒng)的水工模型縮尺影響[J].水利水運(yùn)科學(xué)研究,1985(4):1-8. Xu Qinghua.Scale influence of hydraulic model test for lock filling and emptying system [J].Hydro-Science and Engineering,1985(4):1 -8.

[10] 王作高.船閘設(shè)計(jì)[M].北京:水利電力出版社,1992. Wang Zuogao.Design of Shiplock [M].Beijing:Water Resources and Electric Power Press,1992.

Energy Dissipation Characteristics of Room with Grating Holes of Shiplock Filling and Emptying System

Chen Ming1，Liang Yingchen2，Xuan Guoxiang3，Chen Mingdong1

(1. Key Laboratory of Hydraulic & Waterway Engineering of the Ministry of Education of the People’s Republic of China, Chongqing Jiaotong University, Chongqing 400074, China; 2. Ministry of Transport of the People’s Republic of China, Beijing, 100736, China; 3 Nanjing Hydraulic Research Institute, Nanjing 210029, Jiangsu, China)

As for the lock-head filling system of shiplock composed of short culverts and energy dissipation room with grating holes, its whole 3-D model was established. Through RNGk～εmodel of turbulence flow, the 3-D numerical simulation for dynamic characteristics of filling flow under the normal operation condition were conducted by using dynamic mesh technology for the valve and method. According to the energy dissipation ratio and the amount of free surface turbulence kinetic energy,the energy dissipation effect of room with grating holes was analyzed. According to the turbulence kinetic energy, turbulence energy dissipation ratio and velocity distribution, the energy dissipation mechanism was revealed. The results indicated that more than 90% of thewater flow energy can be consumed by the energy dissipation room with grating holes and the energy consumption mainly occurs in the process of energy exchange in the energy dissipation room.

waterway engineering; shiplock; filling and emptying system; energy dissipation; numerical simulation

10.3969/j.issn.1674-0696.2015.05.17

2013-05-06；

2015-06-10

重慶市教委科學(xué)技術(shù)研究項(xiàng)目(KJ1400322)；重慶市科委自然科學(xué)基金項(xiàng)目(CSTC2011JJA30002)

陳 明(1983—)，男，江西九江人，博士，主要從事通航水力學(xué)方面的研究。E-mail：chenmingjy@126.com。

U641.1

A

1674-0696(2015)05-084-05