蔣筑宇，范召林，劉 波

軸流壓氣機(jī)特性計(jì)算及激波模型研究

蔣筑宇1，范召林1，劉 波2

（1.中國空氣動力研究與發(fā)展中心，四川綿陽621000；2.西北工業(yè)大學(xué)動力與能源學(xué)院，西安710072）

為快速準(zhǔn)確預(yù)估軸流壓氣機(jī)特性和激波損失，基于軸流壓氣機(jī)S2流面流線曲率法，分別采用正激波模型和改進(jìn)的雙激波模型，對某型2級跨聲速風(fēng)扇特性進(jìn)行數(shù)值模擬計(jì)算，得到了1 00%設(shè)計(jì)轉(zhuǎn)速近設(shè)計(jì)點(diǎn)與99.76%設(shè)計(jì)轉(zhuǎn)速近堵塞點(diǎn)的總體性能和氣動參數(shù)，以及95%、1 00%和1 1 0%設(shè)計(jì)轉(zhuǎn)速的特性曲線。通過將計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)進(jìn)行對比，分析研究了各激波損失模型在激波損失預(yù)估和風(fēng)扇/壓氣機(jī)特性計(jì)算方面的差異。分析結(jié)果表明：在跨聲速風(fēng)扇/壓氣機(jī)近設(shè)計(jì)點(diǎn)激波損失和特性參數(shù)的計(jì)算中，正激波模型損失徑向分布計(jì)算結(jié)果接近試驗(yàn)值，總壓比和總效率計(jì)算值分別較試驗(yàn)值約低1.96%和2.54%，模型能夠滿足工程需要。而在近堵塞點(diǎn)，改進(jìn)的雙激波模型總損失計(jì)算值更接近試驗(yàn)值，總壓比計(jì)算值和試驗(yàn)值很吻合，總效率計(jì)算值比試驗(yàn)值約高7.28%。改進(jìn)雙激波模型的不同轉(zhuǎn)速線效率特性曲線也明顯更接近試驗(yàn)值，模型能夠較準(zhǔn)確地預(yù)測遠(yuǎn)離設(shè)計(jì)點(diǎn)激波損失和特性參數(shù)。

軸流壓氣機(jī)；流線曲率法；激波模型；損失預(yù)估；特性計(jì)算

0 引言

全3維數(shù)值模擬技術(shù)已經(jīng)取得了很大進(jìn)展。利用基于雷諾平均的可壓縮非定常Navier-Stokes方程求解葉輪機(jī)械內(nèi)部流動可以得到詳細(xì)的流場信息，特性預(yù)估可信度高；但受湍流模型、轉(zhuǎn)捩模型和網(wǎng)格數(shù)量等限制，模型建立復(fù)雜，計(jì)算量大，主要用于驗(yàn)證計(jì)算[1-3]。而基于葉輪機(jī)械準(zhǔn)3維流動通用理論的S2流面通流設(shè)計(jì)，運(yùn)算效率更高，時(shí)間成本更低，目前仍然是葉輪機(jī)械設(shè)計(jì)體系的核心。Novak提出的流線曲率法[4]是1種數(shù)學(xué)上相對簡單，使用經(jīng)驗(yàn)較豐富，應(yīng)用較成熟，計(jì)算參數(shù)也較準(zhǔn)確的S2流面求解方法[5]。流線曲率法數(shù)學(xué)模型基于無黏假設(shè)，對壓氣機(jī)氣動損失的預(yù)估主要依賴于損失模型。

對于級負(fù)荷不高、相對馬赫數(shù)較低的風(fēng)扇/壓氣機(jī)來說，20世紀(jì)60年代初Miller等發(fā)展的2維通道正激波模型就能較成功地預(yù)測激波損失[6-7]；而對于現(xiàn)代高負(fù)荷壓氣機(jī)，在實(shí)際葉柵流動中，隨著來流馬赫數(shù)的提高，激波和附面層相互作用加強(qiáng)，激波的形狀和結(jié)構(gòu)會發(fā)生很大變化，正激波模型對激波損失的預(yù)估誤差隨之變大，用流線曲率法預(yù)測非設(shè)計(jì)工況激波損失前，必須知道激波的結(jié)構(gòu)和波前馬赫數(shù)才能比較準(zhǔn)確地建立起激波損失模型[8]。Bloch的研究表明，在跨聲速風(fēng)扇轉(zhuǎn)子葉尖，激波結(jié)構(gòu)主要為葉柵入口斜激波和通道正激波，并會隨工作狀態(tài)的變化而變化[9-10]；Biollo的研究也得到了相似的結(jié)果[11]；Boyer在Bloch等研究的基礎(chǔ)上，考慮到來流攻角對激波結(jié)構(gòu)的影響，提出了雙激波模型[10]。

