張 洪，王建華，張清波

（中國(guó)科學(xué)技術(shù)大學(xué)熱科學(xué)與能源工程系，合肥 230027）

0 引 言

在現(xiàn)代燃?xì)馔钙饺~片設(shè)計(jì)中，動(dòng)葉片頂部應(yīng)該盡量接近透平環(huán)形的固定壁面，同時(shí)又必須留有一定間隙，以適應(yīng)葉片因離心力及受熱引起的伸長(zhǎng)和固定壁的熱變形。在葉片壓力面和吸力面之間的壓差驅(qū)動(dòng)下，來(lái)自燃燒室的高溫氣體穿過(guò)這個(gè)頂部間隙，稱(chēng)之為泄漏流。泄漏流將導(dǎo)致以下4個(gè)棘手問(wèn)題：

（1）泄漏流在壓力、吸力面之間的壓差驅(qū)動(dòng)下穿過(guò)葉片間隙后與通道中的主流氣體混和，產(chǎn)生泄漏渦，導(dǎo)致較大的總壓損失；

（2）由于粘性作用，泄漏流在固定環(huán)壁和葉頂表面形成流動(dòng)邊界層，邊界中的氣流速度比轉(zhuǎn)子在葉頂處的線(xiàn)速度小，轉(zhuǎn)子運(yùn)動(dòng)時(shí)受到這部分氣流的阻力作用又造成一定程度的摩擦損失；

（3）泄漏流直接穿過(guò)葉片間隙，沒(méi)有對(duì)葉片做功，因此產(chǎn)生透平級(jí)效率損失。這種損失在一定范圍內(nèi)隨著間隙距離的變化而不同。文獻(xiàn)［1］指出間隙為葉片高度1%時(shí)將產(chǎn)生1%～2%的主流泄漏損失，從而降低1%～3%的透平級(jí)效率；

（4）高溫泄漏流與葉片頂部壁面的對(duì)流換熱增大，使得整個(gè)葉片頂部表面熱負(fù)荷最高，從而導(dǎo)致葉片壽命縮短。

長(zhǎng)期以來(lái)，人們對(duì)葉頂間隙泄漏流特性及高溫泄漏流與葉片頂部、壓力面和吸力面換熱特性做了大量研究，并在此基礎(chǔ)上，改進(jìn)旋轉(zhuǎn)葉片頂部結(jié)構(gòu)設(shè)計(jì)，增加凹槽式迷宮密封及氣膜孔冷卻頂部，這樣既可以減少泄漏流量，又可以冷卻葉片頂部。文獻(xiàn)［2］詳細(xì)回顧了關(guān)于葉片頂部高溫氣體流動(dòng)和傳熱的研究情況。文獻(xiàn)［3］的研究結(jié)果指出葉片頂部的對(duì)流換熱與葉片和固定圍帶之間的相對(duì)運(yùn)動(dòng)無(wú)關(guān)；文獻(xiàn)［4］在較大的相對(duì)速度范圍內(nèi)，其中包括比實(shí)際情況更高的速度下通過(guò)實(shí)驗(yàn)測(cè)量和數(shù)值計(jì)算，結(jié)果都確定了這種無(wú)關(guān)性。

這種現(xiàn)象的本質(zhì)可以解釋為：相對(duì)運(yùn)動(dòng)的影響其實(shí)僅限于一個(gè)附著在圍帶壁面的薄層之內(nèi)，葉頂附近的速度分布實(shí)際上并沒(méi)有因?yàn)橄鄬?duì)運(yùn)動(dòng)而改變。這個(gè)結(jié)論使得對(duì)葉頂泄漏流和傳熱的實(shí)驗(yàn)研究可以在靜態(tài)裝置中進(jìn)行。

在以往的泄漏流特性實(shí)驗(yàn)研究中，不同的測(cè)量技術(shù)被使用。文獻(xiàn)［5］使用激光粒子圖像測(cè)速技術(shù)（particle image velocimetry），觀(guān)察到頂部間隙泄漏流和主流的混合情況。文獻(xiàn)［6－9］利用萘升華技術(shù)，測(cè)量了葉片頂部的局部和平均質(zhì)量傳輸系數(shù)。

