冷卻孔布置對(duì)透平轉(zhuǎn)靜腔室性能影響的數(shù)值研究

張峰,王新軍,李軍

(西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安)

?

冷卻孔布置對(duì)透平轉(zhuǎn)靜腔室性能影響的數(shù)值研究

張峰,王新軍,李軍

(西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安)

采用SST湍流模型數(shù)值研究了透平第二級(jí)轉(zhuǎn)、靜腔室的流動(dòng)與封嚴(yán)特性,分析了冷卻孔布置對(duì)腔室內(nèi)的流動(dòng)、冷卻效率以及主流燃?xì)馊肭痔匦缘挠绊?。研究表?冷卻孔的位置對(duì)上游腔室內(nèi)的流動(dòng)影響較大,對(duì)下游腔室基本無(wú)影響;冷卻孔距離上游越近,上游腔室的旋流比越大,級(jí)間密封進(jìn)口的旋流比越小,密封進(jìn)、出口壓比越小,相應(yīng)的流過(guò)級(jí)間密封的質(zhì)量流量越小,上游輪緣密封的燃?xì)馊肭至吭叫?上游腔室轉(zhuǎn)、靜壁面的冷卻效率就越大;下游輪緣密封出流的旋流比越小,相應(yīng)的主流通道的流動(dòng)損失越大,但冷卻孔位置對(duì)下游腔室內(nèi)的旋流比和壁面冷卻效率的影響很小;隨著冷卻空氣流量的增加,3種冷卻孔布置下上游腔室壁面冷卻效率的差值減小。

轉(zhuǎn)靜腔室;輪緣密封;傳熱;數(shù)值計(jì)算

隨著燃?xì)廨啓C(jī)技術(shù)的發(fā)展,透平通流設(shè)計(jì)已經(jīng)達(dá)到了很高的水平,而通過(guò)透平二次空氣系統(tǒng)的研究可進(jìn)一步改進(jìn)透平通流設(shè)計(jì)。輪盤(pán)轉(zhuǎn)、靜腔室是二次空氣系統(tǒng)中最復(fù)雜的通流單元,因此對(duì)轉(zhuǎn)、靜腔室性能的研究得到了重視。

轉(zhuǎn)、靜腔室的研究集中于預(yù)旋轉(zhuǎn)、靜腔室,這類腔室主要用來(lái)給第一級(jí)動(dòng)葉提供冷卻空氣。Didenko等數(shù)值研究了腔室寬度和預(yù)旋噴嘴位置對(duì)轉(zhuǎn)、靜腔室的預(yù)旋效率、壓力損失系數(shù)和流量系數(shù)的影響[1]。羅翔等數(shù)值研究了某預(yù)旋轉(zhuǎn)、靜腔室的流動(dòng)與換熱特性[2]。

羅羅公司和薩塞克斯大學(xué)熱流中心對(duì)第二級(jí)轉(zhuǎn)、靜腔室進(jìn)行了一系列的實(shí)驗(yàn)和數(shù)值研究,分析了冷卻空氣流量、冷卻孔角度、冷卻孔數(shù)量和導(dǎo)流板對(duì)腔室內(nèi)的流動(dòng)和冷卻效率的影響,研究了腔室壁面的換熱特性和輪緣密封的封嚴(yán)效率[3-7]。Liu采用共軛換熱的方法數(shù)值研究了不同冷卻空氣流量下腔室內(nèi)的流動(dòng)與換熱特性[8]。Amirante采用流熱固耦合的方法研究了腔室內(nèi)的換熱特性[9]。

上述研究多集中于冷卻孔的角度、數(shù)量和導(dǎo)流板布置的影響分析,而本文采用商用軟件ANSYS-CFX12.1數(shù)值研究了冷卻孔軸向位置對(duì)第二級(jí)轉(zhuǎn)、靜腔室內(nèi)的流動(dòng)、冷卻效率和主流燃?xì)馊肭痔匦缘挠绊憽?/p>

1 計(jì)算模型和數(shù)值方法

1.1 計(jì)算模型

