不同出口角度擴(kuò)壓器的內(nèi)部流動及離心壓縮機(jī)級性能數(shù)值研究

趙志偉, 蘇永升(華東理工大學(xué)承壓系統(tǒng)與安全教育部重點實驗室,上海 200237)

不同出口角度擴(kuò)壓器的內(nèi)部流動及離心壓縮機(jī)級性能數(shù)值研究

趙志偉, 蘇永升
(華東理工大學(xué)承壓系統(tǒng)與安全教育部重點實驗室,上海 200237)

以某多級離心壓縮機(jī)首級為研究對象,運用數(shù)值模擬方法研究擴(kuò)壓器出口安裝角對壓縮機(jī)級性能的影響,得出不同出口安裝角情況下壓縮機(jī)單級的性能曲線,并分析擴(kuò)壓器內(nèi)部流動特點和損失機(jī)理,通過壓力恢復(fù)系數(shù)Cp對比不同流量下不同擴(kuò)壓器的擴(kuò)壓效果。結(jié)果表明,當(dāng)擴(kuò)壓器出口安裝角增大時,性能曲線向大流量區(qū)移動,最高效率和壓比先升高后降低;不同出口角度下葉片擴(kuò)壓器的擴(kuò)壓效果由不同工況下的流動特性決定;在大流量下,不同出口角度下擴(kuò)壓器葉背形成分離區(qū)且旋渦位置不同;在小流量下,具有較大出口安裝角的擴(kuò)壓器葉腹率先出現(xiàn)分離區(qū)。

葉片擴(kuò)壓器; 出口安裝角; 性能曲線; 壓力恢復(fù)系數(shù)

離心壓縮機(jī)由于結(jié)構(gòu)緊湊、工藝性好及性能高,在國民經(jīng)濟(jì)中占有重要地位,長期以來很多學(xué)者對離心壓縮機(jī)的氣動特性進(jìn)行了研究[1]。在離心壓縮機(jī)所有靜止件中,擴(kuò)壓器屬于核心部件,對提高壓縮機(jī)效率和級壓比、改善最佳工況點十分重要[2]。針對擴(kuò)壓器的研究已有大量的文獻(xiàn)。周俊安等[3]研究帶葉片擴(kuò)壓器離心壓縮機(jī)模型級內(nèi)流,得出擴(kuò)壓器內(nèi)的流動損失和擴(kuò)壓效果由不同工況下的流動特性決定;席光等[4]通過數(shù)值與實驗研究表明葉片擴(kuò)壓器進(jìn)口安裝角對壓縮機(jī)流動有較大影響;王志恒等[5]以高壓比離心壓縮機(jī)級性能為研究對象,研究了擴(kuò)壓器進(jìn)口安裝角和葉型對級性能的影響。李學(xué)臣等[6]采用數(shù)值模擬對多級離心壓縮機(jī)中間級的流動進(jìn)行了研究,重點考察進(jìn)口段、彎道和回流器對中間級氣動性能的影響;Anish等[7]采用數(shù)值模擬方法分析了不同擴(kuò)壓器、葉輪與擴(kuò)壓器不同間隙情況下離心壓縮機(jī)的性能;Zheng等[8]建立理論分析模型準(zhǔn)確地預(yù)測出擴(kuò)壓器的阻塞流量,證明隨著進(jìn)口預(yù)旋角增大,擴(kuò)壓器的阻塞流量比葉輪阻塞流量下降緩慢。Bousquet等[9]通過數(shù)值模擬質(zhì)量流量減少對葉輪和擴(kuò)壓器進(jìn)口的影響,提出詳細(xì)的失速機(jī)理。

本文以多級離心壓縮機(jī)首級作為研究對象,采用流體計算軟件CFX數(shù)值模擬擴(kuò)壓器在不同出口角度情況下的級性能曲線,對比擴(kuò)壓器壓力恢復(fù)系數(shù)Cp,詳細(xì)分析了葉片擴(kuò)壓器出口安裝角對壓縮機(jī)內(nèi)部流場特點、損失機(jī)理及級性能的影響。

1 幾何模型與網(wǎng)格

1.1 幾何模型