本文利用流線曲率管流計(jì)算方法[12]，結(jié)合損失和落后角模型，并在雙激波模型基礎(chǔ)上改進(jìn)，開發(fā)了1套完整的軸流壓氣機(jī)損失預(yù)估和特性計(jì)算方法。參考NASA公開發(fā)表的試驗(yàn)數(shù)據(jù)[13-14]，利用程序?qū)δ承?級高負(fù)荷跨聲速風(fēng)扇進(jìn)行了數(shù)值模擬，得到了近設(shè)計(jì)點(diǎn)和近堵塞點(diǎn)的總體性能參數(shù)和展向參數(shù)分布，以及全工況的特性曲線。通過與試驗(yàn)結(jié)果對比，重點(diǎn)研究了正激波和改進(jìn)的雙激波模型在軸流壓氣機(jī)損失預(yù)估和特性計(jì)算中的效果和差異。

1 數(shù)值方法

計(jì)算對流場做出了定常、絕熱、無黏和軸對稱的假設(shè)，氣體黏性對級作功能力和流道面積的影響用流道堵塞系數(shù)來改進(jìn)。主要控制方程為計(jì)算站沿徑向的完全徑向平衡方程和連續(xù)方程。

2 攻角和落后角模型

2.1 攻角模型

參考Lieblein和Johnsen等基于NACA-65-A10系列低速葉型試驗(yàn)數(shù)據(jù)建立的經(jīng)驗(yàn)關(guān)系式計(jì)算參考攻角i*[15]，即最小損失攻角imin。利用Herrig提出的有效范圍攻角公式計(jì)算近堵塞攻角ichoke和近失速攻角i stall[16-17]。

2.2 落后角模型

與參考攻角一樣，參考落后角也選為最小損失落后角，參考Lieblein提出的經(jīng)驗(yàn)?zāi)Ｐ蚚15]計(jì)算參考落后角δ*。以Hearsey和Boyer模型分量組合方法為基礎(chǔ)展開[10]，非設(shè)計(jì)點(diǎn)落后角

式中：δva為軸向速度比對落后角的影響；δ3D為壓氣機(jī)中流線在不同展向位置對落后角的影響；δm為馬赫數(shù)變化對落后角的影響，在跨聲速轉(zhuǎn)子流動后，馬赫數(shù)對落后角影響很小，可以忽略；δi為非設(shè)計(jì)狀態(tài)下攻角變化對落后角的影響。

3 總壓損失模型

本文考慮的跨聲速風(fēng)扇/壓氣機(jī)葉片排通道總壓損失主要包括由黏性引起的葉片附面層葉型損失和氣流穿越無黏激波所引起的激波損失。

3.1 葉型損失模型

考慮二次流、端壁損失以及激波與附面層干擾所引起的葉頂及葉根附加損失，采用Hearsey和Boyer的方法計(jì)算附加損失]。

3.2 激波損失模型

3.2.1 正激波模型

正激波模型將葉柵通道激波簡化為1道在通道入口處與通道中線垂直的正激波。當(dāng)來流相對馬赫數(shù)Ma'1＞1時(shí)，激波前有效馬赫數(shù)波Ma'e為

式中：Ma'B為葉型吸力面與激波交界點(diǎn)處的氣流馬赫數(shù)。當(dāng) Ma'1＜1 時(shí)，則

若 Ma'e＜1，則激波損失為 0；若 Ma'e＞1，由正激波關(guān)系式計(jì)算激波損失。

3.2.2 改進(jìn)的雙激波模型

在雙激波模型中，激波結(jié)構(gòu)隨攻角變換有3個(gè)特征形態(tài)點(diǎn)，即近堵塞點(diǎn)（NC），峰值效率點(diǎn)（PE）和近失速點(diǎn)（NS），如圖 1所示。

圖1 雙激波模型

