船閘閘墻長廊道輸水系統閘室三維流場數值模擬研究

彭永勤，彭 濤

(1. 重慶交通大學 西南水運工程科學研究所，重慶 400016；2. 中煤科工集團重慶設計研究院有限公司，重慶 400016)

?

船閘閘墻長廊道輸水系統閘室三維流場數值模擬研究

彭永勤1，彭 濤2

(1. 重慶交通大學 西南水運工程科學研究所，重慶 400016；2. 中煤科工集團重慶設計研究院有限公司，重慶 400016)

采用流場計算分析軟件Fluent，在利用實測數據對數值模型進行驗證后，建立三維數學模型模擬計算了船閘閘墻長廊道輸水系統充水過程中的閘室流場，并對輸水過程中閘室流場垂向演變規律進行了分析。研究表明：在尺度空間一定和各孔射流比較均勻的條件下，增加側支孔的數量，漩渦數量增加，漩渦范圍減小，有利于船體停泊。

航道工程；輸水系統閘室；三維流場；數值模擬

閘墻廊道側支孔輸水系統的型式在船閘工程設計中是靈活多變的[1-4]。目前，沒有一種“最具代表性的閘墻廊道側支孔輸水系統”的閘室水動力學特性能夠涵蓋其他閘墻廊道側支孔輸水系統[5-10]。從理論研究的角度，一種兼具多種閘墻廊道側支孔輸水系統型式共同點的、比較簡單的“概念化閘墻廊道側支孔輸水系統型式”或可以為研究復雜的閘室水動力學問題奠定一個比較通用的基礎。筆者應用流場計算分析軟件Fluent，利用實測數據對數值模型進行驗證，計算模擬船閘閘墻長廊道輸水系統閘室三維流場，分析輸水過程中閘室流場變化規律。

1 網格剖分及邊界條件

本次研究對應物理模型按照1∶1的尺寸比例建立閘室水動力學數學模型。由于數學模型的邊界形狀比較復雜，三維計算網格類型采用四面體非結構化網格。在數學模型計算區域內，側支孔及孔口附近水流流態變化比較復雜，在劃分計算網格時將此區域進行比較細密地劃分，如圖1。計算區域內非結構化網格四面體單元總數為1 304 240個，計算節點數為274 016個，單元格的體積控制在0.63～1 719.85 mm3之間。

圖1 三維數值計算網格

閘室水動力學數學模型的進口邊界類型為“速度-時間”邊界，流速方向垂直進口斷面。根據一般經驗，進口斷面的流速分布對孔口斷面的流速分布影響不大，故在保證流量一致的條件下進口流速取平均值。

進口邊界條件的表達式為：

式中：Q(t)為主廊道流量；s為主廊道斷面面積。

出口運動流體為空氣，采用恒定壓力邊界條件。在數學模型中，由于模型靜止不動，與物模相對應的閘室側壁面、閘底、消力檻壁面、側支孔壁面采用了無滑移壁面條件。閘墻廊道側支孔輸水系統閘室水動力學數學模型的計算區域的輪廓和邊界條件設置如圖2。

圖2 計算區域及邊界條件

2 數學模型驗證

2.1 物理模型試驗

閘室模型尺度為2 000 mm×300 mm ×400 mm(長×寬×高)。閘室兩側輸水主廊道底部與閘室底部高度持平，斷面尺寸為60 mm×60 mm。每支廊道設等間距9個側支孔，側支孔尺寸為20 mm×20 mm，支孔間距為150 mm。閘墻兩側側支孔水平交錯布置，底部與廊道及閘室底部齊平，孔軸線與閘室縱軸線垂直。兩側主廊道由分岔管與上游銜接，分岔管呈對稱布置，以確保上游來流向兩側主廊道分配的流量基本相等。分岔管上游設平板閥門，閥門控制的廊道面積為120 mm×60 mm。平板閥門上游設水箱，溢流水箱內設平水槽，保證閘室充水過程中水位基本不變，以模擬船閘運行過程中的上游水位狀態。

船閘閘墻廊道側支孔輸水系統閘室水動力學試驗物理模型由閘室、側支孔、輸水廊道、岔管、平板閥門、溢流水箱(帶平水槽)、鋼架、水池、輸水管道、水泵及閥門等部分組成，模型如圖3。

圖3 物理模型

閘室瞬態流場采用現階段研究復雜流場問題最先進的工具——圖像粒子測速儀(PIV)進行測定。

2.2 驗證結果

