廣彗冰，王 健*，李澤森，曹博召，劉 超，王康宏

?灌溉水源與輸配水系統(tǒng)?

淤積工況下薄壁三角堰流量特征研究

廣彗冰1，王 健1*，李澤森1，曹博召1，劉 超1，王康宏2

（1.西北農(nóng)林科技大學，陜西 楊凌 712100；2.韓城市水土保持工作站，陜西 韓城 714000）

【目的】探究淤積工況下薄壁三角堰的流量特征?！痉椒ā坷糜嬎懔黧w力學軟件FLUENT 19.0，結(jié)合室內(nèi)試驗，模擬薄壁三角堰渠道從無淤積至淤滿過程中9個淤積工況下的過堰流場，分析流場內(nèi)的流量、流速和水位變化，探討渠道流量特征隨淤積厚度的變化情況。【結(jié)果】模擬流量與理論流量、實測流量數(shù)值吻合，渠道內(nèi)流速分布符合流速分布一般規(guī)律，證明模擬結(jié)果準確可靠；當理論流量為0.221 0 L/s時，隨著堰前淤積的加厚，模擬流量與實測流量呈無規(guī)律的波動變化，波動范圍為理論流量的-0.45%～8.60%；當理論流量為0.221 0 L/s時，對比無淤積和淤積工況下的堰前水面曲線，模擬與實測曲線均在近堰板處出現(xiàn)不同程度的的水位壅高，高度隨淤積厚度的增加而波動上升?！窘Y(jié)論】淤積會使堰前水位壅高且增大薄壁三角堰過流能力，而堰流公式中缺少適應(yīng)這種變化的修正參數(shù)，導(dǎo)致生產(chǎn)中所得理論流量偏小。

薄壁三角堰；淤積；計算流體力學；流量；流速

0 引言

【研究意義】在水資源短缺背景下，精確計量并合理調(diào)配灌區(qū)水資源是提高用水效率的重要手段之一。目前，我國灌區(qū)量水設(shè)備已多達百種。其中，薄壁三角堰能將微小流量變化通過顯著水頭變化體現(xiàn)出來，再結(jié)合堰流公式便可得到較精準的流量數(shù)據(jù)，從而被廣泛應(yīng)用于溝渠流量監(jiān)測當中。同時，三角堰也是多種組合堰中的基礎(chǔ)堰型，分析薄壁三角堰對于改進和發(fā)展復(fù)雜堰型有重要的指導(dǎo)意義。而在實際應(yīng)用中，薄壁三角堰的出現(xiàn)會改變溝渠當中的水流特征，在水土流失區(qū)，堰前易產(chǎn)生泥沙淤積；隨著使用年限的增長，淤積不斷加深，堰頂高度不斷縮短，逐漸偏離堰流公式的測流條件。

【研究進展】渠道泥沙淤積是灌區(qū)存在的普遍問題[1]，許多學者就此開展了廣泛的研究。在無壩引水明渠，徐霖玉等[2]通過物理模擬試驗研究了泥沙淤積的原因，主要為上游來沙量大、渠道底坡較小、缺少有效的排沙系統(tǒng)。高紅艷等[3]在鐙口揚水灌區(qū)展開了研究，認為渠道淤積受地理環(huán)境影響，糙率大、斷面不合理等會加重淤積。在引黃灌區(qū)，泥沙淤積問題更是屢見不鮮，備受關(guān)注[4-6]。同時，Maha[7]研究發(fā)現(xiàn)，隨著沉積深度的增加，流動特征D隨之增加，流量系數(shù)較水平河床增加30%～46%?？傮w來說，渠道淤積在灌區(qū)普遍存在，成因多樣，且會在一定程度上影響測流精度。

