呂偉東 ，劉平昌 ，陳 亮 ,張緒進

（1.重慶交通大學,重慶 400074；2.廣西西江開發(fā)投資集團有限公司，南寧 530022；3.重慶西南水運工程科學研究所,重慶 400016）

船閘閘室明溝消能工消能效果三維數(shù)值模擬

呂偉東1，2，劉平昌1，3，陳 亮1,張緒進1，3

（1.重慶交通大學,重慶 400074；2.廣西西江開發(fā)投資集團有限公司，南寧 530022；3.重慶西南水運工程科學研究所,重慶 400016）

利用多重網(wǎng)格技術及VOF方法，在恒定流情況下，對側墻長廊道單明溝、雙明溝及三明溝消能工充水過程進行三維數(shù)值模擬。分析和比較各型式明溝的剩余比能、動能及流場分布，綜合評判各型式明溝消能工的消能特性。研究結果表明，3種明溝型式比較，三明溝布置消能效果最優(yōu)，單明溝和雙明溝各項消能指標較接近，但雙明溝消能條件下，其布置更為靈活，消能空間更大，水流分布亦較單明溝更為均勻，尤其是在減小閘室內(nèi)橫向水流對船舶的作用力方面較單明溝型式優(yōu)越。

三維數(shù)值模擬；明溝消能工；消能效果

Biography：LV Wei-dong（1968-），male，senior engineer.

對于船閘閘室消能工，早在20世紀30年代，美國開始在部分船閘底部橫支廊道前沿采用單明溝消能工［1］。國內(nèi)從20世紀50年代開始，陸續(xù)在長江三峽船閘、葛洲壩二號等船閘設計、施工中也采用單明溝消能工［2-4］；21世紀初開始，國內(nèi)相關研究單位對閘底廊道雙明溝消能布置進行了物理模型試驗研究，并應用于工程實踐；隨著國內(nèi)一批中低水頭、大尺度閘室船閘側墻廊道工程的建設，由于其輸水廊道布置于閘墻內(nèi)，閘室內(nèi)有足夠的空間布置多明溝消能。因此，需要對其消能特性進行研究。目前采用的研究方法主要以物理模型試驗為主，近幾年通過大量的物理模型試驗和原型觀測成果，使得對輸水時不同消能工的消能機理及水流消能引起的局部水流流態(tài)數(shù)值模擬研究和計算方法已趨一定進展，但系統(tǒng)建立三維數(shù)學模型，模擬計算、評判3種明溝型式的消能特性、消能效果的研究資料較少。因此，以貴港二線船閘為依托，選取其閘室中間段5個支孔為研究對象，利用本數(shù)學模型對各種明溝型式消能工進行三維數(shù)值模擬，系統(tǒng)的分析各明溝消能工消能特性、消能效果，并提出最優(yōu)的明溝型式。為進一步進行物理模型試驗提供方案參考。

1 計算模型及工況

1.1 控制方程

計算模型選用紊流模型中的雷諾時均法RNG k-ε紊流模型。該模型與標準湍流模型的主要區(qū)別：（1）方程中的常數(shù)是用理論推導得出的，并非用實驗方法確定；（2）耗散方程及系數(shù)體現(xiàn)了平均應變率對耗散項的影響，因此對高速射流的各向異性能較好地模擬。因此，本研究采用RNG k-ε紊流模型，其連續(xù)方程、動量方程和k、ε方程可分別表示如下

式中：ρ和μ分別為體積分數(shù)平均的密度和分子粘性系數(shù)；vt為紊流粘性系數(shù)，它可由紊動能k和紊動耗散率ε求出。

以上各張量表達式中，i=1,2,3，即{xi=x,y,z}，{ui=u,v,w}；j為求和下標，方程中通用模型常數(shù)取值為：h0=4.38,b=0.012,Cm=0.085,C2e=1.68,Sk=0.717 9,Se=0.717 9。

對閘室水面模擬，采用VOF模型。引入VOF模型的k-ε紊流模型方程（1）～（4）與單相流的k-ε模型形式是完全相同的。只是密度ρ和μ的具體表達式不同，它們是由體積分數(shù)加權平均值給出，也就是說，ρ和μ是體積分數(shù)的函數(shù)，而不是一個常數(shù)。它們可由下式表示

式中：αw為水的體積分數(shù)；ρw和ρa分別水和氣的密度；μw和μa分別為水和氣的分子粘性系數(shù)。通過對水的體積分數(shù)的αw迭代求解，ρ和μ值都可由式（5）～（6）求出。水的體積分數(shù)αw的控制微分方程為

本文采用有限容積法對上述偏微分方程組進行離散，為反應壓力對速度的影響，壓力—速度耦合求解使用S I M P L E R法。

