李延巖， 孫宏軍

(1.天津大學(xué) 電氣與自動(dòng)化工程學(xué)院，天津 300072；2.過(guò)程檢測(cè)與控制重點(diǎn)實(shí)驗(yàn)室，天津 300072)

?

研究與探討

水平管氣液兩相鈍體繞流穩(wěn)定性的實(shí)驗(yàn)研究*

李延巖1,2， 孫宏軍1,2

(1.天津大學(xué) 電氣與自動(dòng)化工程學(xué)院，天津 300072；2.過(guò)程檢測(cè)與控制重點(diǎn)實(shí)驗(yàn)室，天津 300072)

通過(guò)對(duì)水平管氣液兩相鈍體繞流中的泡狀流和段塞流進(jìn)行實(shí)驗(yàn)研究，從重復(fù)性、斯特勞哈爾數(shù)(St)隨體積含氣率的變化和不同體積含氣率的斯特勞哈爾數(shù)相關(guān)性三個(gè)方面對(duì)渦街穩(wěn)定性進(jìn)行分析，得出體積含氣率在15 %以?xún)?nèi)能形成穩(wěn)定渦街，大于15 %時(shí)，穩(wěn)定渦街受到破壞，但仍有漩渦脫落的結(jié)論。實(shí)驗(yàn)結(jié)果表明：三種方法所得結(jié)果具有一致性，且穩(wěn)定可靠。

氣液兩相； 鈍體繞流； 斯特勞哈爾數(shù)； 重復(fù)性

0　引　言

鈍體繞流現(xiàn)象廣泛存在于工業(yè)現(xiàn)場(chǎng)和實(shí)際生活中。氣液兩相鈍體繞流較常見(jiàn)的情形：火力發(fā)電廠中鍋爐汽包供給蒸汽的計(jì)量、核電站輕水反應(yīng)堆[1]、各種殼管式熱交換器、原油開(kāi)采與輸運(yùn)等。

國(guó)內(nèi)外學(xué)者對(duì)氣液兩相鈍體繞流做了許多有益的研究工作。Hulin J P等人[2]研究了氣液兩相繞流兩鈍體的渦街脫落特性，指出來(lái)流截面含氣率小于10 %時(shí)可得到穩(wěn)定的卡門(mén)渦街。Yokosawa M等人[3]用平板研究?jī)上嗯轄盍鞯奈擦鳎空宫F(xiàn)出兩相流尾流的形狀并討論了其特性。盧家才等人[4，5]研究了矩形截面管內(nèi)氣液兩相沖刷柱體時(shí)，斯特勞哈爾數(shù)(St)數(shù)和脈動(dòng)升力系數(shù)隨雷諾數(shù)和來(lái)流含氣率的變化關(guān)系。戴瑋等人[6]分別對(duì)電導(dǎo)和電容傳感器在氣水和油氣的條件下進(jìn)行兩相流層流標(biāo)定實(shí)驗(yàn)，得出電導(dǎo)和電容傳感器的輸出均和管道內(nèi)含液率有良好的線性關(guān)系。李永光等人[7～9]研究了空氣和水兩相繞流兩種尺寸的T形柱體產(chǎn)生渦街時(shí)氣液兩相斯特勞哈爾數(shù)的變化規(guī)律，獲得了兩相斯特勞哈爾數(shù)、來(lái)流截面含氣率、水流量三者之間的關(guān)系。何泳成等人[10]基于雙截面電阻層析成像技術(shù)進(jìn)行氣液兩相流的相關(guān)速度測(cè)試研究，分析了相關(guān)速度的測(cè)量原理，提出了基于截面含氣率的相關(guān)測(cè)速方法。

本文研究不同體積含氣率β和流速下氣液兩相鈍體繞流的特性。

1　渦街分離原理

流體的繞流運(yùn)動(dòng)會(huì)產(chǎn)生旋渦分離，目前常用發(fā)生體為三角柱結(jié)構(gòu)，如圖1所示，發(fā)生體左側(cè)箭頭為流體運(yùn)動(dòng)方向，u為來(lái)流速度，p為壓力，AB為迎流面,AE和DB為流體截面、直角坐標(biāo)系x，y。

