白　濤

(西安航空學(xué)院 飛行器學(xué)院， 西安　710077)

?

不同加載形式高負(fù)荷低壓渦輪的性能分析

白濤

(西安航空學(xué)院 飛行器學(xué)院， 西安710077)

摘要：通過設(shè)計負(fù)荷能力相同而負(fù)荷分布形式不同的3種葉型來分析在寬廣的雷諾數(shù)范圍內(nèi)負(fù)荷分布對渦輪葉型邊界層發(fā)展的影響規(guī)律。研究結(jié)果表明：后加載葉型對雷諾數(shù)的變化較為敏感，尤其是在低雷諾數(shù)范圍內(nèi)，而前加載和均勻加載則相對不敏感。研究結(jié)果可為設(shè)計適用于不同工況的葉型提供理論指導(dǎo)。

關(guān)鍵詞：負(fù)荷分布；雷諾數(shù)；邊界層；氣動損失

低壓渦輪的質(zhì)量約占到航空發(fā)動機(jī)總質(zhì)量的20%～30%。高負(fù)荷低壓渦輪葉片的設(shè)計是減輕發(fā)動機(jī)質(zhì)量的主要途徑。高負(fù)荷葉片的設(shè)計意味著單級渦輪葉片數(shù)目減少、葉柵稠度降低。典型低壓渦輪工作的Re范圍較低，大約為 50 000～500 000。雷諾數(shù)對低壓渦輪的影響十分顯著[1-4]。尤其當(dāng)?shù)蛪簻u輪工作在低于自?；字Z數(shù)條件下時，雷諾數(shù)的降低會造成邊界層分離和端壁二次流的加劇，從而使得渦輪的效率降低[5]。因此，對高負(fù)荷低壓渦輪工作在低雷諾數(shù)情況下的流動控制變得十分必要。射流、拌線、附面層吹氣等流動控制方法通過加入能量或者改變邊界層的發(fā)展來實現(xiàn)對高負(fù)荷低壓渦輪流動的控制[6-7]，然而這些控制方法會增加葉片和控制機(jī)構(gòu)設(shè)計的難度。

葉片表面負(fù)荷分布決定著邊界層的發(fā)展規(guī)律，因此負(fù)荷分布形式對于渦輪氣動性能的好壞至關(guān)重要。高負(fù)荷葉片的加載形式一般有前加載、均勻加載和后加載3種，不同加載形式下邊界層的發(fā)展不同。Barry等[8]通過數(shù)值模擬和實驗的方法研究了低壓渦輪葉片幾何形狀在跨音流動條件下對渦輪氣動損失的影響，通過保證負(fù)荷能力大小相等、調(diào)整負(fù)荷分布形式得到3種不同加載形式的葉型。研究結(jié)果表明：后加載葉型使得邊界層發(fā)生分離，因此氣動損失明顯增大。Guilot等[9]通過改變?nèi)~片數(shù)、安裝角、尾緣彎折角來改變?nèi)~片的負(fù)荷程度和負(fù)荷分布。結(jié)果表明：負(fù)荷的增大、載荷向后移動以及Re的降低都會使吸力面邊界層增強(qiáng)和葉型氣動損失增加。相比于前加載葉型，后加載葉型的葉型氣動損失更小，但是前加載葉型會使得二次流損失更大。北京航空航天大學(xué)的鄒正平等[11-12]也在此方面做了大量的工作，并且通過設(shè)計后加載葉型來實現(xiàn)減小低壓渦輪氣動損失的目的。但目前在寬廣Re范圍內(nèi)，系統(tǒng)分析不同加載形式對渦輪葉片邊界層發(fā)展規(guī)律的影響的研究開展得較少。

本文首先設(shè)計負(fù)荷能力相同而加載形式不同的3種高負(fù)荷葉型，在寬廣的Re范圍下，分析不同加載形式對Re的敏感性，得出了可供參考的低壓渦輪設(shè)計理論數(shù)據(jù)。

1研究方法及研究對象

