渦輪葉片間隙泄漏流干涉下氣膜冷卻的流動試驗及數(shù)值模擬

摘 要：氣膜冷卻是渦輪葉片的一種重要冷卻手段。為了研究渦輪葉片的間隙泄漏對氣膜冷卻的影響機(jī)制與規(guī)律，采用GE-E3高壓渦輪葉片開展PIV光學(xué)試驗測量和非定常數(shù)值計算研究。分析不同氣膜孔位置、雷諾數(shù)、吹風(fēng)比等主要參數(shù)對氣膜冷卻流動的影響規(guī)律，發(fā)現(xiàn)間隙泄漏流干涉下不同位置氣膜孔射流的流場效果差別較大；主流雷諾數(shù)在一定范圍內(nèi)對氣膜孔射流的影響較?。淮碉L(fēng)比的變化顯著改變了氣膜孔射流的橫向發(fā)展和葉頂氣膜孔射流的強(qiáng)度，使吸力側(cè)泄漏渦存在遠(yuǎn)離葉片表面的趨勢，同時一定程度上抑制了壓力側(cè)的泄漏流發(fā)展。

關(guān)鍵詞：氣膜冷卻；渦輪葉片冷卻；間隙泄漏流；數(shù)值模擬

中圖分類號：TP391.9" 文獻(xiàn)標(biāo)志碼：B" 文章編號：1671-5276（2024）05-0135-05

Experiment and Numerical Simulation of Turbine Blade Film Cooling under Interference of Tip Gap Leakage Flow

Abstract：Film cooling is an important cooling method for turbine blade. In order to study the influence mechanism and law of tip gap leakage flow of turbine blades on film cooling， the GE-E3 high-pressure turbine blade was applied to carry out PIV experimental tests and unsteady numerical calculation investigation. The influence of different film cooling hole position， inlet Reynolds number， blowing ratio and other main parameters on the film cooling flow was analysed. It is found that the flow field effects of the film cooling hole jet at different positions under the interference of the tip gap leakage flow varify greatly; the mainstream Reynolds number has a minor impact on the film cooling jet within a certain range; the variation of blowing ratio significantly changes the lateral development of the film cooling jet and the strength of the film cooling jet at the top of the blade， keeping the leakage vortex on the suction side far away from the blade surface and inhibiting the leakage flow development on the pressure side to a certain extent at the same time.

Keywords：film cooling；cooling for turbine blade；tip leakage flow；numerical simulation

0 引言

航空發(fā)動機(jī)渦輪葉片進(jìn)口溫度在過去幾十年里已經(jīng)升高到1 900K以上，該溫度遠(yuǎn)遠(yuǎn)超過葉片合金材料的耐溫極限。高效的渦輪葉片冷卻是保證渦輪葉片高效安全工作的一個關(guān)鍵技術(shù)。氣膜冷卻是最重要的外部冷卻技術(shù)之一，在渦輪葉片的冷卻設(shè)計中應(yīng)用廣泛。前期研究表明，氣膜冷卻的效率受到眾多因素的影響，包括冷卻流體與高溫來流的動量比值、質(zhì)量比值、溫度比值、主流的壓力梯度、湍流強(qiáng)度等［1-2］。MHETRAS等［3］試驗研究了葉片表面全氣膜冷卻效果，研究了不同氣膜位置、不同氣膜開孔形狀、來流馬赫數(shù)、吹風(fēng)比等參數(shù)對葉片氣膜冷卻的影響，發(fā)現(xiàn)通道中的渦系結(jié)構(gòu)對吸力側(cè)氣膜冷卻效果影響顯著。王昌葉等［4］采用穩(wěn)態(tài)的SST k-ω湍流模型對葉片兩側(cè)氣膜冷卻射流開展了數(shù)值研究，發(fā)現(xiàn)整體上的氣膜覆蓋效率在吸力側(cè)優(yōu)于壓力側(cè)。REZASOLTANI等［5］對某高壓渦輪冷卻開展了實(shí)驗與數(shù)值計算研究工作，采用壓力敏感漆（PSP）技術(shù)和穩(wěn)態(tài)雷諾平均方法（RANS）計算4種葉頂氣膜冷卻的幾何構(gòu)型，并對吹風(fēng)比等重要參數(shù)的影響規(guī)律進(jìn)行了研究。

