苗 輝，朱江楠，馬薏文，劉志剛，張曉東

(1.中國航空發(fā)動機(jī)研究院，北京 101304；2.中國航發(fā)四川燃?xì)鉁u輪研究院，成都 610500)

1 引言

在航空發(fā)動機(jī)實際工作中，燃燒室出口截面(即渦輪進(jìn)口截面)的燃?xì)鉁囟却嬖趪?yán)重的非均勻性，該截面最高溫度會明顯高于平均溫度，差值高達(dá)100～200 K，這種局部的高溫氣流被稱作熱斑[1]。熱斑的存在顯著增加了渦輪的熱負(fù)荷，極大地影響著渦輪的可靠性和使用壽命。充分掌握渦輪進(jìn)口熱斑的影響規(guī)律及影響因素，是渦輪部件冷卻結(jié)構(gòu)精細(xì)化設(shè)計的必要條件。國外學(xué)者已經(jīng)通過數(shù)值與試驗研究對渦輪進(jìn)口熱斑遷移特性開展了探討。其中最早被研究的一個對熱斑有明顯影響的因素是熱斑引入位置與葉片之間的相對位置[2]，其他相關(guān)影響因素還有熱斑的形狀[3-4]、熱斑與主流氣體溫度比[5]、葉柵結(jié)構(gòu)形狀[6-9]以及冷卻流[10-11]等。國內(nèi)有關(guān)熱斑影響的研究起步相對較晚，2000年左右的才開始開展了一些數(shù)值與試驗研究[12-15]。

熱斑引入的周向、徑向位置差異會造成渦輪葉柵內(nèi)溫度場的不同表現(xiàn)。周向位置上，在熱斑完全沖擊第一級靜葉時，熱斑熱流與靜葉尾跡混合后作用在動葉上的熱負(fù)荷明顯變小，第二級導(dǎo)葉上的熱負(fù)荷也可被減輕，但此時付出了第一級導(dǎo)葉熱負(fù)荷激增的代價[16]。徑向位置上，在熱斑往葉尖方向移動時，輪轂附近的熱負(fù)荷逐漸降低，渦輪機(jī)匣附近熱負(fù)荷則上升。綜合考慮理想的熱斑引入位置時，最佳選擇應(yīng)該是將熱斑直接正向沖擊在第一級導(dǎo)葉靠近葉根位置[17]。

分析過往研究發(fā)現(xiàn)，大多工作圍繞單一熱斑開展，但實際情況中，從燃燒室后端引入渦輪的熱斑數(shù)量、形狀、大小、范圍都不盡相同。如對于橢圓形熱斑來說，它與周圍主流有更大的接觸面積，并在徑向上比周向上摻混得更加劇烈[18-19]。熱斑溫比在不同工況下的浮動范圍非常大，對下游部件產(chǎn)生的直接熱負(fù)荷分布不均勻。由于動葉壓力面附近有間歇性的流動分離，一層薄低溫流會出現(xiàn)在熱斑流與壁面之間，當(dāng)熱斑溫比上升時，流動分離的頻率也會隨之增加[20]。

近期的研究越來越多地關(guān)注葉柵彎扭、非軸對稱端壁和葉頂間隙等葉柵結(jié)構(gòu)對熱斑遷移的影響。通過控制靜葉出口流場來調(diào)整熱斑在動葉中的遷移是削弱熱斑產(chǎn)生高熱負(fù)荷區(qū)域的有效方法，葉柵彎扭等結(jié)構(gòu)能夠顯著改變靜葉出口流場分布，降低葉柵氣動損失[21]。正彎葉型可以顯著削弱熱斑向動葉壓力面的徑向遷移，使動葉壓力面和葉頂熱負(fù)荷顯著降低。葉頂間隙泄漏流動也對熱斑的徑向遷移有很大影響，其增強(qiáng)了熱斑向動葉壁面的遷移，使動葉壁面時均溫度增加[22]。使用非軸對稱端壁產(chǎn)生的附加二次流可以有效削弱二次流損失，但對動葉壁面靜溫的影響很小。另外，葉片數(shù)作為渦輪葉柵設(shè)計的一個關(guān)鍵參數(shù)，對渦輪氣熱性能影響很大，導(dǎo)葉數(shù)較多時，帶有熱斑的動葉表面溫度越低，并且存在最優(yōu)的動葉數(shù)使動葉壁面時均熱負(fù)荷最低[23-24]。

