離心壓氣機(jī)相似轉(zhuǎn)速數(shù)值研究

田澤，鐘易成，徐偉祖，沈永涵

(1. 南京航空航天大學(xué) 動(dòng)力與能源學(xué)院，南京 210016； 2. 南京普國科技有限公司，南京 210016)

?

離心壓氣機(jī)相似轉(zhuǎn)速數(shù)值研究

田澤1，鐘易成1，徐偉祖2，沈永涵1

(1. 南京航空航天大學(xué) 動(dòng)力與能源學(xué)院，南京 210016； 2. 南京普國科技有限公司，南京 210016)

摘要：以某型地面燃機(jī)的離心式低壓壓氣機(jī)為對(duì)象進(jìn)行了全三維數(shù)值模擬，分別計(jì)算了折合轉(zhuǎn)速相同但進(jìn)口總溫、總壓及物理轉(zhuǎn)速不同的7種狀態(tài)的壓氣機(jī)工作特性曲線，其中物理轉(zhuǎn)速分別為90%、100%及105%。結(jié)果表明：折合轉(zhuǎn)速相同時(shí)流場(chǎng)的相似程度較高但性能差別較大，隨著進(jìn)口總溫升高、總壓減小雷諾數(shù)降低，壓氣機(jī)增壓比、效率降低。

關(guān)鍵詞：離心式壓氣機(jī)；雷諾數(shù)；壓氣機(jī)性能

0引言

在燃?xì)廨啓C(jī)及其部件的性能研究中相似理論具有重要的地位[1-2]，應(yīng)用相似理論可以把燃?xì)廨啓C(jī)各部件在特定條件下的工作特性推廣到與它相似的情況下。相似理論認(rèn)為當(dāng)雷諾數(shù)進(jìn)入自?；瘏^(qū)即不再是相似準(zhǔn)則。目前國內(nèi)外學(xué)者的研究主要集中于高空低雷諾數(shù)問題[3-5]，對(duì)于地面燃?xì)廨啓C(jī)，在雷諾數(shù)較高的情況下，地面大氣條件、海拔高度或進(jìn)排氣系統(tǒng)損失不同導(dǎo)致發(fā)動(dòng)機(jī)進(jìn)口參數(shù)改變，引起的雷諾數(shù)變化同樣會(huì)對(duì)發(fā)動(dòng)機(jī)各部件的性能產(chǎn)生影響。尤其是低壓壓氣機(jī)進(jìn)口參數(shù)值較小，各參數(shù)通過比值進(jìn)行換算時(shí)偏差較為明顯。本文以某型地面燃機(jī)的低壓離心壓氣機(jī)為研究對(duì)象，采用數(shù)值模擬的方法分別計(jì)算了物理轉(zhuǎn)速為90%、100%及105%，但折合轉(zhuǎn)速均為100%的不同進(jìn)口條件下共7個(gè)狀態(tài)的工作特性，分析了相似轉(zhuǎn)速下離心壓氣機(jī)性能的變化及雷諾數(shù)對(duì)壓氣機(jī)性能的影響。

1物理模型

本文采用商業(yè)NUMECA軟件對(duì)某型地面燃機(jī)低壓級(jí)離心式壓氣機(jī)進(jìn)行數(shù)值模擬，研究了相似轉(zhuǎn)速對(duì)壓氣機(jī)性能的影響。

低壓級(jí)壓氣機(jī)二維結(jié)構(gòu)如圖1所示，幾何模型含6個(gè)部件，分別為進(jìn)口支板、預(yù)旋葉片、工作葉片、徑向擴(kuò)壓器、軸向擴(kuò)壓器和徑向回流器。其中，進(jìn)口支板4個(gè)；預(yù)旋葉片19個(gè)；轉(zhuǎn)子葉片為帶冠葉片，總數(shù)為22；徑向和軸向擴(kuò)壓器，葉片數(shù)分別為17和120；徑向回流器葉片數(shù)為9。

圖1　低壓級(jí)壓氣機(jī)二維結(jié)構(gòu)圖

2數(shù)值計(jì)算方法

高品質(zhì)的網(wǎng)格有利于提高數(shù)值計(jì)算的精度、速度和收斂性。采用IGG/AUTOGRID生成網(wǎng)格,低壓級(jí)壓氣機(jī)網(wǎng)格如圖2所示，網(wǎng)格總量為160萬，其中，轉(zhuǎn)子葉片網(wǎng)格為25萬，如圖3所示。近壁面第一層網(wǎng)格高度為0.006mm，網(wǎng)格最小正交性16.032，最大網(wǎng)格延展比3.074，根據(jù)經(jīng)驗(yàn)滿足湍流模型計(jì)算要求。