圖1 各型式明溝消能閘室網(wǎng)格剖分圖Fig.1 Grid subdivision graph of each open ditch energy dissipation lock chamber

1.2 計算段閘室

以貴港二線船閘閘室中間段5個支孔、取閘室寬度的一半為計算對象，計算閘室段幾何尺寸為4 2.5 m×1 7 m×1 7 m（長×寬×高），數(shù)學模型采用1：2 0比尺，所以計算段閘室尺寸為2.1 2 5 m×0.8 5 m×0.8 5 m（長×寬×高）。為使閘室內(nèi)水流接近于恒定流，在起始水深以下開一個0.6 m×0.2 m泄水閥門。分別建立單明溝、雙明溝及三明溝型式的閘室三維模型，采用六面體和五面體對計算域進行剖分［5－7］。圖1為三明溝消能閘室網(wǎng)格剖分圖，圖2為雙明溝細部布置圖。

1.3 邊界條件處理

對于閘室內(nèi)存在的自由表面，采用三維V O F法進行計算。對于進出口條件，由于閘室內(nèi)水流為恒定流，所以進出口均有流量控制。固壁邊界條件由壁函數(shù)方法給定［8－10］。

圖2 雙明溝布置圖Fig. 2 Arrangement of double open ditches

1.4 計算工況

本文計算工況采用5 m i n開門，閥門開度為n=0.8時水位及流量，其計算模型流量為2 7.5 L/s，計算水深為6 0 c m。

2 計算結果分析

2.1 橫斷面流場分布

圖3顯示了3個消能明溝型式橫向剖面流場分布，可見3個剖面流場明顯不同。單明溝流態(tài)為：水流出支孔后射入閘室，受單明溝內(nèi)低檻阻擋，水股上抬從低檻表面水平射向閘室中部，同時水股分散，受明溝邊墻阻擋，水流上升，在閘室內(nèi)形成逆時間方向旋渦；雙明溝流態(tài)為：水流出支孔后射入閘室，水流受擋板阻擋，部分水流分散上升，部分水流從擋板中孔進入閘室后部，受明溝邊墻阻擋，水流上升，在閘室內(nèi)形成2個逆時間方向旋渦。三明溝流態(tài)為：水流出支孔后射入閘室，水流受第一道擋板阻擋，部分水流分散上升，部分水流從擋板中孔進入第二道明溝，受第二道擋板阻擋，部分水流進一步分散上升，部分水體通過第二道擋板中孔進入閘室中部，受明溝邊墻阻擋，水流上升，在閘室內(nèi)形成2個相反方向旋渦，尺度均較小，第一個旋渦位于第二道明溝之前，為逆時間方向旋渦；第二個旋渦位于閘室后部，為順時間方向旋渦。

圖3 橫斷面流場分布Fig.3 Transverse flow field distribution

圖4 縱向剖面流場分布（Y=0.2 m）Fig. 4 Longitudinal profile flow field distribution（Y= 0.2 m）

2.2 縱剖面流場分布

圖4是Y=0.2 m縱剖面流場二維顯示，可見單明溝受低檻阻擋，明溝底部水流較為分散，雙明溝和三明溝由于擋板透水，流股較為明顯。從流速分布來看，3種型式明溝底部流速均較大，由于三明溝消能更為充分，在Z≥0.2 m區(qū)域內(nèi)，流速基本均勻，且較小。而單明溝和雙明溝消能效果較低，在Z≥0.2 m區(qū)域內(nèi)流速仍較大，且分布不均。Y=0.4 m剖面（圖5）與Y=0.2 m剖面規(guī)律和結論基本相同。

2.3 消能效果比較

對某一水面剖面，剩余比能按式（8）計算，計算結果如表1所示。

由圖6可看出各明溝型式的剩余比能。在無論何種水深下，三明溝的剩余比能均較單明溝和雙明溝要小，雙明溝及單明溝隨著水深的不同剩余比能值互有交替。就剩余比能來看三明溝消能效果最優(yōu)，單明溝和雙明溝指標較接近，消能效果也較為接近。

表1 各水深斷面剩余比能及動能統(tǒng)計表Tab.1 Residual ratio and kinetic energy of each depth section

圖5 縱向剖面流場分布（Y=0.4 m）Fig.5 Longitudinal profile flow field distribution（Y= 0.4 m）

為更好的比較3種型式的消能工，提取各斷面流速按式（9）計算了各斷面動能，計算結果如圖7所示。

從圖7中可見，三明溝消能設施布置各剖面能量均小于單明溝和雙明溝，說明三明溝消能效果最優(yōu)。而單明溝與雙明溝各剖面能量基本相當，而在水面處單明溝能量還略大于雙明溝，說明雙明溝消能效果略好于單明溝。

