王 鵬

（云南創(chuàng)研勘測(cè)設(shè)計(jì)研究院有限公司，云南 昆明 650051）

跌水是常見的渠道工程落差建筑物，在排洪、灌溉、排水等工程中得到了廣泛應(yīng)用[1]。渠道通過坡度較陡地帶時(shí)，為減少深挖方或高填方，通常在落差集中處修建跌水，作為連接渠道上下游的設(shè)施[2，3]。水流流經(jīng)跌水時(shí)由緩流過渡到急流，在重力作用下直接沖擊下游渠道底板，易造成渠道破壞。因此，跌水建筑物一直是水利工作者研究的熱點(diǎn)。

韓冬梅等[4]采用經(jīng)驗(yàn)公式法計(jì)算對(duì)比了跌水、底流2 種消能方式的消能效果，指出在跌水設(shè)計(jì)時(shí)應(yīng)根據(jù)地形條件合理選擇跌水高度及級(jí)數(shù)；譚志偉[5]采用經(jīng)驗(yàn)公式法對(duì)豎井式跌水的進(jìn)口溢流堰、橫洞進(jìn)行了水力計(jì)算，提出了解決橫洞內(nèi)水流流態(tài)不穩(wěn)定、消能效果不佳等問題的措施；黃朝煊[6]基于水力學(xué)理論，采用MATLAB 軟件對(duì)多級(jí)跌水消能進(jìn)行了分析，為多級(jí)跌水消能優(yōu)化設(shè)計(jì)提供了支撐；羅立群等[7]通過1∶20單體水工模型試驗(yàn)觀測(cè)了最大分洪流量下不同體型跌水的流態(tài)及流速，并通過能量方程計(jì)算了跌水的消能率；丁強(qiáng)等[8]采用計(jì)算流體力學(xué)軟件CFX 研究了頂部封閉條件下流量對(duì)跌水窨井內(nèi)流態(tài)的影響；李明達(dá)等[9]采用Fluent軟件結(jié)合流體體積分?jǐn)?shù)法（Volume of Fluid，VOF）與標(biāo)準(zhǔn)k-ε模型對(duì)虹吸井跌水摻氣問題進(jìn)行了數(shù)值模擬，分析了不同流量、不同水位和不同跌水高度時(shí)虹吸井下游的摻氣濃度變化情況；張輝等[10]采用Flow 3D軟件結(jié)合RNGk-ε 紊流模型對(duì)引江濟(jì)淮工程韓橋段跌水的水流壓強(qiáng)及流態(tài)進(jìn)行了數(shù)值模擬，并根據(jù)流速分布情況對(duì)跌水設(shè)計(jì)方案進(jìn)行了優(yōu)化。

本文以某河道多級(jí)跌水為研究對(duì)象，采用計(jì)算流體力學(xué)軟件對(duì)多級(jí)跌水水流流態(tài)、流速、消能率等水力特性進(jìn)行了比較全面的分析研究，可為工程設(shè)計(jì)提供依據(jù)。

1 工程概況

某河道位于低山丘陵區(qū)，總長(zhǎng)4 605.88 m。為提高區(qū)域整體防洪減災(zāi)能力，擬實(shí)施流域水環(huán)境綜合整治工程?？紤]到部分河段底坡較陡，為降低工程挖填方量，設(shè)計(jì)采用多級(jí)跌水連接河道上下游。跌水共2級(jí)，位于河道樁號(hào)2+0.00—2+68.00段，其中一級(jí)跌水位于樁號(hào)2+20.00處、二級(jí)跌水位于樁號(hào)2+42.00處。跌水高3.00 m，寬20.00 m，與河道等寬。為保護(hù)下游河道，在跌水下游設(shè)置消力坎，消力坎高0.5 m。

根據(jù)《防洪標(biāo)準(zhǔn)》（GB 50201-2014）[11]有關(guān)規(guī)定，工程防護(hù)區(qū)等級(jí)為Ⅳ等，防洪標(biāo)準(zhǔn)為30 a一遇，設(shè)計(jì)洪峰流量為216.00 m3/s，設(shè)計(jì)洪水位為743.28 m。

2 數(shù)值模型構(gòu)建

2.1 控制方程

流體運(yùn)動(dòng)控制方程[12]如下。

連續(xù)性方程為：

動(dòng)量方程為：

紊動(dòng)能k方程為：

耗散率ε方程為：

式中：xi，xj表示坐標(biāo)分量（m）；ui，uj表示速度分量（m/s）；i，j分別為1，2，3；t表示時(shí)間（s）；ρ表示密度（kg/m3）；Si表示動(dòng)量方程的廣義源項(xiàng)；μ表示流體黏性系數(shù)（Pa·s）；p表示壓力（Pa）；k表示紊流脈動(dòng)動(dòng)能；ε表示耗散率；μt表示紊流黏性系數(shù)（Pa·s）；常數(shù)σk=1.0，σε= 1.3，C1ε=1.44，C2ε=1.92；GK表示平均速度梯度引起的紊動(dòng)能產(chǎn)生項(xiàng)。

采用有限體積方法[13]離散上述基本方程，以VOF法[14]追蹤水流自由表面。鑒于RNGk-ε 紊流模型模擬復(fù)雜水流流動(dòng)的優(yōu)越性[15]，本文采用RNGk-ε模型封閉基本方程組，對(duì)多級(jí)跌水水力特性進(jìn)行三維數(shù)值模擬研究。

