陳歡歡, 周正貴
(南京航空航天大學能源與動力學院,南京 210016)
航空發(fā)動機朝著高推重比、低耗油率和高可靠性方向發(fā)展,風扇/壓氣機應具有更高的效率、負荷和裕度。轉子葉尖間隙是影響壓氣機性能的主要因素,而超音壓氣機轉子的展弦比較低,同時葉尖區(qū)存在著激波、葉尖泄漏流、端壁附面層與葉片吸力面附面層等的相互作用,葉尖間隙對壓氣機性能的影響也更加嚴重。
外國學者對葉尖間隙流動展開了大量的實驗和數(shù)值研究,Lakshminarayana等[1]研究發(fā)現(xiàn)葉尖泄漏流與通道渦的方向相反,二者相互制約;Inoue等[2]對不同葉尖間隙下的壓氣機轉子展開研究,發(fā)現(xiàn)葉尖泄漏渦的形成點與機匣靜壓值最低點對應,泄漏渦運動軌跡與機匣上的靜壓斜槽相對應;Smith[3]對某低速壓氣機的研究表明:葉尖間隙的增加會導致壓氣機裕度降低;Adamczyk[4]等研究發(fā)現(xiàn)葉尖泄漏渦和通道激波對風扇失速起決定性作用,具體表現(xiàn)為前緣泄漏渦穿過激波,二者相互作用產(chǎn)生的堵塞會導致失速。
中國學者對葉尖泄漏流也進行了相應研究。韓少冰等[5]首次提出在壓氣機轉子葉尖增設小翼的理念,通過系統(tǒng)研究發(fā)現(xiàn)在低速轉子中可以選擇合適的小翼結構達到降低葉尖損失的目的;張凱晨等[6]針對某四級重復大尺度軸流壓氣機進行數(shù)值研究,比較了不同葉尖間隙下設計點與非設計點的流動差別,對葉尖損失的機理進行了詳細分析。敬榮強等[7-8]對某高負荷跨音壓氣機進行三維流動數(shù)值模擬,研究發(fā)現(xiàn)隨著葉尖間隙增大,激波后堵塞區(qū)徑向、軸向尺度增大,壓氣機壓比、效率和流量降低。周正貴等[9]根據(jù)相對運動理念,設計出考慮端壁運動的平面葉柵實驗模型并得出由于端壁由壓力面向吸力面運動會增強葉尖泄漏流流量的結論。楊晰瓊等[10]以三維模型歸納總結了跨音轉子葉尖泄漏渦和激波的相互作用機理,得出葉片吸力面正激波受泄漏渦沖擊干擾被削弱打斷而變形的結論。
對于葉尖間隙的研究,大多數(shù)局限于間隙對轉子性能的影響,對下游靜子以及壓氣機級性能的研究較少,基于此背景,選擇通過數(shù)值研究的方法對不同葉尖間隙下小流量超音壓氣機其壓氣機級、轉子以及下游靜子性能展開研究。
研究對象為由進口導向器、一級轉子和一級靜子構成的超音壓氣機,該超音壓氣機級主要性能參數(shù)如表1所示。
表1 壓氣機設計點性能
對該壓氣機在轉子葉尖間隙分別為0、0.3、0.6、0.9 mm 4種情況下進行三維數(shù)值模擬計算。
對于葉片網(wǎng)格采用O4H型結構,葉尖間隙內(nèi)網(wǎng)格均采用HO型網(wǎng)格,間隙內(nèi)網(wǎng)格層數(shù)為17層,總網(wǎng)格數(shù)為1 867 623,最小正交角為17.09°,最大延展比為2.96,最大長寬比為2 251,各葉片網(wǎng)格數(shù)均超過50×104,由數(shù)值計算結果分析可知,對于單個葉片,當網(wǎng)格數(shù)在50×104以上時,網(wǎng)格數(shù)目對流場影響較小,因此不再進行網(wǎng)格無關性驗證。葉片表面大部分區(qū)域y+均在10以下,只有葉片前緣極小部分區(qū)域大于10,對計算精度影響不大,網(wǎng)格質量基本滿足計算要求。
計算采用單通道定常計算方法,控制方程選用基于三維雷諾平均的Navier-Stokes(N-S)方程,Spalart-Allmaras(S-A)湍流模型,空間離散采用中心差分格式,利用多重網(wǎng)格以及隱式殘差光順方法加速收斂。
進出口邊界條件為:軸向進氣,來流條件為標準大氣,總壓101 325 Pa,總溫為288.15 K,出口條件為基于徑向平衡方程的靜壓分布,轉、靜交界面采用基于周向守恒法定常交界面方法。
計算時保證100%轉速不變,通過調(diào)節(jié)出口壓力得到壓氣機特性曲線,當計算殘差收斂且壓比、效率等參數(shù)穩(wěn)定且進出口流量差控制在0.1%以內(nèi)時,認為計算結果有效,數(shù)值計算發(fā)散的最后一個工況點為喘點,為節(jié)約計算時間,喘點背壓采用二分法確定。
通過不斷改變出口壓力直至壓氣機喘點,得到不同葉尖間隙下該壓氣機特性曲線如圖1、圖2所示。由壓氣機級特性曲線可知,隨著葉尖間隙的增大,該壓氣機所能到到達的最大等熵效率、最大壓比和最大流量降低,壓氣機特性線往左下方移動,壓氣機穩(wěn)定工作范圍降低。
圖1 不同葉尖間隙下流量-效率特性
圖2 不同葉尖間隙下流量-總壓比特性
對于下游靜子,隨葉尖間隙增大,靜子總壓損失特性線上抬,總壓損失增大并且低損失范圍減小,所能達到的最高背壓降低,如圖3所示。轉子葉尖間隙值對下游靜子性能造成較大影響。因此下游靜子性能變化也是不同葉尖間隙下壓氣機性能改變的主要原因。
圖3 不同葉尖間隙下靜子反壓-總壓損失特性
經(jīng)分析比較,選擇不同葉尖間隙下離最高效率點最近的同一反壓(173 kPa)作為設計點反壓,研究表明:當反壓小于173 kPa時,除0.9 mm葉尖間隙外壓氣機流量基本保持不變,轉子處于唯一攻角狀態(tài),即壓氣機流量不隨背壓變化而改變。壓氣機設計點工作狀態(tài)是評判壓氣機性能最重要的指標,下面對不同葉尖間隙下壓氣機設計點性能進行分析。
由表2可知,設計點除葉尖間隙為0外,隨著葉尖間隙增大,壓氣機級流量,總壓比、總效率、級裕度和轉子效率均有不同程度的降低,靜子總壓損失系數(shù)增加,壓氣機性能變差,且葉尖間隙越大這種趨勢越明顯。壓氣機級效率在葉尖間隙為0.3 mm時反而比無葉尖間隙時效率更高的主要原因是該葉尖間隙下轉子效率更高。
表2 不同葉尖間隙下設計點壓氣機性能
表3為不同葉尖間隙下的轉子葉尖總泄漏流流量及相對泄漏量(總泄漏量/總流量),由表3可知,隨著葉尖間隙增大,泄漏流流量、相對泄漏量增大,且泄漏流流量與葉尖間隙尺寸近似成比例關系。
表3 不同葉尖間隙下設計點泄漏流量
圖4為不同葉尖間隙下設計點轉子總壓比、效率沿葉高分布,二者表明葉尖間隙對轉子性能的影響。由圖4(a)可知,與無葉尖間隙相比,有葉尖間隙轉子在80%葉高以上,總壓比明顯下降,并且間隙越大下降幅度越大。由圖4(b)可知,葉尖間隙增加,轉子效率下降,并且葉尖間隙對效率沿葉高的影響區(qū)域明顯大于其對壓比的影響區(qū)域。同時,在葉尖間隙為0.3 mm的80%葉高附近,轉子效率比無葉尖間隙高,這是表2中該間隙下轉子效率更高的原因,在葉尖區(qū)域由于泄漏流動與通道二次流動方向相反[11],適當?shù)男孤┝鲃涌梢栽谝欢ǔ潭壬掀鸬较魅跬ǖ罍u強度,改善葉尖流動狀況的作用。因此對于該壓氣機存在最佳葉尖間隙值使得壓氣機效率達到最高。為分析不同葉尖間隙下設計點靜子性能差異的原因,對下游靜子進口切向氣流角、徑向氣流角及總壓損失沿葉高分布進行分析,如圖5所示。由圖5(a)可知,在70%葉高以上,隨著轉子葉尖間隙增大,靜子進口(轉子出口)切向氣流角增大,因此靜子進氣攻角增大,隨著轉子葉尖間隙增大,葉尖間隙泄漏流增加,葉尖處流動堵塞增大,葉尖附近軸向速度減小,轉子出口切向氣流角增大,造成靜子葉片靠葉尖處前緣載荷增大,如圖6靜子表面靜壓分布所示,而載荷變化主要是靜子前緣吸力面靜壓隨上游葉尖間隙增大而減小造成的。在80%葉高附近,葉尖間隙為0.3 mm時切向氣流角出現(xiàn)局部偏小,這是葉尖間隙下轉子效率偏高的原因。由圖5(b)可知,隨著葉尖間隙增大,轉子通道內(nèi)氣體徑向流動增強,在20%葉高以上靜子進口(轉子出口)徑向氣流角增大。葉尖間隙增大,靜子進口攻角增大,造成60%葉高以上靜子總壓損失增加;當葉尖間隙為0 mm時,在靜子60%葉高以上,總壓損失明顯小于有葉尖間隙,如圖5(c)所示。
圖4 不同葉尖間隙下轉子各參數(shù)沿葉高分布
圖5 不同葉尖間隙下靜子各參數(shù)沿葉高分布
圖6 不同葉尖間隙下下游靜子表面靜壓分布
如圖7所示,轉子葉尖間隙的大小對于進口導向器的影響較小,4種葉尖間隙下其總壓損失沿葉高分布基本重合,可認為下游轉子葉尖間隙值的大小對氣體在進口導向器中的流動幾乎沒有影響。
圖7 不同葉尖間隙下進口導向器總壓損失系數(shù)沿葉高分布
對設計點壓氣機流動特征進行分析可知,葉尖處流動受葉尖間隙值影響最為顯著。圖8為壓氣機90%葉高處S1流面馬赫數(shù)云圖。在無葉尖間隙時,通道內(nèi)激波一斜一正兩道激波構成,當葉尖間隙為0.3 mm,結尾正游波向葉片通道上游移動與斜激波融合,形成λ波。隨著葉尖間隙增大,激波向上游移動且激波角增大。因此隨葉尖間隙增大,轉子葉尖處激波耐反壓能力下降,級穩(wěn)定工作裕度下降(表1)。
圖8 不同葉尖間隙下90%葉高處相對馬赫數(shù)云圖
圖9為轉子葉尖間隙內(nèi)某高度軸向速度云圖,為方便比較,也給出相同高度下無葉尖間隙下速度云圖。黑線包圍部分為軸向速度小于0的區(qū)域,即回流區(qū)。進口斜激波和葉片表面附面層相互作用使得激波后逆壓力梯度增大,因此回流區(qū)主要出現(xiàn)在進口斜激波之后,隨著葉尖間隙的增大進口斜激波逐漸垂直,激波與葉片表面附面層作用位置向葉片前緣移動,且回流區(qū)軸向和周向作用范圍逐漸增大,因此葉尖堵塞增強,裕度降低。在葉尖間隙0.9 mm時,回流區(qū)跨過葉片,葉片前緣發(fā)生溢流轉子接近失速。
圖9 不同葉尖間隙下轉子同高度軸向速度云圖
圖10為不同葉尖間隙下轉子間隙內(nèi)同一高度處相對速度流線。由圖10可知,葉尖間隙泄漏流起始于前緣,在下游卷成泄漏渦,向下游發(fā)展過程中泄漏流與激波相互作用,使得激波形狀發(fā)生輕微改變。一部分泄漏流直接由尾緣流出葉片通道;還有一部分穿過相鄰葉片間隙形成二次泄漏[12]。隨葉尖間隙增大,二次泄漏區(qū)起始位置向葉片前緣移動:在葉尖間隙為0.3 mm時,二次泄漏造成相鄰葉片吸力面后段產(chǎn)生較大的回流區(qū),如圖9(b)所示,隨著葉尖間隙進一步增大,激波向上游移動,與葉尖間隙泄漏流相互作用,使葉片前緣泄漏流向上游偏移,在葉片前段吸力面產(chǎn)生回流區(qū)、二次泄漏產(chǎn)生的回流區(qū)變大,如圖9(c)、圖9(d)所示。
圖10 不同葉尖間隙下間隙內(nèi)相對速度流線
圖11為下游靜子吸力面極限流線,由圖11可知,由于圖5(a)、圖5(b)所示隨著葉尖間隙增大,靜子進口葉尖附近切向和徑向氣流角增大,靜子葉尖角區(qū)分離增大。靜子葉尖角區(qū)分離增大造成流道堵塞不僅會影響靜子工作性能,也會造成上游轉子葉尖激波前移,級穩(wěn)定工作裕度減小。
圖11 不同葉尖間隙下下游靜子吸力面極限流線
以小流量超音帶導葉壓氣機級為研究對象,分析了在轉子葉尖間隙在0、0.3、0.6、0.9 mm情況下的氣動性能與流場結構,得出以下結論。
(1)葉尖間隙增大,壓氣機總壓比、等熵效率和穩(wěn)定工作裕度下降,下游靜子低損失范圍降低;泄漏流流量與間隙尺寸近于呈線性關系。
(2)葉尖處泄漏流與通道渦相互作用,存在最佳葉尖間隙值使得該壓氣機性能達到最高。
(3)葉尖間隙增大,激波向上游移動,與葉尖間隙泄漏流相互作用,使葉片二次泄漏流向上游偏移,在葉片前段吸力面產(chǎn)生回流區(qū),二次泄漏產(chǎn)生的回流區(qū)變大。
(4)隨著葉尖間隙增大,靜子進口葉尖附近切向和徑向氣流角增大,靜子葉尖負荷增加且角區(qū)分離增大,靜子性能下降。
(5)進口導向器內(nèi)流動基本不受下游轉子葉尖間隙值的影響。