高負荷吸附式風扇氣動性能數值模擬

周錚，王掩剛(.中國燃氣渦輪研究院，四川成都60500；.西北工業大學，陜西西安7007)

高負荷吸附式風扇氣動性能數值模擬

周錚1，王掩剛2
(1.中國燃氣渦輪研究院，四川成都610500；2.西北工業大學，陜西西安710072)

采用NUMECA商用軟件，模擬了葉片附面層吸氣技術對高負荷風扇葉片三維流場的影響。重點研究了抽氣模型的建立方法，深入分析了不同開孔位置、尺寸及不同抽吸量對風扇氣動性能的影響。結果表明，在來流條件相同的情況下，轉、靜子選取適當的開孔位置，抽吸流量為總流量的1.5%時，壓比、效率最大可分別提高6.5%和1.6%。

高負荷吸附式風扇；附面層吸氣；抽氣模型；開孔位置；高推重比；氣動性能

1 引言

國內外研究表明：吸附式風扇葉片技術是提高風扇級負荷的一項切實可行的實用性新技術，是減少級數、減輕重量、縮短長度、提高效率的有效途徑。該技術通過在葉片吸力面上強逆壓梯度處開孔，吸除即將發生大分離的附面層，使氣流在葉型后段能完全附著，實現葉片氣流彎角的增加，從而提高風扇級的氣動性能。該技術是開展下一代推重比15～20發動機壓縮部件設計的關鍵技術，而且可用于現役和在研發動機的改進，是提升我國發動機整體設計水平的有效途徑。

本文主要針對高負荷吸附式風扇轉子葉片吸力面開孔進行了三維數值模擬和分析，研究了該方法對高負荷吸附式風扇葉片內部流動模擬的可行性與可靠性，詳細分析了附面層抽吸技術對高負荷跨聲速風扇葉片三維流場及風扇氣動性能的影響。

2 數值模擬方法

采用NUMECA商業軟件，模擬了葉片附面層吸氣技術對高負荷風扇三維流場的影響。采用時間追趕的有限體積法求解三維定常粘性雷諾時均N-S方程。空間離散采用中心差分法，時間離散采用四階龍格-庫塔法，并采用隱式殘差光順和多重網格技術加速收斂。湍流模型采用Spalart-Allmaras模型[1]。

吸附式風扇的計算網格劃分采用多塊網格分區技術[2]。網格拓撲在葉片壁面區域采用O型網格，等距分布；葉片前、后延伸段采用H型網格；機匣壁面與轉子葉片之間選取0.2mm間隙；頂部間隙區采用蝶型網格拓撲結構(O型網格內嵌I型網格，可保證間隙內具有良好的正交性)，并保持頂部間隙幾何形狀不變；葉頂間隙內展向網格均勻分布，機匣附近無需加密。距固體壁面第一層網格滿足0＜Y+≤10，確保能計算出壁面的附面層狀況。計算網格結構共分12個網格塊，總網格節點數約41萬。網格細節分布見圖1。

圖1 三維計算網格圖Fig.1 Three-dimension computationalgrid

3 簡潔抽氣模型建立方法

相關資料顯示[3～5]，抽氣模型的建立過程復雜。因此，針對附面抽吸技術應用于葉片吸力面，及對流場改善的影響，本文提出了建立簡潔抽氣模型的方法和步驟。

本文所選擇的抽氣模型能準確模擬并描述附面層抽吸技術應用于葉片吸力面后，對激波位置、通道流動狀況、附面層分離方面改善的影響。

抽氣模型通過NUMECA軟件的FINE-TURBO模塊建立。步驟如下：

(1)在FINE-TURBO模塊的parameter選項中選擇cooling/bleed模型。特別指出：對于吸附式壓氣機技術，附面層抽取系統可看做一個自身抽吸的系統[6，7]，無需采用吸氣泵，抽吸出的流體從轉子葉片中流出。

(2)根據葉片吸力面流動確定大分離發生處的徑向、軸向及展向位置，根據具體需要建立若干個相匹配的吸氣模型。根據分離情況設置不同抽吸孔數量(即可在一個分離處建立一個或并列的幾個抽吸模型)。此外，考慮到網格生成的難度和質量，采用沿網格線開孔的方法，即縫的寬度為1～2個網格寬，長度則沿徑向取相同流向的若干網格。

