潘炳華，任芳，郭文，陶智

(1.中國燃?xì)鉁u輪研究院，四川成都610500；2.北京航空航天大學(xué)，北京100191)

旋轉(zhuǎn)狀態(tài)下氣膜冷卻效率試驗研究

潘炳華1，任芳1，郭文1，陶智2

(1.中國燃?xì)鉁u輪研究院，四川成都610500；2.北京航空航天大學(xué)，北京100191)

以某型發(fā)動機(jī)高壓渦輪轉(zhuǎn)子葉片吸力面腮區(qū)氣膜孔為研究對象，通過模擬發(fā)動機(jī)狀態(tài)的模型試驗，研究了旋轉(zhuǎn)數(shù)、吹風(fēng)比和主流雷諾數(shù)對氣膜孔冷卻效率的影響。結(jié)果表明，旋轉(zhuǎn)會導(dǎo)致氣膜覆蓋區(qū)域向高半徑方向偏轉(zhuǎn)，且旋轉(zhuǎn)數(shù)越大，偏轉(zhuǎn)角度越大，氣膜冷卻效率越低；同時，旋轉(zhuǎn)會弱化吹風(fēng)比、主流雷諾數(shù)等對氣膜冷氣效率的影響。研究獲得的旋轉(zhuǎn)狀態(tài)下渦輪轉(zhuǎn)子葉片型面典型區(qū)域氣膜冷卻特性的試驗數(shù)據(jù)，可為發(fā)動機(jī)轉(zhuǎn)子葉片冷卻設(shè)計提供參考。

航空發(fā)動機(jī)；渦輪葉片；氣膜冷卻；吹風(fēng)比；旋轉(zhuǎn)數(shù)

1 引言

氣膜冷卻是高負(fù)荷燃?xì)鉁u輪上應(yīng)用廣泛的冷卻方式之一。通過在壁面附近以一定角度射入一股低溫氣流，將高溫燃?xì)馀c壁面隔離，達(dá)到對受熱壁面進(jìn)行冷卻、保護(hù)的作用。國外學(xué)者對此開展了廣泛、深入的研究。如Ito等[1]在不同葉柵上研究了葉片表面曲率對氣膜覆蓋效果的影響；Drost等[2]對壓力面和吸力面的氣膜效率進(jìn)行了對比研究，發(fā)現(xiàn)吸力面的氣膜孔效率高于壓力面；Gritsch等[3]研究了氣膜孔幾何形狀對氣膜效率和孔流量系數(shù)的影響。國內(nèi)高校和科研院所近年來也對氣膜冷卻葉片進(jìn)行了大量的數(shù)值研究。如曾軍等[4]采用源項模擬技術(shù)對氣膜冷卻導(dǎo)向葉片進(jìn)行了數(shù)值模擬，其結(jié)果與試驗數(shù)據(jù)吻合較好；朱惠人等[5]在低雷諾數(shù)k-ε雙方程紊流模型程序的基礎(chǔ)上，引入一種離散孔氣膜冷卻的噴射模型，使該程序能較好地預(yù)測渦輪葉片上有冷氣出流情況下的冷卻效率；吳宏等[6]數(shù)值模擬了旋轉(zhuǎn)狀態(tài)下曲率表面的氣膜冷卻效率，表明不同曲率表面的氣膜冷卻效率不同。但國內(nèi)學(xué)者對氣膜冷卻的試驗研究相對較少，且試驗方案基本局限于靜止?fàn)顟B(tài)，如向安定等[7]進(jìn)行的渦輪葉片表面氣膜冷卻傳熱試驗研究。本文通過模擬發(fā)動機(jī)狀態(tài)的模型試驗，研究了旋轉(zhuǎn)狀態(tài)下旋轉(zhuǎn)數(shù)、吹風(fēng)比和主流雷諾數(shù)對氣膜孔冷卻效率的影響。

2 試驗?zāi)Ｐ?/h2>

試驗?zāi)Ｐ蜑閺澢膱A弧表面，曲率半徑r=200 mm。其彎曲表面并非真實形狀的渦輪葉片，而是采用圓弧段來模擬不同形式的曲面，如圖1所示。試驗件的主流通道部分為縱截面尺寸固定(60 mm×60 mm)的方形通道，其中包含試驗表面的部分為圓弧形，中心角為36°，對應(yīng)測試面圓弧長度約為125 mm，圓弧段兩端分別為進(jìn)、出口直通道。氣膜孔位于圓弧表面的正中間位置，為標(biāo)準(zhǔn)圓柱形，其直徑D為4 mm，流向傾斜角α為30°，側(cè)向傾斜角?為0°，展向傾斜角β為90°。氣膜孔中心到旋轉(zhuǎn)軸心的距離為450 mm。氣膜孔下游圓弧長度s為60 mm，此為試驗有效測試區(qū)域。

