補汽對透平級氣動性能和靜葉汽封轉(zhuǎn)子動力特性影響的數(shù)值模擬

王世柱,李志剛,李軍,冀大偉,肖高繪,葛慶

(1.西安交通大學葉輪機械研究所,710049,西安;2.上海電氣電站設(shè)備有限公司,200240,上海)

?

補汽對透平級氣動性能和靜葉汽封轉(zhuǎn)子動力特性影響的數(shù)值模擬

王世柱1,李志剛1,李軍1,冀大偉2,肖高繪2,葛慶2

(1.西安交通大學葉輪機械研究所,710049,西安;2.上海電氣電站設(shè)備有限公司,200240,上海)

為了研究補汽對汽輪機通流部分氣動性能和轉(zhuǎn)子動力特性的影響,首先采用數(shù)值求解三維RANS方程和k-ε湍流模型的方法,研究了無補汽和2種補汽工況下補汽對透平級氣動性能的影響。在此基礎(chǔ)上,以不同補汽工況獲得的透平級靜葉汽封的進出口氣動參數(shù)為邊界條件,采用多頻橢圓渦動模型和動網(wǎng)格技術(shù)計算了靜葉汽封的非定常氣流激振轉(zhuǎn)子動力特性系數(shù)。研究結(jié)果表明:補汽射流沖擊下游靜葉柵的葉根區(qū)域會造成旋渦流動,導致下游靜葉柵的進口氣動參數(shù)沿周向出現(xiàn)不均勻分布,透平級的氣動性能隨著補汽量的增加而下降;補汽射流改變了的下游靜葉汽封的氣動參數(shù),使得預旋比由0.02逐漸增大為0.06、0.16,進而導致汽封的轉(zhuǎn)子動力特性發(fā)生改變。該結(jié)果可為補汽對機組氣動性能和轉(zhuǎn)子動力特性影響的研究提供參考。

透平級;靜葉汽封;氣動性能;轉(zhuǎn)子動力特性;數(shù)值模擬

目前,超超臨界百萬等級汽輪機一般采用補汽來提高機組的輸出功率,采用抽汽回熱方式提高機組循環(huán)效率[1]。有關(guān)抽汽對汽輪機通流效率的影響,國內(nèi)外學者采用實驗測量和數(shù)值模擬的方法進行了研究。毛研偉等通過數(shù)值研究指出,汽輪機抽汽改變了抽汽口上下游氣動參數(shù),進而改變了汽輪機的通流效率[2]。高俊鵬等計算了抽汽對汽輪機下游透平級氣動性能的影響[3]。Rosic等揭示了抽汽導致抽汽口上下游流場的不均勻性,抽汽量可改變下游靜葉柵進口汽流角和損失的分布,由此提出了改進的抽汽結(jié)構(gòu),以提高汽輪機的通流效率[4]。節(jié)流配氣大功率汽輪機在增加機組出力時采用了補汽結(jié)構(gòu),且通過增加流量來提高機組輸出功率[1]?；粑暮频炔捎脭?shù)值模擬方法研究了補汽結(jié)構(gòu)對高壓缸氣動性能的影響[5]。補汽結(jié)構(gòu)會影響補汽口下游靜葉柵的進氣參數(shù)在周向的分布,可以改變下游靜葉柵汽封處的進口預旋,從而導致下游靜葉柵汽封內(nèi)非定常氣流激振的轉(zhuǎn)子動力特性發(fā)生變化,使得機組運行的穩(wěn)定性受到影響。李志剛等實驗研究了迷宮密封的泄漏特性及其影響機制[6]。Yan等針對汽封的偏心渦動導致的非定常氣流激振特性的變化,發(fā)展了三維數(shù)值方法,以計算分析汽封系統(tǒng)的氣流激振和轉(zhuǎn)子動力特性系數(shù)的變化規(guī)律[7-8]。

本文采用數(shù)值方法研究了補汽對上下游透平級氣動性能的影響,

在此基礎(chǔ)上進一步研究了補汽對

下游靜葉汽封非定常氣流激振轉(zhuǎn)子動力特性系數(shù)的影響規(guī)律,以期為汽輪機的高效安全運行和補汽結(jié)構(gòu)設(shè)計提供參考。

1 計算模型與數(shù)值方法

1.1 計算模型

補汽結(jié)構(gòu)在機組通流的第5級和第6級之間水平布置。補汽結(jié)構(gòu)包括主進氣管路、混合腔室和10個進氣管路。補汽量的大小由調(diào)節(jié)閥在主進氣管路進口處進行調(diào)節(jié)。

