充水速度對空氣閥駝峰管段水力特性的影響

2019-11-28 06:45:12張宏禎李燕輝王玉成

中國農村水利水電 2019年11期

張宏禎，李燕輝，蔣勁，王玉成，羅爽

(1. 甘肅省疏勒河流域水資源局，甘肅玉門 735211；2. 武漢大學水力機械過渡過程教育部重點實驗室，武漢 430072)

0 引言

在長距離輸水工程中，需要沿管線安裝一定數量的空氣閥，使充水過程可排出管內的大量空氣及管道排空過程補進空氣。空氣閥一般安裝在管線局部高點，以避免產生液柱分離和過度的正負壓力波動[1,2]。目前，對空氣閥管道系統水力特性已有較廣泛的研究。在數值模擬方面，王玲等[3]通過改進的牛頓迭代法和直接求解法結合的空氣閥數學模型，分析了單向進氣閥、雙向進排氣和進氣微排等三種空氣閥型對空管充水過程的瞬態特性及水錘防護效果的影響。Lee和Leow[4]提出空氣閥及滯留氣團的改進數值模型，同時對比了空氣閥不同進氣系數和排氣系數對管道停泵壓力瞬變過程的影響。Jung和Karney[5]結合基因算法和粒子群算法，通過瞬態分析對空氣閥尺寸及安裝位置進行優化，并提出了優化的水錘防護方案。Stephenson[6]創建了諾模圖來預估空氣閥快速關閉產生的二次水錘壓力，同時指出了管道流量和水錘壓力的影響參數，包括空氣閥安裝豎管的直徑和長度、空氣閥尺寸和安裝位置、關閥后豎管中的殘余氣體等。在實驗方面，Bergant等[7]通過大型管道實驗臺研究了空氣閥的水力特性，包括單相水過流、動態排氣、動態進氣和同時進排氣等四種工況。Carlos等[8]通過測試空氣閥排氣工況下管道不同位置的流量和壓力，對管道充水過程中的液柱分離和氣團運動數值模型進行了修正。

在實際操作中，空管充水過程應將管內滯留空氣充分排出，否則嚴重的水錘壓力將對管道造成破壞[9]。在安裝空氣閥的基礎上，也要控制好充水速度以免管內空氣從空氣閥過快排出引起二次水錘效應。美國水務協會(AWWA)建議管道充水速度低于0.3 m/s(滿管流速)，該方案比較低效，特別對長距離管道來說非常耗時。以往研究并未涉及大流量充水方案的對比，即探討在較大充水速度下瞬態壓力波動是否在管道能承受的安全范圍內，并分析管道及空氣閥的動態水力特性。另一方面，以往對管道充水排氣過程的數值模擬主要應用一維數值模型[10, 11]，例如特征線法和剛性水柱模型。這些模型實際上是基于相關假設：流動是一維的，并且不考慮流固耦合效應；由剛性水柱代表液相，忽略其可壓縮性，而用氣體多方方程來描述滯留氣團的變化過程，且氣液交界面為垂直于管軸線的剛性面[12, 13]。然而，基于以上假設的一維數值模型在某些情況下不能得到合理的結果，例如在充水速度較小時水柱在重力作用下將首先沿管內氣團的下方流動而將氣體往局部高點擠壓，此時兩相交界面與管道垂直及其剛性的假定就不再有效[14]，另外剛性水柱模型也會高估初始時刻的充水流速[15]。為克服上述種種問題，本文應用格子玻爾茲曼方法(LBM)對裝有空氣閥的駝峰管道充水排氣過程進行CFD數值模擬，建立簡化的駝峰管段分析模型來代表實際管道的局部高點(空氣閥安裝位置)，分析較大充水速度(0.3～0.8 m/s)對系統水力特性的影響。研究結果可為實際輸水管道充水操作方案提供參考價值。

1 分析建模

1.1 物理模型及工況參數

建立空氣閥及簡化的駝峰管段幾何模型如圖1所示，其中，DN50的單體式空氣閥由不銹鋼浮球(重179 g)、閥座及閥體組成。管道內徑為100 mm，充水前空氣閥為開啟狀態，管內充滿空氣。下游閥門在充水過程中保持關閉，在建模中設為壁面邊界。另外，在數值模擬中將空氣閥浮球設置為自由剛體，其運動由重力及浮力的合力控制。開始充水時，流進管內的水體對空氣擠壓，通過空氣閥排出管外。當水體充滿管道時，浮球在浮力的作用下上升并堵住空氣閥出口，完成充水排氣過程。

圖1 空氣閥及駝峰管段幾何模型Fig.1 Geometric models of air valve and hump pipe

常規操作規程建議管道充水速度低于0.3 m/s，而在本次數值模擬中分別考慮了0.3、0.4、0.5、0.6、0.7、0.8 m/s等6種充水方案，對比分析較大的充水速度對空氣閥管段動態水力特性的影響，特別是殘余氣團體積分數和二次水錘壓力。而管道進口和空氣閥出口對應的邊界條件則分別設為速度進口和壓力出口(大氣壓)。

1.2 數值方法

對管道充水過程的氣液兩相流動進行CFD數值模擬可采用傳統的數值方法，如歐拉模型、混合模型和VOF模型[16, 17]等。而本次研究采用了無網格的格子玻爾茲曼方法(LBM)，相關文獻[18-20]已從理論和實驗等方面驗證了其對兩相流模擬的適用性。與傳統的CFD方法求解Navier-Stokes方程不同，LBM方法通過離散的玻爾茲曼方程來描述介觀尺度下的流體運動，可提供清晰的物理概念和簡單的算法。此外，由于粒子間的相互作用，LBM模擬的兩相流動中相分離可自然產生，不需要對兩相界面進行動態重構和跟蹤。

