黎 軍，薛偉鵬

(中國航發(fā)四川燃?xì)鉁u輪研究院，成都610500)

1 引言

由于葉片排間的相對運(yùn)動，轉(zhuǎn)靜葉片排間需有一定的軸向間隙避免碰磨。同時，為防止流道中高溫燃?xì)馔ㄟ^葉片排間的軸向間隙進(jìn)入盤腔內(nèi)部，需在渦輪盤腔內(nèi)引入一股封嚴(yán)冷氣，并通過葉片排間間隙進(jìn)入渦輪級內(nèi)。該股封嚴(yán)冷氣與通道中的主流流動和二次流動發(fā)生摻混，導(dǎo)致渦輪性能惡化。因此封嚴(yán)冷氣對渦輪性能的影響引起國內(nèi)外學(xué)者的廣泛關(guān)注。Dénos等[1]對一單級軸流渦輪在不同封嚴(yán)冷氣流量下進(jìn)行的試驗(yàn)研究發(fā)現(xiàn)，1.49%的封嚴(yán)冷氣導(dǎo)致導(dǎo)葉出口7%的靜壓升，相應(yīng)的導(dǎo)葉出口馬赫數(shù)降低，從而使渦輪性能明顯改變，渦輪功率降低2.8個百分點(diǎn)。張晶輝等[2]的非定常數(shù)值研究表明，封嚴(yán)冷氣與上游導(dǎo)葉尾跡的相互作用引起轉(zhuǎn)子通道內(nèi)熵增，當(dāng)封嚴(yán)冷氣量占主流流量的1.37%時，渦輪效率降低2.1%。由此可見，對于高性能渦輪，該股封嚴(yán)冷氣的影響不容忽視。

目前封嚴(yán)冷氣對渦輪主流影響的數(shù)值研究，主要有設(shè)定源項(xiàng)邊界模擬封嚴(yán)流動和建立封嚴(yán)模型并與主流進(jìn)行兩股流動耦合分析兩種方法。Turner等[3]采用源項(xiàng)方法通過給定軸向和徑向動量模擬低壓導(dǎo)向器緣板的封嚴(yán)泄漏流動，封嚴(yán)泄漏流的流量、溫度和旋向通過一個單獨(dú)的迷宮封嚴(yán)分析程序計算獲得，結(jié)果表明封嚴(yán)泄漏對端壁區(qū)的溫度剖面有很大影響。Mirzamoghadam等[4-5]對級間交界面分別位于封嚴(yán)腔上游和下游的定常與非定常計算結(jié)果進(jìn)行了對比，認(rèn)為導(dǎo)葉尾跡對封嚴(yán)流動的影響比動葉前緣勢場的影響更明顯。相比交界面位于封嚴(yán)腔上游，采用級間交界面位于封嚴(yán)腔下游的定常計算，能夠與非定常計算吻合得更好，并在文中提出了最小封嚴(yán)流量模型。Chilla等[6]對帶真實(shí)級間封嚴(yán)結(jié)構(gòu)的渦輪級進(jìn)行了定常和非定常數(shù)值模擬，認(rèn)為封嚴(yán)冷氣和主流之間的相互作用受到封嚴(yán)冷氣與主流之間速度虧損的嚴(yán)重影響，同時也受到由動葉產(chǎn)生的非均勻周向壓力場影響。周楊等[7]采用簡化結(jié)構(gòu)的封嚴(yán)模型與主流耦合分析，就不同封嚴(yán)流量對高壓渦輪二次流動進(jìn)行了研究。胡松林等[8]研究了封嚴(yán)冷氣噴射角度對主流的影響。張偉昊[9]、孔祥林[10]等對封嚴(yán)結(jié)構(gòu)的結(jié)構(gòu)形式和結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化，了解了其對主流的影響。

