賀　象，銀越千，黃生勤，馬宏偉，魏　巍

（1.中國航空機(jī)械動力研究所，湖南株洲412002；2.北京航空航天大學(xué)能源與動力工程學(xué)院，北京100191）

低速軸流壓氣機(jī)旋轉(zhuǎn)失速演化機(jī)制研究

賀象1，銀越千1，黃生勤1，馬宏偉2，魏巍2

（1.中國航空機(jī)械動力研究所，湖南株洲412002；2.北京航空航天大學(xué)能源與動力工程學(xué)院，北京100191）

為深入探討軸流壓氣機(jī)旋轉(zhuǎn)失速機(jī)理，利用全環(huán)非定常數(shù)值模擬在出口邊界添加節(jié)流閥模型模擬壓氣機(jī)節(jié)流過程，通過試驗(yàn)測量與數(shù)值模擬結(jié)果對比分析，澄清了低速軸流壓氣機(jī)失速先兆及失速團(tuán)的產(chǎn)生和演化過程。研究結(jié)果表明：壓氣機(jī)在轉(zhuǎn)子葉片安裝角存在幾何偏差的情況下，spike型失速先兆是由2階模態(tài)的整階擾動演化而來，并最終直接形成2個(gè)大尺寸全葉高的失速團(tuán)；與軸對稱轉(zhuǎn)子模型相比，考慮存在幾何偏差的轉(zhuǎn)子模型的數(shù)值模擬結(jié)果更為接近試驗(yàn)測量結(jié)果。

軸流壓氣機(jī)；全環(huán)非定常數(shù)值模擬；節(jié)流閥模型；失速先兆；失速團(tuán)；航空發(fā)動機(jī)

0　引言

旋轉(zhuǎn)失速會造成壓氣機(jī)性能急劇下降，嚴(yán)重時(shí)可能對壓氣機(jī)造成機(jī)械損傷，對發(fā)動機(jī)安全具有極大危害。因此，弄清其產(chǎn)生機(jī)理和發(fā)展過程一直是重要的研究方向。

20世紀(jì)90年代的大量研究表明，失速團(tuán)形成之前存在失速先兆，一般可分為2類：模態(tài)波先兆［1］和突尖波先兆［2-3］。此外，還存在高頻失速、整階擾動等先兆擾動形式［3］。整階擾動的周向傳播速度和方向與壓氣機(jī)旋轉(zhuǎn)速度相同，Day等［3］認(rèn)為該擾動的形成與壓氣機(jī)幾何偏差（幾何不對稱性）存在很大關(guān)系。

壓氣機(jī)和渦輪在實(shí)際加工和裝配過程中難免出現(xiàn)偏差，一般可分為2類：1類是固定在絕對坐標(biāo)系下的機(jī)匣圓度、機(jī)匣與轉(zhuǎn)子同心度偏差造成的葉尖間隙不均勻，其不但會造成壓氣機(jī)性能降低，在近失速點(diǎn)還會產(chǎn)生短尺度非定常擾動，并可進(jìn)一步發(fā)展成spike失速先兆［4-6］；另1類是固定在相對坐標(biāo)系下的轉(zhuǎn)子葉型的幾何偏差，該情況形成了Day等［3］所發(fā)現(xiàn)的整階擾動。其在很多壓氣機(jī)上均有不同程度的體現(xiàn)，不僅出現(xiàn)在近失速點(diǎn)［6-7］，還可能出現(xiàn)在遠(yuǎn)離穩(wěn)定邊界的設(shè)計(jì)點(diǎn)［8］，而且在整機(jī)環(huán)境的壓縮部件試驗(yàn)中也觀察到這種擾動形式［9］。