二級(jí)軸流實(shí)驗(yàn)透平原型[10]由薩塞克斯大學(xué)熱流中心提供,實(shí)驗(yàn)透平包含了39個(gè)靜葉和78個(gè)動(dòng)葉,為了簡(jiǎn)化計(jì)算,選取了整體結(jié)構(gòu)的1/39作為本文的研究對(duì)象。計(jì)算模型如圖1a所示,其中包括渦輪主流通道和第二級(jí)轉(zhuǎn)、靜腔室。第二級(jí)轉(zhuǎn)、靜腔室分上游腔室和下游腔室,上游腔室由第一級(jí)動(dòng)葉輪盤(pán)的下游壁面和第二級(jí)靜葉葉根的上游壁面圍成,而下游腔室由第二級(jí)靜葉葉根的下游壁面和第二級(jí)動(dòng)葉輪盤(pán)的上游壁面圍成。轉(zhuǎn)、靜壁面之間必然存在間隙,因此在2個(gè)腔室頂部與主流通道連接處布置了輪緣密封結(jié)構(gòu)。2個(gè)輪緣密封結(jié)構(gòu)分別位于第二級(jí)靜葉的進(jìn)、出口位置,因此這2處的主流氣體之間存在壓差。當(dāng)沒(méi)有為腔室提供冷卻空氣時(shí),在壓降的作用下主流燃?xì)鈴纳嫌屋喚壝芊馓幦肭智皇?并由級(jí)間密封進(jìn)入下游腔室,最后從下游輪緣密封處流入主流燃?xì)庵?。此時(shí),主流燃?xì)馊肭值馁|(zhì)量流量主要由上、下游輪緣密封處的壓降和級(jí)間密封結(jié)構(gòu)決定。為了防止主流燃?xì)馊肭?必須為腔室提供冷卻空氣,而腔室內(nèi)氣流的流動(dòng)情況受到了2股氣流相互作用的影響。模型中靜葉流道為靜止域,動(dòng)葉流道為旋轉(zhuǎn)域,2個(gè)域之間由交界面連接。冷卻孔沿軸向的布置如圖1b所示,分別定義為上游冷卻孔、中間冷卻孔和下游冷卻孔,它們距上游腔室旋轉(zhuǎn)壁面的距離與腔室寬度的比值x/s分別為0.15、0.58和0.86。

(a)腔室及主流通道幾何結(jié)構(gòu)與網(wǎng)格示意

(b)冷卻孔軸向分布圖1 計(jì)算模型示意圖

采用商用軟件ICEM進(jìn)行了結(jié)構(gòu)化網(wǎng)格劃分,其中葉片周圍區(qū)域采用O網(wǎng)格,腔室倒圓角部分采用Y網(wǎng)格,其他位置采用H網(wǎng)格,對(duì)輪緣密封、級(jí)間密封及一些流動(dòng)和結(jié)構(gòu)復(fù)雜區(qū)域的網(wǎng)格進(jìn)行了加密處理。

1.2 數(shù)值方法及驗(yàn)證

采用商用軟件ANSYS-CFX12.1數(shù)值求解了帶有主流通道的第二級(jí)轉(zhuǎn)、靜腔室的三維定常黏性雷諾時(shí)均N-S方程,并進(jìn)行了網(wǎng)格無(wú)關(guān)性和湍流模型驗(yàn)證。

圖2 不同網(wǎng)格數(shù)下腔室內(nèi)虛線所示位置的氣流旋流比沿徑向的變化

冷卻空氣質(zhì)量流量mc=50 g/s時(shí),3種網(wǎng)格數(shù)下腔室內(nèi)虛線所示位置的氣流旋流比β沿徑向的變化如圖2所示?？梢钥闯?在低半徑處250萬(wàn)網(wǎng)格的預(yù)測(cè)值與420萬(wàn)網(wǎng)格的預(yù)測(cè)值非常接近,且最大相對(duì)變化量低于1%,因此在后續(xù)計(jì)算中網(wǎng)格數(shù)為250萬(wàn)。

圖3給出了mc=50 g/s時(shí)轉(zhuǎn)、靜腔室中監(jiān)控點(diǎn)處的溫度分布,計(jì)算條件與文獻(xiàn)[11]的實(shí)驗(yàn)條件相同。由于未考慮葉片及葉根的導(dǎo)熱,3種湍流模型的數(shù)值預(yù)測(cè)值均略低于實(shí)驗(yàn)值。上游腔室監(jiān)控點(diǎn)處的預(yù)測(cè)值與實(shí)驗(yàn)值比較接近,而下游腔室的預(yù)測(cè)值存在一定差異,SST湍流模型的計(jì)算結(jié)果與實(shí)驗(yàn)值吻合最好。因此,后續(xù)計(jì)算均采用SST湍流模型。