由于閘室充水過程的水流為非恒定流，通過物理模型試驗和數學模型計算得到的數據量較大。根據速度矢量的時間和空間分布，整體比較數值模擬結果和試驗結果的工作是冗余的。本次采用抽取特定時刻和空間的樣本比較的方法，抽取物理模型和數學模型樣本的時刻為5 s，特征斷面與PIV測量的斷面相同，即距模型閘底高1 cm處的水平斷面。得出左側3#和9#支孔的流場對比結果如圖4、圖5。

圖4 3#支孔流場

圖5 9#支孔流場圖

由圖4、圖5可知，采用重正化群k-ε紊流模型模擬的閘墻長廊道輸水系統充水階段閘室三維流場與閘室水動力學物理模型試驗的測量結果吻合良好，證明所建立的數學模型是準確的，可以用于模擬實際的閘室充水過程，且其計算結果將作為閘室水動力特性分析的依據。

3 結果分析

采用側支孔輸水系統布置型式的船閘在充水過程中閘室水流形態最復雜的區域集中在閘室底部側支孔布置段。筆者圍繞上述時間、空間區域，對船閘充水過程的水動力學特性進行分析。

3.1 閘室底部特征斷面的流態

選取距離閘室底面1 cm(即z=1 cm)的水平截面作為閘室底部的水流形態的特征斷面。側支孔的高度為2 cm，所選取的特征面處于側支孔豎向中面位置，可以較清晰地反應底部側支孔射流的形態。特征斷面上的流速分布如圖6。

圖6 閘室特征斷面不同時刻的流態(z=1 cm)

在輸水初始時刻t=0.5 s〔圖6(a)〕，各側支孔射流的流量沿主廊道從上游到下游呈現逐漸減小的趨勢，射流的中心線沿側支孔軸線延伸?？拷嫌蔚膫戎Э咨淞髟谀┒私诲e，下游側支孔射流處于互不干擾的狀態。2 s時刻及以后〔圖6(b)、圖6(c)〕，閘室兩側側支孔射流處于充分的交錯擺動狀態，擺動幅度較大，劇烈程度隨時間不斷變化。

3.2 閘室底部流態的垂向演變

一般來說，閘室充水之前閘室底部與船舶底部之間有一定的富裕水深。在閘室設計中保持一定的富裕水深主要是為了避免閘室底部流速較大、流態復雜的水體大面積地與船體接觸而導致船體受到過大的水流沖擊力。因此，在閘室中與船體接觸的主要是靠近閘室自由水面的上部水體。閘室底部水體的流態必然影響上部水體。閘室底部流態的垂向演變反映了閘室底部水體的動力因素如何傳遞至上部水體，是研究閘室水動力學特性和船舶停泊條件重要內容之一。閘室底面以上6 cm(即z= 6 cm)處水平斷面流態的數值模擬結果如圖7。

圖7 閘室特征斷面不同時刻的流態(z=6 cm)

距閘室底面6 cm處的水平斷面相對1 cm處的水平斷面，流速絕對值降低，絕大部分區域在0.15 m/s以下。閘室兩側壁面區域的流速值相對中間區域的流速值較高，是射流末端遇對面邊壁阻擋后，轉而沿邊壁上升形成的。斷面流態呈多漩渦排列狀，沒有明顯射流流速分布狀態。相鄰漩渦方向相反，各漩渦范圍相對閘底斷面上對應位置漩渦的范圍擴大。通過分析上述特征斷面流態可知，閘室底部側支孔射流形成的流態隨著水深增加，射流流態逐漸消失，漩渦范圍擴大，流速降低，存在射流末端沿邊壁上升現象。直觀推斷，漩渦以順時針和逆時針依次排列、流速降低和射流主流區消失等現象，對船體在閘室上部水體內受力均勻是有利的。

4 結 語

閘室底部流態的垂向演變的根本原因是水流與閘室邊壁的摩擦作用和內摩擦作用。影響水流與邊壁的摩擦作用的主要因素為貼壁射流區域的大小、射流流速和射流方向等。在特定來流條件下，上述因素與主廊道和側支孔設計型式、邊壁形狀緊密聯系?？梢酝茢?，通過調整側支孔和閘室底面的形狀，增強射流與邊壁的內摩擦在一定條件下對于形成更有利于船舶停泊的上部水體是有益的。

