鄧清華,何偉,張亮,彭傲然,趙卓斌,豐鎮(zhèn)平

(1.西安交通大學(xué)葉輪機械研究所,710049,西安;2.西安交通大學(xué)陜西省葉輪機械及動力裝備工程實驗室,710049,西安;3.上海核工程研究設(shè)計院有限公司,200233,上海)

均勻穩(wěn)定以及不存在明顯的流態(tài)畸變是能源、化工等工藝流程中旋轉(zhuǎn)機械、流量測量裝置對來流狀態(tài)的基本要求。非均勻的來流會導(dǎo)致泵產(chǎn)生振動、氣蝕、揚程降低等問題。當(dāng)流量計上游流動速度剖面不對稱時,至少要布置50倍管徑距離的直管段來削弱流態(tài)畸變的影響[1]。然而在實際管路布置條件下,彎管、三通、孔板、擴張管等結(jié)構(gòu)是非常常見的,流體流過這些結(jié)構(gòu)則必然出現(xiàn)不對稱速度分布和復(fù)雜二次流[2-3]。為了將上述流態(tài)畸變削弱到滿足泵或流量計要求的范圍,同時避免占用大量物理空間和增加系統(tǒng)成本,一般采用流動調(diào)整器裝置而不是鋪設(shè)幾十倍乃至上百倍管徑長的直管段。按照結(jié)構(gòu)的不同,流動調(diào)整器分為葉片式、管束式、格柵式、孔板式、翅片式等[4-6]。

