王春利，喬渭陽

(1.江西洪都航空工業(yè)集團(tuán)有限責(zé)任公司，江西南昌330024；2.西北工業(yè)大學(xué)動力與能源學(xué)院，陜西西安710072)

畸變進(jìn)氣對兩級風(fēng)扇穩(wěn)定性影響的數(shù)值模擬

王春利1，喬渭陽2

(1.江西洪都航空工業(yè)集團(tuán)有限責(zé)任公司，江西南昌330024；2.西北工業(yè)大學(xué)動力與能源學(xué)院，陜西西安710072)

基于多級軸流壓氣機(jī)的逐級特性，建立了一種預(yù)估多級軸流壓氣機(jī)在均勻進(jìn)氣和周向畸變進(jìn)氣條件下的喘振邊界的一維數(shù)值模擬方法。根據(jù)動態(tài)壓縮系統(tǒng)模型，對一臺兩級風(fēng)扇的喘振邊界進(jìn)行了數(shù)值預(yù)測，與基于李亞普諾夫理論的線化一維模型和試驗(yàn)結(jié)果的比較表明，該模型能較為準(zhǔn)確地預(yù)估多級軸流壓氣機(jī)的喘振邊界。對畸變進(jìn)氣條件下兩級風(fēng)扇穩(wěn)定性進(jìn)行的詳細(xì)數(shù)值分析表明：進(jìn)氣總壓畸變在流動過程中會生成總溫畸變并伴隨著總壓畸變的衰減，進(jìn)氣總溫畸變則會生成總壓畸變并伴隨著總溫畸變的衰減；反向總溫總壓組合畸變進(jìn)氣時，畸變衰減快穩(wěn)定裕度損失小，而正向總溫總壓組合畸變進(jìn)氣時，畸變衰減慢穩(wěn)定裕度損失大。

航空發(fā)動機(jī)；軸流壓氣機(jī)；一維模型；畸變；氣動穩(wěn)定性；數(shù)值模擬

1 引言

軍用航空燃?xì)鉁u輪發(fā)動機(jī)在研制和使用過程中，經(jīng)常遇到失穩(wěn)引起的發(fā)動機(jī)壓縮系統(tǒng)不穩(wěn)定工作問題。壓縮系統(tǒng)作為航空燃?xì)鉁u輪發(fā)動機(jī)的一個重要組成部分，正常工作條件下，為發(fā)動機(jī)提供合適的質(zhì)量流量和增壓比，最好能保持氣動穩(wěn)定性。然而，實(shí)際飛行過程中，特別是進(jìn)行機(jī)動飛行和武器發(fā)射操作時，壓縮系統(tǒng)通常要受到進(jìn)口流場波動和其他瞬態(tài)變化影響，這些都可能導(dǎo)致風(fēng)扇/壓氣機(jī)的惡劣工作狀態(tài)甚至失穩(wěn)(如旋轉(zhuǎn)失速或喘振)。

由于經(jīng)濟(jì)和技術(shù)條件的限制，通過試驗(yàn)測試所有氣流擾動對壓氣機(jī)的影響不現(xiàn)實(shí)，于是人們通過有效的數(shù)學(xué)模型來獲得。鑒于目前三維模型發(fā)展得還不夠成熟，基于準(zhǔn)一維、時間相關(guān)模型方程的發(fā)動機(jī)壓縮系統(tǒng)穩(wěn)定工作理論預(yù)測模型，即一維模型，也稱stage-by-stage模型，仍具有重要的應(yīng)用價值。

國外，美國阿諾德工程發(fā)展中心的Davis等利用平行壓氣機(jī)理論，采用直接數(shù)值求解的時間推進(jìn)計算方法，對畸變進(jìn)氣條件下的壓縮系統(tǒng)穩(wěn)定性進(jìn)行了較全面的研究[1～3]。而俄羅斯中央航空研究院，則采用基于李亞普諾夫穩(wěn)定性理論的計算方法，通過線化控制方程的方法，發(fā)展了針對畸變進(jìn)氣條件下的全臺發(fā)動機(jī)穩(wěn)定性分析程序。本課題組已采用這種方法進(jìn)行過相關(guān)的數(shù)值計算[4]。

