基于PAREADJ方法的航空翼型的氣動力分析研究

胡志明 樊艷紅

摘 要：提出了一種用于航空翼型的氣動力優化的PAREADJ方法，采用該方法可以優化得到氣動性能更佳的航空翼型，為機翼的設計提供理論方法。首先，由函數Rosen Brock的求解，說明PAREADJ方法的正確性；然后，以NACA-63-215翼型的反設計為例，以采用RANS方法計算的該翼型的氣動力結果作為目標及約束，驗證該方法的可行性和有效性。

關鍵詞：航空翼型；氣動力； PAREADJ

中圖分類號：V211.1 文獻標識碼：A 文章編號：2095-7394（2018）02-0036-04

航空翼型作為飛機機翼剖面，是構成飛機外形的基本要素，其氣動力特性與優化是機翼的基礎和核心技術。其氣動外形優化設計方法主要包括全局優化方法和局部優化設計方法等。國內在該方面的研究成果[1-9]，他們求解歐拉方程或者N-S方程得到流場信息，求解Ad joint方程得到梯度，根據流場信息和梯度方向進行優化，得到優化結果。劉俊等[10]進行了基于Kriging模型的翼型反設計研究。楊梅花、夏露等[11]進行了基于粒子群和人工蜂群混合算法的氣動優化設計研究。本文提出一種基于Pareto優化設計方法和Ad joint優化方法的新型混合優化方法——PAREADJ方法，并通過Rosen Brock函數和 NACA0012航空翼型的反設計對 PAREADJ方法進行驗證分析。

1 PAREADJ方法

PAREADJ方法是在基于Pareto優化設計方法和Adjoint優化方法兩種優化方法的基礎上提出的新型混合優化方法，它綜合了Pareto優化設計方法的全局性和Adjoint優化方法的局部性的優點，其流程圖如圖1。

2 算例

例1 Rosen Brock算例

該算例的優化模型如下：

[min：f（x）=k=1ny-1100（xk+1-x2k）2+（xk-1）2，s.t. -2xk2，k=1，2，…，nv。]

該算例的解析解是[x?=（1，…，1），f（x?）=0]。圖2和圖3分別表示[nv=2]時目標函數的軌跡及收斂歷程，其中，第1、3、8、12、16迭代步數的解（0.05，0）、（0.394 82，0.140 14）、（0.709 14，0.490 91）、 （0.911 28，0.828 91）、（1.000 00，1.000 00）；該算例從理論上說明：PAREADJ算法優于Pareto方法與 Adjoint方法，收斂速度更快。

例2 優化反設計

15%相對厚度航空翼型的優化反設計

下面采用N-S方程數值模擬翼型粘性繞流，對15%相對厚度的NACA-63-215翼型進行優化反設計，計算網格如圖4所示。采用Hicks-Henne參數化方法，共64個設計變量。這里，Ma=0.165，α=5.0°， Re=3.0×106。設計過程中升力系數保持不減，阻力系數保持不增。采用PAREADJ算法進行該翼型的優化反設計研究。

圖5給出了采用PAREADJ算法設計結果與目標值的比較（（a）幾何外形；（b）壓力分布），可見，設計的翼型和壓力分布均與目標值吻合。圖6表示目標函數的收斂歷程，從圖中可以看出：PAREADJ算法優于Pareto方法與 Adjoint方法，收斂速度更快。

3 結論

由以上算例，得到結論如下：

（1）通過Rosen Brock算例，驗證了該文的算法是正確的。

（2）該算法可以用于翼型的優化設計。

（3）該算法優于Pareto方法和Adjoint方法。

參考文獻：

[1] 唐智禮，黃明恪.基于控制理論的Euler 方程翼型減阻優化設計[J].空氣動力學學報，2001，19（3）：262-270.

[2] 周鑄，陳作斌.基于N-S 方程的翼型氣動優化設計 [J].空氣動力學學報，2002，20（2）：141-149.

[3] 楊旭東，喬志德，朱兵.基于控制理論和NS方程的氣動設計方法研究[J].空氣動力學學報，2005，23（1）：46-51.

[4] 楊旭東，喬志德，朱兵.氣動/幾何約束條件下翼型優化設計的最優控制理論方法[J].計算物理，2006，23（1）：66-72.

[5] 熊俊濤，喬志德，楊旭東，等.基于黏性伴隨方法的跨聲速機翼氣動優化設計[J].航空學報，2007，28（2）：281-285.

[6] 左英桃，高正紅，詹浩.基于N-S 方程和離散共軛方法的氣動設計方法研究[J].空氣動力學學報，2009，27（1）：67-72.

[7] 王龍，宋文萍，楊旭東.一種基于伴隨方程攻角不變的多目標優化減阻方法[J].西北工業大學學報，2012，30（1）：68-72.

[8] HAN Zhong-Hua. Improving Adjoint-Based Aerodynamic Optimization via Gradient-Enhanced Kriging[C].50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition，2012-0670，2012.

[9] 樊艷紅，宋文萍，韓忠華. 提高基于Adjoint方法翼型優化設計魯棒性的研究[J]. 西北工業大學學報， 2013，31（4）： 547-555.

[10] 劉俊，宋文萍，韓忠華，等. Kriging模型在翼型反設計中的應用研究[J].空氣動力學學報，2014，32（4）： 518-526.

[11] 楊梅花，夏露， 張欣，等.基于粒子群和人工蜂群混合算法的氣動優化設計[J].航空工程進展， 2017，8（2）： 182-189.

The Research of Aerodynamic Analysis of Aviation Airfoil Based on the PAREADJ Method

HU Zhi-ming， FAN Yan-hong

（School of Applied Mechanics， Taiyuan University of Science and Technology， Taiyuan 030024， China）

Abstract： The PAREADJ method of aviation airfoil aerodynamic optimization is put forward， with which we can optimize by aviation airfoil for better aerodynamic performance. In this way， we can provide accurate and reliable geometric shape and aerodynamic data for the design of wing. First of all， take the RosenBrock analytic function for example， we show that the optimization method is correct. Then， for the inverse design of NACA0012 airfoil， the aerodynamic result by the RANS method is as the target and constraints， the feasibility and effectiveness of the method is verified.

Key words： aviation airfoil； aerodynamic；PAREADJ

責任編輯 祁秀春