魏文禮， 邵世鵬， 劉玉玲

(西安理工大學 水利水電學院，陜西 西安 710048)

梯形斷面明渠丁壩繞流水力特性三維大渦模擬

魏文禮， 邵世鵬， 劉玉玲

(西安理工大學 水利水電學院，陜西 西安 710048)

對梯型斷面明渠非淹沒式單一丁壩繞流水力特性進行三維數(shù)值模擬，為丁壩的設計和施工提供理論基礎支持。數(shù)值模擬采用兩相流混合模型，并分別選取大渦模型(LES)和RNGk-ε湍流模型封閉兩相流時均方程。速度與壓力的耦合使用半隱式SIMPLE算法，模擬自由水面采用了VOF法。通過不同截面流線圖和流速矢量圖的比較得出大渦模型能更好地捕獲水流瞬時流動特性，動態(tài)再現(xiàn)二次流動結構；并將其模擬的丁壩后回流區(qū)域邊界線與試驗值對比，兩者吻合良好，表明大渦模型能夠很好地模擬明渠丁壩繞流的水力特性分布規(guī)律。

丁壩繞流； 數(shù)值模擬； 大渦模型； RNGk-ε模型； 水力特性

丁壩是廣泛使用的一種水工建筑物，其主要作用是保護河岸免受水流沖刷、束窄河道提高河道通航能力等。修建丁壩后水流被導向河流中心，所以河流主流區(qū)的流速明顯提高[1]。早期的研究者主要采用實驗的方法對丁壩進行研究，但隨著計算機技術和計算流體力學的逐漸成熟，采用數(shù)值模擬的方法對丁壩進行研究已被越來越多的研究者所接受。夏云峰等[2]將二維的水流運動方程與SIMPLER方法結合，得到了丁壩水流流態(tài)及水面線變化規(guī)律，且與實測值的吻合度較高；鄧紹云等[3]采用數(shù)值模擬的方法對一非淹沒丁壩的繞流水力特性進行了研究，并對水流流速場和壓力的分布規(guī)律進行了分析，其符合實際的丁壩繞流水力特性，說明數(shù)值模擬結果達到了捕捉丁壩周圍水流流態(tài)的目的；周宜林等[4]采用大渦模擬方法研究了不同挑流形式下的丁壩水力特性，從模擬結果可以看出，不同挑流形式對水流影響較小，但是相比較而言，當丁壩采用下挑形式時可以更好地保護壩頭；李冰凍等[5]通過研究丁壩試驗和數(shù)值模擬兩方面的成果，分析了丁壩附近區(qū)域復雜的水流特征，通過比較得出數(shù)值計算能夠反映丁壩擋水所引起的水流流態(tài)的變化過程；Akahori R等[6]用大渦模型對丁壩附近二次流進行了數(shù)值模擬；Mayerle R[7]和Jia Y等[8]采用數(shù)值模擬的方法對丁壩附近水流的流動規(guī)律進行了研究；白靜等[9]將數(shù)值模擬和PIV測量實驗結合在一起，對不同丁壩間距下丁壩群附近的水力特性進行了研究，得出丁壩的長度與丁壩之間距離的比值L/D對湍動強度，渦量分布和流場結構等有很大的影響；蔣昌波等[10]采用簡化的二維大渦模型，對一非淹沒丁壩群的丁壩繞流水力特性進行了研究，并將計算的結果與實驗結果相比較，結果吻合良好，說明采用簡化的大渦模型可以很好地反映丁壩繞流中的流場、流線及渦的變化等水力特性；辛永政等[11]采用不同的紊流數(shù)值模型對明渠丁壩繞流進行了數(shù)值模擬并與試驗值進行了比較，結果表明標準模型可以比較準確地模擬明渠中丁壩繞流的情況。雷亞等[12]采用FLUENT軟件對一非淹沒丁壩進行了模擬研究，自由水面的捕捉采用VOF法；得出丁壩后回流區(qū)范圍會因丁壩類型而變化；曹曉萌等[13]運用數(shù)值模擬的方法運用不同劃分準則對一非淹沒丁壩進行了研究，并將計算結果與實驗值進行了對比，結果吻合良好，說明數(shù)值模擬達到了反映實際丁壩繞流水力特性的目的。

