過堰水流數(shù)值模擬研究

□康迎賓 □姚園園

（華北水利水電大學(xué)）

過堰水流數(shù)值模擬研究

□康迎賓 □姚園園

（華北水利水電大學(xué)）

文章運用FLUENT軟件，以某城市景觀河道的水流為研究對象，對小型實用堰的過堰水流進(jìn)行了數(shù)值模擬研究。對過堰水流進(jìn)行了數(shù)值模擬，采用FLUENT計算軟件中的RNGk-ε模型（即重整化群k-ε模型）建立了三維水-氣兩相流模型，由于過堰水流水氣交界面的波動較大，故選擇VOF方法對水流的水氣交界面進(jìn)行跟蹤模擬?？刂品匠痰碾x散方法則采用有限體積法，并采用迎風(fēng)差分格式進(jìn)行空間離散。

數(shù)值模擬；過堰水流；有限體積法；VOF方法

1 引言

隨著水生態(tài)文明城市的創(chuàng)建，城市河道規(guī)劃及城市水系景觀設(shè)計等在民生水利建設(shè)工程中越來越多地受到人們的的關(guān)注。在設(shè)計景觀河道時，鑒于對河道的防洪、排澇、景觀以及生態(tài)等功能的考慮，經(jīng)常采用復(fù)式斷面作為河道整治斷面形式的最佳選擇。但是降雨的時空分布不均，水文的洪枯流量大小直接影響著河道的水景觀。為保持河道的水景觀，經(jīng)常在河道上修建小型實用堰。文章針對這種小型景觀河道工程受試驗經(jīng)費與試驗周期等限制，大多只采用水力計算的問題，提出了采用數(shù)值模擬代替物理模型試驗的方法對水力設(shè)計方案加以驗證。

2 紊流模型控制方程

文章模擬所采用的RNG k-ε模型的控制方程見公式（1）與（2）。

式中：ρ表示密度；t為時間；k為紊動能；ε為耗散率；α、β均為修正系數(shù)；u為速度矢量；ui為速度分量；xi為坐標(biāo)分量；μ為動力粘滯系數(shù)；μt為動力粘度；Gk為紊動能的產(chǎn)生項；C1ε、C2ε為經(jīng)驗常數(shù)；

3 計算模型及網(wǎng)格劃分

3.1 計算模型區(qū)域

算例1：在模型計算區(qū)域中，取l=5 m長的堰前段，堰高P1=1 m，進(jìn)口形狀為90°直角進(jìn)口，堰厚δ=3 m，槽深H=3.50 m，計算區(qū)域總長L=18 m。水流的流態(tài)為自由出流，堰前水深h與進(jìn)口與流速uint詳見表1。

算例2：在算例2中，選取距離折線型堰l=5 m處的斷面作為進(jìn)口斷面，堰高P1=1 m，堰厚δ=3 m，折線型堰的下游坡a:b=1:1，槽深H=3.50 m，計算區(qū)域總長不變，取L=18 m。水流的流態(tài)也為自由出流，堰前水深h與流速uint見表1。

表1 各算例中的堰前水深表

3.2 網(wǎng)格劃分

對計算區(qū)域中的模型進(jìn)行網(wǎng)格劃分，網(wǎng)格類型采用非均勻結(jié)構(gòu)化（Map）網(wǎng)格，網(wǎng)格單元采用三維的六面體（Hex）。網(wǎng)格劃分單元個數(shù)、面（線）數(shù)與節(jié)點個數(shù)如表2所示。

表2 模型網(wǎng)格劃分信息表

在對網(wǎng)格進(jìn)行劃分時，為避免網(wǎng)格間距的大小對網(wǎng)格質(zhì)量的影響，本章在對算例1與算例2中的模型進(jìn)行網(wǎng)格劃分時，均采用同一尺寸大小的網(wǎng)格間距。

4 指定邊界條件

文章采用VOF法對自由表面進(jìn)行處理，所以將進(jìn)口邊界分成兩部分，即空氣進(jìn)口與水流進(jìn)口。其中，位于水面之上的區(qū)域指定為空氣進(jìn)口邊界，對其進(jìn)行壓力邊界條件的設(shè)置，設(shè)定值為1個工程大氣壓；位于水面以下的區(qū)域指定為水流進(jìn)口邊界，設(shè)定為流速進(jìn)口邊界條件。自由表面指定為壓力邊界，這是由于該邊界處的壓力值已知，但是，未知進(jìn)出自由表面的空氣通量，故將此設(shè)定為壓力邊界條件。計算區(qū)域中的水流為自由出流，因此設(shè)定出口為壓力出口邊界條件，其出口壓力參數(shù)值設(shè)置為大氣壓強值。邊壁面設(shè)置為無滑移邊界，在近壁粘性流動區(qū)，則采用標(biāo)準(zhǔn)的壁函數(shù)法進(jìn)行設(shè)置。其余邊界默認(rèn)為固體壁面。

5 求解參數(shù)設(shè)定

