垂直下水管中氣液兩相流動特性分析

姚昱婷， 孫 陽,2， 林欣然， 程 玖

(1. 河海大學(xué) 港口海岸與近海工程學(xué)院， 江蘇 南京 210098； 2.南京水利科學(xué)研究院水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室， 江蘇 南京 210098)

1 研究背景

針對目前我國低水頭水力能利用不足的問題，提出了一種對環(huán)境影響小且更有效的解決方案：山區(qū)低落差河流環(huán)境下液氣泵系統(tǒng)[1]。該系統(tǒng)將低落差河流的水力能轉(zhuǎn)化為高壓氣能，再使用產(chǎn)生的高壓氣泵水，可以解決山區(qū)居民用水困難的現(xiàn)狀。本文的研究對象是該系統(tǒng)第一部分(利用低水頭的河流水力能將自然空氣轉(zhuǎn)化為高壓氣體)垂直下水管中的氣液兩相流的流動特性。

氣液兩相流現(xiàn)象普遍存在于自然界及人類的各種活動中[2-6]。國內(nèi)外對氣液兩相流的研究，主要基于垂直上升管道、傾斜管道和水平管道，對垂直下水管道的研究基本沒有開展。這些研究中，除Almabrok等[7]對直徑為101.6 mm的垂直下水管道上的氣液兩相流流型進(jìn)行了一項(xiàng)試驗(yàn)研究，許多研究人員研究了兩相上升流和下降流的流型，這些研究大多集中在直徑小于100mm的小型管道上[8-12]，沒有針對直徑100 mm以上大直徑圓管道的研究。受液氣黏滯性、表面張力的影響，不同直徑的圓管道可能產(chǎn)生完全不同的運(yùn)動特性，因而大直徑圓管道中氣液兩相流的研究具有重要研究和工程實(shí)際借鑒價值。本文基于均相流解析模型來分析垂直下水管道的氣液兩相流動特性，有助于更進(jìn)一步的得到該系統(tǒng)各影響因素與運(yùn)轉(zhuǎn)性能之間的規(guī)律，尤其是能夠進(jìn)一步探討包括液相體積流量在內(nèi)的關(guān)鍵因素對系統(tǒng)的復(fù)雜影響。

2 試驗(yàn)裝置

試驗(yàn)?zāi)Ｐ脱b置如圖1所示，主要由3部分組成：垂直下水管、水氣分離罐和出水管道。在實(shí)際工程中，這3個部分均需埋設(shè)在地下。本試驗(yàn)是在0℃～ 30℃的常溫和常壓下進(jìn)行的。在設(shè)備通電開啟后，水將從進(jìn)水管中流入垂直下水管，由于管道內(nèi)存在負(fù)壓，空氣會自發(fā)地混合到水中，在垂直下水管道中形成兩相流。氣相在垂直下水管道中被不斷壓縮，并在下落過程中與水充分混合，最終在下水管道末端形成高壓氣體。兩相流在水氣分離罐中發(fā)生氣相與液相的分離，氣相通過出氣管道排出，而液相通過出水口排出。試驗(yàn)中采用渦街流量計對液相體積流量Ql進(jìn)行測量，渦街流量計安裝在出水管末端。圖1中，dg為進(jìn)氣口管道直徑，d為垂直下水管直徑，d′為出水管直徑，Δh為水流落差，h′為出水管道高度。所采用的試驗(yàn)方案如表1所示。

圖1 試驗(yàn)裝置

3 均相流解析模型

為了簡化分析，基于假設(shè)：

(1) 管中的氣液兩相是充分混合且均勻分布的，忽略管道中液體顆粒和氣團(tuán)的橫向運(yùn)動，只考慮縱向運(yùn)動；

(2) 氣液兩相的實(shí)際速度是相同的，即氣液兩相流無相對滑移；

(3) 由于垂直下水管道底部的壓力以及流體與管道內(nèi)壁之間的摩擦勻被忽略，因此，主流是做自由落體運(yùn)動的。

根據(jù)垂直落水管中高壓氣形成的機(jī)理，可以建立一維流動的均相流模型[1]：

(1)

式中：β為容積含氣率；Vl為液相折算速度，m/s；L為垂直下水管道的長度,m。

表1 試驗(yàn)方案

4 試驗(yàn)結(jié)果與分析

4.1 β和h′的關(guān)系

取公式(1)關(guān)于L的一階偏導(dǎo)數(shù)，得到公式(2)：

(2)

根據(jù)公式(2)，β隨著L的增長而增長，且增長趨勢逐漸趨緩。L和Δh之間的數(shù)學(xué)關(guān)系是：

L=Δh+h′

(3)

落差Δh若為恒定值，根據(jù)公式(3)，則反映到β和h′的關(guān)系為：β隨h′的增大而增大。

圖2為4種試驗(yàn)工況不同流量下β與h′的實(shí)際關(guān)系。如圖2所示，顏色越淺，代表β越高，在h′方向顏色變化不明顯，可見h′對β的影響較小，與上述基于模型的預(yù)測存在一定的差別。按照模型的預(yù)測，β隨出水管高度h′的增大而增大。有研究發(fā)現(xiàn)，水流在排水立管道中流動，達(dá)到一定長度后，流速維持一個定值，不再遞增，這個長度稱為終限長度，流速稱為終限流速[13]。這項(xiàng)發(fā)現(xiàn)與本文的假想是一致的。h′的變化主要使垂直下水管道的后段部分長度發(fā)生變化，在垂直下水管道的后段部分，形成穩(wěn)定均勻的兩相流，兩相流的速度基本保持不變，因此幾乎不影響β。

圖2 4個試驗(yàn)工況不同流量下β與h′的實(shí)際關(guān)系圖

4.2 β和Ql的關(guān)系