1.2 數(shù)值方法

圖1為補汽腔室和前后兩級葉柵通道的通流模型和計算網(wǎng)格。葉柵通道為多塊結(jié)構(gòu)化網(wǎng)格,由IGG/AutoGrid軟件生成;補汽結(jié)構(gòu)為非結(jié)構(gòu)網(wǎng)格,由ANSYS-ICEM生成。采用六面體網(wǎng)格和三棱柱網(wǎng)格對壁面進行局部加密處理,以保證k-ε湍流模型的y+>11.06。表1給出了計算區(qū)域的網(wǎng)格數(shù)。

圖1 補汽腔室和前后兩級葉柵通道的通流模型和計算網(wǎng)格

計算域補汽結(jié)構(gòu)第5級靜葉第5級動葉第6級靜葉第6級動葉第5級靜葉汽封第5級動葉汽封第6級靜葉汽封第6級動葉汽封計算區(qū)域網(wǎng)格數(shù)/10417703868312029293841344643683683435731382

采用商用CFD軟件ANSYS-CFX數(shù)值求解了三維RANS方程,求解方法為時間推進法,對流項的離散采用高精度格式(二階),工質(zhì)采用真實水蒸氣模型中的Steam.3vl,該模型是基于IAPWS IF97工業(yè)標準的。湍流模型為標準k-ε兩方程模型。進口給定總溫、總壓和速度方向,出口給定平均靜壓。為了保證動、靜葉交接面的速度場連續(xù),動、靜葉交接面的處理采用“Stage”方法,其中動葉通道、葉頂密封、隔板密封區(qū)域設(shè)為旋轉(zhuǎn)域,靜止域與旋轉(zhuǎn)域之間的交接面設(shè)為“Stage”。

圖2為用于非定常轉(zhuǎn)子動力特性計算的補汽管道下游第6級靜葉汽封剖面圖,表2給出了第6級靜葉汽封的幾何尺寸。圖3為靜葉汽封質(zhì)量流量與周向節(jié)點數(shù)的關(guān)系。可以看出,當周向節(jié)點數(shù)由360增加到480時,計算得到的質(zhì)量流量差別很小,對此迷宮密封而言,周向節(jié)點數(shù)為480時已滿足網(wǎng)格無關(guān)性要求。

圖2 第6級靜葉汽封剖面圖

表2 第6級靜葉汽封幾何參數(shù)

圖3 靜葉汽封質(zhì)量流量與周向節(jié)點數(shù)的關(guān)系

(1)

-Fx=(Kxx+jΩCxx)Dx+(Kxy+jΩCxy)Dy

(2)

-Fy=(Kyy+jΩCyy)Dy+(Kyx+jΩCyx)Dx

(3)

式中:Dx和Dy為轉(zhuǎn)子在x和y兩個方向渦動位移的頻域信號。

多頻橢圓渦動模型假設(shè)轉(zhuǎn)子在單頻下做橢圓渦動(見圖4),利用FFT將時域信號轉(zhuǎn)化為頻域信號,可以從一次計算結(jié)果中提取出多種渦動頻率下的轉(zhuǎn)子動力特性系數(shù),由此大大加快了計算效率。采用ANSYS-CFX數(shù)值求解了非定常RANS方程和k-ε湍流模型,數(shù)值模擬了不同預旋比和渦動頻率下的密封內(nèi)非定常流場,計算了在不同預旋比和頻率下密封的轉(zhuǎn)子動力特性系數(shù)。

迷宮密封的實驗研究[10]表明:迷宮密封一般具有負的阻尼系數(shù),不利于削弱振動,影響轉(zhuǎn)子穩(wěn)定;環(huán)形腔室內(nèi)的周向旋流使得交叉剛度增大,轉(zhuǎn)子振動趨勢增大。交叉剛度是轉(zhuǎn)子非同步低頻渦動激振力的來源,直接阻尼則表示系統(tǒng)阻尼對轉(zhuǎn)子渦動的抑制作用,2個參數(shù)對氣流激振力的影響較大。因此,定義有效阻尼以考慮二者的綜合作用,即

Ceff=Cxx-Kxy/Ω

(4)

式中:Cxx為直接阻尼;Kxy為交叉剛度。

圖4 汽封轉(zhuǎn)子非穩(wěn)態(tài)橢圓渦動模型[8]