圖1　渦街分離原理示意圖

流體流經(jīng)AB迎流面， x方向流速降為0，y方向流速先增大后減小。B處的實(shí)際流速方向?yàn)閗。實(shí)際流體有粘性，在斜邊AC和BC上產(chǎn)生邊界層，邊界層外流體作減速增壓，邊界層內(nèi)流體出現(xiàn)倒流現(xiàn)象。到達(dá)一定程度產(chǎn)生邊界層分離現(xiàn)象，在后端形成兩列漩渦，即卡門(mén)渦街。渦街頻率f與來(lái)流流速u(mài)的關(guān)系為

(1)

式中f為渦街頻率，Hz；St為斯特勞哈爾數(shù)；u為管道來(lái)流速度，m/s；d為發(fā)生體迎流面寬度，m；D為管道直徑，m。

2　實(shí)驗(yàn)系統(tǒng)

在天津大學(xué)過(guò)程檢測(cè)與控制重點(diǎn)實(shí)驗(yàn)室三相流實(shí)驗(yàn)裝置DN50管段上進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)系統(tǒng)如圖2所示。

圖2　實(shí)驗(yàn)系統(tǒng)裝置圖

兩相介質(zhì)為水和空氣?？諝饨?jīng)過(guò)濾器過(guò)濾、壓縮機(jī)加壓和冷干機(jī)降濕后送入儲(chǔ)氣緩沖罐2，儲(chǔ)氣罐和計(jì)量管段之間安裝穩(wěn)壓閥，保證氣路壓力穩(wěn)定。調(diào)節(jié)控制閥可調(diào)節(jié)氣流量的大小。從引射器噴出的氣液混合流經(jīng)過(guò)直管段充分發(fā)展后，進(jìn)入實(shí)驗(yàn)管段。實(shí)驗(yàn)條件：水為連續(xù)相，氣為離散相，溫度為20 ℃，液相流量點(diǎn)為：7,8,9,10 m3/h；體積含氣率：0 %～20 %。實(shí)驗(yàn)流型：β>3 %時(shí)，為段塞流；β<4 %時(shí)，為泡狀流。采用有機(jī)玻璃透明管道，在發(fā)生體后1.5D(D為水平管道內(nèi)徑)處安裝壓電探頭檢測(cè)渦街信號(hào)，在發(fā)生體前5D處取壓，以檢測(cè)流場(chǎng)壓力變化。渦街信號(hào)和壓力信號(hào)采用NI—6009板卡進(jìn)行采集并通過(guò)LabVIEW軟件實(shí)現(xiàn)實(shí)時(shí)數(shù)據(jù)記錄。

3　實(shí)驗(yàn)結(jié)果與分析

3.1含氣率對(duì)渦街信號(hào)的影響

對(duì)于氣液兩相流來(lái)說(shuō)，離散相的加入必會(huì)對(duì)渦街的產(chǎn)生和發(fā)展造成影響，且隨著離散相的增多，這種影響會(huì)不斷加強(qiáng)。下面以體積含氣率β為特征量，分析空氣的加入給渦街信號(hào)帶來(lái)的影響。

由圖3知，β=0 %時(shí)渦街信號(hào)穩(wěn)定性最好且信號(hào)幅值最大。當(dāng)β≤15 %時(shí)，隨著β增加，渦街信號(hào)變得不規(guī)則，穩(wěn)定性逐漸變差，有很大衰減，但仍能采集到較穩(wěn)定的渦街信號(hào)。當(dāng)β>15 %時(shí)，隨著β增加，渦街信號(hào)產(chǎn)生很大波動(dòng)，穩(wěn)定性極差，仍有旋渦脫落，但已檢測(cè)不出穩(wěn)定渦街信號(hào)。當(dāng)β>20 %時(shí)，采集的信號(hào)基本上是干擾信號(hào)?？傻?，β=15 %為渦街信號(hào)穩(wěn)定與不穩(wěn)定的分界點(diǎn)。

圖3　不同體積含氣率下的渦街信號(hào)圖