圖1 試驗?zāi)Ｐ褪疽鈭DFig.1 Experimental model

試驗?zāi)Ｐ筒捎脤?dǎo)熱系數(shù)較小的夾布膠木材料加工而成。在進(jìn)行氣膜冷卻效率試驗時，為盡量準(zhǔn)確模擬絕熱壁面邊界條件，在試驗?zāi)Ｐ蛢?nèi)部采用增加空氣隔熱層及背面采用增加絕熱材料的方法，以減小試驗件導(dǎo)熱損失，降低絕熱壁溫測試誤差。

3 試驗狀態(tài)

以空氣作為冷卻工質(zhì)。選取D為特征尺度，利用試驗方法研究旋轉(zhuǎn)狀態(tài)下各主要影響因素(如吹風(fēng)比M、主流雷諾數(shù)ReD、旋轉(zhuǎn)數(shù)Rt等)對彎曲葉片氣膜冷卻效率η的影響規(guī)律。主要研究內(nèi)容為：

(1)確定M在0.2～2.0范圍時，研究M分別約為0.4、0.8、1.2和1.6工況下η的變化規(guī)律；

(2)確定ReD在1 800～5 000范圍時，研究對應(yīng)主流流量分別為200 kg/h、300 kg/h及400 kg/h時η的變化規(guī)律；

(3)確定Rt在0～0.020 0范圍時，研究對應(yīng)旋轉(zhuǎn)速度ω分別為0 r/min、500 r/min及800 r/min時η的變化規(guī)律。

4 試驗結(jié)果

4.1 吹風(fēng)比的影響

圖2、圖3分別給出了靜止和旋轉(zhuǎn)狀態(tài)下，不同M時模型表面η的二維等值線圖，圖中所示為ReD= 4 797時的試驗結(jié)果。如圖2所示，靜止?fàn)顟B(tài)時，氣膜冷卻軌跡沿主流方向中心線(s軸)呈對稱分布；并且隨著流向無量綱距離(s/D)的增加，各個M下η均單調(diào)下降。原因為隨著流動的發(fā)展，冷卻氣體逐漸耗散，其冷卻能力逐漸降低，導(dǎo)致η下降。從圖中還可明顯看出，隨著M的增大，氣膜軌跡覆蓋范圍明顯減小。這是因為M增大使得冷卻氣體的質(zhì)量流量相應(yīng)增加，其穿透主流的能力也逐漸增強(qiáng)，導(dǎo)致氣膜貼附壁面的能力下降，分離趨勢增強(qiáng)，從而對η產(chǎn)生較大影響。

如圖3所示，旋轉(zhuǎn)狀態(tài)時，在模型表面上旋轉(zhuǎn)使得氣膜沿徑向(z軸方向)發(fā)生偏轉(zhuǎn)，且偏向于高半徑(+z)方向，這是射流受到旋轉(zhuǎn)哥氏力和主流綜合作用的結(jié)果。理論分析可知，在試驗?zāi)Ｐ偷耐贡砻嫔?對應(yīng)渦輪葉片的吸力面)，哥氏力指向高半徑方向，使得冷氣在噴出氣膜孔后有向該方向偏轉(zhuǎn)的趨勢；同時，由于壁面曲率的影響，主流速度方向與旋轉(zhuǎn)角速度方向存在一定夾角，使得主流流體也因旋轉(zhuǎn)哥氏力的存在明顯影響z方向的分速度，從而影響射流軌跡。從數(shù)值上看，隨著M的增加，旋轉(zhuǎn)狀態(tài)下的η略有下降，但下降趨勢沒有靜止?fàn)顟B(tài)時的明顯。

4.2 旋轉(zhuǎn)數(shù)的影響

圖4給出了ReD=4 797、Rt=0.015 9條件下η的二維等值線圖。對比圖2(a)、圖3(a)和圖4(a)，可清楚看到旋轉(zhuǎn)對氣膜覆蓋區(qū)的影響。靜止?fàn)顟B(tài)下冷氣在氣膜孔下游的覆蓋區(qū)域分布基本對稱，且隨著s/D的增加，η逐漸降低。受旋轉(zhuǎn)哥氏力和主流流場的綜合作用，旋轉(zhuǎn)狀態(tài)時氣膜明顯向高半徑方向偏轉(zhuǎn)，且隨著Rt的增加，氣膜偏轉(zhuǎn)更加明顯。