葉片式流動調(diào)整器的典型代表是Etoile式流動調(diào)整器[6],其采用八個周向分布的葉片將流動空間分為均勻的扇形部分[7],優(yōu)點是易加工、結(jié)構(gòu)簡單、壓力損失低,但整流效果一般。Laws等研究發(fā)現(xiàn),較短的長度和空心設(shè)計能夠進(jìn)一步提高葉片式流動調(diào)整器的整流性能[8]。管束式流動調(diào)整器由布置在主管道內(nèi)的19根以上的軸向平行管束構(gòu)成[9],加工工藝簡單,但流動阻力較大。格柵式流動調(diào)整器的典型代表是AMCA流動調(diào)整器[10],其原理與管束式流動調(diào)整器類似,但其分割密度大于管束式流動調(diào)整器,整流性能更佳,缺點是其流動損失也更大,而且裝卸、檢修較為煩瑣[10]?？装迨搅鲃诱{(diào)整器是在垂直于管道的板件上布置對稱的多組不同開孔直徑的整流裝置[11],典型的有Laws調(diào)整器[12]、NEL調(diào)整器[13]、Zanker流動調(diào)整器[14]等?？装迨搅鲃诱{(diào)整器具有加工工藝簡單、整流性能強、壓力損失大的特點。

翅片式流動調(diào)整器是1991年由美國Vortab公司的Smith等提出的一種內(nèi)置翅片產(chǎn)生旋渦擾亂原始速度分布并利用二次流摻混形成均勻流場的整流裝置[15-16]。根據(jù)其專利介紹,該結(jié)構(gòu)能有效應(yīng)用于氣體、液體的整流中,可以以極低的壓力損失為代價在較短的距離內(nèi)實現(xiàn)均勻、無旋渦的速度分布[15]。目前,翅片式流動調(diào)整器按照結(jié)構(gòu)主要分為直管式和彎管式兩種,直管式調(diào)整器的翅片沿周向均勻布置,而彎管式則僅在彎管外側(cè)沿流向布置,其基本原理和結(jié)構(gòu)如圖1所示。由圖1可見,當(dāng)流體流經(jīng)管內(nèi)翅片時,在翅片迎風(fēng)面與背風(fēng)面的壓差作用下會產(chǎn)生強烈的剪切力[17],與航空領(lǐng)域的翼尖渦類似,這種剪切力會形成發(fā)夾渦并進(jìn)一步演變?yōu)閷D(zhuǎn)渦,大大促進(jìn)邊界層內(nèi)外的質(zhì)量與能量交換[18]。發(fā)卡渦結(jié)構(gòu)是湍流邊界層中的基礎(chǔ)結(jié)構(gòu),由Theodorsen于1952年首次發(fā)現(xiàn)[19]。Gretta等則率先利用該原理將翅片作為渦生成器增強流體的摻混[20]。Elavarasan等利用平面激光誘導(dǎo)熒光(PLIF)技術(shù)觀測了平板梯形翅片產(chǎn)生的發(fā)卡渦及其破裂后的對轉(zhuǎn)渦結(jié)構(gòu)[21]。Yang等采用粒子成像測速(PIV)技術(shù)進(jìn)行了類似的研究,同樣觀測到了反向渦對[19]。Lin研究了翅片式渦生成器的高度、傾斜角等幾何參數(shù)對其邊界層分離控制的影響[22]。盡管關(guān)于翅片式渦生成器的研究比較廣泛,但大多集中在高效換熱領(lǐng)域[17],而且多與平板結(jié)構(gòu)搭配,其應(yīng)用于彎頭內(nèi)作為整流裝置的相關(guān)研究鮮見公開報道。在核電站、艦船等場景內(nèi),水路管網(wǎng)十分復(fù)雜,泵或流量計上游短距離內(nèi)極有可能布置有彎頭結(jié)構(gòu),往往需要在彎頭內(nèi)進(jìn)行整流工作。因此,研究翅片式彎管流動調(diào)整器具有十分重要的理論研究意義和工程應(yīng)用價值。

圖1 翅片式流動調(diào)整器整流原理示意圖[15]Fig.1 Rectification mechanism of tab-type flow conditioner[15]

本文以翅片式彎頭流動調(diào)整器為研究對象,采用3種梯度方向的速度分布作為進(jìn)口邊界條件,系統(tǒng)研究了該結(jié)構(gòu)相對光滑彎管在不同雷諾數(shù)下的整流能力,包括壓力均勻度、速度均勻度和流向角均勻度等,深入分析了其流場分布特征,揭示了其整流機理,總結(jié)了其壓力損失特性。希望通過本文研究揭示翅片在彎管內(nèi)的整流特性,為研制高效的彎管流動調(diào)整器提供參考。

1 數(shù)值方法

1.1 計算域設(shè)置

采用商用軟件ANSYS CFX進(jìn)行定常數(shù)值模擬,帶翅片計算域布置參照文獻(xiàn)[15],主要結(jié)構(gòu)包括抗渦肋、二次抗渦肋、上游翅片、下游翅片等,如圖2所示。其中上游翅片與抗渦肋為一體式結(jié)構(gòu),下游翅片則為離散的翅片組。計算域主要尺寸見表1,圓管內(nèi)徑D為154.08 mm,上游進(jìn)口段長度Lu=1D,下游發(fā)展段長度Ld=10D。對照計算域為相同尺寸的光滑管。上下游翅片長度L均為0.5D,寬度W均為0.4D,收縮角α均為15°,抬升角β分別為50°、40°,其他幾何參數(shù)詳見文獻(xiàn)[15]。

(a)三維圖

表1 計算域主要幾何參數(shù)

1.2 邊界條件

基于CAP1400核電機組的管路流動參數(shù),表2列出了本文研究的邊界條件。模擬工質(zhì)為25 ℃水,進(jìn)口速度與6×105～1.2×106的雷諾數(shù)相匹配,出口壓力為0.18 MPa。進(jìn)口設(shè)置了3種均勻線性速度梯度,用于模擬最惡劣的彎管進(jìn)口不均勻速度分布,如圖3所示。圖中,變量VR是無量綱速度,定義為當(dāng)?shù)厮俣扰c進(jìn)口平均速度之比。相關(guān)算例按照梯度方向簡寫為PY、PZ、NZ,同時,光滑管結(jié)構(gòu)和流動調(diào)整器結(jié)構(gòu)分別簡寫為SP和FC,如SP-PY代表進(jìn)口為+Y方向速度梯度的光滑管算例,而FC-NZ代表進(jìn)口為-Z方向速度梯度的流動調(diào)整器算例。值得說明的是,由于+Y與-Y方向的進(jìn)口線性速度分布在下游造成的流動狀態(tài)是沿XZ平面對稱的,為了避免冗余,本文沒有開展-Y方向速度梯度的相關(guān)研究。

表2 邊界條件參數(shù)Table 2 Parameters of boundary conditions

(a)SP-PY

1.3 湍流模型驗證