在數(shù)值計(jì)算方面關(guān)于葉片間隙區(qū)域的傳熱傳質(zhì)特性也有許多研究報(bào)告。文獻(xiàn)［10］采用標(biāo)準(zhǔn)k－ω湍流模型，模擬了帶有頂部凹槽的透平葉片的三維流場(chǎng)和傳熱特性。文獻(xiàn)［11］利用商業(yè)軟件Fluent，用不同湍流模型進(jìn)行數(shù)值計(jì)算，并指出：雷諾應(yīng)力模型下的結(jié)果和實(shí)驗(yàn)數(shù)據(jù)比較吻合。文獻(xiàn)［12］比較了7種不同的湍流模型和不同的網(wǎng)格劃分方法，最后指出k－ω模型的結(jié)果和實(shí)驗(yàn)結(jié)果最為接近，同時(shí)他們還分析了葉頂形狀對(duì)流動(dòng)和傳熱特性的影響。

實(shí)驗(yàn)在靜止通道內(nèi)進(jìn)行，利用PIV技術(shù)獲得不同截面上的速度矢量分布圖，同時(shí)采用CFX12.0計(jì)算5種湍流模型下的流場(chǎng)分布，并將實(shí)驗(yàn)結(jié)果與數(shù)值結(jié)果進(jìn)行比較，由此確定與實(shí)際流動(dòng)趨勢(shì)最吻合的湍流模型。為下一步進(jìn)行流－固－熱耦合計(jì)算頂部傳熱特性提供依據(jù)。

1 實(shí)驗(yàn)裝置

1.1 實(shí)驗(yàn)管道和葉片模型

為了有利于捕捉內(nèi)部流場(chǎng)粒子圖像，實(shí)驗(yàn)管道用透射率非常好的有機(jī)玻璃做成。管道進(jìn)口截面為164mm×72mm，并在進(jìn)口段放置一蜂窩器，以降低來(lái)流湍流度。蜂窩器由8mm×8mm的方孔組成，孔間距為3mm，長(zhǎng)度為40mm。試驗(yàn)件是經(jīng)過(guò)發(fā)黑處理并放大3倍的GE－E3不銹鋼葉片，高度為60mm，葉片頂部開(kāi)有深度為4.5mm的凹槽。頂部沿葉片型面中軸線(xiàn)按7.5mm等間距分布11個(gè)氣膜孔，壓力側(cè)凸臺(tái)上按等間距分布5個(gè)氣膜孔，每個(gè)孔直徑1.5mm（如圖1所示）。為了滿(mǎn)足周期性條件，試驗(yàn)件兩側(cè)均用有機(jī)玻璃做成了葉片形狀，并與通道連接（如圖2所示），葉片的間距為87mm，葉片與來(lái)流間的夾角為58°，葉片與出口段的夾角為24°。中間葉片與管道上表面有一個(gè)間隙，以模擬真實(shí)燃?xì)馔钙饺~片頂部間隙。葉片上方的管道處開(kāi)了一個(gè)多邊形的孔，在孔處加一個(gè)活動(dòng)蓋板，這樣可以方便改變?cè)囼?yàn)件與壁面的間隙。本實(shí)驗(yàn)因?yàn)橐臄z葉片頂部間隙處的流場(chǎng)，將間隙值設(shè)置為12mm。本次研究不涉及傳熱和物性相關(guān)方面的問(wèn)題，主要針對(duì)流場(chǎng)特性，因此，在保證與真實(shí)條件相同吹風(fēng)比的條件下，依據(jù)雷諾相似定律，在常溫（293K）下進(jìn)行的水流模擬實(shí)驗(yàn)可以反應(yīng)相同雷諾數(shù)下的氣體流動(dòng)特性。具體實(shí)驗(yàn)參數(shù)如表1所示。

圖1 葉片外形示意圖Fig.1 Schematic of blade shape

圖2 實(shí)驗(yàn)管道截面示意圖Fig.2 Schematic of test section

