平衡孔對軸流式透平氣動性能影響的 試驗和數(shù)值研究

鐘主海， 孫 奇

(東方汽輪機(jī)有限公司， 四川德陽 618000)

對于典型的沖動式透平級而言，平衡孔的存在使得葉輪兩側(cè)腔室相通，推力減小。后期在考慮平衡孔對效率的影響時，主要是從定性角度出發(fā)，認(rèn)為上游隔板漏氣量完全通過平衡孔流入下游腔室，不灌入主流，這可以降低動葉片根部的二次流損失。孫弼等[1]將邊界層計算引入到透平子午通流計算中，對某汽輪機(jī)高、低壓缸的計算結(jié)果表明，透平級靜、動葉根部軸向間隙內(nèi)的汽流經(jīng)葉輪上的平衡孔流入下一級，可以減少端部的二次流損失，改善汽輪機(jī)氣動性能。Moroz等[2]用數(shù)值方法詳細(xì)研究了帶平衡孔的葉輪腔室中的流動情況。計算結(jié)果表明，當(dāng)工質(zhì)密度相對較大時，腔室中出現(xiàn)了沿周向方向不穩(wěn)定的旋渦。Luo等[3]用數(shù)值方法詳細(xì)研究了帶平衡孔的葉輪腔室中的流動情況。計算結(jié)果表明，與無平衡孔相比，帶平衡孔方案效率降低1.81%。劉網(wǎng)扣等[4]對某沖動式透平級進(jìn)行了平衡孔面積與效率影響的研究，其研究模型為單級透平，結(jié)果表明：葉輪上開平衡孔可以明顯提高效率，與無平衡孔相比，隨著平衡孔面積的增大，效率逐漸提高，最佳平衡孔對應(yīng)效率與無平衡孔時效率相差約0.8%。Slma[5]用數(shù)值方法對單級透平模型進(jìn)行了平衡孔和徑向汽封齒對效率和軸向推力影響的研究，結(jié)果表明：徑向汽封齒對效率基本沒有影響，無平衡孔時效率比有平衡孔時低約0.3%。李慶等[6]對無平衡孔和有平衡孔的通流結(jié)構(gòu)進(jìn)行數(shù)值模擬，計算表明：有平衡孔的通流效率要低于無平衡孔的通流效率，在平衡孔數(shù)量固定的情況下，隨著平衡孔直徑的增大，效率逐漸降低，在平衡孔直徑固定的情況下，隨著平衡孔數(shù)量的增加，效率和軸向推力都呈現(xiàn)逐漸降低的趨勢。劉網(wǎng)扣[7]對沖動式機(jī)組級間漏氣對效率的影響進(jìn)行了試驗研究，結(jié)果表明：靜、動葉間存在少量流出主流的漏氣量時，透平能獲得較高的效率。肖曉玲等[8]對不同平衡孔孔徑進(jìn)行不同隔板汽封間隙的試驗研究，結(jié)果表明：隨著平衡孔孔徑及面積的增大，效率明顯提高，當(dāng)平衡孔面積增加到一定程度后，效率開始緩慢下降，即存在最佳平衡孔面積值，使得效率達(dá)到最佳。

從以上關(guān)于平衡孔對效率影響的研究情況來看，結(jié)論基本都是平衡孔可以緩解推力，但其存在相比無平衡孔時會影響效率。因此，在當(dāng)前純通流部分效率提升難度日益增加的情況下，妥善設(shè)計葉輪腔室可以使透平效率進(jìn)一步提升。那么，如何根據(jù)實(shí)際的透平級特性參數(shù)設(shè)計最佳的平衡孔面積，以達(dá)到最高效率，是當(dāng)前的主要研究課題之一。從前文所述的文獻(xiàn)情況來看，許多結(jié)論均是級間葉輪抽吸流量占主流比例需要達(dá)到1%～2%量級才能對應(yīng)最高效率，當(dāng)然這與上游隔板漏氣量占主流比例相關(guān)。一般對于典型的高壓缸而言，隔板漏氣量占主流比例在1%的量級，抽吸量超過漏氣量1倍或以上的設(shè)計無疑需要更加慎重，因為這部分氣流直接從主流出來，并未在動葉片中轉(zhuǎn)化做功，更重要的是，漏過平衡孔的氣流由于大大超過下游隔板的通流能力，有相當(dāng)部分氣流會在動葉片下游根部進(jìn)入主流，從而擾亂下游靜葉片的根部流動，降低下游級次的效率。