近年來，采用3維黏性多通道甚至全環(huán)的非定常數(shù)值模擬也可以捕捉到壓氣機(jī)失速先兆以及旋轉(zhuǎn)失速的發(fā)展過程。Hoying［10］首次捕捉到突變型失速先兆擾動，指出泄漏渦是否溢出相鄰葉片前緣可作為壓氣機(jī)是否失速的判段依據(jù)，Vo等［11］進(jìn)一步完善了該判斷準(zhǔn)則。隨后，Gourdian等［12］和Chen等［13］分別針對低速和高速軸流壓氣機(jī)進(jìn)行全環(huán)非定常數(shù)值模擬，均捕捉到短尺度擾動（spike先兆）的形成以及失速團(tuán)的演化過程，但是計(jì)算結(jié)果與試驗(yàn)結(jié)果仍存在較大差異。Gourdian指出［12］，在將來的研究中應(yīng)采用接近真實(shí)情況的幾何模型（進(jìn)、出段長度和幾何偏差等），才能得到更準(zhǔn)確的模擬結(jié)果。Choi等［14］改變了某個(gè)轉(zhuǎn)子葉片的安裝角并進(jìn)行節(jié)流失速全環(huán)非定常數(shù)值模擬，研究表明，小失速團(tuán)最先產(chǎn)生于該葉片附近的尖區(qū)；Pullan等［15］也調(diào)整了某個(gè)轉(zhuǎn)子葉片安裝角，數(shù)值模擬結(jié)果捕捉到了整階擾動，但是很難界定失速先兆是否由整階擾動發(fā)展而來。因此，存在幾何偏差條件下的壓氣機(jī)整階擾動與失速先兆之間關(guān)聯(lián)性仍存在許多值得考究的地方，需要進(jìn)一步深入研究。

本文結(jié)合試驗(yàn)與數(shù)值模擬方法，對比分析研究了轉(zhuǎn)子葉片安裝角均勻和非均勻2種情況下低速軸流壓氣機(jī)旋轉(zhuǎn)失速先兆及失速團(tuán)的產(chǎn)生和演化機(jī)制。

1　試驗(yàn)裝置及測量方案

試驗(yàn)裝置為北京航空航天大學(xué)低速大尺寸軸流壓氣機(jī)，機(jī)匣直徑1 m，輪轂比0.6，葉片數(shù)17，葉尖間隙3mm（1.5%葉高）。試驗(yàn)在折合轉(zhuǎn)速n=900 r/min下進(jìn)行。

壓氣機(jī)流量管周向布置4個(gè)靜壓測點(diǎn)測量進(jìn)口平均靜壓，根據(jù)大氣總壓計(jì)算進(jìn)口速度V，從而得到流量系數(shù)Φ=V/Utip。其中Utip為葉尖切向速度。在轉(zhuǎn)子上、下游機(jī)匣壁面周向分別布置4個(gè)測點(diǎn)測量平均靜壓，計(jì)算壓氣機(jī)靜壓升系數(shù)ψ=△Ps/0.5，從而測量壓氣機(jī)靜壓升特性。

在轉(zhuǎn)子上游10%弦長機(jī)匣壁面，周向均布8個(gè)動態(tài)壓力測點(diǎn)；在轉(zhuǎn)子葉頂機(jī)匣壁面弦向布置14個(gè)測點(diǎn)。數(shù)據(jù)采集系統(tǒng)多路并行連續(xù)采集動態(tài)壓力信號和鎖相信號，采樣速率為20 kHz。試驗(yàn)測量記錄多個(gè)穩(wěn)定工況點(diǎn)，隨后在近失速點(diǎn)快速關(guān)閉節(jié)流閥，壓氣機(jī)進(jìn)入完全失速狀態(tài)，記錄失速過程動態(tài)壓力信號。

2　數(shù)值模擬方案

在壓氣機(jī)分解過程中對轉(zhuǎn)子幾何參數(shù)進(jìn)行計(jì)量，發(fā)現(xiàn)該壓氣機(jī)轉(zhuǎn)子存在明顯的幾何偏差［16］，其中由安裝角不均勻產(chǎn)生的影響最大。本文壓氣機(jī)數(shù)值計(jì)算模型包括轉(zhuǎn)子葉片安裝角周向均勻和非均勻2種情況。對于非均勻情況，依據(jù)轉(zhuǎn)子實(shí)際計(jì)量的安裝角分布，在轉(zhuǎn)子的4個(gè)周向位置共7個(gè)葉片調(diào)整了安裝角大小如圖1所示，得到轉(zhuǎn)子葉片安裝角周向分布的波數(shù)為4。采用ICEM生成的六面體結(jié)構(gòu)網(wǎng)格如圖2所示，計(jì)算模型總網(wǎng)格數(shù)為865萬。其中進(jìn)、出口計(jì)算域長度分別為0.5、1倍流道外徑，網(wǎng)格分別為200萬和250萬，設(shè)置為靜止域。轉(zhuǎn)動域計(jì)算網(wǎng)格拓?fù)錇镠OH，單通道網(wǎng)格數(shù)約為24萬，整圈網(wǎng)格總數(shù)為415萬。

