邱亞松，白俊強，黃 琳，朱 軍，陳迎春，李亞林，周 濤

（1.西北工業(yè)大學(xué) 航空學(xué)院，陜西 西安710072；2.上海飛機設(shè)計研究院，上海200232）

0 引 言

當代大多數(shù)的民用干線飛機均采用翼吊發(fā)動機的布局。翼吊布局具有發(fā)動機短艙進氣口不容易受機身或機翼尾流的干擾（發(fā)動機來流直接，不受飛機其他部件干擾，進氣流場好）、發(fā)動機維護方便、更輕的機翼和機身結(jié)構(gòu)，以及整個飛機的重量分布更合理等優(yōu)點［1］。當然，翼吊布局形式也有缺點，其中比較嚴重的一點是：翼吊發(fā)動機短艙會顯著惡化增升裝置的氣動效率［2］。而針對該問題的最佳解決方案就是在發(fā)動機短艙上加裝渦流片，進行流動控制。

國外對于翼吊短艙對全機增升裝置氣動性能影響及短艙渦流片設(shè)計有很好的研究與應(yīng)用，并已有豐富的工程應(yīng)用經(jīng)驗。但鮮有從機理上闡述問題本質(zhì)的文獻。早在20世紀50年代的DC－8發(fā)展項目中，為了弄清楚流動機理，道格拉斯公司對有短艙掛架、無短艙掛架構(gòu)型進行了大量風洞試驗，收集了大量的數(shù)據(jù)并進行了許多分析。在60年代末期至70年代初期設(shè)計的DC－10上率先采用了短艙渦流片［3］。近期，為了發(fā)展自己的大飛機，日本JAXA（日本宇宙航空研究開發(fā)機構(gòu)）進行了一系列的高升力研究計劃，其研究的重點內(nèi)容之一就是翼吊發(fā)動機短艙對于全機增升裝置氣動性能的影響及短艙渦流片位置的優(yōu)化。其在研究過程中大量使用了CFD技術(shù)及眾多先進的風洞試驗技術(shù)。從其發(fā)表的文獻可以看出，日本對于該方面的研究已經(jīng)相當深入［4－5］。歐洲各個相關(guān)研究機構(gòu)也在近期組織了兩期名為歐洲高升力的研究計劃（EUROLIFT），參與的研究機構(gòu)包括DLR（德國宇航研究院）、NLR（荷蘭國家宇航局）、ONERA（法國國家宇航局）、FOI（瑞典國防研究署）等。同樣，其研究的重點內(nèi)容之一就是翼吊發(fā)動機短艙對增升裝置氣動性能影響及短艙渦流片設(shè)計。并且其規(guī)模遠遠超過了JAXA的研究計劃，投入的人力、物力都相當大，采用的手段也更為先進。從公開發(fā)表的文獻來看，其取得的成果更為顯著［6－10］。國內(nèi)在這方面還未見有公開發(fā)表的研究文獻。在發(fā)展我國的大飛機的背景下，這項研究顯得極為迫切。

本文運用數(shù)值模擬的方法，對加掛翼吊發(fā)動機短艙的增升裝置流動控制進行研究。首先，運用本文的計算方法對某風洞試驗構(gòu)型進行數(shù)值模擬，并將計算結(jié)果與風洞數(shù)據(jù)進行對比，確定方法的可行性。之后，對比計算了有無翼吊發(fā)動機短艙的兩個增升裝置構(gòu)型的氣動特性，從流動機理上分析了翼吊發(fā)動機短艙惡化增升裝置氣動特性的原因。在此基礎(chǔ)上，在發(fā)動機短艙上進行了渦流片安裝位置的優(yōu)化。通過流場顯示及流動機理分析，解釋了短艙渦流片改善增升裝置氣動性能的機理并闡明了渦流設(shè)計的準則。本文設(shè)計的翼吊發(fā)動機短艙渦流片經(jīng)過了風洞試驗驗證。

1 計算構(gòu)型

