周俊安,劉立軍,肖萍,陳彥峰,郭江峰,陳雪江

?

帶葉片擴壓器離心壓縮機模型級內(nèi)流研究

周俊安1,劉立軍1,肖萍2,陳彥峰2,郭江峰1,陳雪江1

離心式壓縮機廣泛應用于冶金、石油化工、天然氣輸送、制冷以及動力工業(yè)等多個行業(yè),長期以來國內(nèi)外專家學者對離心壓縮機進行了研究和探索[1-4],其中模型級相似設計是目前工程中應用最為廣泛的一種設計方法。針對模型級的研發(fā)和研究工作已有大量的文獻,例如:聞蘇平等分析了馬赫數(shù)對離心壓縮機小流量系數(shù)模型級性能的影響[5];張楚華等研究了半開式離心壓縮機葉輪的內(nèi)部流場、氣動性能隨葉頂間隙的變化規(guī)律[6];劉瑞韜等研究了壓縮丙烯時模型級內(nèi)部及彎道回流器內(nèi)部流動規(guī)律[7];么立新對小流量模型級進行了結(jié)構(gòu)設計優(yōu)化[8]。

擴壓器作為離心壓縮機的核心靜止部件,對整個模型級的影響十分顯著。Liu等采用數(shù)值模擬的方法分析了葉輪和擴壓器之間的相互作用對葉輪尾部及擴壓器內(nèi)部氣流的影響[9];Anish等用數(shù)值模擬的方法分析了4種不同擴壓器下離心壓縮機的性能[10];聞蘇平等數(shù)值研究了某離心風機葉輪和擴壓器耦合的非定常流動[11];李凱等采用實驗和數(shù)值研究相結(jié)合的方法對某離心風機的擴壓器進行了研究[12]。

翼型葉片擴壓器由于造價高的緣故,所以多用于航空航天領域[13],具有結(jié)構(gòu)緊湊、尺寸小和高轉(zhuǎn)速的特點,而且運行在高馬赫數(shù)下。一般離心壓縮機模型級研究采用的結(jié)構(gòu)多是葉輪、擴壓器和蝸殼[14],但是同時帶有彎道和回流器的整個模型級的實驗與數(shù)值模擬的對比研究并不多見。

本文對某研發(fā)的工業(yè)用帶翼型葉片擴壓器離心壓縮機模型級的內(nèi)部流動進行了數(shù)值模擬,同時開展了模型級外特性的實驗研究,分析了該模型級內(nèi)部流動特點、損失機理及其對級性能的影響。

1 模型級及其計算模型

本文離心壓縮機模型級葉輪和葉片擴壓器的基本參數(shù)如表1和表2所示,其中葉輪轉(zhuǎn)速為11 515 r/min,進口總溫為30 ℃,進口總壓為101 325 Pa,設計工況下質(zhì)量流量為4.05 kg/s,介質(zhì)為空氣。葉片擴壓器的葉片形式為翼型。

圖1為模型級的計算網(wǎng)格。根據(jù)葉輪葉片的特點,對壁面附近和葉片進、出口處的網(wǎng)格進行了加密處理。進口邊界條件為給定總溫、總壓和速度方向,出口邊界條件為給定質(zhì)量流量。旋轉(zhuǎn)葉輪和葉片擴壓器的耦合計算采用滑移運動網(wǎng)格技術(shù)。在模型級數(shù)值模擬試算階段采用了3套計算網(wǎng)格,其中單通道葉輪內(nèi)部的網(wǎng)格數(shù)分別為206 000、309 500和460 000,模型級其他部件內(nèi)部的計算網(wǎng)格數(shù)按照相應比例分別進行調(diào)整。根據(jù)網(wǎng)格無關性原則和計算效率原則,最后確定采用第二套網(wǎng)格用于詳細的計算分析(單通道葉輪的計算網(wǎng)格數(shù)為309 500)。計算中取壁面網(wǎng)格y+值為30,壁面邊界層效應計算采用壁面函數(shù)法。

表1 模型級葉輪的幾何參數(shù)

表2 模型級葉片擴壓器的幾何參數(shù)

(a)壓縮機級網(wǎng)格

(b)單通道網(wǎng)格圖1 模型級計算網(wǎng)格

2 翼型葉片擴壓器模型級外特性

模型級的外特性是由研發(fā)模型級的內(nèi)部流動特性決定的。為了研究模型級的外特性及其影響因素,在對模型級內(nèi)流進行詳細分析之前,首先采用數(shù)值模擬和實驗的方法對模型級的外特性進行了對比研究,以驗證數(shù)值模型和相關算法的準確性和可靠性。

