低雷諾數下翼型后緣變化的氣動特性研究

吳書山, 周洲, 甘文彪, 許曉平

(西北工業(yè)大學 無人機特種技術重點實驗室, 陜西 西安 710065)

低雷諾數下翼型后緣變化的氣動特性研究

吳書山, 周洲, 甘文彪, 許曉平

(西北工業(yè)大學 無人機特種技術重點實驗室, 陜西 西安 710065)

以低雷諾數下高空長航時無人機翼型為背景,研究了翼型后緣變化的氣動特性。以低雷諾數流動中氣動性能良好的SD7062翼型為研究對象,采用求解二維雷諾平均N-S方程的有限體積法,選用k-kl-w湍流轉捩模型,給出了數值方法與實驗數據的對比驗證。對翼型后緣變化的流場特性進行了細致分析,總結了后緣變化對翼型氣動特性影響的變化規(guī)律,為低雷諾數流動中翼型設計提供參考。

低雷諾數; 翼型后緣變化; 分離氣泡; 流動轉捩; 氣動特性

引言

高空長航時無人機大都采用大展弦比機翼,其翼型設計是氣動設計的主要問題。在超高空(>20 km)大氣密度極其稀薄,故低雷諾數流動問題十分突出。在低雷諾數復雜流動中有效提升翼型的氣動特性,對高空長航時無人機至關重要,采用后緣加厚和后緣轉角等后緣變化的方法是一種有效手段。目前針對后緣變化對翼型氣動特性影響的研究相對較少,冉景洪等[1-2]對翼型在低雷諾數流動中相對厚度和彎度的影響進行了研究;劉杰平等[3]研究了后緣加厚方式對典型風力機翼型氣動性能的影響。他們的工作都只是進行氣動特性的分析,沒有涉及流場的分析和后緣轉角的影響。

本文選取在低雷諾數流動中氣動性能良好的SD7062翼型為研究對象,對翼型后緣變化的流場特性進行了細致分析,總結了后緣變化對翼型氣動特性影響的變化規(guī)律,為低雷諾數流動中翼型的設計提供參考。

1 數值方法

本文計算采用求解二維雷諾平均N-S方程的有限體積方法,湍流模型采用Transitionk-kl-w[4]模型。在低雷諾數流動中時間推進方案難以收斂,因此采用預處理矩陣減輕方程的“剛性”。

1.1 湍流模型

分離氣泡是低雷諾數下翼型表面流動的重要特征之一,其生成機理與氣流的轉捩特性有很大關系,故采用轉捩模型可以較好地捕捉到這一流動特征。本文采用k-kl-w湍流轉捩模型,此方程包括三個輸運方程[4]:

(1)

(2)

(3)

1.2 離散格式與邊界條件

控制方程中對流通量項采用二階精度Roe迎風格式離散;粘性通量項采用中心差分格式離散;時間離散采用隱式LU-SGS方法。翼型表面滿足無滑移邊界條件,遠場采用壓力遠場。

2 計算模型與計算條件

本文選取在低雷諾數流動中氣動性能良好的SD7062翼型為研究對象,以提升翼型的氣動特性為研究目的,采用后緣加厚和后緣轉角的后緣變化方法進行翼型后緣變化的氣動特性研究。具體方法為:后緣加厚采用吸力面和壓力面同時變化,加厚厚度值從后緣到前緣某一位置處呈線性遞減直至為零的方法;后緣轉角以旋轉位置處中弧線上的點為旋轉中心,從后緣至旋轉位置處吸力面和壓力面同時繞旋轉中心旋轉某一角度,偏轉時以角度向下為正,旋轉位置處采用光滑過渡。由于修改模型較多,為了進行區(qū)分,采用以下命名方法:

(1)0.26-0.4%:表示在距離翼型前緣0.26c位置處開始加厚,后緣加厚厚度值0.4%c;

(2)0.26-0.4%-0.25-2°:表示在0.26-0.4%翼型的基礎上從0.25c位置處開始后緣向下偏轉2°。具體計算模型如圖1所示。

圖1 計算模型示意圖

