姜世杰,李志剛,李軍

(西安交通大學能源與動力工程學院,710049,西安)

燃氣渦輪進口溫度和壓力的升高在提高能量轉(zhuǎn)換效率、降低熱耗率的同時,推動了燃氣渦輪向大功率方向發(fā)展。渦輪動葉葉頂間隙的高溫泄漏燃氣與動葉主流通道渦和二次流渦相互摻混,降低了動葉氣動效率,同時導致動葉頂部產(chǎn)生高傳熱系數(shù),具有高熱負荷的動葉頂部傳熱及冷卻問題導致葉頂燒蝕是渦輪動葉失效的主要原因[1]。動葉葉頂和機匣之間存在相對運動,高溫泄漏燃氣在旋轉(zhuǎn)作用下受到離心力、科氏力的影響,葉頂間隙的復雜流動形態(tài)為凹槽葉頂構(gòu)型設(shè)計帶來了挑戰(zhàn)。目前,為減少葉頂間隙泄漏損失、降低葉頂?shù)臒嶝摵?從而提高渦輪動葉的氣熱性能,Bunker等指出,構(gòu)型設(shè)計是提高葉頂氣熱性能的關(guān)鍵技術(shù)之一[2-3]。

科研人員對渦輪動葉葉頂間隙泄漏流和傳熱特性做了大量研究。Kwak等通過實驗測量方法對6種靜止狀態(tài)下不同肩壁結(jié)構(gòu)葉頂壁面換熱系數(shù)進行了研究,結(jié)果表明凹槽狀葉頂明顯降低了葉頂壁面?zhèn)鳠嵯禂?shù),不同葉頂肩壁結(jié)構(gòu)幾何排列改變了泄漏流動路徑,導致了不同傳熱系數(shù)分布[4]。鐘兢軍等基于直列葉柵實驗驗證了葉頂小翼結(jié)構(gòu)能減小泄漏量并提升葉柵氣動性能[5]。張敏等采用數(shù)值方法研究了1.5級動葉葉頂小翼結(jié)構(gòu)對渦輪級氣動性能的影響,結(jié)果表明,葉頂小翼結(jié)構(gòu)能抑制泄漏流與主流的摻混,進而減少泄漏損失[6]。

在葉頂傳熱特性研究方面,Azad等對凹槽狀葉頂泄漏流動和傳熱特性進行了研究,結(jié)果表明,凹槽狀結(jié)構(gòu)比平葉頂?shù)男孤┝坑酗@著的減少并降低了葉頂平均傳熱系數(shù)[7]。楊佃亮等數(shù)值研究了凹槽狀葉頂肩壁位置對泄漏流動和傳熱特性的影響,結(jié)果表明吸力側(cè)肩壁的凹槽葉頂具有最低的平均傳熱系數(shù)[8]。Coull等研究了帶小翼凹槽狀葉頂?shù)膫鳠岷蜌鈩有阅?結(jié)果發(fā)現(xiàn),相比于平葉頂結(jié)構(gòu),帶小翼的葉頂結(jié)構(gòu)可提升46%氣動性能,同時降低了葉頂平均傳熱系數(shù)[9]。Zou等提出了葉頂小翼結(jié)構(gòu)并進行了氣動性能和傳熱特性的研究,指出了葉頂小翼結(jié)構(gòu)的設(shè)計原則[10]。黃琰等數(shù)值研究了靜止狀態(tài)下葉頂壓力側(cè)小翼對其氣熱性能的影響,指出壓力側(cè)小翼在靜止狀態(tài)能減少葉頂泄漏流量并降低葉頂傳熱系數(shù)[11]。Zhong等指出了小翼結(jié)構(gòu)在亞聲速和跨聲速流動條件下均能改善葉頂傳熱性能[12-13]。Zhou等發(fā)展了提高葉頂氣熱性能的小翼結(jié)構(gòu)設(shè)計方法,提高了葉頂?shù)臍鉄嵝阅躘14]。