因此，要準(zhǔn)確設(shè)計葉輪腔室中的平衡孔以及評估平衡孔對葉輪腔室中流場的影響情況，需要定性研究葉輪腔室中的流動情況。目前大部分的研究都集中在數(shù)值模擬。從國內(nèi)外研究情況來看，我們發(fā)現(xiàn)葉輪腔室中的流動實(shí)際上是高度三維非定常黏性流動，因此，簡單的一維簡化理論很難直接得到通用且高效的工程設(shè)計方法[9]。本文的研究從工程快速準(zhǔn)確的設(shè)計原則出發(fā)，首先通過多級空氣透平試驗得到在試驗?zāi)Ｐ蜅l件下的級宏觀效率與平衡孔面積的變化規(guī)律，然后通過CFD計算模型來模擬試驗?zāi)Ｐ偷牧鲃?，試圖對比CFD宏觀效率與試驗效率變化規(guī)律，如果兩者在趨勢上吻合，則進(jìn)一步分析流場以得到效率變化的原因。

1 試驗條件與計算模型

1.1 試驗條件

本次試驗的目的主要是從試驗角度定量評估平衡孔對效率的影響大小，從工程應(yīng)用的角度出發(fā)，確定優(yōu)化平衡孔設(shè)計對效率提升的空間；另一方面則是驗證和校驗CFD方法，為后續(xù)CFD的研究提供基礎(chǔ)數(shù)據(jù)，為今后工程設(shè)計和優(yōu)化提供支撐。

試驗主要從宏觀性能出發(fā)，通過測量首級靜葉片前和末級動葉片后的氣流參數(shù)、流量以及輪周功率來得到多級透平的宏觀效率。本次試驗在東方汽輪機(jī)有限公司多級空氣透平實(shí)驗臺上進(jìn)行，多級空氣透平實(shí)驗臺如圖1所示。透平為3級沖動式機(jī)組模型，試驗實(shí)物模型如圖2所示。為了方便更換平衡孔，將平衡孔加工為內(nèi)螺紋，另外加工平衡螺塞，平衡孔尺寸分別為5 mm×16 mm(長度×直徑，下同)、5 mm×20 mm、5 mm×25 mm、5 mm×30 mm、5 mm×40 mm、5 mm×50 mm、5 mm×60 mm，每次更換試驗件只需將上半隔板拆開更換平衡螺塞，轉(zhuǎn)子、下半隔板等均不變，以保證其他結(jié)構(gòu)參數(shù)完全一致。

圖1 多級空氣透平實(shí)驗臺

圖2 三級透平實(shí)物模型

圖3為試驗?zāi)Ｐ褪疽鈭D。如圖3所示，在首級靜葉前有總溫和總壓探針，以及壁面靜壓測量孔；在末級動葉片下游的排氣導(dǎo)流板上，根部和頂部都具有壁面靜壓孔；在各葉片排之間也設(shè)有靜壓測量孔。同時，采用流量計和水力測功器分別測量試驗條件下的總流量和總輸出功率。多級空氣透平的試驗效率則可以由式(1)[10]計算得到：

(1)

式中：T0,in為入口總溫，K；p0,in為入口總壓，Pa；pex為出口靜壓，Pa；qm為質(zhì)量流量，kg/s；Tq為扭矩，N·m；ω為轉(zhuǎn)速,r/min；k為空氣等熵指數(shù)；cp為比定壓熱容，kJ/(kg·K)。

圖3 試驗?zāi)Ｐ褪疽鈭D