在多相LBM模型中，在空間位置x和時間t下的總粒子密度分布函數如下[21]：

(1)

(2)

式中：k代表液相或氣相；x+eiVt代表x附近最近點的坐標。

(3)

在碰撞算子的組成中，第一項表示由非平衡態發展到局部平衡狀態的松弛過程，其線性化形式為：

(4)

(5)

式中：Fk(x)是由兩相相互作用而在各相產生的等效力。

另外，在碰撞步中強制實現了各相質量守恒和總動量守恒[22]：

(6)

式中：u為局部速度；ρu為局部總動量；液相密度ρw和氣相密度ρa之和組成局部總密度。

2 結果與討論

2.1 瞬態流動特性

各充水方案下水流初次到達駝峰頂部時刻的整體流速分布對比如圖2所示，水柱推動管內空氣通過空氣閥往外界排放，且充水速度越大，排氣速度越大。然而，只有液面附近的空氣才會受到充水水柱的擠壓和推動往空氣閥加速排放，而無論充水速度如何變化，下游管段內的氣體均不受上游水柱的影響而保持靜止或低速流動狀態。這是因為空氣閥的過流面積小流速大，形成局部低壓，產生強烈的分流作用，所以對下游管內空氣的擾動較小。另外，對于駝峰管段充水的前半程，當充水速度達到某個臨界值后，管內流速分布類似，如0.7與0.8 m/s充水方案對比所示。

圖2 水流初次到達駝峰管段頂部時的流速分布Fig.2 Velocity distributions at the instant of water column initially reaching the top end of the pipe based on different filling flow discharges

當水流充滿管道時，空氣閥浮球在浮力的驅動下上升并關閉閥門出口，而關閥時刻及關閥速度與充水速度大小有關。不同充水方案對應的浮球豎向位移變化如圖3所示，顯然，隨著充水速度增加，空氣閥關閥時刻提前。需要注意的是，在某些充水速度下，如0.4與0.5 m/s，會出現浮球竄動回彈現象，這是因為關閥階段空氣閥內水體入流不連續所致。特別是0.5 m/s充水方案，當水柱初次經過駝峰頂部時產生較大的液面波動，部分水流進入空氣閥使浮球竄動至頂部，會出現短暫的堵氣現象。就排氣平穩性而言，0.5 m/s的充水速度對本次研究的駝峰管段來說是不利的充水方案。

圖3 充水過程空氣閥浮球豎向位移變化Fig.3 Time histories of vertical displacement of the valve float in filling process based on different filling flow discharges

對應不同的充水方案，當空氣閥完全關閉后，管內存在不同程度的殘余氣團，氣團分布及體積率數值分別如圖4和圖5所示。在推薦的充水速度0.3 m/s工況中，不產生殘余氣團，這是因為充水水柱在重力作用下能完全沿管內氣體底部緩慢推進，通過空氣閥將氣體全部排出。當充水速度增大，殘余氣團量增加，且在關閥瞬間主要集中在下游管段，而大流量充水時殘余氣團則有向駝峰頂部發展的趨勢。在圖5中，隨著充水速度增大，關閥殘余氣團的體積分數將近似拋物線地增加。

圖4 關閥時刻殘余氣團分布Fig.4 Residual air distribution at the instant of air valve closure based on different filling flow discharges

圖5 關閥時刻殘余氣團體積分數Fig.5 Residual air volume fraction at the instant of valve closure

空氣閥排氣速度變化過程對比如圖6所示。各充水方案的前半程均為穩定排氣階段(圖中箭頭所指)，排氣速度在一恒定值附近小幅波動，且隨充水速度增加而整體上升。這是因為水柱液面越過駝峰頂部前，液面后方至空氣閥間的氣體流動阻力小，氣體受水柱推動而非強烈壓縮排出，排氣速度與液面推進速度接近。而后半程因水流在下游管段與管內氣體混流擾動，排氣速度波動明顯。整體來說，較大充水速度時關閥殘余氣團量增加，則對應的排氣總量減少，所以后半程排氣速度波動相比小充水速度方案要小。與圖3中0.5 m/s充水工況的浮球劇烈竄動情況對應，該方案的排氣速度波動也最大。對比各充水工況，除0.3 m/s充水速度(常規方案)外，排氣速度峰值均出現在空氣閥關閥時刻附近。整體上，隨充水速度增大，排氣速度峰值呈下降趨勢。而在0.5 m/s充水工況中因空氣閥發生短暫堵氣，所以浮球跌落后出現瞬間極大的排氣速度，其峰值在所有工況中最大。

圖6 空氣閥排氣速度變化過程Fig.6 Time histories of air release velocity through the air valve outlet based on different filling flow discharges

2.2 二次水錘壓力

二次水錘壓力是空氣閥快速關閉產生的壓力陡升，與系統參數有關[23]，包括空氣閥安裝位置、空氣閥尺寸、管道長度、管道縱剖面分布和波速等，而本次研究主要探究充水速度對二次水錘的影響。空氣閥及下游閥門位置的二次水錘壓力變化過程分別如圖7(a)和7(b)所示，兩者有相似的變化規律，即隨著充水速度增加，空氣閥關閉引起的二次水錘壓力呈線性增加趨勢。充水速度過快，管內氣體來不及完全排出，水柱對殘余氣團的強烈壓縮使壓力急劇增加。另外，下游閥門處的水錘正壓整體比空氣閥位置略大，這是因為前者有更大的水柱體量，沖擊能量也更大。