方晨蕾，葛曉偉，王 芳，董詩媛

（1.高郵市水利綜合服務(wù)中心，江蘇 高郵 225600；2.高郵市水利局臨澤水利站，江蘇 高郵 225600）

水閘是建在河道、渠道及水庫、湖泊岸邊具有擋水和泄水功能的低水頭水工建筑物，對攔洪、擋潮、沖沙、取水或下游流量調(diào)節(jié)具有重要作用[1-3]。在水閘的水力參數(shù)設(shè)計中，流量設(shè)計十分重要，合理控制水閘開度不僅可以使水閘下游流量均勻分布，還可以防止兩岸和河床的水土流失[4，5]。

孫海超[6]在理論上分析推導(dǎo)了起推水位偏差對水面線的影響范圍，推算了閘下水位流量關(guān)系；迮振榮等[7]分析了影響平原河網(wǎng)地區(qū)感潮河段水閘下泄流量的成因，并通過成因分析研究了提高水閘流量控制方法；周全等[8]介紹了在線監(jiān)測、水文水動力學(xué)模型、水力學(xué)及多元回歸法過流公式3 種水閘流量計算方法，并計算了上海市浦東新區(qū)典型閘門流量過程；邢端生[9]以粵東某水閘為例，在分析流域產(chǎn)匯流特性基礎(chǔ)上，分別用廣東省綜合單位法和平均排除法計算山區(qū)和平原面積的設(shè)計流量；孟健[10]以水閘不同啟閉情況下的實測過閘流量為依據(jù)，結(jié)合水閘自身特性與水力學(xué)相關(guān)理論，探求可操作性的過閘流量推求方案；劉鵬等[11]對黃埔涌水閘正常和特殊設(shè)計工況過流流量進行了模型試驗研究，分析了幾種工況下閘門運行時開度與流量之間的關(guān)系。

本文通過水閘開度控制流量大小，采用Fluent數(shù)值模擬方法，引入能量耗散和流動阻力，分析對比了不同閘門開度下RNG k-ε、RNG k-ω 和RSM 湍流模型對流體速度分布和能量系數(shù)的影響，并通過室內(nèi)水閘試驗進行驗證對比，研究結(jié)果可為相關(guān)工程提供參考。

1 室內(nèi)試驗裝置

試驗在國內(nèi)某大學(xué)水力實驗室水槽上進行，水槽為玻璃墻和鋼底，尺寸為8 m×30 cm×50 cm（長×寬×深），水閘位于距離上游端4 m 處。流量由入口管上的閘門進行調(diào)節(jié)，通過放置在渠尾的V 形缺口堰進行測量。下游尾水深度由水槽尾部的閘門固定，有機玻璃門矩形門寬5 mm，采用精度為±0.5 mm的點式水位計測量水位。使用Nortek 的Vectrino 聲學(xué)多普勒測速儀測量流速，輸出采樣頻率為25 Hz，測量記錄時間為40 s，通過分析平均瞬時數(shù)據(jù)計算湍流特性。試驗裝置，如圖1所示。

圖1 試驗裝置

2 數(shù)值模擬

2.1 計算方案

數(shù)值分析采用非穩(wěn)態(tài)條件下雷諾平均Navier-Stockes 方程求解，可在70 s 的模擬時間內(nèi)求出一個恒定的解，相比穩(wěn)態(tài)計算，更具有時效性。計算時，采用SIMPLE 算法獲得壓力場，并采用PRESTO 框架將壓力離散化，其他方程采用二階離散格式，參數(shù)由Fluent Default 值給出。此外，湍流模型采用了3 種計算公式即RNG k-ε、RNG k-ω 和RSM，同時為了跟蹤自由表面，使用了基于VOF 公式的流體體積模型。在整個計算域中，考慮了水和空氣兩相狀態(tài)，流體的性質(zhì)也被考慮在內(nèi)，RANS 和連續(xù)性方程采用流體特性的體積分數(shù)平均值求解。

淹沒流二維RANS模擬的域和邊界條件，如圖2所示。試驗在淹沒流中進行，總水頭H0為0.3 m，閘門開度W為0.03～0.24 cm，下游水位h2為0.20～0.27 cm，相對閘門開度a=W/H0，淹沒度s=h2/H0，x為閘門距離上下游的距離。

圖2 淹沒流二維RANS模擬的域和邊界條件

2.2 數(shù)值網(wǎng)格

