小容積流量下汽輪機(jī)調(diào)節(jié)級動葉三維流場與級后工質(zhì)溫度的研究*

李 勇 王涵北

(東北電力大學(xué)能源與動力工程學(xué)院)

小容積流量下汽輪機(jī)調(diào)節(jié)級動葉三維流場與級后工質(zhì)溫度的研究*

李 勇**王涵北

(東北電力大學(xué)能源與動力工程學(xué)院)

以某300MW汽輪機(jī)調(diào)節(jié)級為研究對象，介紹了比例系數(shù)法和數(shù)值模擬法的基本原理。對不同相對容積流量比和不同工況下的調(diào)節(jié)級動葉三維流場進(jìn)行了分析。將數(shù)值模擬方法得到的調(diào)節(jié)級級后工質(zhì)溫度與比例系數(shù)法、速度三角形法的計(jì)算結(jié)果進(jìn)行對比，得到不同計(jì)算方法間的相對誤差，并對不同相對容積流量下級后工質(zhì)溫度的計(jì)算方法進(jìn)行了探討。

汽輪機(jī) 動葉三維流場 級后工質(zhì)溫度 相對容積流量比

隨著風(fēng)能、太陽能等清潔能源在電力能源系統(tǒng)中所占比例的不斷增加，要求參與調(diào)峰的大容量火電機(jī)組必須通過降低負(fù)荷，甚至是啟停調(diào)峰的方式來適應(yīng)新能源發(fā)電的介入。因此，火電機(jī)組的快速啟停工作受到國內(nèi)外學(xué)者的廣泛關(guān)注。在確定火電機(jī)組啟停速度和壽命損耗中，汽輪機(jī)調(diào)節(jié)級后的溫度變化率對優(yōu)化火電機(jī)組啟停具有重要影響。為確定主蒸汽流量較小時的級后工質(zhì)參數(shù)，文獻(xiàn)[1～4]提出當(dāng)汽輪機(jī)負(fù)荷低于80%額定負(fù)荷時，采用比例系數(shù)法進(jìn)行參數(shù)計(jì)算，然而這種估算方法并未考慮汽輪機(jī)沖轉(zhuǎn)流量變化帶來的影響，必然引起調(diào)節(jié)級級后溫度變化率產(chǎn)生較大誤差。

隨著CFD技術(shù)的興起與日漸成熟，數(shù)值模擬方法在汽輪機(jī)調(diào)節(jié)級和小容積流量研究中獲得了廣泛應(yīng)用。文獻(xiàn)[5～7]通過模擬調(diào)節(jié)級三維流場，分析了工質(zhì)經(jīng)過調(diào)節(jié)級堵塞區(qū)域時通道中靜壓的變化情況。謝金偉等通過建立完整的調(diào)節(jié)閥至調(diào)節(jié)級段模型，對不同工況下該段的三維流動與壓力損失特性進(jìn)行了研究[8]。文獻(xiàn)[9～12]對汽輪機(jī)調(diào)節(jié)級氣動性能和氣流激振突變方面進(jìn)行了研究，提出了優(yōu)化調(diào)節(jié)級結(jié)構(gòu)的可行方案。文獻(xiàn)[13,14]通過對低負(fù)荷運(yùn)行的汽輪機(jī)末幾級進(jìn)行數(shù)值模擬，得到當(dāng)容積流量小到一定程度時末級動葉會做負(fù)功的結(jié)論。張仲彬等對多級流動進(jìn)行了數(shù)值分析，結(jié)果表明，小容積流量工況時上級流場的分離流動會導(dǎo)致下級流場內(nèi)分離區(qū)域增多[15]。

目前國內(nèi)外學(xué)者普遍關(guān)注的是小容積流量工況下汽輪機(jī)的末級流動狀態(tài)，而調(diào)節(jié)級小容積流量方面的研究卻鮮有報道。在此，筆者對不同相對容積流量比和不同工況下的調(diào)節(jié)級動葉三維流場進(jìn)行分析，并將采用數(shù)值模擬方法得到的調(diào)節(jié)級級后工質(zhì)溫度與比例系數(shù)法、速度三角形法計(jì)算結(jié)果進(jìn)行對比，得到不同計(jì)算方法間的相對誤差，以指導(dǎo)不同容積流量下級后工質(zhì)溫度的計(jì)算。

1 比例系數(shù)法