堰流研究手段通常有公式推算法和直接測定法，但由于理論計算具有諸多限制條件，野外流速測量存在較多風險與困難，計算流體力學隨著科技發(fā)展越來越受到青睞。許多研究者采用數(shù)值模擬[8-10]來研究堰流，證明在FLUENT軟件中采用合適的湍流模型與氣液兩相流模型能夠較準確的解出泄流流場的各種水力參數(shù)，是一種行之有效的研究方法。柳雙環(huán)等[11]利用RNG k-ε模型分析紊流，基于VOF法追蹤自由表面，對小型U形渠道三角剖面堰進行了數(shù)值模擬，對比發(fā)現(xiàn)流量、水位模擬值與實際值相對誤差較小，證實所用方法可以快速可靠地模擬水流狀況。徐君玲等[12]采用帶自由表面的k-ε模型，并使用自定義UDF函數(shù)處理入口條件，對整個流場進行數(shù)值模擬計算，結(jié)果表明模擬結(jié)果與實際數(shù)據(jù)吻合良好，說明所選模型能夠準確模擬矩形薄壁量水堰的流場特性。魏文禮等[13]通過對比理論推導(dǎo)所得的梯型薄壁堰堰流公式與利用RNG k-ε模型和VOF模型解出的流量計算式，發(fā)現(xiàn)二者流量值與堰流特性結(jié)果基本一致，證明模擬方法準確可靠。此外，數(shù)值模擬方法應(yīng)用于其他模型[14-16]的情況均為本研究數(shù)值模擬模型的選取積累了經(jīng)驗，并提供了可行性保障。

【切入點】當前，大量研究集中在渠道淤積的成因、減淤措施和無淤積時堰流特征的分析，而針對淤積狀態(tài)下薄壁三角堰測流問題的研究較少?！緮M解決的關(guān)鍵問題】因此，本研究采用數(shù)值模擬與室內(nèi)試驗相結(jié)合的方法，對薄壁三角堰渠道從無淤積至淤滿過程中9個淤積工況下的過堰流場進行模擬，定量分析淤積對流量觀測產(chǎn)生的影響，以及淤積下堰流特征的變化，為灌溉渠系和河溝流量精確測定提供依據(jù)，促進灌區(qū)水資源管理科學化。

1 試驗設(shè)計

1.1 試驗方法

試驗在西北農(nóng)林科技大學水工實驗室內(nèi)進行，利用特制水槽進行流量試驗。具體裝置由水泵、DK800-6玻璃轉(zhuǎn)子流量計、水槽與水流收集箱構(gòu)成，如圖1所示。水槽高0.3 m，長0.9 m，由薄鋼板制成，分為平水柵、輸水渠道以及薄壁三角堰3部分。其中，平水柵長0.1 m，內(nèi)接進水管；輸水渠道長0.8 m，內(nèi)部以堰板為起點，向上游每隔0.05 m劃分一個過水斷面，共計11個，分別命名為1～11號，如圖2所示。

圖1 試驗裝置示意

圖2 過水斷面標號示意

在現(xiàn)實生產(chǎn)中，薄壁三角堰的尺寸根據(jù)當?shù)厍赖膶捝钋闆r進行適應(yīng)性設(shè)計，其中直角堰的應(yīng)用最為廣泛。因此，試驗設(shè)計堰板尺寸為高0.3 m，寬0.3 m，厚度2 mm，缺口形狀為一個斜邊邊長0.2 m的等腰直角三角形，居堰板上邊緣中央。渠道內(nèi)的淤積厚度受使用時長、當?shù)厮亮魇顩r等因素的影響，各渠道內(nèi)不等。因此，試驗設(shè)計渠道內(nèi)淤積變化過程為自渠道無淤積開始，淤積厚度以2.5 cm為步長逐級加深直至渠道淤滿，即劃分淤積厚度為0.0、2.5、5.0、7.5、10.0、12.5、15.0、17.5、20.0 cm共9個淤積工況，堰板尺寸與淤積工況劃分情況如圖3所示。

圖3 淤積劃分示意

試驗利用攪拌土和水泥定型成尺寸合適的塊體，浸潤后置于水槽中來模擬淤積。具體操作時先根據(jù)設(shè)定水深在槽中放水，當水深達到設(shè)定值且液面穩(wěn)定后，由玻璃轉(zhuǎn)子流量計讀取流量，精確為0.01 L/s；水深利用標尺為0.40 m的水位測針讀取，估讀至0.000 1 m，每個斷面測量3次，取平均值作為最終水深。試驗標號和過程如圖4、圖5所示。

圖4 試驗標號

圖5 試驗過程

1.2 數(shù)據(jù)處理

根據(jù)《水工建筑物與堰槽測流規(guī)范（SL 537—2011）》[17]（以下簡稱《規(guī)范》），當薄壁三角堰的堰口角在20°～100°之間時，根據(jù)水位計算過堰流量的公式為：