5.1 求解器參數(shù)設(shè)定

文章模擬涉及到多相流計算。選擇Volume Fraction項，在Number of phase中選擇水-氣兩相流動，其他保留默認(rèn)設(shè)置；選擇本次模擬所要采用的紊流模型k-epsilon；對模型的材料進(jìn)行定義，設(shè)置水的密度為1 000 kg/m3，動力粘度為0.001 Pa·s，并設(shè)置基本相與第二相。然后，設(shè)置操作環(huán)境為大氣壓強101 325 Pa，重力加速度g為9.81 m2/s。

在對邊界條件中紊流參數(shù)進(jìn)行設(shè)置時，其紊動參數(shù)k與耗散率ε可按公式（3）與（4）計算。

式中：uint是水流進(jìn)口流速；Hint為進(jìn)口平均水深。

5.2 求解控制參數(shù)設(shè)定

壓力-速度耦合選擇PISO算法，PISO算法在計算時能夠進(jìn)行相鄰校正與偏斜校正，其Iteration項采用默認(rèn)值1。

文章在流場初始化對話框中的Compute From項中選擇velocity-inlet；在迭代計算設(shè)置對話框中，Time Step Size項設(shè)置為0.10 s，在需要求解的時間步數(shù)項設(shè)置為1 000，其他設(shè)置項為默認(rèn)值。

6 求解結(jié)果分析

文章模擬寬頂堰在枯水時期的水流流態(tài)及水力特性，計算工況分兩種。分別對相同流量下的寬頂堰的過堰水流進(jìn)行了三維紊流數(shù)值模擬。計算結(jié)果均采用進(jìn)出口流量差≤0.50%時，認(rèn)為計算完成。

水流在未流向?qū)嵱醚邥r，流速幾乎處于均勻狀態(tài)。當(dāng)水流經(jīng)過寬頂堰或折線型堰時，在進(jìn)口前端處，水流過水?dāng)嗝骈_始發(fā)生收縮，水流的流速也相應(yīng)地增大。當(dāng)水流流過堰頂后，從流速矢量圖中可以觀察到堰后水流的最大流速均處在水面以下范圍內(nèi)。

水流在流經(jīng)寬頂堰后，會在堰后下游面處產(chǎn)生順時針方向的分離漩渦，并且不同的流量下的分離漩渦的尺度大小也不一樣，流量增大時，水流流速值增大，分離漩渦的尺度也增大。該分離漩渦的產(chǎn)生，主要是由于堰體上游水流流速與堰后水流流速平行接觸后在堰體尾部附近，水流發(fā)生分離，引起剪切層流動，進(jìn)而又由于剪切層流動的不穩(wěn)定，使剪切層在不同水流流場中自發(fā)地產(chǎn)生大小尺度不同的分離漩渦。該分離漩渦內(nèi)有關(guān)紊動強度與時均流速梯度都非常大，并且大量的紊動能將耗散在漩渦處。

還可分析得出分離漩渦處的流速分布特點。流速在漩渦內(nèi)部的中心偏右上角范圍內(nèi)，出現(xiàn)流速較小值；而在分離漩渦的外邊緣，流速較小值的出現(xiàn)往往會在堰后腳處，并且隨著堰后分離漩渦尺寸的減小而減小。分離漩渦處流速的分布趨勢為由外向內(nèi)，流速在減小。

水流經(jīng)過寬頂堰時，水流會產(chǎn)生分離漩渦，而水流經(jīng)過折線型堰時在堰下游面處不產(chǎn)生漩渦，并且由圖中折線型堰下游坡面處，速度等值線較密集，可得出水流流經(jīng)折線型堰時，水流受下游坡邊壁影響。另外，下游坡頂角處的速度等值線的斜向上彎曲度會隨著流量Q值的增大，范圍也相應(yīng)地加大。

文章采用RNG k-ε紊流模型，有限體積法與VOF方法，分別模擬了相同流量下，直角進(jìn)口寬頂堰與折線型堰在自由出流時的水流流態(tài)與水力特性。對寬頂堰的過堰水流與折線型堰的過堰水流進(jìn)行對分析，得出水流在經(jīng)過這兩種堰體時，水面都會發(fā)生兩次下降。但水流在經(jīng)過寬頂堰時，在堰后會產(chǎn)生順時針方向的分離漩渦，而水流經(jīng)過折線型堰時，則在堰后不存在分離式漩渦，并且水流緊貼堰的下游坡流出。

通過對直角進(jìn)口的寬頂堰與折線型堰的過堰水流的三維數(shù)值模擬結(jié)果分析，得出同等條件下折線型堰的過堰水流要比寬頂堰的過堰水流的流態(tài)較好，考慮選取折線型實用堰，與實際方案堰型選擇類型相符。其研究成果可為同類小型河道工程設(shè)計，提供一定的參考價值和借鑒作用。

編輯：劉長垠 韋詩佳

TV132

A

1673-8853（2017）03-0072-02

2017-1-5