比例系數(shù)法中首先根據(jù)熱力特性說明書查出汽輪機(jī)額定負(fù)荷下的主蒸汽溫度、壓力，再進(jìn)行計(jì)算：

(1)

(2)

式中p0、p01——額定工況和變工況下的主蒸汽壓力，MPa；

pi、pi1——額定工況和變工況下某一級特征面的蒸汽壓力，MPa；

t0、t01——額定工況和變工況下的主蒸汽溫度，℃；

ti、ti1——額定工況和變工況下某一級特征面的蒸汽溫度，℃。

由式(1)、(2)可知，比例系數(shù)法存在兩方面問題。首先，對于汽輪機(jī)沖轉(zhuǎn)階段而言，主蒸汽參數(shù)基本不變，因此該方法計(jì)算出的調(diào)節(jié)級級后工質(zhì)參數(shù)也不隨時間改變。而事實(shí)上，汽輪機(jī)沖轉(zhuǎn)時需要調(diào)節(jié)閥全開來保證汽缸和轉(zhuǎn)子受熱均勻，因此多數(shù)電廠采用高壓主汽閥控制汽輪機(jī)轉(zhuǎn)速以便進(jìn)行低速檢查和中速暖機(jī)，這就使主汽閥開度在沖轉(zhuǎn)階段隨時間變化，進(jìn)而引起調(diào)節(jié)級入口、出口工質(zhì)參數(shù)隨時間改變。其次，式(1)、(2)中沒有體現(xiàn)沖轉(zhuǎn)流量改變對調(diào)節(jié)級級后工質(zhì)參數(shù)的影響。

2 物理模型和數(shù)值方法

2.1 物理模型

筆者對某300MW汽輪機(jī)調(diào)節(jié)級葉片及其流域進(jìn)行建模，葉片幾何尺寸見表1。

表1 某300MW汽輪機(jī)調(diào)節(jié)級葉片幾何尺寸

采用Gambit軟件建模時，由表1可知，調(diào)節(jié)級動、靜葉片數(shù)不同，同一工況下單位靜葉流道與單位動葉流道的質(zhì)量流量不同，因此不能采用單葉片模型進(jìn)行研究，而建立整周模型所需網(wǎng)格數(shù)量大，容易導(dǎo)致計(jì)算機(jī)內(nèi)存超限并延長計(jì)算時間，經(jīng)觀察發(fā)現(xiàn)，5個靜葉與4個動葉構(gòu)成的流域具有相同質(zhì)量流量，遂建立如圖1所示的物理模型。筆者選用的調(diào)節(jié)級模型部分進(jìn)汽度高達(dá)0.96，建模時可忽略其影響。

在調(diào)節(jié)閥全開的情況下，過熱蒸汽通過汽輪機(jī)蒸汽室進(jìn)入調(diào)節(jié)級靜葉，經(jīng)噴嘴加速、動葉做功后流入調(diào)節(jié)級汽室進(jìn)行充分混合。為保證建立的

圖1 某300MW汽輪機(jī)調(diào)節(jié)級物理模型

模型與實(shí)際情況相似，同時考慮到網(wǎng)格劃分的需要，設(shè)置動葉入口截面距動葉前緣點(diǎn)0.2mm，動葉出口截面距動葉尾緣點(diǎn)0.2mm，并在動葉后建立軸向距離等于真實(shí)長度的調(diào)節(jié)級汽室模型，在保證動葉出口流體充分發(fā)展的同時使流動更加真實(shí)。汽室部分在動葉生成流域的基礎(chǔ)上切割而成，不但提高了動葉與汽室連接處的順滑性，也為劃分高質(zhì)量網(wǎng)格提供了保證。葉柵通道網(wǎng)格(結(jié)構(gòu)化網(wǎng)格)劃分情況如圖2所示。

圖2 葉柵通道網(wǎng)格劃分情況

2.2 網(wǎng)格無關(guān)性討論

為兼顧計(jì)算速度和準(zhǔn)確度，筆者對模型進(jìn)行了網(wǎng)格無關(guān)性討論。在保證計(jì)算方法、邊界條件不變的前提下，分別對網(wǎng)格數(shù)為184萬、292萬、346萬、422萬的模型進(jìn)行數(shù)值模擬。結(jié)果表明，當(dāng)網(wǎng)格數(shù)大于346萬后，隨著網(wǎng)格數(shù)的增加模擬結(jié)果幾乎不變，故綜合考慮后取346萬網(wǎng)格的方案。