當(dāng)來流攻角i等于最小損失攻角imin時(shí)，通道入口斜激波波后馬赫數(shù)即為通道正激波波前馬赫數(shù)。由斜激波關(guān)系式計(jì)算氣流通過斜激波損失。Boyer提出可將激波聲速點(diǎn)激波角εsonic作為激波平均傾斜角，這也是目前的普遍處理方法[7,18]，εsonic由Bloch-Moeckel脫體激波模型方法求解[9]

εsonic僅為馬赫數(shù)的函數(shù)，采用εsonic作為激波角可以排除葉型前緣半徑的影響，但這樣會使斜激波波后馬赫數(shù)降為1。對此，本文在正激波損失計(jì)算中引入波前修正馬赫數(shù)Ma2e

若 Ma2e＞1，則計(jì)算正激波損失；若 Ma2e＜1，則不考慮正激波損失。

隨著來流攻角i減小，通道激波向后移動，當(dāng)攻角等于近堵塞攻角ichoke時(shí)，正激波后移至通道出口，將斜激波和正激波之間的流動假設(shè)為拉瓦爾噴管內(nèi)的流動，由K?nig提出的方法求出正激波前馬赫數(shù)[19]。在接近堵塞狀態(tài)時(shí)激波損失激增，對此，在i介于ichoke和imin之間時(shí)，激波損失利用函數(shù) F（（i-imin）4）插值求得，如圖2所示。

圖2 攻角小于i min時(shí)激波損失函數(shù)

F（（i-imin）4）沿 i的變化趨勢接近實(shí)際激波損失，但與試驗(yàn)結(jié)果的對比表明，在接近堵塞工況時(shí)損失計(jì)算值偏小，致使效率計(jì)算值偏大。對此，本文引入?yún)?shù)l（l＞1）將，損失函數(shù)改進(jìn)為

改進(jìn)函數(shù)仍然為i-imin的4次多項(xiàng)式函數(shù)。如圖2所示，改進(jìn)后，在接近堵塞狀態(tài)時(shí)，激波損失更大；而在攻角較大時(shí)，損失基本與原函數(shù)值相同。

隨著來流攻角i增大，通道激波向前移動，通道入口斜激波增強(qiáng)，將斜激波看作無限層排列的脫體弓型激波[10]。忽略來流馬赫數(shù)通過每道斜激波的變化，氣流經(jīng)過所有激波后壓比為[20]

按斜激波理論計(jì)算脫體激波后馬赫數(shù)，激波斜角ε 由函數(shù) G（（i-imin）4）求得。

當(dāng)攻角i＞inorm時(shí)，葉柵通道激波變?yōu)槿肟诿擉w正激波，激波損失計(jì)算方法與正激波模型一致。inorm由經(jīng)驗(yàn)公式求得

4 計(jì)算結(jié)果及分析

本文研究對象為某2級跨聲速風(fēng)扇/壓氣機(jī)[13-14]，在設(shè)計(jì)狀態(tài)下的總體性能參數(shù)為：總壓比2.8，絕熱效率83.9%，質(zhì)量流量83.5 kg/s。設(shè)計(jì)轉(zhuǎn)速為10720 r/min；進(jìn)口條件為標(biāo)準(zhǔn)大氣條件，軸向進(jìn)氣；靜子出口軸向出氣；第1級轉(zhuǎn)子有減振凸臺，葉型為多圓弧葉型。

4.1 設(shè)計(jì)點(diǎn)程序驗(yàn)證計(jì)算

采用文獻(xiàn)[13]給出的風(fēng)扇設(shè)計(jì)點(diǎn)各級損失和落后角，利用本文程序?qū)υO(shè)計(jì)點(diǎn)進(jìn)行數(shù)值模擬，將計(jì)算得到的總體性能參數(shù)與設(shè)計(jì)值進(jìn)行對比，以驗(yàn)證程序的有效性。設(shè)計(jì)點(diǎn)各級累積增壓比和絕熱效率的設(shè)計(jì)值和數(shù)值模擬結(jié)果對比見表1、2。

表1 設(shè)計(jì)點(diǎn)累積增壓比設(shè)計(jì)值與計(jì)算值

表2 設(shè)計(jì)點(diǎn)累積絕熱效率設(shè)計(jì)值與計(jì)算值 %

從表中可見，設(shè)計(jì)點(diǎn)總體氣動性能計(jì)算值跟設(shè)計(jì)值很接近?？傇鰤罕扔?jì)算值跟設(shè)計(jì)值吻合很好，總絕熱效率計(jì)算值比設(shè)計(jì)值略高，本文程序計(jì)算精度較高。

