賀登輝, 陳森林, 白博峰

(1.西安理工大學(xué)省部共建西北旱區(qū)生態(tài)水利國(guó)家重點(diǎn)實(shí)驗(yàn)室,陜西西安 710048;2.西安交通大學(xué)動(dòng)力工程多相流國(guó)家重點(diǎn)實(shí)驗(yàn)室,陜西西安 710049)

氣液分層流是濕氣多相流(濕天然氣、濕蒸汽等)中的一種重要的流動(dòng)形態(tài),廣泛存在于工農(nóng)業(yè)生產(chǎn)過(guò)程中。氣液分層流的測(cè)量非常重要。目前雖然有許多種方法可以測(cè)量氣液分層流流量,但應(yīng)用最多的仍然是差壓法[1-2]。作為一種新型的差壓式流量計(jì),V錐流量計(jì)因其具有信號(hào)穩(wěn)定、壓損低、量程比寬、所需直管段短等優(yōu)點(diǎn)[3-7],近年來(lái)在多相流測(cè)量領(lǐng)域受到了越來(lái)越多關(guān)注。對(duì)V錐流量計(jì)測(cè)量濕氣時(shí)的差壓響應(yīng)特性研究表明, V錐流量計(jì)濕氣測(cè)量特性主要受氣、液相含量、壓力及節(jié)流比等因素影響[3, 7-11]。Stewart等[8]采用Steven文丘里關(guān)聯(lián)式[10]的建立方法,得到對(duì)不同節(jié)流比V錐流量計(jì)的測(cè)量關(guān)聯(lián)式;Steven等[3,12]對(duì)管徑為101.6和152.4 mm,節(jié)流比為0.75的V錐流量計(jì),改進(jìn)了de Leeuw關(guān)聯(lián)式測(cè)量模型;徐英等[13]通過(guò)對(duì)V錐流量計(jì)測(cè)量濕氣時(shí)產(chǎn)生的“過(guò)讀”進(jìn)行修正,建立了一套濕氣測(cè)量模型;愛(ài)默生公司[14]基于V錐流量計(jì),開(kāi)發(fā)了Roxar濕氣流量計(jì),并成功應(yīng)用于實(shí)際油氣井氣液流量測(cè)量中,但該公司未公布其V錐流量計(jì)具體采用何種修正模型。由于V錐流量計(jì)在濕氣測(cè)量時(shí)比文丘里管對(duì)液相測(cè)量精度的依賴更低[9,15],這對(duì)于提高其氣相測(cè)量精度十分有利。目前對(duì)于V錐流量計(jì)關(guān)注的重點(diǎn)多是不同結(jié)構(gòu)參數(shù)及流動(dòng)參數(shù)下V錐流量計(jì)的差壓響應(yīng)特性,但V錐節(jié)流裝置內(nèi)的氣液兩相流動(dòng)特性及其對(duì)差壓測(cè)量的影響機(jī)制知之較少。筆者針對(duì)實(shí)際濕氣中常見(jiàn)的氣液分層流采用數(shù)值模擬方法對(duì)V錐節(jié)流裝置內(nèi)氣液分層流動(dòng)特性進(jìn)行研究。

1 數(shù)值模擬方法及模型

1.1 數(shù)值模擬方案

V錐節(jié)流元件的結(jié)構(gòu)如圖1所示。節(jié)流元件由前、后錐角分別為φ和θ的兩個(gè)V形錐體組成,并且由支撐桿固定在管道上;高壓取壓口位于V錐元件上游,低壓取壓口位于后錐體的頂點(diǎn)處,穿過(guò)錐體由支撐桿引出管外。管道內(nèi)徑為50 mm,V錐的節(jié)流比為0.55,前、后錐角分別為45°和135°。

圖1 V錐節(jié)流裝置Fig.1 V-Cone throttle device and structure