目前,葉頂小翼氣熱性能研究中,小翼結(jié)構(gòu)設(shè)計采用平面小翼結(jié)構(gòu)來提升葉頂性能,但這些小翼結(jié)構(gòu)沒有使用凹槽設(shè)計。本文基于GE-E3渦輪級動葉凹槽狀葉頂,在平面小翼的基礎(chǔ)上設(shè)計了壓力側(cè)、吸力側(cè)和兩側(cè)設(shè)計凹槽狀小翼3種葉頂小翼結(jié)構(gòu),并采用數(shù)值方法對比分析凹槽狀小翼結(jié)構(gòu)對葉頂泄漏流動和傳熱特性的影響。本文研究工作將為提高渦輪葉頂?shù)臍鉄嵝阅茉O(shè)計提供參考。

1 數(shù)值方法和驗證

(a)GE-E3渦輪級計算域

(b)GE-E3動葉葉頂計算網(wǎng)格圖1 GE-E3渦輪級和動葉葉頂計算網(wǎng)格

GE-E3渦輪級結(jié)構(gòu)[15]和動葉凹槽狀葉頂計算網(wǎng)格如圖1所示,渦輪級主要幾何參數(shù)如表1所示。葉頂間隙厚度和肩壁厚度為動葉葉高的1%,凹槽深度為動葉葉高的2%。采用ANSYS-ICEM生成結(jié)構(gòu)化網(wǎng)格,計算域包含靜葉和動葉流道,流道兩側(cè)采用周期性邊界條件;為防止回流造成影響,動葉出口設(shè)在動葉尾緣1.8倍軸向弦長處。為保證網(wǎng)格質(zhì)量,葉片壁面周圍和葉頂間隙內(nèi)均采用O型網(wǎng)格,壁面網(wǎng)格進行加密處理以保證壁面y+<1.0,數(shù)值計算采用二階迎風格式,動靜交界面為混合平面方法,當計算殘差達到10-4以下時可認為計算收斂。

表1 渦輪級幾何參數(shù)[15]

數(shù)值計算邊界條件[15]如表2所示,渦輪級靜葉進口給定總溫、總壓,動葉出口給定靜壓,壁面采用無滑移等溫條件。傳熱系數(shù)定義為

h=q/(Tw-Tad)=(q1-q2)/(T1-T2)

(1)

式中:q為壁面的傳熱量;Tw為壁面溫度;Tad為靠近壁面的流體溫度。采用文獻[10]的方法,根據(jù)壁面溫度條件下的傳熱量之差來計算傳熱系數(shù),q1、q2是在壁面溫度為T1、T2條件下的壁面?zhèn)鳠崃?T1=490 K,T2=510 K。

表2 渦輪級數(shù)值計算邊界條件[15]

4種湍流模型下葉頂傳熱系數(shù)與實驗數(shù)據(jù)[15]云圖的對比如圖2所示,可知標準k-ω湍流模型的傳熱系數(shù)分布與實驗數(shù)據(jù)吻合。Ledezma等驗證了采用標準k-ω湍流模型研究小翼葉頂氣熱性能的有效性[16-18]。

(a)實驗數(shù)據(jù)[15]

(b)標準k-ε模型 (c)標準k-ω模型

(d)RNGk-ε模型 (e)SSTk-ω模型
圖2 不同湍流模型下葉頂傳熱系數(shù)等值線云圖對比

表3 計算網(wǎng)格無關(guān)性驗證

(a)壓力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)2)

(b)吸力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)3)

(c)壓力側(cè)和吸力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)4)圖3 3種凹槽狀小翼結(jié)構(gòu)尺寸和計算網(wǎng)格

2 葉頂小翼結(jié)構(gòu)