設(shè)計了3種加載形式的葉型。在設(shè)計過程中保證3種加載形式的載荷大小是相同的。載荷大小用Zweifel載荷系數(shù)來進(jìn)行衡量。本文主要通過改變安裝角和尾緣彎折角來實現(xiàn)不同的加載形式，其中：在均勻加載葉型的基礎(chǔ)上，通過增大尾緣彎折角、減小安裝角來得到后加載葉型；相反，通過減小尾緣彎折角、增大安裝角得到前加載葉型。在造型過程中，對前尾緣楔形角也進(jìn)行了適當(dāng)?shù)恼{(diào)整。由于高負(fù)荷低壓渦輪的需求，3種葉型的Zweifel載荷系數(shù)均為1.25。葉型設(shè)計的幾何參數(shù)如表1所示，葉型幾何截面如圖1所示。

表1　葉型幾何參數(shù)

圖1　 葉型示意圖

采用商用軟件CFX13.0求解三維定常黏性雷諾平均N-S方程，數(shù)值方法采用時間追趕的有限體積法，空間離散采用二階迎風(fēng)格式，時間離散應(yīng)用二階后差歐拉格式。選用SST湍流模型和γ-θ轉(zhuǎn)捩模型。數(shù)值模擬單層網(wǎng)格數(shù)取為12萬，近壁處的Y+均小于1，近壁處的延展比在1.2左右。計算邊界條件為：進(jìn)口給定總溫、總壓、氣流角；出口給定背壓。

2結(jié)果分析

圖2～5為不同Re下不同加載形式葉片表面的負(fù)荷分布。3種葉型的出口馬赫數(shù)均為0.6。前加載葉型在前緣后迅速加速到吸力峰，然后開始減速，直到尾緣附近，即前加載葉型的逆壓力梯度區(qū)域很長；而均勻加載葉型則是氣流緩慢加速，然后緩慢減速，即葉身前后部分的負(fù)荷是近似相同的；后加載葉型則是從前緣后開始加速，一直加速到葉身后半部分，然后減速。相比前加載葉型，后加載葉型的擴(kuò)壓段較短，但是其擴(kuò)散因子較大。當(dāng)前加載出口Re=2.5×104時(即低雷諾數(shù)狀態(tài))黏性力作用增強(qiáng)，因此對于擴(kuò)壓段而言，邊界層很快變厚，從而使得邊界層很容易發(fā)生分離。由于本文的葉型均為高負(fù)荷，后加載葉型的設(shè)計使得葉型后半段的負(fù)荷更大，擴(kuò)壓段的強(qiáng)逆壓力梯度使得葉型大概在65%軸向位置處發(fā)生分離，邊界層直到尾緣處沒有發(fā)生再附，但是在負(fù)荷分布上并沒有出現(xiàn)Mayle理論的開放式壓力平臺，這可能是由于分離泡中已經(jīng)出現(xiàn)了回流或者是后加載葉型在低Re下的適用性還有待研究。而前加載葉型雖然其擴(kuò)壓段較長，但擴(kuò)散因子較低，因而即使在低Re下，邊界層的分離泡也保持為很小。均勻加載葉型的變化則介于二者中間，邊界層在靠近尾緣位置處發(fā)生了分離，分離點相對于后加載葉型向下游方向移動，并且分離泡的長度也變短。

隨著Re的增大，滯止點后氣流的加速程度增大，因此使得葉型的負(fù)荷向前移動。由吸力峰的位置和峰值點可以看出：當(dāng)Re從2.5×104增大到5×104時，前加載葉型的吸力峰位置從25.3%移動到24.6%處，吸力峰的峰值點也略有所增加，因此隨著雷諾數(shù)的增大，前加載葉型在吸力峰后發(fā)生分離，但是由于分離泡和高Re的作用，邊界層可能發(fā)生了轉(zhuǎn)捩，因此分離泡發(fā)生了再附。對于均勻加載和后加載葉型，隨著Re的增大，葉身后半段的分離泡長度和分離高度都有所減弱。由于葉型高負(fù)荷的設(shè)計，使得在所研究的雷諾數(shù)范圍內(nèi)邊界層都有分離情況發(fā)生。

圖2　3種葉型葉片表面負(fù)荷分布(Re=2.5×104)

圖3　3種葉型葉片表面負(fù)荷分布(Re=5×104)

圖4　3種葉型葉片表面負(fù)荷分布(Re=10×104)

圖5　3種葉型葉片表面負(fù)荷分布(Re=15×104)

圖6　3種葉型葉片表面負(fù)荷分布(Re=20×104)

圖7～9為不同加載形式葉片吸力面形狀因子分布和動量厚度分布。形狀因子的變化可以反映邊界層的發(fā)展情況，而動量厚度則可以表征邊界層的氣動損失。下面通過這2個參數(shù)的分析來得出不同加載形式在不同Re下邊界層的變化特征。

首先分析Re=2.5×104時不同葉型邊界層的發(fā)展情況。在前緣滯止點后氣流劇烈加速，以及前緣吸力峰后的逆壓力梯度使得形狀因子迅速增大，但是依然沒有達(dá)到層流的分離值。

圖7　3種葉型吸力面形狀因子分布和

圖8　不同葉型吸力面形狀因子分布和

圖9　3種葉型吸力面形狀因子分布和

由于3種葉型前緣楔形角的差異導(dǎo)致前緣附近形狀因子也呈現(xiàn)差異，其中后加載葉型前緣楔形角較大，因此前緣附近邊界層的變化較均勻加載和前加載葉型都較為平緩。在葉身順壓梯度的作用下，邊界層形狀因子落回到2.5左右。氣流在葉身處分別經(jīng)歷了加速段和擴(kuò)壓段的減速過程。在葉身處3種葉型的邊界層發(fā)展呈現(xiàn)出了很大的差異，其中前加載葉型使邊界層經(jīng)過較短的順壓力梯度后迅速發(fā)展到逆壓力梯度，因此形狀因子開始增大的位置也很靠近前緣，形狀因子的峰值點出現(xiàn)在葉背弧長的40%左右，最大值為4.5，此時邊界層發(fā)生了分離，但是分離泡較弱(分離泡很小，在流線圖上沒有示意出來)。對于后加載葉型，邊界層在前緣后首先在較長的順壓力梯度下發(fā)展，因此邊界層發(fā)展較前加載葉型更為緩慢，直到逆壓力梯度的出現(xiàn)，邊界層發(fā)展才迅速加快，其厚度急劇增大。由后加載葉型形狀因子的分布可以看出：形狀因子開始增大的地方相比前加載葉型更加靠后；在70%弧長位置處形狀因子達(dá)到最大，即此處邊界層的分離高度達(dá)到最大，邊界層的分離大概發(fā)生在60%弧長位置處。這是因為在分離泡中發(fā)生了分離流的轉(zhuǎn)捩，轉(zhuǎn)捩區(qū)的長度較短，而轉(zhuǎn)捩后由于邊界層逆壓力梯度較大，因此邊界層始終沒有發(fā)生再附，也可以認(rèn)為湍流邊界層一直沒有發(fā)生重建。均勻加載葉型形狀因子的變化則介于前加載和后加載葉型之間，在80%弧長處邊界層的值達(dá)到最大，但峰值點小于后加載葉型，一直到尾緣處，形狀因子的值都保持在5左右，此時的邊界層為開放式分離，分離的區(qū)域已經(jīng)達(dá)到很大，邊界層的外緣不容易判斷，可能使得形狀因子的計算不夠準(zhǔn)確。由形狀因子的大小可以看出：邊界層亦發(fā)生分離，但分離的長度和強(qiáng)度均小于后加載葉型，同時邊界層在分離流中發(fā)生了轉(zhuǎn)捩。前加載葉型的轉(zhuǎn)捩長度更長，而均勻加載葉型轉(zhuǎn)捩完成后邊界層直到尾緣處都沒有來得及再附。

由動量厚度的大小可以得出類似對應(yīng)的規(guī)律：在60%弧長以前，3種葉型的邊界層均沒有發(fā)生分離。動量厚度是邊界層氣動損失大小的衡量，而對于附著的邊界層，氣動損失的大小同邊界層外緣速度即主流速度呈正比。在前加載葉型前面部分劇烈的加速使得邊界層的損失大于后加載和均勻加載葉型，但是在60%弧長后邊界層發(fā)生分離，分離使得邊界層的氣動損失迅速增大。由圖7可以看出：后加載葉型的損失迅速增大，增大的趨勢遠(yuǎn)大于前加載和均勻加載葉型，而前加載葉型保持為最小，均勻加載葉型居中。值得注意的是：在30%～60%弧長范圍內(nèi)，前加載葉型的氣動損失都要高于前加載和均勻加載葉型，這是由于前加載葉型提前轉(zhuǎn)捩造成的。

通過對動量厚度分布的分析可以得出：無論是在低雷諾數(shù)還是在高雷諾數(shù)的情況下，在前緣后的順壓力梯度區(qū)域，由于邊界的發(fā)展較為緩慢，因此3種葉型的動量厚度差異很小。但是前加載葉型的設(shè)計使得邊界層更易于提前發(fā)生轉(zhuǎn)捩，因此在順壓梯度后，前加載葉型的動量厚度均保持最大，隨著雷諾數(shù)的增大，轉(zhuǎn)捩的位置進(jìn)一步提前，因此前加載葉型動量厚度分布較大的區(qū)域增大。當(dāng)雷諾數(shù)為2.5×104時，只有在發(fā)生分離的區(qū)域內(nèi)，均勻加載葉型吸力面的氣動損失才會大于前加載葉型，而隨著雷諾數(shù)增大到10×104時，均勻加載葉型的分離損失進(jìn)一步減少，而前加載葉型由于更加提前的轉(zhuǎn)捩造成的湍流損失則增大，因此在這個吸力面的邊界層內(nèi)，前加載葉型的氣動損失都大于均勻加載葉型。對于后加載葉型，其在葉身前半部分的氣動損失都保持為最小，而在擴(kuò)壓段內(nèi)，分離使得邊界層內(nèi)的氣動損失急劇增大，因此其損失在3種葉型中保持為最大。由形狀因子的分布還可以得出：前緣對雷諾數(shù)也體現(xiàn)出了敏感性，當(dāng)雷諾數(shù)為10×104時，均勻加載和前加載葉型前緣附近均出現(xiàn)了分離泡。

由于在低雷諾數(shù)下，后加載葉型擴(kuò)壓段較大的分離區(qū)域使得邊界層的損失非常大，因此在Re=2.5×104時，后加載葉型的氣動損失很大，遠(yuǎn)遠(yuǎn)大于均勻加載和前加載葉型。而隨著Re的增大，邊界層黏性作用減小，因此3種葉型的氣動損失均有所降低。對于后加載葉型，由于其分離強(qiáng)度隨著雷諾數(shù)的增大變化十分敏感，其擴(kuò)壓段由開式的強(qiáng)分離泡變成閉式的短分離泡，因此對于后加載葉型而言，其不同Re氣動損失的差異主要是由于層流的分離程度造成的，體現(xiàn)出了對雷諾數(shù)的敏感性。而對于前加載葉型，隨著Re的增大，其氣動損失亦降低，但其降低的程度較后加載和均勻加載葉型均較為緩慢，這是因為雖然Re的增大使得整個邊界層內(nèi)的速度分布變得更加飽滿，因此邊界層內(nèi)整體的氣動損失會降低，但是隨著雷諾數(shù)增大，前加載葉型會使得邊界層轉(zhuǎn)捩的位置提前，從而使得邊界層的氣動損失增大，當(dāng)Re達(dá)到很高時，這種損失會和分離泡(隨著Re的增大，分離泡再減弱)帶來的氣動損失相當(dāng)。從損失的分布可以看出：當(dāng)Re>5×104時，前加載葉型的氣動損失均大于均勻加載葉型，而當(dāng)Re達(dá)到20×104時，前加載葉型的氣動損失和后加載葉型的氣動損失相當(dāng)。

圖10　3種葉型葉柵通道中的氣動損失系數(shù)

3結(jié)論

通過分析3種不同加載形式葉型在不同雷諾數(shù)下邊界層發(fā)展規(guī)律的影響得出：就葉型氣動損失而言，后加載葉型對雷諾數(shù)的變化非常敏感，其在整個雷諾數(shù)范圍內(nèi)的氣動損失都較大；而隨著雷諾數(shù)的增大，前加載葉型的高氣動損失特性開始顯現(xiàn)出來。因此，在低雷諾數(shù)下可以選擇前加載和均勻加載葉型來降低葉型氣動損失，而在高雷諾數(shù)下，前加載葉型則不是很好的選擇。

參考文獻(xiàn)：

[1]鄒正平,葉建,劉火星,等.低壓渦輪內(nèi)部流動及其氣動設(shè)計研究進(jìn)展[J].力學(xué)進(jìn)展,2007,37(4):551-562.

[2]鄒正平,周琨,王鵬,等.大涵道比渦扇發(fā)動機(jī)渦輪內(nèi)部流動機(jī)理及氣動設(shè)計技術(shù)研究進(jìn)展[J].航空制造技術(shù),2012 (13):49-54.

[3]DENTONJD.Lossmechanismsinturbomachines[J].ASMEJournalofturbomachinery.1993,115(4):621-656.

[4]HODSONHP,HOWELLRJ.Theroleoftransitioninhigh-liftlow-pressureturbinesforaeroengine[J].progressinaerospacescience,2005,41:419-514.

[5]楊琳,馮濤,鄒正平.低雷諾數(shù)渦輪內(nèi)部流場分析[J].北京航空航天大學(xué)學(xué)報,2006,31(11):1195-1197.

[6]CHENPP,QIAOWY,HOUWT,etal.Studyonhigh-loadingturbineairfoildesignbasedonflowcontrol[J].Journalofaerospacepower,2010,25(11):2528-2538.

[7]ZHANGXF,VeraHODSONM.Separationandtransitioncontrolonanaft-loadedultra-high-liftLPturbinebladeatlowReynoldsnumber:low-speedinvestigation[J].JTurbomach,2005,128(3)：517-527.

[8]BARRYJ,BROWN,MALANDRAA.Effectofturbineairfoilshapeonaerodynamiclossesforturbineairfoilsoperatingundertransoniccondition[C]//proceedingofASMEturbineExpro.USA:[s.n.],2011.

[9]GUILOTS,NGWF.EffectofloadinglevelanddistributiononLPTlosses.[C]//Turbomachinery:AxialFlowTurbineAerodynamics.USA:[s.n.],2008.

[10]ZORICT,POPOVICT,SJOLANDERSA.ComparativeinvestigationofthreehighlyloadedLPturbineairfoils:Part2-Measuredprofileandsecondarylossesatoff-designincidence[C]//ASMETurboExpo2007.USA:[s.n.],2007.

[11]馮濤,程洪貴,楊琳,等.邊界層特性對雷諾數(shù)變化的敏感性分析[J].推進(jìn)技術(shù),2005,26(4):329-334.

[12]楊琳,劉火星,鄒正平.低雷諾數(shù)條件下低壓渦輪氣動設(shè)計[J].工程熱物理學(xué)報,2005,26(2):229-230.

(責(zé)任編輯劉舸)

Performance Analysis of High Loaded Low PressureTurbinewithDifferentLoadingDistribution

BAI Tao

(School of Aircraft, Xi’an Aeronautical University, Xi’an 710077, China)

Abstract:TheeffectlawofloadingdistributiononboundarylayerdevelopmentatwiderangeReynoldsnumberwasstudiedthroughdesignthreeturbinebladeswithdifferentloaddistributionbutsameloadlevel.Theresearchindicatedthatafter-loadedairfoilismoresensitivetoReynoldsnumberespeciallywhenatthelowReynoldsnumberswhilethefrontandmidloadedairfoilislesssensitive.Theresearchcanprovidetheoryguidanceforthedesignofairfoilwhichworkedatdifferentworkcondition.

Keywords:loaddistribution;Reynoldsnumber;boundarylayer;aerodynamicloss

收稿日期：2016-01-18

基金項目：國家自然科學(xué)基金青年基金資助項目(51406003)

作者簡介：白濤(1988—)，女，陜西榆林 人，碩士研究生，主要從事動力工程及工程熱物理研究。

doi:10.3969/j.issn.1674-8425(z).2016.06.011

中圖分類號：V231.3

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

文章編號：1674-8425(2016)06-0064-06

引用格式：白濤.不同加載形式高負(fù)荷低壓渦輪的性能分析[J].重慶理工大學(xué)學(xué)報(自然科學(xué))，2016(6):64-69.

Citationformat：BAITao.PerformanceAnalysisofHighLoadedLowPressureTurbinewithDifferentLoadingDistribution[J].JournalofChongqingUniversityofTechnology(NaturalScience)，2016(6):64-69.