3.2渦街穩(wěn)定性分析

3.2.1渦街信號(hào)重復(fù)性分析

氣液兩相流中，隨著體積含氣率β不斷增大，渦街有一定耗散，且逐漸不穩(wěn)定，實(shí)驗(yàn)的可重復(fù)性就變差。重復(fù)性σ表征測(cè)量結(jié)果的一致性，選 σ作為判斷氣液兩相流穩(wěn)定性的依據(jù)，其公式如下

(2)

σ=[σi]max

(3)

式中fij為第i個(gè)實(shí)驗(yàn)點(diǎn)第j次檢測(cè)的渦街頻率值；fi為第i個(gè)實(shí)驗(yàn)點(diǎn)的渦街平均頻率值；σi為第i個(gè)實(shí)驗(yàn)點(diǎn)渦街頻率值的重復(fù)性。

由圖4可知，當(dāng)體積含氣率β≤15 %時(shí)，每個(gè)流量下各實(shí)驗(yàn)點(diǎn)的渦街頻率重復(fù)性在5 %以?xún)?nèi)，重復(fù)性較好，流場(chǎng)比較穩(wěn)定；當(dāng)體積含氣率β>15 %時(shí)，渦街頻率重復(fù)性基本上大于5 %，重復(fù)性變差，流場(chǎng)變得不穩(wěn)定。

圖4　重復(fù)性隨體積含氣率變化圖

3.2.2氣液兩相流斯特勞哈爾數(shù)分析

斯特勞哈爾數(shù)是分析氣液兩相流動(dòng)穩(wěn)定性的主要參數(shù)。在單相流中，斯特勞哈爾數(shù)(St)在一定雷諾數(shù)(Re)范圍內(nèi)為常數(shù)，在氣液兩相流中，St由式(4)計(jì)算可得

(4)

式中

(5)

式中f為氣液兩相渦街頻率，Hz； w為渦街發(fā)生體體征尺寸，m；u為氣液兩相混合物經(jīng)過(guò)發(fā)生體時(shí)的平均流速，m/s；A為發(fā)生體兩側(cè)有效截面積，如圖5中陰影部分所示，m2；Qvg為液相體積流量，m3/h；Qvl為氣相體積流量，m3/h；

圖5　發(fā)生體前部截面圖

由圖6知，β≤15 %時(shí)，各流量點(diǎn)下St隨β 增大而減小，相關(guān)性較強(qiáng)； β>15 %時(shí)，St變化雜亂無(wú)章，沒(méi)有明顯規(guī)律，可推斷每個(gè)流量下β≤15 %時(shí)有穩(wěn)定渦街產(chǎn)生。

圖6　斯特勞哈爾數(shù)隨體積含氣率變化散點(diǎn)圖

3.2.3斯特勞哈爾數(shù)St相關(guān)系數(shù)分析

由圖7知,對(duì)各流量點(diǎn)分β≤15 %,15 %<β≤20 %兩種情形的St隨β變化進(jìn)行線性擬合，對(duì)比相關(guān)系數(shù)和標(biāo)準(zhǔn)差。以8m3/h為例，當(dāng)β≤15 %時(shí)，近似成線性關(guān)系，擬合關(guān)系式為y=-0.073 6×10-4x + 0.165 2，相關(guān)系數(shù)R為-0.889 2，標(biāo)準(zhǔn)差為0.001 53；當(dāng)15 %<β≤20 %時(shí)，R為-0.152 8，標(biāo)準(zhǔn)差0.004 14，明顯不相關(guān)，且誤差很大。

圖7　St隨β變化圖

7,9,10 m3/h流量點(diǎn)的St隨β變化擬合結(jié)果如表1所示：β≤15 %時(shí)，St隨β增大而線性減小，線性相關(guān)性好，說(shuō)明渦街較穩(wěn)定；15 %<β≤20 %時(shí)，St變化與β變化相關(guān)性不強(qiáng)，偏差較大，說(shuō)明此β下無(wú)穩(wěn)定渦街生成。所以渦街從穩(wěn)定到不穩(wěn)定是一個(gè)漸變的過(guò)程。

