彎管角度對充水閥過流能力影響的數值計算研究

舒崚峰，劉 婧，方 杰，李勝兵，姚志峰，3

(1.中國電建集團華東勘測設計研究院有限公司，浙江杭州311122；2.中國農業大學水利與土木工程學院，北京100083；3.中國農業大學北京市供水管網系統安全與節能工程技術研究中心，北京100083)

0 引 言

水力機組進水口閘門一般均為靜水開啟，開啟前需要利用庫水向閘門后至機組前的隧洞或者壓力管道充水平壓，充水閥是目前國內最常用的充水平壓方式。受水輪發電機組漏水、上游水位、充水閥流量系數等因素的綜合作用，充水平壓過程中，當閘門前后水位差達到某一值后充水量和漏水量會實現平衡，所以實際閘門開啟時存在1～5 m的水位差，但隨著機組設計水頭、額定容量的不斷加大，已有若干工程出現了8 m的水位差還無法實現充水量和漏水量平衡，使進水口閘門無法開啟并影響發電的情況。

為保證進水閘門在規定水位差內實現正常開啟，需對影響充水閥過流能力的因素進行深入分析。由于以往設計充水閥對于大部分電站均能滿足要求，為此國內對于閘門充水閥的研究并不多：夏維紅[1]通過公式分析和水力試驗分析了充水閥的空穴特性，劉素英[2]從工程實際出發介紹了充水閥的設計要點，肖段龍[3]等經理論計算與分析推導了充水閥流量和開度的相互關系，龔洪波[4]等對快速閘門充水方式進行優化改進，鄧亞新[5]對已運行多年的充水閥改造要點進行了總結，陶云冬[6]結合抽水蓄能電站充水特點揭示該類充水閥設計及運用中需要注意的要點。

本文依托某水利工程，利用數值計算方法(CFD)對充水閥在不同出水彎管角度和不同開度下的過流能力進行了分析。

1 研究對象與數值計算方案

1.1 充水閥基本參數

本文選取的充水閥出口彎管角度為工程上常用為90°和105°。閥門剖面圖如圖1所示，閥門開度由閥體向上提升的行程表示：圖中閥體行程為0 mm，處于全關位置；閥體向上行程250 mm為全開位置。

圖1 充水閥結構

1.2 模型建立與網格劃分

1.2.1模型簡化

考慮上游側來流通過充水閥流至下游側的距離短，其對閘門上方流場影響范圍很有限。為此，在控制總體計算量和計算精度的前提下，將閥門井流體域限定在閘門上方3 m處。此外，充水閥上游流場流速低，水力損失不顯著，故對充水閥上游側模型做一定簡化，最終簡化計算模型如圖2所示。

圖2 計算模型

1.2.2網格劃分

計算采用非結構化網格，對主體流場計算區域(閥體附近位置)做網格加密，對于非主要流體計算區域采取一定的疏化處理。網格單元總數約為380 000，網格節點總數約為65 000，網格劃分情況如圖3所示。

圖3 網格劃分

1.3 計算設置

通過充水閥的水流主要來自閘門上游，設置上游流體進口為20 kPa的總壓進口，并將閥門井流體域限定在閘門上方3m處。充水閥過流量相對于閘門下游側的引水隧洞體積非常小，故不考慮引水隧洞內流場。出口定義在閥門彎管出口，采用0 kPa靜壓出口，并設置流量監測點。最終邊界條件如圖4所示。

圖4 計算邊界

本次數值計算采用k-ε湍流模型，該湍流模型具有穩定、簡單、經濟的特點[7]，對于一般工程均能獲得較為滿意的計算精度。為充分了解彎管角度對不同充水閥過流能力的影響，分別選取25、75、100、150、250 mm共5個開度進行計算。

2 計算結果與分析

2.1 收斂性分析

收斂情況用最后100步計算結果的峰-峰值與平均值的比值反映。整體收斂情況顯示：大開度下，計算結果較為穩定；在小開度下，計算結果具有較高振蕩性。收斂振蕩性與實際流場紊亂程度呈現關聯性。

圖6 速度矢量場

圖7 總壓分布

2.2 流場分析

圖5 全局流場情況(100 mm開度)