表1 實(shí)驗(yàn)參數(shù)Table 1 Experiment parameters

1.2 PIV系統(tǒng)

實(shí)驗(yàn)采用的粒子測(cè)速（PIV）系統(tǒng)由Nd：YAG激光器、片光發(fā)生器、海信CCD攝像系統(tǒng)以及配套的圖像處理軟件（Dantec Flowmap 1500）組成。Nd：YAG激光器發(fā)射紅外線(xiàn)，并可延伸至532nm的綠色可見(jiàn)光波長(zhǎng)段，激光脈沖間隔可控制在300ns～300ms之間。實(shí)驗(yàn)用示蹤粒子為平均直徑10μm的中空球，由光敏材料制造，球與水的質(zhì)量比為1.05，可以跟隨流體運(yùn)動(dòng)。這些中空球既能較為準(zhǔn)確地顯示流場(chǎng)運(yùn)動(dòng)特性，又可以清晰地被CCD相機(jī)捕捉。實(shí)驗(yàn)時(shí)，激光器觸發(fā)的激光脈沖間隔選取為3ms，CCD像素為1024pixel×768pixel，每一種測(cè)量狀態(tài)連續(xù)拍攝150幅圖片，然后用系統(tǒng)軟件分析得到所拍攝截面上的二維速度分布。

2 數(shù)值計(jì)算模型及方法

采用商業(yè)軟件ANSYS CFX 12.0進(jìn)行數(shù)值計(jì)算，邊界處理和模型尺寸與實(shí)驗(yàn)條件一致。進(jìn)口段長(zhǎng)度取一倍葉片軸向弦長(zhǎng)，出口段延伸至足夠長(zhǎng)以防止出現(xiàn)回流。

計(jì)算的邊界條件與實(shí)驗(yàn)條件一致，流動(dòng)介質(zhì)為293K的水，進(jìn)口條件給定流量入口，流率為1.0586kg／s，給定壓力出口條件為常壓，頂部氣膜孔出口速度為0.3m／s，雷諾數(shù)為1.3×104。計(jì)算收斂殘差設(shè)定為小于10－6。數(shù)值計(jì)算中，先后使用了k－ε、SST、SSG、RNGk－ε、k－ω5種湍流模型，并將其計(jì)算結(jié)果與實(shí)驗(yàn)測(cè)量數(shù)據(jù)進(jìn)行比較。利用ICEMCFD工具劃分網(wǎng)格。為了考察網(wǎng)格對(duì)計(jì)算結(jié)果的影響，選用5種不同疏密的網(wǎng)格分別進(jìn)行計(jì)算，選取的計(jì)算參量是頂部間隙泄漏質(zhì)量流量。從表2可以看出，隨著網(wǎng)格的進(jìn)一步加密，間隙流質(zhì)量流量變化很小，最大僅僅為1.3%，所以可認(rèn)為此處的計(jì)算結(jié)果與網(wǎng)格無(wú)關(guān)。計(jì)算時(shí)選取的網(wǎng)格單元總數(shù)為596萬(wàn)。

表2 網(wǎng)格無(wú)關(guān)性驗(yàn)證Table 2 Validation of mesh independent

計(jì)算區(qū)域分為通道區(qū)和葉頂區(qū)兩部分，兩部分的網(wǎng)格無(wú)關(guān)性驗(yàn)證是同時(shí)進(jìn)行的，在驗(yàn)證網(wǎng)格與計(jì)算結(jié)果無(wú)關(guān)后，通道區(qū)選用結(jié)構(gòu)化六面體網(wǎng)格，網(wǎng)格數(shù)為198萬(wàn)；葉頂區(qū)及氣膜孔內(nèi)選用非結(jié)構(gòu)化四面體網(wǎng)格，葉頂近壁面網(wǎng)格采用三棱柱網(wǎng)格進(jìn)行加密，網(wǎng)格數(shù)為398萬(wàn)；總網(wǎng)格數(shù)約為596萬(wàn)，壁面第一層網(wǎng)格的y＋值控制在1之內(nèi)。圖3展示了計(jì)算區(qū)域的部分網(wǎng)格。