國內(nèi)對此也進(jìn)行了相關(guān)研究，但研究得還很不夠。有學(xué)者基于Davis等發(fā)展的stage-by-stage模型，計算了均勻進(jìn)氣和各種畸變進(jìn)氣(包括組合畸變進(jìn)氣)條件下軸流壓氣機(jī)的穩(wěn)定裕度損失[5，6]，但只比較了不同畸變下的損失大小，并沒有分析畸變在壓縮系統(tǒng)內(nèi)部的發(fā)展變化和穩(wěn)定裕度損失機(jī)理，且其中文獻(xiàn)[6]所用的失穩(wěn)判據(jù)還不太合適。

本文參照Davis等的研究方法，將葉輪機(jī)對氣流的作用作為源項(xiàng)來考慮，應(yīng)用時間推進(jìn)求解初值問題的方法，對航空燃?xì)鉁u輪發(fā)動機(jī)的壓縮系統(tǒng)進(jìn)行數(shù)值模擬；并利用平行壓氣機(jī)理論，對畸變進(jìn)氣條件下的壓縮系統(tǒng)穩(wěn)定性進(jìn)行研究。還通過具體算例，深入研究了畸變在壓縮系統(tǒng)內(nèi)部的非定常變化和組合畸變進(jìn)氣下的穩(wěn)定裕度損失機(jī)理。

2 模型方程

如圖1(a)所示，發(fā)動機(jī)空氣壓縮系統(tǒng)由進(jìn)氣道和多級壓氣機(jī)組成。控制體的建立如圖1(b)所示，總控制體由壓氣機(jī)和進(jìn)出氣管道組成，包括氣流抽吸影響。其中Fx是軸向力分布，包括葉片和涵道壁面對氣流的作用力。外部對控制體提供的能量包括傳入的熱量Q和軸功SW，從控制體抽吸的氣體流量為WB。進(jìn)氣道入口邊界條件給定總壓和總溫，出口邊界條件可給定靜壓、流量或出口馬赫數(shù)。為更準(zhǔn)確地體現(xiàn)壓縮系統(tǒng)對氣流的作用，總控制體被劃分為一系列單元控制體，如圖1(c)所示。在壓氣機(jī)區(qū)域，一個小控制體包括一級靜子和轉(zhuǎn)子。進(jìn)氣道和排氣管道同樣被分為一些小控制體。圖中還示出了單元控制體上的進(jìn)出口參數(shù)。

圖1 壓縮系統(tǒng)物理模型和控制體概念Fig.1 Physical compression system model and control volume concepts

對每個單元控制體應(yīng)用質(zhì)量、動量和能量守恒定律，得到控制方程組。連續(xù)方程、動量方程和能量方程依次為：

式中：ρ為氣體密度，A為控制體截面積，W為質(zhì)量流率，IMP=Wu+psA為沖量函數(shù)，為葉片和機(jī)匣環(huán)壁面對氣流的軸向作用力，E為氣體單位體積的內(nèi)能，H為氣體總焓，HB為抽吸部分的熱焓。

對于軸向力、軸功等動量和能量方程中的源項(xiàng)，由壓氣機(jī)的準(zhǔn)穩(wěn)態(tài)級特性確定。級特性可由試驗(yàn)數(shù)據(jù)直接取得或采用級疊加法計算獲得，本文算例中采用計算獲得[7]。對于無葉通道，方程組中的力和功分別表示為：

對于逐級葉片通道，力和功可分別表示為：

式中：f為摩擦系數(shù)，D為水力直徑，ηp為絕熱效率，TR為總溫比。

3 計算方法

守恒型控制方程組的有限差分形式可表示為：

對上述氣動方程組，采用二階精度的MacCor?mack預(yù)估-校正顯式差分格式求解。計算過程中應(yīng)滿足CFL條件，即：

因?yàn)楸灸Ｐ图仁且粋€初值問題也是一個邊值問題，邊界條件處理不當(dāng)會給數(shù)值求解帶來麻煩，所以采用特征分析法處理邊界條件。在進(jìn)口應(yīng)用特征線差分格式，內(nèi)部控制體用MacCormack格式；出口邊界在非堵塞流動情況下用特征線法，在堵塞流動情況下用等熵噴管邊界條件。

將進(jìn)出口段控制體取為等截面管道，若不考慮摩擦和熱交換，根據(jù)一維非定常歐拉氣動方程組和相容性關(guān)系，可給出特征方程和相容性方程：