基于ANSYS Fluent平臺(tái),對(duì)空氣、水流經(jīng)V錐節(jié)流裝置的流動(dòng)特性進(jìn)行數(shù)值模擬。由于來(lái)流是分層流,主要包括光滑分層流、波狀分層流及滾動(dòng)波狀流,不考慮分層流內(nèi)液滴的運(yùn)動(dòng),采用VOF模型捕捉氣液界面[16]。如圖2(a)所示,所模擬的V錐裝置上、下游直管段分別為5D和9D,這與試驗(yàn)中的測(cè)壓點(diǎn)布置保持一致。采用速度入口邊界條件,分為兩部分,其中氣相由中間的圓形區(qū)域進(jìn)入管道,液相由環(huán)形區(qū)域進(jìn)入。對(duì)于本文中研究工況,氣液兩相經(jīng)過(guò)5D的長(zhǎng)度,在進(jìn)入V錐節(jié)流段時(shí)可達(dá)到氣液分層流態(tài);出口采用壓力出口邊界條件,壁面無(wú)滑移;模擬的工況進(jìn)出口邊界條件均為試驗(yàn)值。湍流模型采用RNGk-ε模型[17];采用結(jié)構(gòu)化六面體網(wǎng)格對(duì)模擬區(qū)域進(jìn)行劃分,并在錐體附近進(jìn)行網(wǎng)格加密處理(圖2(b))。

圖2 幾何模型及網(wǎng)格劃分Fig.2 Geometry model and mesh

1.2 數(shù)學(xué)模型

由于涉及的數(shù)學(xué)模型眾多,僅對(duì)幾個(gè)主要的模型進(jìn)行介紹,其他模型可以參考ANSYS Fluent幫助文獻(xiàn)[18]。

1.2.1 連續(xù)性方程

采用VOF模型模擬氣、液兩相流時(shí),兩相流體使用同一個(gè)方程組,每一相的體積分?jǐn)?shù)在整個(gè)計(jì)算域內(nèi)被追蹤。通過(guò)求解其中一相體積分?jǐn)?shù)的連續(xù)方程來(lái)跟蹤氣、液兩相之間的界面。第q相的連續(xù)性方程為

(1)

由于不考慮相變,傳質(zhì)及源項(xiàng)均取為零。

氣、液兩相的體積分?jǐn)?shù)之和滿足

αg+αl=1.

(2)

連續(xù)性方程采用隱式時(shí)間離散格式進(jìn)行求解,單元內(nèi)的物性參數(shù)根據(jù)氣、液兩相體積分?jǐn)?shù)加權(quán)平均得到。氣、液混合物的平均密度ρ和黏度μ分別為

ρ=αgρg+(1-αg)ρl,

(3)

μ=αgμg+(1-αg)μl.

(4)

式中，μg和μl分別為氣相和液相的動(dòng)力黏度，Pa·s。

1.2.2 動(dòng)量方程

動(dòng)量方程通過(guò)氣、液混合密度ρ和動(dòng)力黏度μ與體積分?jǐn)?shù)相聯(lián)系,在整個(gè)計(jì)算區(qū)域內(nèi)求解同一動(dòng)量方程,得到的速度場(chǎng)也被兩相共同使用。動(dòng)量方程為

ρg+Fvol.

(5)

其中

式中,p為壓力,Pa;Fvol為作用于氣、液間相界面的體積力,N/m3;κg為曲率,1/m。

1.2.3 能量方程

由于V錐節(jié)流的影響,流體溫度在錐體附近可能會(huì)發(fā)生變化,因此,考慮了溫度對(duì)能量傳遞的影響。氣液兩相共用能量方程為

(6)

式中,keff為有效熱導(dǎo)率,由氣液兩相共用,W/(m·K);T為溫度,K;E為能量,J。

VOF模型中取兩相質(zhì)量加權(quán)平均值計(jì)算能量及溫度,對(duì)于能量有

(7)

式中,Eq為按單相比熱容和共用的溫度計(jì)算的每一項(xiàng)能量,J。

1.3 模型驗(yàn)證

1.3.1 網(wǎng)格獨(dú)立性驗(yàn)證

采用3組不同的網(wǎng)格數(shù)進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證(圖3,其中Usg為表觀氣速,m/s;Usl為表觀液速,m/s),網(wǎng)格單元數(shù)依次增大約1.5倍,如表1所示。由圖3可知,網(wǎng)格2與網(wǎng)格3計(jì)算結(jié)果相比,沿流動(dòng)方向(X軸方向)的靜壓力分布變化很小,可認(rèn)為網(wǎng)格2能夠滿足計(jì)算需求。