高壓渦輪在實(shí)際應(yīng)用中葉片葉尖與機(jī)匣之間通常預(yù)留出很小的間隙，以適應(yīng)葉片旋轉(zhuǎn)，防止葉片與機(jī)匣之間的剮蹭。渦輪葉片壓力側(cè)與吸力側(cè)存在較大的壓力梯度，從而在間隙中形成顯著的葉尖泄漏流，會導(dǎo)致較大的氣動損失。葉尖的泄漏流動流場結(jié)構(gòu)復(fù)雜，呈現(xiàn)強(qiáng)烈的三維非定常特性。前期研究發(fā)現(xiàn)，渦輪葉片葉尖附近冷卻設(shè)計的難度較大，主要原因是該區(qū)域的流場結(jié)構(gòu)復(fù)雜，存在間隙泄漏流、氣膜冷卻、通道渦等［6］。因此，針對渦輪葉尖間隙區(qū)域研究氣膜冷卻對于優(yōu)化葉尖冷卻設(shè)計具有重要的意義。

本文在前期研究基礎(chǔ)上采用先進(jìn)的PIV光學(xué)測量手段和高精度的非定常數(shù)值模擬方法對葉尖泄漏流干涉下的氣膜冷卻流場結(jié)構(gòu)開展試驗和數(shù)值計算，并對關(guān)鍵參數(shù)的影響規(guī)律進(jìn)行研究，為渦輪葉尖冷卻設(shè)計提供一定的數(shù)據(jù)支撐。

1 實(shí)驗系統(tǒng)及參數(shù)

實(shí)驗以經(jīng)典的GE-E3高壓渦輪葉片為研究對象，實(shí)驗中氣膜孔的位置分別位于壓力側(cè)、吸力側(cè)和葉頂，如圖1所示。葉片軸向弦長為86.1mm，柵距為91.5mm，葉高為119.5mm，間隙高度為4.5mm。本研究中考慮粒子圖像測速系統(tǒng)（PIV）激光片光厚度以及保障PIV中大量的示蹤粒子可以通過葉尖間隙，選用了3.7%的高度作為葉尖間隙值。

實(shí)驗中分別在葉片壓力側(cè)和吸力側(cè)98%葉高處及葉頂中弧線上設(shè)置氣膜冷卻孔，按等弧長設(shè)置7個氣膜孔（自葉片流向前向后編號為1—7）。葉片內(nèi)部設(shè)有按比例縮放的冷氣進(jìn)氣腔以起到實(shí)驗中穩(wěn)定氣流的作用，保證氣膜孔的出流均勻。

實(shí)驗系統(tǒng)如圖2所示，各個子系統(tǒng)主要由供氣、測量、實(shí)驗段、PIV等組成。其中測量系統(tǒng)包含流量測量、溫度測量、壓力測量等子系統(tǒng)。本研究的PIV系統(tǒng)使用的激光光源為Nd YAG雙脈沖式激光器，粒子發(fā)生器為LaVision氣溶膠發(fā)生器，示蹤粒子為DEHS（癸二酸二辛脂）粒子。

本文重點(diǎn)關(guān)注氣膜冷卻中主次流摻混的流動特性，分別研究流場對氣膜孔位置、主流雷諾數(shù)、吹風(fēng)比3個參數(shù)的變化規(guī)律，試驗工況如表1所示。

2 數(shù)值計算模型

2.1 非定常湍流模型

本文研究中采用的非定常高精度湍流模型為VLES方法（very-large eddy simulation），該模型的計算可靠性和精度在前期研究中得到了廣泛的驗證［7］?；谂ｎD流體基本假設(shè)求解Favre過濾的控制方程。亞格子尺度湍流應(yīng)力的?；赩LES k-ω模型進(jìn)行求解。通過最核心的分辨率控制函數(shù)Fr對湍流黏性系數(shù)μt進(jìn)行重新?；?。分辨率控制函數(shù)Fr的形式為

式中Lc、Li和Lk的物理意義分別是截斷長度尺度、積分長度尺度和最小kolmogorov長度尺度，具體形式為：

式（2）中Cx為重要的模型參數(shù)，在本文研究中取值為0.61。

VLES k-ω模型中湍動能k和比耗散率ω的輸運(yùn)方程以及湍流黏性系數(shù)表達(dá)式為：

式中：d表示物質(zhì)導(dǎo)數(shù)；sub表示亞格子尺度的湍流量。

具體數(shù)值計算中，根據(jù)局部3個湍流尺度之間的關(guān)系，F(xiàn)r的取值在0 ～ 1.0之間自動地平滑取值，由此確定湍流的?；潭?。

2.2 計算模型和網(wǎng)格劃分