本文在常用平面小翼結(jié)構(gòu)基礎(chǔ)上設(shè)計使用凹槽狀小翼結(jié)構(gòu),3種凹槽狀小翼的幾何尺寸和計算網(wǎng)格如圖3所示,圖中小翼寬度為3 mm,為葉片厚度的0.4倍,小翼上肩壁厚度、高度均和中間肩壁相同。沒有小翼結(jié)構(gòu)的凹槽狀葉頂稱為結(jié)構(gòu)1,壓力側(cè)凹槽狀小翼葉頂稱為結(jié)構(gòu)2,吸力側(cè)凹槽狀小翼葉頂稱為結(jié)構(gòu)3,壓力側(cè)和吸力側(cè)具有凹槽狀小翼葉頂稱為結(jié)構(gòu)4。動葉凹槽狀小翼葉頂?shù)闹饕獛缀螀?shù)如表4所示。

表4 凹槽狀小翼結(jié)構(gòu)幾何尺寸

3 結(jié)果分析

3.1 傳熱特性

凹槽狀小翼結(jié)構(gòu)葉頂壁面?zhèn)鳠嵯禂?shù)等值線云圖和壁面極限流線如圖4所示。由圖4可知:葉頂結(jié)構(gòu)1的凹槽底部再附線周圍出現(xiàn)一條平行于壓力側(cè)肩壁的高傳熱系數(shù)區(qū);凹槽前緣位置由葉頂泄漏流沖擊作用產(chǎn)生另一個高傳熱系數(shù)區(qū);壓力側(cè)和吸力側(cè)凹槽狀小翼結(jié)構(gòu)2、3的凹槽底部高傳熱系數(shù)區(qū)得到有效減弱;壓力側(cè)和吸力側(cè)兩側(cè)凹槽狀小翼結(jié)構(gòu)4的凹槽底部高傳熱系數(shù)區(qū)進一步降低,葉頂前緣附近的高傳熱系數(shù)區(qū)基本消失。

(a)結(jié)構(gòu)1 (b)結(jié)構(gòu)2

(c)結(jié)構(gòu)3 (d)結(jié)構(gòu)4圖4 葉頂壁面?zhèn)鳠嵯禂?shù)等值線云圖和極限流線

4種凹槽狀小翼結(jié)構(gòu)的葉頂平均傳熱系數(shù)的對比如圖5所示,可知相比于結(jié)構(gòu)1,結(jié)構(gòu)2、3、4的葉頂表面平均傳熱系數(shù)分別降低了12.2%、17.1%、19.8%。結(jié)果表明,葉頂吸力側(cè)凹槽狀小翼結(jié)構(gòu)的傳熱系數(shù)比壓力側(cè)凹槽狀小翼結(jié)構(gòu)更小,在葉頂壓力、吸力兩側(cè)同時設(shè)計凹槽狀小翼結(jié)構(gòu)能最大程度降低葉頂熱負荷。

3.2 葉頂間隙內(nèi)流動結(jié)構(gòu)

葉頂間隙泄漏流動形態(tài)決定凹槽壁面?zhèn)鳠嵯禂?shù)分布,4種凹槽狀小翼結(jié)構(gòu)的葉頂間隙泄漏流形態(tài)如圖6所示。對于葉頂結(jié)構(gòu)1,間隙泄漏流在凹槽壓力側(cè)形成從葉頂前緣伸向葉頂尾緣的壓力側(cè)角渦和刮削渦,泄漏流從壓力側(cè)角渦和刮削渦之間通過;在葉頂60%軸向弦長下游部分,由于刮削渦逐漸偏離到吸力側(cè)肩壁,凹槽內(nèi)逐漸由壓力側(cè)角渦主導流動,泄漏流從凹槽上方直接流過葉頂間隙,不再流入凹槽內(nèi)部。3種凹槽狀小翼結(jié)構(gòu)中,泄漏流從壓力側(cè)角渦和刮削渦之間通過,但壓力側(cè)角渦、刮削渦的強度和位置產(chǎn)生了變化,影響凹槽內(nèi)流動形態(tài)。相比于凹槽狀葉頂結(jié)構(gòu)1,葉頂結(jié)構(gòu)2的壓力側(cè)凹槽內(nèi)部形成從凹槽前緣流向凹槽尾緣的壓力側(cè)凹槽渦,該渦在壓力側(cè)凹槽內(nèi)產(chǎn)生堵塞作用,阻止泄漏流進入壓力側(cè)凹槽;泄漏流在掠過壓力側(cè)凹槽上方間隙后流入中間凹槽,壓力側(cè)角渦強度減小,刮削渦起始位置遠離葉片前緣。在葉頂結(jié)構(gòu)3中,由于吸力側(cè)凹槽的影響,壓力側(cè)角渦強度進一步減弱,刮削渦更加靠近壓力側(cè)肩壁并最后從靠近凹槽尾緣位置流出凹槽。吸力側(cè)凹槽內(nèi)流體主要由兩部分組成:從前緣位置直接流入凹槽形成吸力側(cè)凹槽渦;從中間凹槽流出的泄漏流。這兩部分流體在凹槽50%軸向弦長下游位置處相互作用,并形成與中間凹槽相似的迷宮密封流動結(jié)構(gòu)。葉頂結(jié)構(gòu)4的葉頂流動特征兼具結(jié)構(gòu)2、結(jié)構(gòu)3的特點,刮削渦在靠近壓力側(cè)肩壁的同時起始位置遠離葉片前緣,壓力側(cè)角渦的強度減弱。