圖3為垂直下水管道不同模擬直徑的液相體積流量數(shù)據(jù)范圍。圖3中的偏離點(diǎn)代表各次試驗(yàn)中的極大或極小工況，由此可知本裝置所采用的Ql在75～300 m3/h這個區(qū)間。因此，基于均相流模型，在Ql的區(qū)間內(nèi)分別對直徑為160、200、250 mm的垂直下水管道的幾組試驗(yàn)數(shù)據(jù)進(jìn)行簡單的模擬計算，即將Δh、d、Ql這3項(xiàng)參數(shù)代入均相流模型得到圖4。而實(shí)際的試驗(yàn)結(jié)果如圖5和6。

分析圖5和6，可以發(fā)現(xiàn)在不同試驗(yàn)工況下，β與Ql的曲線的變化趨勢均為由上升段或者水平段逐漸轉(zhuǎn)入下降段，并且存在一個轉(zhuǎn)折點(diǎn)，這些轉(zhuǎn)折點(diǎn)近似位于同一直線上。

圖3 液相體積流量數(shù)據(jù)范圍圖

圖4 β與Ql均相流模型預(yù)測圖

圖5 β與Ql實(shí)際關(guān)系圖(出水管道高度h′=5 m)

圖6 β與Ql實(shí)際關(guān)系圖(出水管道高度h′= 6 m)

根據(jù)圖5、6，本文假設(shè)在下降段區(qū)間內(nèi)，容積含氣率β與液相體積流量Ql近似呈線性關(guān)系。在不同Δh以及不同d的工況下，這段曲線變化幅度不大，液相體積流量Ql是影響容積含氣率β的主要因素?；谏鲜鲈囼?yàn)結(jié)果，下降段直線擬合結(jié)果為：

β=-7.969×10-4Ql+0.2779

(4)

如公式(4)所示，下降段中液相體積流量Ql與容積含氣率β為近似線性關(guān)系。轉(zhuǎn)折點(diǎn)在下降曲線上所處的位置受到Δh、h′ 和d這3個因素的影響，不同Δh、不同h′和不同d會影響進(jìn)入下降段轉(zhuǎn)折點(diǎn)的位置，但是這條擬合直線的位置幾乎不發(fā)生變動。經(jīng)過進(jìn)一步的觀測發(fā)現(xiàn)，較細(xì)的垂直下水管(Φ160 mm)的β與Ql關(guān)系曲線首先達(dá)到轉(zhuǎn)折點(diǎn)，較粗的垂直下水管(Φ250 mm)的β與Ql關(guān)系曲線后達(dá)到轉(zhuǎn)折點(diǎn)；較長的垂直下水管的β與Ql關(guān)系曲線首先達(dá)到轉(zhuǎn)折點(diǎn)，較短的垂直下水管的β與Ql關(guān)系曲線后達(dá)到轉(zhuǎn)折點(diǎn)。根據(jù)這條擬合出的下降曲線可以估算出一定流量的水通過本系統(tǒng)時，理想情況下最多可以產(chǎn)生的高壓氣量Qg：

Qg=βQl=(-7.969×10-4Ql+0.2779)Ql

(5)

即可認(rèn)為，只要本裝置各部分結(jié)構(gòu)安排合理，就能產(chǎn)生流量為Qg的高壓氣。Qg=Qg(Ql)是一條拋物型曲線，拋物線必然存在的下降段，顯然與實(shí)際Qg=Qg(Ql)呈現(xiàn)不斷上升的變化趨勢是不吻合的，但可以反映一定的趨勢。根據(jù)這樣一條拋物線，當(dāng)液相體積流量Ql達(dá)到174 m3/h時，Qg基本不再發(fā)生變化。

同時，在實(shí)際工程中可以根據(jù)河流的流量，選擇合理直徑的垂直下降管道。在一般情況下，盡量采用較小直徑的垂直下水管，其能量轉(zhuǎn)換效率較高；在天然河流的流量特別大的工況下，采用較大直徑，能量轉(zhuǎn)換效率不會受到影響，并且由于較細(xì)的垂直下水管道的過流能力有所限制，大流量的天然河流需要采用較大直徑的垂直下水管。根據(jù)圖5和6的數(shù)據(jù)，在實(shí)際工程中，對應(yīng)于特定流量的河流，垂直下水管道直徑d的選取應(yīng)當(dāng)遵循表2。需要注意的是，表2適用于Δh=2～7 m的工況。對于流量小于150 m3/h的情況，現(xiàn)有的試驗(yàn)數(shù)據(jù)尚不能反映，應(yīng)當(dāng)選擇更小直徑的垂直下水管。

表2 垂直下水管直徑d建議取值表

5 結(jié) 論

(1) 基于多工況試驗(yàn)數(shù)據(jù)，對垂直下水管中的氣液兩相流運(yùn)動特性進(jìn)行了分析，能夠反映垂直下水管中的氣液混合程度的關(guān)鍵參數(shù)β受出水管高度h′的影響較??；β隨著液相體積流量Ql的變化趨勢均為由上升段或者水平段逐漸轉(zhuǎn)入下降段，并且存在一個轉(zhuǎn)折點(diǎn)，這些轉(zhuǎn)折點(diǎn)近似位于同一直線上。據(jù)所得試驗(yàn)結(jié)果，為垂直下水管道直徑的選取提供了建議。

(2) 通過試驗(yàn)結(jié)果與基于一定假設(shè)提出的均相流模型的解析解對比，進(jìn)一步分析并論證了該模型的可靠性及存在的不足之處。

(3) 還需開展更多工況的試驗(yàn)，進(jìn)一步分析均相流模型存在誤差的原因以及為小流量情況下垂直下水管直徑的選擇提供試驗(yàn)依據(jù)。