計算域采用單葉片通道，如圖3所示，進(jìn)口距離為2倍的弦長，出口為1倍弦長。在葉片表面、間隙內(nèi)部、冷氣通道壁面進(jìn)行網(wǎng)格局部加密。通過多套網(wǎng)格計算結(jié)果對比，發(fā)現(xiàn)網(wǎng)格量達(dá)到320萬之后對計算結(jié)果流場的影響較小，由此確定最終采用320萬的網(wǎng)格進(jìn)行計算。數(shù)值計算參數(shù)設(shè)置與試驗葉片參數(shù)保持完全相同，如表1所示。

計算邊界條件是兩側(cè)為周期性邊界條件，主流來流設(shè)置為速度進(jìn)口，出口設(shè)置為壓力出口，次流采用的邊界條件為質(zhì)量流量進(jìn)口。計算中采用SIMPLEC壓力-速度耦合算法，動量方程采用限定的中心差分格式，時間離散為二階隱式格式。

3 結(jié)果與討論

3.1 實(shí)驗結(jié)果與數(shù)值計算結(jié)果對比

圖4所示為實(shí)驗中PIV拍攝到的兩個葉片截面處通過葉尖區(qū)域的示蹤粒子圖片。可以看出，流體在壓差作用下通過間隙區(qū)域從壓力側(cè)向吸力側(cè)泄漏，從而形成泄漏流，且隨著流動下游發(fā)展。圖中結(jié)果清晰顯示間隙泄漏流與主流存在復(fù)雜的相互作用，流動的非定常特性顯著。氣膜孔位于吸力面和壓力面時，間隙泄漏流特性呈現(xiàn)顯著的差異。

圖5所示是吸力側(cè)開設(shè)氣膜孔、雷諾數(shù)為50 000且吹風(fēng)比為1.0工況下某截面的PIV測量速度分布云圖以及數(shù)值計算結(jié)果。結(jié)果顯示，實(shí)驗測量和數(shù)值計算均揭示出了氣膜孔產(chǎn)生的射流對間隙泄漏流的影響，是氣膜射流使得泄漏渦有所抬升。圖5中虛線框標(biāo)注區(qū)域為氣膜射流在間隙泄漏流的作用下形成的低速滯止區(qū)。該工況下，由于氣膜射流的動量相對于主流的動量較小，在主流與泄漏流的綜合作用下，射流的動量衰減較快，穿透到主流中的距離也較短。氣膜射流迅速衰減后被泄漏渦裹挾。數(shù)值計算結(jié)果與實(shí)驗結(jié)果宏觀上比較吻合，但數(shù)值計算的泄漏渦較為飽滿，且流動表現(xiàn)出顯著的脈動特征。表明實(shí)驗和數(shù)值計算均較好地復(fù)現(xiàn)了間隙泄漏渦與氣膜孔射流之間的相互作用結(jié)果。

3.2 氣膜孔位于吸力面工況

圖6顯示的是吹風(fēng)比M=1.0時不同主流雷諾數(shù)影響下吸力面的氣膜孔射流與間隙泄漏流的干涉作用流場。增大雷諾數(shù)導(dǎo)致渦核有所增大，但結(jié)果并未明顯表明會對葉片吸力面的換熱帶來改善。次流射流干擾了泄漏渦發(fā)展，泄漏渦流線隨著雷諾數(shù)增加向葉頂方向發(fā)展，并由此影響間隙內(nèi)的流動。

圖7給出了雷諾數(shù)Re=50 000時吹風(fēng)比影響下的氣膜孔射流與泄漏流干涉下的流場結(jié)果。吹風(fēng)比的增大意味著次流的動量增大，從而對泄漏流作用更強(qiáng)，使得泄漏渦的大小有所增大且流線有遠(yuǎn)離葉片向葉頂偏轉(zhuǎn)的趨勢，進(jìn)而抑制了間隙流過來的氣流。

3.3 氣膜孔位于壓力面工況

圖8給出了吹風(fēng)比M=1.0時不同雷諾數(shù)下壓力面設(shè)置氣膜孔的流場結(jié)構(gòu)。結(jié)果顯示壓力面主流流速低于吸力側(cè)截面速度。壓力側(cè)的氣膜孔射流由于受到對向泄漏流的裹挾偏轉(zhuǎn)嚴(yán)重。由此泄漏流在間隙區(qū)域內(nèi)流線更加彎曲，在壁面附近形成回流區(qū)。增大雷諾數(shù)導(dǎo)致泄漏渦流線加劇偏轉(zhuǎn)，抑制形成回流區(qū)。