現(xiàn)有文獻(xiàn)中關(guān)于翅片式流動調(diào)整器的實驗研究數(shù)據(jù)較少,本文以水為工質(zhì),對雙扭管(兩個連續(xù)彎管)下游布置Zanker式流動調(diào)整器的實驗[23]進(jìn)行湍流模型驗證,其流體域如圖4所示。采用定常數(shù)值模擬常用的k-ε、RNGk-ε、k-ω和SSTk-ω模型進(jìn)行湍流模型驗證,驗證結(jié)果如圖5a所示。對比工況的雷諾數(shù)為6.0×105,對比數(shù)據(jù)為無量綱軸向速度,其中Vb為充分發(fā)展段的流動速度,r為流動調(diào)整器下游3D截面處相對管道軸心的高度,以管道軸心為界,下側(cè)為負(fù),上側(cè)為正?？梢钥闯?相對于實驗結(jié)果,數(shù)值模擬得到的速度分布曲線要更平滑一些;4種湍流模型的計算結(jié)果中,SSTk-ω模型的結(jié)果與實驗結(jié)果最接近,最大偏差為13.27%,平均偏差2.71%。由于壁面附近速度梯度較高,呈現(xiàn)較強的各項異性特征,因此采用各向同性渦黏性假設(shè)的雷諾時均方法求解邊界層時會存在一定的誤差。SSTk-ω湍流模型得到的誤差在可以接受的范圍內(nèi),因此,本文采用SSTk-ω模型進(jìn)行數(shù)值仿真研究。

圖4 湍流模型驗證計算域Fig.4 Computational domain for turbulence model validation

(a)速度分布預(yù)測結(jié)果與實驗數(shù)據(jù)對比(Re=6.0×105)

為了檢驗所選湍流模型對壓力損失的預(yù)測精度,對比了不同雷諾數(shù)下數(shù)值模擬所得壓力損失系數(shù)與實驗結(jié)果[23]的差異見圖5b。定義壓力損失系數(shù)

(1)

式中:ΔP為流動調(diào)整器前后壓力差;ρ為流體密度。由圖5b可以看出,數(shù)值模擬的壓力損失系數(shù)略高于實驗結(jié)果,最大偏差3.72%,平均偏差1.88%,相對誤差較小,說明所用數(shù)值方法可以較準(zhǔn)確地預(yù)測流動調(diào)整器的壓力損失系數(shù)。

1.4 網(wǎng)格無關(guān)性驗證

由于翅片式流動調(diào)整器包含多個方向的翅片,結(jié)構(gòu)較為復(fù)雜,因此對計算域采用非結(jié)構(gòu)化網(wǎng)格剖分。近壁面第一層網(wǎng)格厚度0.002 mm,增長率為1.2,邊界層層數(shù)為30。經(jīng)計算,壁面Y+為1.0以下,符合所采用的SST湍流模型要求。壁面網(wǎng)格和橫截面網(wǎng)格如圖6所示。本節(jié)采用5組精細(xì)度的網(wǎng)格進(jìn)行網(wǎng)格無關(guān)性驗證,按照壁面和主流等比例加密的原則進(jìn)行加密,網(wǎng)格數(shù)量分別設(shè)置為200萬、300萬、400萬、500萬、600萬。驗證結(jié)構(gòu)為FC-PY,進(jìn)口雷諾數(shù)為9×106,對比數(shù)據(jù)為彎管下游5個位置處的速度均勻度,驗證結(jié)果如圖7所示。

(a)壁面網(wǎng)格

從圖7可以看到,在網(wǎng)格數(shù)量超過400萬后,彎頭下游5個截面數(shù)據(jù)點的速度均勻度相對誤差很小,500萬與600萬網(wǎng)格所得到的預(yù)測結(jié)果的最大相對誤差為0.4%。為了平衡計算精度要求與計算資源消耗,本文所有算例均采用500萬網(wǎng)格對應(yīng)的加密設(shè)置,即最大體網(wǎng)格限定10 mm,翅片表面加密至1.5 mm。

圖7 網(wǎng)格數(shù)量對彎管下游速度均勻度的影響Fig.7 Effect of grid number on velocity uniformity

2 流場均勻度評價方法