圖1　轉(zhuǎn)子周向葉片安裝角分布

圖2　低速大尺寸軸流壓氣機(jī)數(shù)值模擬計(jì)算網(wǎng)格

采用CFX12.0求解3維黏性雷諾平均N-S方程，模擬壓氣機(jī)全環(huán)非定常流場。空間離散采用2階格式，湍流模型采用κ-ε。標(biāo)準(zhǔn)壁面函數(shù)處理近壁區(qū)域流動，壁面最大y+約為20，滿足壁面函數(shù)的要求。進(jìn)口邊界條件給定標(biāo)準(zhǔn)大氣壓、軸向進(jìn)氣，湍流度為5%。轉(zhuǎn)、靜摻混面選擇Transient Rotor Stator。物理時(shí)間步的選擇保證轉(zhuǎn)子每轉(zhuǎn)過1個(gè)通道的步數(shù)均為20。出口邊界條件設(shè)置如下：

（1）穩(wěn)定工況點(diǎn)非定常數(shù)值模擬

出口邊界條件給定平均靜壓。

收斂準(zhǔn)則：當(dāng)監(jiān)測的轉(zhuǎn)子進(jìn)口壓力信號具有明顯的周期性，并且在500個(gè)物理時(shí)間步內(nèi)，出口流量脈動不大于0.5%，即認(rèn)為此工況點(diǎn)非定常數(shù)值收斂。

（2）節(jié)流失速過程的數(shù)值模擬

出口平均靜壓 Ps，out根據(jù)節(jié)流閥模型［17］確定，定義為

式中：Pt為參考大氣壓；密度ρ為常數(shù)；k0為常數(shù)；k1為閥門開度，隨計(jì)算時(shí)間步的增大而減小。

3　壓氣機(jī)靜壓升特性

對比安裝角非均勻和均勻情況數(shù)值模擬獲得的壓氣機(jī)特性如圖3所示。從圖中可見，轉(zhuǎn)子安裝角非均勻引起靜壓升系數(shù)減小，壓氣機(jī)性能衰減，而且數(shù)值模擬捕捉到的非均勻情況失速點(diǎn)流量比均勻情況的更大。對比數(shù)值模擬與試驗(yàn)測得的特性曲線，非均勻情況數(shù)值模擬捕捉到的失速點(diǎn)流量與試驗(yàn)測得的一致，而與均勻情況數(shù)值模擬捕捉到的失速點(diǎn)流量偏差較大。表明采用更接近真實(shí)壓氣機(jī)幾何情況的計(jì)算模型進(jìn)行非定常數(shù)值模擬，所獲得的節(jié)流失速特性更接近試驗(yàn)測量結(jié)果。

圖3　靜壓升特性

4　壓氣機(jī)節(jié)流失速過程

在轉(zhuǎn)子上游10%弦長機(jī)匣壁面的傳感器測得的節(jié)流失速過程動態(tài)壓力信號分布和信號小波變換［16］能量譜分別如圖4（a）、（b）所示。對于圖4（b），橫坐標(biāo)為時(shí)間，用轉(zhuǎn)子轉(zhuǎn)過的轉(zhuǎn)數(shù)表示，縱坐標(biāo)為無量綱頻率，1 frot為轉(zhuǎn)頻，17為葉片通過頻率（BPF），譜圖的幅值代表不同頻率分量在不同時(shí)刻下的能量高低。在第70轉(zhuǎn)之前，除BPF之外，2 frot整階擾動［16］最強(qiáng)，4 frot整階擾動次之。

圖4　壓氣機(jī)節(jié)流失速過程無量綱壓力信號及其小波譜