2.3 數(shù)值方法

采用Fluent軟件進(jìn)行計(jì)算，選擇適用于參數(shù)研究的標(biāo)準(zhǔn)k-ε湍流模型，應(yīng)用Simplec算法和二階迎風(fēng)離散格式，采用旋轉(zhuǎn)周期性邊界條件，動葉流域及其固體壁面均繞x軸以給定轉(zhuǎn)速旋轉(zhuǎn)，靜葉、汽室流域及其固體壁面保持靜止。靜葉、動葉、調(diào)節(jié)級汽室三者交界面采用滑移網(wǎng)格技術(shù)進(jìn)行數(shù)據(jù)交換。入口邊界條件為壓力入口，出口邊界條件為質(zhì)量流量出口(mass-flow inlet中的outward normals)，過熱蒸汽采用相較于理想氣體更接近真實(shí)水蒸氣性質(zhì)的Redlich-Kwong氣體模型，能量方程收斂精度為10-6，速度方程收斂精度為10-3。

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

3.1 動葉三維流場的分析

由文獻(xiàn)[15]可知，當(dāng)本級靜葉入口工質(zhì)均勻時，容積流量減小很難引起靜葉流域發(fā)生明顯的流動分離，故主要分析動葉三維流場的情況。為了排除汽輪機(jī)轉(zhuǎn)速對流場的影響，將動葉及其流域轉(zhuǎn)速均設(shè)為3kr/min。相對容積流量比k為：

k=G1v1/(Gv)

(3)

式中G、G1——閥門全開和變工況下調(diào)節(jié)級動葉出口蒸汽質(zhì)量流量，kg/s；

v、v1——閥門全開和變工況下調(diào)節(jié)級動葉出口蒸汽比容，m3/kg。

以閥門全開工況為基準(zhǔn)，圖3為不同k下同一動葉流道50%葉高截面處的流線圖，圖4為k為0.290時同一工況下不同葉高截面處的流線圖。

圖3 不同k下同一動葉流道50%葉高截面處的流線圖

圖4 同一工況下不同葉高截面處的流線圖

由圖3可知，隨著k值的減小，調(diào)節(jié)級動葉壓力面前緣點(diǎn)附近會出現(xiàn)流動分離，分離區(qū)與主流之間有一條明顯的分離線。出現(xiàn)流動分離的原因是容積流量較小時調(diào)節(jié)級靜葉出口汽流絕對速度較小，由于動葉高速旋轉(zhuǎn)，導(dǎo)致汽流以較大負(fù)攻角進(jìn)入動葉流域并被吸力面分為兩支，其中一支在旋轉(zhuǎn)葉片的作用下沿本流道繼續(xù)向下游運(yùn)動，另一分支流入相鄰流道并在動葉壓力面附近產(chǎn)生漩渦。當(dāng)k>0.026時，壓力面漩渦軸向影響范圍隨k值的減小而增大，這是因?yàn)榇藭r負(fù)攻角較大，而且容積流量相對較高，工質(zhì)經(jīng)堵塞段加速后具有較大的軸向分速度，使再附點(diǎn)逐漸遠(yuǎn)離分離點(diǎn)；當(dāng)k<0.026時，動葉出口吸力面附近出現(xiàn)了流動分離且分離區(qū)隨k值的減小而增大，此分離區(qū)的出現(xiàn)使動葉內(nèi)工質(zhì)橫向流動明顯加劇，加之此時容積流量過小，汽流軸向分速度較低，使壓力面回流區(qū)軸向影響范圍稍有減小。

由圖4可知，同一工況下不同葉高截面處的流動狀態(tài)存在差異，壓力面流動分離現(xiàn)象首先在動葉根部前緣點(diǎn)附近產(chǎn)生，且當(dāng)k=0.290時回流區(qū)范圍隨葉高的增加而減小直至消失。為了進(jìn)一步說明壓力面發(fā)生流動分離時的動葉三維流場情況，引入圖5所示的不同k值時動葉壓力面和下端壁極限流線圖。

a． k=0.290

b． k=0.330