流場均勻度評價采用3個參數(shù),即壓力均勻度、速度均勻度和流向角均勻度。數(shù)據(jù)采集方法為彎頭下游1D～5D以1D距離為間隔設(shè)置截面,在截面上以點陣的方式收集當(dāng)?shù)貕毫?、速度和流向角。?shù)據(jù)點分布示意如圖8所示,共8條輻線,每條輻線上有10個數(shù)據(jù)點,加上原點共81個采樣點。

圖8 截面數(shù)據(jù)點分布示意圖Fig.8 Distribution of data sampling points

壓力均勻度和速度均勻度的計算公式[24]為

(2)

(3)

其中Vj為第j個采樣點的壓力或速度,n為采樣點個數(shù)。通過比較不同工況下的C值來評判截面上流場的均勻度,C值越大,流場均勻度越高。

根據(jù)GB-T 17611—1998,旋渦角θ為橫截面某給定點的局部速度與管道軸線之間的夾角?？紤]到無量綱處理要求分母不為0,定義旋渦角的余角為流向角θ,其計算方法為

(4)

式中:Vw代表給定點速度在管道軸線方向上的分量;V代表給定點的速度。

流向角均勻度的定義為

(5)

式中:θj為第j個采樣點的流向角。Cθ越大,流場均勻度越高。

3 結(jié)果與分析

本文算例共兩種幾何結(jié)構(gòu),即光滑管與翅片式彎頭流動調(diào)整器,進(jìn)口速度共5種,用于匹配5種雷諾數(shù)(6×105、7.5×105、9×105、1.05×106和1.2×106),進(jìn)口速度梯度方向共3種。計算結(jié)果顯示雷諾數(shù)對于兩種結(jié)構(gòu)各自的流場特性、均勻度等沒有明顯影響,為避免結(jié)果繁雜,除了3.3節(jié)關(guān)于壓力損失分析外,其他部分僅展示進(jìn)口雷諾數(shù)為9×105的結(jié)果。

3.1 不同進(jìn)口速度分布下的整流性能對比

壓力均勻度影響管路的穩(wěn)定情況,速度均勻度體現(xiàn)了流體的流量分布情況,流向角均勻度體現(xiàn)了流體流動的旋渦尺度大小。通過這3個參數(shù)可以系統(tǒng)地驗證流動調(diào)整器在不同邊界條件下的整流性能及適用范圍。

圖9是流動調(diào)整器對彎頭下游壓力均勻度的影響。可以看到,光滑管彎頭下游的壓力均勻度非常高,3種進(jìn)口邊界條件中,只有PY方向的速度分布產(chǎn)生了略低的壓力均勻度,且在5D范圍內(nèi)呈現(xiàn)持續(xù)下降的趨勢。相比之下,增加流動調(diào)整器后,壓力均勻度均在研究范圍內(nèi)呈上升趨勢,數(shù)值上與光滑管彎頭產(chǎn)生的壓力均勻度在同一區(qū)間,均保持在98%以上?？梢哉f,在壓力均勻度方面,增加翅片式彎頭流動調(diào)整器沒有顯著的影響,但其產(chǎn)生的均勻度完全符合測量要求。

圖9 流動調(diào)整器對壓力均勻度的影響(Re=9×105) Fig.9 Effect of flow conditioner on pressure uniformity(Re=9×105)

圖10是流動調(diào)整器對速度均勻度的影響?？梢钥吹?在3種邊界條件下,流動調(diào)整器均能在彎頭下游3D范圍內(nèi)將速度均勻度提升至85%以上,在5D范圍內(nèi)提升至90%以上。在下游3D處,流動調(diào)整器產(chǎn)生的速度均勻度相對于光滑彎頭在PY、PZ、NZ條件下分別提升了31%、21%、36%。其中,PY、PZ條件下的光滑彎頭下游速度均勻度提升較快,而NZ條件下的速度均勻度始終在65%以下,這體現(xiàn)了彎頭整流相對于直管整流的特殊性,即進(jìn)口速度梯度方向的改變會對下游流動產(chǎn)生很大的影響。同時,也體現(xiàn)了翅片式流動調(diào)整器對于彎頭進(jìn)口不均勻來流的整流能力是顯著而全面的。

圖10 流動調(diào)整器對速度均勻度的影響(Re=9×105) Fig.10 Effect of flow conditioner on velocity uniformity(Re=9×105)