從公開文獻(xiàn)看，大部分采用封嚴(yán)流體域與主流域進(jìn)行流動耦合分析的方法來開展研究工作，以期獲得更豐富的封嚴(yán)冷氣與主流相互作用的流動現(xiàn)象，了解其相互影響機(jī)理。但鑒于源項(xiàng)方法的便利性，在需要開展大量迭代計算的工程應(yīng)用領(lǐng)域，目前源項(xiàng)方法還是得到了廣泛采用。然而由于源項(xiàng)方法給定的封嚴(yán)冷氣進(jìn)口邊界對渦輪性能評估結(jié)果有很大的影響，如何通過源項(xiàng)方法準(zhǔn)確模擬封嚴(yán)流動的研究還開展得較少，為此本文重點(diǎn)對比不同源項(xiàng)進(jìn)口邊界(氣流角和所在計算域)對渦輪性能的影響。在此基礎(chǔ)上，采用封嚴(yán)冷氣流與主流耦合分析方法，獲取封嚴(yán)冷氣在真實(shí)盤腔結(jié)構(gòu)環(huán)境下的流動特點(diǎn)，據(jù)此對源項(xiàng)邊界進(jìn)行修正，以提高工程方法的計算精度，并對封嚴(yán)冷氣對渦輪性能的影響機(jī)理進(jìn)行深入分析。

2 源項(xiàng)模擬結(jié)果對比

計算模型選用某高負(fù)荷跨聲渦輪的葉型和流道，計算網(wǎng)格采用TurboGrid 13.0劃分，網(wǎng)格節(jié)點(diǎn)總數(shù)為95萬。計算了渦輪級性能評估中常用的四種不同封嚴(yán)冷氣邊界條件和一組無封嚴(yán)冷氣的結(jié)果。封嚴(yán)冷氣源項(xiàng)進(jìn)口分別設(shè)定在靜子計算域和轉(zhuǎn)子計算域上，徑向氣流角給定為11.3°，周向氣流角分別給定為軸向和與當(dāng)?shù)亓飨蛳喈?dāng)?shù)那邢驓饬鹘?，進(jìn)口流量給定為渦輪進(jìn)口流量的0.6%。詳細(xì)邊界信息和計算結(jié)果見表1。

表1 不同封嚴(yán)冷氣進(jìn)口邊界的計算結(jié)果Table 1 The results of different boundary conditions of the rim seal flow

對比表1結(jié)果可見，封嚴(yán)冷氣在靜子計算域上的計算效率低于轉(zhuǎn)子計算域。這一方面是由于靜子計算域上引入封嚴(yán)冷氣時，封嚴(yán)冷氣與主流之間的速度差異導(dǎo)致在級間交界面上的周向平均引起摻混損失增加，從而效率降低；另一方面，在轉(zhuǎn)子計算域上給定封嚴(yán)冷氣的源項(xiàng)進(jìn)口時，實(shí)際上已經(jīng)給封嚴(yán)冷氣加載了一個與轉(zhuǎn)速相關(guān)的切向動量，而該動量在進(jìn)入葉片通道后參與了做功，因此相比封嚴(yán)冷氣在靜子計算域的情況功率略有增加，從而效率有所提高。對比封嚴(yán)冷氣軸向和切向進(jìn)入的情況可見，當(dāng)封嚴(yán)冷氣邊界給定為具有切向速度分量時，計算效率相比給定軸向速度邊界的算例略有增加。這是因?yàn)樵谙嗤鈬?yán)冷氣縫面積下，給定切向進(jìn)氣和流量，封嚴(yán)冷氣可獲得較高的流速，與主流流速差較小。同時，其與主流的流動方向的差異減小，也將減少該封嚴(yán)冷氣流進(jìn)入主流道后與主流的摻混損失，并更好地參與做功。綜上所述，不同的級間封嚴(yán)冷氣邊界設(shè)置，將對渦輪級的計算效率和反力度等性能參數(shù)產(chǎn)生影響，不同算例之間計算效率的最大差異為0.6個百分點(diǎn)。