本文共涉及到六個增升裝置構(gòu)型的計算，均為前緣縫翼加單縫福勒襟翼的形式，其幾何模型分別如圖1所示。其中構(gòu)型A為風洞標模，也就是用來驗證本文數(shù)值方法的構(gòu)型。構(gòu)型B與構(gòu)型A類似，為無翼吊發(fā)動機短艙的增升裝置翼身組合體構(gòu)型。構(gòu)型C為在構(gòu)型B基礎(chǔ)上加裝翼吊發(fā)動機短艙后的構(gòu)型。構(gòu)型D1、D2、D3，分別為在構(gòu)型C的翼吊短艙三個不同位置安裝了渦流片的增升裝置構(gòu)型。需要注意的是，在構(gòu)型C及D1、D2、D3上機翼前緣與掛架相接的地方，由于短艙掛架的存在，前緣縫翼在此被切除了一段。

圖1 增升裝置幾何構(gòu)型Fig.1 Geometrical configurations of high－lift system

2 網(wǎng)格劃分與計算方法

2.1 網(wǎng)格劃分策略

由于對復(fù)雜外形生成結(jié)構(gòu)網(wǎng)格比較困難，有些情況下幾乎無法生成。所以對于類似帶翼吊發(fā)動機短艙的三維增升裝置幾何構(gòu)型，普遍采用非結(jié)構(gòu)網(wǎng)格。但在單個網(wǎng)格單元最大尺度相同的情況下，非結(jié)構(gòu)網(wǎng)格的耗散比結(jié)構(gòu)網(wǎng)格要大得多，尤其是對于有強剪切流動的情況［7］。在近期日本及歐洲公布的高升力研究計劃中，數(shù)值模擬幾乎都采用了非結(jié)構(gòu)網(wǎng)格，氣動力模擬結(jié)果與試驗結(jié)果相比，絕對值誤差很大。但是構(gòu)型局部變化之后的氣動力變化量捕捉較為準確［4，6－8］。為了更為準確地捕捉流場特征，本文對于圖1中的六個構(gòu)型均采用多塊結(jié)構(gòu)網(wǎng)格進行流場計算。構(gòu)型A、B的網(wǎng)格單元數(shù)為1050萬，構(gòu)型C的網(wǎng)格單元數(shù)為1850萬、構(gòu)型D1～D3的網(wǎng)格單元數(shù)為1970萬。其中對于D1～D3三個最為復(fù)雜的構(gòu)型，如果單獨進行網(wǎng)格劃分約需40個人工時。本文借鑒了文獻［11］中的非結(jié)構(gòu)網(wǎng)格局部網(wǎng)格自動更新技術(shù)，發(fā)展了多塊結(jié)構(gòu)網(wǎng)格局部網(wǎng)格更新技術(shù)。構(gòu)型D1～D3的網(wǎng)格均在構(gòu)型C的網(wǎng)格基礎(chǔ)上進行局部更新。除短艙渦流片局部網(wǎng)格變化之外，其余部分的網(wǎng)格均不改變。運用此技術(shù)后，劃分構(gòu)型D1～D3的網(wǎng)格只需約半個人工時。圖2展示了構(gòu)型A、C及D3的局部網(wǎng)格。

2.2 數(shù)值方法

本文采用的控制方程是三維可壓縮非定常N－S方程的積分形式。其在直角坐標系下的表達式為：

上式中各個變量的具體定義見參考文獻［12］。采用格心格式的有限體積法對控制方程進行離散。設(shè)網(wǎng)格單元i，其控制體記為Ωi，定義為Ωi內(nèi)的各物理量Qi的平均值，并假設(shè)在網(wǎng)格的格心Qi＝，則方程（1）在Ωi內(nèi)可以近似地寫成如下形式：

其中?Ωi表示Ωi的邊界。上式在單元Ωi上離散為：

這里Vi是單元體Ωi的體積，Ri和分別是單元體Ωi的無粘通量和粘性通量，