2.2 網(wǎng)格及邊界條件

網(wǎng)格是模型計(jì)算進(jìn)行離散的基礎(chǔ)，是影響數(shù)值模擬精度的關(guān)鍵因素。本文采用AutoCAD 軟件建立河道及跌水的三維模型，采用結(jié)構(gòu)化網(wǎng)格對(duì)模型進(jìn)行網(wǎng)格劃分。建立三維模型時(shí)，本著與實(shí)際工程一致的原則[16]，根據(jù)河道及跌水的實(shí)際尺寸進(jìn)行建模，模型包括河道、多級(jí)跌水、消力坎等。河道及跌水三維模型，如圖1所示。

圖1 河道及跌水三維模型

經(jīng)多次試算，確定結(jié)構(gòu)化網(wǎng)格尺寸為Δx=Δy=Δz=0.10 m，其中順?biāo)鞣较驗(yàn)閤方向、重力方向?yàn)閦方向，總共960 萬網(wǎng)格。為避免計(jì)算中發(fā)生水流封頂現(xiàn)象影響模擬精度，河道最大模擬高程為743.80 m，高于設(shè)計(jì)洪水位743.28 m。

模型入口采用壓力入口，通過設(shè)計(jì)洪水位給定，水的體積分?jǐn)?shù)設(shè)為1；出口設(shè)置在第二道消力坎下游20.00 m 處，以確保出口水流平順，出口邊界條件設(shè)置為自由出流；河道頂部給定為壓力入口邊界條件，水的體積分?jǐn)?shù)設(shè)為0；底部及兩側(cè)邊壁采用固體壁面邊界條件。

3 數(shù)值模擬結(jié)果分析

3.1 水流流態(tài)

設(shè)計(jì)洪水位下河道及跌水水流流態(tài)，如圖2所示。

圖2 河道及跌水水流流態(tài)

從圖2 可以看出，一級(jí)跌水上游段水流流態(tài)比較平順，在一級(jí)跌水處水流由緩流轉(zhuǎn)變?yōu)榧绷?，之后快速?zèng)_泄向第一道消力坎，在消力坎阻擋下產(chǎn)生明顯的水躍。水躍越過第一道消力坎后經(jīng)二級(jí)跌水繼續(xù)向下游流動(dòng)，在第二道消力坎前劇烈翻滾，急速消散能量。越過第二道消力坎后，水流逐漸趨于平順。

值得注意的是，由于跌水在平面上呈弧形，因此河道內(nèi)水流流程不同，靠近兩側(cè)岸坡的水流流程較短，自跌水處跌落時(shí)間早于弧頂處水流。在一級(jí)跌水下游，兩側(cè)水流跌落后產(chǎn)生朝向?qū)Π兜臎_擊波，沖擊波在河道中間位置交匯后形成較為明顯的波峰疊加效應(yīng)，疊加后水深約2.82 m。數(shù)值模擬很好地模擬了沖擊波及水流疊加現(xiàn)象。

3.2 水流流速

沿河道中心線位置，繪制順?biāo)鞣较蛄魉偾€，如圖3所示。

圖3 河道沿程流速變化曲線

從圖3 可以看出，在一級(jí)跌水上游（樁號(hào)2+0—20.00），隨著水深逐漸減小，水流流速逐漸增大；在二級(jí)跌水上游（樁號(hào)2+20.00—40.90），水流流速快速上升后趨于穩(wěn)定；在第二道消力坎前后，由于水流波峰疊加效應(yīng)，消力坎頂部水深稍高，流速出現(xiàn)下降。經(jīng)過二級(jí)跌水后，水流流速快速增大，至樁號(hào)2+58.44 處，達(dá)到最大值約11.94 m/s；水流經(jīng)過第二道消力坎后流速趨于平穩(wěn)，平均流速約6.19 m/s。數(shù)值模擬很好地呈現(xiàn)出河道水流在重力、跌水及消力坎作用下流速的變化情況。

3.3 消能率

消能率η是反映跌水消能效果的指標(biāo)[17]，其定義為：

式中：E1，E2分別為跌水上，下游斷面總水頭（m）；Z1，Z2分別為跌水上，下游水面相對(duì)于基準(zhǔn)面的落差（m）；v1，v2分別為跌水上，下游斷面平均流速（m/s）；α1，α2分別為對(duì)應(yīng)的流速系數(shù)。

經(jīng)計(jì)算，設(shè)計(jì)洪水位下多級(jí)跌水消能率達(dá)62.83%，有效消減了水流能量。

4 結(jié)論

采用計(jì)算流體力學(xué)軟件對(duì)河道多級(jí)跌水水流流態(tài)、流速等進(jìn)行了數(shù)值模擬，并計(jì)算得到了多級(jí)跌水消能率。一級(jí)跌水上游水流流態(tài)較平穩(wěn)，流速呈逐漸增大趨勢(shì)；越過一級(jí)跌水后水流翻滾劇烈，流速相應(yīng)地波動(dòng)較大，多級(jí)跌水消能效果較好。數(shù)值模擬能很好地呈現(xiàn)出多級(jí)跌水水流流態(tài)、流速等，可用于工程設(shè)計(jì)方案分析驗(yàn)證。由于項(xiàng)目時(shí)間緊迫，未進(jìn)行水工模型試驗(yàn)研究，后繼將采用模型試驗(yàn)進(jìn)行深入研究。