(b)壓力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)2)

(c)吸力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)3)

(d)壓力側(cè)和吸力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)4)圖6 4種凹槽狀葉頂間隙泄漏流

凹槽狀葉頂結(jié)構(gòu)1的間隙泄漏渦系結(jié)構(gòu)和典型截面二維流動圖譜如圖7所示,渦系結(jié)構(gòu)采用Q判據(jù)進行繪制。凹槽狀葉頂間隙的流動結(jié)構(gòu)主要由壓力側(cè)角渦和刮削渦主導,壓力側(cè)角渦的強度強于刮削渦。在Ⅰ區(qū)域截面,壓力側(cè)角渦和吸力側(cè)角渦主導流動特征,刮削渦的影響基本可以忽略,泄漏流從葉頂肩壁的兩側(cè)流入凹槽內(nèi)并沖擊凹槽底部產(chǎn)生高傳熱區(qū);在Ⅱ區(qū)域截面,葉頂間隙內(nèi)刮削渦逐漸發(fā)展,葉頂泄漏流在壓力側(cè)角渦、刮削渦和吸力側(cè)角渦的共同作用下流動軌跡呈S狀,并沖擊到凹槽底部形成高傳熱系數(shù)條紋;在Ⅲ區(qū)域截面,壓力側(cè)角渦逐漸成為凹槽內(nèi)的主要渦系結(jié)構(gòu),刮削渦逐漸偏離出凹槽并移動到吸力側(cè)肩壁上方,只有少量泄漏流進入凹槽,凹槽底部沒有高換熱系數(shù)區(qū)形成;在Ⅳ區(qū)域截面,壓力側(cè)角渦完全主導葉頂間隙的流動,泄漏流直接從凹槽上方流過葉頂區(qū)域,不再流入凹槽內(nèi)部。

(a)無小翼結(jié)構(gòu)凹槽狀葉頂間隙渦系結(jié)構(gòu)

(b)間隙二維流動圖譜圖7 無小翼結(jié)構(gòu)的凹槽狀葉頂結(jié)構(gòu)1間隙渦系結(jié)構(gòu)和典型截面的二維流動圖譜

(a)結(jié)構(gòu)2渦系結(jié)構(gòu)

(b)結(jié)構(gòu)2二維流動圖譜

(c)結(jié)構(gòu)3渦系結(jié)構(gòu)

(d)結(jié)構(gòu)3二維流動圖譜

(e)結(jié)構(gòu)4渦系結(jié)構(gòu)

(f)結(jié)構(gòu)4二維流動圖譜圖8 3種凹槽狀小翼結(jié)構(gòu)葉頂間隙渦系結(jié)構(gòu)和典型截面的二維流動圖譜