全局流線圖顯示(見圖5)，從上游側進入閘門上方水體的水流，呈螺旋式流向充水閥進口，但水體流速非常低。水力損失主要集中在閥體附近及彎管區域。速度矢量場顯示(見圖6)，小開度下，閥門出口的水流以類似射流的方式進入彎管，在內、外側壁處伴有漩渦形成；大開度下，閥門內流場順暢，閥門出口流速均勻；在彎管段，離心作用導致二次流的形成，彎管外側處流體降速增壓，內側處流體升速降壓。隨著開度增大，旋渦范圍逐漸減小，這增大了彎管處的過流面積，減小了能量損耗。此外，對比90°與105°彎管出口流場，大開度下，105°彎管出口流場分布更加均勻、流暢；小開度下，兩種結構紊亂程度均非常高。

圖7顯示了總壓的變化情況。總壓變化情況與速度場的變化相符合。大開度下：內側壁出口附近存在部分能量耗散；彎管段主流區能量在速度能與壓力能之間轉換，幾乎無損失。隨著開度減小，內側壁出口能量損耗區增大，且主流區受到漩渦區域的影響，速度能與壓力能的轉換過程中能量也被消耗，主流區也無法再恢復到原先的總壓能。

圖8顯示了壓力梯度的變化情況。壓力梯度的大小反應了能量轉換過程中速度能與壓力能變化最為劇烈的位置。從圖中可以發現，閥體附近一直為壓力能變化最劇烈的地方；對于彎管段，開度較大時，彎管內側壁能量轉換劇烈，而在小開度下，彎管外側壁有較大的能量轉換。

2.3 過流特性分析

2.3.1阻力系數

閥門的過流能力可由其阻力系數反應，流量與阻力系數的關系如下：

圖8 壓力梯度

(1)

式中，Q為通過閥門的流量，m3/h；H表示閥門前后水頭壓差，m；K為阻力系數(流量系數倒數)，m0.5/(m3/h)。

計算區域內的水力損失主要集中在閥體附近，故按式(1)計算阻力系數時，可將計算域進出口壓差作為閥門前后壓差。根據計算得到90°和105°兩種彎管出口下的充水閥流量-開度曲線(見圖9)和阻力系數-開度曲線(見圖10)，可以發現：0～100 mm開度范圍內，增加閥門開度能顯著提高閥門過流能力，降低阻力系數，過流量增加幅度約為5.82 kg·s-1/mm；150～250 mm開度范圍內，阻力系數趨向恒定，增大開度對提升過流能力的幫助有限，過流量增加幅度約為2.0 kg·s-1/mm；此外，0～50 mm開度范圍內，改變彎管角度對充水閥過流能力幾乎無影響；當開度超過75 mm時，105°出口彎管比90°出口彎管更能有效提高充水閥的過流能力，最大提高約16.4%。

圖9 流量—開度曲線計算結果

圖10 阻力系數—開度曲線計算結果

2.3.2水力損失

流速平均化對比情況(見圖11)顯示：大開度下，質量平均和面積平均的計算結果差別較小；小開度下，兩者計算結果差異較大。由于小開度下的有效過流面積小，閥門出口過流斷面流態已不符合過流斷面上各點流速均勻分布的假設，故采用斷面平均流速無法真實反映損失特性。為此，均按照質量平均考慮水力損失。最終計算得到不同開度下的閥門水損情況如表1所示，計算結果可為充實閥的設計和選擇提供參考。

圖11 不同平均化流速對比

3 結 論

對某水力工程進水口閘門的DN500充水閥流場進行了數值模擬，研究了開度及出口彎管角度對充水閥過流能力和流場特性的影響情況，得到主要結論如下：

表1 充水閥水損

注：壓力用相應的水柱高度表示，進口質量平均總壓均為2.04 m。

(1)流場隨閥門開度變化規律為：0～100 mm開度下，彎管在內、外側壁形成漩渦，減小過流面積，增大能量損耗，開度增大能快速降低漩渦影響范圍，增大過流能力；100 ～250 mm開度下，流體順暢通過閥門，彎管出口流速較均勻，增大開度對流場的改善作用有限。此外，當出口未出現回流現象時，105°彎管出口流場分布更加均勻、流暢；而在小開度下，兩種結構紊亂程度均非常高。

(2)閥門開度增大能夠提高閥門過流能力、降低阻力系數。開度小于100 mm時，過流能力提升效果顯著，每毫米開度的平均流量增長率可達5.82 kg/s；開度大于150 mm時候，過流能力提升效果降低，每毫米開度的平均流量增長率為2.0 kg/s。

(3)開度為0～75 mm時，改變閥門出口彎管角度對閥門過流能力影響不大；開度為75～250 mm時，105°彎管相比90°彎管能有效提高閥門的過流能力，整個開度范圍內的平均過流量增加14.1%，最大過流量增加16.39%。