4.2 近設(shè)計(jì)點(diǎn)特性計(jì)算

分別結(jié)合2種激波模型和落后角模型，利用程序?qū)︼L(fēng)扇近設(shè)計(jì)點(diǎn)進(jìn)行數(shù)值模擬計(jì)算，近設(shè)計(jì)點(diǎn)轉(zhuǎn)速為100%設(shè)計(jì)轉(zhuǎn)速，流量為83.574 kg/s。近設(shè)計(jì)點(diǎn)各級累積增壓比和絕熱效率的設(shè)計(jì)值和數(shù)值模擬結(jié)果對比見表 3、4。

表3 近設(shè)計(jì)點(diǎn)累積增壓比試驗(yàn)值與計(jì)算值

表4 近設(shè)計(jì)點(diǎn)絕熱效率試驗(yàn)值與計(jì)算值 %

從表中數(shù)據(jù)可見，采用正激波模型計(jì)算得到的增壓比和絕熱效率均小于試驗(yàn)值，而雙激波模型計(jì)算結(jié)果均大于試驗(yàn)值。正激波模型總壓比的計(jì)算結(jié)果較試驗(yàn)值約低1.96%，總效率的計(jì)算結(jié)果較試驗(yàn)值約低2.54%。而雙激波模型總壓比的計(jì)算結(jié)果較試驗(yàn)值約高3.42%，總效率的計(jì)算結(jié)果較試驗(yàn)值約高0.53%。

各級轉(zhuǎn)、靜子損失試驗(yàn)值和計(jì)算值沿展向分布如圖3、4所示。

圖3 第1級損失試驗(yàn)值與計(jì)算值

圖4 第2級損失試驗(yàn)值與設(shè)計(jì)值

從圖3中可見，第1級轉(zhuǎn)子損失試驗(yàn)值沿展向總體呈增長趨勢，受凸臺的影響，轉(zhuǎn)子葉中位置損失出現(xiàn)突增。第1級靜子葉中部分損失試驗(yàn)值較小，葉尖和葉根部分較大。轉(zhuǎn)子損失計(jì)算值大致反映出試驗(yàn)值變化趨勢，靜子損失計(jì)算值與試驗(yàn)值也比較接近。從圖4中可見，第2級轉(zhuǎn)子和靜子損失試驗(yàn)值均表現(xiàn)為葉展中部相對較小，葉尖和葉根部分相對較大。轉(zhuǎn)子葉根部分損失計(jì)算值較試驗(yàn)值偏低，而葉中到葉尖部分計(jì)算值與試驗(yàn)值較接近。靜子損失計(jì)算值和試驗(yàn)值吻合很好。在2級轉(zhuǎn)子全展向范圍內(nèi)，正激波模型損失計(jì)算結(jié)果基本高于雙激波模型的，且沿展向差值增大。

各級絕熱效率試驗(yàn)值和計(jì)算值沿展向分布如圖5、6 所示。

圖5 第1級絕熱效率試驗(yàn)值與計(jì)算值

圖6 第2級絕熱效率試驗(yàn)值與計(jì)算值

從圖5中可見，第1級絕熱效率試驗(yàn)值總體上沿展向呈下降趨勢，由于轉(zhuǎn)子凸臺使轉(zhuǎn)子中部損失突增，葉展中部效率因此降低。絕熱效率計(jì)算值較好反映了試驗(yàn)值變化趨勢。從圖6中可見，第2級絕熱效率試驗(yàn)值葉展中部相對較高，葉尖和葉根部相對較低。效率計(jì)算值與試驗(yàn)值也較符合。對比損失圖和絕熱效率圖可見，各級損失試驗(yàn)值沿展向變化趨勢與效率試驗(yàn)值基本相反，說明級損失大小直接影響了級效率高低。由于正激波模型損失計(jì)算值更大，所以在2級全展向范圍內(nèi)，其效率計(jì)算結(jié)果基本低于雙激波模型的，且沿展向差值增大。

各級轉(zhuǎn)子出口相對氣流角和靜子出口氣流角試驗(yàn)值和計(jì)算值沿展向分布如圖7、8所示。

圖7 轉(zhuǎn)子出口相對氣流角試驗(yàn)值與計(jì)算值

圖8 靜子出口氣流角試驗(yàn)值與設(shè)計(jì)值