當(dāng)流動為亞聲速時，入口給定總壓p*和總溫T*，其他參數(shù)輔之特征、相容性方程求解；出口截面可給定靜壓p，再加上特征、相容性方程求解。

本研究發(fā)展的計算進(jìn)氣畸變對發(fā)動機(jī)穩(wěn)定性影響的程序，使用了經(jīng)典的平行壓氣機(jī)理論模型[8]。

失穩(wěn)判斷：利用這種動態(tài)壓縮系統(tǒng)模型的特性，根據(jù)文獻(xiàn)[1]中方法進(jìn)行逼喘，給定壓氣機(jī)出口靜壓升。計算時，為使預(yù)測結(jié)果足夠精確，靜壓的時間變化率應(yīng)充分慢。

4 算例分析

為檢驗(yàn)所建模型的實(shí)際應(yīng)用效果，對文獻(xiàn)[9]中的一臺兩級風(fēng)扇進(jìn)行數(shù)值計算。其壓縮系統(tǒng)劃分為風(fēng)扇進(jìn)氣管道(劃分為三個控制體)、兩個風(fēng)扇級控制體和一個風(fēng)扇排氣管道控制體。

4.1 均勻進(jìn)氣下的程序驗(yàn)證

為考核本文方法、計算程序的可靠性和合理性，首先在均勻進(jìn)氣時，將本模型計算結(jié)果與基于李亞普諾夫穩(wěn)定性理論的線化模型[4]結(jié)果和試驗(yàn)數(shù)據(jù)進(jìn)行比較，結(jié)果如圖2所示?？梢姡?0%～100%換算轉(zhuǎn)速范圍內(nèi)，兩個模型所得的特性線和穩(wěn)定邊界與試驗(yàn)結(jié)果都相當(dāng)接近，其中線化模型預(yù)測結(jié)果偏于保守。故可認(rèn)為，本研究所用的級特性計算方法和穩(wěn)定性判別方法可行。

圖2 計算模型的試驗(yàn)驗(yàn)證Fig.2 Model validation with experimental results

4.2 進(jìn)氣總壓畸變和總溫畸變下的分析

算例中研究畸變進(jìn)氣時將風(fēng)扇分為兩個相等的180°扇形子區(qū)，總壓畸變區(qū)為低壓區(qū)，總溫畸變區(qū)為高溫區(qū)。

圖3為進(jìn)口總壓畸變度和總溫畸變度都為3%時的喘振邊界比較?？梢?，相同畸變度下，總溫畸變比總壓畸變對壓縮系統(tǒng)的影響更大，穩(wěn)定裕度的減小也更明顯，這與文獻(xiàn)[1]、[5]的結(jié)果一致。

圖4示出了100%換算轉(zhuǎn)速下3%進(jìn)氣壓力畸變時風(fēng)扇內(nèi)總溫和總壓的非定常變化過程?？梢?，從0.10 s開始進(jìn)口給定總壓畸變，到0.15 s時流動重新趨于穩(wěn)定，顯然這一過程中葉柵通道內(nèi)生成了總溫畸變，并伴隨著總壓畸變在下游的大幅衰減。計算數(shù)據(jù)還表明，總壓畸變的高壓區(qū)是生成總溫畸變的低溫區(qū)，低壓區(qū)是高溫區(qū)。圖5為不同轉(zhuǎn)速下總溫生成和總壓衰減的快慢對比分析?？梢?，進(jìn)口總壓畸變時，風(fēng)扇轉(zhuǎn)速對畸變的影響不是很大。

圖3 壓力畸變和溫度畸變下的穩(wěn)定邊界比較Fig.3 Comparison of the stability limit between pressure and temperature distortion

圖4 進(jìn)口總壓畸變條件下的溫度和壓力非穩(wěn)態(tài)圖譜Fig.4 The unsteady process after inlet total pressure distortion

圖5不同發(fā)動機(jī)轉(zhuǎn)速條件下進(jìn)口總壓畸變產(chǎn)生的下游壓力和溫度畸變比較Fig.5 Comparison of the distortion index with inlet pressure distortion under different speeds