圖5a中k=0.290時的工況與文獻(xiàn)[16]類似，此時在動葉壓力面與下端壁附近均發(fā)生了流動分離。k=0.290時動葉壓力面前緣點(diǎn)附近沿葉高方向形成一條分離線，靠近分離線的下游流體在40%左右軸向弦長內(nèi)形成了大面積回流區(qū)且漩渦占據(jù)了整個葉高，導(dǎo)致近壁區(qū)域出現(xiàn)強(qiáng)烈的回流和徑向流動。40%軸向弦長及其后區(qū)域徑向流動逐步減弱，但整體流線呈斜下方趨勢。此外，k=0.290時下端壁附近也出現(xiàn)了明顯的流動分離，分離區(qū)內(nèi)反向流體沿斜上方運(yùn)動并向分離線靠攏。對比發(fā)現(xiàn)，k=0.330時壓力面與下端壁附近均未出現(xiàn)回流，流動效果較好。綜上所述，k=0.290為調(diào)節(jié)級動葉流域出現(xiàn)流動分離現(xiàn)象的起始點(diǎn)。

3.2 級后工質(zhì)溫度的分析

速度三角形法和比例系數(shù)法都是計(jì)算級后工

質(zhì)參數(shù)的傳統(tǒng)方法，由于筆者主要介紹數(shù)值模擬方法，因此只將新方法與傳統(tǒng)方法所得結(jié)果進(jìn)行比較。假設(shè)此時主汽閥和調(diào)節(jié)閥全開，則3種方法計(jì)算結(jié)果之間的相對誤差如圖6所示。

圖6 不同級后工質(zhì)溫度計(jì)算方法間

由圖6可知，在k>0.026范圍內(nèi)數(shù)值模擬法與速度三角形法計(jì)算結(jié)果的相對誤差小于2%，且基本保持不變，這是因?yàn)檎{(diào)節(jié)級徑高比較大，流場沿半徑方向雖有變化但變化程度較小，平均直徑處汽流參數(shù)仍適用于變工況的計(jì)算，雖然k≤0.290時動葉壓力面附近已出現(xiàn)流動分離，但此時仍可用速度三角形法估算級后工質(zhì)溫度；當(dāng)k≤0.026時，由于動葉出口吸力面附近出現(xiàn)了流動分離，使平均直徑處汽流出汽角相比于設(shè)計(jì)工況發(fā)生了較大變化，導(dǎo)致兩種方法計(jì)算結(jié)果的相對誤差呈大幅上升趨勢，此時速度三角形法得到的級后工質(zhì)溫度存在較大偏差，應(yīng)采用數(shù)值模擬方法進(jìn)行計(jì)算。比例系數(shù)法計(jì)算結(jié)果與數(shù)值模擬結(jié)果的相對誤差在k=0.330時已大于2%，且隨k值的減小而增大，在k<0.110范圍內(nèi)增幅更為明顯，可見比例系數(shù)法所得結(jié)果在容積流量較小時并不理想。綜上所述，為得到準(zhǔn)確的調(diào)節(jié)級級后工質(zhì)溫度，當(dāng)k>0.026時可采用速度三角形法進(jìn)行計(jì)算，當(dāng)k≤0.026時應(yīng)采用數(shù)值模擬方法。

4 結(jié)論

4.1對動葉三維流場的分析表明，當(dāng)k=0.290時，由于入口汽流負(fù)攻角較大，下端壁和壓力面附近首次出現(xiàn)流動分離。當(dāng)0.026

4.2對數(shù)值模擬結(jié)果與速度三角形法計(jì)算結(jié)果的分析表明，當(dāng)k>0.026時，兩者的相對誤差小于2%，雖然在0.026

4.3對數(shù)值模擬結(jié)果與比例系數(shù)法計(jì)算結(jié)果的分析表明，當(dāng)k≤0.330時，兩者相對誤差大于2%，可見容積流量較小時，采用比例系數(shù)法所得結(jié)果存在較大誤差，相比之下，數(shù)值模擬結(jié)果準(zhǔn)確度更高。

[1] 陳堅(jiān)紅,阮偉,任浩仁,等.冷態(tài)啟動溫升率對汽輪機(jī)轉(zhuǎn)子壽命影響的研究[J].動力工程,1997,17(6):11～15.

[2] 袁鵬飛,盛德仁,陳堅(jiān)紅,等.汽輪機(jī)轉(zhuǎn)子熱疲勞壽命損耗監(jiān)測面的有限元分析[J].電站系統(tǒng)工程,2004,20(2):43～44.

[3] 付忠廣,杜志鋒,王夢嬌.660MW超超臨界汽輪機(jī)轉(zhuǎn)子熱應(yīng)力的研究[J].現(xiàn)代電力,2012,29(1):60～64.

[4] 袁利軍.300MW機(jī)組中壓缸啟動熱應(yīng)力計(jì)算與啟動優(yōu)化[D].北京:華北電力大學(xué),2005.

[5] 訾宏達(dá),方韌,徐星仲,等.汽輪機(jī)部分進(jìn)汽調(diào)節(jié)級三維流場的非定常數(shù)值模擬[J].動力工程,2007,27(5):721～724.

[6]Lampart P,Szymaniak M,Rzadkowski R.Unsteady Load of Partial Admission Control Stage Rotor of a Lager Power Steam Turbine[C].ASME Turbo Expo 2004:Power for Land,Sea and Air.Glasgow:American Society of Mechanical Engineers, 2004:237～246.

[7] Hills N J.Whole Turbine CFD Modelling[C].ASME Turbo Expo 2007:Power for Land,Sea and Air.Montreal:American Society of Mechanical Engineers,2007:817～824.

[8] 謝金偉,王新軍,廖高良,等.汽輪機(jī)調(diào)節(jié)閥-調(diào)節(jié)級段三維流動與壓力損失數(shù)值研究[J].中國電機(jī)工程學(xué)報,2014,34(20): 3376～3383.

[9] Hushmandi N B,Hu J,Fridh J,et al.Numerical Study of Unsteady Flow Phenomena in a Partial Admission Axial Steam Yurbine[C].ASME Turbo Expo 2008:Power for Land,Sea and Air.Berlin:American Society of Mechanical Engineers,2008: 713～722.

[10] He L.Computation of Unsteady Flow through Steam Turbine Blade Row at Partial Adimission[J].Proceeding of the Institution of Mechanical Engineers Part A:Journal of Power and Energy,1997,211(3):197～205.

[11] 陽虹,楊建道,李軍,等.汽輪機(jī)調(diào)節(jié)級氣動性能分析和結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].熱力透平,2011,40(1):5～10.