整階擾動向失速先兆演變過程如圖5所示。其中P1信號上虛線圈M標(biāo)示了a和b 2個(gè)2 frot整階擾動，其傳播速度與轉(zhuǎn)速相同，與圖4（b）中M所標(biāo)出的2 frot對應(yīng)。在第70轉(zhuǎn)左右，整階的2 frot擾動向1.5 frot擾動演化，該過程實(shí)際是a和b 2個(gè)擾動周向傳播速度降低、尺度和強(qiáng)度增大、向失速先兆擾動演化的過程。在第73轉(zhuǎn)時(shí)刻其傳播速度約為82%。再經(jīng)歷10多轉(zhuǎn)的時(shí)間，2個(gè)先兆擾動尺度和強(qiáng)度明顯增大，最終形成2個(gè)大尺度失速團(tuán)，如圖6所示。其傳播速度為65%轉(zhuǎn)速，壓氣機(jī)進(jìn)入完全失速狀態(tài)。由此可知，該壓氣機(jī)失速先兆波是由2 frot（模態(tài)數(shù)為2）發(fā)展而來的，并最終形成了2個(gè)失速團(tuán)。

圖5　整階擾動向失速先兆演變過程

圖6　失速先兆和失速團(tuán)的周向傳播

在轉(zhuǎn)子安裝角非均勻情況壓氣機(jī)節(jié)流失速過程數(shù)值模擬監(jiān)測的壓力信號及其小波譜分別如圖7（a）、（b）所示。虛擬探針布置位置與試驗(yàn)相同。壓氣機(jī)先兆擾動產(chǎn)生之前，除BPF和1 frot擾動強(qiáng)度較大之外，4 frot整階擾動強(qiáng)度也較大。對于圖7（b），字母a和b對應(yīng)4 frot整階擾動，這與壓氣機(jī)轉(zhuǎn)子計(jì)算模型在周向4個(gè)位置調(diào)整了安裝角存在直接關(guān)系（安裝角周向分布的波數(shù)為4）。其傳播速度與轉(zhuǎn)速相同，傳播特征與近失速點(diǎn)的情況相同［16］，且隨壓氣機(jī)節(jié)流擾動強(qiáng)度升高逐漸增大。而在接下來的1轉(zhuǎn)中，表示能量較高斑點(diǎn)的字母c向低頻方向傾斜，表明4 frot整階擾動向失速先兆波演化，其尺度增大、傳播速度降低，到第19.5轉(zhuǎn)左右，周向傳播速度下降至83%轉(zhuǎn)速。在第20轉(zhuǎn)左右時(shí)刻，C5、C7這2個(gè)周向位置的原始信號上能夠清晰地觀察到失速先兆波，分別用黑色和紅色虛線箭頭標(biāo)識了這2個(gè)先兆擾動的周向傳播路徑，如圖8所示。在C1信號上，N所標(biāo)示的先兆擾動與圖7中的擾動相對應(yīng)。隨后，2個(gè)先兆擾動尺度進(jìn)一步增大，傳播速度進(jìn)一步降低，演化成2個(gè)小失速團(tuán)。在失速團(tuán)形成初期，其周向傳播速度約為75%轉(zhuǎn)速，進(jìn)入完全失速狀態(tài)之后，形成2個(gè)大尺度的失速團(tuán)，周向傳播速度降至67%轉(zhuǎn)速。

圖7　轉(zhuǎn)子安裝角非均勻情況，數(shù)值模擬監(jiān)測節(jié)流失速過程的壓力信號及其小波譜

圖8　轉(zhuǎn)子安裝角非均勻情況，數(shù)值模擬監(jiān)測的失速先兆周向傳播

轉(zhuǎn)子安裝角均勻情況節(jié)流失速過程數(shù)值模擬監(jiān)測的壓力信號如圖9所示，虛擬探針布置位置與試驗(yàn)相同。在第8.8轉(zhuǎn)左右之前，原始壓力信號及其小波譜上均觀察不到整階擾動。之后壓氣機(jī)周向多個(gè)通道突然同時(shí)產(chǎn)生失速先兆擾動，并合并成3個(gè)小失速團(tuán)，如圖9中A、B和C所標(biāo)示傳播路徑。在第12.5轉(zhuǎn)左右時(shí)，失速團(tuán)A被失速團(tuán)C吞并，最終形成2個(gè)失速團(tuán)，傳播速度為63%轉(zhuǎn)速，壓氣機(jī)進(jìn)入完全失速狀態(tài)。