需要說明的是，出口靜壓由測量的根、頂部靜壓值平均之后得到。

表1給出了試驗?zāi)Ｐ拖鄳?yīng)的氣動和幾何數(shù)據(jù)。需要說明的是，方案中3級透平的根直徑、隔板汽封直徑和間隙等均按照相等設(shè)計。

表1 試驗主要?dú)鈩雍蛶缀螖?shù)據(jù)

1.2 計算模型

1.2.1 幾何模型

圖4為上述試驗?zāi)Ｐ湍臣塁FD計算模型的示意圖。其中，整個計算域包含靜止域(靜葉片流道域、靜葉片根部汽封域)和旋轉(zhuǎn)域(動葉片流道域、葉輪腔室域)。靜葉片流道域和動葉片流道域之間、靜葉片根部汽封域和葉輪腔室域之間都分別設(shè)置混合平面；動葉片流道域和葉輪腔室域則通過直接粘接起來。

圖4 單級CFD模型示意圖

1.2.2 網(wǎng)格劃分

圖5為本次試驗?zāi)Ｐ湍臣塁FD計算的網(wǎng)格示意圖。其中，單級靜、動葉片流道域的網(wǎng)格數(shù)約為30萬和27萬，葉輪腔室域網(wǎng)格數(shù)為28萬。第1層網(wǎng)格尺度的大小保證合理的壁面Y+值，本次模型Y+值控制在20以下，分布云圖如圖6所示。本次數(shù)值模擬僅考慮單級模型，其原因為：(1) 出于工程實(shí)際應(yīng)用的考慮，希望能夠快速評估；(2) 前文已經(jīng)提到，本次設(shè)計的3級透平方案，其根部反動度均按照相等的原則設(shè)計，且各級隔板汽封直徑、間隙等均相同，各級輪盤直徑、寬度和平衡孔大小、數(shù)目均相同，因此，從定性上講，各級輪盤腔室的流動應(yīng)該是基本相同的。

(a) 葉片域網(wǎng)格

圖6 Y+值云圖

1.2.3 控制方程求解及邊界條件

CFD計算采用全三維N-S方程和SST湍流模型，計算中對控制方程和邊界條件進(jìn)行有限元離散；動量、能量、湍動能以及湍流耗散率的離散格式均采用高分辨率的二階迎風(fēng)格式。計算工質(zhì)為空氣，入口給定均勻的總溫約80 ℃和總壓約145 kPa的氣流，進(jìn)氣方向為軸向；動葉片排氣口1給定平均靜壓約127 kPa；葉輪腔室出口2給定質(zhì)量流量，大小與上游隔板汽封漏氣量相等。

2 平衡孔漏氣對級性能的影響

2.1 平衡孔面積對效率影響的試驗結(jié)果

對于每個不同平衡孔直徑的方案，均測量了不同轉(zhuǎn)速下的性能參數(shù)，在變轉(zhuǎn)速的同時，透平級入口和出口的邊界參數(shù)基本保持不變，通過將轉(zhuǎn)速和級等熵焓降的等效速度無量綱化，可以得到效率與速比關(guān)系的曲線，本文中的速比均指按照圓周速度無量綱化后的速比。圖7為3種速比條件下，3級整體效率隨平衡孔直徑變化的情況，分別包含設(shè)計速比條件、高速比條件(為設(shè)計速比+0.1)和低速比條件(為設(shè)計速比-0.1)。從圖7可以看出，3種情況下，設(shè)計速比方案效率最高，速比變化量相同時，速比增大情況比速比減小情況的效率要高，這主要源于速比增大時級內(nèi)動葉片入口將出現(xiàn)負(fù)攻角，型面損失增大不顯著，且氣流折轉(zhuǎn)角減小，二次流損失變化不明顯。當(dāng)速比減小時，動葉片入口呈現(xiàn)正攻角，且氣流折轉(zhuǎn)角增大，二次流和型面損失均迅速增大。隨平衡孔直徑的增大，效率呈現(xiàn)先增大后減小的趨勢，即存在最佳平衡孔。當(dāng)前試驗?zāi)Ｐ蜅l件下，最佳平衡孔直徑為25～30 mm。由于葉輪腔室中的平衡孔流量與平衡孔總面積直接相關(guān)，而平衡孔總面積除了與直徑有關(guān)外，還與孔的個數(shù)有關(guān)，因此圖7并不能直觀表示效率受平衡孔面積的影響程度。圖8為設(shè)計速比條件下，效率隨平衡孔面積比的變化情況(面積比是指平衡孔總面積與上游隔板漏氣總面積的比值)。從圖8可以看出，本試驗方案的最佳平衡孔面積比為2～3，當(dāng)平衡孔面積比小于2時，效率隨平衡孔面積比的變化而變化得較快；當(dāng)平衡孔面積比大于3時，效率則隨平衡孔面積比的變化而變化得相對較慢。需要特別指出的是，與前文中的文獻(xiàn)[2]調(diào)研結(jié)果不同，本次試驗結(jié)果中，可以發(fā)現(xiàn)當(dāng)平衡孔面積比大于最佳值時，效率仍然會有顯著下降。