式中：為流量（m3/s）；D為流量系數(shù)，隨水頭稍有變化，可視為常數(shù)，本研究中堰頂高與渠道行進寬度的比值在0.08～0.67范圍內(nèi)，設(shè)計有效水深與堰頂高的比值在0.15～1.20范圍內(nèi)，通過查閱《規(guī)范》中流量系數(shù)圖可知，D值可取為0.60；為堰頂角（°）；g為重力加速度，取值為9.81 m/s2；e為有效水位（m）；為實測堰上水頭（m）；h為黏滯力和表面張力綜合影響的校正值，根據(jù)《規(guī)范》可知，在=90°時，h為0.000 85 m。

2 數(shù)值模擬

2.1 模型建立與網(wǎng)格劃分

使用ICEM 19.0構(gòu)建渠道的三維模型并劃分網(wǎng)格。軸、軸、軸分別表示展向、垂向、流向。三維模型主要由上游渠道、薄壁三角堰以及下游渠道組成。無淤積情況下的計算域為展向0.3 m，垂向0.4 m，流向1.4 m，建模結(jié)果如圖6所示。淤積厚度分別為2.5、5.0、7.5、10.0、12.5、15.0、17.5、20.0 cm時，計算域?qū)Ρ葻o淤積時上游渠道渠底抬高，即自坐標軸零點起減少展向0.3 m、流向0.9 m，垂向等于淤積厚度的矩形區(qū)域，淤積為10.0 cm時的建模結(jié)果如圖7所示。

圖6 無淤積情況下渠道的三維模型

圖7 淤積10.0 cm情況下渠道的三維模型

計算域的網(wǎng)格劃分采用以六面體為基本單元的結(jié)構(gòu)化網(wǎng)格[18]，對堰口的三角形區(qū)域進行Y型切分。驗證網(wǎng)格無關(guān)性時，在整個計算域使用統(tǒng)一的網(wǎng)格大小，依次選取0.005、0.010、0.015、0.020、0.025 m五種尺寸進行網(wǎng)格劃分，利用Patch功能測試網(wǎng)格大小對模擬結(jié)果的影響，采取水的體積分數(shù)為0.5時的結(jié)果作為水面線的位置，模擬結(jié)果如表1所示。

表1 各網(wǎng)格尺寸的模擬結(jié)果

由表1可知，水位模擬結(jié)果的精度隨網(wǎng)格尺寸的增加相對誤差逐漸增大，當網(wǎng)格尺寸為0.005 m時，相對誤差僅有0.67%。當相對誤差小于3%時，認為模擬誤差可忽略。由于薄壁三角堰主要利用V型缺口進行測流，需重點關(guān)注過堰水流的水面高度與流量變化情況，因此綜合考慮精度與時長，在計算域=0.00～0.30 m，=0.00～0.30 m，=0.00～1.05 m區(qū)域及堰板處，選取尺寸0.005 m進行網(wǎng)格劃分，而后網(wǎng)格尺寸以1.1的增長倍率向四周增大至0.010 m。無淤積情況下的網(wǎng)格總數(shù)約為152萬，劃分情況如圖8所示。

圖8 網(wǎng)格劃分

2.2 邊界條件設(shè)置

求解器中設(shè)置瞬態(tài)時間模型，操作條件設(shè)置中保持默認的101 325 Pa大氣壓強，再添加垂向重力加速度值9.81 m/s2。邊界條件設(shè)定中將空氣入口定義為只有空氣進入的壓強進口（PRESSURE INLET）；水流入口定義為只有水流進入的流速入口（VELOCITY INLET），隨著淤積的加深，增大入口流速以保證輸入流量基本穩(wěn)定，流速具體值根據(jù)堰流公式與明渠流公式求得，數(shù)值如表2所示。頂部是與大氣聯(lián)通的開放面，定義為只有空氣溢出的壓力出口（PRESSURE OUTLET）；渠道出口定義為無回流的無壓出流（PRESSURE OUTLET）。堰板、渠底與邊壁定義為無滑移的固體邊壁（WALL），采用標準壁面函數(shù)處理水流黏性底層[19]。

表2 流速分配情況

2.3 數(shù)學模型設(shè)定

本研究選用RNG k-ε模型分析流體介質(zhì)之間動量和能量的變化，RNG k-ε模型是在標準k-ε模型的基礎(chǔ)上對紊動能耗散率方程進行了修正，在其源項中加入了一個系數(shù)來反映主流時均應(yīng)變率，可以更好處理對流線彎曲程度較大的流動，其模型運輸方程詳見文獻[20]。