離心壓縮機(jī)首級由進(jìn)口段、后彎閉式葉輪、葉片擴(kuò)壓器、彎道和回流器組成,葉輪葉片和擴(kuò)壓器葉片均采用等厚度型,葉輪及擴(kuò)壓器流道結(jié)構(gòu)如圖1所示,設(shè)計參數(shù)如表1和表2所示,其中葉輪轉(zhuǎn)速為8 000 r/min,進(jìn)口總溫度為20 ℃,進(jìn)口總壓為0.101 MPa,介質(zhì)為空氣,相對寬度為葉片出口寬度與葉輪外徑之比,β2A為葉輪葉片出口安裝角,α3A和α4A分別為擴(kuò)壓器進(jìn)口、出口安裝角。

圖1 葉輪及擴(kuò)壓器流道結(jié)構(gòu)

Outerdiameter/mmNumbersofbladeRelativewidthβ2A/(°)600200．07345

表2 葉片擴(kuò)壓器設(shè)計參數(shù)

1.2 計算網(wǎng)格及無關(guān)性驗證

為全面考察擴(kuò)壓器出口安裝角對上下游部件流動的影響,采用整級數(shù)值模擬,所選計算區(qū)域包括進(jìn)口段、葉輪、葉片擴(kuò)壓器、彎道和回流器。采用ICEM前處理軟件,生成全場四面體非結(jié)構(gòu)化網(wǎng)格。通過選用網(wǎng)格拓?fù)淠K實現(xiàn)網(wǎng)格拓?fù)涞淖詣庸芾?、網(wǎng)格優(yōu)化、網(wǎng)格塊的附著和關(guān)聯(lián)等功能,提高網(wǎng)格的質(zhì)量和效率。同時加密靠近壁面的網(wǎng)格以滿足湍流模型對y+的要求。本文采用全流道計算,計算域及網(wǎng)格如圖2所示。

在數(shù)值模擬計算時,所有模型級的網(wǎng)格數(shù)均滿足無關(guān)性要求。圖3示出了在質(zhì)量流量為7.25 kg/s時擴(kuò)壓器Ⅱ的網(wǎng)格無關(guān)性驗證。由圖3可知級效率(η)與壓比(ε)變化不明顯,認(rèn)為網(wǎng)格達(dá)到無關(guān)性?？紤]到求解精度和計算速度,采用表3所示的擴(kuò)壓器網(wǎng)格計算分析,整級的網(wǎng)格數(shù)大于2.7×106。

圖2 計算域及網(wǎng)格

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

No．NumbersofgridNodesⅠ1175821204597Ⅱ1178060204853Ⅲ1179426205006Ⅳ1179040204852

2 數(shù)學(xué)模型與邊界條件

2.1 計算模型

應(yīng)用流體計算軟件CFX來求解三維可壓縮雷諾平均守恒型N-S方程。湍流模型選取κ-ε兩方程模型[10],壁面邊界層效應(yīng)采用壁面函數(shù)法。為保證求解準(zhǔn)確,采用高階求解模式。數(shù)值模擬控制方程可表示為方程(1)通用形式:

(1)

式中:ρ為密度;V為體積;t為時間;φ為通用變量;Γ為廣義擴(kuò)散系數(shù);S為廣義源項。

2.2 邊界條件及收斂標(biāo)準(zhǔn)

本文計算過程分為定常計算和非定常計算兩個階段。為有效減少流動不充分所帶來的誤差,進(jìn)口邊界設(shè)置為進(jìn)口段入口,出口邊界為回流器出口。定常計算的計算域入口均為給定總溫、總壓和軸向速度方向;出口為質(zhì)量流量;固體壁面均為絕熱、無滑移條件。非定常計算以定常計算結(jié)果為初始流場,柯朗數(shù)CFL取2,時間步長設(shè)置為2.3×10-4s,內(nèi)迭代設(shè)置為20步。計算殘差下降到10-4數(shù)量級以下,當(dāng)進(jìn)出口流量誤差小于0.1%,其他監(jiān)控參數(shù)(效率和壓比)保持恒定時,即認(rèn)為計算收斂。

3 數(shù)值模擬結(jié)果分析

3.1 性能曲線