圖7 不同平衡孔直徑下效率的變化情況

圖8 效率隨平衡孔面積比的變化情況(設(shè)計速比)

2.2 隔板汽封間隙對效率影響的試驗結(jié)果

從定性認(rèn)識角度出發(fā)，平衡孔能夠改善效率的原因在于其能將上游隔板漏氣量部分或全部抽吸到葉輪下游腔室中，使得上游隔板泄漏氣流不會灌入級間的主流中，而對動葉片入口根部流場進(jìn)行干擾。當(dāng)平衡孔面積太大時，其將過多的主流氣流抽入到葉輪腔室中，一方面會減小動葉片的做功能力，另一方面，由于下游隔板汽封流動面積有限，故有部分通過上游平衡孔的氣流重新灌入到動葉片下游的主流中，對下游靜葉片根部入口流場進(jìn)行干擾。為此，補(bǔ)充試驗來驗證上述分析。補(bǔ)充試驗內(nèi)容是將各級隔板汽封間隙從設(shè)計值0.6 mm補(bǔ)充加工到1.0 mm，并測量平衡孔尺寸為5 mm×25 mm和5 mm×40 mm 2種方案下的效率。圖9為2種方案的對比情況。從圖9可以看出，平衡孔尺寸為5 mm×25 mm時，將隔板汽封間隙從0.6 mm增大到1.0 mm，效率降低約0.4%；平衡孔尺寸為5 mm×40 mm時，將隔板汽封間隙從0.6 mm增大到1.0 mm，效率僅僅下降0.15%左右，2種方案效率變化差異從側(cè)面證明了平衡孔面積需要匹配上游隔板漏氣面積設(shè)計，而兩者之間的具體面積比則與隔板汽封類型、齒數(shù)和驅(qū)動平衡孔流動的葉輪兩側(cè)壓差，或動葉片的根部反動度等參數(shù)直接相關(guān)。

(a) 平衡孔直徑為25 mm

2.3 平衡孔面積對效率影響的數(shù)值計算與試驗結(jié)果分析