現(xiàn)有對丁壩繞流水力特性的研究都是在矩形斷面水槽中進行的，而天然河道斷面大多不是矩形的，天然河道斷面更接近梯形斷面，因此有必要對梯形斷面水槽中丁壩繞流的水力特性進行數(shù)值模擬，并通過數(shù)值模擬比較選出能更好反映實際河道中丁壩繞流水力特性的數(shù)學模型。

1 數(shù)學模型

利用空間濾波函數(shù)對三維湍流N-S方程進行濾波后，得大渦模型的控制方程[14-15]為：

(1)

(2)

這表明小尺度變量對大尺度變量有作用，必須進行建模。

LES湍流模擬方法假定湍流運動可分為大尺度渦和小尺度渦，兩者之間的分離對大渦的演變不會有大的影響，并且小渦受水流的幾何形狀和邊界條件的影響較小。

Smagorinsky模型是基于渦粘性假設的最常用的亞格子雷諾應力模型。

LES直接由一組濾波控制方程解決大尺度場，亞格子尺度是各向同性的，主要用來表示耗散，所以可以很容易地模擬其對網(wǎng)格尺度的影響[16]。對于這種模型，亞格子雷諾應力假定與應變速率張量成比例，即：

(3)

以 RNGk-ε模型封閉的描述氣液兩相流的時均控制方程見文獻[17]。

自由水面的捕捉采用VOF法，其基本原理見文獻[18-19]。

2 數(shù)值模擬及驗證

2.1 計算區(qū)域及網(wǎng)格劃分

本研究要進行數(shù)值模擬研究的水槽長為L=8.0 m，底寬B=1.5 m，高H=0.2 m，邊坡系數(shù)為1.33。

數(shù)值模型參數(shù)：丁壩為方頭丁壩，布設在水槽的右側，丁壩迎水面軸線位于距計算入口3.0 m的地方，丁壩下底軸線長為D=0.2 m，丁壩高度為0.2 m，丁壩厚度為d=0.008 m，水槽水深為h=0.15 m。進口流量恒定為Q=0.05 m3/s。

計算區(qū)域體型圖如圖1(a)所示，圖1(b)為水槽橫斷面尺寸示意圖。

2.2 網(wǎng)格劃分

初始計算網(wǎng)格由GAMBIT程序生成。

計算時采用結構化網(wǎng)格進行劃分，沿水流方向計算網(wǎng)格劃分為漸變網(wǎng)格，在丁壩處加密，網(wǎng)格總數(shù)為36 398個，見圖2。

2.3 邊界條件

水流進口設為速度進口，速度為0.196 m/s；空氣進口也設為速度進口，速度為0.196 m/s。

水流出口邊界設定為壓力出口， 初始水位為0.15 m；空氣出口邊界也為壓力出口，相對壓強為0。水與空氣接觸面的相對壓強設定為0。

水槽和丁壩固體壁面上的邊界條件符合流速在壁面上的不分離和無滑移條件。

模型的離散采用有限體積法；邊界的擬合采用非規(guī)則網(wǎng)格。自由水面按VOF法確定。

2.4 結果分析和討論

2.4.1 水平方向流線圖與流速矢量圖分析

不同水深平面處的流線圖如圖3、4所示。

由圖3、4可以看出，由于丁壩的存在，部分水流繞過丁壩流向下游，而部分水流則沿丁壩上游壩面向水槽底部下潛，然后繞過丁壩壩頭流向下游，在丁壩附近形成主流區(qū)和回流區(qū)；而且在丁壩后都形成了一個大的二次流，沿順時針方向流動，這符合實際繞壩水流的運動規(guī)律。從流線圖可以看出,由RNGk-ε模擬的渦規(guī)則均勻，能夠反映出強流線彎曲；而大渦模型把湍流分成大尺度和小尺度兩部分運動，能夠更好地模擬丁壩后大尺度渦的基本特征，模擬的流線也更符合實際回流的流態(tài)。大渦模型模擬的丁壩后回流區(qū)域邊界線長度比RNGk-ε模擬的小，而寬度較大；大渦模型模擬的主流區(qū)比RNGk-ε模擬的窄，說明大渦模擬能更好地反映丁壩后回流區(qū)對主流區(qū)的影響。大渦模型模擬的結果捕捉到了丁壩前后與邊坡拐角處的小漩渦；而RNGk-ε模型只捕捉到了丁壩前與邊坡拐角處的小漩渦。說明大渦模型能夠捕獲水流瞬時流動特性，動態(tài)再現(xiàn)二次流動結構。同時由圖3、4也可以看出隨著水深的增加丁壩后回流區(qū)域逐漸變小。