圖4、圖5分別示出了不同質(zhì)量流量(Qm)下級壓比和多變效率曲線。由于當(dāng)質(zhì)量流量小于5.75 kg/s時,整級模擬不能收斂,所以數(shù)值模擬從質(zhì)量流量為5.75 kg/s開始計算;當(dāng)質(zhì)量流量大于5.75 kg/s并小于7.00 kg/s時,定常計算難以收斂,此時以定常計算結(jié)果為初始流場,采用非定常計算。非定常計算結(jié)果在本文中采用時均值處理[11]。

圖4 壓比曲線

圖5 效率曲線

如圖4、圖5所示,隨著擴(kuò)壓器出口安裝角度增大,級性能曲線向大流量方向偏移,雖最高效率和壓比先升高后降低，但它們的最高值提升并不明顯,且最大工況點變化不大。在設(shè)計質(zhì)量流量7.25 kg/s附近工況點,帶有4種不同出口安裝角的擴(kuò)壓器的壓縮機(jī)級壓比、效率接近;在質(zhì)量流量為6.50 kg/s時,擴(kuò)壓器Ⅳ (α4A=36.0°)的離心壓縮機(jī)級壓比、效率明顯下降;在質(zhì)量流量為8.15 kg/s時,擴(kuò)壓器Ⅰ (α4A=29.0°)的級壓比、效率都較低。通過性能曲線對比可以得出:當(dāng)擴(kuò)壓器進(jìn)口安裝角度一定時,選取合適的出口安裝角度既可以增大工作范圍,又可以改善壓縮機(jī)的工作效率。

3.2 擴(kuò)壓器內(nèi)部流動分析

3.2.1 概述 本文采用3種工況,即質(zhì)量流量分別取6.50、7.25、8.15 kg/s,分析擴(kuò)壓器內(nèi)部的流動情況和壓力損失。

3.2.2 擴(kuò)壓器等半徑截面上Cp分析 擴(kuò)壓器的作用是將動能轉(zhuǎn)化為靜壓力能。當(dāng)考察擴(kuò)壓器性能時,采用靜壓恢復(fù)系數(shù)Cp[12]也是一種判定方法。定義靜壓恢復(fù)系數(shù)為

(2)

其中:Ps3,Psx分別為擴(kuò)壓器葉片前緣進(jìn)口、等截面的靜壓平均值;Pt3為擴(kuò)壓器葉片前緣進(jìn)口總壓的平均值。

圖6～圖8示出了不同工況下擴(kuò)壓器靜壓恢復(fù)系數(shù)分布。

如圖6,當(dāng)質(zhì)量流量為6.50 kg/s時,不同出口安裝角度的擴(kuò)壓器在擴(kuò)壓器出口處都有較高的靜壓恢復(fù)系數(shù),但擴(kuò)壓器Ⅳ的靜壓恢復(fù)系數(shù)相比于其他擴(kuò)壓器較低,其次較低的是擴(kuò)壓器Ⅰ,擴(kuò)壓器Ⅱ和擴(kuò)壓器Ⅲ的靜壓恢復(fù)系數(shù)相近。

圖6 Qm=6.50 kg/s時擴(kuò)壓器靜壓恢復(fù)系數(shù)

當(dāng)質(zhì)量流量為7.25 kg/s時(如圖7),擴(kuò)壓器Ⅱ、Ⅲ、Ⅳ的靜壓恢復(fù)系數(shù)良好,而擴(kuò)壓器Ⅰ的靜壓恢復(fù)系數(shù)相比其他擴(kuò)壓器較低。此時,隨著擴(kuò)壓器出口安裝角度的增大,擴(kuò)壓器靜壓恢復(fù)系數(shù)增大。

當(dāng)質(zhì)量流量為8.15 kg/s時(如圖8),靜壓恢復(fù)系數(shù)出現(xiàn)較大差異,擴(kuò)壓器流道中靜壓恢復(fù)系數(shù)為負(fù)值。且隨著出口安裝角的增大，極小值向出口處偏移。在擴(kuò)壓器出口處，擴(kuò)壓器Ⅰ的靜壓恢復(fù)系數(shù)為負(fù)值,失去擴(kuò)壓器功能；擴(kuò)壓器Ⅱ的靜壓恢復(fù)系數(shù)為零,也失去擴(kuò)壓功能,擴(kuò)壓器Ⅲ、Ⅳ的靜壓恢復(fù)系數(shù)為正值,有部分?jǐn)U壓功能。