射流主流區之間的相互摩擦是導致底面和上層面漩渦存在的主要因素。底面漩渦在向上延伸過程中，并沒有變得雜亂無章，而是有次序地排列。如果船體受順時針和逆時針漩渦擾動比較平均，則在垂向上所受的力矩則較小。漩渦的數量與交錯布置的側支孔的數量是對應的，在尺度空間一定和各孔射流比較均勻的條件下，增加側支孔的數量，漩渦數量增加，漩渦范圍減小，有利于船體停泊。

[1] 黎賢訪.船閘閘墻廊道側支孔輸水系統水力學研究[D].南京:南京水利科學研究院,2008. Li Xianfang.Hydraulics of Culvert Side Orifice Filling and Emptying System for Navigation Locks[D].Nanjing:Nanjing Hydraulic Research Institute,2008.

[2] 楊朝東.船閘復雜輸水系統輸水數學模型及應用研究[D].南京:南京水利科學研究院,1997. Yang Chaodong.A Mathematical Model of Unsteady Flow in Complicated Filling and Emptying System and Its Application[D].Nanjing:Nanjing Hydraulic Research Institute,1997.

[3] 張瑞凱,宣國祥.帶調節水池船閘水力計算數學模型[J].水利水運科學研究,1995(1):33-41. Zhang Ruikai,Xuan Guoxiang.Numerical model of hydraulic calculation for lock with regulation ponds[J].Hydro-Science and Engineering,1995(1):33-41.

[4] 吳時強,宣國祥,張瑞凱.三峽船閘末級輸水長廊道水力特性數值模擬[R].南京:南京水利科學研究院,1996. Wu Shiqiang,Xuan Guoxiang,Zhang Ruikai.Numerical Simulation on Hydraulic Property of Final Stage Long-Culvert Filling of Three Gorges Ship Lock[R].Nanjing:Nanjing Hydraulic Research Institute,1996.

[5] 楊忠超,陳明棟,楊斌,等.高水頭船閘閘室消能工消能效果實驗研究[J].重慶交通大學學報:自然科學版,2010,29(3):461-465. Yang Zhongchao,Chen Mingdong,Yang Bin,et al.Experimental study on energy dissipation effect of chamber energy dissipater in superhigh-head lock[J].Journal of Chongqing Jiaotong University:Natural Science,2010,29(3):461-465.

[6] 陳明,宣國祥,陳明棟.船閘輸水系統水動力學研究綜述[J].重慶交通大學學報:自然科學版,2013,32(1):157-160. Chen Ming,Xuan Guoxiang,Chen Mingdong.Review of hydrodynamics of lock filling and emptying system[J].Journal of Chongqing Jiaotong University:Natural Science,2013,32(1):157-160.

[7] Deardorff J W.A numerical study of three-dimensional turbulent channel flow at large reynolds numbers[J].Journal of Fluid Mechanics,1970,41(2):453-480.

[8] Harlow F H,Welch J E.Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface[J].Physics of Fluids,1965,8(12):2182-2189.

[9] Richard L S.Modeling hydrodynamic forces on vessels during navigation lock operations[C]// Hydraulic Measurements and Experimental Methods 2002.U.S.A.:American Society of Civil Engineers,2004:1-9.

[10] Richard L S.Mooring model coefficients for barge tows in a navigation lock[J].Journal of Waterway,Port,Coastal and Ocean Engineering, 2003,129(5):233-236.

Three Dimensional Flow Fields Numerical Simulation of Long CorridorFilling-Emptying System Chamber of Ship Lock Wall

Peng Yongqin1, Peng Tao2

(1. Southwestern Hydro Engineering Research Institute for Water Way, Chongqing Jiaotong University,Chongqing 400016, China;2. China Coal Technology & Engineering Group, Chongqing Design & Research Institute, Chongqing 400016, China)

After using the measured data to validate numerical models, the flow field analysis software Fluent was adopted to build a three-dimensional mathematical model, for calculating the chamber flow field in the process of filling water, and analyzing evolution of vertical flow field in the filling process. The study shows that when scale space is certain and each hole jets uniformly, if the number of collateral hole increases, the number of swirl also increases, but the swirl range is reduced, which is in favor of the hull moored.

waterway engineering; filling-emptying system chamber; three dimensional flow fields; numerical simulation

10.3969/j.issn.1674-0696.2015.03.15

2013-10-30；

2014-03-12

彭永勤(1984—)，女，重慶人，工程師，碩士，主要從事通航論證方面的研究。E-mail: 174588703@qq.com。

U612

A

1674-0696(2015)03-072-04