由圖5、6可以看出在丁壩前水流流速矢量分布均勻，在靠近丁壩上游壩面處流速減慢水流雍向壩面后轉向，流速矢量在壩頭處和主流區(qū)分布密集，流速變化劇烈；在丁壩前后流速矢量分布稀疏，速度較小。同時可以看出在丁壩后是一個死角，水流速度很小，渦也很弱。計算結果也表明隨著渦的擴散，在丁壩下游一段距離后水流流速逐漸恢復均勻。

2.4.2 垂直方向流線圖與流速矢量圖分析

由圖7可見，在丁壩前RNGk-ε模型和大渦模型模擬的流線圖都很均勻，但在靠近丁壩上游壩面處大渦模型模擬的流線圖能更好地反映出丁壩前水流的壅水現(xiàn)象。在丁壩后RNGk-ε模型捕捉到了水流在垂直截面上的回流，流線分布均勻，但只是渦的局部；大渦模型模擬的回流比較完整，得到的渦沒有RNGk-ε模型模擬的規(guī)則，但更符合實際水流的流動，可以看出隨著渦的擴散，在近底面處水流流向恢復到與進口水流相同的方向，能更好地反映丁壩后水流在垂直截面上的流動狀態(tài)。由于水平方向的渦和垂直方向的渦的相互作用，形成典型的三維渦結構，因此在丁壩附近的水流是強三維紊流。

比較圖7,圖8可以看出丁壩對水槽中心斷面水流有一定的影響，但在丁壩附近的二次流現(xiàn)象消失了，水流方向相同且均勻，說明由于丁壩的作用不同垂直截面上的水流流線分布方式不一致，即不同垂直截面上水流流態(tài)不一致。

比較大渦模型和RNGk-ε模型模擬的丁壩壩軸線中心縱向斷面矢量圖(圖9(a),(b))可以看出，大渦模型比RNGk-ε模型模擬的丁壩前水流的上升趨勢明顯，能更好地反映丁壩前的壅水現(xiàn)象。大渦模型模擬的丁壩后水流分布表現(xiàn)出明顯的各向異性，有二次流產(chǎn)生，這更符合實際繞壩水流的流動特性，說明大渦模型能更好地捕獲水流瞬時流動特性，動態(tài)再現(xiàn)二次流動結構。

由圖10可以看出，在水槽中心斷面矢量均勻，二次流現(xiàn)象消失。

2.4.3 回流區(qū)域邊界線比較

由以上流線圖分析可知，大渦模型能更好地反映水流的瞬時流態(tài)和二次流的水流結構，所以將大渦模型模擬的水槽底面丁壩后的回流區(qū)域邊界線與試驗值相比較。由圖11可以看出計算模擬的回流長度與試驗值[20]吻合良好，但回流寬度有一定的誤差。主要原因是，在試驗時丁壩是圓頭的，而在數(shù)值模擬時考慮到丁壩很薄(0.008 m)，為方便畫網(wǎng)格，將丁壩近似畫為方頭，這樣導致了模擬中的繞壩水流沒有試驗中的平順，產(chǎn)生了一定的誤差。

3 結 論

本研究采用兩相流混合模型，并分別選取大渦模型和RNGk-ε湍流模型封閉兩相流時均方程，對梯型斷面明渠非淹沒式單一丁壩繞流水力特性進行三維數(shù)值模擬，由以上的分析和比較可以得出以下結論。

1) 分析比較由大渦模型和RNGk-ε模型模擬的流線圖和流速矢量圖得出，大渦模型能更好地捕獲丁壩繞流的瞬時流動特性和流動結構。

2) 大渦模型模擬的丁壩后回流區(qū)域邊界線與試驗得出的結果吻合良好，說明大渦模型可以很好地預測丁壩附近水流的回流區(qū)域。