前文已經(jīng)提到，平衡孔面積對效率的影響比平衡孔直徑對效率的影響的描述更為合理，而在隔板漏氣面積一定時，平衡孔面積比與效率關(guān)聯(lián)則更顯得有物理意義。圖10為設(shè)計速比條件下，CFD數(shù)值計算結(jié)果和試驗結(jié)果的對比情況。從圖10可以看出，CFD數(shù)值計算結(jié)果與試驗結(jié)果基本上有相同的最佳平衡孔面積比(對應(yīng)最高效率)，為2～3；CFD數(shù)值計算結(jié)果和試驗曲線基本上具有類似的效率-面積比變化趨勢，即最低效率點(diǎn)均在無平衡孔時，如果平衡孔的面積繼續(xù)增大，則效率還會有一定程度的降低。CFD數(shù)值計算結(jié)果和試驗結(jié)果同樣存在一定的差異，主要體現(xiàn)在兩方面：(1) 在最高效率點(diǎn)后，試驗曲線的效率變化程度大于CFD數(shù)值計算結(jié)果，尤其是在平衡孔面積比為7之后的區(qū)域；(2) CFD數(shù)值計算結(jié)果的最高效率與最低效率差異為0.6%左右，而試驗結(jié)果差異為0.8%左右，CFD數(shù)值計算結(jié)果差異小于試驗結(jié)果。上述2方面差異是由于本文中的單級模型簡化，但計算結(jié)果仍然存在諸多不確定性，如定常算法、湍流模型的影響等。如前所述，本文的目的不在于定量地數(shù)值研究流場和機(jī)理，而是評估簡化的CFD方法是否能夠評估出合理的平衡孔面積設(shè)計范圍。通過與試驗結(jié)果的對比，可以發(fā)現(xiàn)本節(jié)中簡化的CFD模型能夠評估出相對準(zhǔn)確的平衡孔面積-效率變化趨勢，因此，可以通過CFD計算結(jié)果進(jìn)一步判斷效率變化的原因。

圖10 平衡孔面積比對效率影響的CFD數(shù)值計算結(jié)果與試驗結(jié)果的對比

圖11則是不同平衡孔直徑下動葉片上下游軸向平面處的熵分布云圖(熵變化范圍0.96～1 kJ/(kg·K))，圖中的上下游位置與圖4中的平面1和平面2一一對應(yīng)。從圖11可以進(jìn)一步看出：當(dāng)無平衡孔或平衡孔直徑較小時，在動葉片上游平面根部處出現(xiàn)了高熵增流體，這股流體經(jīng)過動葉片根部的二次流進(jìn)一步發(fā)展后，使得動葉片下游出現(xiàn)了更大的高熵增區(qū)域；當(dāng)平衡孔直徑大于20 mm，動葉片上游根部區(qū)域的熵增分布云圖基本相同，且范圍較小，這主要是由于上游靜葉片端部二次流發(fā)展。因此，僅從動葉片的損失來看，進(jìn)一步增大平衡孔尺寸將不能再降低動葉片損失。但過多的抽吸氣流會經(jīng)過下游葉輪腔室重新灌入主流，影響下游靜葉流動，因此效率會降低。圖12給出了平衡孔直徑為7 mm、25 mm、50 mm時，葉輪腔室內(nèi)的流線和壓力分布圖。由圖12可以看出：隨著平衡孔面積的增大，平衡孔兩側(cè)壓差逐漸減小，這是導(dǎo)致平衡孔流量增大程度減緩的主要因素；平衡孔出口處的流量和流速都較大時，氣流會在出口處沖擊到壁面上，形成旋渦，具有一定的阻塞作用；平衡孔面積變化時，接近主流區(qū)域的葉輪腔室壓力變化很小，因此平衡孔面積對級根部反動度影響較小。

(a) 平衡孔直徑7 mm

3 結(jié) 論

(1) 平衡孔的設(shè)計可以有效提高效率，相比無平衡孔時提高的量級可以達(dá)到0.8%左右。對于當(dāng)前的試驗?zāi)Ｐ投?，最佳平衡孔面積比在2～3。

(2) CFD數(shù)值計算和試驗曲線基本上具有類似的效率-面積比變化趨勢，即最低效率點(diǎn)都在無平衡孔時，當(dāng)平衡孔的面積比大于最佳值時，繼續(xù)增大平衡孔面積，則效率還會存在一定程度的降低。

(3) CFD數(shù)值計算的流場分析證明了葉輪微抽吸時，可以改善動葉片上游端部邊界層狀態(tài)，從而降低動葉片端部二次流損失，提高效率；但進(jìn)一步提高葉輪抽吸流量對改善動葉片端部二次流損失并無作用，相反會導(dǎo)致下游葉輪腔室的氣流重新灌入主流，影響下游靜葉片的端部二次流，從而降低效率。