3種凹槽狀小翼結(jié)構(gòu)葉頂間隙渦系結(jié)構(gòu)和典型截面二維流動圖譜如圖8所示,再附線附近的高傳熱區(qū)系數(shù)強弱主要取決于壓力側(cè)角渦強度,起始位置主要取決于刮削渦。相比于凹槽狀葉頂結(jié)構(gòu)1,葉頂結(jié)構(gòu)2的刮削渦起始位置遠離葉片前緣,吸力側(cè)角渦擴大,區(qū)域Ⅰ擴大;在Ⅱ區(qū)域截面,壓力側(cè)角渦的馬赫數(shù)減小,凹槽底面伴隨再附線產(chǎn)生的高傳熱系數(shù)條紋強度隨之降低;區(qū)域Ⅲ、Ⅳ的位置和范圍基本沒有變化。對于葉頂結(jié)構(gòu)3,由于吸力側(cè)凹槽內(nèi)吸力側(cè)凹槽渦的形成,吸力側(cè)角渦的發(fā)展受到抑制,進而導致區(qū)域Ⅰ減小;在Ⅱ區(qū)域截面,壓力側(cè)角渦的馬赫數(shù)相比結(jié)構(gòu)2更低,凹槽底部的高傳熱系數(shù)條紋強度更低;區(qū)域Ⅱ在葉頂尾緣方向范圍擴大,區(qū)域Ⅲ、Ⅳ的范圍被壓縮。對于葉頂結(jié)構(gòu)4,受到壓力側(cè)、吸力側(cè)兩側(cè)凹槽狀小翼的共同影響,葉頂流動特征集合結(jié)構(gòu)2、結(jié)構(gòu)3的特點,區(qū)域Ⅰ的范圍與結(jié)構(gòu)2類似,相比結(jié)構(gòu)1、結(jié)構(gòu)3有所增加;區(qū)域Ⅱ的范圍與結(jié)構(gòu)3類似,在葉頂尾緣方向擴大;在區(qū)域Ⅱ內(nèi),壓力側(cè)角渦的馬赫數(shù)相比結(jié)構(gòu)1、2、3更低,降低了凹槽底面的高傳熱系數(shù)條紋區(qū)域。

3.3 氣動性能

采用葉柵通道總壓損失系數(shù)來表征葉頂?shù)臍鈩有阅?總壓損失系數(shù)為

(2)

式中:Pt,1為動葉進口的面積平均相對總壓;Pt為當?shù)氐南鄬倝?ρ為當?shù)亓黧w密度;V為旋轉(zhuǎn)域內(nèi)的平均速度。

(a)無小翼結(jié)構(gòu)的凹槽狀葉頂結(jié)構(gòu)(結(jié)構(gòu)1)

(b)壓力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)2)

(c)吸力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)3)

(d)壓力側(cè)和吸力側(cè)凹槽狀小翼結(jié)構(gòu)(結(jié)構(gòu)4)圖9 吸力側(cè)葉柵流道內(nèi)總壓損失系數(shù)等值線云圖

凹槽狀小翼結(jié)構(gòu)葉頂吸力側(cè)葉柵流道總壓損失系數(shù)等值線云圖如圖9所示。對葉頂結(jié)構(gòu)1、2,高總壓損失系數(shù)區(qū)域主要分為兩個部分:第1個區(qū)域位于60%～80%葉高位置,由葉片通道渦產(chǎn)生;第2個區(qū)域位于80%～100%葉高位置,由葉頂?shù)男孤u產(chǎn)生。對葉頂結(jié)構(gòu)3、4,由于吸力側(cè)凹槽狀小翼的影響,流道內(nèi)只有一處高總壓損失系數(shù)區(qū)域,位于吸力側(cè)小翼的下方,由葉頂?shù)男孤u產(chǎn)生。吸力側(cè)葉柵通道渦強度大大減弱,和泄漏渦摻混在一起,與文獻[13-14]得到的結(jié)論一致。