其中N表示網(wǎng)格單元i的面的總數(shù)，nk，i是網(wǎng)格單元i的第k個面的外法線的單位矢量，ΔSk，i是這個面的面積，F(xiàn)k，i（Q）和Gk，i（Q）是此面上的F（Q）和G（Q）的平均值，在數(shù)值計算中，用面心處的F（Q）和G（Q）的值近似代替面上的平均值。空間離散格式為二階迎風Roe格式。時間推進方式采用LU－SGS隱式時間推進算法。該算法不像常用的隱式算法需要大量的矩陣計算和存儲空間，有利于提高計算效率。同時可以證明，對于LU－SGS方法，由于作積分裂時引入交叉項，因此在時間方向不到一階精度，但是這并不影響格式顯式部分在空間的精度，并且它在定常計算中有較好的穩(wěn)健性和收斂性，因此只在定常計算中使用它是合適的［13］。動力方程中的雷諾應(yīng)力通過k－ωSST兩方程湍流模型和自動壁面處理計算。k－ωSST湍流模型能夠較準確地模擬分離流、漩渦流等復(fù)雜流動現(xiàn)象［14］。

圖2 增升裝置網(wǎng)格Fig.2 Mesh of high－lift systems

3 計算方法可靠性驗證

增升裝置流動特性復(fù)雜，準確的氣動力計算及流場特征捕捉較為困難。為了驗證本文所采用數(shù)值方法的準確性，本文首先計算了圖1中的風洞標模A的氣動特性。在風洞試驗雷諾數(shù)、馬赫數(shù)下，本文計算所得的構(gòu)型A升力特性與風洞試驗數(shù)據(jù)對比如圖3所示。從圖中可以看出，本文計算的升力特性與試驗結(jié)果吻合良好，尤其是失速迎角及最大升力系數(shù)的捕捉相當準確。本文計算失速迎角比試驗值只小0.33°，而最大升力系數(shù)僅小0.00502。線性段的升力系數(shù)最大差值也不超過0.1。這表明本文針對增升裝置采用的數(shù)值計算方法是合理、可信的。該計算結(jié)果相比較日本及歐洲高升力研究計劃中公布的CFD 計算數(shù)據(jù)與試驗數(shù)據(jù)比較 情 況［4，6－8］而言準確 性相當高。初步分析，認為主要由于本文采用了根據(jù)高升力系統(tǒng)流動特征進行拓撲劃分的多塊結(jié)構(gòu)網(wǎng)格，尤其在各個部件尾跡重疊區(qū)垂直物面方向上的網(wǎng)格尺度很小。而國外的高升力研究計劃中主要采用了耗散性大的非結(jié)構(gòu)網(wǎng)格。根據(jù)非結(jié)構(gòu)網(wǎng)格的性質(zhì)，從其表面網(wǎng)格尺度可以推斷其空間網(wǎng)格在各個部件尾跡重疊區(qū)垂直物面方向上的尺度較大。當然具體原因與湍流模型、數(shù)值計算格式等等都有關(guān)系，具體需進一步分析。

圖3 構(gòu)型A升力特性計算值與試驗值對比Fig.3 Lift characteristic of model A comparison for prediction and experiment

4 數(shù)值模擬結(jié)果

運用前述的數(shù)值模擬方法，本文首先對比計算了構(gòu)型B、C在自由來流馬赫數(shù)M∞＝0.2，以機翼平均氣動弦長為參考長度的雷諾數(shù)Re＝2.0×107下的大迎角升力特性，結(jié)果對比如圖4所示。從圖中可以看出：在構(gòu)型C上加裝發(fā)動機短艙和掛架后全機升力特性急劇惡化，失速迎角減小了2°，最大升力系數(shù)減小了0.52。這比文獻［2］中介紹的波音737換裝大尺寸發(fā)動機帶來的增升裝置氣動性能惡化更為嚴重?？梢娨淼醢l(fā)動機短艙對全機增升裝置的升力特性有嚴重的負面影響。

圖4 構(gòu)型B、構(gòu)型C大迎角升力特性計算對比Fig.4 Predictive lift characteristic comparison for model B and model C at high angles of attack

隨后本文對在構(gòu)型C發(fā)動機短艙上加裝渦流片的三個構(gòu)型D1、D2、D3的氣動特性進行了計算。比較構(gòu)型C而言，在構(gòu)型C的失速迎角下，D1、D2、D3的升力系數(shù)變化量分別為＋0.24、＋0.08、－0.21?？梢娫诎l(fā)動機短艙加裝渦流片進行流動控制的方法對于改善由翼吊發(fā)動機短艙引起的增升裝置氣動特性惡化相當有效，但其前提是渦流片的安裝位置必須合適。