3 結(jié)果與討論

3.1 頂部泄漏渦

為了觀(guān)測(cè)泄漏渦的產(chǎn)生和發(fā)展，實(shí)驗(yàn)選取經(jīng)過(guò)最靠近尾緣的3個(gè)氣膜孔中心且垂直于葉片壁面的截面為測(cè)量面，如圖4所示。

圖4 二維流場(chǎng)拍攝面位置示意圖Fig.4 Schematic of measurement location

為降低測(cè)量結(jié)果的波動(dòng)性，測(cè)量時(shí)每個(gè)截面連續(xù)拍攝150張圖片，處理時(shí)將這150組數(shù)據(jù)做平均，得到準(zhǔn)穩(wěn)態(tài)速度矢量分布。

計(jì)算時(shí)選用了5種不同的湍流模型，最終經(jīng)過(guò)與實(shí)驗(yàn)數(shù)據(jù)的比較發(fā)現(xiàn)：RNGk－ε模型得到的計(jì)算結(jié)果與實(shí)驗(yàn)值最為接近。圖5展示了3個(gè)測(cè)量面上的實(shí)驗(yàn)與計(jì)算得到的速度場(chǎng)分布。其中，左圖為PIV系統(tǒng)捕捉到的泄漏渦，右側(cè)為用RNG k－ε模型計(jì)算得到的數(shù)值模擬結(jié)果。

由圖5可以看出：（1）越接近葉片尾緣，泄漏渦越偏離吸力面，而且渦的范圍越大；（2）無(wú)論是泄漏渦的位置還是速度值，數(shù)值模擬給出的結(jié)果與實(shí)驗(yàn)值都很接近。泄漏渦的產(chǎn)生是由于流體在葉片頂部接近吸力面處具有較大的速度梯度，于是在接近吸力面外側(cè)形成了旋渦。同時(shí)由于頂部泄漏流具有主流方向的分速度，旋渦會(huì)因?yàn)橹髁鞯挠绊懚粩嗟仄x吸力面。

3.2 頂部氣膜孔射流

為了觀(guān)察氣膜孔射流對(duì)頂部泄漏流的影響，在圖4所示的截面2位置，將CCD相機(jī)重新對(duì)焦，把拍攝范圍縮小在頂部間隙區(qū)域，當(dāng)氣膜孔注射流以平均流速0.3m／s進(jìn)入頂部凹槽時(shí)，拍攝了150幅流場(chǎng)矢量圖，在進(jìn)行數(shù)字圖像平均處理后，得到圖6所示速度分布。圖中所示速度矢量只指示速度的方向，其長(zhǎng)短不代表速度大小，速度大小值由上方的顏色標(biāo)尺標(biāo)定。

圖6 頂部氣膜孔對(duì)間隙流的影響（單位：mm）Fig.6 Effect of tip film hole to tip leakage（unit：mm）

由于射流垂直葉頂噴出，為了更直觀(guān)地展示通過(guò)氣膜孔射流對(duì)泄漏流的阻擋作用，圖6展示了二維速度場(chǎng)測(cè)量結(jié)果。從圖6可以看到葉片頂部氣膜孔射流對(duì)泄漏流產(chǎn)生的阻擋作用，即：葉頂氣膜孔附近垂直向上的射流速度分量較大，遠(yuǎn)離氣膜孔這個(gè)速度分量逐漸衰減，而這種向上的射流直接阻礙了橫穿間隙的泄漏流。圖7展示了三維速度數(shù)值計(jì)算結(jié)果，這種阻擋作用在三維合速度分布圖中更加清晰。由圖7可見(jiàn)，垂直向上的射流擠壓了泄漏流通道，在氣膜孔中心的最大速度射流甚至與泄漏流方向相反，因此阻擋效果更加明顯。

圖7 頂部間隙流的數(shù)值模擬結(jié)果Fig.7 Simulation results of tip leakage flow

為了定量地展示頂部氣膜孔射流對(duì)泄漏流的抑制作用，通過(guò)數(shù)值手段分別計(jì)算了不同氣膜孔分布下的泄漏流流量，如表3所示：