參考美國伊利諾伊大學(UIUC)低湍流亞聲速風洞進行SD7062翼型風洞實驗時的雷諾數范圍[5]。為了盡可能真實地模擬飛機巡航時所處的流場環(huán)境,計算選取與實驗條件一致的Re=299 500狀態(tài)。計算結果與實驗結果的對比如圖2所示。

圖2 升力系數和阻力系數隨迎角的變化曲線

由圖2可知,上述數值方法模擬所得的計算結果與實驗值吻合得很好,說明這一數值計算方法適用于本類物理問題的求解。

3 基本流場特性分析

低雷諾數下翼型表面流動情況復雜,其典型特征是存在分離氣泡和流動轉捩。為了深入研究翼型后緣變化的氣動特性,需要對翼型后緣變化的基本流場進行分析。選取4.17°迎角時的SD7062翼型、后緣加厚的0.26-0.4%翼型和后緣轉角的0.26-0.4%-0.25-2°翼型進行綜合對比研究。圖3為三種翼型表面摩擦阻力系數(Cf)分布。

圖3 翼型表面摩阻系數圖

由圖3可知，在翼型上表面0.3c位置處Cf出現一段平臺區(qū),分離氣泡開始出現,由于此處粘性和速度梯度都很小,因此Cf接近于零且基本不變;由平臺區(qū)向后,氣泡發(fā)展迅速,產生了分離渦,逆壓梯度急劇增大,出現了Cf負值區(qū);此后流動發(fā)生轉捩,粘性急劇增大,壓力梯度逐漸變化直至轉為順壓,流動發(fā)生了湍流再附;三種翼型中后緣加厚使得Cf的最低點提高,后緣偏轉使得平臺區(qū)和轉捩位置提前。

圖4為翼型上表面速度矢量圖和湍動能圖。速度矢量圖表征了上表面流動速度和壓力梯度的變化,清晰地反映了流動速度由層流型變化到分離回流型并最終過渡到湍流形態(tài)的過程。湍動能圖反映了上表面由層流轉捩到湍流的流動發(fā)展特征,在回流速度型出現之后,湍動能變化逐漸變得迅速。綜合來看,分離的基本要素是粘性和壓力梯度,速度型深刻地反映了壓力梯度的變化;轉捩的一種重要形式是分離流轉捩,湍動能的變化恰好體現了這種轉捩形式的發(fā)展情況。

圖4 翼型上表面速度矢量圖和湍動能圖

總體來說,通過數值模擬,翼型表面的分離氣泡和分離流轉捩特性得到了捕捉和細致分析。

4 后緣變化后的氣動特性分析

由上面的流場特性分析可知,后緣加厚和后緣轉角對低雷諾數下翼型表面流動的分離、轉捩等特性產生了非常明顯的變化,這種變化在翼型氣動力方面必然會引起很大的不同,所以有必要對后緣變化后的翼型氣動特性進行分析。

4.1 后緣加厚對翼型氣動特性的影響分析

分離氣泡隨著迎角的增大逐漸向前緣移動,為了使后緣加厚對分離氣泡產生直接的影響,同時又不引起原翼型氣動力特性的較大變化,選取靠近前緣的0.26c處改變厚度,厚度值取0.4%c和1.0%c,其氣動力特性如圖5所示。

圖5 后緣不同加厚厚度的氣動特性

從圖中可以看出,隨著厚度的增加升力線斜率和阻力系數都逐步提高,其中1.0%c厚度升力和阻力提高最多。在極曲線中,最大升阻比和最大升阻比所對應的升力系數略微有些下降。總體而言,變化不是特別明顯。由此可以看出,后緣小幅度加厚時不會帶來翼型氣動特性的明顯變化,這樣既可以保證飛機在制造過程中所需要的后緣厚度,同時又不會使得其氣動性能損失太多。

4.2 后緣轉角對翼型氣動特性的影響分析

4.2.1原翼型的基礎上偏轉不同角度