采用VOF多相流模型迭代追蹤水氣交界面[21]，并設(shè)定空氣與水的表面張力系數(shù)為0.072 N/m。VOF的基本思想是定義一個表示計算區(qū)域內(nèi)流體體積與計算區(qū)域總體積的相對比例的體積率函數(shù)W=(x,t)=1。對于某一具體計算單元而言，當W=(x,t)=1時，表示單元完全被液體充滿；當W=(x,t)=0表示一個空單元；若W=(x,t)=1∈(0,1)，則表示液體只充滿單元的一部分。描述W的控制方程詳見文獻[21]。

通過有限體積法對控制方程進行離散，使用SIMPLEC算法耦合壓力項和速度項，采用Body Force Weighted選項，對動量等的離散均采用二階迎風格式，由此獲得更高的計算精度[22]。流場初始化時，設(shè)定上游渠道在流動開始前已存在0.3×0.2×0.9（展向×垂向×流向）m3的水，其余流場空間充滿氣體。計算步長取0.005 s，設(shè)置進口處與堰口處的流量監(jiān)測，二者流量差小于3%時，認為計算收斂，水流呈穩(wěn)定狀態(tài)[12]。

3 結(jié)果與分析

3.1 模型驗證

設(shè)定上游堰前水深分別為0.030 0、0.032 0、0.034 0、0.036 0、0.038 0 m，運算求解得相應(yīng)模擬堰前水面曲線與模擬流量。通過模擬水面曲線可得有效水深，帶入堰流公式得到理論流量。模擬流量與理論流量結(jié)果見表3。由表3可知，隨著有效水深的增大，模擬流量與理論流量均呈指數(shù)增大趨勢，且在同一堰頂水頭下，模擬流量均略大于理論流量，但平均相對誤差僅為0.66%，最大相對誤差為1.09%，二者基本吻合。

表3 模擬流量與理論流量對比

由于控制模擬水深與實測水深完全一致具有一定難度，因此，二者的驗證采取對比其流量-水深的線性擬合方程完成，具體結(jié)果如圖9所示。模擬流量趨勢線的斜率為21.775，實測流量趨勢線的斜率為20.923，二者相近，表明變化趨勢基本一致，且2條趨勢線截距差為0.051，在位置上也十分靠近。隨流量的增大，二者數(shù)值誤差不斷減小，證明模擬流量與實測流量結(jié)果吻合良好。通過理論流量與實測流量的雙重驗證，證明本研究針對薄壁三角堰渠道采用的數(shù)值模擬方法合理，可用于進一步的計算分析。

圖9 模擬流量與實測流量關(guān)系

3.2 過堰水流的流速分布規(guī)律

取薄壁三角堰渠道內(nèi)的數(shù)值模擬結(jié)果，自入口開始，在上游渠道選取5個典型斷面（圖10），即在-平面=0.15 m、=0.30 m、=0.45 m、=0.60 m、=0.75 m處做以速度為變量的contour圖（圖11）。同時取-平面=0.90 m處的流速分布圖，得到堰板處的流速分布圖（圖12）。

圖10 橫斷面位置

由圖11可知，水流流速沿軸方向?qū)ΨQ分布，最大值出現(xiàn)在水氣交界面位置，且流速隨水流從邊壁向中心靠近的過程中逐漸增大。隨著水流在行進方向上的不斷推移，邊壁緩流區(qū)的范圍也不斷縮小，這是由于模擬流動開始前，上游渠道設(shè)定已存在0.3×0.2×0.9 m3的水，而后流速受邊壁影響逐漸調(diào)整為更接近實際情況的分布狀態(tài)，符合渠道水流行進規(guī)律。由圖12可知，水流在經(jīng)過堰板時仍舊表現(xiàn)為水體自由表面處流速最大。此外，因水流受到堰板的阻擋與導(dǎo)流，產(chǎn)生一種豎窄作用使得中心部分水流流速增大，呈現(xiàn)出以堰口為中心，流速環(huán)狀向外逐漸減小的變化過程。

圖11 X-Y平面（展向）典型橫斷面流速分布

圖12 堰板處流速分布

3.3 不同淤積工況下的流量變化