圖9　轉(zhuǎn)子安裝角均勻情況，數(shù)值模擬監(jiān)測的失速先兆和失速團(tuán)周向傳播

經(jīng)比較可知，在安裝角非均勻情況下，數(shù)值模擬獲得的節(jié)流失速過程與試驗(yàn)測量結(jié)果類似，2個(gè)失速先兆是由整階擾動演化而來，并最終直接形成2個(gè)以約65%轉(zhuǎn)速沿周向傳播的大尺度失速團(tuán)。而在轉(zhuǎn)子安裝角均勻情況下，周向存在多個(gè)失速先兆，雖然同樣最終形成周向傳播速度為63%轉(zhuǎn)速的2個(gè)大失速團(tuán)，但失速先兆和失速團(tuán)的產(chǎn)生和演化過程與試驗(yàn)結(jié)果偏差較大。因此，考慮了轉(zhuǎn)子葉片幾何偏差的壓氣機(jī)節(jié)流失速數(shù)值模擬獲得的結(jié)果顯然更為接近試驗(yàn)測量結(jié)果。

5　葉尖泄漏流動與失速先兆的關(guān)系

對比分析數(shù)值模擬和試驗(yàn)測量的壓氣機(jī)節(jié)流失速過程的靜壓信號之后，進(jìn)一步剖析失速先兆產(chǎn)生和發(fā)展過程中轉(zhuǎn)子葉尖流場的演化過程。

在轉(zhuǎn)子通道前緣附近熵增分布中，梯度較大的位置即高熵增區(qū)與主流之間的分界線，可表示泄漏渦影響區(qū)與主流之間的分界線［11］。98%葉高的壓氣機(jī)葉頂葉片前緣發(fā)生溢流導(dǎo)致失速先兆如圖10所示。從圖中可見，在失速先兆產(chǎn)生時(shí)，安裝角均勻情況下，周向至少4個(gè)通道的泄漏流溢出相鄰葉片前緣，而在安裝角非均勻情況下，周向僅有2個(gè)位置存在泄漏流溢出相鄰葉片前緣的現(xiàn)象（8、13號葉片），其它通道泄漏流與主流的交界線并沒有超出額線。泄漏流溢出相鄰葉片前緣的流動現(xiàn)象符合Hoying、Vo等［10-11］提出的突變型失速先兆產(chǎn)生條件，可見數(shù)值模擬獲得的失速起始擾動為spike失速先兆。

圖10　98%葉高的壓氣機(jī)葉頂葉片前緣發(fā)生溢流導(dǎo)致失速先兆

試驗(yàn)測量結(jié)果表明壓氣機(jī)失速先兆是由具有模態(tài)特征的2 frot演化而來，試驗(yàn)所觀察到的失速先兆屬于spike還是modal先兆，需要進(jìn)一步對比分析。

在失速先兆產(chǎn)生和演化過程中的轉(zhuǎn)子葉頂流動如圖11所示。從圖11（a）中可見，A1為整階擾動，其特征是2個(gè)相鄰?fù)ǖ漓o壓明顯比兩側(cè)通道高，即與時(shí)序信號中所觀察到的壓力突增相對應(yīng)。

圖11　失速先兆產(chǎn)生和演化過程中的轉(zhuǎn)子葉頂流動

從圖11（b）中可見，隨壓氣機(jī)節(jié)流，整階擾動向失速先兆演化。此時(shí)A1位置的2個(gè)通道壓升也比較高，與圖11（a）中的類似。需要注意的是，在A1右側(cè)的B1所標(biāo)示的2個(gè)通道的前緣位置表示泄漏渦的低靜壓區(qū)延伸到相鄰葉片的壓力面，整個(gè)通道前緣附近的靜壓均非常低。即失速先兆的特征為高壓區(qū)A1和低壓區(qū)B1組合流場結(jié)構(gòu)。