圖11為流動調(diào)整器對流向角均勻度的影響。可以看出,在PZ和NZ條件下,光滑彎頭下游流向角均勻度全程高于90%,而同樣條件下的翅片式流動調(diào)整器反而降低了彎頭下游的均勻度。在PY條件下,兩種結(jié)構(gòu)在彎頭下游產(chǎn)生的流向角均勻度均在80%附近,流動調(diào)整器略微提升了均勻度。這說明,在多數(shù)情況下,翅片式流動調(diào)整器對于削弱彎頭下游的旋渦強度沒有明顯的作用,若要在短距離內(nèi)實現(xiàn)無旋渦流動則需在下游直管段搭配管束式整流器等其他裝置。

圖11 流動調(diào)整器對流向角均勻度的影響(Re=9×105)Fig.11 Effect of flow conditioner on flow angle uniformity(Re=9×105)

3.2 流場特性分析

圖12、圖13分別展示的是PZ條件下兩種彎管內(nèi)Q準(zhǔn)則等值面和下游渦量分布,其中ωz為渦量。Q是流體仿真中用于捕捉渦核位置的常用參數(shù),其定義為

(6)

圖12 光滑彎管內(nèi)的Q準(zhǔn)則等值面和下游渦量分布(SP-PZ,Re=9×105)Fig.12 Q-criterion iso-surface and vorticity distribution in smooth pipe(SP-PZ,Re=9×105)

圖13 翅片式流動調(diào)整器內(nèi)的Q準(zhǔn)則等值面和下游渦量分布(FC-PZ,Re=9×105)Fig.13 Q-criterion iso-surface and vorticity distribution in tab-type flow conditioner(FC-PZ,Re=9×105)

在均勻來流下,彎管下游一般會出現(xiàn)對轉(zhuǎn)形態(tài)的“迪恩渦”。然而,在PZ條件下,彎頭內(nèi)側(cè)速度高,在彎頭下游誘發(fā)了兩組流向渦對,體現(xiàn)了不均勻來流對流動結(jié)構(gòu)的影響。增設(shè)翅片后,其上游翅片前后面間的壓力差會在翅片邊緣誘導(dǎo)產(chǎn)生指向彎頭內(nèi)側(cè)的對轉(zhuǎn)渦對,在彎頭下游1D距離內(nèi)與其下游翅片誘發(fā)的對轉(zhuǎn)渦對相互摻混融合,產(chǎn)生新的流向渦。在這一過程中,原始線性速度分布被擾亂,在彎頭下游1D后逐漸形成新的速度分布。

圖14為光滑管內(nèi)3種進(jìn)口條件下的流線分布。在PY條件下,高速流體和低速流體在進(jìn)入彎頭后在壓力梯度的驅(qū)使下形成大尺度的流向渦,并在彎頭下游始終保持該形態(tài),可以看到,管壁邊界層的剪切和流體之間的動量交換在一定程度上削弱了這種速度差,但仍然需要較長的直管段才能完全消除。在PZ條件下,彎頭內(nèi)側(cè)的高速流體和彎頭外側(cè)的低速流體同樣會在壓力梯度的作用下產(chǎn)生動量交換,并由此形成復(fù)雜的小尺度流向渦,其速度不均勻程度在這一過程中得到削弱。在NZ條件下,原始速度差產(chǎn)生的壓力梯度可以在一定程度上抵消流經(jīng)彎頭時產(chǎn)生的徑向壓力梯度,因此在經(jīng)過彎頭后,原始速度分布得以保留,沒有產(chǎn)生宏觀尺度的流向渦。

(a)SP-PY

圖15是不同進(jìn)口條件下翅片式流動調(diào)整器內(nèi)的流線分布(Re=9×105)。在PY條件下,彎頭內(nèi)的單螺旋流向渦被翅片誘導(dǎo)渦攪拌干擾,在彎頭下游1D位置前充分摻混,有效降低了速度差,但其“以渦擾渦”的原理使其無法削弱流向渦的強度。在PZ條件下,彎頭外側(cè)的低速流體在流經(jīng)上游翅片時產(chǎn)生了指向彎頭內(nèi)側(cè)的誘導(dǎo)渦對,但由于其動量和強度較低,無法進(jìn)行高低速流體間的充分摻混,需要下游翅片誘導(dǎo)渦的補充才能完成動量交換,同時會附帶產(chǎn)生低強度的流向渦。在NZ條件下,彎頭外側(cè)的高速流體在流經(jīng)上游翅片時產(chǎn)生了高強度的誘導(dǎo)渦對,沖擊到彎頭內(nèi)側(cè)壁面,與該位置處的低速流體進(jìn)行充分動量交換。在經(jīng)過彎頭后,依然有部分流體的動量難以被翅片有效削弱。