5 流動機理分析

5.1 翼吊短艙惡化增升裝置氣動性能的機理分析

圖5為構(gòu)型B、C在構(gòu)型C的失速攻角下的空間低速流動區(qū)對比圖。從圖中可以明顯看出，比較構(gòu)型B而言，構(gòu)型C主翼及襟翼上方在空間范圍內(nèi)出現(xiàn)了很大范圍的低速流動區(qū)。該低速流動區(qū)從前緣縫翼切口處開始，往后逐漸擴大，在內(nèi)襟翼上方幾乎擴展到整個內(nèi)翼段。該低速區(qū)的存在使得構(gòu)型C主翼及襟翼上表面壓力升高，從而使得升力降低。故該低速區(qū)的存在是導(dǎo)致構(gòu)型C氣動性能惡化的主要原因。

構(gòu)型C與構(gòu)型B的主要區(qū)別在于構(gòu)型C機翼下方加掛了發(fā)動機短艙，并且由于短艙掛架的存在，其前緣縫翼也被切除了一段。翼吊的發(fā)動機短艙相對于機翼來說有一個較大的前伸量。在中、大迎角下，流過短艙后面的主翼上表面的氣流首先必須繞過短艙。繞短艙流動與繞圓柱體的流動相似，在迎角不大時，繞圓柱的流動很容易發(fā)生分離。所以在迎角不大時，繞過短艙的流動極易發(fā)生分離，而分離區(qū)的低能量氣流在自由來流的帶動下，會流向主翼上表面，并與其附面層混合。其結(jié)果會使主翼上表面的附面層較正常情況下增厚很多，從而使主翼上表面出現(xiàn)圖5所示的大范圍低速流動區(qū)域。圖6顯示了構(gòu)型C在失速迎角時的表面流向摩擦力系數(shù)小于零的區(qū)域，也就是發(fā)生物面分離的的區(qū)域。從圖中可以明顯看出：發(fā)動機短艙前部上方及后部內(nèi)側(cè)上方發(fā)生了大面積的物面分離。

圖5 構(gòu)型B、C空間低速流動區(qū)域?qū)Ρ菷ig.5 Low－speed flow space bound comparison for model B and model C

圖6 構(gòu)型C物面分離區(qū)Fig.6 The surface separation area of model C

前緣縫翼在掛架處被切除了一段也對構(gòu)型C的氣動性能有一定的負面影響。根據(jù)A.M.O.Smith提出的在多段翼繞流中縫道的5個主要有利影響［15－16］，繞前緣縫翼流動所產(chǎn)生的環(huán)量使后段翼面前緣的流動減速，這就是縫翼效應(yīng)。縫翼效應(yīng)會減弱后段翼面前緣吸力峰強度，使壓力恢復(fù)要求減小，從而推遲物面上分離的產(chǎn)生，提高失速迎角。所以前緣縫翼切除一段之后，相應(yīng)有一段主翼的失速迎角就會減小。進而會影響整個增升裝置的失速迎角和最大升力系數(shù)。

5.2 短艙渦流片流動控制機理分析

圖7為構(gòu)型D1、D2、D3在構(gòu)型C失速攻角下主翼上方空間截面上的馬赫數(shù)云圖。從圖中可以看出：構(gòu)型D1中的渦流片產(chǎn)生了一個極強的漩渦，該強漩渦向原構(gòu)型C主翼上方的低速區(qū)內(nèi)注入了大量能量，從而攪動了該低速區(qū)的流動，使得低速區(qū)的范圍大大減小。故構(gòu)型D1中的渦流片流動控制效果較好。構(gòu)型D2中的渦流片也產(chǎn)生了一個較強的漩渦，但由于渦流片安裝位置比較靠前，漩渦到達主翼上方后已經(jīng)被耗散得較弱。在主翼后方該漩渦已經(jīng)弱化、破裂，沒有明顯的渦核存在。所以構(gòu)型D2中主翼上方低速區(qū)較構(gòu)型C而言范圍有所減小，但效果明顯不如構(gòu)型D1中的渦流片好。構(gòu)型D3中的渦流片產(chǎn)生的漩渦較弱，并且渦的軌跡沒有向翼根方向發(fā)生偏移，而是一直在短艙掛架的外側(cè)沿流向發(fā)展。而從圖5可以看出，構(gòu)型C主翼上方的低速區(qū)從前向后逐漸擴張主要是向翼根方向擴展，該低速區(qū)主要位于短艙掛架的內(nèi)側(cè)。故構(gòu)型D3中渦流片產(chǎn)生的弱漩渦沒能將自由流的能量注入到主翼上方的低速區(qū)。所以即便構(gòu)型D3中的渦流片能產(chǎn)生像構(gòu)型D1中那樣的強漩渦，也不能改善增升裝置的氣動性能。不僅如此，該漩渦的出現(xiàn)還增加了對自由來流的耗散，從而使得構(gòu)型D3的氣動性能較構(gòu)型C更差。

