明渠急流分流結構數值模擬

燕春林 何孟狄 余 振 汪德金

（中國市政工程西南設計研究總院 四川 成都 610082）

河道或渠道采取分流措施已廣泛的應用于各種水利工程中，諸如河道防洪、城市防洪、農業灌溉等方面。明渠分流常見布置是在明渠側墻部位設置分流渠道，具體體型則由實際效用決定。在明渠分流處，存在有一定的流線分離，從而出現局部回流區和束窄現象。Neary and Odgaard研究了低佛汝德數時非等寬明渠直交分流區的流場結構、分流能力和分流點處的水深。A.S.Ramamurthy等[2]和Jianchun Huang等[3]采用數值模擬研究了分流渠和主渠交角為90°時的分流區的流場結構。Bahareh Pirzadeh,Hamid Shamloo采用二維和三維數值模擬相結合的手段，模擬了分流明渠的水流結構。Amruthur S.Ramamurthy[5]等采用試驗研究了分流明渠流態和分流渠進口上游主渠進佛汝德數對分流量的影響。

以上研究的分流明渠底板坡度變化較小，分流渠和主渠坡度基本相同，僅對不同分流渠角度的流場和分流量等水力特性指標進行了詳細研究。本文以分流渠與主渠交角為40°，且分流渠底板為急流底坡的三維異型分流結構。

研究對象，并將數值模擬結果與模型試驗結果相比較以驗證數值模擬的有效性。

1 計算模型及邊界條件

1.1 數學模型和計算方法

本文采用引入具有自由面追蹤功能的VOF法和標準kε雙方程模型來對中閘室出口空腔進行數值模擬，其基本方程為：

式中：ρ和μ分別為密度和粘性系數，其表達式為：ρ＝αwρw＋（1－αw）ρa，μ＝αwμw＋（1－αw）μa，其中，αw為網格上水的體積分數，0≤αw≤1，ρw、μw和ρa、μa分別指水和空氣的密度、粘性系數；p為修正的壓力；μt為紊流粘性系數，可由紊動能k和紊動耗散率ε求出其中，Cμ為經驗常數，fi為流體受到的慣性力，本文僅有重力。σk和σε分別為k和ε的紊流普朗特數，以上各式中的常數取值見表1。

表1 控制方程中的常數值Table1 The constant value of control equation

方程離散采用有限體積法，對流項采用二階迎風格式，擴散項采用中心差分格式，壓力速度耦合采用PISO算法。每個網格上的體積分數變化采用求解下式完成：

其中，t為時間，ui和 xi分別為速度分量和坐標分量。

1.2 邊界條件

數值模擬體型如圖1所示，由主管道、分流結構、分流管道和主渠收縮段組成。分流結構分為分流中墩、1#分流道和2#分流道分流管道。分流結構進口處設分流墩，其前沿導向角（前沿尖角射線與主渠軸線的夾角）為40°，分流墩段底板為三維異型結構，由多段弧組成?；緟禐椋褐髑罊M斷面為3.1×2.9m （寬×高），長33m；1#分流管道進口橫斷面為5.42×4m （寬×高）；2#分流管道進口橫斷面為5.17×4m （寬×高）；主管道束窄處斷面為1.3×4m（寬×高）；分流管下游主管道典型橫斷面為4.0×4.0m（寬×高），長41.5m。該結構的主要特點是：分流結構進口不是簡單設置于主渠側墻處，而是將進口設置于主渠內并與主渠收縮段相結合，分流渠底板坡降較大，使水流呈急流形式，以達到增加分流能力的目的。網格劃分采用全局非結構化網格，最小計算尺寸為0.1m。進口采用速度邊界條件，空氣進口和出口均采用壓力邊界條件，固壁采用無滑移邊界條件，粘性底層采用標準壁面函數法處理。

1.3 計算工況

根據模型試驗結果，并與模型試驗相對應，其數值模擬的計算工況如表2所示。

圖1 計算體型Fig.1 Numerical modle

表2 數值模擬工況Table2 Performance of numerical simulation

2 成果分析

2.1 流態及流場分布