從圖7中可見，2級轉(zhuǎn)子出口相對氣流角試驗(yàn)值沿展向逐漸增大，2級轉(zhuǎn)子出口相對氣流角計(jì)算值和試驗(yàn)值比較接近，轉(zhuǎn)子根部計(jì)算值比試驗(yàn)值偏大。從圖8中可見，第1級靜子出口氣流角試驗(yàn)值沿展向變化較曲折，第2級靜子試驗(yàn)結(jié)果則基本呈軸向出氣。第1級靜子葉展中部出口氣流角計(jì)算值和試驗(yàn)值比較接近，葉根和葉尖部分差異較大。第2級靜子出氣角計(jì)算值比試驗(yàn)值小。2級轉(zhuǎn)子根部雙激波模型出口相對氣流角計(jì)算結(jié)果比正激波的大。2級靜子各展向位置，雙激波模型出口氣流角計(jì)算值比正激波模型的高 0.5°左右。

綜合來看，在近設(shè)計(jì)點(diǎn)損失預(yù)估和特性計(jì)算上，2種模型沒有明顯的優(yōu)劣區(qū)別，正激波模型計(jì)算結(jié)果能夠滿足工程需求，采用該模型更簡單有效。

4.3 近堵塞點(diǎn)特性計(jì)算

利用程序?qū)氯c(diǎn)進(jìn)行數(shù)值模擬，近堵塞點(diǎn)轉(zhuǎn)速為99.76%設(shè)計(jì)轉(zhuǎn)速，流量為83.866 kg/s。近堵塞點(diǎn)級累積增壓比和絕熱效率試驗(yàn)值和數(shù)值模擬計(jì)算結(jié)果對比見表5、6。

表5 近堵塞點(diǎn)累積增壓比試驗(yàn)值與計(jì)算值

表6 近堵塞點(diǎn)累積絕熱效率試驗(yàn)值與計(jì)算值 %

從表5、6中可見，相對于近設(shè)計(jì)點(diǎn)，近堵塞點(diǎn)總壓比和總效率試驗(yàn)結(jié)果有很大跌落，而單級性能的下降體現(xiàn)在第2級。第1級壓比和效率計(jì)算值較試驗(yàn)值均偏小。雙激波模型第1級壓比和效率計(jì)算結(jié)果較正激波模型的略低。正激波模型總壓比計(jì)算值比試驗(yàn)值約高1.63%，總效率計(jì)算值比試驗(yàn)值約高14.18%。而雙激波模型總壓比計(jì)算值跟試驗(yàn)值很吻合，總效率計(jì)算值比試驗(yàn)值約高7.28%。雙激波模型總壓比和總效率計(jì)算結(jié)果更接近試驗(yàn)值。

近堵塞點(diǎn)各級轉(zhuǎn)、靜子損失試驗(yàn)值和計(jì)算值沿展向分布如圖9、10所示。

圖9 第1級損失試驗(yàn)值與計(jì)算值

圖10 第2級損失試驗(yàn)值與計(jì)算值

從圖9中可見，第1級轉(zhuǎn)、靜子損失試驗(yàn)值沿展向變化趨勢與近設(shè)計(jì)點(diǎn)趨勢比較接近。第1級轉(zhuǎn)子尖部2種模型損失計(jì)算值很接近，而根部雙激波模型轉(zhuǎn)子損失計(jì)算值較試驗(yàn)值偏大，正激波模型轉(zhuǎn)子損失計(jì)算值跟試驗(yàn)值更接近。第1級靜子損失計(jì)算值跟試驗(yàn)值比較接近。從圖10中可見，第2級轉(zhuǎn)子損失試驗(yàn)值在葉尖、葉中和葉根局部有跌落，靜子葉根損失試驗(yàn)值有提高很大，葉尖損失試驗(yàn)值稍有增大。第2級轉(zhuǎn)子損失計(jì)算值沿展向變化不大，各展向位置雙激波模型計(jì)算值比正激波大約0.05，更接近試驗(yàn)值。第2級靜子損失計(jì)算值相比試驗(yàn)值有很大差異。

各級絕熱效率試驗(yàn)值和計(jì)算值沿展向分布如圖11、12所示。

圖11 第1級絕熱效率試驗(yàn)值與計(jì)算值