綜合圖6～圖8可見,在小流量和設(shè)計流量時,擴(kuò)壓器靜壓恢復(fù)系數(shù)在前半段上升較快,后半段變平緩,說明擴(kuò)壓器的擴(kuò)壓功能主要集中在前半段;在大流量時,出口安裝角度小的擴(kuò)壓器率先失去擴(kuò)壓功能。

圖7 Qm=7.25 kg/s時擴(kuò)壓器靜壓恢復(fù)系數(shù)

圖8 Qm=8.15 kg/s時擴(kuò)壓器靜壓恢復(fù)系數(shù)

3.2.3 擴(kuò)壓器50%葉高截面處靜壓和流線分布 在穩(wěn)態(tài)流動時,流線的形狀和位置不隨時間變化,采用流線能直觀反映流動情況,同時靜壓云圖可直觀反映擴(kuò)壓器內(nèi)靜壓力分布。下文以擴(kuò)壓器Ⅰ、Ⅱ(與Ⅲ相近)、Ⅳ作為研究對象,分析了擴(kuò)壓器50%葉高截面處流體流動情況和壓力分布。當(dāng)質(zhì)量流量為6.50 kg/s時,由于計算結(jié)果為非定常計算所得,故小流量采取時均處理。

當(dāng)質(zhì)量流量為6.50 kg/s時(如圖9),擴(kuò)壓器葉片前緣沖角i=6°,擴(kuò)壓器內(nèi)流動出現(xiàn)差異。擴(kuò)壓器Ⅰ并未產(chǎn)生旋渦區(qū),而擴(kuò)壓器Ⅱ葉腹尾緣產(chǎn)生較小旋渦區(qū),但擴(kuò)壓器Ⅱ有較高的靜壓;擴(kuò)壓器Ⅳ葉腹尾緣處產(chǎn)生較大分離和回流,形成氣流旋渦,同時部分流道由于旋渦的存在靜壓分布不均勻。

當(dāng)質(zhì)量流量為7.25 kg/s時(如圖10),擴(kuò)壓器葉片前緣沖角i≈0°,氣流分布均勻且速度變化平緩,沒有明顯的氣流分離或回流,流動情況良好。此時擴(kuò)壓器Ⅱ、Ⅳ有較高的靜壓。

圖9 Qm=6.50 kg/s時50%葉高截面處靜壓和流線時均圖

圖10 Qm=7.25 kg/s時50%葉高截面處靜壓和流線圖

當(dāng)流量為8.15 kg/s時(如圖11),擴(kuò)壓器葉片前緣沖角i=-8°,擴(kuò)壓器Ⅰ、Ⅱ、Ⅳ葉背都產(chǎn)生明顯旋渦,而擴(kuò)壓器Ⅰ旋渦主要出現(xiàn)在葉背前緣,且在喉部流動變?yōu)槭諗啃再|(zhì);擴(kuò)壓器Ⅳ產(chǎn)生的旋渦位于葉背靠近出口處且旋渦偏小,擴(kuò)壓器Ⅱ介于兩者之間。擴(kuò)壓器靜壓都明顯下降,同時在葉片前緣處產(chǎn)生較大的沖擊損失。

圖11 Qm=8.15 kg/s時50%葉高截面處靜壓和流線圖

3.2.4 擴(kuò)壓器等半徑截面的總壓分布 擴(kuò)壓器內(nèi)的總壓分布反映內(nèi)部流動損失的分布,且等半徑截面的總壓分布能全面地反映出擴(kuò)壓器周向與徑向的流動損失分布。圖12～圖14示出了不同工況下擴(kuò)壓器Ⅰ、Ⅱ、Ⅳ沿氣體流向在半徑截面上的總壓分布,其量綱為一半徑分別為X1=0.2、X2=0.5和X3=0.8。圖中S、H、SS、PS分別表示擴(kuò)壓器的前蓋板、后蓋板、吸力面和壓力面。