設(shè)定堰前有效水深約為0.030 0 m，取模擬流量、實測流量與理論流量隨淤積厚度的變化情況繪制于圖13。隨淤積厚度的增長，模擬流量與實測流量都出現(xiàn)了無序變化，而理論流量恒定為0.221 0 L/s，未隨淤積的變化產(chǎn)生適應(yīng)性改變，且均小于模擬流量，75%情況下小于實測流量。

圖13 各淤積工況下的流量變化

將模擬流量、實測流量與理論流量的相對誤差列于表4。淤積厚度為0.0 cm時，模擬流量、實測流量與理論流量的相對誤差分別為0.23%、0.45%，數(shù)值較為吻合；當淤積厚度為2.5～20.0 cm時，不同淤積工況下模擬流量與理論流量出現(xiàn)最大為5.88%的相對誤差，且平均相對誤差較表3中無淤積時兩者的平均相對誤差由0.66%增長至4.10%；不同淤積工況下實測流量與理論流量出現(xiàn)最大為8.60%的相對誤差，且其平均相對誤差為3.73%?？梢?，淤積會使薄壁三角堰理論流量值偏小，存在約-0.45%～8.60%的誤差。

3.4 過堰水流的水位變化特征

當理論流量為0.221 0 L/s時，取各淤積工況下的數(shù)值模擬結(jié)果，做-平面=0.15 m的切片處理，在此面上生成水的體積分數(shù)為0.5的等值線圖，得到各淤積工況下過堰水流的模擬水位變化情況。將各淤積工況下的實測堰前水位與模擬堰前水位繪制于坐標軸中（圖14）。各淤積工況下模擬堰前水面曲線走勢較為一致，均表現(xiàn)為初始水面平直，靠近堰板時出現(xiàn)不同程度的水位壅高，高度隨淤積厚度的增加波動上升。而后受堰板垂向收縮的作用，以一條光滑的降落曲線跌落。實測堰前水面曲線由于人工測量誤差，存在些許波動，但水位也在靠近堰板時出現(xiàn)不同程度的抬升，與模擬堰前水面曲線的特征基本一致。證明淤積會導(dǎo)致渠道水流在靠近薄壁三角堰時產(chǎn)生水位壅高，從而影響渠道過流能力，導(dǎo)致理論流量偏小。

表4 理論流量為0.221 0 L/s時不同淤積工況下實測流量、模擬流量與理論流量的相對誤差

圖14 理論流量為0.221 0 L/s時各淤積工況下的模擬與實測堰前水面曲線

4 討論

數(shù)值模擬方法的選取是影響試驗結(jié)果準確性的重要因素。根據(jù)前人研究經(jīng)驗，采用RNG k-ε模型與VOF模型模擬薄壁三角堰渠道水流，流量模擬結(jié)果較為可靠，實測流量略大于模擬流量的原因可能為實際測流渠道中存在一定的沿程損失，而模擬中對壁面做了簡化處理，沿程損失可以忽略不計。具體分析如下：總水頭由位置水頭、流速水頭和沿程損失構(gòu)成，模擬水流沿程損失忽略不計，堰板前兩者總水頭一致，若滿足水面高度相同的條件，則入水口處模擬平均流速大于實測平均流速，根據(jù)明渠流公式，可知實測流量略大于模擬流量。

此外，流速分布符合渠道水流行進規(guī)律，證明水力學特征模擬結(jié)果無誤。但由于流動開始前在上游渠道設(shè)定已存在0.3×0.2×0.9 m3的水，導(dǎo)致入口斷面流速分布呈現(xiàn)為上下差異較大且過渡段較短的分布情況，在流速逐漸調(diào)整為符合實際情況的分布后，從能量的角度分析，入口斷面處流速動能會略小于調(diào)整后流速分布動能，忽略短距離內(nèi)的沿程損失，即總水頭不變的情況下，入口處水深大于調(diào)整后水深，表現(xiàn)為隨著水流的行進，水面出現(xiàn)一定程度的降低。因此，在條件允許的情況下，最好采用UDF方法在入口處將流速設(shè)置為符合實際情況的不均勻分布，以避免這樣的情況。