從圖11（c）中可見，由數(shù)值模擬獲得的spike先兆流動（1～3葉片間A1、B1通道），與圖11（b）中A1、B1流動結(jié)構(gòu)相似。由于A1通道中間存在較大尺度的低速旋渦結(jié)構(gòu)（速度矢量圖），堵塞轉(zhuǎn)子通道導(dǎo)致靜壓增大（靜壓云圖），而相鄰的B1通道存在向相鄰葉片前緣溢出的流動現(xiàn)象（速度矢量圖），即形成spike失速先兆。由于存在溢流現(xiàn)象，B1通道流動速度并不低，因而靜壓相對較低（靜壓云圖）。從圖8可知，失速先兆沿著周向傳播速度為75%轉(zhuǎn)速，從相對坐標(biāo)系下看，在接下來的某個(gè)時(shí)刻，溢出的流動與3號葉片泄漏渦相互作用形成大尺度旋渦，然后向4號葉片壓力面移動，并溢出4號葉片前緣，即失速先兆從3號葉片前緣（B1通道）傳播至4號葉片。在該過程中，B1通道的旋渦尺度進(jìn)一步增大，通道堵塞也增大，形成類似A1通道流動狀態(tài)，因而B1通道靜壓增大。同時(shí)，A1通道中低速旋渦流出通道，流動恢復(fù)正常。

從圖11（d）中可見，隨著壓氣機(jī)節(jié)流，失速先兆在周向傳播的過程中，隨著尺度的不斷增大，其傳播速度也會降低，趨于形成小失速團(tuán)。A1、B1所示的失速先兆的影響范圍達(dá)到3～4個(gè)通道（前緣低靜壓區(qū)），這與童志庭［18］所測得“spike”失速先兆流動結(jié)構(gòu)相同。數(shù)值模擬同樣捕捉到類似的流動結(jié)構(gòu)，如圖11（e）所示。

由以上試驗(yàn)和數(shù)值模擬結(jié)果對比分析可以確定，該壓氣機(jī)失速先兆形式為spike先兆，是由模態(tài)數(shù)為2的2 frot整階擾動演化而來。

6　完全失速狀態(tài)葉頂流場

在完全失速狀態(tài)下試驗(yàn)結(jié)果與數(shù)值模擬結(jié)果的對比以及熵增分布和速度矢量分布如圖12所示。圖中可以明顯觀察到周向存在的2個(gè)失速團(tuán)。

圖12　在完全失速狀態(tài)下數(shù)值模擬與試驗(yàn)結(jié)果的對比

從圖12（c）中可見，位于葉排上游的低靜壓區(qū)實(shí)際為順時(shí)針方向旋轉(zhuǎn)的回流旋渦。旋渦左側(cè)氣流逆著進(jìn)氣方向，因而滯止形成高靜壓區(qū)，左側(cè)主流偏轉(zhuǎn)，對于轉(zhuǎn)子葉片從相對坐標(biāo)系來看減小了其攻角；而在旋渦右側(cè)，主流趨向于向右側(cè)偏轉(zhuǎn)，從相對坐標(biāo)系來看增大了轉(zhuǎn)子葉片的攻角；從而，右側(cè)轉(zhuǎn)子通道分離增大，左側(cè)分離減小，分離區(qū)相對于葉片排向右側(cè)傳播，即形成相對轉(zhuǎn)子葉片按照轉(zhuǎn)速方向相反的方向轉(zhuǎn)動的旋轉(zhuǎn)失速現(xiàn)象。王洪偉［19］以及Lepicovsky［20］利用熱線風(fēng)速儀測量了失速團(tuán)內(nèi)部速度矢量，獲得了類似的位于失速團(tuán)內(nèi)的大尺度回流旋渦流動結(jié)構(gòu)。

從圖12（d）中可見失速團(tuán)的影響范圍，每個(gè)失速團(tuán)占據(jù)了約5個(gè)轉(zhuǎn)子通道的寬度。失速區(qū)之外的通道流動比較穩(wěn)定，如圖12（b）、（d）所示。由于失速區(qū)的流動堵塞，造成壓氣機(jī)周向流量重新分配，失速區(qū)流量迅速減小，其余通道流量有所增大流動趨于穩(wěn)定，從而使壓氣機(jī)適應(yīng)過小的流量狀態(tài)，而非試圖把流量的減小平均分配到所有的葉片通道。