標準模型級性能實驗采用了JB/T3165—1999《離心和軸流式鼓風機和壓縮機熱力性能試驗方法》中的開式進、出氣實驗方法,進口管道裝有端面噴嘴以測量流量,出氣管道裝有閘閥以調(diào)節(jié)流量。圖2為模型級實驗簡圖,測量截面為模型級的進、出口截面。

圖2 模型級實驗簡圖

圖3為通過數(shù)值模擬和實驗分別獲得的模型級無因次性能曲線。由圖3可見,數(shù)值模擬預測的模型級性能曲線與實驗獲得的結(jié)果一致,由數(shù)值模擬獲得的級效率和總能頭系數(shù)較實驗值略大,級效率計算值與實驗值的平均相對偏差為3.5%。

由圖3還可見,在工況范圍內(nèi),整個模型級的性能較好,效率性能曲線非常平坦,且平均效率高達82%,表明所研發(fā)的模型級具有優(yōu)良性能。對比計算值和實驗值表明,本文采用的數(shù)值模型和獲得的模擬結(jié)果是可靠的,能夠比較準確地預測離心壓縮機模型級的內(nèi)部流動和性能。

圖3 模型級無因次性能曲線

3 模型級內(nèi)部流動分析

對3種典型流量工況下模型級內(nèi)部的流動和總壓損失進行了分析,這3種工況分別為:小流量工況(質(zhì)量流量為2.475 kg/s)、設計流量工況(質(zhì)量流量為4.05 kg/s)、大流量工況(質(zhì)量流量為5.40 kg/s)。

3.1 設計工況下模型級內(nèi)部通流速度發(fā)展

圖4為設計工況下葉輪和擴壓器出口截面上的通流速度(Cm/u2)分布。由圖4a可見,在葉輪出口通流速度并不均勻,在輪蓋與葉片吸力面的角區(qū)有明顯的尾跡區(qū)存在,而最大速度分布在輪蓋與壓力面的角區(qū)以及輪盤與吸力面的角區(qū)。由圖4b可見,對于擴壓器,主流區(qū)域的流動速度逐漸均勻化,同時在前、后壁面上形成了明顯的壁面邊界層,且在擴壓器出口截面最為明顯。

(a)葉輪出口截面

(b)葉片擴壓器出口截面圖4 設計流量下葉輪和擴壓器出口截面內(nèi)的通流速度分布

3.2 不同工況下擴壓器子午面內(nèi)速度分布

圖5為不同流量工況下不同跨盤蓋截面上擴壓器流道子午面的速度分布。其中,小流量工況取靠近葉腹的截面,設計工況取中間截面,大流量工況取靠近葉背的截面,靠近葉腹的截面、中間截面和靠近葉背的截面分別指流道內(nèi)在周向距葉腹20%、50%和80%流道寬度的位置。

小流量下,氣流主要沖擊葉片擴壓器葉背,葉腹處易出現(xiàn)旋渦(見圖5a),在靠近葉腹的截面上,整個子午面內(nèi)幾乎都有回流存在,而其他子午面內(nèi)的流動狀況良好。在設計流量下,所有子午面內(nèi)的流動速度分布都比較均勻(見圖5b)。大流量下,氣流主要沖擊葉片擴壓器葉腹,葉背處易形成旋渦區(qū)(見圖5c),在靠近葉背子午面內(nèi)的大部分區(qū)域都有回流形成。

(a)小流量 (b)設計流量 (c)大流量圖5 擴壓器子午面內(nèi)速度分布

3.3 不同工況下跨葉片截面內(nèi)流動跡線分布

跡線表示氣體運動的軌跡,它是將氣體的流動速度在整個區(qū)域內(nèi)隨時間進行積分后得到的氣體運動的軌跡線。了解了氣體的運動軌跡,就可以知道氣流在擴壓器中的流動情況。圖6～圖8分別顯示了不同工況下不同跨葉片截面內(nèi)(5%葉高、50%葉高、95%葉高)氣流跡線分布。

小流量下(見圖6),葉片擴壓器葉腹附近旋渦的大小和位置在不同葉高處并不相同。在靠近輪盤處,旋渦區(qū)范圍最大,幾乎影響了整個流道;在中間葉高截面內(nèi),旋渦局限在葉腹尾緣;在靠近輪蓋處,旋渦緊貼整個葉腹。

(a)5%葉高 (b)50%葉高 (c)95%葉高圖6 小流量下氣流跡線分布