圖3 網(wǎng)格數(shù)量變化時(shí)壓力模擬結(jié)果對(duì)比Fig.3 Wall pressure profiles under different meshes

網(wǎng)格節(jié)點(diǎn)數(shù)單元數(shù) 網(wǎng)格1963199936360 網(wǎng)格214879811452880網(wǎng)格322140012168992

1.3.2 試驗(yàn)數(shù)據(jù)驗(yàn)證

采用網(wǎng)格2方案進(jìn)行模擬,得到沿流動(dòng)方向靜壓力分布,與試驗(yàn)結(jié)果(試驗(yàn)裝置和方法見(jiàn)參考文獻(xiàn)[4])進(jìn)行對(duì)比 (圖4)。由圖4可知,兩者吻合良好,最大相對(duì)誤差小于5.0%,表明所建立的數(shù)值方法可靠,能夠反映出氣液兩相流在V錐節(jié)流裝置內(nèi)的壓力分布特性。

圖4 數(shù)值模擬與試驗(yàn)得到的沿流動(dòng)方向靜壓對(duì)比Fig.4 Comparison of wall pressure profile between numerical simulation and experimental results

2 數(shù)值結(jié)果分析

2.1 錐尾尾渦特性

2.1.1 尾渦基本特征

當(dāng)來(lái)流為單相氣體時(shí),如圖5(p=0.2 MPa,Usg=5.71 m·s-1,Usl=0)所示,單相氣體流經(jīng)V錐,在喉部形成環(huán)形射流,射流撞擊壁面,發(fā)生反彈,并在喉部下游一定距離處射流速度達(dá)到最大;在射流剪切作用下,V錐錐尾區(qū)域壓力梯度變化較大,錐尾下游出現(xiàn)了尾流,形成了尾流漩渦(尾渦),尾渦渦心速度較低。流動(dòng)速度在尾渦下游端點(diǎn)附近達(dá)到最小,沿流動(dòng)方向速度接近于0處的位置即為尾渦的下游端點(diǎn)。尾渦尺度較大,長(zhǎng)度Lvortex約為錐體直徑d的2～3倍;由于V錐節(jié)流裝置邊壁收縮,尾渦靠近管壁處的流動(dòng)方向與主流相同,中心區(qū)的流動(dòng)方向與主流相反(圖5(b)和(c))。

當(dāng)來(lái)流為氣液分層流時(shí),V錐下游尾渦特征呈現(xiàn)出上部大下部小的形態(tài),靠近管道下壁面處的尾渦長(zhǎng)度與單相氣體中的相比,其長(zhǎng)度大大縮短。例如,圖6(p=0.2 MPa,Usg=5.71 m·s-1,Usl=0.057 m·s-1)工況中的尾渦長(zhǎng)度與圖5的工況相比,縮短了約40%。其主要影響機(jī)制為:①流體物性的影響。單相氣體中加入液相后導(dǎo)致流體的混合密度和黏度等物性均發(fā)生變化,而密度增大幅度遠(yuǎn)遠(yuǎn)大于黏度,再加上喉部射流速度相對(duì)于單相氣體急劇增大,從而導(dǎo)致流體雷諾數(shù)增大,加劇了氣液動(dòng)量交換,能量損失增大,尾渦長(zhǎng)度縮短;②V錐喉部射流的影響。射流使氣液撞擊壁面,液體發(fā)生反彈和破碎,少量液體被卷吸進(jìn)尾渦中,消耗了尾渦的能量。同時(shí)隨液量增大,反彈的液體高度增大[4],從而阻礙了尾渦發(fā)展。

圖5 單相氣體條件下的流動(dòng)Fig.5 Streamline and velocity vector distribution under single-phase gas flow

圖6 氣液兩相流條件下的流動(dòng)Fig.6 Streamline and velocity vector distribution under annular flow

