張少博,史偉,劉火星
(1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院航空發(fā)動(dòng)機(jī)氣動(dòng)熱力國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100191;2.先進(jìn)航空發(fā)動(dòng)機(jī)協(xié)同創(chuàng)新中心,北京100191)
脈沖式噴氣調(diào)節(jié)渦輪流量的數(shù)值研究
張少博1,史偉1,劉火星2
(1.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院航空發(fā)動(dòng)機(jī)氣動(dòng)熱力國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100191;2.先進(jìn)航空發(fā)動(dòng)機(jī)協(xié)同創(chuàng)新中心,北京100191)
采用三維數(shù)值模擬方法,對(duì)在導(dǎo)葉表面射入脈沖式氣流的渦輪流場(chǎng)進(jìn)行計(jì)算,對(duì)比了脈沖式射流與靜態(tài)射流對(duì)渦輪流量的影響,分析了脈沖頻率、脈沖幅值對(duì)渦輪流量調(diào)節(jié)的影響規(guī)律。結(jié)果表明:在保證與靜態(tài)射流同等工況調(diào)節(jié)效果的同時(shí),脈沖式射流對(duì)渦輪的效率有所提升,低頻率的脈沖噴氣有更寬廣的渦輪流量調(diào)節(jié)范圍,較高幅值的脈沖噴氣更有利于渦輪流量調(diào)節(jié)。
變循環(huán)發(fā)動(dòng)機(jī);渦輪;流量調(diào)節(jié);脈沖頻率;脈沖幅值;三維數(shù)值模擬
超聲速運(yùn)輸和多任務(wù)戰(zhàn)斗機(jī)對(duì)發(fā)動(dòng)機(jī)提出了兩方面的要求:兼具高馬赫數(shù)飛行狀態(tài)和超聲速巡航,以提高飛機(jī)生存能力;低油耗的經(jīng)濟(jì)工作狀態(tài)以增強(qiáng)飛機(jī)遠(yuǎn)程奔襲能力。變循環(huán)發(fā)動(dòng)機(jī)可改變自身的空氣流量和單位推力以適應(yīng)不同的飛行狀態(tài),在整個(gè)飛行包線內(nèi)都能獲得較好的性能,是未來(lái)發(fā)動(dòng)機(jī)的發(fā)展趨勢(shì)和方向[1-3]。為適應(yīng)發(fā)動(dòng)機(jī)工作過(guò)程中循環(huán)模式的改變,需要采取調(diào)節(jié)手段來(lái)調(diào)整渦輪部件的工作狀態(tài)。
渦輪工況調(diào)節(jié)技術(shù)主要包括幾何調(diào)節(jié)技術(shù)和氣動(dòng)調(diào)節(jié)技術(shù),但其實(shí)質(zhì)都是通過(guò)改變渦輪喉道的流通面積,來(lái)改變渦輪的流通能力。幾何調(diào)節(jié)由于結(jié)構(gòu)復(fù)雜及額外增重等弊端,降低了發(fā)動(dòng)機(jī)的可靠性和推重比[4]。氣動(dòng)調(diào)節(jié)則因其具有結(jié)構(gòu)簡(jiǎn)單、噴氣可以利用渦輪導(dǎo)葉內(nèi)部原有的冷卻結(jié)構(gòu)和空氣系統(tǒng)等優(yōu)點(diǎn)[5],近年來(lái)越來(lái)越受到人們的重視。北京航空航天大學(xué)的楊旸[6]和付超[7]研究了葉片表面噴氣對(duì)渦輪流量的調(diào)節(jié),閆晨等[8]研究則是利用端區(qū)射流來(lái)調(diào)節(jié)渦輪流量,但以上研究的都是靜態(tài)噴氣對(duì)渦輪流量的影響。Saracoglu等[9]采用能量分析方法研究了高壓渦輪導(dǎo)葉尾緣的脈沖式冷氣射流,發(fā)現(xiàn)相對(duì)于不噴氣狀態(tài),脈沖式射流使得尾緣損失顯著降低。因此,借鑒尾緣脈沖式冷氣射流的研究方案,對(duì)渦輪流量的氣動(dòng)調(diào)節(jié)技術(shù)進(jìn)行深入研究,為未來(lái)變循環(huán)發(fā)動(dòng)機(jī)的設(shè)計(jì)提供理論支持,具有重要的工程應(yīng)用意義。
基于上述認(rèn)識(shí),本文采用三維數(shù)值模擬方法,研究葉片表面脈沖式噴氣對(duì)渦輪性能的影響,探討其對(duì)渦輪流量的調(diào)節(jié)規(guī)律。
2.1研究對(duì)象
根據(jù)總體設(shè)計(jì)要求和研究需求,選取某1+1對(duì)轉(zhuǎn)渦輪的低壓渦輪作為研究對(duì)象。通過(guò)在葉片表面開槽的方式,向渦輪氣流通道內(nèi)部引入一股高速氣流,與主流相互作用。為將取自高壓壓氣機(jī)的氣流以很高的速度引入射流通道,選定拉法爾噴管作為噴氣的氣路裝置,如圖1所示。經(jīng)優(yōu)化研究后,選定噴氣位置為94.38%軸向弦長(zhǎng)處,噴氣方向?yàn)榇怪比~片表面。
圖1 帶超聲速噴氣通道的葉型Fig.1 Airfoil with supersonic airflow passage
2.2計(jì)算域和網(wǎng)格
數(shù)值模擬采用商業(yè)軟件CFX12.1進(jìn)行,模擬方法為時(shí)間追趕的有限體積法,并利用多重網(wǎng)格技術(shù)加速收斂??臻g離散采用二階迎風(fēng)格式,時(shí)間離散采用二階后差歐拉格式。湍流模型采用SST模型。計(jì)算網(wǎng)格如圖2、圖3所示,高壓動(dòng)葉葉尖間隙內(nèi)部給定9個(gè)網(wǎng)格節(jié)點(diǎn)。環(huán)繞葉片表面均使用O型網(wǎng)格,導(dǎo)葉進(jìn)口段和動(dòng)葉出口段均增加一個(gè)H型網(wǎng)格與之相連。對(duì)于帶有噴氣通道的導(dǎo)葉,參照DENTON程序處理網(wǎng)格的思路,將網(wǎng)格分成多個(gè)部分進(jìn)行拓?fù)?,且?dǎo)葉的壓力面和吸力面不再作為一個(gè)整體,而是分開在兩個(gè)周期面上,這樣可避免在葉片周圍采用O-block;在前緣、尾緣及噴氣孔所在位置,將整體結(jié)構(gòu)分成三個(gè)部分,并簡(jiǎn)化每一部分的拓?fù)洹?/p>
圖2 低壓渦輪網(wǎng)格Fig.2 Low pressure turbine grid
圖3 超聲速噴氣通道網(wǎng)格放大圖Fig.3 Grid of supersonic airflow passage
數(shù)值模擬前,選擇了4套網(wǎng)格進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證,網(wǎng)格分布如表1所示,模擬結(jié)果如圖4所示。可見,隨著網(wǎng)格數(shù)的增加,計(jì)算的導(dǎo)葉進(jìn)口流量不斷減小。Case 3和Case 4計(jì)算的導(dǎo)葉進(jìn)口流量已趨于穩(wěn)定且比較接近,因此下文所有算例選擇Case 3所示網(wǎng)格數(shù)及網(wǎng)格節(jié)點(diǎn)分布進(jìn)行計(jì)算。
表1 網(wǎng)格數(shù)分布Table 1 Grid number distribution
2.3算例選取
Saracoglu研究脈沖式冷卻射流時(shí),選取了頻率為200 Hz、幅值為0.8倍平均值的正弦波形。基于此,本文選擇噴氣孔進(jìn)口流量時(shí)均值為0.214 kg/s的正弦式周期振蕩的脈沖m=[1+mrsin(2πft)]作為研究對(duì)象。其中,m為脈沖流量,為時(shí)均流量,mr為相對(duì)幅值,f為脈沖頻率,如圖5所示。選取噴氣孔位置附近作為監(jiān)測(cè)點(diǎn),選擇100 Hz、300 Hz、600 Hz三個(gè)頻率值探究頻率對(duì)渦輪流量的影響規(guī)律;為探討幅值的影響,分別選取0.05 kg/s、0.10 kg/s、0.15 kg/s三個(gè)幅值(相對(duì)幅值分別為23.36%、46.73%、70.09%)進(jìn)行討論。由于研究的重點(diǎn)是噴氣孔的脈沖射流效應(yīng),故在非定常計(jì)算中沒(méi)有使用Domain Scaling方法,主要是通過(guò)在射流孔進(jìn)口給定周期性脈動(dòng)來(lái)設(shè)置非定常特性。對(duì)于一個(gè)特定頻率的非定常計(jì)算,將一個(gè)周期分成200份作為時(shí)間步長(zhǎng),計(jì)算6個(gè)周期,對(duì)后面兩個(gè)周期的信息進(jìn)行統(tǒng)計(jì)處理。除此之外,非定常計(jì)算采用的數(shù)值方法與定常計(jì)算的數(shù)值方法一致。
圖5 頻率與幅值的分布關(guān)系Fig.5 Distribution relationship of frequency and amplitude
圖6 脈沖式射流收斂過(guò)程Fig.6 Convergence process of impulsive airflow injection
以300 Hz為例,圖6示出了非定常計(jì)算檢測(cè)點(diǎn)的收斂過(guò)程,圖中橫坐標(biāo)為無(wú)量綱化時(shí)間,縱坐標(biāo)為按照不噴氣時(shí)導(dǎo)葉進(jìn)口流量參數(shù)進(jìn)行的無(wú)量綱化處理??梢?,檢測(cè)到的導(dǎo)葉進(jìn)口流量和監(jiān)測(cè)點(diǎn)的靜壓、速度呈現(xiàn)出很好的周期性,說(shuō)明脈沖式噴氣計(jì)算結(jié)果收斂較好。另外,導(dǎo)葉進(jìn)口流量與噴氣孔出口檢測(cè)點(diǎn)速度的變化并不同步,其峰值點(diǎn)有一定相位差,表現(xiàn)出一定的滯后性,這主要是由于非定常計(jì)算里的容腔效應(yīng)所致。如果把渦輪通道看作一個(gè)大容腔,從葉片表面一條狹縫周期性地往容腔里噴氣,壓力可快速變化,而速度的變化會(huì)表現(xiàn)得較為緩慢。從噴氣孔出口檢測(cè)點(diǎn)速度變化看,一個(gè)周期內(nèi)70%的時(shí)間無(wú)量綱速度均為0.6,而另外30%的時(shí)間速度有一個(gè)急劇的增減,最高值甚至達(dá)到1.8。這可能是容腔導(dǎo)致的氣體累積,進(jìn)而出現(xiàn)速度躍升。另一方面,沿著葉片表面進(jìn)行三維脈沖噴氣,沿葉高方向會(huì)產(chǎn)生不均勻性,也會(huì)導(dǎo)致氣流累積。
表2給出了定常噴氣和三個(gè)不同頻率脈動(dòng)噴氣對(duì)渦輪流量的影響。以定常結(jié)果為基準(zhǔn)可以看出,與靜態(tài)噴氣相比,脈沖式噴氣渦輪導(dǎo)葉進(jìn)口流量及流量變化率略微減小,而噴管出口的馬赫數(shù)明顯降低,有利于降低導(dǎo)葉尾緣復(fù)雜激波系的強(qiáng)度,從而使尾緣處的損失減小,渦輪效率有所增加。因此,脈沖式射流在保證同等工況調(diào)節(jié)效果的同時(shí),還提高了渦輪效率。
表2 頻率變化時(shí)射流總參數(shù)對(duì)比Table 2 Comparison of injection airflow parameters with different frequency
通常,噴管出口馬赫數(shù)在一定范圍內(nèi)越大,射流深度越大,形成的低速摻混區(qū)就越大,其調(diào)節(jié)效果就越好。由上表可以看出,噴氣流量相同的情況下,較小的脈動(dòng)頻率對(duì)應(yīng)的噴管出口馬赫數(shù)較低。這表明要達(dá)到相同的流量調(diào)節(jié)效果,采用低頻率的脈沖式噴氣有利于減少冷氣的使用量,這就使進(jìn)一步增加噴氣流量有了更大的可能,進(jìn)而增加渦輪流量調(diào)節(jié)范圍。
為表征葉片負(fù)荷大小及分布,定義葉片表面壓力系數(shù)Cp=(pt1-ps)/(pt1-ps1),其中pt1代表葉柵進(jìn)口總壓,ps1表示葉柵進(jìn)口靜壓,ps代表葉片表面靜壓。圖7、圖8給出了定常噴氣和三個(gè)不同頻率脈動(dòng)噴氣時(shí)葉片表面的壓力系數(shù)分布,其中脈沖式噴氣時(shí)采用的都是時(shí)均參數(shù),Z表示軸向位置,Cax表示軸向弦長(zhǎng)。圖中表明,四種情況下葉片表面的壓力系數(shù)幾乎一致,這說(shuō)明脈沖式射流基本上不會(huì)引起葉片表面負(fù)荷變化,僅在射流孔位置附近表現(xiàn)出一定的差異性。
圖7 葉中截面壓力系數(shù)分布Fig.7 Pressure coefficient distribution of mid-span
圖8 葉片尾緣負(fù)荷放大圖Fig.8 Trailing edge loading distribution
圖9 不同時(shí)刻葉中截面馬赫數(shù)分布云圖Fig.9 Mach number distribution contours of mid-span at different time
圖9給出了100 Hz頻率下不同時(shí)刻葉中截面的馬赫數(shù)分布云圖。可見,隨著時(shí)間的推移,噴氣量的脈動(dòng)僅會(huì)影響尾緣后的低速區(qū)大小,而對(duì)其他部分的流場(chǎng)影響不大。
表3給出了三個(gè)不同幅值脈沖噴氣對(duì)渦輪流量的影響。以定常結(jié)果為基準(zhǔn)可以看出,隨著脈動(dòng)幅值的逐漸增大,流量變化率逐漸減小,對(duì)渦輪流量的控制有一定程度下降;噴管出口馬赫數(shù)明顯減小,利于渦輪流量調(diào)節(jié)范圍的增加;同時(shí),渦輪效率明顯提高。相對(duì)而言,與噴管出口馬赫數(shù)下降的幅度相比,流量變化率減小不明顯,因此較高的脈動(dòng)幅值對(duì)渦輪流量調(diào)節(jié)更有利。
表3 幅值變化時(shí)射流參數(shù)對(duì)比Table 3 Comparison of injection airflow parameters with different amplitude
圖10 不同幅值時(shí)葉中截面時(shí)均馬赫數(shù)云圖分布Fig.10 Mach number distribution contours of mid-span with different amplitude
圖10給出了不同幅值時(shí)葉中截面時(shí)均馬赫數(shù)云圖分布和聲速線位置,圖中黑線表示馬赫數(shù)為1的等值線,即喉道位置。對(duì)比幾幅圖可知,增加脈動(dòng)幅值可有效減少尾跡區(qū)域及摻混區(qū)的強(qiáng)度,而射流深度也有所減小。氣動(dòng)調(diào)節(jié)主要是通過(guò)射流堵塞渦輪氣流通道,脈沖式噴氣的幅值對(duì)喉道面積減少量影響很小,但由射流和主流形成的低速區(qū)卻隨著脈動(dòng)幅值的增加而減小。這說(shuō)明脈沖式射流在控制導(dǎo)葉進(jìn)口流量方面作用不大,高幅值的脈沖射流能使損失明顯降低。
(1)脈沖式噴氣與靜態(tài)噴氣相比,在保證同等工況調(diào)節(jié)效果的同時(shí),對(duì)渦輪的效率有所提升。
(2)低頻率的脈沖噴氣有更寬廣的渦輪流量調(diào)節(jié)范圍,同時(shí)渦輪的效率還有所提升;較高幅值的脈沖噴氣更有利于渦輪流量調(diào)節(jié)。
致謝:感謝北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院鄒正平教授,在本課題研究中給予的指導(dǎo)和幫助。
[1]胡軍.ADVENT,渦輪噴氣發(fā)動(dòng)機(jī)的新里程碑[R].北京:中國(guó)航空工業(yè)發(fā)展研究中心航空技術(shù)所,2008.
[2]方昌德.變循環(huán)發(fā)動(dòng)機(jī)[J].燃?xì)鉁u輪試驗(yàn)與研究,2004,17(3):1—5.
[3]劉大響,金捷.21世紀(jì)世界航空動(dòng)力技術(shù)發(fā)展趨勢(shì)與展望[J].中國(guó)工程科學(xué),2004,9(6):1—8.
[4]胡松巖.變幾何渦輪及其設(shè)計(jì)特點(diǎn)[J].航空發(fā)動(dòng)機(jī),1996,3(24):30—36.
[5]Latimer R J.Variable flow turbines[R].AGARD CP-205,1977.
[6]楊旸.氣動(dòng)方式調(diào)節(jié)渦輪流量[D].北京:北京航空航天大學(xué),2009.
[7]付超.高性能渦輪氣動(dòng)設(shè)計(jì)若干關(guān)鍵技術(shù)研究[D].北京:北京航空航天大學(xué),2010.
[8]閆晨,葉建,付超.利用端區(qū)射流調(diào)節(jié)渦輪流量的數(shù)值研究[J].燃?xì)鉁u輪試驗(yàn)與研究,2011,24(2):9—12.
[9]Saracoglu B H,Paniagua G,Salvadori S.Energy analysis of pulsating coolant ejection[R].ASME GT2014-25868,2014.
Numerical investigation of turbine mass-flow rate regulation using impulsive airflow injection
ZHANG Shao-bo1,SHI Wei1,LIU Huo-xing2
(1.National Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics,School of Energy and Power Engineering,Beijing University of Aeronautics and Astronautics,Beijing 100191,China;2.Collaborative Innovation Center of Advanced Aero-Engine,Beijing 100191,China)
The turbine flow field with guide vane surface using impulsive airflow injection was calculated with three-dimensional numerical simulation method.The influence of impulsive airflow injection on turbine mass-flow rate was compared with the influence of steady airflow injection,and the effects of impulsive frequency and impulsive amplitude on turbine mass-flow rate were analyzed.The results indicate that the impulsive airflow injection can improve turbine's efficiency when the modulation effects of impulsive airflow injection are the same as the effects of steady airflow injection;low-frequency impulsive airflow injection can bring a wider modulation range of turbine mass-flow rate;large-amplitude impulsive airflow injection is beneficial to the modulation of turbine mass-flow rate.
variable cycle engine;turbine;mass-flow rate modulation;impulsive frequency;impulsive amplitude;three-dimensional numerical simulation
V231.3
A
1672-2620(2015)06-0026-05
2015-03-18;
2015-06-23
張少博(1991-),男,河南安陽(yáng)人,碩士研究生,從事葉輪機(jī)內(nèi)部流動(dòng)機(jī)理研究。