表1　各流量點(diǎn)在兩種情形下St的相關(guān)系數(shù)

4　結(jié)　論

本文從重復(fù)性、氣液兩相斯特勞哈爾數(shù)(St)隨體積含氣率β的變化和斯特勞哈爾數(shù)(St)的相關(guān)性三個(gè)方面對(duì)渦街穩(wěn)定性進(jìn)行分析，得出相同流量下，當(dāng)β≤15 %時(shí)，有穩(wěn)定渦街產(chǎn)生;當(dāng)β>15 %時(shí)，依然存在旋渦脫落，但由于氣相對(duì)旋渦耗散的作用，渦街穩(wěn)定性逐漸被破壞直至無(wú)法形成。

[1]Schroeder H J.Chemistry in nuclear-power stations with light water-reactors[J].Chemiker-zeitung, 1984, 108(2):55-64.

[2]Hulin J P，F(xiàn)ierfort C，Condol R.Experimental study of vortox emission behind bluff obstacles in a gas-liquid vertical two-phase flow[J].International Journal of Multiphase Folw，1982，8(5)：475-490.

[3]Yokosawa M, Kozawa Y, Inoue A.Studies on two-phase cross flow, part II:Transition Reynolds number and drag and drag coefficient[J].Int’ J of Multiphase Flow, 1986, 12(2):169-184.

[4]Lu Jiacai，Xie Zhengwu，Lin Zonghu.Study on the fluctuating lift acting on a circular cylinder under two-phase cross-flow [J].Chinese Journal of Applied Mechanics，2000，17(1)：76-80.

[5]Lu Jiacai，Xie Zhengwu，Wang Yangpeng.Pressure distribution around circular cylinder for gas-liquid two-phase cross flow[J].Journal of Xi’an Jiao Tong University，1999，33(1)：55-58.

[6]戴瑋，譚超，董鋒.水平管道油氣水三相流含水率測(cè)量[J].傳感器與微系統(tǒng)，2015，34(1):135-138.

[7]李永光，林宗虎.垂直上升氣液兩相流中三角形柱體兩相斯托拉赫數(shù)的研究[J].水動(dòng)力學(xué)研究和進(jìn)展，1999，14(4):438-443.

[8]李永光，林宗虎.氣液兩相渦街穩(wěn)定性的研究[J].力學(xué)學(xué)報(bào)，1998，30(2):138-144.

[9]李永光，林宗虎.氣液兩相繞流柱體時(shí)的渦街特性的研究[J].應(yīng)用力學(xué)學(xué)報(bào)，1994,11(1):37-40.

[10] 何泳成，王化祥，張立峰，等.基于截面含氣率的氣/液兩相相關(guān)速度測(cè)量[J].傳感器與微系統(tǒng)，2009，28(10):103-106.

Experimental research on stability of gas-liquid two-phase flow around bluff body in horizontal pipe*

LI Yan-yan1,2， SUN Hong-jun1,2

(1.School of Electrical Engineering and Automation, Tianjin University,Tianjin 300072,China;2.Key Laboratory of Process Measurement and Control, Tianjin 300072，China)

Experimental researches are carried out on bubble flow and slug flow of flow around bluff body in horizontal gas-liquid two-phase flow.Stability of vortex street is analyzed based on repeatability, variation trend ofStand correlation research ofStwith changing volume void fraction.It is concluded that the flow field can produce a stable vortex street when volume void fraction is within 15 % but the stable vortex street is destroyed and it can still produce vortex street when volume void fraction is greater than 15 %.Experimental studies show that the results of three methods are consistent and have good stability.

gas-liquid two-phase flow； flow around bluff body； Strouhal number； repeatability

10.13873/J.1000—9787(2016)09—0005—04

2015—11—11

國(guó)家自然科學(xué)基金青年科學(xué)基金資助項(xiàng)目(50906061)

TH 814

A

1000—9787(2016)09—0005—04

李延巖(1990-)，男，河南焦作人，碩士研究生,主要研究方向?yàn)閭鞲衅骷夹g(shù)與信號(hào)處理、流體參數(shù)檢測(cè)與數(shù)值仿真。