尾渦形成和發(fā)展對(duì)沿流動(dòng)方向的靜壓分布產(chǎn)生影響。圖7為單相氣體和氣液兩相流條件下壁面靜壓沿流動(dòng)方向分布(p=0.2 MPa)。由圖7可知,流體流經(jīng)V錐,壓力逐漸降低,在錐尾下游某處壓力達(dá)到最低,之后逐漸恢復(fù);其恢復(fù)長(zhǎng)度受尾渦的影響,尾渦越長(zhǎng),則壓力恢復(fù)所需的距離越長(zhǎng);隨著尾渦影響的減弱,下游壓力也逐漸恢復(fù),經(jīng)過(guò)一定長(zhǎng)度的過(guò)渡區(qū)域,當(dāng)流線與主流方向趨于平行時(shí),流動(dòng)即可恢復(fù)。圖7中的氣液兩相流尾渦長(zhǎng)度與單相氣體中的基本相等,因此所需的壓力恢復(fù)長(zhǎng)度也相差不大。

圖7 流動(dòng)對(duì)管壁靜壓的影響Fig.7 Effect of vortex on wall pressure

2.1.2 影響因素

單相氣體中引入液相,使得錐后尾渦發(fā)生變化。圖8為表觀液速對(duì)尾渦的影響(刻度尺單位為mm,上圖為XOY截面,下圖為XOZ截面)。由于尾渦是立體的,為了更全面的認(rèn)識(shí)尾渦,給出了其在XOY和XOZ兩個(gè)截面上的分布。單相氣體中,尾渦沿軸線基本對(duì)稱(圖8(a));加入液相后,尾渦開(kāi)始呈現(xiàn)上部大下部小的形態(tài),且隨著表觀液速增加,尾渦長(zhǎng)度發(fā)生變化,且尾渦下游端點(diǎn)越來(lái)越向管道下部?jī)A斜。這意味著管道上部流場(chǎng)的穩(wěn)定區(qū)域增大,尾渦對(duì)下游的影響也減弱。尾渦形狀的變化主要是由于分層流來(lái)流時(shí),液相主要分布在管道下部,因此下部流速較低,V錐喉部射流作用較弱,導(dǎo)致了尾渦向管道下部?jī)A斜,而XOZ截面上尾渦分布的對(duì)稱性基本不受表觀液速影響。此外,氣、液表觀流速基本相同的情況下,尾渦基本不受來(lái)流入口壓力的影響。在本文氣液分層流范圍內(nèi),流體溫度變化較小(小于1 ℃),因此溫度對(duì)測(cè)量的影響亦較小。

圖8 液相對(duì)尾流漩渦的影響Fig.8 Effect of liquid velocity on vortex

2.1.3 尾渦對(duì)下游流動(dòng)影響

尾渦影響下游流體速度及壓力分布特性。圖9為來(lái)流經(jīng)過(guò)V錐之后,上下游截面速度沿Y軸方向分布,其中0D表示V錐錐尾取壓口所在的管道橫截面?？梢园l(fā)現(xiàn),無(wú)論在單相氣體還是氣液分層流條件下,錐尾取壓口附近(r/R=0,0D截面)流速基本均為0,說(shuō)明錐尾取壓口附近的流動(dòng)較為穩(wěn)定,有利于低壓測(cè)量;當(dāng)來(lái)流為單相氣體時(shí)(圖9(a)),上游-1D處流動(dòng)較為規(guī)則;從下游1D截面上的流速分布可知,從管道中心軸線r/R=0沿Y軸方向到r/R=1范圍內(nèi),流速先減小后增大,這是由于該范圍位于錐尾尾渦區(qū)域之內(nèi)(圖8),尾渦的影響使得流速分布發(fā)生變化。當(dāng)來(lái)流為分層流時(shí),位于尾渦區(qū)內(nèi)的截面上也存在類(lèi)似的速度分布(圖9(b)～(d))。對(duì)于分層流,管道上下速度并不對(duì)稱,上部速度大于下部速度,且表觀液速越大,其不對(duì)稱性越明顯;速度分布亦與圖8的尾渦分布特征相吻合。