設計流量下(見圖7),氣流沿著流道中心線流動,未產(chǎn)生明顯的旋渦,說明此時擴壓器中的流動損失較小,運行效率較高。

(a)5%葉高 (b)50%葉高 (c)95%葉高圖7 設計流量下氣流跡線分布

(a)5%葉高 (b)50%葉高 (c)95%葉高圖8 大流量下氣流跡線分布

大流量下(見圖8),葉背附近的旋渦由輪盤到輪蓋截面逐漸向出口移動。在靠近輪盤處,旋渦幾乎覆蓋了整個葉背;中間葉高處,旋渦向擴壓器出口方向移動;在靠近輪蓋處,旋渦位于擴壓器出口位置。

4 擴壓器內(nèi)總壓和馬赫數(shù)分布

4.1 不同工況下擴壓器子午面內(nèi)總壓分布

擴壓器內(nèi)的總壓分布反映了擴壓器內(nèi)部的流動損失。圖9～圖11為不同流量工況下模型級葉片擴壓器內(nèi)不同跨盤蓋子午面上(靠近葉背、中間截面、靠近葉腹)的總壓(P0/P0,in)分布。

小流量下(見圖9),靠近葉背處的總壓損失較大,中間子午面和靠近葉腹子午面總壓損失較小,而在擴壓器進口處總壓均是輪蓋側(cè)高于輪盤側(cè),總壓損失也是輪蓋側(cè)高于輪盤側(cè)。

(a)靠近葉背 (b)中間子午面 (c)靠近葉腹圖9 小流量下擴壓器子午面內(nèi)總壓分布

設計流量下(見圖10),沿子午線方向的總壓梯度相對較小,因此總壓損失也就最小。

(a)靠近葉背 (b)中間子午面 (c)靠近葉腹圖10 設計流量下擴壓器子午面內(nèi)總壓分布

大流量下(見圖11),靠近葉背的總壓在擴壓器進口處開始急劇減小,輪盤處總壓減小較輪蓋處更快,而在靠近擴壓器出口區(qū)域總壓趨于均勻;在中間子午面內(nèi),總壓從進口到出口的梯度比較大;在靠近葉腹處,氣流不斷地向葉腹靠攏,總壓不但沒有損失,反而有局部上升的趨勢,輪盤處最明顯。

(a)靠近葉背 (b)中間子午面 (c)靠近葉腹圖11 大流量下擴壓器子午面內(nèi)總壓分布

對比圖9～圖11可見,大流量下在靠近葉背處和中間子午面內(nèi)的總壓損失要遠大于小流量及設計流量下的總壓損失。比較而言,小流量下的總壓損失并不十分明顯,這表明大流量下在葉背附近形成的旋渦對總壓損失的影響較大,而小流量下在葉腹處形成的旋渦對總壓損失的影響并不顯著。

4.2 不同工況下擴壓器子午面內(nèi)馬赫數(shù)分布

擴壓器中氣流的絕對馬赫數(shù)是衡量氣流在其中將動能轉(zhuǎn)化為壓力能程度(即擴壓)的一個指標,表征氣流是否得到有效滯止。圖12～圖14為不同工況下擴壓器內(nèi)不同跨盤蓋子午面(靠近葉背、中間子午面、靠近葉腹)上的絕對馬赫數(shù)分布。

小流量下(見圖12),在靠近葉背處氣流的滯止效果最好。中間子午面內(nèi)靠近輪盤側(cè)沿子午線方向的馬赫數(shù)梯度較小,氣流的滯止效果并不明顯,這是由于靠近輪盤處葉腹附近的旋渦幾乎影響到整個流道的緣故。在葉腹附近氣流的絕對馬赫數(shù)普遍高于其他子午流面,說明受葉腹附近旋渦的影響,葉腹的流動沒有得到有效滯止。

(a)靠近葉背 (b)中間子午面 (c)靠近葉腹圖12 小流量下擴壓器子午面內(nèi)絕對馬赫數(shù)分布

(a)靠近葉背 (b)中間子午面 (c)靠近葉腹圖13 設計流量下擴壓器子午面內(nèi)絕對馬赫數(shù)分布

(a)靠近葉背 (b)中間子午面 (c)靠近葉腹圖14 大流量下擴壓器子午面內(nèi)絕對馬赫數(shù)分布

設計流量下(見圖13),氣流的絕對馬赫數(shù)沿軸向和子午線方向分布都比較均勻且在擴壓器前半段的梯度較大,說明在擴壓器中的各個方向上氣流都得到了有效的滯止,擴壓情況良好。