圖2　低壓級(jí)壓氣機(jī)通道網(wǎng)格

圖3　轉(zhuǎn)子葉片計(jì)算網(wǎng)格

采用FINE/TURBO進(jìn)行數(shù)值計(jì)算，求解三維雷諾平均N-S方程，數(shù)值方法采用時(shí)間追趕的有限體積法，空間離散采用二階精度中心差分格式，時(shí)間離散應(yīng)用二階后差歐拉格式，湍流模型采用S-A (spalart-allmaras)單方程模型，流動(dòng)介質(zhì)為理想氣體，使用多重網(wǎng)格法、局部時(shí)間步長和殘差光順等多種措施加速收斂。

根據(jù)相似理論，壓氣機(jī)特性曲線可以在任意進(jìn)口溫度下測(cè)取，通過相似理論進(jìn)行換算應(yīng)該得到同樣的結(jié)果。折合參數(shù)[2]就是應(yīng)用相似理論，把不同大氣條件下所得參數(shù)換算成海平面國際標(biāo)準(zhǔn)大氣條件下相似狀態(tài)的對(duì)應(yīng)參數(shù)。

折合轉(zhuǎn)速：

折合流量：

表1為不同狀態(tài)下低壓級(jí)壓氣機(jī)通道計(jì)算域的進(jìn)口邊界參數(shù)，其中編號(hào)NO.4_100per是在設(shè)計(jì)轉(zhuǎn)速下，模擬壓氣機(jī)標(biāo)準(zhǔn)大氣狀態(tài)下的工作特性；調(diào)整進(jìn)口總溫使物理轉(zhuǎn)速為90%、105%的折合轉(zhuǎn)速均為100%，根據(jù)相似理論總壓對(duì)結(jié)果沒有影響，故相同物理轉(zhuǎn)速下分別給定了不同進(jìn)口總壓。

表1　不同狀態(tài)下壓氣機(jī)參數(shù)

3結(jié)果與分析

3.1　氣動(dòng)性能

運(yùn)用NUMECA軟件包對(duì)上述7種狀態(tài)的低壓壓氣機(jī)進(jìn)行了數(shù)值模擬，得到了壓氣機(jī)壓比、效率隨折合流量變化的工作特性。

圖4及圖5給出了不同狀態(tài)下壓氣機(jī)的流量-壓比特性及流量-效率特性曲線。由圖可以看出，在同一折合轉(zhuǎn)速下，不同狀態(tài)的壓氣機(jī)工作特性具有明顯差異，說明本文的離心式壓氣機(jī)特性受進(jìn)口總溫、總壓影響比較顯著。隨著壓比的減小，壓氣機(jī)通流能力減弱，折合流量減小，當(dāng)壓比減小到約3.5以下時(shí)，流量發(fā)生堵塞，壓比繼續(xù)減小流量基本不再變化，且不同狀態(tài)下堵塞流量均不同；隨著折合流量的增大壓氣機(jī)效率減小，當(dāng)流量堵塞時(shí)效率急劇下降。

圖4　壓氣機(jī)流量壓比特性

圖5　壓氣機(jī)流量效率特性

為更直觀的研究壓氣機(jī)內(nèi)部性能的變化，圖6給出了折合流量相同時(shí)不同狀態(tài)的壓比、效率隨雷諾數(shù)的變化。根據(jù)相似理論，壓氣機(jī)在某一工作狀態(tài)與標(biāo)準(zhǔn)大氣條件下工作的折合轉(zhuǎn)速與折合流量相等時(shí)，這兩個(gè)狀態(tài)就相似，壓氣機(jī)其他各性能相似參數(shù)也應(yīng)該相等。雷諾數(shù)的定義[6]基于離心壓氣機(jī)轉(zhuǎn)子葉片出口周向速度及轉(zhuǎn)子葉片出口高度。