(a)FC-PY

圖16和圖17分別給出了彎頭下游不同位置截面上Z方向渦量分量和無量綱速度的分布情況。渦量分布方面,對于光滑管來說,在PY條件下,彎頭下游存在明顯的Z方向單螺旋流向渦,在研究范圍內(nèi)沒有削弱。在PZ條件下,彎頭內(nèi)外側(cè)流體的充分摻混形成了有對稱特征的兩組對轉(zhuǎn)渦,且維持了較遠(yuǎn)距離。在NZ條件下,原始速度分布在經(jīng)過彎頭后得以保留,沒有劇烈的動量交換過程,因此沒有產(chǎn)生大尺度流向渦。對翅片式流動調(diào)整器來說,在1D位置處可以看到下游翅片形成的多組誘導(dǎo)對轉(zhuǎn)渦,然而其位置和形態(tài)受到主流速度分布的影響。在PY條件下,單螺旋流向渦的周向動量使下游翅片誘導(dǎo)渦對沿順時針方向偏離了翅片位置,經(jīng)過動量交換后,下游仍存在單螺旋渦。在PZ條件下,經(jīng)過整流后彎頭下游僅存在一組流向?qū)D(zhuǎn)渦。在NZ條件下,上游翅片誘導(dǎo)渦對彎頭內(nèi)側(cè)的強烈沖擊使得下游翅片誘導(dǎo)渦對向彎頭外側(cè)方向偏離。經(jīng)過摻混后,下游仍存在一對大尺度“迪恩渦”。

(a)光滑管

(a)光滑管

無量綱速度VR分布方面,對光滑管來說,在PY條件下,高速流體沖向管壁并沿壁面呈周向運動,以單螺旋渦包裹內(nèi)部的低速渦核。在PZ條件下,彎頭內(nèi)側(cè)高速流體沿壁面以對轉(zhuǎn)渦的形態(tài)從兩側(cè)包裹內(nèi)部低速流體,并進(jìn)一步演變?yōu)閮山M對轉(zhuǎn)渦。在NZ條件下,由于沒有形成流向渦,高低速流體在彎頭下游保持進(jìn)口速度分布,缺乏劇烈的動量交換過程。對于彎頭來說,來流速度梯度方向的影響會與彎頭方向的影響產(chǎn)生耦合疊加,因而不均勻來流產(chǎn)生的不均勻度情況要比在直管更復(fù)雜多變,更難以消除。對翅片式流動調(diào)整器來說,在PY條件下,單螺旋流向渦產(chǎn)生的速度差在下游翅片誘導(dǎo)渦的有效摻混下大幅降低,彎頭下游3D距離后的速度分布比較均勻。在PZ條件下,管壁附近的流體在經(jīng)過整流后相對于管道中心的流體有一定的速度差,且在彎頭下游3D～5D距離呈現(xiàn)快速降低趨勢。在NZ條件下,彎頭外側(cè)的高速流體在流經(jīng)上下游翅片后依然具有相當(dāng)?shù)膭恿?并在下游彎頭外側(cè)翅片誘導(dǎo)渦的裹挾下在管道中心形成橫向高速帶。盡管不能在5D距離內(nèi)完全消除速度差,但經(jīng)過整流后的速度差在各個進(jìn)口條件下、各個下游位置處均明顯優(yōu)于光滑管。

3.3 不同進(jìn)口雷諾數(shù)下的壓力損失對比

流動調(diào)整器在改變流體速度分布的同時會難以避免地造成流動損失,圖18為在所研究的雷諾數(shù)范圍內(nèi),翅片式彎頭流動調(diào)整器產(chǎn)生的壓力損失。根據(jù)文獻(xiàn)[25],不均勻來流的影響需要200D直管才能徹底消除。因此,圖中還對比了充分發(fā)展流動在200D直管內(nèi)產(chǎn)生的壓力損失,其計算公式為