圖6 剩余比能比較圖F i g.6 C o m p a r i s o n o f r e s i d u a l r a t i o

圖7 各水深斷面動能比較圖Fig.7 Comparison of kinetic energy for each depth section

4 結論

（1）經(jīng)三維數(shù)值模擬閘室內(nèi)各型式明溝消能工消能效果比較，其結果表明，三明溝消能設施較單明溝、雙明溝，流場分布更均勻，且剩余比能及動能等指標均較優(yōu)。說明三明溝較單、雙明溝消能效果更好，對閘室內(nèi)水流的二次調(diào)整作用更明顯；而雙明溝與單明溝消能設施比較，雖然其閘室內(nèi)縱橫向流場分布、剩余比能指標較接近，但其斷面動能在水面處單明溝能量略大于雙明溝，說明雙明溝消能效果略好于單明溝。

（2）明溝型式的消能機理主要是在出流擋板末端消能，孔口出流的突擴生成大量的旋渦以及水流的相互劇烈摻混，明溝及閘室中部的水體進行質量、動量和能量的交換，由于內(nèi)外摩擦作用，在水流動能轉變?yōu)閯菽艿倪^程中，流場急劇改變，致使水流內(nèi)部相對運動加劇、碰撞和摩擦作用亦加強，消耗了大部分能量。

［1］U S Army Corps Engineers．Hydraulic Design of Nay Lgation Locks［R］．Washington：U S Army Corps.，1975.［2

［2］JTJ306-2001,船閘設計規(guī)范［S］.

［3］須清華．葛洲壩工程二號船閘輸水系統(tǒng)水工模型終結試驗報告［R］．南京：南京水利科學研究院，1981．

［4］宗慕偉.三峽雙線連續(xù)五級船閘總體運行水力學試驗研究［R］.南京：南京水利科學研究院，1995.

［5］劉平昌,彭永勤,王召兵.貴港二線船閘側墻廊道輸水系統(tǒng)水力學試驗研究報告［R］.重慶：重慶西南水運工程科學研究所,2010.

［6］張緒進.貴港航運樞紐二線船閘輸水系統(tǒng)水力學試驗研究［J］.水利水運工程學報,2012（8）：34-38.

ZHANG X J. Research of Hydraulic Experiment About Second-line Ship lock′s Filling-emptying system of Guigang Navigation Project ［J］.Hydro-science and engineering, 2012（8）:34-38.

［7］張緒進，宣國祥，劉平昌.船閘輸水過程消能機理及船舶系纜力計算方法研究總報告［R］.重慶：重慶交通大學，2 0 1 3.

［8］宣國祥，李君.船閘輸水系統(tǒng)水力學創(chuàng)新技術研究與實踐［R］.南京：南京水利科學研究院,2009.

［9］楊忠超.高水頭船閘閘室明溝消能效果及機理分析［J］.水運工程,2009（12）：169-173.

YANG Z C.Energy dissipation effect and mechanism of open ditch energy dissipater in high-head lock［J］. Port & WaterwayEngineering,2009（12）：169-173.

［10］楊忠超.高水頭船閘閥門段體型優(yōu)化三維數(shù)值模擬［J］.水利水電科技進展，2010（4）：10-14.

YANG Z C. Three dimensional numerical simulation of figure optimization of valves of high head Iock［J］.Advances Science andTechnology of Water Resources，2010（4）:10-14.

Three-dimensional numerical simulation of energy dissipation effect in energy dissipater of lock chamber open ditch

LV Wei-dong1,2,LIU Ping-chang1,3,CHEN Liang1,ZHANG Xu-jin1,3
（1.Chongqing Jiaotong University,Chongqing400074,China;2.Guangxi Xijiang Development&Investment Group Co.Ltd.,Nanning530022,China;3.Southwestern Research Institute of Water Transport Engineering,Chongqing400016,China）

By using the multi-grid technique and VOF （Volume of Fluid）method,the local mathematical model of lock chamber was developed.The energy dissipation process of the open ditch,double ditches and three ditches were simulated by three-dimensional model in the steady flow circumstances.Then the distribution of residual energy,kinetic energy and flow field were analyzed,and the energy dissipation effects of all types of energy dissipaters were compared in order to point out the optimal open ditch.The results show that various criteria on energy dissipation of open ditch are close to those of double ditches,and the energy dissipation effect of three ditches is optimized,but the energy dissipater of double ditches has much larger space.

three-dimensional numerical simulation;energy dissipater;energy dissipation effect

U 641；O 242.1

A

1005-8443（2013）06-0508-05

2013-03-07；

2013-04-15

交通運輸部科技項目（20118323501580）

呂偉東（1968-），男，廣西陸川人，高級工程師，主要從事內(nèi)河航運的開發(fā)、運營及研究工作。