在24核80 GB內(nèi)存的工作站上,對帶有補汽結(jié)構(gòu)的全流場進行了定常計算,收斂周期為60～80 h(不同的邊界條件和初始條件會引起收斂時間出現(xiàn)差異)。非定常計算時汽封入口給定總溫、總壓,出口給定靜壓,具體計算邊界條件見表2。計算中,當連續(xù)方程、動量方程、能量方程和湍流方程的均方根殘差小于10-5,相鄰周期轉(zhuǎn)子面受力Fx和Fy周期性振蕩,且相鄰周期受力波動小于0.2%時,認為非定常計算收斂。

2 結(jié)果分析

2.1 補汽對透平級氣動性能的影響

表3給出了不同補汽率下第5、第6級相對總總等熵效率(以第5級0補汽的總總等熵效率為1)的比較。總總等熵效率

(5)

表3 不同補汽率下第5、第6級相對

圖5為補汽率=23%時通過第6級靜葉的72個葉片通道的質(zhì)量流量分布?？梢钥闯?未補汽側(cè)(垂直方向)質(zhì)量流量分布均勻,通過各個葉柵通道的質(zhì)量流量基本相等;補汽側(cè)(水平方向)的工質(zhì)分布沿周向差異較大,呈鋸齒狀,這種周向不均勻可能會引起一部分葉柵通道內(nèi)無足夠的蒸汽流過,葉片做功能力不能充分發(fā)揮,另一部分葉柵通道內(nèi)流過的工質(zhì)過多,導致壅塞。

內(nèi)圓周數(shù)字為葉片通道編號圖5 第6級靜葉質(zhì)量流量的周向分布

圖6為補汽口下游第6級靜葉進口速度分布云圖。靜葉進口速度沿周向也存在不均勻性。補汽管道噴出的蒸汽會撞擊葉柵通道的葉根面,形成局部滯止區(qū)域。補汽管附近的速度矢量分布雜亂,且有很多旋渦存在,隨著工質(zhì)向未補汽區(qū)域的轉(zhuǎn)移,流動逐漸趨于平穩(wěn)。

圖6 第6級靜葉進口速度分布

圖7為第6級靜葉進口截面(0-0)、靜葉與動葉之間截面(1-1)和第6級動葉出口截面(2-2)的溫度分布。比較不同補汽工況可以發(fā)現(xiàn):隨著補汽量的增加,第6級葉柵進口處溫度的周向不均勻性逐漸增大,每個與補汽管臨近的葉柵前方溫度較高,由此導致一部分葉柵通道焓降大,另一部分進口溫度低的葉柵通道焓降小,這種周向分布的不均勻性隨著工質(zhì)在葉柵通道內(nèi)的流動而逐漸減弱。定義工質(zhì)周向速度與軸向速度的比值為氣流角的正切值。

(a)補汽率=0%

(b)補汽率=9%

(c)補汽率=23%圖7 補汽口下游第6級透平級3個截面的溫度分布云圖

圖8為第6級靜葉進口氣流角的分布,圖中給出了幾處氣流角惡化的分布?？梢钥闯?受補汽的影響,補汽管處的工質(zhì)流速的周向分量較大,部分工質(zhì)甚至直接沿周向流動,這會導致工質(zhì)垂直沖擊葉片表面,造成葉柵損失增大、效率降低。

圖8 第6級靜葉進口汽流角分布云圖

2.2 補汽對隔板汽封轉(zhuǎn)子動力特性的影響

表4給出了不同補汽獲得的第6級靜葉汽封的進出口邊界條件(周向質(zhì)量平均)。補汽率的增加對密封進口預旋比影響較大,而進口預旋比λ影響密封腔室內(nèi)的周向流動,引起密封轉(zhuǎn)子動力特性發(fā)生變化。

表4 第6級靜葉汽封進出口邊界條件

預旋比

(6)

圖9為3種補汽率下第6級靜葉汽封的直接剛度Kxx、Kyy隨轉(zhuǎn)子渦動頻率的變化?？梢钥闯?Kxx與Kyy近似相等;3種補汽率下直接剛度均為正值,且受渦動頻率影響不大;過載補汽對直接剛度有明顯影響,隨著補汽量的增加,直接剛度不斷增大。

圖9 第6級靜葉汽封的直接剛度隨渦動頻率的變化

圖10為3種補汽率下第6級靜葉汽封的交叉剛度Kxy、Kyx隨轉(zhuǎn)子渦動頻率的變化。可以看出,每種工況下交叉剛度的絕對值近似相等,但值的符號發(fā)生了轉(zhuǎn)變。