由圖6可以看出各狀態(tài)雷諾數(shù)均比較高，從NO.1到NO.7雷諾數(shù)5×105增加到7.6×105，增加超過50%，說明不同狀態(tài)參數(shù)對(duì)雷諾數(shù)影響比較嚴(yán)重。雷諾數(shù)與壓比、效率呈現(xiàn)出近似線性關(guān)系，雷諾數(shù)越低壓氣機(jī)壓比越小，損失越大，效率越低。當(dāng)進(jìn)口總壓保持不變(NO.1、4、7)時(shí)，隨著物理轉(zhuǎn)速的增加及進(jìn)口總溫的升高雷諾數(shù)減小，壓氣機(jī)的增壓比、效率均降低。當(dāng)物理轉(zhuǎn)速及進(jìn)口總溫保持不變(NO.1、2、3及NO.5、6、7)時(shí)，隨著進(jìn)口總壓的減小雷諾數(shù)降低，壓氣機(jī)的增壓比、效率均降低。

圖6　壓比、效率隨雷諾數(shù)變化

為更好并且更方便地對(duì)比不同條件下壓氣機(jī)內(nèi)部流動(dòng)變化，本文選取編號(hào)分別為NO.2_90per、NO.4_100per、NO.6_105per的3個(gè)狀態(tài)(因?yàn)榇?個(gè)狀態(tài)的折合流量與物理流量相同)，比較流量相同的計(jì)算點(diǎn)，分析其他參數(shù)的變化。具體參數(shù)如表2所示，可以看出隨著物理轉(zhuǎn)速由90%增加到105%，以標(biāo)準(zhǔn)狀態(tài)為基準(zhǔn)，雷諾數(shù)由+19.87%減小到-8.35%，，壓氣機(jī)增壓比由+4.92%減小至-2.66%，效率由+3.76%降至-2.03%。雷諾數(shù)共降低28.225 %，增壓比及效率分別降低了7.581%、5.787%。

表2　不同狀態(tài)等折合流量壓氣機(jī)性能參數(shù)

表3給出了3種狀態(tài)壓氣機(jī)各部件通道內(nèi)出口總壓與進(jìn)口總壓的比值?？梢钥闯鲭S著雷諾數(shù)的減小，進(jìn)口支板、預(yù)旋葉片、徑向擴(kuò)壓器、軸向擴(kuò)壓器和徑向回流器的總壓恢復(fù)系數(shù)都減小。其中進(jìn)口支板及預(yù)旋葉片的變化都比較小，隨雷諾數(shù)的減小轉(zhuǎn)子葉片產(chǎn)生的增壓比由2.688%減小到-1.784%，而由于轉(zhuǎn)子葉片出口氣流速度較高，位于轉(zhuǎn)子葉片后方的徑向擴(kuò)壓器和軸向擴(kuò)壓器總壓損失也遠(yuǎn)高于其他非旋轉(zhuǎn)部件。

表3　不同狀態(tài)等折合流量各部件總壓比

3.2　流場(chǎng)相似性

為了更詳細(xì)地分析雷諾數(shù)對(duì)壓氣機(jī)性能的影響，研究其對(duì)壓氣機(jī)性能的影響機(jī)理，圖7給出了3種狀態(tài)下轉(zhuǎn)子葉片S2流面流場(chǎng)的相對(duì)馬赫數(shù)分布云圖?？梢钥闯?種流場(chǎng)的馬赫數(shù)分布比較相似，葉片葉根與葉尖都存在一個(gè)較明顯的低速區(qū)，馬赫數(shù)較小。葉根的低速區(qū)主要是葉輪進(jìn)口處的擴(kuò)壓作用造成的，低速區(qū)范圍較小，沒有導(dǎo)致?lián)p失明顯增加，但是由于此處低能流體的存在，不利于控制二次流的發(fā)展，在葉尖可看到明顯的馬赫數(shù)變化，說明此處存在泄露流。

圖7　轉(zhuǎn)子通道相對(duì)Ma分布

對(duì)比不同狀態(tài)的流場(chǎng)分布，可以看出隨著雷諾數(shù)的降低，物理轉(zhuǎn)速增加，馬赫數(shù)分布呈現(xiàn)“兩極分化”的趨勢(shì)，通道內(nèi)的高馬赫數(shù)范圍增加，但同時(shí)低馬赫數(shù)范圍也增加。葉片前緣處的馬赫數(shù)增大，導(dǎo)致分離區(qū)的范圍增大；葉尖低速區(qū)逐漸增大，同時(shí)低速區(qū)的起始位置逐漸向前移動(dòng)，這主要是由于離心壓氣機(jī)轉(zhuǎn)子葉片受離心力影響，氣流從葉根向葉尖遷移流動(dòng)，流體低能區(qū)的由葉根向葉尖的遷移流動(dòng)也逐漸增強(qiáng)，低能區(qū)的堆積導(dǎo)致葉尖低速區(qū)加強(qiáng)，低馬赫數(shù)區(qū)的增大，導(dǎo)致了工作葉片的增壓比減小，這也是離心葉輪總壓損失重要原因。