在4種葉頂結(jié)構(gòu)中,吸力側(cè)的泄漏渦均明顯存在且位置基本相似:靠近葉片吸力面,與機匣保持平行,在葉高方向沒有位移。4種葉頂結(jié)構(gòu)通道渦的位置和強度具有較大變化:在葉頂結(jié)構(gòu)1、2中,由于沒有吸力側(cè)小翼的影響,流體從間隙分離并形成通道渦,在向下游的流動過程中逐漸向葉根方向移動;在葉頂結(jié)構(gòu)3、4中,由于吸力側(cè)凹槽狀小翼結(jié)構(gòu)的影響,泄漏渦的位置與葉片吸力面之間存在一定距離,當流體從間隙分離時,會受到泄漏渦卷吸作用,流道內(nèi)的通道渦并不明顯。

為定量分析4種凹槽狀小翼葉頂結(jié)構(gòu)的氣動性能,在動葉尾緣下游10%軸向弦長位置給出截面總壓損失系數(shù)等值線云圖和葉頂間隙泄漏流三維流線,如圖10所示。對于葉頂結(jié)構(gòu)1、2,葉頂泄漏渦引起的氣動損失小于通道渦引起的氣動損失,在葉頂結(jié)構(gòu)2中,部分泄漏流直接從動葉尾緣處流出,增加了泄漏渦產(chǎn)生的氣動損失;對于葉頂結(jié)構(gòu)3、4,其通道渦和泄漏渦相互摻混,總壓損失要小于葉頂結(jié)構(gòu)1、2的。

(a)結(jié)構(gòu)1 (b)結(jié)構(gòu)2

(c)結(jié)構(gòu)3 (d)結(jié)構(gòu)4圖10 4種凹槽狀葉頂結(jié)構(gòu)的葉柵出口截面總壓損失系數(shù)等值線云圖和葉頂間隙泄漏流線

圖11 4種凹槽狀葉頂結(jié)構(gòu)葉柵平均總壓損失系數(shù)沿葉高的分布

圖12 4種凹槽狀葉頂結(jié)構(gòu)葉柵平均總壓損失系數(shù)比較

4 結(jié) 論

采用數(shù)值方法研究了壓力側(cè)、吸力側(cè)和兩側(cè)凹槽狀小翼葉頂結(jié)構(gòu)對其傳熱特性和葉柵流道總壓損失的影響,得到如下結(jié)論。

(1)壓力側(cè)凹槽狀小翼結(jié)構(gòu)2、吸力側(cè)凹槽狀小翼結(jié)構(gòu)3和兩側(cè)凹槽狀小翼結(jié)構(gòu)4的葉頂表面平均傳熱系數(shù),相比于無小翼凹槽狀葉頂結(jié)構(gòu)1的平均傳熱系數(shù)分別降低12.2%、17.1%、19.8%;在葉頂?shù)膬蓚?cè)均設(shè)置凹槽小翼能最大程度降低葉頂?shù)臒嶝摵伞?/p>

(2)在葉頂凹槽內(nèi)壓力側(cè)角渦和刮削渦的作用下,泄漏流在凹槽內(nèi)呈S形流動,在凹槽底面沖擊形成高傳熱系數(shù)條紋;在葉頂兩側(cè)設(shè)置凹槽小翼時,泄漏流的流速降低,壓力側(cè)角渦和刮削渦的強度減弱,從而降低了凹槽底面的傳熱系數(shù)。

(3)對比葉頂結(jié)構(gòu)1,壓力側(cè)凹槽狀小翼結(jié)構(gòu)2的平均總壓損失系數(shù)增加8.5%,吸力側(cè)凹槽狀小翼結(jié)構(gòu)3和兩側(cè)凹槽狀小翼結(jié)構(gòu)4的平均總壓損失系數(shù)分別減少8.5%、2.5%;吸力側(cè)凹槽狀小翼結(jié)構(gòu)3能有效減少氣動損失。綜合葉頂傳熱特性和葉柵流道總壓損失,吸力側(cè)凹槽狀小翼結(jié)構(gòu)具有最佳的氣熱性能。