圖1給出了不同級間封嚴(yán)冷氣位置和方向邊界下的封嚴(yán)冷氣流線對比，從圖中可看出不同邊界設(shè)置對動葉根部二次流動的影響。封嚴(yán)冷氣邊界加載在靜子計算域中時，由于經(jīng)過交界面的周向平均，封嚴(yán)冷氣對導(dǎo)葉出口速度矢量產(chǎn)生的不均勻影響將被抹平，因此傳遞到動葉進(jìn)口的不均勻信息減弱，使得動葉端壁附面層中流體受到封嚴(yán)冷氣擾動較小，動葉進(jìn)口附面層形成的二次流發(fā)展較緩慢，由此將減小動葉根部通道渦的影響范圍。封嚴(yán)冷氣邊界加載在轉(zhuǎn)子計算域中時，封嚴(yán)冷氣流更容易受到通道二次流動的影響，更早地脫離動葉緣板并卷入通道渦，影響更大范圍的葉高。當(dāng)封嚴(yán)冷氣邊界給定切向分量時，封嚴(yán)冷氣流與主流速度矢量差異減小，更不容易被動葉根部二次流動所卷吸，產(chǎn)生的通道渦及其影響范圍也相應(yīng)減小。

通過以上對比可看出，采用不同的邊界設(shè)置模擬級間封嚴(yán)冷氣將對渦輪級的效率、反力度等性能參數(shù)產(chǎn)生影響，也導(dǎo)致動葉通道中流動模式有一定差異。為明確哪一種進(jìn)氣邊界條件設(shè)置能夠更好地模擬渦輪級間封嚴(yán)冷氣對渦輪級的影響，下文將對主流與和封嚴(yán)冷氣流路開展耦合分析，以了解真實(shí)渦輪級間封嚴(yán)冷氣對主流的影響模式，并對源項(xiàng)方法進(jìn)行修正，給出封嚴(yán)冷氣源項(xiàng)模擬邊界設(shè)置的建議。

3 耦合分析結(jié)果

耦合分析算例是將渦輪主流計算模型與預(yù)旋噴嘴、盤前封嚴(yán)冷氣流路等多股流路進(jìn)行耦合分析，以準(zhǔn)確獲得封嚴(yán)冷氣的流動邊界。主流計算網(wǎng)格與上述算例相同，封嚴(yán)冷氣流路計算網(wǎng)格采用Work?bench劃分非結(jié)構(gòu)化網(wǎng)格。圖2給出了耦合分析的計算網(wǎng)格，網(wǎng)格節(jié)點(diǎn)總數(shù)為229萬。

圖3給出了通過級間間隙進(jìn)入動葉通道的三維流線，可以看出與上述源項(xiàng)方法給定的邊界有很大不同。由于封嚴(yán)冷氣的絕對速度很低，相對于高速旋轉(zhuǎn)的動葉，封嚴(yán)冷氣以很大的負(fù)攻角沖擊到動葉吸力面，隨后沿葉片向徑向流動卷起通道渦，這一過程甚至在動葉前緣就開始。相比之下，之前的源項(xiàng)計算結(jié)果，其封嚴(yán)冷氣大部分緊貼輪轂，直至葉片吸力面最大厚度位置才開始徑向流動。

圖4對比了無封嚴(yán)冷氣和帶封嚴(yán)冷氣流路的耦合分析計算得到的動葉根部吸力面角區(qū)的二次流動。無封嚴(yán)冷氣時，在圖4(a)所示的第三個觀察截面才在角區(qū)發(fā)現(xiàn)明顯的通道渦。而封嚴(yán)冷氣流路與主流耦合計算結(jié)果中，封嚴(yán)冷氣流沖擊到吸力面上，很快在動葉前緣形成通道渦并不斷擴(kuò)展，在第二個觀察截面已發(fā)現(xiàn)明顯的徑向流動。

圖5為動葉葉根載荷分布對比圖?？梢?，在吸力面50%軸向弦長之前的區(qū)域，兩組計算的載荷分布差異很大。帶封嚴(yán)冷氣流路的耦合計算結(jié)果大約在動葉30%軸向弦長位置出現(xiàn)了一個壓力峰值。這是由于封嚴(yán)冷氣沖擊到吸力面并滯止，導(dǎo)致沖擊區(qū)域壓力升高，葉根負(fù)荷降低。封嚴(yán)冷氣對吸力面根部的影響范圍在0～55%軸向弦長之間的葉片表面，直至喉部附近才重新加速至與無封嚴(yán)冷氣時的流動一致。這一影響將降低葉片所受的氣動力，進(jìn)而影響葉片輸出功率和渦輪級效率。