圖3 湍流模型的數(shù)值預(yù)測(cè)值與實(shí)驗(yàn)值對(duì)比

1.3 參數(shù)定義

腔室壁面的冷卻效率

(1)

式中:T0,inlet為主流通道進(jìn)口總溫;Tw為壁面溫度;Tc為冷卻工質(zhì)進(jìn)口總溫。

氣流在腔室內(nèi)流動(dòng)的旋流比

β=Vφ/ωr

(2)

式中:Vφ為氣流的周向速度;ω為輪盤(pán)轉(zhuǎn)速;r為半徑。

冷卻空氣及上游腔室輪緣密封、級(jí)間迷宮型密封處的質(zhì)量流量與主流燃?xì)赓|(zhì)量流量的比值分別為Mc、Mrs、Mls,即

(3)

(4)

(5)

式中:mh為主流燃?xì)獾馁|(zhì)量流量;mrs為流過(guò)上游輪緣密封的質(zhì)量流量;mls為流過(guò)級(jí)間迷宮型密封的質(zhì)量流量。

迷宮型密封進(jìn)、出口壓比

式中:P0,ls為迷宮型密封進(jìn)口總壓;Pls,out為迷宮型密封出口靜壓。

1.4 邊界條件

表1給出了計(jì)算時(shí)采用的邊界條件。計(jì)算工質(zhì)為理想氣體,密度為溫度的函數(shù),動(dòng)力黏性系數(shù)采用Sutherland定義;主流通道進(jìn)口給定總溫T0,inlet、總壓P0,inlet和進(jìn)口湍動(dòng)度IT,出口給定靜壓Poutlet;冷卻空氣進(jìn)口給定總溫Tc,冷卻空氣質(zhì)量流量與主流通道質(zhì)量流量的比值Mc(冷卻空氣來(lái)自壓氣機(jī)抽氣,不同的抽汽位置和沿程密封結(jié)構(gòu)等都會(huì)改變冷卻空氣進(jìn)口總壓,從而改變冷卻空氣質(zhì)量流量,一定的進(jìn)口總壓對(duì)應(yīng)著一定的冷卻空氣質(zhì)量流量,此處考慮3種冷卻空氣質(zhì)量流量的影響)和進(jìn)口湍動(dòng)度ITc;動(dòng)葉和旋轉(zhuǎn)輪盤(pán)的轉(zhuǎn)速為ω;轉(zhuǎn)、靜域之間采用交界面連接,周向給定周期性邊界條件,其他壁面為絕熱無(wú)滑移壁面。當(dāng)連續(xù)方程、動(dòng)量方程、能量方程、湍流方程的均方根殘差小于10-6,且位于上游輪緣密封處的監(jiān)控點(diǎn)監(jiān)測(cè)的溫度和質(zhì)量流量曲線平穩(wěn)時(shí)計(jì)算收斂。

表1 邊界條件

2 計(jì)算結(jié)果及分析

2.1 轉(zhuǎn)、靜腔室內(nèi)的流動(dòng)分析

圖4為3種冷卻孔布置下上游腔室內(nèi)氣流流動(dòng)對(duì)比。可以看出,當(dāng)Mc=0.61%時(shí),少量冷卻空氣進(jìn)入腔室。在旋轉(zhuǎn)輪盤(pán)離心力的作用下,入侵的主流燃?xì)庵苯恿飨蜢o止壁面,形成一個(gè)順時(shí)針?lè)较蛐D(zhuǎn)的渦,同時(shí)一個(gè)逆時(shí)針旋轉(zhuǎn)的渦出現(xiàn)在腔室右下角。來(lái)自上游冷卻孔的冷卻空氣受到旋轉(zhuǎn)輪盤(pán)離心力的作用,沿轉(zhuǎn)盤(pán)壁面徑向流出,在腔室頂部與入侵燃?xì)鈸交?見(jiàn)圖4a;來(lái)自中間和下游冷卻孔的冷卻空氣由于動(dòng)量較小,且受到入侵燃?xì)獾淖钃?所以冷卻射流僅僅滲入到靜葉的葉根底部,且被入侵燃?xì)鈮褐圃谇皇业挠蚁陆?其中來(lái)自中間冷卻孔的冷卻空氣進(jìn)入腔室后形成了一個(gè)順時(shí)針旋轉(zhuǎn)的渦,然后向級(jí)間密封處流動(dòng),見(jiàn)圖4b,而來(lái)自下游冷卻孔的冷卻空氣受到葉根底部靶面的阻擋,直接向級(jí)間密封處流動(dòng),見(jiàn)圖4c。