3) 由水平方向和垂直方向流線圖可以看出，由于水平方向和垂直方向的渦相互作用，形成典型的三維渦結構使得丁壩附近水流為強三維紊流。

4) 由于丁壩的作用在丁壩后會形成一個大的順時針流向的二次流，在丁壩前后與底面的拐角處會形成小的角渦，在丁壩后會形成一個死角，水流速度很慢，渦也很弱，容易造成泥沙沉積。

[1]Wu B S, Wang G Q, Ma J M, et al. Case study: river training and its effects on fluvial processes in the lower Yellow River[J]. Hydraul English ASCE, 2005, 131(2): 85-96.

[2]夏云峰，孫梅秀，李昌華. 用水深平均 紊流模型計算淹沒丁壩[J].水利水運科學研究,1993,(2):109-118.

Xia Yunfeng, Sun Meixiu, Li Changhua.2-D depth averaged turbulence model for flow around submerged spur dike[J]. Journal of Nanjing Hydraulic Research Institute,1993,(2):109-118.

[3]鄧紹云，王義剛，邱清華. 非淹沒丁壩三維繞流數(shù)值模擬[J].水運工程,2010,(6):90-94.

Deng Shaoyun, Wang Yigang, Qiu Qinghua, 3D-circumferential flow value simulation to an un-submerged spur dike[J].Port & Waterway Engineering, 2010,(6):90-94.

[4]周宜林.淹沒丁壩附近三維水流運動大渦模擬[J]長江科學學院院報,2001,18(5):28-31.

Zhou Yilin. Division mechanism and flow scale investigation of circulation zone around spur dike[J]. Yangtze River Scientific Research Institute, 2001, 18(5):28-31.

[5]李冰凍，李嘉，李克鋒. 丁壩水流的水槽實驗及數(shù)值模擬研究[J].水動力學研究與進展,2013,28(2):176-183.

Li Bingdong, Li Jia, Li Kefeng. Flume experiment and numerical simulation on flow around a spur dike[J]. Chinese Journal of Hydrodynamics,2013,28(2):176-183.

[6]Akahori R，Schmeeck E M．Numerical analysis of secondary-flow around a spur dike using a three-dimensional free water surface LES model: Parker & García． River，Coastal and Estuarine Morphodynamics, RCEM 2005[C]. Taylor & Francis Group，London, 2006: 921-930.

[7]Mayerle R，Wang S S Y，Toro F M. Verification of a three-dimensional numerical model simulation of flow in the vicinity of spur-dikes[J]. Journal of Hydraulic Research, 1995, 33(2):243-256.

[8]Jia Y，Scott S，Xu Y，et al. Three-dimensional numerical simulation and analysis of flows around a submerged weir in a channel bend way[J]. Journal of Hydraulic Engineering, 2005, 131(8):682-693.

[9]白靜，方紅衛(wèi)，何國建. 非淹沒丁壩繞流的三維大渦模擬研究[J]. 力學學報,2013,45(2):151-156.

Bai Jing, Fang Hongwei, He Guojian. Study of non-submerged groin turbulence flow in a shallow open channel by LES[J].Chinese Journal of Theoretical and Applied Mechanics, 2013,45(2):151-156.

[10]蔣昌波，呂昕，楊宜章. 丁壩繞流的二維大渦模擬數(shù)值模擬[J].長沙交通學院學報，1999,15(3):68-72.

Jiang Changbo, Lü Xin, Yang Yizhang. Large eddy simulation of 2D flow around a spur dike numerical simulation[J]. Journal of Changsha Communications University, 1999,15(3):68-72.

[11]辛永政，張莉莉. 3種紊流數(shù)值模型模擬明渠丁壩三維水流初探[J]. 貴州水力發(fā)電, 2011,(5): 56-58.

Xin Yongzheng, Zhang Lili, Three kinds of turbulent flow numerical model for the simulation of three-dimensional flow in open channel[J].Guizhou Water Power, 2011，(5):56-58.

[12]雷亞，劉杰. 基于FLUENT的丁壩回流區(qū)數(shù)值模擬[J].水利科技與經(jīng)濟，2015,21(2):4-6.