從圖11中可見，第1級絕熱效率試驗(yàn)值與近設(shè)計(jì)點(diǎn)試驗(yàn)值比較接近。第1級葉尖處兩種模型效率計(jì)算結(jié)果基本一致，正激波模型在葉根部分效率計(jì)算值與試驗(yàn)值很接近，雙激波模型由于葉根部分有較大的損失計(jì)算值，葉根效率計(jì)算值相對較低。從圖12中可見，第2級效率試驗(yàn)值沿展向變化較曲折，且在葉根處很低。第2級效率計(jì)算值比試驗(yàn)值高很多，沿展向變化較平滑。雙激波模型計(jì)算值比正激波模型低5%左右，更接近試驗(yàn)結(jié)果。

圖12 第2級絕熱效率試驗(yàn)值與計(jì)算值

各級轉(zhuǎn)子出口相對氣流角和靜子出口氣流角試驗(yàn)值和計(jì)算值沿展向分布如圖13、14所示。

圖13 轉(zhuǎn)子出口相對氣流角試驗(yàn)值與計(jì)算值

圖14 靜子出口氣流角試驗(yàn)值與計(jì)算值

從圖13中可見，第2級轉(zhuǎn)子出口相對氣流角試驗(yàn)值與近設(shè)計(jì)點(diǎn)試驗(yàn)值很接近。計(jì)算值比較接近試驗(yàn)值，在葉根部分比試驗(yàn)值偏大。從圖14中可見，第1級靜子出口氣流角試驗(yàn)值沿展向變化趨勢與近設(shè)計(jì)點(diǎn)趨勢仍比較接近，第2級靜子出口氣流角試驗(yàn)值在葉尖處接近0°，而在葉中到葉根部分基本為-10°。第1級靜子葉展中部出口氣流角計(jì)算值跟試驗(yàn)值很貼近，在葉尖和葉根處差別較大。第2級靜子出口氣流角計(jì)算結(jié)果基本比試驗(yàn)結(jié)果大。

綜合來看，利用雙激波模型計(jì)算近堵塞工況，損失計(jì)算結(jié)果比正激波模型的大，總壓比和總效率結(jié)果更接近試驗(yàn)值。

4.4 全工況特性計(jì)算

分別對95%、100%、110%設(shè)計(jì)轉(zhuǎn)速下工況進(jìn)行了計(jì)算，得到了特性曲線，并與試驗(yàn)曲線進(jìn)行了比較，如圖15、16所示。

圖15 增壓比全工況特性試驗(yàn)值與計(jì)算值

圖16 絕熱效率全工況特性試驗(yàn)值與計(jì)算值

從圖15中可見，在等轉(zhuǎn)速下較低流量下，雙激波模型壓比計(jì)算值比正激波模型高。在95%、100%設(shè)計(jì)轉(zhuǎn)速線近堵塞點(diǎn)處2種模型計(jì)算值很接近，但110%設(shè)計(jì)轉(zhuǎn)速線近堵塞點(diǎn)處雙激波模型計(jì)算值與試驗(yàn)值吻合更好。雙激波模型100%設(shè)計(jì)轉(zhuǎn)速線更接近試驗(yàn)曲線。從圖16中可見，雙激波模型等轉(zhuǎn)速效率計(jì)算值隨流量變換更陡峭，而正激波模型計(jì)算值變化比較平緩。各轉(zhuǎn)速線近堵塞點(diǎn)處，雙激波模型計(jì)算值有較大跌落，比較符合試驗(yàn)結(jié)果，而正激波模型計(jì)算值下降不大。雙激波模型100%設(shè)計(jì)轉(zhuǎn)速線也更接近試驗(yàn)曲線。

綜合來看，雙激波模型的各轉(zhuǎn)速特性線更接近試驗(yàn)結(jié)果，尤其在絕熱效率特性圖上，正激波模型效率計(jì)算值隨流量變化不夠明顯，而雙激波模型效率變化趨勢跟試驗(yàn)值符合較好。

5 結(jié)論

本文利用流線曲率法結(jié)合正激波和改進(jìn)的雙激波模型對某型2級風(fēng)扇進(jìn)行了數(shù)值模擬，得到了近設(shè)計(jì)點(diǎn)、近堵塞點(diǎn)和全工況的特性，結(jié)論如下：

（1）本文發(fā)展的性能預(yù)估方法能有效計(jì)算跨聲速軸流風(fēng)扇/壓氣機(jī)性能和氣動參數(shù)沿展向分布的規(guī)律，對設(shè)計(jì)點(diǎn)特性計(jì)算精度較高。