圖6示出了100%換算轉(zhuǎn)速下3%進(jìn)氣溫度畸變時風(fēng)扇內(nèi)總溫和總壓的非定常變化過程?？梢?，從0.10 s開始進(jìn)口給定總壓畸變，到0.15 s時流動重新趨于穩(wěn)定，很明顯這一過程中葉柵通道內(nèi)生成了總壓畸變，并伴隨著總溫畸變在下游的衰減。這里也要指出，總溫畸變的高溫區(qū)是生成總壓畸變的低壓區(qū)，低溫區(qū)是高壓區(qū)。圖7為不同轉(zhuǎn)速下總壓生成和總溫衰減的快慢對比分析?？梢姡M(jìn)口總溫畸變時風(fēng)扇轉(zhuǎn)速越高，總溫畸變衰減越快，而總壓畸變生成就越慢。

圖8示出了三種不同大小進(jìn)口壓力畸變所生成溫度畸變的比較?？梢?，總壓畸變越大，生成的總溫畸變越大；隨著總壓畸變的增大，所產(chǎn)生的總溫畸變也成比例增大，兩者間似乎存在定量關(guān)聯(lián)關(guān)系。

4.3 組合畸變下的分析

以對進(jìn)氣單獨(dú)總壓畸變和單獨(dú)總溫畸變的分析為基礎(chǔ)，研究進(jìn)口總壓總溫組合畸變對兩級風(fēng)扇的影響。圖9為反向組合畸變(壓力溫度畸變處于兩個不同的子區(qū)間)和正向組合畸變(壓力溫度畸變處于兩個相同的子區(qū)間)時穩(wěn)定裕度下降對比，其中總溫畸變和總壓畸變的強(qiáng)度都是3%。可見，反向組合畸變對系統(tǒng)穩(wěn)定性的影響很小，正向組合畸變穩(wěn)定裕度的損失比反向組合畸變的大得多，這與文獻(xiàn)[1]、[5]、[6]的結(jié)果相符。若與前文的進(jìn)口總壓和進(jìn)口總溫畸變一起比較，則反向組合畸變引起的損失最小，正向組合畸變引起的損失最大。這可理解為，由總壓畸變生成的總溫畸變，低壓區(qū)是高溫區(qū)，高壓區(qū)是低溫區(qū)；與之對應(yīng)，由總溫畸變引起的總壓畸變，高溫區(qū)為低壓區(qū)，低溫區(qū)為高壓區(qū)。如圖10所示，反向組合時由于這種生成畸變抵消的影響，壓力、溫度畸變都衰減得特別快，快于進(jìn)口單項(xiàng)總壓或總溫畸變；而正向組合時生成畸變有加強(qiáng)作用，所以兩種畸變都衰減得格外慢，慢于進(jìn)口單項(xiàng)總壓或總溫畸變。這就解釋了為什么反向組合時穩(wěn)定裕度損失最小、正向組合時損失最大。

圖6 進(jìn)口總溫畸變條件下的溫度和壓力非穩(wěn)態(tài)圖譜Fig.6 The unsteady process after inlet total temperature distortion

圖7 不同發(fā)動機(jī)轉(zhuǎn)速條件下進(jìn)口總溫畸變產(chǎn)生的下游壓力和溫度畸變比較Fig.7 Comparison of the distortion index with inlet temperature distortion under different speeds

圖8 不同進(jìn)口總壓畸變產(chǎn)生的下游溫度畸變Fig.8 Temperature distortion produced by different inlet pressure distortion intensity

圖9 反向總壓總溫畸變和正向總壓總溫畸變下的穩(wěn)定邊界比較Fig.9 Comparison of the stability limit between opposing and overlap combined distortion

圖10 反向總壓總溫畸變和重疊總壓總溫畸變產(chǎn)生的下游壓力和溫度畸變Fig.10 Comparison of the distortion index between opposing and overlap combined distortion

5 結(jié)論

(1)通過與試驗(yàn)數(shù)據(jù)和基于李亞普諾夫理論線化模型的比較，驗(yàn)證了所建立的用于分析多級軸流壓氣機(jī)瞬態(tài)性能及穩(wěn)定性的逐級預(yù)測模型的可靠性。

(2)進(jìn)口總壓畸變時，在氣流流動過程中會生成總溫畸變，且生成畸變的高溫區(qū)在總壓畸變的低壓區(qū)，低溫區(qū)在總壓畸變的高壓區(qū)；與之對應(yīng)，進(jìn)口總溫畸變同樣會生成總壓畸變，且生成畸變的高壓區(qū)在總溫畸變的低溫區(qū)，低壓區(qū)在總溫畸變的高溫區(qū)。