Mc=0.61%(a)上游冷卻孔 (b)中間冷卻孔 (c)下游冷卻孔

Mc=1.22%(d)上游冷卻孔 (e)中間冷卻孔 (f)下游冷卻孔

Mc=1.84%(g)上游冷卻孔 (h)中間冷卻孔 (i)下游冷卻孔圖4 上游腔室的流線分布

隨著冷卻空氣質(zhì)量流量的增加,入侵的燃?xì)饬髁繙p小,當(dāng)Mc=1.22%時(shí),少量燃?xì)馊肭智皇?而當(dāng)Mc=1.84%時(shí),冷卻空氣從上游輪緣密封處溢出,見(jiàn)圖5,因此冷卻空氣對(duì)腔室內(nèi)氣流流動(dòng)的影響增強(qiáng)。比較圖4a、4d、4g可以發(fā)現(xiàn),在上游布置了冷卻孔,腔室內(nèi)氣流流動(dòng)情況基本不受冷卻空氣質(zhì)量流量的影響。這是因?yàn)樯嫌卫鋮s孔緊鄰旋轉(zhuǎn)輪盤(pán)壁面,使得來(lái)自上游冷卻孔的冷卻射流受到壁面旋轉(zhuǎn)離心力的作用較大,氣流始終沿壁面徑向出流。對(duì)于中間冷卻孔,當(dāng)Mc=1.22%時(shí),燃?xì)馊肭至髁繙p小,使得入侵氣流對(duì)冷卻射流的阻擋作用減弱,冷卻射流的動(dòng)量增加,因此來(lái)自中間冷卻孔的冷卻空氣能夠向腔室頂部滲入,形成一個(gè)逆時(shí)針旋轉(zhuǎn)的渦,此外,受轉(zhuǎn)盤(pán)離心力的作用,到達(dá)轉(zhuǎn)盤(pán)的冷卻射流在腔室中心位置形成了一個(gè)順時(shí)針旋轉(zhuǎn)的渦,見(jiàn)圖4e;當(dāng)Mc=1.84%時(shí),冷卻射流的動(dòng)量進(jìn)一步增大,使得位于腔室底部的渦向腔室頂部膨脹,見(jiàn)圖4h。對(duì)于下游冷卻孔,當(dāng)Mc=1.22%時(shí),冷卻射流的動(dòng)量增大,入侵燃?xì)獾淖钃踝饔脺p弱,因此來(lái)自下游冷卻孔的冷卻射流在沖擊靜葉葉根底部的靶面之后,能夠沿軸向向腔室上游滲入,由此形成了一個(gè)逆時(shí)針旋轉(zhuǎn)的渦,見(jiàn)圖4f;當(dāng)Mc=1.84%時(shí),冷卻射流的動(dòng)量進(jìn)一步增強(qiáng),使得渦向腔室上游和頂部膨脹,見(jiàn)圖4i。

圖5 上、下游腔室輪緣密封處凈流量隨冷卻空氣流量的變化

上游腔室的氣流最終通過(guò)級(jí)間密封流向下游腔室,氣流在級(jí)間密封中流動(dòng)時(shí),密封結(jié)構(gòu)能夠消除上游腔室氣流流動(dòng)的差異,使得氣流在下游腔室的流動(dòng)情況相同,如圖6所示。因此,冷卻孔的位置對(duì)下游腔室的流動(dòng)情況基本無(wú)影響,下游腔室的氣流流動(dòng)主要由一個(gè)逆時(shí)針旋轉(zhuǎn)的渦來(lái)控制。