對于渦輪入口處溫度場畸變的影響，熱斑特性參數(shù)和葉型結(jié)構(gòu)的合理選擇可在一定程度上削弱熱斑產(chǎn)生的高熱負(fù)荷區(qū)域，但目前可行性最高且應(yīng)用較多的方法是采用先進(jìn)的氣膜冷卻技術(shù)。對比其他冷卻方法來講，氣膜冷卻結(jié)構(gòu)效率更高，對冷卻氣流速度的限制更少，便于有針對性地減小熱斑負(fù)面影響[25]。氣膜冷卻對熱斑的削弱作用在不同位置差異很大，前緣冷卻氣流主要削弱輪轂區(qū)域的高溫區(qū)，壓力面冷氣對熱斑的影響甚微，而吸力面上的氣膜幾乎可消除熱斑影響。

綜上，國內(nèi)外針對熱斑問題開展的工作，鮮有涉及熱斑壓力與主流壓力之比(本文定義為熱斑壓力比)的影響規(guī)律的研究。在實際發(fā)動機(jī)工作條件下，由于燃燒室火焰筒二次孔橫向射流造成的速度不均勻性會導(dǎo)致渦輪進(jìn)口截面的總壓出現(xiàn)不均勻性。這種不均勻性雖然不大(如±3%)，但其與溫度不均勻性同時存在，對渦輪部件會有明顯的影響。本文針對圓形熱斑對渦輪葉柵的影響進(jìn)行試驗和數(shù)值研究，以便深化對熱斑效應(yīng)的認(rèn)識，為更精細(xì)化地設(shè)計渦輪冷卻系統(tǒng)、合理減少冷卻氣量提供依據(jù)。

2 試驗介紹

2.1 試驗系統(tǒng)

試驗系統(tǒng)內(nèi)部結(jié)構(gòu)見圖1，主要由進(jìn)口延長段、轉(zhuǎn)接段、湍流熱斑段、試驗段等組成。轉(zhuǎn)接段的作用是將延長段中的流體引入流道；湍流熱斑段的作用是通過擾流棒和次流引氣管來模擬試驗所需的進(jìn)口湍流度和熱斑；試驗段的作用是固定試驗葉柵和冷氣集氣腔，并通過其上安裝的測試受感部采集葉柵進(jìn)出口數(shù)據(jù)。需說明的是，雖然真實發(fā)動機(jī)的熱斑溫度本身并不均勻，試驗中次流引氣管以均勻溫度的方式來模擬熱斑條件時，保持試驗中的熱斑溫度比接近于實際發(fā)動機(jī)燃燒室出口截面熱斑平均溫度與整個截面平均溫度之比，且熱斑面積相似。詳細(xì)介紹見文獻(xiàn)[26]。

圖1 試驗系統(tǒng)示意圖Fig.1 Schematic views of test setup

2.2 試驗件

試驗中共有5 個渦輪葉片形成的4 個葉柵通道。5個葉片中只有1個葉片為測量葉片，其余為實心陪襯葉片；葉片材料為環(huán)氧樹脂，采用C3X模型，有7 排氣膜孔，其中前緣3 排，壓力面和吸力面各2排；葉片表面有為鋪設(shè)熱偶絲而切出的槽道。渦輪葉片模型和試驗件照片如圖2所示。

圖2 試驗葉片(左-示意圖；右-照片)Fig.2 Test turbine vane(Left:schematic view;Right:photo)

2.3 試驗條件

本試驗研究的目的是為校核數(shù)值研究方案提供基礎(chǔ)數(shù)據(jù)。受試驗資源限制，試驗條件如表1 所示。試驗中，在試驗條件穩(wěn)定后錄取葉片表面各點熱電偶溫度數(shù)據(jù)。

表1 試驗條件Table 1 Test Conditions

3 數(shù)值計算方法

計算域如圖3 所示，取一個渦輪葉片周期為計算對象，沿兩側(cè)渦輪通道中心線形成的曲面為邊界，并取為周期性邊界。熱斑管模型與試驗條件的一致。由于葉片材料為低導(dǎo)熱率的環(huán)氧樹脂，壁面導(dǎo)熱可以忽略，故計算域中只有流體部分。主流和熱斑從葉柵前面進(jìn)入葉柵通道；冷卻氣流從葉片下部進(jìn)入葉片內(nèi)部，經(jīng)由氣膜孔進(jìn)入主流。

圖3 計算域Fig.3 Computational domain

采用網(wǎng)格劃分軟件生成的非結(jié)構(gòu)化四面體網(wǎng)格，在保證計算精度的條件下，具有快速分網(wǎng)的特點。對葉片表面附近網(wǎng)格進(jìn)行加密，并附著邊界層。網(wǎng)格無關(guān)性驗證所用總網(wǎng)格量分別是580 萬、1 100萬和1 600萬，以前兩套和后兩套網(wǎng)格方案計算的葉片平均溫度相差分別為2.1%和0.3%，故選擇1 100萬網(wǎng)格方案。