當(dāng)質(zhì)量流量為6.50 kg/s時(如圖12),由于氣流進(jìn)氣角度偏小,氣流向葉背靠攏,在靠近葉背處產(chǎn)生較大的總壓,且輪蓋側(cè)高于輪盤側(cè),而靠近葉腹的總壓急劇減小。擴(kuò)壓器Ⅳ沿周向總壓梯度變化比較大且在X3截面上靠近葉腹的總壓低,產(chǎn)生流動損失,其次是擴(kuò)壓器Ⅱ、Ⅰ。

當(dāng)質(zhì)量流量為7.25 kg/s時(如圖13),擴(kuò)壓器各截面內(nèi)總壓梯度變化小,但也有差異;擴(kuò)壓器Ⅰ較其他二者總壓梯度變化大,擴(kuò)壓器Ⅱ在流道內(nèi)流動尤其平穩(wěn);擴(kuò)壓器Ⅱ、Ⅳ在X3截面的壓力已分布均勻。擴(kuò)壓器Ⅰ、Ⅱ、Ⅳ未產(chǎn)生總壓明顯減小的區(qū)域,說明此工況下擴(kuò)壓器內(nèi)的流動損失小。

當(dāng)質(zhì)量流量為8.15 kg/s時(如圖14),由于氣流進(jìn)氣角度偏大,氣流不斷向葉腹靠攏,在葉腹靠近輪蓋側(cè)總壓局部上升;而在葉背處總壓急劇減小產(chǎn)生分離,流動損失較大,總壓沿周向和徑向分布都不均勻。擴(kuò)壓器Ⅰ在X1截面葉背處存在明顯低壓區(qū),而擴(kuò)壓器Ⅱ在X2截面葉背處存在明顯低壓區(qū),擴(kuò)壓器Ⅳ在X2截面存在低壓區(qū)且隨著截面半徑的增大低壓區(qū)范圍擴(kuò)大。

圖13 Qm=7.25 kg/s時擴(kuò)壓器Ⅰ、Ⅱ、Ⅳ沿流向各半徑截面上的總壓分布圖

圖14 Qm=8.15 kg/s時擴(kuò)壓器Ⅰ、Ⅱ、Ⅳ沿流向各等半徑截面上的總壓分布圖

綜合圖12～14可知,小流量(6.50 kg/s)和設(shè)計流量(7.25 kg/s)下，擴(kuò)壓器Ⅰ、Ⅱ、Ⅳ總壓損失并不明顯,表明葉腹產(chǎn)生的旋渦對總壓影響并不大;而大流量(8.15 kg/s)下總壓損失較前兩種工況大,表明分離損失產(chǎn)生的旋渦對總壓損失的影響較大。

4 結(jié) 論

(1) 數(shù)值研究表明,當(dāng)擴(kuò)壓器葉片出口安裝角度增大時,級性能曲線向大流量方向移動,級最高效率和壓比先升高后降低;對某一流量應(yīng)存在最佳的出口安裝角度,使級性能最優(yōu)。

(2) 擴(kuò)壓器靜壓恢復(fù)系數(shù)在擴(kuò)壓器前半段上升較快,后半段變平緩。小流量下擴(kuò)壓器的擴(kuò)壓效果隨出口角度的增大先增強(qiáng)后減弱;設(shè)計流量下擴(kuò)壓器擴(kuò)壓效果隨出口角度的增大而增強(qiáng);大流量下擴(kuò)壓器中間段不起擴(kuò)壓效果,喪失部分?jǐn)U壓功能,且出口安裝角度小的擴(kuò)壓器率先失去擴(kuò)壓功能。研究表明不同出口安裝角葉片擴(kuò)壓器的擴(kuò)壓效果由不同工況下的流動特性決定。

(3) 就單級流動特性而言,小流量下出口安裝角度大的擴(kuò)壓器葉腹率先出現(xiàn)分離區(qū),葉腹出現(xiàn)的旋渦對總壓損失影響較小;大流量下不同出口角度的擴(kuò)壓器其葉背形成分離區(qū)且旋渦的位置不同,葉背出現(xiàn)的旋渦對總壓損失影響較大。