大流量下(見圖14),氣流不但沒有得到滯止,還因為流道的阻塞而加速,葉腹附近尤為嚴重?？梢钥吹?在這種情況下,葉片擴壓器整體擴壓已經(jīng)不再有效。

5 結(jié) 論

本文對研發(fā)的帶翼型葉片擴壓器三元葉輪離心壓縮機模型級的內(nèi)部流動和外特性進行了數(shù)值模擬和實驗研究,得到以下結(jié)論。

(1)所研發(fā)的離心壓縮機模型級具有較寬的穩(wěn)定運行范圍,氣動性能曲線的數(shù)值結(jié)果與實驗值吻合較好。

(2)就模型級內(nèi)部流動特性而言,小流量下葉片擴壓器葉腹出現(xiàn)了旋渦區(qū)且在5%葉高處的旋渦范圍最大;設計流量下葉輪出口輪蓋與葉片吸力面的角區(qū)依然存在明顯的尾跡區(qū);大流量下葉片擴壓器葉背形成了旋渦區(qū)。

(3)葉片擴壓器內(nèi)的流動損失和擴壓效果由不同工況下的流動特性決定。小流量下氣流在葉片擴壓器葉腹沒有得到有效滯止;大流量下(5.4 kg/s),氣流不但沒有得到滯止,還因為流道的阻塞效應使氣流加速,葉腹附近尤為嚴重,同時靠近葉背處和中間子午面內(nèi)的總壓損失遠大于小流量及設計流量下的總壓損失。

[1] KATANIS T. Use of arbitrary quasi-orthogonals for calculating flow distribution in meridional plane of a turbomachine, NASA-TN-D 2546 [R]. Washington DC, USA: NASA, 1964.

[2] HATHAWAY M D, CHRISS R M, STRAZISAR A J, et al. Laser anemometer measurement of the three-dimensional rotor flow field in the NASA low-speed centrifugal compressor, NASA 3527 [R]. Washington DC, USA: NASA, 1995.

[3] 饒江, 李京平. 一種離心壓氣機子午流道的優(yōu)化成形方法 [J]. 航空動力學報, 1992, 7(4): 325-328. RAO Jiang, LI Jingping. A method for optimizing the meridional passage in a centrifugal compressor [J]. Journal of Aerospace Power, 1992, 7(4): 325-328.

[4] 樊會元, 席光, 王尚錦. 離心壓縮機二元和三元葉輪內(nèi)流分析 [J]. 風機技術(shù), 1999(1): 5-8. FAN Huiyuan, XI Guang, WANG Shangjin. Internal flow analysis between 2-D & 3-D impellers of centrifugal compressor [J]. Compressor Blower & Fan Technology, 1999(1): 5-8.

[5] 聞蘇平, 胡小文, 茹重智, 等. 馬赫數(shù)對離心壓縮機高壓端小流量系數(shù)基本級性能影響的研究 [J]. 化工設備與管道, 2011, 48(z2): 16-21. WEN Suping, HU Xiaowen, RU Zhongzhi, et al. Influence of Ma on model stage performance for small flow centrifugal compressor [J]. Process Equipment & Piping, 2011, 48(z2): 16-21.

[6] 張楚華, 王寶潼, 欒輝寶, 等. 葉頂間隙對離心葉輪氣動性能影響研究 [J]. 流體機械, 2006, 34(12): 13-16. ZHANG Chuhua, WANG Baotong, LUAN Huibao, et al. Study on influence of tip clearance on aerodynamic performance of centrifugal impeller [J]. Fluid Machinery, 2006, 34(12): 13-16.

[7] 劉瑞韜, 徐忠. 離心壓氣機實際氣體三維粘性流場分析 [J]. 西安交通大學學報, 2004, 8(3): 317-321. LIU Ruitao, XU Zhong. Numerical investigation to the flow of real gas in a centrifugal compressor [J]. Journal of Xi’an Jiaotong University, 2004, 8(3): 317-321.

[8] 么立新. 小流量模型級內(nèi)部流動數(shù)值實驗與實驗研究 [D]. 大連: 大連理工大學, 2006.

[9] LIU Y W, LIU B J. Numerical investigation of impeller-diffuser interaction in a centrifugal compressor [J]. Journal of Aerospace Power, 2009, 24(12): 695-702.

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