圖9給出了雷諾數(shù)Re=50 000時吹風(fēng)比影響下射流與主流的干涉流場。吹風(fēng)比增大后導(dǎo)致泄漏流在間隙區(qū)域更加穩(wěn)定，但在葉頂壁面附近的回流區(qū)有增大的趨勢。增大吹風(fēng)比意味著動量增大，從而對泄漏流起到抑制作用，減少泄漏流損失。

3.4 氣膜孔位于葉頂工況

圖10給出了吹風(fēng)比M=2.0時葉頂氣膜出流與泄漏流的干涉流場。間隙內(nèi)的射流導(dǎo)致泄漏流向遠(yuǎn)離葉頂方向偏轉(zhuǎn)，相應(yīng)地次流更加貼近葉頂壁面。葉頂射流由于會明顯抑制泄漏流成為控制泄漏量損失的方法之一。

圖11給出了Re=60 000時吹風(fēng)比影響下葉頂噴氣與泄漏流干涉流場。葉頂噴氣作用下葉頂壁面位置會形成較小的回流區(qū)。吹風(fēng)比增大導(dǎo)致對泄漏流阻滯作用增強(qiáng)，從而泄漏流在間隙中遠(yuǎn)離葉頂表面。吹風(fēng)比增加還導(dǎo)致泄漏渦有所增大，形狀由扁長向圓形發(fā)展，從而使得影響吸力面換熱的區(qū)域有所減小。

4 結(jié)語

采用PIV光學(xué)手段和非定常數(shù)值模擬方法對渦輪葉片不同位置開氣膜孔射流與間隙泄漏渦相互作用開展了研究，重點(diǎn)分析了雷諾數(shù)、吹風(fēng)比對主次流相互干涉的影響規(guī)律。得出以下主要結(jié)論。

1）間隙泄漏流影響下，不同位置的氣膜孔出流效果差異較大。吸力面射流導(dǎo)致泄漏渦流線明顯抬升，射流受到泄漏渦干擾在縱向發(fā)展上受到抑制，使得流動更加不穩(wěn)定；壓力面射流受到泄漏流的對向裹挾，從而一定程度上抑制了泄漏流；葉頂射流與泄漏流干擾下將更加貼近葉頂表面，同時泄漏流遠(yuǎn)離葉頂表面。

2）雷諾數(shù)對射流影響相對較小，但使得吸力面渦核有所增大，壓力面流入間隙泄漏流偏轉(zhuǎn)，進(jìn)而抬升葉頂泄漏流；改變吹風(fēng)比使得吸力面和壓力面射流橫向發(fā)展得到一定程度的增強(qiáng)，使吸力面泄漏渦有所遠(yuǎn)離葉片，對壓力面泄漏流起到了抑制作用。葉頂射流的增大抑制了間隙中泄漏流的發(fā)展。

3）實(shí)驗結(jié)果與非定常數(shù)值計算結(jié)果比較吻合，表明數(shù)值計算方法精度較高。所獲得的PIV實(shí)驗和數(shù)值計算數(shù)據(jù)對于理解和優(yōu)化間隙泄漏流及氣膜射流設(shè)計具有一定的參考價值。

參考文獻(xiàn)：

［1］ HAN J C，DUTTA S，EKKAD S. Gas turbine heat transfer and cooling technology［M］. New York：CRC Press，2012.

［2］ 張強(qiáng)，吉洪湖，梁萌，等. 異型擾流-狹縫氣膜層板結(jié)構(gòu)的冷卻特性分析［J］. 機(jī)械制造與自動化，2022，51（4）：44-47.

［3］ MHETRAS S，HAN J C，RUDOLPH R. Effect of flow parameter variations on full coverage film-cooling effectiveness for a gas turbine blade［J］. Journal of Turbomachinery，2012，134（1）：011004-011013.

［4］ 王昌葉，劉釗，葉綠，等. 葉片沖擊/氣膜復(fù)合冷卻的數(shù)值研究［J］. 工程熱物理學(xué)報，2015，36（1）：55-59.

［5］ REZASOLTANI M，LU K，SCHOBEIRI M T，et al. A combined experimental and numerical study of the turbine blade tip film cooling effectiveness under rotation condition［J］. Journal of Turbomachinery，2015，137（5）：051009-051020.

［6］ SUNDEN B，XIE G N. Gas turbine blade tip heat transfer and cooling：a literature survey［J］. Heat Transfer Engineering，2010，31（7）：527-554.

［7］ HAN XS，KRAJNOVIC＇ S. An efficient very large eddy simulation model for simulation of turbulent flow［J］. International Journal for Numerical Methods in Fluids，2013，71（11）：1341-1360.