(3)選擇抽吸孔的開孔形式。開孔形式有：抽吸孔(line)、抽氣點(point)、抽氣縫(slot)三種，本文選擇line形式。line形式在幾何學和流動特性概念下，表示為一系列抽吸孔均勻分布為一條線。這種開孔形式適用范圍較廣，且開孔位置自動均勻分布，有利于準確捕捉流場及開孔位置處的細微流動變化。

(4)選擇抽吸孔所在固面，根據實際數值模擬對象，定義抽吸孔數量(1～n)、抽吸孔直徑及抽吸流量(一般為總流量的0.5%～3.0%)[8]。抽吸孔幾何位置可由輸入具體抽吸位置的X、Y、Z坐標確定。這種方法準確，即使葉片彎扭度很大也可準確選取抽吸處，尤其適用于分離嚴重的角區。也可選擇圓柱坐標系或網格點坐標給定，或在幾何體上手動選取(選取初始點和結束點位置即可)，這種方法選取抽吸位置直觀，適用于彎扭度不大的靜子葉片。也可幾種方法結合，共同確定抽吸位置。

4 抽氣位置選取原則及方法

文獻[9]的研究表明，吸附式風扇葉片抽吸孔位置的選取在研究中最為關鍵。本算例選擇在轉子葉片吸力面開孔，開孔位置由軸向位置和徑向位置共同確定。

在葉片吸力面開孔，須特別注意激波流動狀況，必須選擇在激波之后的壓力恢復區并沿附面層流動方向開孔。

抽吸孔位置的設置通過NUMECA軟件的FINE-TURBO模塊實現。該開孔模型可設置吸氣模型結構、抽吸孔數量及吸氣流量等參數。若同一位置處抽吸量較大，可同時開設兩排及以上的抽吸孔，以便在抽吸孔直徑一定的情況下，保證抽吸流量和吸氣效果。開孔效果見圖2。

圖2 轉子葉片開孔位置示意圖Fig.2 Rotor suction slot location

圖3 75%葉高擬S1流面相對馬赫數等值線對比圖Fig.3 Distribution of S1 stream surface relativeMacHnumberat75%span witHdifferentaspirated flow rates

表1 抽吸孔參數列表Table 1 Details of suction slotdistribution

5 全三維數值模擬結果與分析

以某單級高負荷吸附式風扇為研究對象，利用NUMECA數值模擬軟件，研究不同轉、靜子葉片位置及不同吸氣量對風扇氣動性能的影響。

5.1 吸氣位置及吸氣量計算分析

吸氣位置的選取方法：利用三維分析軟件對未進行抽吸處理的風扇進行數值模擬。對于超聲葉型，根據激波與邊界層相互作用情況，選取合適的吸氣位置；對于亞聲葉型，則根據可能存在的附面層分離情況，選取合適的吸氣位置。

本算例選擇的開孔抽吸位置：轉子葉片40%～95%葉高、5%弦長位置，轉子葉片全葉高、80%弦長位置，靜子葉片全葉高、30%弦長位置。

計算中轉子葉片抽吸孔參數設置見表1。

多組對比計算結果表明，開孔直徑必須隨吸氣量的增加而增大，否則計算無法正常進行。大致規律為：吸氣流量為0.01～0.03 kg/s，孔直徑1mm；吸氣流量為0.03～0.05 kg/s，孔直徑2mm；吸氣流量為0.06～0.10 kg/s，孔直徑3 mm；吸氣流量為0.12～0.20 kg/s，孔直徑4mm。

5.2 抽吸前后內部流場細節對比分析

按上述吸氣位置和吸氣量，對該風扇進行設計點的全三維分析，并與無吸氣狀態進行對比。圖3、圖4分別為75%和95%葉高吸氣前后的擬S1流面相對馬赫數等值線對比圖，圖5為吸氣前后的擬S2流面馬赫數對比云圖，圖6為轉子吸力面極限流線對比圖。

從圖3、圖4中可看出，在無吸氣條件下，存在弓形激波和通道激波，95%葉高處葉型尾緣出現較大分離流動，分離起始位置位于激波后。抽吸掉即將發生大分離的低能高熵流體，較好地改善了通道流通狀況，通道激波強度明顯減弱，在75%葉高處改善效應最明顯。從圖5和圖6中看，由于抽吸孔1的抽吸，使得前緣弓形激波位置向下游漂移，損失增大，高馬赫數區增大。抽吸孔2的抽氣效果明顯，轉子葉片吸力面馬赫數分布更均勻，吸氣位置前附面層內流體動能增加，有助于避免低能流體過早分離。

圖4 95%葉高擬S1流面相對馬赫數等值線對比圖Fig.4 Distribution of S1 stream surface relativeMacHnumber at95%span witHdifferentaspirated flow rates