圖10 第6級靜葉汽封交叉剛度隨渦動頻率的變化

隨著補汽率的增加,Kxy由未補汽時的正值逐漸減小為負值,而Kyx由未補汽時的負值變?yōu)榱苏?。對于?shù)值計算給定的轉(zhuǎn)子渦動軌跡,Kxy<0、Kyx>0時,流體激振力對轉(zhuǎn)子產(chǎn)生的阻尼為正,有利于轉(zhuǎn)子的穩(wěn)定;相反,當Kxy>0、Kyx<0時,阻尼為負,不利于轉(zhuǎn)子的穩(wěn)定[9]。

圖11為3種補汽率下第6級靜葉汽封的直接阻尼Cxx、Cyy隨轉(zhuǎn)子渦動頻率的變化。可以看出:3種工況下Cxx=Cyy,且隨渦動頻率的增加不斷增大;在90 Hz附近,直接阻尼的值由負變正,正阻尼有利于轉(zhuǎn)子的穩(wěn)定;在所計算的頻域內(nèi),補汽量越大,直接阻尼越小,從阻尼方面考慮,補汽降低了轉(zhuǎn)子的穩(wěn)定性;補汽引起的直接阻尼之間的差異隨著渦動頻率的增加而減小。

圖11 第6級靜葉汽封直接阻尼隨渦動頻率的變化

圖12為3種補汽率下第6級靜葉汽封的交叉阻尼Cxy、Cyx隨轉(zhuǎn)子渦動頻率的變化。可以看出:3種工況下的交叉阻尼系數(shù)近似相等;交叉阻尼的絕對值隨渦動頻率的增加而不斷減小;隨著補汽量的增加,交叉阻尼的符號發(fā)生了變化,Cxy從未補汽時的正值減小為負值,Cyx從未補汽時的負值增大為正值;補汽引起的交叉阻尼之間的差異隨著渦動頻率的增加而減小。

圖12 第6級靜葉汽封交叉阻尼隨渦動頻率的變化

圖13為3種補汽率下第6級靜葉汽封的有效阻尼Ce隨轉(zhuǎn)子渦動頻率的變化。對于低頻區(qū),補汽量越大,同頻率下的有效阻尼越小,即補汽會降低轉(zhuǎn)子運行的穩(wěn)定性;隨著渦動頻率的增加,補汽引起的有效阻尼之間的差異逐漸減小;在90 Hz時,有效阻尼由負變正,正阻尼有利于轉(zhuǎn)子的穩(wěn)定。

圖13 第6級靜葉汽封有效阻尼隨渦動頻率的變化

3 結(jié) 論

本文采用數(shù)值模擬的方法研究了補汽對大型汽輪機上下游透平級氣動性能和下游靜葉汽封非定常氣流激振轉(zhuǎn)子動力特性的影響,獲得以下結(jié)論。

補汽量的增加導致補汽結(jié)構(gòu)上下游透平級的氣動效率下降,特別是下游透平級的氣動效率下降值大于上游透平級。補汽導致下游透平級靜葉進口參數(shù)周向分布不均勻,造成下游透平級部分靜葉進口氣流角攻角變化較大,引起較大的葉柵攻角損失,由此降低了氣動效率。

補汽改變了機組通流部分補汽結(jié)構(gòu)下游靜葉汽封進口的氣動參數(shù)和預旋比。補汽量的增加導致下游靜葉汽封的進口預旋比增加,引起下游靜葉汽封的直接剛度系數(shù)增大。低頻區(qū)內(nèi),補汽量增加導致下游靜葉汽封的有效阻尼減小,補汽量變化所引起的下游靜葉汽封有效阻尼的差異隨著渦動頻率的增加而減小。

[1] QUINKERTZ R, THIEMANN T, GIERSE K. Validation of advanced steam turbine technology: a case study of an ultra super critical steam turbine power plant, GT2011-45816 [R]. New York, USA: ASME, 2011.

[2] 毛研偉, 竺曉程, 周代偉, 等. 超超臨界汽輪機抽汽系統(tǒng)數(shù)值研究 [J]. 汽輪機技術(shù), 2012, 54(2): 95-98. MAO Yanwei, ZHU Xiaocheng, ZHOU Daiwei, et al. Numerical investigation on extraction system in ultra-supercritical steam turbine [J]. Turbine Technology, 2012, 54(2): 95-98.