[1]孫濤,王毅等.變厚度葉片對離心壓縮機(jī)結(jié)構(gòu)強(qiáng)度和性能影響[J].大連理工大學(xué)學(xué)報,2012,52(1):29-34.

[2]徐玉梅,劉丹丹,謝俊良,等.寬窄擴(kuò)壓器失速機(jī)理模擬分析[J].輕工機(jī)械,2015,33(3):33-37.

[3]周俊安,劉立軍,肖萍,等.帶葉片擴(kuò)壓器離心壓縮機(jī)模型級內(nèi)流研究[J].西安交通大學(xué)學(xué)報,2015,49(9):1-6.

[4]席光,周莉,丁海萍,等.葉片擴(kuò)壓器進(jìn)口安裝角對離心壓縮機(jī)性能影響的數(shù)值與實驗研究[J].工程熱物理學(xué)報,2006,27(1):61-64.

[5]李學(xué)臣,席光,李心偉,等.離心壓縮機(jī)中間級內(nèi)流動數(shù)值研究[J].工程熱物理學(xué)報,2010,31(6):955-958.

[6]王志恒,張鯤羽,王宏亮,等.葉片擴(kuò)壓器對離心壓氣機(jī)整機(jī)性能影響的實驗及數(shù)值研究[J].工程熱物理學(xué)報,2015,36(10):84-88.

[7]ANISH S,SITARAM N.Computational investigation of impeller-diffuser interaction in a centrifugal compressor with different type of diffusers [J].Proceeding of the Institution of Mechanical Engineers:Part A.Journal of Power and Energy,2009,223(A2):167-178.

[8]ZHENG Xinqian,HUANG Qiangqiang,LIU Anxiong.Loss mechanisms and flow control for improved efficiency of a centrifugal compressor at high inlet prewhirl[J].Journal of Turbomachinery,2016,138(10):101011-1-11.

[9]BOUSQUET Y,BINDER N.Numerical simulation of stall inception mechanisms in a centrifugal compressor with vaned diffuser[J].Journal of Turbomachinery,2016,138(12):121005-1-9.

[10]陶文銓.數(shù)值傳熱學(xué)[M].第2版.西安:西安交通大學(xué)出版社,2011:347-353.

[11]BENICHOU E,TRéBINJAC I.Comparison of steady and unsteady flows in a transonic radial vaned diffuser[J].Journal of Turbomachinery,2016,138(12):121002-1-10.

[12]FILIPENCO V G,DENIZ S,JOHNSTON J M,etal.Effects of inlet flow field conditions on the performance of centrifugal compressor diffusers:Part 1.Discrete-passage diffuser [J].Journal of Turbomachinery-Transactions of the Asme,2000,122 (1):1-10.

Numerical Study of the Different Outlet Angles of Diffuser and Performance of Centrifugal Compressor

ZHAO Zhi-wei, SU Yong-sheng

(Key Laboratory of Pressurized Systems and Safety,Ministry of Education,East China University of Science and Technology,Shanghai 200237,China)

Based on the first stage of centrifugal compressor ,the effects of outlet angles of diffuser on the stage performance were investigated by the numerical simulation.The performance curves at stage with different outlet angle were measured.The internal flow and the total pressure loss on the stage were analyzed,and by comparison of pressure recovery coefficient of different vaned diffusers to judge diffusion effect.The results show that the increase of outlet angle makes the performance curve move toward the large flow rate side,and the peak stage efficiency and pressure ratio increase firstly and then decrease.The diffusion effect of vaned diffuser depends upon flow condition during different operations.When the flow mass is larger than the design level,the separation zone generates near the shroud of diffuser,and the different vortex regions enlarge with the outlet angle of diffuser.When the flow mass is smaller than the design level,the separation zone is firstly emerged in large outlet angle at vaned diffuser.

vane diffuser; outlet angles; performance map; pressure recovery coefficient

1006-3080(2017)02-0266-07

10.14135/j.cnki.1006-3080.2017.02.019

2016-06-21

趙志偉(1989-),男,山西人,碩士生,從事流體機(jī)械方向研究。E-mail:zhaozhiwei012@126.com

蘇永升,E-mail:yssu@ecust.edu.cn

TH452

A