對于較大吹風(fēng)比的情況(M=1.2)，如圖2(c)、圖3(c)和圖4(b)所示，η的分布規(guī)律與M=0.4時的結(jié)果較為相似，同樣可看到氣膜在旋轉(zhuǎn)狀態(tài)下的偏轉(zhuǎn)，且轉(zhuǎn)速越高偏轉(zhuǎn)越明顯。

圖2 靜止?fàn)顟B(tài)下不同吹風(fēng)比時氣膜冷卻效率的二維等值線圖(Rt=0，ReD=4 797)Fig.2 2D contours of film cooling efficiency under stationary status with different blowing ratio(Rt=0，ReD=4 797)

圖3 旋轉(zhuǎn)狀態(tài)下不同吹風(fēng)比時氣膜冷卻效率的二維等值線圖(Rt=0.01，ReD=4 797)Fig.3 2D contours of film cooling efficiency under rotational status with different blowing ratio(Rt=0.01，ReD=4 797)

圖4 Rt=0.015 9時氣膜冷卻效率的二維等值線圖(ReD=4 797)Fig.4 2D contours of film cooling efficiency (ReD=4 797，Rt=0.015 9)

圖5 為不同旋轉(zhuǎn)數(shù)時側(cè)向平均的η隨s/D的分布曲線?？梢?，當(dāng)M=0.4時，靜止?fàn)顟B(tài)下η的下降趨勢較緩，但η值相對較高，證明此時氣膜可較好地貼附壁面，冷氣與主流的摻混較弱。隨著Rt的增加，如圖5(a)中的Rt=0.010 0、0.015 9，η從氣膜孔出口處便開始急劇下降，這是因為旋轉(zhuǎn)使得氣膜附壁能力下降所致。此外，由于旋轉(zhuǎn)加劇了主流與冷氣的摻混，降低了冷氣的冷卻能力，加劇了氣膜冷卻效率的降低。

對比圖5可見，當(dāng)M增大到1.2時，Rt對η的影響與M=0.4時相似。對于靜止?fàn)顟B(tài)，η下降較平緩，且數(shù)值上高于旋轉(zhuǎn)狀態(tài)(Rt=0.015 9的0～1范圍內(nèi)除外)；旋轉(zhuǎn)狀態(tài)下，η隨s/D明顯降低。這說明M相同時，旋轉(zhuǎn)狀態(tài)下的氣膜較靜止時更易發(fā)生分離。

4.3 主流雷諾數(shù)的影響D

圖6給出了M=0.4，Rt=0.010 0時不同ReD下η的二維等值線圖。結(jié)合圖3(a)可看出，相同旋轉(zhuǎn)數(shù)和吹風(fēng)比條件下，ReD的改變對氣膜偏轉(zhuǎn)沒有產(chǎn)生較大影響，但卻使得主流流體的流動狀態(tài)發(fā)生了變化，因此會在一定程度上影響η的分布。

圖5 不同旋轉(zhuǎn)數(shù)時氣膜冷卻效率隨s/D的分布(ReD=4 797)Fig.5 Average film cooling efficiency vs.s/Dwith different rotation number(ReD=4 797)

圖7 給出了不同ReD時η隨s/D的分布曲線。從圖中看，M=0.4時，ReD的升高使得η先上升后大幅降低。這是因為在M不變時，ReD的升高增強(qiáng)了主流對射流的抑制，使得冷卻氣體能更好地貼附壁面，增大了氣膜冷卻效率；而當(dāng)ReD進(jìn)一步升高時，主流的擾動能力越來越強(qiáng)，增強(qiáng)了主流與射流的摻混，降低了射流的冷卻能力。

M=0.8時，ReD的影響有所下降。當(dāng)ReD升高時，η先是變化不大，而后降低。不過ReD的這一影響在氣膜孔附近已變得不太明顯，只在遠(yuǎn)離氣膜孔的下游區(qū)可清楚觀察到。這是因為隨著M的增大，射流的能力也有所增加，因此在ReD較小時射流對于主流特性的變化不太敏感，而當(dāng)ReD進(jìn)一步增大時才對射流產(chǎn)生較大影響。在遠(yuǎn)離氣膜孔位置，由于射流的耗散，其保持自身穩(wěn)定性的能力也逐漸降低，因此更易受到主流流動的影響。

圖6 不同主流雷諾數(shù)時氣膜冷卻效率的二維等值線圖(M=0.4，Rt=0.010 0) Fig.6 2D contours of film cooling efficiency with different mainstream Reynolds number(M=0.4，Rt=0.010 0)

圖7 不同主流雷諾數(shù)時氣膜冷卻效率隨s/D的分布(Rt=0.010 0)Fig.7 Average film cooling efficiency vs.s/Dwith different mainstream Reynolds number(Rt=0.010 0)