為使后緣轉角時對翼型表面的分離氣泡產生直接影響又小幅度改變翼型彎度,偏轉位置選取0.25c位置處,角度范圍為-4°～6°,以SD7062翼型為研究對象進行偏轉,不同偏轉角時氣動力的差別如圖6所示。

圖6 SD7062翼型偏轉不同角度的氣動特性

從圖中可以看出，隨著轉角的增大,升力系數曲線整體向上平移,阻力系數在負迎角時下降正迎角時增加;最大升阻比和最大升阻比所對應的升力系數都隨轉角的增大而逐步提高。

4.2.2后緣加厚的基礎上偏轉不同角度

以0.26-0.4%翼型為研究對象,參照前面的方法,偏轉位置選取0.25c位置處,角度范圍為-4°～6°,氣動力差別如圖7所示。

從圖中可以看出，在加厚的基礎上進行偏轉時,隨著轉角的增大,升力系數曲線整體向上平移,阻力系數在負迎角時下降、正迎角時增加;最大升阻比和最大升阻比所對應的升力系數都隨轉角的增大而逐步提高。

圖7 0.26-0.4%翼型偏轉不同角度的氣動特性

5 結論

通過本文研究得到以下結論:

(1)低雷諾數流動中,后緣加厚可以提高翼型升力線斜率,增大阻力系數;小幅度加厚時既可以保證飛機在制造過程中所需要的后緣厚度,同時又不會使氣動性能損失太多。

(2)后緣向下轉角時,升力系數曲線整體向上平移,阻力系數在負迎角時下降、正迎角時增加;最大升阻比和最大升阻比所對應的升力系數隨轉角的增大而逐步提高。

(3)合理的后緣變化方式可以使翼型在設計點處的氣動特性得到有效提升。

[1] 冉景洪,劉子強,白鵬.相對厚度對低雷諾數流動中翼型動態(tài)氣動力特性的影響[J].空氣動力學學報,2008,26(2):178-185.

[2] 冉景洪,劉子強,白鵬.相對彎度對低雷諾數流動中翼型動態(tài)氣動力特性的影響[J].計算力學學報,2010,27(1):88-94.

[3] 劉杰平,陳培,張衛(wèi)平.后緣加厚方式對典型風力機翼型氣動性能的影響[J].太陽能學報,2009,30(8):1092-1096.

[4] Walters D K,Cokljat D.A three-equation eddy-viscosity model for Reynolds-averaged Navier-Stokes simulations of transitional flows [J].Journal of Fluids Engineering,2008,130(1):1-14.

[5] Lyon C A,Broeren A P,Giguere P,et al.Summary of low-speed airfoil data [M].SoarTech Publication,Virginia Beach,Virginia,USA,1997.

(編輯：崔立峰)

FlowcharacteristicsofchangingthetrailingedgeatlowReynoldsnumber

WU Shu-shan, ZHOU Zhou, GAN Wen-biao, XU Xiao-ping

(National Key Laboratory of Science and Technology on UAV, NWPU, Xi’an 710065, China)

Based on high altitude and long endurance (HALE) UAVs applications, numerical computations of airfoil were conduced to analyze the aerodynamic characteristics of changing the trailing edge at low Reynolds number. The SD7062 airfoil, with good aerodynamic performance at the low Reynolds number, was selected. Finite volume method and transitionk-kl-wturbulence model were used to solve the 2D Reynolds-averaged Navier-Stokes equations. Verification was done between the numerical results and experimental results. This paper conducted a detailed analysis of flow characteristics of the trailing edge changes, and summarized the variations of the aerodynamic characteristics to provide useful references for low Reynolds number airfoil design.

low Reynolds number; airfoil trailing edge change; separation bubble; transition; aerodynamic characteristics

V211.3

A

1002-0853(2012)06-0494-04

2012-03-15;

2012-07-17; < class="emphasis_bold">網絡出版時間

時間:2012-11-23 14∶03

吳書山(1986-)，男，山東臨沂人，碩士研究生，主要研究方向為飛行器設計、計算流體力學；

周洲(1966-)，女，湖南長沙人，教授，博士生導師，主要研究方向為無人機總體、氣動布局設計。