Mc=0.61%(a)上游冷卻孔 (b)中間冷卻孔 (c)下游冷卻孔圖6 下游腔室的流線分布

由分析可知,腔室中氣流的流動(dòng)受到了入侵燃?xì)夂屠鋮s空氣的影響。入侵燃?xì)鈦?lái)自第一級(jí)動(dòng)葉出口,由于動(dòng)葉出口氣流角與轉(zhuǎn)子旋轉(zhuǎn)方向相反,因此入侵燃?xì)獾男鞅群苄?而冷卻孔位于第一級(jí)動(dòng)葉輪盤(pán)與第二級(jí)動(dòng)葉輪盤(pán)之間的驅(qū)動(dòng)臂上,冷卻空氣在冷卻孔中隨輪盤(pán)旋轉(zhuǎn),使得進(jìn)入腔室的冷卻空氣旋流比很大,因此氣流旋流比能夠定性地反映入侵燃?xì)夂屠鋮s空氣的相對(duì)大小。圖7為轉(zhuǎn)、靜腔室內(nèi)氣流旋流比分布云圖?？梢钥闯?隨著冷卻空氣質(zhì)量流量的增加,腔室內(nèi)的旋流比增大。

Mc=0.61%(a)上游冷卻孔 (b)中間冷卻孔 (c)下游冷卻孔

Mc=1.22%(d)上游冷卻孔 (e)中間冷卻孔 (f)下游冷卻孔

Mc=1.84%(g)上游冷卻孔 (h)中間冷卻孔 (i)下游冷卻孔圖7 腔室內(nèi)氣流旋流比分布云圖

2.2 轉(zhuǎn)、靜腔室內(nèi)冷卻效率和密封效率分析

圖8為3種冷卻孔布置下級(jí)間密封進(jìn)口處的氣流旋流比隨冷卻空氣質(zhì)量流量的變化?？梢钥闯?隨著冷卻空氣質(zhì)量流量的增加,氣流的旋流比增大。冷卻孔距離下游越近,氣流的旋流比越大。這是因?yàn)槔鋮s孔越靠近下游,直接進(jìn)入級(jí)間密封的冷卻空氣量就越大,來(lái)自上游冷卻孔的冷卻空氣與入侵的主流燃?xì)鈸交?使得氣流周向速度減小,氣流徑向流過(guò)靜止壁面后才進(jìn)入級(jí)間密封,靜止壁面對(duì)氣流的黏性阻力的作用又使氣流周向速度進(jìn)一步減小,因此來(lái)自上游冷卻孔的冷卻空氣在級(jí)間密封進(jìn)口處的氣流旋流比最小。

圖8 級(jí)間密封進(jìn)口旋流比隨冷卻空氣質(zhì)量流量的變化

圖9為3種冷卻孔布置下級(jí)間密封進(jìn)、出口壓比隨冷卻空氣質(zhì)量流量的變化?？梢钥闯?上游布置冷卻孔后級(jí)間密封進(jìn)、出口壓比最小。這是因?yàn)閬?lái)自上游冷卻孔的冷卻空氣全部與入侵燃?xì)鈸交旌螽a(chǎn)生的壓力損失較大,而來(lái)自中間冷卻孔的部分冷卻空氣直接流向級(jí)間密封,相應(yīng)的冷卻空氣與燃?xì)鈸交旌螽a(chǎn)生的壓力損失較小,因此中間布置冷卻孔的級(jí)間密封進(jìn)、出口壓比大于上游布置冷卻孔;下游冷卻孔的位置正對(duì)靜葉葉根的底部靶面,冷卻射流受葉根底部靶面的限制被阻塞在級(jí)間密封進(jìn)口處,使得進(jìn)口總壓升高,因此下游布置冷卻孔后級(jí)間密封進(jìn)、出口壓比最大。

圖9 級(jí)間密封進(jìn)、出口壓比隨冷卻空氣質(zhì)量流量的變化

