內部橫流對外部氣膜冷卻特性的影響

駱劍霞，朱惠人，劉存良，賈廣森

（西北工業大學動力與能源學院，陜西 西安710072）

0 引言

在現代航空發動機中，渦輪前溫度已高于葉片材料的熔點，因而必須采用高效的冷卻技術。大多數渦輪葉片冷卻方式均為內外綜合冷卻，包括復雜的內部冷卻通道與外部的氣膜孔結構。目前，一般研究大多都是針對單獨內部冷卻或單獨外部冷卻展開，兩者分別有著大量的研究成果。但是關于內外冷卻結構之間相互影響的研究比較少。掌握內外流之間的相互影響機理，協調優化內外冷卻結構，才能設計出更高效的葉片冷卻方式。

內外流之間的相互影響可以分為2部分內容［1－11］，一是關于氣膜孔出流對內部流動與換熱影響的研究，二是關于內部流動對外部氣膜冷卻特性影響的研究。在此，分析對比了內部橫流對外部氣膜冷卻特性的影響。

1 研究方法

1．1 實驗系統

實驗系統結構如圖1所示。在系統中，二次流與主流分別由離心風機提供。

主流流入實驗臺穩定段和收縮段A后，均勻地通過快速加熱器，在加熱器內得到均勻一致的加熱?？焖偌訜崞髦筮€有一段收縮段B，以確保均勻穩定的主流流入實驗測量段。實驗測量段橫截面尺寸為220 mm×80 mm。為了消除加熱器對主流邊界層帶來的影響，在實驗測量段入口處一側的測量板上布置有一刃狀劈縫。在氣膜孔中心上游140 mm處測得主流湍流度小于1%。

實驗通道中布置有5個氣膜孔，如圖1所示。氣膜孔直徑D為10 mm，氣膜孔長徑比L／D＝3，氣膜孔孔間距P／D＝4，氣膜孔流向傾角α為30°。

二次流流量由通道前后2個孔板流量計控制。二次流通道截面尺寸為70 mm×70 mm。

實驗條件：主流進口平均速度約為17 m／s，基于氣膜孔直徑的進口雷諾數ReD＝10 000；二次流進口平均速度約為23 m／s，基于二次流通道水力直徑的進口雷諾數ReDh為50 000和100 000。

實驗中氣膜孔吹風比為0．5，1．0和2．0。吹風比定義為：

ρc為二次流密度；Uc為二次流孔內平均速度；ρm為主流密度；Um為主流進口平均流速。

坐標軸原點定義為中間氣膜孔出口中心，X軸為主流流動方向，Y軸為二次流通道橫流流動方向，Z軸為高度方向。

圖1 實驗系統

1．2 測量方法及誤差分析

實驗使用窄帶液晶測量氣膜孔下游壁面上的溫度，多次測量求解出表面換熱系數和冷卻效率。熱色液晶測量技術已被廣泛應用，關于其原理和方法，參考文獻［12］有詳細的介紹和分析。參考文獻［13］分析了使用瞬態液晶同時測量換熱系數和氣膜冷卻效率實驗中的不確定度，參考其研究方法，計算得出實驗中冷卻效率的不確定度為6．15%～10%，換熱系數的不確定度為5%～7．14%。

2 結果與分析

2．1 氣膜冷卻效率

各表面展向平均氣膜冷卻效率分布曲線如圖2所示。實驗測量所得不同條件下表面冷卻效率η分布云圖如圖3所示。圖中坐標X／D和Y／D為無量綱距離，其中，X，Y為該位置處橫向、展向與氣膜孔中心距離，D為氣膜孔直徑。

圖2 氣膜孔下游表面展向平均冷卻效率分布曲線

圖3 氣膜孔下游表面冷卻效率分布云圖

吹風比為0．5時，在不同的內部橫流條件下，冷卻效率分布均呈現為：高冷卻效率區域向＋Y側偏移，沿流動方向冷卻效率逐漸下降。而且內部橫流越大，偏斜越嚴重，氣膜展向分布越窄，沿流動方向下降越迅速。內部橫流的存在，使得進口氣流具有橫流方向速度分量。在小吹風比下，氣膜孔出流抽吸作用較弱，對橫流影響的削弱作用也較弱，橫流使得出口射流向橫流方向（＋Y方向）偏斜，從而影響了氣膜分布，橫流越大，這部分影響越明顯。由圖2可以看出，橫流ReDh＝100 000與ReDh＝50 000相比，展向平均冷卻效率下降了約20%。