（2）不同激波模型對激波損失預(yù)估有較大差異，損失計(jì)算值會影響到增壓比和絕熱效率的預(yù)估，尤其是絕熱效率受損失計(jì)算結(jié)果影響很大。

（3）對于跨聲速風(fēng)扇/壓氣機(jī)，在近設(shè)計(jì)點(diǎn)激波損失預(yù)估和特性計(jì)算上，正激波模型基本滿足工程需求。而改進(jìn)的雙激波模型在非設(shè)計(jì)點(diǎn)損失預(yù)估和特性計(jì)算上結(jié)果更準(zhǔn)確。

要更準(zhǔn)確地計(jì)算非設(shè)計(jì)點(diǎn)的損失和特性參數(shù)，還需要研究更有效的二次流損失和葉型損失模型，需要在今后工作中進(jìn)一步完善。

[1]胡江峰,竺曉程,杜朝輝.一種跨聲速軸流壓氣機(jī)性能預(yù)測的數(shù)值方法[J].航空動力學(xué)報(bào),2011,26（1）:122-127.HU Jiangfeng,ZHU Xiaocheng,DU Zhaohui.Numerical method of transonic axial compressor performance prediction [J].Journal of Aerospace Power,2011,26（1）:122-127.（in Chinese）

[2]劉波,馬乃行,楊小東,等.適應(yīng)較大葉型彎角范圍的軸流壓氣機(jī)落后角模型[J].航空動力學(xué)報(bào),2014,29（8）:1824-1831.LIU Bo,MA Naixing,YANG Xiaodong,et al.Deviation angle models suitable for wider range of blade profile camber in axial-flow compressor[J].Journal of Aerospace Power,2014,29（8）:1824-1831.（in Chinese）

[3]祝啟鵬,高麗敏,李瑞宇,等.跨聲速多級軸流壓氣機(jī)特性預(yù)估及分析[J].推進(jìn)技術(shù),2014,35（10）:1342-1348.ZHU Qipeng,GAO Limin,LI Ruiyu,et al.Performance prediction and analysis of multi-stage transonic axial compressors[J].Journal of Propulsion Technology,2014,35（10）:1342-1348.（in Chinese）

[4]Novak R A.Streamline Curvature Computing Procedure for Fluid-Flow Problems[J].ASME Journal of Engineering for Power,1967,89（4）:478-490.

[5]袁寧,張振家,顧中華,等.渦噴發(fā)動機(jī)壓氣機(jī)三種S2流面計(jì)算程序的比較[J].推進(jìn)技術(shù),1998,19（1）:50-56.YUAN Ning,ZHANG Zhenjia,GU Zhonghua,et al.A comparison of three kinds of calculation program of S2 stream surface in the compressor of aero-engine[J].Journal of Propulsion Technology，1998,19（1）:50-56.（in Chinese）

[6]Miller G R,Lewis G W,Hartmann M J.Shock losses in transonic compressor blade rows[J].ASME Journal of Engineering for Power,1961,83（2）:235-242.

[7]趙勇,胡駿,鄭大勇.用雙激波模型計(jì)算風(fēng)扇/壓氣機(jī)非設(shè)計(jì)點(diǎn)的性能[J].南京航空航天大學(xué)學(xué)報(bào),2006,38（5）:545-550.ZHAO Yong,HU Jun,ZHENG Dayong.Dual-shock loss model for off-design performance calculation of fan/compressor[J].Journal of Nanjing University of Aeronautics&Astronautics,2006,38（5）:545-550.（in Chinese）

[8]LI Bo,GU Chunwei,LI Xiaotang,et al.Development and application of a through flow method for high-loaded axial flow compressors[J].Science China Technology Sciences,2016,59（1）:93-108.

[9]Bloch G S.Flow Losses in supersonic compressor cascades[D].Virginia:Virginia Polytechnic Institute and State University,1996.

[10]Boyer K M.An Improved streamline curvature approach for off-design analysis of transonic compression systems[D].Virginia:Virginia Polytechnic Institute and State University,2001.

[11]Biollo R,Benini E.Recent advances in transonic axial compressor aerodynamics[J].Progress in Aerospace Science,2013,56:1-18.