數值模擬結果表明：各計算工況下，該體型參數均能滿足泄流能力，管道內流態基本穩定，隨著流量增加，1#和2#分流道內部沿側壁有翻轉水流，分流墩下游局部有較大水躍，紊流程度劇烈，主渠收縮段后有水力坡降會陡增。當流量小于20m3/s時，管道內水流平順，分流墩后主渠收縮段對水流束窄能力較小，水流進入分流管道彎段以后，受彎道離心力作用而向外偏轉，在水流進入直段前以弱水躍銜接。當流量大于30m3/s時，分流墩后主渠收縮段對水流束窄能力增強，分流管道流量相應增加，導致分流段流速增加，分流管道偏轉水流和分流墩下游水躍強度均增加。當流量大于65m3/s時，分流墩后主渠收縮段對水流束窄能力明顯增強，水流在分流段段進口形成回流和平面旋流，但受下游水位影響，分流墩后呈淹沒水躍（圖2～圖4）。

圖2 分流結構流態（Q=10m3/s）Fig.2 The flow pattern of diversion structure(Q=10m3/s)

圖3 分流結構流態（Q=30m3/s）Fig.3 The flow pattern of diversion structure(Q=30m3/s)

圖4 分流結構流態（Q=50m3/s）Fig.4 The flow pattern of diversion structure(Q=50m3/s)

2.2 水面線分布

圖5 計算結果水面線（主渠）Fig.5 Water surface profile of numerical simulation

水面線分布是衡量分流結構體型是否合理的一個重要水力特性指標，通過與模型試驗測量水面線相比較，從而驗證數值模擬成果。在流量小于20m3/s時，主渠收縮段對水流的束窄能力有限，主渠內水面線均平順；流量大于30m3/s后，主渠收縮段對水流的束窄能力逐漸增強，束窄斷后水力坡降也逐漸增大。對于分流結構段，在流量Q≤10m3/s時，水面線分布均勻，當流量Q≥20m3/s后，因流速增大并受離心力影響而向左側墻偏轉，流道內橫斷面水面線分布不均勻。在流量Q≤5m3/s時，分流墩后以弱水躍銜接，當流量Q≥20m3/s后，分流墩后水躍強度逐漸增大，當流量Q≥50m3/s后，隨著下游水位升高，水躍會呈淹沒狀態（如圖5、圖6）。數值模擬結果與模型試驗結果對比表明，二者水面線結構符合良好（如圖7、圖8）

圖6 計算結果水面線(分流渠)Fig.6 Water surface profile of numerical simulation

圖7 計算值與模型值對比（主渠道、工況6）Fig.7 Comparison of stream line between numerical and measured(main channel,performance 6)

圖8 計算值與模型值對比（分流管道、工況6）Fig.8 Comparison of stream line between numerical and measured(diversion channel,performance 6)

2.3 分流比

分流比是衡量分流結構分流能力大小的一個重要指標。比較分析表明，數值模擬成果與模型試驗成果符合良好（如表3）。表3說明，各工況下，分流結構分流后兩渠道流量比值在1.8∶1～3.3∶1之間，說明分流結構具有較強的分流能力，該體型能滿足分流要求。

表3 各工況下的分流比（分流渠道流量：分流后主渠道流量）Table 3 The split ration under different operation condition(discharge of main channel:discharge of diversion channel)

3 結論

急流分流結構是一種較為新型的分流結構，其分流能力和分流道體型是影響水力學指標的關鍵因素。本文通過數值模擬模型試驗研究急流分流結構并與模型試驗成果相比較，二者在流場、水面線和分流比方面均具有較好的符合性，說明數值模擬結果有效。該急流分流結構能分走主渠2/3以上的流量，說明具有較強的分流能力。

［1］Neary,V.S.and Odgaard,A.J.,(1993).Three-Dimensi-onal Flow Structure at Open Channel Diversions[J].Hydr.Engrg.,ASCE,119(11)∶1223-1230

［2］A.S.Ramamurthy,Junying Qu,Diep Vo,and Chao Zhai.3-D Simulation of Dividing Flows in 90°Rectangular Closed Conduits[J].J. Fluids Eng.September 2006,128(5)∶1126-1129.

［3］Jianchun Huang,Larry J.Weber,and Yong G.Lai.Three-Dimensional Numerical Study of Flows in Ope-n-Channel Junctions[J]. J.Hydr.Engrg.2002,128(3),pp.268-280

［4］Amruthur S.Ramamurthy,Duc Minh Tran,Luis B.Carballada. Dividing Flow in Open Channels[J].Journal of Hydraulic Engineering, Vol.116,No.3,March 1990,pp.449-455.

［5］Bahareh Pirzadeh,Hamid Shamloo(2007).Numerical investigation of Velocity Field in Dividing Open-Channel Flow [A].12th WSEAS Int. Conf.on APPLIED MATHEMATICS,Cairo,Egypt,December 29-31,pp. 194-198.