通過分析同一設(shè)定水深，不同淤積工況下理論流量、模擬流量與實測流量的結(jié)果發(fā)現(xiàn)，隨著渠道內(nèi)淤積厚度的增加，理論流量數(shù)值恒定，而模擬流量和實測流量出現(xiàn)無規(guī)律的波動變化，其中模擬流量均表現(xiàn)為偏大，實測流量在75%的淤積情況下偏大?？紤]到人工測量會存在一定的誤差，因此有一定把握可以認為，淤積的存在會使薄壁三角堰過流能力增大。其原因可能是：當淤積存在且厚度小于等于堰頂高時，淤積對渠底起到了等同于抬升的作用，則在堰上水深不變時，上游渠道的過流區(qū)域縮小，水的流速加快，導(dǎo)致堰口處過堰流量增大。在理論公式中，雖然D值會根據(jù)堰頂高與渠道行進寬度的比值和有效水深與堰頂高的比值的變化而變化，但以本文選取的情況為例，其變化程度無法適應(yīng)淤積造成的影響，使得測流結(jié)果偏小。在之后的研究中，可以就淤積情況下流量公式的修正系數(shù)進行更深入的研究。

當理論流量相同時，比較各淤積工況下模擬與實測堰前水面曲線，發(fā)現(xiàn)薄壁三角堰上游渠道遠堰板處水面線平直，近堰板處出現(xiàn)一定程度的水位壅高。分析其主要原因為：淤積擠占了堰板前兩側(cè)與底部受阻擋和導(dǎo)流的水流的空間，加之流量的增大使得從堰口集中涌出的水流變得擁擠，為了緩解這種擁擠，因而發(fā)生了堰前水位的升高。此外，模擬流量與實測流量基本一致的情況下，實測水流的動能大于模擬水流的動能，因而實測水面線略低于模擬水面線。

5 結(jié)論

實測流量數(shù)據(jù)與模擬流速分布規(guī)律表明，采用計算流體力學軟件FLUENT19.0中的RNG k-ε模型與VOF模型模擬薄壁三角堰渠道水流是準確可靠的。分析試驗中所得的不同淤積工況下的模擬流量和實測流量，以及對應(yīng)的堰前水面曲線得知，淤積會增大薄壁三角堰過流能力，并使堰前產(chǎn)生水位壅高，而堰流公式中缺少適應(yīng)這種變化的修正參數(shù)，導(dǎo)致生產(chǎn)中所得理論流量偏小。

[1] 史紅玲, 戴清, 袁玉平, 等. 引黃灌區(qū)泥沙處理措施及提水設(shè)施的減淤作用[J]. 泥沙研究, 2003(3): 12-16.

SHI Hongling, DAI Qing, YUAN Yuping, et al. Effects of pumping on sediment treatment in irrigation districts along the Yellow river[J]. Sediment Research, 2003(3): 12-16.

[2] 徐霖玉, 路明, 李冉, 等. 無壩引水明渠淤積成因物理模擬試驗研究[J]. 水電能源科學, 2015, 33(1): 108-110.

XU Linyu, LU Ming, LI Ran, et al. Physical simulation test study of deposition causes of diversion open channel without dam[J]. Water Resources and Power, 2015, 33(1): 108-110.

[3] 高紅艷, 魏占民. 鐙口揚水灌區(qū)輸水干渠來沙特性解析[J]. 灌溉排水學報, 2018, 37(7): 99-105

GAO Hongyan, WEI Zhanmin. Analysis of the factors affecting sedimentation in canals ditching Yellow river water for irrigation[J]. Journal of Irrigation and Drainage, 2018, 37(7): 99-105.

[4] 趙志華, 吳文勇, 王佳盛, 等. 引黃灌渠泥沙遷移特性與渠道挾沙力模型試驗研究[J]. 灌溉排水學報, 2019, 38(10): 63-71.

ZHAO Zhihua, WU Wenyong, WANG Jiasheng, et al. Experimental study on sediment transport characteristics and channel sediment carrying capacity model of Yellow river irrigation canal[J]. Journal of Irrigation and Drainage, 2019, 38(10): 63-71.

[5] 史立紅, 王云輝, 王思勤, 等. 山東小開河引黃灌區(qū)高質(zhì)量發(fā)展的對策措施[J]. 中國農(nóng)村水利水電, 2022(4): 166-169, 175.

SHI Lihong, WANG Yunhui, WANG Siqin, et al. High quality development in Shandong Xiaokaihe irrigation area[J]. China Rural Water and Hydropower, 2022(4): 166-169, 175.

[6] 王婧, 韓啟彪, 那巍. 山西柳林灌區(qū)引黃水的泥沙特性分析[J]. 灌溉排水學報, 2018, 37(S2): 140-143.