(7)

式中:l0為管內(nèi)流動距離;f為摩擦系數(shù),參照下式確定。

(8)

從圖18可以看出,隨著進(jìn)口雷諾數(shù)的增加,各種工況的壓力損失都有所提高,其中,PZ條件下的翅片式流動調(diào)整器產(chǎn)生的壓力損失最低,其他兩種進(jìn)口速度分布產(chǎn)生的壓力損失情況幾乎完全相同。一方面,由于翅片在管道內(nèi)的投影面積較大,其引起的壓差損失要比誘導(dǎo)渦引起的摻混損失大很多。在PY和NZ條件下,高速流體會在彎頭內(nèi)沖擊翅片,造成翅片正面較大的滯止壓力,而在PZ條件下,沖擊翅片的是低速流體,在翅片正面產(chǎn)生的滯止壓力較低,其引起的誘導(dǎo)渦進(jìn)一步阻止了高速流體向彎頭外側(cè)翅片流動,所以PZ條件下的壓力損失明顯低于其他兩種條件。另一方面,盡管PY條件產(chǎn)生的壓差阻力低于NZ條件,但其二次流結(jié)構(gòu)更復(fù)雜更不均勻,因而其摻混損失要高于NZ條件,總壓力損失與NZ條件接近。在研究范圍內(nèi)的各種工況下,翅片式彎管流動調(diào)整器產(chǎn)生的流動損失都遠(yuǎn)低于充分發(fā)展的流動在200D的直管內(nèi)所造成的壓力損失。

圖18 彎頭流動調(diào)整器與普通直管壓損比較Fig.18 Comparison of pressure losses between elbow flow conditioner and straight pipe

表3橫向?qū)Ρ攘顺崞綇濐^流動調(diào)整器與前文提到的其他類型典型流動調(diào)整器的壓力損失系數(shù)。需要說明的是,流動調(diào)整器的壓力損失系數(shù)隨雷諾數(shù)變化略有變化,但一般而言幅度不大,因此表中翅片式流動調(diào)整器的壓力損失系數(shù)為各個雷諾數(shù)下的壓損系數(shù)平均值。其他流動調(diào)整器的壓力損失系數(shù)來自文獻(xiàn)[5]和文獻(xiàn)[23]。從表3可以看出,翅片式彎頭流動調(diào)整器的壓力損失系數(shù)隨進(jìn)口條件的變化有較大幅度的變化,為0.36～1.17,但基本與管束式流動調(diào)整器相當(dāng),高于簡單的八葉片式流動調(diào)整器,低于大流阻的孔板式流動調(diào)整器。

表3 各種流動調(diào)整器無量綱壓力損失比較Table 3 Comparison of pressure losses among different flow conditioners

4 結(jié) 論

(1)翅片式彎頭流動調(diào)整器的主要作用機理在于利用翅片前后的壓差在翅片邊緣產(chǎn)生對轉(zhuǎn)的誘導(dǎo)渦,使部分彎頭外側(cè)的流體沖擊彎頭內(nèi)側(cè),增強動量交換,擾亂來流速度分布,并利用下游翅片誘導(dǎo)渦加強摻混,降低流體間的速度差。

(2)翅片式彎頭流動調(diào)整器對彎頭下游壓力均勻度影響不大,對速度均勻度有21%～36%的提高,對流向角均勻度的影響隨進(jìn)口速度梯度方向的變化而變化,一般呈負(fù)面效果。

(3)在泵或流量計前存在彎頭的情況下,來流速度梯度方向的影響會與彎頭壓力梯度的影響產(chǎn)生耦合疊加,因而來流產(chǎn)生的不均勻度情況要比在直管內(nèi)更復(fù)雜多變,更難以消除。

(4)翅片式彎頭流動調(diào)整器適用于削弱不均勻來流對彎頭下游流動結(jié)構(gòu)的影響,能夠在下游短距離內(nèi)構(gòu)建均勻的速度分布,但難以消除直管段內(nèi)的流向渦。

(5)翅片式彎頭流動調(diào)整器的壓力損失系數(shù)隨進(jìn)口條件的變化在0.36～1.17間變化,基本與管束式流動調(diào)整器相當(dāng),高于簡單的八葉片式流動調(diào)整器,低于大流阻的孔板式流動調(diào)整器。