吹風比為1時，氣膜孔出流抽吸作用增強，因而橫流作用相對減弱。橫流的作用使得氣膜孔出口＋Y側射流高度高于－Y側，因而氣膜孔下游－Y側壁面的冷卻效率較高。內部橫流越大，射流與主流的摻混作用越強，展向分布越寬，尤其是在靠近氣膜孔出口處，展向平均冷卻效率提高了50%，沿著流動方向該作用越來越弱。

吹風比為2時，氣膜孔射流的Z向分量較大，射流易沖出壁面。在ReDh＝50 000時，射流沖出壁面，冷卻效率很低。在ReDh＝100 000時，橫流作用增強了氣膜孔出口射流與主流的摻混作用，這使得一部分冷氣仍能貼附壁面，冷卻效率相對較高。與小橫流條件對比，大橫流條件下展向平均冷卻效率提高了約180%。

在各橫流條件下，隨著吹風比的增加，冷卻效率均表現為先增加后下降。

2．2 換熱系數

實驗測量所得不同條件下表面換熱系數（hf／h0）分布云圖如圖4所示。各表面展向平均換熱系數比分布曲線如圖5所示。

圖4 氣膜孔下游表面展向平均換熱系數比分布云圖

圖5 氣膜孔下游表面展向平均換熱系數比分布曲線

由圖4和圖5可以看出，不同的橫流雷諾數使得氣膜孔下游換熱系數分布有較大的區別。總體說來，橫流的作用使得換熱系數分布呈現非對稱性，橫流Re數越大，非對稱性越明顯。

根據參考文獻［14］，無橫流條件下，圓柱型氣膜孔出口射流與主流相互作用，形成一對對稱分布的外轉對轉渦。對轉渦在射流兩側的孔間區域可以形成類似于沖擊的流動結構，該結構使得相應區域的邊界層厚度減薄，從而增強該區域的對流換熱強度，因而強換熱區基本位于旋轉渦的正下方。同時，在2對轉渦內部中心區域，形成了與沖擊相反的流動結構，這種流動結構會減小壁面附近的速度梯度，從而使得邊界層厚度增加，削弱對流換熱強度，這種現象只有在射流速度較低或相鄰對轉渦比較接近的時候，才會對換熱產生明顯影響。

吹風比為0．5，ReDh＝50 000條件下，高換熱區域向橫流流動方向偏斜，在X／D＞8后，分成2個區域，其中，＋Y側區域換熱系數迅速下降，而另一區域換熱系數下降則相對較緩。這說明此條件下氣膜孔出口射流與主流相互摻混，形成的2個對轉漩渦強度不對等。橫流ReDh＝100 000條件下，高換熱系數區域偏斜非常嚴重，而且下游也沒有分成2個區域，說明這種情況下對轉渦的強度非常不對等。

吹風比為1，ReDh＝50 000條件下，高換熱系數區域沒有向＋Y側偏斜，下游高換熱系數也分成2個區域。2個區域中間的低換熱系數區域非常明顯，越往下游越低。橫流ReDh＝100 000條件下，高換熱區域仍只有1個，而且這部分換熱系數非常高。

吹風比為2，ReDh＝50 000條件下，換熱系數分布基本對稱，說明在該工況下橫流作用不明顯（氣膜孔抽吸作用相對橫流作用很強，橫流作用因而顯現不出來）。但是ReDh＝100 000條件下，出現了不對稱的高換熱系數分布，而且－Y側換熱系數值非常高。

ReDh＝50 000條件下，氣膜孔出口下游均出現了2個高換熱系數區域，隨著吹風比的增加，橫流的作用相對減弱，2個高換熱系數區域差別減小。ReDh＝100 000條件下，吹風比為0．5和1時，均只有1個高換熱系數區域（此時橫流的影響相對較強），而吹風比為2時才出現了2個不對稱的高換熱系數區域。