表3 不同氣膜孔分布下的泄漏流流量Table 3 Leakage flow rate with different film holes distribution

從表3可以發(fā)現(xiàn)，氣膜孔射流最多能減少6.12%的主流泄漏，說(shuō)明通過(guò)氣膜孔射流來(lái)抑制主流泄漏是一種行之有效的手段。

4 結(jié) 論

用PIV實(shí)驗(yàn)測(cè)量和CFX數(shù)值模擬兩種方法研究了透平葉片頂部泄漏流運(yùn)動(dòng)特性。通過(guò)實(shí)驗(yàn)和計(jì)算結(jié)果比較分析，得出以下結(jié)論：

（1）從所得速度分布圖可以看出泄漏渦的產(chǎn)生和發(fā)展趨勢(shì)：越靠近葉片尾緣，吸力面上產(chǎn)生的泄漏渦越偏離吸力面，而且渦的范圍越來(lái)越大；

（2）通過(guò)與實(shí)驗(yàn)數(shù)據(jù)的對(duì)比發(fā)現(xiàn)，所采用RNG k－ε模型計(jì)算結(jié)果與實(shí)驗(yàn)捕捉的泄漏流特性在速度和變化趨勢(shì)上都很接近，因此RNG k－ε模型可用來(lái)預(yù)估泄漏流特性；

（3）不論是實(shí)驗(yàn)還是數(shù)值模擬，在葉片頂部氣膜孔出口，都能清楚地看到由氣膜孔出來(lái)的射流對(duì)泄漏流的阻擋作用；

（4）該實(shí)驗(yàn)條件下，由于氣膜孔的作用，最大能降低6.12%的泄漏流量損失。

［1］ BOOTH T C，DODGE P R，HEPWORTH H K.Rotor－tip leakage：part 1－basic methodology［J］.ASME J.Eng.Power，1982，104：154－161.

［2］ BUNKER R S.A review of turbine blade tip heat transfer［J］.Annals of the New York Academy of Sciences，2001，934（1）：64－79.

［3］ MAYLE R E，METZGER D E.Heat transfer at the tip of an unshrouded turbine blade［C］／／Proceedings of the 7thinternational heat transfer conference，1982，3：87－92.

［4］ CHYU M K，MOON H K，METZGER D E.Heat transfer in the tip region of a rotor blade simulator［R］.NASA Lewis Research Center，Turbine Engine Hot Section Technology，1986：177－192.

［5］ WANG J H，LIU Y L，WANG X C.Characteristics of tip leakage flow of the turbine blade with cutback squealer and coolant injection［R］.ASME GT－2010－22566，2010.

［6］ PAPA M，GOLDSTEIN R J，GORI F.Effects of tip geometry and tip clearance on the mass／heat transfer from a large－scale gas turbine blade［J］.ASME J.Turbomach.，2003，125（1）：90－96.

［7］ JIN P，GOLDSTEIN R J.Local mass／heat transfer on a turbine blade tip［J］.Int.J.Rotating Mach.，2003，9（2）：81－95.

［8］ JIN P，GOLDSTEIN R J.Local mass／heat transfer on turbine blade near－tip surfaces［J］.ASME J.Turbomach.，2003，125：521－528.

［9］ SRINIVASAN V，GOLDSTEIN R J.Effect of endwall motion on blade tip heat transfer［J］.ASME J.Turbomach.，2003，125：267－273.

［10］AMERI A A，STEINTHORSSON E，Rigby D L.Effect of squealer tip on rotor heat transfer and efficiency［J］.ASME J.Turbomach.，1998，120：753－759.

［11］YANG H，ACHARYA S，EKKAD S V.Numerical simulation of flow and heat transfer past a turbine blade with a squealer－tip［R］.ASME GT－2002－30193，2002.

［12］BRIAN M T，TANG P P，DAVID R H，et al.Computational modelling of tip heat transfer to a superscale model of an unshrouded gas turbine blade［R］.ASME GT－2008－51212，2008.