轉(zhuǎn)子出口熵增如圖13所示。從圖中可見，在完全失速狀態(tài)下，壓氣機(jī)失速團(tuán)影響范圍擴(kuò)展至葉根，顯然這是1種全葉高失速。另外，數(shù)值模擬獲得的失速團(tuán)在周向并不是對稱分布，試驗(yàn)測量結(jié)果也顯示了類似的情況。

圖13　完全失速狀態(tài)下轉(zhuǎn)子出口熵增

7　結(jié)論

本文針對低速大尺寸軸流壓氣機(jī)，結(jié)合試驗(yàn)和全環(huán)非定常數(shù)值模擬方法，詳細(xì)地研究了在有幾何偏差條件下壓氣機(jī)失速先兆以及旋轉(zhuǎn)失速演化過程，得出以下結(jié)論：

（1）壓氣機(jī)在轉(zhuǎn)子周向非軸對稱情況下，會產(chǎn)生1種頻率與轉(zhuǎn)頻成整數(shù)倍、傳播速度與轉(zhuǎn)速相同的整階擾動；

（2）數(shù)值模擬結(jié)果表明，整階擾動階數(shù)與轉(zhuǎn)子幾何參數(shù)周向非均勻分布的波數(shù)一致；

（3）隨壓氣機(jī)節(jié)流，2階模態(tài)的整階擾動誘導(dǎo)轉(zhuǎn)子葉尖周向2個(gè)位置形成spike失速先兆，并最終演化成2個(gè)失速團(tuán)；

（4）采用更接近真實(shí)壓氣機(jī)幾何情況計(jì)算模型進(jìn)行非定常數(shù)值模擬，所獲得壓氣機(jī)特性以及節(jié)流失速過程更接近試驗(yàn)測量結(jié)果。

［1］McDougall N M，Cumpsty N A，Hynes T P.Stall inception in axial compressors［J］.Journal of Turbomachinery，1988，112（1）：116-125.

［2］Day I J.Stall inception in axial flow compressors［J］.Journal of Turbomachinery，1993，115（1）：1-9.

［3］Day I J，Breuer T，Escuret J，et al.Stall inception and the prospects for active control in four high-speed compressors［J］.Journal of Turbomachinery，1999，121（1）：18-27.

［4］Young A，Day I，Pullan G.Stall warning by blade pressure signature analysis［R］.ASME 2011-GT-45850.

［5］Cameron J D，Bennington M A，Ross M H，et al.Effects of steady tip clearance asymmetry and rotor whirl on stall inception in an axial compressor［R］.ASME 2007-GT-28278.

［6］LIN Feng，CHEN Jingyi，LI Meilin.Wavelet analysis of rotor-tip disturbances in an axial-flow compressor［J］.Journal of Propulsion and Power，2004，20（2）：319-333

［7］LIN Feng，CHEN Jingyi，LI Meilin.Experimental investigation of unsteady rotor tip flows in a high speed compressor throttled to stall［R］. ASME 2002-GT-30360.

［8］Hyung-Soo L，Young-Cheon L，Shin-Hyoung K，et al.Experimental analysis about the magnitude of the shaft frequency growth near stall in the axial compressor［J］.Journal of Mechanical Science and Technology，2011，25（5）：1317-1324.

［9］Hoss B，Leinhos D，F(xiàn)ottner L.Stall inception in the compressor system of a turbo fan engine［J］.Journal of Turbomachinery，1998，120（1）：32-44.

［10］Hoying D A，Tan C S，Vo H D，et al.Role of blade passage flow structures in axial compressor rotating stall inception［J］.Journal of Turbomachinery，1999，121（4）：735-742.

［11］Vo H D，Tan C S，Greitzer E M.Criteria for spike initiated rotating stall［J］.Journal of Turbomachinery，2005，130（1）：155-165.