WANG Jing, HAN qibiao, NA Wei. Characteristics of mud and sand in the water of the Yellow river in Liulin irrigation district of Shanxi province[J]. Journal of Irrigation and Drainage, 2018, 37(S2): 140-143.

[7] MAHA Rashad Fahmy. Effect of sediment deposition on the efficiency of Fayoum weir[J]. Flow Measurement and Instrumentation, 2015, 46: 133-138.

[8] 陳大宏, 陳娓. 溢流堰水流的三維模擬[J]. 武漢大學學報(工學版), 2005, 38(5): 54-58, 64.

CHEN Dahong, CHEN Wei. Three-dimensional simulation of flow over weirs[J]. Engineering Journal of Wuhan University, 2005, 38(5): 54-58, 64.

[9] 李玲, 陳永燦, 李永紅. 三維VOF模型及其在溢洪道水流計算中的應(yīng)用[J]. 水力發(fā)電學報, 2007, 26(2): 83-87.

LI Ling, CHEN Yongcan, LI Yonghong. Three-dimensional VOF model and its application to the water flow calculation in the spillway[J]. Journal of Hydroelectric Engineering, 2007, 26(2): 83-87.

[10] 李然, 李洪, 李嘉, 等. 氣液兩相流理論在明渠水氣界面計算中的應(yīng)用[J]. 水動力學研究與進展, 2002(1): 77-83.

LI Ran, LI Hong, LI Jia, et al. Application of gas-liquid two-phase theory for the water surface calculation in open channels[J]. Chinese Journal of Hydrodynamics, 2002(1): 77-83.

[11] 柳雙環(huán), 馬孝義, 陳超飛. U形渠道三角剖面堰水力特性的數(shù)值模擬[J]. 節(jié)水灌溉, 2018(4): 37-41.

LIU Shuanghuan, MA Xiaoyi, CHEN Chaofei. Numerical simulation of hydraulic characteristics of triangular profile weirs in U-shape channel[J]. Water Saving Irrigation, 2018(4): 37-41.

[12] 徐玲君, 陳剛, 李國棟, 等. 薄壁堰泄流能力的數(shù)值模型計算及模擬自由水面的評價[J]. 西北農(nóng)林科技大學學報(自然科學版), 2010, 38(1): 225-230.

XU Lingjun, CHEN Gang, LI Guodong, et al. Numerical simulation about rectangle sharp-crested weir and evaluation of free-surface[J]. Journal of Northwest Agriculture and Forest University (Natural Science Edition), 2010, 38(1): 225-230.

[13] 魏文禮, 李盼盼, 白朝偉, 等.梯形薄壁堰流量規(guī)律數(shù)值模擬研究[J]. 應(yīng)用力學學報, 2016, 33(4): 659-664, 740.

WEI Wenli, LI Panpan, BAI Zhaowei, et al. Numerical simulation of discharge law for a trapezoid thin-walled weir[J]. Chinese Journal of Applied Mechanics, 2016, 33(4): 659-664, 740.

[14] BIRTE Domnik, SHIVA Prasad Pudasaini, ROLF Katzenbach, et al. Coupling of full two-dimensional and depth-averaged models for granular flows[J]. Journal of Non-Newtonian Fluid Mechanics, 2013, 201(11): 56-68.

[15] SHARAREH Mahmodinia, MITRA Javan, AFSHIN Eghbalzadeh. The flow field and free surface pattern of the submerged side weir with different lengths[J]. Arabian Journal for Science and Engineering, 2014, 39(6): 4 461-4 472.

[16] 魏文禮, 張澤偉, 白朝偉, 等. 雙曲線型薄壁堰泄流規(guī)律的數(shù)值模擬[J].西北農(nóng)林科技大學學報(自然科學版), 2016, 44(8): 213-218.

WEI Wenli, ZHANG Zewei, BAI Zhaowei, et al. Numerical simulation of discharge of a hyperbolic weir[J]. Journal of Northwest Agriculture and Forest University (Natural Science Edition), 2016, 44(8): 213-218.

[17] 水工建筑物與堰槽測流規(guī)范. SL537—2011[S].

Code for Measurement of Discharge by Hydraulics Structures Weirs and Flumes. SL537—2011[S].

[18] 朱亞磊. 平坦V形堰過流特性數(shù)值模擬研究[D]. 楊凌: 西北農(nóng)林科技大學, 2010.