圖7 構(gòu)型D1、D2、D3主翼上方截面馬赫數(shù)云圖Fig.7 Mach number counters on sections above main wing of model D1，D2，D3

圖8為構(gòu)型D1、D2、D3在構(gòu)型C的失速攻角下主翼上方的空間低速流動區(qū)對比圖。顯然構(gòu)型D1主翼上方的低速流動區(qū)比較構(gòu)型C而言明顯減小，并且渦流片產(chǎn)生的強漩渦在主翼后方仍很明顯。而構(gòu)型D2中的低速區(qū)比較而言也有所減小，但從圖中看不出明顯的漩渦區(qū)域。構(gòu)型D3中的低速流動區(qū)域較構(gòu)型C而言不僅沿流向尺寸擴大，而且展向方向上也擴展到了掛架的外側(cè)。說明，渦流片產(chǎn)生的漩渦在主翼前方就弱化、破裂，形成了流動雜亂無章的低速流動區(qū)。

圖9顯示了構(gòu)型D1、D2、D3在構(gòu)型C的失速攻角下表面流向摩擦力系數(shù)小于零的區(qū)域，也就是發(fā)生物面分離的區(qū)域。顯然相對構(gòu)型C而言，構(gòu)型D1、D2都只在短艙上表面發(fā)生了物面分離。并且構(gòu)型D1的短艙物面分離區(qū)有所減小。而構(gòu)型D3的短艙分離區(qū)有所減小，但主翼后方出現(xiàn)了較大范圍的物面分離。這又從另一個方面說明了構(gòu)型D3上的渦流片流動控制起到了與期望相反的效果。

圖8 構(gòu)型D1、D2、D3空間低速流動區(qū)域?qū)Ρ菷ig.8 Low－speed flow space bound comparison for model D1，D2and D3

綜上所述，在對構(gòu)型C進行流動控制以改善翼 吊發(fā)動機短艙引起的增升裝置氣動性能惡化的三種渦流片嘗試中，構(gòu)型D1上安裝的渦流片效果最好，構(gòu)型D3上安裝的渦流片效果最差。所以構(gòu)型D1上的渦流片被選為風洞試驗驗證對象。風洞試驗結(jié)果表明：在構(gòu)型C上安裝構(gòu)型D1上的渦流片之后，升力系數(shù)的增量與本文計算結(jié)果吻合良好。這進一步證明了本文計算方法的可靠性。

圖9 構(gòu)型D1、D2、D3物面分離區(qū)比較Fig.9 The surface separation area comparison for model D1，D2and D3

6 結(jié) 論

（1）翼吊發(fā)動機短艙對于增升裝置的氣動性能有嚴重的負面影響。主要原因在于：在大迎角時繞過發(fā)動機短艙的氣流發(fā)生分離，分離后的氣流被自由流帶到主翼上方形成一個低速流動區(qū)從而使得增升裝置升力下降、失速迎角減小。同時，由于發(fā)動機短艙掛架使得增升裝置的前緣縫翼被切除了一段，進而使得對應(yīng)一段主翼失去了縫翼效應(yīng)，這也會使得增升裝置的氣動性能惡化。

（2）在發(fā)動機短艙安裝渦流片，采用流動控制的方法改善翼吊發(fā)動機短艙引起的增升裝置氣動性能惡化相當有效。其基本原理是：渦流片在大迎角時產(chǎn)生的強漩渦能將自由流的能量注入到主翼上方的低速流動區(qū)，減小該區(qū)域的范圍。但是渦流片的位置必須進行優(yōu)化，在合適的位置安裝渦流片能起到正面效果，而在不適當?shù)奈恢冒惭b渦流片則會進一步惡化增升裝置的氣動性能。