［12］Gourdain N，Burguburu S，Leboeuf F，et al.Simulation of rotating stall in a whole stage of an axial compressor［J］.Computers&Fluids，2010，39（3）：1644-1655.

［13］Chen J P，Hathawary M D，Herrick G P.Prestall behavior of a transonic axial compressor stage via time-accurate numerical simulation［J］.Journal of Turbomachinery，2008，130（4）：353-368.

［14］Choi M，Vahdati M.Numerical strategies for capturing rotating stall in fan［J］.Proc.IMechE Part A：J.Power and Energy，2011，225（5）：655-664.

［15］Pullan G.，Young A，Day I，et al.Origins and structure of spike-type rotating stall［R］.ASME 2012-GT-68707.

［16］HE Xiang，MA Hongwei，ZHANG Jun，et al.Wavelet analysis of the shaft order perturbation and stall inception in an axial compressor［J］. Journal of Thermal Science，2013，22（3）：223-228.

［17］蔣康濤.低速軸流壓氣機(jī)旋轉(zhuǎn)失速的數(shù)值模擬研究［D］.北京：中國科學(xué)院研究生院，2004. JIANG Kangtao.Numerical investigation on rotating stall in a low-speed axial compressor［D］.Beijing：Graduate School of the Chinese Academy of Sciences，2004.（In Chinese）

［18］童志庭.軸流壓氣機(jī)中葉尖泄漏渦、失速先兆、葉尖微噴氣非定常關(guān)聯(lián)性的試驗(yàn)研究［D］.北京：中國科學(xué)院研究生院，2006. TONG Zhiting.The interactive unsteady mechanism between tip leakage vortex，stall inception and micro tip injection in low-speed axial compressor［D］.Beijing：Graduate School of the Chinese Academy of Sciences，2006.（In Chinese）

［19］王洪偉，宋志強(qiáng)，王治.單斜熱線探針測量壓氣機(jī)旋轉(zhuǎn)失速流場［J］.航空動力學(xué)報(bào)，2011，26（2）：249-254. WANG Hongwei，SONG Zhiqiang，WANG Zhi.Measurement of rotating stall in an axial compressor using a single slanted hot-wire probe［J］.Journal of Aerospace Power，2011，26（2）：249-254.（In Chinese）

［20］Lepicovsky J，Braunscheidel E P.Measurement of flow pattern within arotatingstallcellinanaxialcompressor［R］.ASME 2006-GT-91209.

（編輯：栗樞）

Investigation on the Rotating Stall in a Low-Speed Axial Compressor

HE Xiang1，YIN Yue-qian1，HUANG Sheng-qin1，MA Hong-wei2，WEI Wei2
（1.AVIC Aviation Powerplant Research Institute，Zhuzhou Hunan，412002，China；2.School of Energy and Power Engineering，Beihang University，Beijing 100191，China）

In order to further investigate the axial flow compressor rotating stall mechanism，the compressor stall process was simulated via the methods of full-annulus unsteady numerical simulation in which a throttle model was imposed on the outlet boundary.With the comparative analysis of the experimental and numerical results，the generation and evolution of the stall inception and the stall cell in a low-speed axial-compressor was clarified.The results indicate that two spike inceptions are induced by the shaft order perturbations of the second order and evolve into two full-span stall cells in this compressor in the case of some deviation in the rotor blade stagger angle. Comparing with the symmetric rotor model，the numerical results of the rotor model with geometry deviation agreed well with the experimental results.

axial compressor；full-annulus unsteady simulation；throttle model；stall inception；stall cell；aeroengine

V 231.3

Adoi：10.13477/j.cnki.aeroengine.2016.05.014

2016-03-27

賀象（1984），男，工程師，主要從事壓氣機(jī)氣動設(shè)計(jì)工作；E-mail：41213014@qq.com。

引用格式：賀象，銀越千，黃生勤，等.低速軸流壓氣機(jī)旋轉(zhuǎn)失速演化機(jī)制研究［J］.航空發(fā)動機(jī)，2016，42（5）：81-87.HEXiang，YINYueqian，HUANG Shengqin，etal.Investigationontherotatingstallinalow-speedaxialcompressor［J］.Aeroengine，2016，42（5）：81-87.