[12]閆轉(zhuǎn)運(yùn),成金鑫,陳江.多級軸流壓氣機(jī)通流造型一體化設(shè)計(jì)研究[J].工程熱物理學(xué)報(bào),2016,37（6）:1218-1224.YAN Zhuanyun, CHENG Jinxin, CHEN Jiang. Integrated through-flow and blade shape design of multi-stage axial flow compressor[J].Journal of Engineering Thermophysics,2016,37（6）:1218-1224.（in Chinese）

[13]Messenger H E,Kennedy E E.Two-stage fan I.aerodynamic and mechanical design[R].NASA-CR-120859.

[14]Ruggeri R S,Benser W A.Performance of a highly loaded two-stage axial-flow fan[R].NASA-TM-X-3076.

[15]Johnsen I A,Bullock R O.Aerodynamic design of axial-flow compressors[R].NASA-SP-36.

[16]Herrig L J,Emery J C,Erwin J R.Systematic two-dimensional cascade tests of NACA 65-Series compressors[R].NACA-RM-L51G31.

[17]鄧熙,劉波,馬乃行.高亞聲速大彎角軸流壓氣機(jī)平面葉柵損失模型研究[J].推進(jìn)技術(shù),2015,36（9）:1302-1308.DENG Xi,LIU Bo,MA Naixing.Investigation of loss model applicable to large range of high subsonic cascade in axial-flow compressor[J].Journal of Propulsion Technology,2015,36（9）:1303-1308.（in Chinese）

[18]杜文海,吳虎,阮建剛.跨聲速軸流壓氣機(jī)特性計(jì)算用激波損失模型[J].推進(jìn)技術(shù),2007,28（4）:378-382.DU Wenhai,WU Hu,RUAN Jiangang.A shock loss model for the transonic axial compressor[J].Journal of Propulsion Technology,2007,28（4）:378-382.（in Chinese）

[19]K?nig W M,Hennecke D K,Fottner L.Improved Blade Profile loss and deviation angle models for advanced transonic compressor bladings:Part II-A model for supersonic flow[J].ASME Journal of Turbomachinery,1996,117（1）:81-87.

[20]Klapproth J F.Approximate relative-total-pressure Losses of an Infinite cascade of supersonic blades with finite leading-edge thickness[R].NACA-RM-E9L21.

Study on Performance Simulation and Shock Model for Axial-Flow Compressor

JIANG Zhu-yu1,FAN Zhao-lin1,LIU Bo2
（1.China Aerodynamics Research and Development Center,Mianyang Sichuan 621000,China;2.School of Power and Energy,Northwestern Polytechnical University,Xi'an 710072,China）

In order to quickly and accurately predict performance and shock loss for axial-compressor,based on axial-flow compressor S2 flow surface streamline curvature method,simulation of the performance of a two-stage transonic fan was performed involving normal shock model and improved dual-shock model.Thus the overall performance and aerodynamic parameters on near-design and near-choke point,separately at 100%and 99.76%design rotation speed,and characteristic curves at 95%,110%,and 110%design speed were obtained.The differences in prediction of shock loss and calculation of fan/compressor performance between the two shock loss models have been analyzed by comparing the simulation results and experiment data.The results show that on the near-design point,with the normal shock model,calculation result of radial distribution of loss is close to test result,and results of total pressure ratio and total efficiency are proximately lower than test result by 1.96%and 2.54%separately.Thus,the normal shock model can satisfy engineering need.While on the near-choke point,with improved dual-shock shock model,calculation result of total loss is closer to test result,result of total pressure ratio is consistent with test result,and result of total efficiency is about 7.28%higher than test result.And calculated characteristic curves are more similar to test curves.Therefore the use of improved dual-shock loss model can provide a more accurate calculation of shock loss and performance parameters on far off-design point.

axial-flow compressor;streamline curvature calculation;shock model;loss prediction;performance simulation

V 231.3

A

1 0.1 3477/j.cnki.aeroengine.201 7.02.01 4

2016-11-10 基金項(xiàng)目：國家自然科學(xué)基金（51676162）資助

蔣筑宇（1991），男，在讀碩士研究生，研究方向?yàn)槿~輪機(jī)氣動熱力學(xué)；E-mail：kongjianglieren@163.com。

蔣筑宇，范召林，劉波.軸流壓氣機(jī)特性計(jì)算及激波模型研究[J].航空發(fā)動機(jī)，2017，43（2）：48-55.JIANG Zhuyu,FAN Zhaolin,LIU Bo.Study on performance simulation and shock mdelforaxial-flow compressor[J].Aeroengine，2017，43（2）：48-55.

（編輯：李華文）