圖6 轉子葉片吸力面表面極限流線分布Fig.6 Distribution ofstream linewitHdifferentaspirated flow rates

吸氣位置后，激波結構的改變使得分離區后移并減小，通道擴壓能力得到恢復，同時也改善了靜子葉片排的流場分布，提高了整個級的壓比和效率，但轉子葉片角區出現了低速區。因此，可在角區開孔抽吸，進一步改善流動狀況。

5.3 抽吸前后總性能對比分析

為進一步驗證和研究吸附式風扇的性能，對相對換算轉速1.00、0.95、0.90的各工況點，針對相同引氣量進行了總性能參數計算。各相對換算轉速下的特性線如圖7所示。

圖7 抽吸前后吸附式風扇特性曲線Fig.7 Characteristicsof fan witHaspiration

從特性圖看出，在設計點工況，吸氣后風扇效率和壓比均大幅提高，效率最大提高1.6%，壓比最大提高6.5%；在非設計點工況，吸氣后風扇壓比和效率較無吸氣情況有較大提升，同時裕度也有一定改善。

6 結論

(1)本文建立的抽氣模型簡潔、合理，易于調整開孔位置及開孔數量，能準確模擬計算對象流場內部流動狀況，可滿足吸附式風扇數值模擬要求。

(2)在風扇轉、靜子葉片吸力面開孔吸除適量的低能高熵流體，可有效控制葉片尾緣附面層的大分離流動，加強流道通流能力，提高葉片負荷、風扇級效率和壓比。

(3)在設計點工況，在轉子葉片吸力面80%弦長、靜子葉片吸力面30%弦長、全葉高開孔，且吸氣流量為總流量的1.5%時，吸附效果最為顯著，壓比和效率最大可分別提高6.5%和1.6%；在非設計點工況，吸氣對風扇性能的提升，具有相同的效果。

[1]陶文銓.數值傳熱學[M].西安：西安交通大學出版，2001.

[2]盧新根.軸流壓氣機內部流動失穩及其被動控制策略研究[D].西安：西北工業大學，2007

[3]張華良，王松濤，王仲奇.采用壁面吸氣改善葉柵性能的數值模擬[J].動力工程，2006，26(6)：795—798.

[4]周楊，鄒正平，劉火星，等.邊界層吹吸氣對高負荷擴壓葉柵性能的影響[J].推進技術，2007，28(6)：647—652.

[5]宋彥萍，陳浮，趙桂杰，等.附面層吸除對大轉角壓氣機葉柵氣動性能影響的數值研究[J].航空動力學報，2005，20(4)：561—562.

[6]Merchant A A，Drela M，Jack L.etal.A Family of Designs for Aspirated Compressor[R].ASME 98-GT-196，1998.

[7]Merchant A A，Drela M，Jack L.et al.Aerodynamic De?sign and Analysis ofa HigHPressure Ratio Aspirated Com?pressor stage[R].ASMEGT2000-619，2000.

[8]王松濤，潛紀儒，馮國泰，等.壁面吸氣抑制分離減少流動損失的研究[J].工程熱物理學報，2006，27(1)：48—50.

[9]王掩剛，程榮輝，蘭發祥，等.吸附式葉柵抽吸流與激波相干性研究[J].燃氣渦輪試驗與研究，2008，21(2)：15—18.

Numerical Simulation of Aerodynamic Performance for High-Loading Aspirated Fan

ZHOU Zheng，WANG Yan-gang
(1.China Gas Turbine Establishment，Chengdu 610500，China；2.Northweastern Polytechnical University，Xi’an 710072，China)

The numerical simulation of a high-loading aspirated fan was introduced.NUMECA software was adopted to simulate the influence of high-loading fan witHboundary layer suction on the three-dimen?sional flow field.Themethod of building the suction model and the influence of suction location and suction flow rate on the performance of fan were investigated in details.The results showed thatwitHthe same flow and appropriate suction locations for stator and rotor,the pressure ratio and adiabatic efficiency could be im?proved to 6.5%and 1.6%respectivelywitH1.5%suction ratio.

high-loading aspirated fan；boundary layer suction；suctionmodel；suction slot location；higHthrust-to-weight ratio；aerodynamic performance

V231.3

A

1672-2620(2012)04-0016-04

2012-02-02；

2012-08-28

航空基金(2011ZB24002)

周錚(1982-)，女，重慶人，工程師，碩士，主要從事航空發動機壓氣機氣動設計和研究。