級(jí)間迷宮型密封的質(zhì)量流量主要受密封進(jìn)、出口壓比和密封結(jié)構(gòu)的影響。當(dāng)密封結(jié)構(gòu)不變時(shí),進(jìn)、出口壓比越大,流過(guò)級(jí)間密封的質(zhì)量流量越大,因此下游冷卻孔的級(jí)間密封質(zhì)量流量最大,中間冷卻孔次之,上游冷卻孔最小,如圖10所示。流過(guò)級(jí)間密封的氣流由冷卻空氣和入侵燃?xì)饨M成。當(dāng)冷卻空氣供給量一定時(shí),增加的級(jí)間密封質(zhì)量流量由入侵燃?xì)馓钛a(bǔ),因此冷卻孔距離下游越近,燃?xì)獾娜肭至吭酱?如圖11所示。

圖10 3種冷卻孔布置下級(jí)間密封質(zhì)量流量對(duì)比

圖11 3種冷卻孔布置下燃?xì)馊肭至髁繉?duì)比

圖12 3種冷卻孔布置的上游腔室壁面平均冷卻效率對(duì)比

圖13 3種冷卻孔布置的下游腔室壁面平均冷卻效率對(duì)比

圖12和圖13分別給出了上游腔室和下游腔室的轉(zhuǎn)、靜壁面平均冷卻效率的對(duì)比?？梢钥闯?冷卻空氣流量越大,腔室壁面的冷卻效率越大;旋轉(zhuǎn)壁面的冷卻效率大于靜止壁面,但下游腔室轉(zhuǎn)、靜壁面的冷卻效率相差很小;3種冷卻孔布置對(duì)上游腔室壁面的冷卻效率影響明顯,而對(duì)下游腔室基本無(wú)影響。當(dāng)Mc=0.61%、1.22%時(shí),冷卻孔距離上游越近,上游腔室壁面的冷卻效率越大。這是因?yàn)閬?lái)自上游冷卻孔的冷卻空氣受到轉(zhuǎn)盤(pán)的泵效應(yīng)影響比較顯著,使得冷卻射流沿輪盤(pán)壁面徑向出流,最終沿靜止壁面回流,這部分冷卻空氣能夠減少入侵燃?xì)馀c壁面的直接接觸,降低燃?xì)鉁囟?并對(duì)輪盤(pán)和葉根壁面進(jìn)行冷卻保護(hù),且由上游冷卻孔布置引起的燃?xì)馊肭至枯^小,因此冷卻效率較高;來(lái)自中間和下游冷卻孔的氣流被壓制在腔室右下角和腔室底部,使得輪盤(pán)和葉根的大部分壁面直接暴露在入侵燃?xì)庵?。此?2種冷卻孔布置的燃?xì)馊肭至枯^大,且部分冷卻空氣直接流向級(jí)間密封,沒(méi)有起到冷卻上游腔室壁面的作用,因此冷卻效率較低。當(dāng)Mc=1.84%時(shí),冷卻空氣充滿了整個(gè)上游腔室,并從上游輪緣密封處溢出,此時(shí)3種冷卻孔布置下的壁面冷卻效率基本相同,但中間冷卻孔布置的壁面冷卻效率稍高。這是因?yàn)閬?lái)自下游冷卻孔的氣流直接流向級(jí)間密封的質(zhì)量流量較大,而來(lái)自上游冷卻孔的氣流受到的轉(zhuǎn)盤(pán)黏性加熱作用明顯。此外,隨著冷卻空氣質(zhì)量流量的增加,3種冷卻孔布置下上游腔室壁面冷卻效率的差值減小。

腔室內(nèi)的氣流最終從下游輪緣密封處流出且進(jìn)入主流與燃?xì)鈸交?。下游輪緣密封位于第二?jí)靜葉出口,靜葉出口氣流角與轉(zhuǎn)子旋轉(zhuǎn)方向相同,因此氣流的旋流比較大,大于腔室出流的旋流比。圖14給出了下游腔室出口處氣流旋流比隨冷卻空氣質(zhì)量流量的變化。可以看出,冷卻孔距離上游越近,腔室出流旋流比越小。Demargne等的研究指出,腔室出流和通道主流的旋流比相差越大,摻混產(chǎn)生的氣動(dòng)損失就越大[12]。因此,上游布置冷卻孔后流動(dòng)損失最大,中間布置冷卻孔后流動(dòng)損失次之,下游布置冷卻孔后流動(dòng)損失最小。

圖14 下游腔室出口處氣流旋流比隨冷卻空氣流量的變化

3 結(jié) 論