表2給出了上述不同源項(xiàng)計算給定的封嚴(yán)冷氣流動參數(shù)與帶封嚴(yán)冷氣流路的耦合計算結(jié)果的差異。帶封嚴(yán)冷氣流路耦合計算的封嚴(yán)冷氣相對氣流角為-58°，該截面動葉葉根的進(jìn)口結(jié)構(gòu)角為43°，即封嚴(yán)冷氣以101°的負(fù)攻角沖擊到動葉根部，從而導(dǎo)致上述很強(qiáng)的二次流動。根據(jù)該計算結(jié)果，繪制封嚴(yán)冷氣的速度三角形(圖6)?？梢?，目前通常采用的源項(xiàng)計算速度邊界與真實(shí)流動情況有所差異，均不能反應(yīng)真實(shí)封嚴(yán)冷氣流動對主流的影響。

聯(lián)系表1中的渦輪效率計算結(jié)果，從以上速度三角形可以看出，隨著封嚴(yán)冷氣相對氣流角從-72°到-58°、-8°、0°、45°變化，渦輪效率依次升高。由此可見，封嚴(yán)冷氣相對于動葉的流向決定了其能產(chǎn)生的影響大小。如果能較好地組織該股封嚴(yán)冷氣的流動，就可利用該封嚴(yán)冷氣輸出較大功率并抑制端壁二次流動，降低流動損失，提高渦輪效率，同時對動葉緣板進(jìn)行有效冷卻等。

4 對源項(xiàng)方法的修正對比

根據(jù)上述帶級間封嚴(yán)冷氣流路的耦合分析結(jié)果，將該組計算得到的封嚴(yán)冷氣進(jìn)入主流道的方向和流量作為邊界條件進(jìn)行帶級間封嚴(yán)冷氣的源項(xiàng)計算，并對比兩者的差異。從前文耦合分析結(jié)果可看出，封嚴(yán)冷氣進(jìn)入主流后，其與主流之間存在嚴(yán)重的剪切作用，并迅速加強(qiáng)端壁二次流動。如果在靜子域上加載封嚴(yán)冷氣，則定常計算所采用的混合面平均處理方法將無法向下游準(zhǔn)確傳遞兩者的速度差異，而該速度差異又是準(zhǔn)確模擬封嚴(yán)冷氣影響的關(guān)鍵。因此，為更準(zhǔn)確模擬封嚴(yán)冷氣與主流的相互影響，應(yīng)將封嚴(yán)冷氣加載在轉(zhuǎn)子計算域上，并且參照耦合計算結(jié)果給定較大的負(fù)攻角。表3為修正后的源項(xiàng)計算與耦合分析方法的邊界和結(jié)果對比。可見，兩者的封嚴(yán)冷氣縫出口馬赫數(shù)很接近，這也表明修正后的源項(xiàng)方法對封嚴(yán)冷氣模擬的有效性。兩組計算的效率、反力度都很接近。故可認(rèn)為，只要較為準(zhǔn)確地給定了封嚴(yán)冷氣源項(xiàng)進(jìn)口邊界，就可以較為準(zhǔn)確地模擬封嚴(yán)冷氣對渦輪性能的影響。

表2 源項(xiàng)方法與耦合分析方法封嚴(yán)縫位置流動參數(shù)對比Table 2 The contrast of flow parameter of rim seal using the source term method and coupling analysis method

表3 修正后的源項(xiàng)邊界與耦合分析方法對比Table 3 The comparison between the corrected source term and coupling flow analysis

圖7給出了兩種方法得到的封嚴(yán)冷氣流線及葉片表面極限流線對比?？梢?，當(dāng)封嚴(yán)冷氣的進(jìn)口邊界條件與帶級間封嚴(yán)冷氣流路的耦合分析一致時，模擬結(jié)果與耦合分析結(jié)果很接近。

圖8給出了兩種方法計算得出的葉片根部表面壓力分布對比?？梢?，兩種方法得到的葉片表面載荷分布整體趨勢、吸力面壓力峰值位置一致，僅在吸力面壓力峰值上略有差異。