圖8給出了徑向擴(kuò)壓器50%高度絕對(duì)馬赫數(shù)云圖和流線圖，可以看出3種流場(chǎng)的馬赫數(shù)分布及流線均比較相近。在葉背均出現(xiàn)了明顯的低速區(qū)，氣流發(fā)生了分離。隨著雷諾數(shù)的減小，轉(zhuǎn)子葉片出口Ma減小，同時(shí)氣流角發(fā)生改變，攻角增大，在擴(kuò)壓器葉片葉背分離區(qū)起始位置前移、范圍增大，從而導(dǎo)致了總壓損失增大，通道的通流能力減弱。

圖8　徑向擴(kuò)壓器絕對(duì)馬赫數(shù)及流線圖

4結(jié)論

本文以某型地面燃機(jī)離心式低壓壓氣機(jī)為對(duì)象，采用數(shù)值模擬的方法分析了7種狀態(tài)的部件特性，對(duì)工作特性曲線及壓氣機(jī)流場(chǎng)進(jìn)行了分析，得出的結(jié)論如下：

1) 對(duì)于本文所研究的離心壓氣機(jī)，折合轉(zhuǎn)速與折合流量均相同時(shí)，從NO.2_90per到NO.6_105per狀態(tài)雷諾數(shù)共降低28.225 %，增壓比降低了7.581%，效率降低了5.787%。

2) 折合轉(zhuǎn)速相等時(shí)，隨著進(jìn)口總溫的升高、總壓的降低雷諾數(shù)減小，壓氣機(jī)的增壓比、效率均減小。

3) 本文所研究的離心壓氣機(jī)在折合轉(zhuǎn)速相等時(shí)，通道內(nèi)的流場(chǎng)相似程度比較高，流場(chǎng)存在一定分離現(xiàn)象，但相似換算參數(shù)與標(biāo)準(zhǔn)狀態(tài)存在較大差別。

參考文獻(xiàn):

[1] 鄒滋祥. 相似理論在葉輪機(jī)械模型研究中的應(yīng)用[M]. 北京：科學(xué)出版社, 1984.

[2] 朱行健，王雪瑜. 燃?xì)廨啓C(jī)工作原理及性能[M]. 北京：科學(xué)出版社, 1992.12.

[3] Wassell A B. Reynolds Number effects in Axial Compressors[J ] . ASME Journal of Engineering f or Power, April,1968.

[4] 鄒正平，葉建，劉火星,等. 低雷諾數(shù)低壓渦輪內(nèi)部流動(dòng)及其氣動(dòng)設(shè)計(jì)技術(shù)研究進(jìn)展[J].力學(xué)進(jìn)展，2007,37(4).

[5] 喬渭陽，王占學(xué)，伊進(jìn)寶. 低雷諾數(shù)渦輪流動(dòng)損失控制技術(shù)研究[J].推進(jìn)技術(shù).2005.

[6] CASEY MV. The effects of Reynolds number on the efficiency of centrifugal compressor stages. ASME Journal of Engineering for Gas Turbines and Power . 1985.

Numerical Analysis of Centrifugal Compressor in Similar Speed

TIAN Ze1, ZHONG Yi-cheng1, XU Wei-zu2,SHEN Yong-han1

(1. College. of Energy and Power, Nanjing University. of Aeronautics and Astronautics, Nanjing 210016, China；

2. The Limited Company of Nanjing Puguo Technology, Nanjing 210016, China)

Abstract:A type of low pressure centrifugal compressor of a ground gas turbine is used to do its three dimensional numerical simulation and calculate the characteristic curres of the compressor in the seven different cases of inlet total temperature, total pressure, and different physical speed of 90%, 100% and 105%, with same conversion speed. The results show that in the same conversion speed the flow of the compressor is quite similar, but its performance is different largely, and with the rising of the inlet total temperature and the reducing the inlet total pressure, Reynolds number, the compressor pressure ratio and efficiency are reduced.

Keywords:centrifugal compressor; Reynolds number; Compressor performance

中圖分類號(hào)：V231.3

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1671-5276(2015)02-0059-04

作者簡介：田澤(1989-)，男，河北石家莊人，碩士研究生，研究方向?yàn)楹娇瞻l(fā)動(dòng)機(jī)總體及葉輪機(jī)氣動(dòng)設(shè)計(jì)。