Lei Ya, Liu Jie. Simulation of the backflow zone around spur dike based on FLUENT[J]. Water Conservancy Science and Technology and Economy, 2015,21(2):4-6.

[13]曹曉萌，顧正華.3種非淹沒雙體丁壩作用尺度劃分準則及比較[J].浙江大學學報:工學版，2015,49(2):200-207.

Cao Xiaomeng,Gu Zhenghua. Three classification criteria and their comparison impact scale between double non-submerged spur dikes[J].Journal of Zhejiang University(Engineering Science),2015,49(2):200-207.

[14]魏文禮，戴會超. 紊流模型理論及工程應用[M]. 西安:陜西科學技術出版社,2006.

[15]戴會超，王玲玲，魏文禮，等. 線性波中二維射流的大渦模擬[J]. 華中科技大學學報, 2006,34(9):45-47.

Dai Huichao, Wang Lingling, Wei Wenli, et al. Large-eddy simulation of the interaction between jets and linear waves[J].Journal of Huazhong University of Science and Technology, 2006,34(9):45-47.

[16]Ferziger J H. Large eddy simulation: an introduction and perspective[M]. New York: Springer-Verlag, 1996.

[17]戎貴文，魏文禮，劉玉玲，等. 涌潮作用下丁壩附近水流運動特性的數(shù)值模擬研究[J].水利學報,2012,43(3):296-301.

Rong Guiwen, Wei Wenli, Liu Yuling, et al. Study on flow characteristics near spur dikes under tidal bore[J]. Journal of Hydraulic, 2012,43(3):296-301.

[18]Shen Y M, NG C O, Zheng Y H. Simulation of wave propagation over a submerged bar using the VOF method with a two-equation κ-ε turbulence modeling[J]. Ocean Engineering, 2004, 31(1): 87-95.

[19]袁麗蓉，沈永明，鄭永紅. 用VOF方法模擬橫流下窄縫紊動射流[J]. 海洋學報,2005,(4): 155-160.

Yuan Lirong,Shen Yongming,Zheng Yonghong. Simulation on a slit-type turbulent jet in a cross-flow using the VOF method[J]. Act A Oceanologica Sinica,2005,(4): 155-160.

[20]馬騰云. 岸坡對丁壩回流影響的實驗研究[D].南京：河海大學，1988.

Ma Tengyun. Experimental investigations of shore effecting on the size of groyne's recirculating flow zone[D]. Nanjing: Hohai University,1988.

(責任編輯 王衛(wèi)勛)

3D numerical simulation of hydraulic behaviors of spur dike flow in a trapezoid open channel

WEI Wenli, SHAO Shipeng, LIU Yuling

(Faculty of Water Resources and Hydroelectric Engineering,Xi’an University of Technology，Xi’an 710048，China)

The 3-D numerical computation of flow hydraulic behaviors around non-submerged spur-dikes in a trapezoid open channel is carried out to provide theoretical basis support for design and construction of spur-dikes. The gas-liquid two-phase large eddy simulation and the gas-liquid two-phase mixture model with 3-D RNGk-εturbulence model are used in the simulation. The SIMPLE algorithm is used to solve velocity and pressure. The VOF method is used to simulate the free surface. The comparisons of the streamlines at different section simulated by the large eddy simulation and RNGk-εmodel show that the large eddy simulation can better capture the instantaneous flow behaviors and can simulate dynamic secondary flow structures. The comparisons of the biggest recirculation boundary line between the computed and the measured data are in a good agreement, which indicates that the large eddy simulation can well simulate the hydraulic characteristics of spur dike flow in trapezoid open channel.

spur dike flow; numerical simulation; large eddy simulation; RNGk-εmodel; hydraulic characteristics

1006-4710(2015)04-0385-06

2014-12-20

國家自然科學基金資助項目(51178391);陜西省科學研究計劃基金資助項目(2014K15-03-05)；陜西省西北旱區(qū)生態(tài)水利工程重點實驗室開放基金項目(106-221223);中央財政支持地方高校發(fā)展專項資金特色重點學科項目(106-00X101) 。

魏文禮，男，博士，教授，研究方向為環(huán)境水力學、水污染控制理論與技術。E-mail:wei_wenli@126.com。

TV131.4

A