[11]聞蘇平, 張楚華, 李景銀. 旋轉(zhuǎn)葉輪和葉片擴壓器耦合的非定常流動計算 [J]. 西安交通大學學報, 2004, 38(7): 754-757. WEN Suping, ZHANG Chuhua, LI Jingyin. Calculation of unsteady flow for interaction between rotating impeller and vaned diffuser [J]. Journal of Xi’an Jiaotong University, 2004, 38(7): 754-757.

[12]李凱, 曹淑珍, 林梅, 等. 直板型葉片擴壓器流場的實驗測量與數(shù)值研究 [J]. 西安交通大學學報, 2006, 40(3): 279-283. LI Kai, CAO Shuzhen, LIN Mei, et al. Experimental and numerical study on flows in straight plate diffuser of centrifugal fan [J]. Journal of Xi’an Jiaotong University, 2006, 40(3): 279-283.

[13]STEGLICH T, KITZINGER J, SEUME J R, et al. Improved diffuser/volute combinations for centrifugal compressors [J]. ASME Journal of Turbomachinery, 2008, 130(1): 011014.

[14]TREBINJAC I, KULISA P, BULOT N, et al. Effect of unsteadiness on the performance of a transonic centrifugal compressor stage [J]. ASME Journal of Turbomachinery, 2009, 131(4): 041011.

(編輯 苗凌)

(1.西安交通大學能源與動力工程學院,710049,西安;2.中國石油寧夏石化分公司機動設備處,750026,銀川)

通過對在研工業(yè)用帶翼型葉片擴壓器離心壓縮機模型級的內(nèi)部流動進行數(shù)值模擬,以及對模型級外特性進行了實驗研究,分析了模型級內(nèi)部流動特點、損失機理及其對級性能的影響。研究發(fā)現(xiàn):該離心壓縮機模型級性能較好,工況范圍內(nèi)效率曲線非常平坦,穩(wěn)定運行范圍較寬,整級效率高達82%;設計流量下葉輪出口輪蓋與葉片吸力面的角區(qū)有明顯的尾跡區(qū)存在,小流量下進入葉片擴壓器內(nèi)的氣流在葉腹的流動沒有得到有效滯止,大流量下氣流以沖擊葉片擴壓器葉腹為主,同時在葉背上出現(xiàn)了較大的旋渦區(qū);不同工況下葉片擴壓器內(nèi)的流動損失和擴壓效果由不同工況下的流動特性決定。數(shù)值模擬和實驗研究獲得的模型級外特性曲線一致,表明所采用的數(shù)值模型能夠準確預測模型級內(nèi)部的流動特性。

離心壓縮機;葉片擴壓器;模型級;數(shù)值模擬;實驗研究

Internal Flow at Developed Model Stage in Centrifugal Compressor with Vaned Diffuser

ZHOU Jun’an1,LIU Lijun1,XIAO Ping2,CHEN Yanfeng2,GUO Jiangfeng1,CHEN Xuejiang1

(1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China; 2. Power-Mechanical Department, Ningxia Petrochemical Company of Petro China, Yinchuan 750026, China)

The internal flow at model stage in centrifugal compressor with 3D centrifugal impeller and vaned diffuser is investigated by numerical simulation, the external performance at model stage is determined by experiments. The internal flow at model stage and the total pressure loss, as well as their influences on the stage performance, are analyzed. The results show that the centrifugal compressor at the model stage well performs with a very flat efficiency curve and the average efficiency gets as high as 82%. However, a wake region can be obviously observed at the outlet of the impeller near the shroud and the suction surface of the impeller when it runs under design flow rate condition. When the flow rate is greatly lower than the design level, the flow region near the pressure surface of the diffuser vane cannot keep effectively stagnant, leading to an extra loss of energy. When the flow rate is greatly higher than the design level, the airflow notably impinges on the pressure surface of the diffuser vane, and a large vortex region generates near the suction surface of the diffuser vane. The vaned diffuser performance depends upon the internal flow conditions during different operations. The predicted performance curve at this stage coincides well with that obtained by experiments.

centrifugal compressor; vaned diffuser; model stage; numerical simulation; experimental investigation

2014-12-22。 作者簡介:周俊安(1990—),男,碩士生;劉立軍(通信作者),男,教授。 基金項目:國家重點基礎研究發(fā)展計劃資助項目(2012CB026001);中央高校基本科研業(yè)務費專項資金資助項目(2010jdgz08)。

時間:2015-06-17

http:∥www.cnki.net/kcms/detail/61.1069.T.20150617.0902.006.html

10.7652/xjtuxb201509006

TH452

A

0253-987X(2015)09-0030-06