從換熱系數比展向平均曲線可以看出，橫流雷諾數增加，下游表面換熱系數增強；吹風比增大，表面換熱系數增強。

2．3 流量系數

不同橫流雷諾數條件下氣膜孔的流量系數結果如圖6所示。由圖6可以看出，隨著吹風比的增大，流量系數增加；而隨著橫流Re的增大，流量系數下降。內部橫流的影響主要是改變了氣膜孔進口速度分布，但是氣流在氣膜孔內流動會在一定程度上削弱橫流的影響，吹風比增大，氣膜孔抽吸作用增強，內部橫流的影響相對越弱。

圖6 氣膜孔流量系數分布曲線

3 結束語

研究了二次流通道內部橫流對外部氣膜冷卻特性的影響。通過以上分析，可以得到如下結論：

a．內部橫流使得氣膜孔下游表面冷卻效率和換熱系數分布產生偏斜。

b．吹風比為0．5時，小橫流條件下氣膜冷卻效率高于大橫流條件；隨著吹風比的增加，增強內部橫流反而增強了下游氣膜冷卻效率，尤其是在吹風比為2時，展向平均冷卻效率提高了約180%。

c．橫流雷諾數增加，下游表面換熱系數增強；吹風比增大，表面換熱系數增強。

d．橫流增強，氣膜孔流量系數降低。吹風比增大，流量系數增大。

［1］ 神家銳，Ireland P T，Wang Z，等．收斂通道內氣膜孔及其與矩形肋復合的強化傳熱［J］．工程熱物理學報，1992，13（2）：185－190．

［2］ Ekkad S V，Hung Y，Han J C．Detailed heat transfer distributions in two－pass square channels with rib turbulators and bleed holes［J］．International Journal of Heat and Mass Transfer，1998，41（23）：3781－3791．

［3］ Cukurel B，Selcan C，Arts T．Film cooling extraction effects on the areo－thermal characteristics of rib roughened cooling channel flow［R］．ASME Paper GT2012－68680，2012．

［4］ Kim K M，Park S H，Jeon Y H，et al．Heat／mass transfer characteristics in angled ribbed channels with various bleed ratios and rotation numbers［R］．ASME Paper GT2007－27166，2007．

［5］ Jeon Y H，Park S H，Kim K M，et al．Effects of bleed flow on heat／mass transfer in a rotating rib－roughened channel［R］．ASME Paper No．GT2006－91122，2006．

［6］ Kissel H P，Weigand B，von Wolfersdorf J，et al．An

experimental and numerical investigation of the effect of cooling channel crossflow on film cooling performance［R］．ASME Paper No．GT2007－27102，2007．

［7］ 彭 威，姜培學．冷卻通道內流動對氣膜冷卻影響的大渦模擬研究［J］．工程熱物理學報，2010，31（8）：1359－1362．

［8］ 白江濤，朱惠人，劉存良．內通道交錯橫流對氣膜冷卻效率的影響［J］．航空動力學報，2008，23（8）：1353－1358．

［9］ Gritsch M，Schulz A，Wittig S．Effect of crossflows on the discharge coefficient of film cooling holes with varying angles of inclination and orientation［J］．Journal of Turbomachinery，2001，123（4）：781－787．

［10］ Gritsch M，Saumweber C，Schulz A，et al．Effect of internal coolant crossflow orientation on the discharge coefficient of shaped film－cooling holes［J］．Journal of Turbomachinery，2000，122（1）：146－152．

［11］ Bunker R S，Bailey J C．Film cooling discharge coefficient measurements in a turbulated passage with internal crossflow ［J］．Journal of Turbomachinery，2001，123（4）：774－780．

［12］ Irelang P T，Jones T V．Liquid crystal measurements

of heat transfer and surface shear stress［J］．Measurement Science and Technology，2000，11（7）：969－986．

［13］ 白江濤，朱惠人，劉存良．雙參數傳熱實驗的液晶瞬態測量不確定度分析［J］．航空動力學報，2009，24（9）：1945－1951．

［14］ Walters D K，Leylek J H．A detailed analysis of filmcooling physics Part 1：streamwise injection with cylindrical holes［R］．ASME Paper 97－GT－269，1997．