數(shù)值算法選用CFD商業(yè)軟件進(jìn)行求解，計算采用Faver 平均的N-S 方程。對流項離散采用二階精度迎風(fēng)格式，擴(kuò)散項為中心差分格式，時間迭代采用穩(wěn)態(tài)的基于密度的隱式耦合算法。湍流模型選用SSTk-ω模型。邊界條件為壓力入口/壓力出口條件，葉片表面及計算域底面/頂面為絕熱壁面。

4 結(jié)果與討論

4.1 試驗結(jié)果和計算結(jié)果的對比

針對試驗工況進(jìn)行計算，邊界條件按照試驗條件設(shè)置，對應(yīng)的葉片表面溫度場如圖4 所示。可看出，熱斑直接沖擊在葉片前緣，雖然在前緣附近有3排氣膜孔，但氣膜孔未能給前緣提供良好的保護(hù)，形成了高溫區(qū)。另外，熱斑正對葉片前緣，但熱斑氣流對葉片吸力面的影響更大。在壓力面上經(jīng)過氣膜孔后熱斑氣流已沒有明顯影響，而吸力面上在氣膜孔間覆蓋的邊緣區(qū)形成了多個溫度較高的狹長區(qū)域。

圖4 試驗條件(熱斑壓力比1.028，冷氣流量比4%)下的溫度場云圖Fig.4 Numerical temperature contour of vane under test condition with pressure ratio of 1.028 and cooling air ratio of 4%

圖5 為葉柵通道中截面氣流的溫度與速度云圖。可看出，熱斑氣流受位勢效應(yīng)的影響，向吸力面發(fā)生較為明顯的偏轉(zhuǎn)。速度最高區(qū)域出現(xiàn)在吸力面鰓區(qū)，壓力面靠近葉片前緣區(qū)域則出現(xiàn)低速區(qū)。

圖5 葉片中截面流體域的溫度場與速度場Fig.5 Temperature and velocity contour at the middle section plane of the vane

圖6示出了葉片中截面表面溫度計算結(jié)果與試驗結(jié)果的對比?？煽闯?，中截面上計算結(jié)果有數(shù)個溫度陡降區(qū)，這是因為氣膜射流的注入使下游臨近區(qū)域溫度得到明顯降低。而試驗中葉片表面測點要避開氣膜孔區(qū)域，故試驗結(jié)果未有溫度突然下降。經(jīng)分析，壁面溫度的試驗結(jié)果和對應(yīng)計算結(jié)果的平均誤差為1.05%，認(rèn)為計算方案足夠精確。

圖6 試驗和數(shù)值計算的中截面溫度Fig.6 Experimental and numerical calculation results of middle section plane of the vane

4.2 熱斑壓力比的影響

由于試驗中未考慮設(shè)置總壓可變，因此采用數(shù)值計算方法研究熱斑壓力比的影響。試驗中熱斑壓力比為1.040，但熱斑氣流經(jīng)過彎管會產(chǎn)生壓力損失，葉柵入口處熱斑壓力比約為1.028。數(shù)值計算中設(shè)置熱斑壓力比變化范圍在相差±2%(即熱斑壓力比范圍為0.980～1.020)，其中設(shè)定主流壓力不變，熱斑進(jìn)口壓力根據(jù)熱斑壓力比計算；冷氣流量比選擇2%和4%；其余參數(shù)與表1中的相同。圖7示出了不同熱斑壓力比條件下葉片表面的平均溫度。可看出，隨著熱斑壓力比的增加，熱斑使葉片表面的平均溫度逐漸增加。其中，壓力面平均溫度受熱斑壓力比變化的影響較小，吸力面則較大，尤其是熱斑壓力比較高時。

圖7 不同壓力比條件下葉片壓力面和吸力面的平均溫度Fig.7 Average temperature of pressure side and suction side of vane under different perssure ratios

圖8 示出了熱斑壓力比分別為0.980、1.000、1.010、1.020，冷氣流量比為4%時葉片表面的溫度場。從圖8及圖4可看出，隨著熱斑壓力比的增加，溫度較高的區(qū)域逐漸擴(kuò)大。這是因為熱斑壓力比增加，引起熱斑射流速度和流量的增加，高溫?zé)岚呱淞骺梢愿采w葉片表面更多的面積。另外，對于此類熱斑正對葉片前緣的計算條件，當(dāng)熱斑壓力比較小(0.980、1.000)時，熱斑對帶氣膜射流的葉片無明顯影響；而隨著熱斑壓力比的增加(壓力比為1.010、1.020、1.028)，熱斑逐漸對葉片表面產(chǎn)生影響，造成前緣和吸力面出現(xiàn)高溫區(qū)。