[3] 高俊鵬. 汽輪機回熱抽汽對下游流場影響的數(shù)值研究 [J]. 汽輪機技術(shù), 2010, 52(1): 49-50, 76. GAO Junpeng. Numerical investigation of the effect of regeneration extraction steam on downstream flow field [J]. Turbine Technology, 2010, 52(1): 49-50, 76.

[4] ROSIC B, MAZZONI C M, BIGNELL Z. Aerodynamic analysis of steam turbine feed-heating steam extraction [J]. ASME Journal of Engineering for Gas Turbines and Power, 2014, 136: 112602.

[5] 霍文浩, 晏鑫, 李軍, 等, 超超臨界汽輪機高壓缸旁路冷卻系統(tǒng)流動特性的研究 [J]. 西安交通大學學報, 2013, 47(3): 42-47. HUO Wenhao, YAN Xin, LI Jun, et al. Numerical investigations on flow characteristics of bypass cooling system in high pressure cylinder of ultra-supercritical steam turbines [J]. Journal of Xi’an Jiaotong University, 2013, 47(3): 42-47.

[6] 李志剛, 郞驥, 李軍, 等. 迷宮密封泄漏特性的試驗研究 [J]. 西安交通大學學報, 2011, 45(3): 48-52. LI Zhigang, LANG Ji, LI Jun, et al. Experimental investigations on the leakage flow characteristics of labyrinth seals [J]. Journal of Xi’an Jiaotong University, 2011, 45(3): 48-52.

[7] YAN Xin, LI Jun, FENG Zhenping. Investigations on the rotordynamic characteristics of a hole-pattern seal using transient CFD and periodic circular orbit model [J]. ASME Journal of Vibration and Acoustics, 2011, 133(4): 041007.

[8] LI Zhigang, LI Jun, YAN Xin. Multiple frequencies elliptical whirling orbit model and transient RANS solution approach to rotordynamic coefficients of annual gas seals prediction [J]. ASME Journal of Vibration and Acoustics, 2013, 135(3): 031005.

[9] CHILDS D W. Turbomachinery rotordynamics: phenomena, modeling and analysis [M]. New York, USA: Wiley, 1993: 292.

[10]VANCE J, ZIERER J, CONWAY E M. Effect of straight-through labyrinth seals on rotordynamics [C]∥Proceedings of the ASME Vibration and Noise Conference. New York, USA: ASME, 1993: 159-172.

(編輯 苗凌)

Numerical Investigations for Effect of Supplementary Steam on Aerodynamic Performance of Turbine Stage and Rotordynamics of Stator Seal

WANG Shizhu1, LI Zhigang1, LI Jun1, JI Dawei2, XIAO Gaohui2, GE Qing2

(1. Institute of Turbomachinery, Xi’an Jiaotong University, Xi’an 710049, China; 2. Shanghai Electric Power Generation Co., Ltd., Shanghai 200240, China)

To investigate the influence of supplementary steam on the aerodynamic performance of turbine through flow and rotordynamic characteristics of diaphragm seal, three-dimensional Reynolds-Averaged Navier-Stokes (RANS) andk-εturbulent model are used to analyze the aerodynamic performance of turbine stage without supplementary steam and with two kinds of supplementary steam. The inlet and outlet boundary conditions of the turbine stage stator seal at three kinds of supplementary steam operating conditions were specified according to flow field calculations. The unsteady fluid excitation rotordynamic characteristics of the stator seal were numerical calculated using multi-frequency elliptical whirling orbit model and dynamic grids. The numerical results show that the supplementary steam jets impacts on the hub regions of the downstream stator blade and results in the vortex flow. This flow behavior leads to the non-uniform inlet aerodynamic parameters at the entrance of the stator blade. The aerodynamic performance decreases with the increase of the mass flow rate of the supplementary steam. The supplementary steam jet changes the inlet preswirl and aerodynamic parameters of the downstream stator seal. The unsteady fluid excitation force and rotordynamic coefficients of the stator seal also change with the variation of the inlet flow parameters.

turbine stage; stator seal; aerodynamic performance; rotordynamics; numerical simulation

2014-09-22

王世柱(1990—),男,碩士生;李軍(通信作者),男,教授,博士生導師。

國家自然科學基金資助項目(51376144)。

時間:2015-02-27

10.7652/xjtuxb201505009

TK474.7

A

0253-987X(2015)05-0056-06

網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20150227.0845.001.html