（3）本文中設(shè)計的效果最好的短艙渦流片經(jīng)過了風洞試驗的驗證。安裝渦流片之后的效果與本文計算結(jié)論吻合良好。

［1］黨鐵紅.翼吊布局民用飛機發(fā)動機安裝設(shè)計［J］.民用飛機設(shè)計與研究，2008，2：8－14.（DANG Tie－h(huán)ong.The engine installation of civil air－craft with wing mounted engine configuration［J］.CivialAircraftDesignandResarch，2008，2：8－14.）

［2］方寶瑞主編.飛機氣動布局設(shè)計［M］.北京：航空工業(yè)出版社，1997：1173－1175，432－433.（FANG Bao－rui.Aerodynamic configuration design of airplane［M］.Beijing：Aeronautic Industry Press，1997：1173－1175，432－433.）

［3］VAN DAM C P.The aerodynamic design of multi－element high－lift systems for transport airplanes［J］.ProgressinAerospaces，2002，38：101－144.

［4］YUZURU YOKOKAWA，MITSUHIRO MURAYAMA，TAKESHI ITO，et al.Experiment and CFD of a high lift configuration civil transport aircraft model［R］.AIAA 2006－3452.

［5］YUZURU YOKOKAWA，MITSUHIRO MURAYAMA，MASAHIRO KANAZAKI，et al.Investigation and improvment of high－lift aerodynamic performence in lowspeed wind tunnel testing［R］.AIAA 2008－350.

［6］RUDNIK R，F(xiàn)RHR H，GEYR V.The European highlift project EUROLIFTII objectives，approach，and structure［R］.AIAA 2007－4296.

［7］FRHR H，GERY V，SCHADE N，et al.CFD prediction of maximum lift effects on realistic high－lift－commercial－aircraft－configuration within the European project EUROLIFT II［R］.AIAA 2007－4299.

［8］JOCHEN WILD，JOEL BREZILLON，OLIVIER AMOIGNON，et al.Advance high－lift design by numercal methods and wind tunnel verification within European project EOROLIFTII［R］.AIAA 2007－4300.

［9］RUDNIK R.Stall behaviour of the EUROLIFT high lift configurations［R］.AIAA 2008－836.

［10］HARALD QUIX，MATTHIAS SCHULZ，JüRGEN QUEST，et al.Low speed high lift validation tests within the European project EUROLIFT II［R］.AIAA 2007－4298.

［11］YASUSHI ITO，MITSUHIRO MURAYAMA，KA－ZUOMI YAMAMOTO，et al.Efficient computational fluid dynamics evaluation of small－device locations with automatic local remeshing［J］.AIAAJournal，2009 ，47（5）：1270－1276.

［12］陳科甲，白俊強，朱軍.發(fā)動機短艙對翼身組合體跨音速氣動特性影響研究［J］.航空計算技術(shù)，2010，40（1）：63－66.（CHEN Kejia，BAI Junqiang，ZHU Jun.Research on transonic aerodynamic performance impact of the engine nacelle［J］.AeronauticalComputingTechnique，2010，40（1）：63－66.）

［13］楊明智，袁先旭，謝昱飛，等.前體渦非對稱分離機理及前緣吹氣控制研究［J］.空氣動力學(xué)學(xué)報，2009，27（2）：186－192.（YANG Mingzhi，YUAN Xianxu，XIE Yufei，et al.Numerical research on the asymmetry mechanism of fore－body vortex and active control［J］.ActaAerodynamicaSinica，2009，27（2）：186－192.）

［14］MENTER F R.Zonal two－equation K－W turbulence model for aerodynamic flows［R］.AIAA 1993－2906.

［15］SMITH AMO.Aerodynamics of high－lift airfoil systems.fluid dynamics of aircraft stalling［R］.AGARD CP 102，November 1972，p.10－1－10－27.

［16］SMITH AMO.High－lift aerodynamics［J］.JAircraft，1975，12（6）：501－30.