圖8 冷氣流量比為4%時不同壓力比條件下的葉片溫度場云圖Fig.8 Temperature contour of vane under different pressure ratio with the cooling air ratio of 4%

圖9 示出了熱斑壓力比分別為0.980、1.000 和1.020，計算域在50%葉高水平面上的溫度場云圖?？煽闯觯m然熱斑射流正對葉片前緣，但受位勢效應(yīng)影響，熱斑射流向吸力面一側(cè)偏轉(zhuǎn)。熱斑在渦輪通道中向下游流動時，溫度不斷下降。熱斑壓力比越小，熱斑溫度衰減越快，方向偏轉(zhuǎn)也越明顯。其中，當(dāng)熱斑壓力比為0.980時，熱斑氣流未及前緣區(qū)域，其溫度便衰減至和主流相當(dāng)?shù)臏囟人健?/p>

圖9 冷氣流量比為4%時不同壓力比條件下葉片中截面流體域溫度場Fig.9 Temperature contour at the middle section plane of the vane under different pressure ratio with the cooling air ratio of 4%

從上述分析可看出，熱斑壓力比的微小變化對熱斑特征會產(chǎn)生明顯的影響。在前人的研究中未考慮這一因素，也許是造成研究結(jié)果偏差較大的原因之一。由于實際發(fā)動機(jī)中，渦輪前截面的溫度不均勻性(熱斑)與總壓不均勻性是相互獨立的，且各自的影響規(guī)律又較為復(fù)雜，用于指導(dǎo)渦輪冷卻結(jié)構(gòu)設(shè)計時，可按照最高壓力比校核熱斑的影響。

4.3 冷氣流量比的影響

圖10 示出了冷氣流量比為1%～5%，熱斑壓力比為1.020，其余參數(shù)與表1 試驗條件相同下，葉片壓力面和吸力面的平均溫度?？煽闯?，隨著冷氣流量比增加，葉片表面溫度顯著下降。

圖10 不同冷氣流量比條件下葉片壓力面和吸力面的平均溫度Fig.10 Average temperature of pressure side and suction side of vane under different cooling air ratios

圖11示出了熱斑壓力比1.020、冷氣流量比2%條件下葉片表面的溫度場?？梢姡c冷氣流量比為4%的圖8(d)相比，冷氣流量比越小，氣膜對葉片表面的覆蓋越差，葉片前緣和吸力面的高溫區(qū)域面積越大。此外，熱斑壓力比較高(如1.020)時，即使冷氣流量比高達(dá)4%，本文選用的氣膜冷卻結(jié)構(gòu)也未能很好地阻止熱斑對葉片前緣和吸力面的加熱。為此，在實際渦輪設(shè)計中，可以根據(jù)熱斑的規(guī)律對冷卻結(jié)構(gòu)進(jìn)行針對性的設(shè)計或優(yōu)化。

圖11 壓力比1.020時冷氣流量比2%條件下的葉片表面溫度場Fig.11 Temperature contour of vane with pressure ratio of 1.020 and cooling air ratio of 2%

5 結(jié)論

(1)在壓力比1.028、冷氣流量比4%條件下，數(shù)值計算結(jié)果和試驗結(jié)果吻合較好，各測點試驗與計算的壁面溫度結(jié)果的平均誤差為1.05%。

(2)熱斑壓力比的微小變化對熱斑特征會產(chǎn)生明顯的影響。在計算條件下，熱斑壓力比小于1.000 時，熱斑對葉片表面熱負(fù)荷幾乎無影響；但隨著熱斑壓力比的增加，熱斑對葉片表面的加熱效果越來越明顯，尤其是在葉片前緣和吸力面。

(3)受位勢效應(yīng)影響，熱斑射流向吸力面一側(cè)偏轉(zhuǎn)。熱斑在渦輪通道中向下游流動時，溫度不斷下降。熱斑壓力比越小，熱斑溫度衰減越快，方向偏轉(zhuǎn)也越明顯。當(dāng)熱斑壓力比為0.980 時，熱斑氣流未及前緣區(qū)域，其溫度便衰減至和主流相當(dāng)?shù)臏囟人健?/p>

(4)考慮到實際發(fā)動機(jī)渦輪前截面的溫度不均勻性(熱斑)與總壓不均勻性是相互獨立的，為應(yīng)對熱斑效應(yīng)給葉片帶來的額外熱負(fù)荷，在渦輪冷卻結(jié)構(gòu)設(shè)計時，可按照最高壓力比來校核熱斑的影響。