ZHU Yalei. Numerical simulation study on behavior of flow in flat V-weir[D]. Yangling: Northwest Agriculture and Forest University, 2010.

[19] 張敏. 弧底梯形短喉道量水槽水力性能數(shù)值模擬[D]. 楊凌: 西北農(nóng)林科技大學, 2018.

ZHANG Min. Numerical simulation on hydraulic performance of short-throat flume in arc-based trapezoidal[D]. Yangling: Northwest Agriculture and Forest University, 2018.

[20] 鐘新銘, 馬孝義, 楊珮珮. 弧底梯形渠道無喉道量水槽水位流量關(guān)系數(shù)值模擬[J]. 水土保持研究, 2016, 23(2): 340-345.

ZHONG Xinming, MA Xiaoyi, YANG Peipei. Numerical simulation on discharge relation of arc-based trapezoidal canaland cut-throat flow measurement flume[J]. Research of Soil and Water Conservation, 2016, 23(2): 340-345.

[21] SHEN Yongming, NG C O, ZHENG Yonghong. Simulation of wave propagation over a submerged bar using the VOF method with a two-equation k-ε turbulence modeling[J]. Ocean Engineering, 2004, 31(1): 87-95.

[22] 藍敏樂, 譚博仁, 許東兵, 等. 管型混合澄清槽內(nèi)的液-液兩相流的數(shù)值模擬[J]. 化工學報, 2021, 72(4): 1 965-1 974.

LAN Minle, TAN Boren, XU Dongbing, et al. Numerical simulation of liquid-liquid two-phase flow in tubular mixer-settler[J]. The Chemical Industry and Engineering Society of China Journal, 2021, 72(4): 1 965-1 974.

Performance of Triangular Crested Weir under Impact of Sediment Deposition

GUANG Huibing1, WANG Jian1*, LI Zesen1, CAO Bozhao1, LIU Chao1, WANG Kanghong2

(1. Northwest Agriculture and Forest University, Yangling 712100, China; 2. Hancheng District Soil and Water Conservation Workstation, Hancheng 714000, China)

【Objective】Weir is a device to measure water flow in open channels. Its accuracy and robustness depend on many factors. In this paper, we studied the performance of the triangular crested weir as impacted by sediment deposition in open channels.【Method】The study was based on computational fluid dynamic modelling. We simulated water flow under nine conditions ranging from no sediment deposition to full sediment deposition using the software FLUENT 19.0. The simulated results were verified against experimental measurements.【Result】①The simulation results were reliable and robust as comparison with experimental data showed that the simulated water flow rate and velocity distribution in the channels agreed well with those experimentally measured. ②When the theoretical flow rate was 0.221 0 L/s, the simulated and measured flow rates fluctuated erratically as the thickness of the sediment deposition increased, and the simulated flow water differed from the theoretical one with the difference varying from -0.45% to 8.60%. ③When the theoretical flow water was 0.221 0 L/s, comparison of the free-surface profiles with and without sediment deposition showed that the simulated and measured profiles occurred at different degrees of water level accumulation near the weir plate; the height of water level increased erratically with the increase in sediment deposition.【Conclusion】Comparison with experimental data showed that the simulation results are accurate and reliable. The water level in front of the weir rises because of the sediment deposition, and the flow capacity of the weir increases as a result. However, as the wire traditional formula does not have a correction term to account for sediment deposition, the flow rate estimated is smaller than the real ones.

triangular sharp-crested weir; deposition; computational fluid dynamics; flow rate; velocity

廣彗冰, 王健, 李澤森, 等. 淤積工況下薄壁三角堰流量特征研究[J]. 灌溉排水學報, 2022, 41(10): 139-146.

GUANG Huibing, WANG Jian, LI Zesen, et al. Performance of Triangular Crested Weir under Impact of Sediment Deposition[J]. Journal of Irrigation and Drainage, 2022, 41(10): 139-146.

S277.7

A

10.13522/j.cnki.ggps.2021651

1672 - 3317（2022）10 - 0139 - 08

2022-01-03

國家自然科學基金項目（41771308）

廣彗冰（1999-），女。碩士研究生，主要從事土壤侵蝕研究。E-mail: ghice@nwsuaf.edu.cn.

王?。?973-），男。教授，博士生導(dǎo)師，主要從事土壤侵蝕與流域管理研究。E-mail: wangjian@nwsuaf.edu.cn

責任編輯：韓 洋