M=1.2和1.6時，ReD對氣膜孔下游的影響較小，但對遠(yuǎn)離氣膜孔位置的冷卻效率的影響十分顯著。這是因為隨著M的增大，冷氣的穿透能力也在增加，且在氣膜孔下游某個位置發(fā)生分離，此時主流擾動能力的增強(qiáng)會進(jìn)一步增強(qiáng)射流與主流的摻混，因此η明顯降低。

5 結(jié)論

(1)冷卻氣體噴出氣膜孔后，在向下游的發(fā)展過程中逐漸耗散，主流與冷氣的摻混逐漸增強(qiáng)，氣膜冷卻效率隨流向無量綱弧長逐漸降低，氣膜覆蓋范圍逐漸減小。

(2)吹風(fēng)比是影響氣膜冷卻效率的一個重要因素。隨著吹風(fēng)比的增大，在靜止?fàn)顟B(tài)下，逐漸增強(qiáng)的冷氣動量及穿透能力使得氣膜孔的冷卻效率逐漸降低；在旋轉(zhuǎn)狀態(tài)下，吹風(fēng)比對氣膜冷卻效率的影響相對靜止?fàn)顟B(tài)明顯減小，且隨著轉(zhuǎn)速的升高，吹風(fēng)比的影響進(jìn)一步降低。

(3)旋轉(zhuǎn)數(shù)會對氣膜冷卻效率的分布產(chǎn)生兩個主要影響：一是使得氣膜覆蓋區(qū)域向高半徑方向偏轉(zhuǎn)，且旋轉(zhuǎn)數(shù)越大，偏轉(zhuǎn)角度越大；二是旋轉(zhuǎn)數(shù)增加使得氣膜冷卻效率逐漸降低。

(4)旋轉(zhuǎn)狀態(tài)下主流雷諾數(shù)的改變對于氣膜軌跡偏轉(zhuǎn)沒有影響，但較高的主流雷諾數(shù)會使得空氣冷卻時的氣膜冷卻效率明顯降低。

[1]Ito S，Goldstein R J，Eckert E R.Film Cooling of a Gas Turbine Blade[J].Journal of Engineering for Power，1978，100：476—481.

[2]Drost U，Bolcs A，Hoffs A.Utilization of the Transient Liq?uid Crystal Technique for Film Cooling Effectiveness and Heat Transfer Investigations on a Flat Plate and a Turbine Airfoil[R].ASME 97-GT-26，1997.

[3]Gritsch M，Colban W，Schar H.Effect of Hole Geometry on the Thermal Performance of Fan-Shaped Film Cooling Holes[J].Journal of Turbomachinary，2005，127(4)：718—725.

[4]曾軍，王彬，康涌.氣膜冷卻渦輪導(dǎo)向葉片流場數(shù)值模擬[J].燃?xì)鉁u輪試驗與研究，2006，19(4)：16—19.

[5]朱惠人，劉松齡，余志紅.渦輪葉片氣膜冷卻的數(shù)值模擬[J].航空學(xué)報，1999，20(5)：416—420.

[6]吳宏，孟恒輝，陶智.旋轉(zhuǎn)對曲率表面氣膜冷卻效率影響的數(shù)值研究[J].航空學(xué)報，2009，30(9)：1624—1629.

[7]向安定，羅小強(qiáng)，朱惠人.渦輪葉片表面氣膜冷卻的傳熱試驗研究[J].航空動力學(xué)報，2002，17(5)：577—581.

Experimental Research of Film-Hole Cooling Efficiency in Rotation Status

PAN Bing-hua1，REN Fang1，GUO Wen1，TAO Zhi2
(
1.China Gas Turbine Establishment，Chengdu 610500，China；2.Beijing University of Aeronautics and Astronautics，Beijing 100191，China)

Experimental investigations were carried out to study the influence of rotation number,blowing ratio and mainstream Reynolds number on film cooling efficiency by simulating engine working condition. The results show that rotation will make the film coverage region lean to higher radial direction,the deflec?tion angle will become larger and the film cooling efficiency lower as the rotation number gets greater.At the same time,rotation will weaken the effect of blowing ratio and mainstream Reynolds number on film cooling efficiency.The test data of film-hole cooling characteristics at the typical regions of HPT surface gained through the study could be referential to the experiment database of the engine blade cooling design. Key words:aero-engine；turbine blade；film cooling；blowing ratio；rotation number

V231.3

A

1672-2620(2013)01-0030-05

2012-05-14；

2012-09-12

潘炳華(1976-)，男，重慶合川人，高級工程師，碩士，主要從事航空發(fā)動機(jī)空氣系統(tǒng)與熱分析研究。