在同一截面上,隨著表觀液速增大,沿Y軸方向的上部流速越來(lái)越大,而下部流速越來(lái)越小,如圖10所示。由圖7可知,流體流過(guò)V錐之后,壓力恢復(fù)位置與尾渦的影響密切相關(guān)。當(dāng)表觀液速增大時(shí),V錐差壓也隨之增大,同時(shí)下游恢復(fù)壓力也越小,壓力損失越大。對(duì)于分層流來(lái)流,下游壓力恢復(fù)位置受表觀液速影響較小(圖11),這與尾渦對(duì)下游流動(dòng)的影響規(guī)律相一致。

圖9 不同橫截面上沿Y軸方向的速度分布Fig.9 Velocity along Y axis direction under different cross sections

圖10 V錐下游3D和6D截面上速度沿Y軸的分布隨表觀液速的變化Fig. 10 Velocity along Y axis direction with different superficial liquid velocity at 3D and 6D downstream V-cone

圖11 不同表觀液速下壁面靜壓力沿流動(dòng)方向X的分布Fig.11 Wall static pressure along flow direction X under different superficial liquid velocity

2.2 氣液相分布變化及其對(duì)測(cè)量影響

對(duì)V錐上游1D、V錐前端、錐尾取壓口所在的管道截面及下游3D和6D處管道橫截面上的液相分布進(jìn)行分析。當(dāng)來(lái)流為分層流時(shí)(圖12(a)Usg=5.64 m·s-1,Usl=0.007 m·s-1;(b)Usg=5.63 m·s-1,Usl=0.028 m·s-1;(c)Usg=5.71 m·s-1,Usl=0.057 m·s-1;(d)Usg=5.79 m·s-1,Usl=0.11 m·s-1),氣相中液相極少,流經(jīng)V錐后,尾渦卷吸少量液體,同時(shí)在喉部高速射流作用下,管道下部液膜變薄,在V錐下游一定距離處才逐漸恢復(fù);喉部射流的影響還導(dǎo)致液體發(fā)生飛濺,并有少量液體進(jìn)入氣相中。試驗(yàn)中還發(fā)現(xiàn),氣液兩相流流經(jīng)V錐,液膜在后錐體表面吸附,而在尾渦的影響下,液膜并沒(méi)有完全覆蓋后錐體表面,而是存在一個(gè)液膜邊界,使得取壓口附近基本沒(méi)有液體分布(圖13)。液相在V錐表面的這種分布是由“科恩達(dá)效應(yīng)”(Coanda effect)引起的。尾部取壓口附近基本沒(méi)有液相分布,從而使得低壓取壓孔處于一個(gè)相對(duì)穩(wěn)定的流動(dòng)環(huán)境中,十分有利于氣液兩相差壓測(cè)量。

圖12 液相體積含率分布Fig.12 Water volume fraction

圖13 試驗(yàn)觀察到的錐后液體分布Fig.13 Liquid distribution observed in experiment

3 結(jié) 論

(1)V錐錐尾下游形成大尺度的尾流漩渦。分層流來(lái)流條件下,由于液相主要分布在管道下部,因此下部流速較低,同時(shí)射流強(qiáng)度較弱,導(dǎo)致尾渦呈現(xiàn)上部大下部小的分布形態(tài),表觀液速越大,尾渦下游端點(diǎn)越向管道下部?jī)A斜。氣液兩相流的平均密度與黏度均比單相氣體中的高,以及V錐喉部環(huán)形射流撞擊管壁產(chǎn)生的反彈液滴擴(kuò)散現(xiàn)象,是導(dǎo)致V錐后尾渦特征發(fā)生變化的主要原因。

(2)V錐下游壓力恢復(fù)長(zhǎng)度和橫截面上的速度分布與尾渦長(zhǎng)度和形狀變化密切相關(guān)。對(duì)于分層流,管道上下部速度并不對(duì)稱,上部速度大于下部速度,其不對(duì)稱性隨表觀液速增大而增加。

(3)錐尾取壓口附近的流動(dòng)較為穩(wěn)定,且錐尾取壓孔附近液相含量極低,使得低壓取壓環(huán)境較為穩(wěn)定,進(jìn)而有利于差壓測(cè)量。