本文數(shù)值研究了透平第二級(jí)轉(zhuǎn)、靜腔室的流動(dòng)和封嚴(yán)特性,分析了不同位置冷卻孔對(duì)腔室內(nèi)的流動(dòng)、冷卻效率及主流燃?xì)馊肭痔匦缘挠绊懸?guī)律,結(jié)論如下。

(1)腔室內(nèi)氣流的流動(dòng)受冷卻空氣和入侵燃?xì)獾挠绊?冷卻孔的位置對(duì)上游腔室內(nèi)的流動(dòng)影響較大,對(duì)下游腔室基本無(wú)影響。

(2)冷卻孔距離上游越近,上游腔室的旋流比越大,級(jí)間密封進(jìn)口的旋流比越小,密封進(jìn)、出口壓比越小,相應(yīng)的流過(guò)級(jí)間密封的質(zhì)量流量越小,上游輪緣密封的燃?xì)馊肭至烤驮叫?上游腔室轉(zhuǎn)、靜壁面的冷卻效率越大,隨著冷卻空氣流量的增加,3種冷卻孔布置下壁面冷卻效率差值減小;下游腔室的旋流比和壁面的冷卻效率基本不受冷卻孔位置的影響。

(3)冷卻孔距離上游越近,下游輪緣密封出流的旋流比越小,相應(yīng)的主流通道的流動(dòng)損失越大。

(4)綜合考慮冷卻孔位置影響可知,冷卻孔靠近上游布置時(shí)腔室的性能最好。

[1] DIDENKO R A, KARELIN D V, IEVLEV D G, et al. Pre-swirl cooling air delivery system performance study [C]∥ASME 2012 Turbine Technical Conference and Exposition. New York, USA: ASME, 2012: 1921-1932.

[2] 羅翔, 馮野, 徐國(guó)強(qiáng), 等. 直接供氣預(yù)旋轉(zhuǎn)靜系流動(dòng)和換熱數(shù)值模擬 [J]. 航空動(dòng)力學(xué)報(bào), 2012, 27(10): 2188-2193. LUO Xiang, FENG Ye, XU Guoqiang, et al. Numerical simulation of flow and heat transfer performances in a direct transfer pre-swirl system [J]. Journal of Aerospace Power, 2012, 27(10): 2188-2193.

[3] DIXON J A, BRUNTON I L, SCANLON T J, et al. Turbine stator well heat transfer and cooling flow optimisation [C]∥ASME Turbo Expo 2006: Power for Land, Sea, and Air. New York, USA: ASME, 2006: 1375-1383.

[4] DIXON J A, VALENCIA A G, COREN D, et al. Main annulus gas path interactions: turbine stator well heat transfer [C]∥ASME 2012 Turbine Technical Conference and Exposition.New York,USA: ASME, 2012: 2013-2024.

[5] VALENCIA A G, DIXON J A, DA SOGHE R, et al. An investigation into numerical analysis alternatives for predicting re-ingestion in turbine disc rim cavities [C]∥ASME 2012 Turbine Technical Conference and Exposition. New York,USA: ASME, 2012: 2025-2035.

[6] VALENCIA A G, DIXON J A, GUARDINI A, et al. Heat transfer in turbine hub cavities adjacent to the main gas path including FE-CFD coupled thermal analysis [C]∥ASME 2011 Turbine Technical Conference and Exposition. New York, USA: ASME, 2011: 833-843.

[7] EASTWOOD D. Investigation of rim seal exchange and coolant re-ingestion in rotor stator cavities using gas concentration techniques [D]. Brighton, UK: University of Sussex, 2014.

[8] LIU H, AN Y, ZOU Z. Aerothermal analysis of a turbine with rim seal cavity [C]∥ASME 2014 Turbine Technical Conference and Exposition. New York, USA: ASME, 2014: V05AT11A006.

[9] AMIRANTE D, HILLS N J, BARNES C J. Thermo-mechanical finite element analysis/computational fluid dynamics coupling of an interstage seal cavity using torsional spring analogy [J]. ASME Journal of Turbomachinery, 2012, 134(5): 051015.

[10]STEFANIS V. Investigation of flow and heat transfer in stator well cavities of a two-stage axial turbine [D]. Brighton, UK: University of Sussex, 2007.