(3)進(jìn)氣總壓畸變條件下，風(fēng)扇轉(zhuǎn)速對畸變的生成和衰減影響不是很大，隨著總壓畸變強(qiáng)度的增大，生成的總溫畸變成比例增大；進(jìn)口總溫畸變條件下，轉(zhuǎn)速的影響更明顯，轉(zhuǎn)速越高，生成壓力畸變越慢，溫度畸變衰減就越快。

(4)反向組合畸變時由于生成畸變抵消的影響，壓力、溫度畸變都衰減得非?？?；與之相反，正向組合畸變時生成畸變有加強(qiáng)作用，所以兩種畸變都衰減得特別慢。這解釋了反向組合時穩(wěn)定裕度損失最小、正向組合時損失最大的原因。

(5)研究畸變在壓縮系統(tǒng)內(nèi)部的發(fā)展變化，有助于理解和預(yù)測各種畸變進(jìn)氣條件下穩(wěn)定裕度的損失。

[1]Hale A A，Davis Jr M W.DYNamic Turbine Engine Com?pressor Code:DYNTECC-TheoryandCapabilities[R].AIAA 92-3190，1992.

[2]Gorrell S E，Davis Jr M W.Application of a Dynamic Com?pression System Model to a Low-Aspect Ratio Fan:Casing Treatment and Distortion[R].AIAA 93-1871，1993.

[3]Shahrokhi K A，Davis Jr M W.Application of a Modified Dynamic Compression System Model to a Low-Aspect Ra?tio Fan:Effects and Inlet Distortion[R].AIAA 95-0301，1995.

[4]喬渭陽，蔡元虎，陳玉春，等.進(jìn)氣畸變對壓縮系統(tǒng)穩(wěn)定性影響的數(shù)值模擬[J].推進(jìn)技術(shù)，2001，22(4)：307—310.

[5]吳虎，廉小純，陳輔群，等.畸變進(jìn)氣下壓縮系統(tǒng)穩(wěn)定性分析的通用方法[J].推進(jìn)技術(shù)，1997，18(5)：62—64.

[6]董金鐘，王苗苗.總溫總壓組合畸變對壓氣機(jī)穩(wěn)定性影響的數(shù)值模擬研究[J].航空動力學(xué)報，2003，18(6)：788—793.

[7]Steinke R J.A Computer Code for Predicting Multi-Stage Axial-Flow Compressor Performance by a Meanline Stage Stacking Method[R].NASA TP-2020，1982.

[8]劉大響，葉培梁，胡駿，等.航空燃?xì)鉁u扇發(fā)動機(jī)穩(wěn)定性設(shè)計與評定技術(shù)[M].北京：航空工業(yè)出版社，2004.

[9]Urasek D C，Gorrell E T，Cunnan W S.Performance of Two-Stage Fan Having Low-Aspect Ratio First-Stage Ro?tor Blading[R].NASA TP-1493，1979.

Numerical Simulation of Influence by Inlet Distortion on Aerodynamic Stability in a Two Stage Fan

WANG Chun-li1，QIAO Wei-yang2
(1.Jiangxi Hongdu Aviation Industry Group Corporation Limited，Nanchang 330024，China；2.School of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China)

A quasi one-dimensional,stage-by-stage axial compression system mathematical model for sim?ulating surge line with clean or distortion inlet flow has been constructed based on compressor stage charac?teristics.The surge margin of a two stage fan was predicted by simulation according to the dynamic compres?sion system model.Compared with the experimental data and the linear model based on Liapunov’s theory, it is apparent that this model can predict the axial flow compressor’s stability margin correctly.According to the detailed numerical results for the stability of the two stage fan,inlet total pressure distortion will lead to the generation of total temperature distortion and decrease in the flowage,inlet total temperature distor?tion will cause total pressure distortion generation and reduce at the same time.Further more,inlet oppos?ing combined distortion of pressure and temperature cause drastic distortion decrease and moderate stabili?ty margin loss,while inlet overlaped combined distortion of pressure and temperature cause comparatively moderate distortion decrease and drastic stability margin loss.

aero-engine；axial flow compressor；stage-by-stage model；distortion；aerodynamic stability；numerical simulation

V231.3

：A

：1672-2620(2014)04-0001-06

2013-11-29；

：2014-07-28

王春利(1982-)，男，湖北紅安人，工程師，碩士，研究方向?yàn)橥七M(jìn)系統(tǒng)氣動熱力學(xué)。