5 結(jié)論

采用源項(xiàng)方法對四種常用封嚴(yán)冷氣邊界和無封嚴(yán)冷氣的渦輪級全三維流動模型進(jìn)行了對比計算，了解了邊界設(shè)置對計算結(jié)果的影響。同時，進(jìn)行了級間封嚴(yán)冷氣流路與主流流路耦合的渦輪級數(shù)值模擬，以考慮真實(shí)封嚴(yán)冷氣流動造成的影響。將耦合分析結(jié)果與源項(xiàng)方法計算結(jié)果進(jìn)行對比，并據(jù)此對源項(xiàng)方法進(jìn)行修正。結(jié)果表明：

(1)不同封嚴(yán)冷氣邊界條件對計算結(jié)果有一定的影響，計算效率最大偏差0.6個百分點(diǎn)；

(2)常用的封嚴(yán)冷氣邊界設(shè)置計算結(jié)果與帶封嚴(yán)冷氣流路和主流耦合計算的結(jié)果差異較大，特別是在流動模式上；

(3)封嚴(yán)冷氣在進(jìn)入主流道后，由于其軸向、切向速度低，對動葉葉根形成101°負(fù)攻角，大大加強(qiáng)了通道渦強(qiáng)度，影響了渦輪效率、葉片表面流動以及葉片表面壓力分布；

(4)封嚴(yán)冷氣對渦輪效率的影響取決于封嚴(yán)冷氣速度三角形與端區(qū)主流速度三角形的差異，兩者差異越小對渦輪性能的不利影響越小；

(5)在轉(zhuǎn)子計算域上以較大負(fù)功角形式給定封嚴(yán)冷氣的源項(xiàng)進(jìn)口邊界的模擬，更能真實(shí)反映封嚴(yán)冷氣的實(shí)際流動行為，與封嚴(yán)冷氣流路和主流耦合分析結(jié)果更為接近，可作為快速迭代設(shè)計時的封嚴(yán)冷氣進(jìn)口邊界。

[1]Dénos R，Paniagua G.Influence of the hub endwall cavity flow on the time-averaged and time-resolved aero-thermo?dynamicsofan axialHP turbine stage[R].ASME GT-2002-30185，2002.

[2]張晶輝，馬宏偉.輪緣封嚴(yán)氣體對渦輪轉(zhuǎn)子性能影響的非定常數(shù)值研究[J].推進(jìn)技術(shù)，2014，35(4)：470—478.

[3]Turner M G，Vitt P H，Topp D A，et al.Multistage simula?tions of the GE90 turbine[R].NASA/CR—1999-209311，1999.

[4]Mirzamoghadam A V，Heitland G，Morris M C，et al.3D CFD ingestion evaluation of a high pressure turbine rim seal disk cavity[R].ASME GT2008-50531，2008.

[5]Mirzamoghadam A V，Heitland G，Hosseini K M.The ef?fect of annulus performance parameters on rotor-stator cavity sealing flow[R].ASME GT2009-59380，2009

[6]Chilla M，Hodson H，Newman D.Unsteady interaction be?tween annulus and turbine rim seal flows[J].ASME GT2012-69089，2012.

[7]周 楊，牛為民，鄒正平，等.輪轂封嚴(yán)氣體對高壓渦輪二次流動的影響[J].推進(jìn)技術(shù)，2006，27(6)：515—520.

[8]胡松林，劉火星，明 磊.渦輪封嚴(yán)氣體噴射角度變化對主流的影響研究[J].燃?xì)鉁u輪試驗(yàn)與研究，2009，22(4)：25—29.

[9]張偉昊，周 穎，鄒正平，等.輪轂緣板進(jìn)口圓角對渦輪氣動性能的影響[J].熱能動力工程，2011，26(2)：140—146.

[10]孔祥林，陶加銀，馮增國，等.結(jié)構(gòu)參數(shù)對輪緣密封封嚴(yán)特性影響的數(shù)值研究[J].動力工程學(xué)報，2014，34(4)：280—285.