[12] 李建蘭,郭緒宏,余辰,等.汽輪機(jī)調(diào)節(jié)級汽流激振突變模型研究[J].中國電機(jī)工程學(xué)報,2013,33(11):39～46.

[13] 楊建道.空冷50MW機(jī)組低壓缸變工況三維流場數(shù)值分析[J].熱力透平,2007,36(1):53～56.

[14] Senoo S,Segawa K,Hamatake H,et al.Computation for Unsteady Compressible Flows in a Multi-Stage Steam Turbine with Steam Properties at Low Load Operations[J].Journal of Engineering for Gas Turbines and Power,2011,133(10): 711～722.

[15] 張仲彬,田艷靜,曹麗華,等.小容積流量工況下汽輪機(jī)級內(nèi)流場的數(shù)值分析[J].化工機(jī)械,2013,40(1):94～97.

[16] 朱光宇,俞茂錚,戴麗萍,等.汽輪機(jī)調(diào)節(jié)級動葉柵大負(fù)沖角工況下的三維分離流動研究[J].動力工程學(xué)報,2004,24(2): 189～194.

Studyof3DFlowFieldsinGoverningStageandSteamTemperatureafterStageofSteamTurbineatSmallVolumetricFlow

LI Yong, WANG Han-bei

(SchoolofEnergyandPowerEngineering,NortheastDianliUniversity,Jilin132012,China)

Taking the governing stage of a 300MW steam turbine as the object of study, the basic principle of proportional coefficient method and numerical simulation method was described; and the 3D flow fields in moving blades at governing stage with different volume flow ratios and at various working conditions were analyzed. Through comparing the steam temperature after governing stage and from the numerical simulation method with the results from the proportional coefficient method and velocity triangle method, the relative error between various calculation methods can be reached. These different calculation methods for steam temperature at different volumetric flows were discussed.

steam turbine, 3D flow field in moving blade, steam temperature after stage, relative volume flow ratio

*國家自然科學(xué)基金項(xiàng)目(51376041)。

**李 勇，男，1964年1月生，教授。吉林省吉林市，132012。

TQ053

A

0254-6094(2016)04-0489-06

2015-09-21)