[11]ANDREINI A, DA SOGHE R, FACCHINI B. Turbine stator well CFD studies: effects of coolant supply geometry on cavity sealing performance [J]. ASME Journal of Turbomachinery, 2011, 133(2): 021008.

[12]DEMARGNE A A J, LONGLEY J P. The aerodynamic interaction of stator shroud leakage and mainstream flows in compressors, 2000-GT-0570 [R]. New York, USA: ASME, 2000.

[本刊相關(guān)文獻(xiàn)鏈接]

張峰,王新軍,李軍.球凹平板沖擊冷卻性能的數(shù)值研究及結(jié)構(gòu)改進(jìn).2016,50(1):124-130.[doi:10.7652/xjtuxb201601 019]

杜長(zhǎng)河,范小軍,李亮,等.噴嘴長(zhǎng)寬比和雷諾數(shù)對(duì)旋流冷卻流動(dòng)與傳熱特性的影響.2015,49(12):124-129.[doi:10.7652/xjtuxb201512020]

張峰,王新軍,李軍.帶肋U型通道中的汽霧/空氣流動(dòng)與換熱數(shù)值研究.2015,49(9):52-57.[doi:10.7652/xjtuxb201509 010]

李森,李亮,杜長(zhǎng)河.球窩/球凸結(jié)構(gòu)下的U型通道蒸汽冷卻性能.2015,49(9):63-69.[doi:10.7652/xjtuxb201509012]

王新軍,周子杰,宋釗,等.核電站凝汽器的壓力瞬態(tài)變化特性.2015,49(7):6-10.[doi:10.7652/xjtuxb201507002]

謝金偉,王新軍,周駿飛.空心靜葉汽膜孔排吹掃除濕的數(shù)值研究.2015,49(7):61-66.[doi:10.7652/xjtuxb201507011]

費(fèi)繼友,曾俊雄,高鐵瑜,等.渦旋核心分布對(duì)斜置肋片蒸汽冷卻通道換熱特性的影響.2015,49(5):7-15.[doi:10.7652/xjtuxb201505002]

霍文浩,李軍,鐘剛云,等.冷卻結(jié)構(gòu)對(duì)中壓透平級(jí)蒸汽冷卻性能的影響.2015,49(5):36-42.[doi:10.7652/xjtuxb201505 006]

鐘兵,洪偉,蘇巖,等.點(diǎn)火時(shí)刻對(duì)怠速工況缸內(nèi)直噴汽油機(jī)微粒排放特性的影響.2015,49(3):32-37.[doi:10.7652/xjtuxb201503006]

(編輯 苗凌)

Numerical Investigation for Effect of Coolant Inlet Arrangements on Performance of Turbine Rotor-Stator Cavity

ZHANG Feng,WANG Xinjun,LI Jun

(School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

The flow and sealing performance of secondary rotor-stator cavity was numerically investigated with SST turbulence model. The effect of coolant inlet position on cavity flow, cooling performance and mainstream ingestion phenomenon were analyzed. The results show that the influence of coolant inlet position on the flow in upstream cavity is significant, while that in downstream cavity is slight. As the coolant inlet position moves toward upstream, the swirl ratio in the upstream cavity increases, the swirl ratio at the inlet of interstage labyrinth seal decreases, the pressure ratio and mass flow through interstage labyrinth seal decrease, the mass flow of mainstream ingestion on the upstream rim seal decreases, and meanwhile the cooling performance of the upstream cavity walls increases; whilst the swirl ratio on the downstream rim seal decreases, correspondingly the flow losses of mainstream flow increase, while the swirl ratio and walls cooling performance in the downstream cavity slightly depends on the coolant inlet position. As the mass flow rate of cooling air increases, the difference in cooling performance of the upstream cavity walls set in different coolant inlet positions decreases.

rotor-stator cavity; rim seal; heat transfer; numerical calculation

10.7652/xjtuxb201603009

2015-06-19。 作者簡(jiǎn)介:張峰(1993—),男,博士生;王新軍(通信作者),男,副教授。

時(shí)間:2015-12-10

http:∥www.cnki.net/kcms/detail/61.1069.T.20151210.1128.004.html

TK263

:A

:0253-987X(2016)03-0055-07