數(shù)值模型二維結(jié)構(gòu)網(wǎng)格通過Gambit 生成，尺寸為6 m×0.4 m，閘門距離下游邊界3 m 處為水平橫坐標(biāo)的原點。相對閘門開度a為0.1～0.8，創(chuàng)建了8 個不同的網(wǎng)格。網(wǎng)格在墻邊界（底部和水閘）附近和自由表面區(qū)域被細化。構(gòu)建網(wǎng)格是為了驗證閘門到墻的無量綱距離在12＜y+＜250的范圍內(nèi)。

2.3 邊界條件

邊界條件和初始值，如圖2 所示。水和空氣在上游端以2 種速度條件分別注入，空氣速度不敏感且等于零。假設(shè)水流速度呈指數(shù)分布，其公式為：

式中：U(y)為水流速度（m/s）；γ為形狀系數(shù)，根據(jù)從明渠中觀察到的標(biāo)準(zhǔn)值取0.1；U0為進水口處的平均流速（m/s）；y為水流垂直位置（m）；h0為河道深度（m）。

第一次計算采用Garbrecht 流量方程確定流量和平均速度，以細化自由表面的網(wǎng)格，并在下游端固定靜水壓力出口。

3 結(jié)果分析

3.1 剖面速度分析

閘門小開口（a=0.2 和s=0.8，x=-5，-2，2，5，10）時的流速分布，如圖3所示。

圖3 閘門小開口時的流速分布

從圖3可以看出，與室內(nèi)試驗結(jié)果對比，RSM 和RNG k-ε湍流模型能夠很好地模擬出水流實際速度分布情況，而標(biāo)準(zhǔn)RNG k-ω 模型通常適用于射流模擬，在收縮射流上的模擬效果最差。

閘門大開口（a=0.7 和s=0.9，x=-1，-0.5，1，2，3）時的流速分布，如圖4所示。

圖4 閘門大開口時的流速分布

從圖4 可以看出，在閘門大開度下，除閘門和收縮段之間存在射流外，所有RANS 模擬結(jié)果與試驗一致。RSM 和RNG k-ε 湍流模型的結(jié)果非常相似，而RNG k-ω 模型計算出的結(jié)果與試驗有很大偏差。在閘門上游，由于存在二維模擬未模擬出的再循環(huán)區(qū)，所以模擬出的縱向速度略低。模擬結(jié)果表明，靠近收縮截面的混合區(qū)域的流速略有增加，這與僅在x=1 處的試驗相同，其誤差可能是由于湍流模型和橫向流速分布不均勻造成的。盡管如此，從試驗和數(shù)值方法推導(dǎo)出的收縮系數(shù)非常接近，對于所有淹沒比而言，小開口閘門的水流可以認為是射流，而大開口閘門下的流速分布在收縮斷面中更不均勻。

3.2 能量修正系數(shù)分析

通過對流速的估計結(jié)果可以計算動能和動量系數(shù)。淹沒條件s=0.9 時能量修正系數(shù)α變化規(guī)律，如圖5所示。

圖5 s=0.9時能量修正系數(shù)α變化規(guī)律

從圖5可以看出，α的最小值對應(yīng)于縱向流速分布（邊界層除外）最均勻處，同時與閘門開度有明顯相關(guān)性，其大約為閘門開度的2～4倍。在下游部分，收縮段的α值在1.02～1.03，動量系數(shù)的平均值為1.02（閘門上游）和1.015（閘門下游）。上游水池和收縮段之間的水頭損失約為總水頭的1%，主要是由于邊界層中的高流速梯度引起的。

4 結(jié)論

本文為分析不同水閘開度對水流特性的影響，采用Fluent 數(shù)值模擬方法，引入能量耗散和流動阻力，分析對比了不同閘門開度下RNG k-ε、RNG k-ω和RSM 湍流模型對流體速度分布和能量系數(shù)的影響，并通過室內(nèi)水閘試驗進行驗證對比，結(jié)論如下。

（1）當(dāng)閘門為小開口時，RSM 和RNG k-ε 湍流模型能夠很好地模擬出水流實際速度分布情況，而標(biāo)準(zhǔn)RNG k-ω 模型通常適用于射流模擬，在收縮射流上的模擬效果最差。

（2）在閘門大開度下，除閘門和收縮段之間存在射流外，所有RANS 模擬結(jié)果與試驗一致。RSM 和RNG k-ε 湍流模型的結(jié)果非常相似，而RNG k-ω 模型計算出的結(jié)果與試驗有很大偏差。

（3）對于所有淹沒比，小開口閘門的水流可以認為是射流，而大開口閘門下的流速分布在收縮斷面中更不均勻。