吸氣式高超聲速飛行器大迎角氣動特性分析

羅文莉， 李道春， 向錦武

北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院， 北京 100191

吸氣式高超聲速飛行器大迎角氣動特性分析

羅文莉， 李道春， 向錦武*

北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院， 北京 100191

吸氣式高超聲速飛行器在飛行過程中受到大氣紊流等外部干擾的作用時，飛行姿態(tài)很可能會出現(xiàn)大迎角情況。針對大迎角飛行時飛行器可能出現(xiàn)的氣動問題，對一種典型吸氣式高超聲速飛行器的流場進(jìn)行了數(shù)值模擬。以雷諾平均Navier-Stokes(RANS)方程為控制方程，采用標(biāo)準(zhǔn)k-ε湍流模型求解，得到其流場特征和氣動特性。重點針對大迎角情況，分別對整機(jī)氣動特性、進(jìn)氣道性能和全動尾翼氣動性能進(jìn)行了分析，并結(jié)合流場特征作出解釋。結(jié)果表明，機(jī)身和發(fā)動機(jī)之間存在氣動/推進(jìn)耦合現(xiàn)象。大迎角下飛行器的氣動參數(shù)表現(xiàn)出非線性特性，升阻比減小，整機(jī)縱向表現(xiàn)為靜不穩(wěn)定，且不穩(wěn)定性隨迎角增大而增大；進(jìn)氣道性能在大迎角下降低，從而導(dǎo)致發(fā)動機(jī)推力下降，不利于發(fā)動機(jī)的正常工作，但卻適當(dāng)降低了整機(jī)的縱向靜不穩(wěn)定度；全動尾翼操縱效率降低從而使得配平難度增大。

吸氣式； 高超聲速； 數(shù)值模擬； 大迎角； 進(jìn)氣道； 尾翼

高超聲速巡航飛行器可用作天地往返運輸系統(tǒng)和高速導(dǎo)彈，具有巨大的軍事價值和潛在的經(jīng)濟(jì)價值。根據(jù)推進(jìn)系統(tǒng)的不同，可將其分為火箭動力高超聲速飛行器(RHV)和吸氣式高超聲速飛行器(AHV)。其中AHV由于使用了高比沖、輕載荷的超燃沖壓發(fā)動機(jī)，利用空氣中的氧氣作為氧化劑，因而具有結(jié)構(gòu)輕、成本低等優(yōu)異性能。近年來各航空航天大國都將AHV作為發(fā)展重點，對其進(jìn)行了廣泛研究[1-2]，其中最具代表性的是美國的Hyper-X計劃[3]。

高超聲速流動的復(fù)雜性給飛行器設(shè)計帶來了極大的困難，氣動問題就是其中的一個重點難題。目前國內(nèi)外已開展了大量針對AHV氣動特性的實驗研究。Holland等[4]使用美國國家航空航天局(NASA)Langley 研究中心風(fēng)洞試驗測量出了一種X-43飛行器縮小模型在進(jìn)氣道閉合時的氣動力參數(shù)，F(xiàn)an等[5]利用風(fēng)洞試驗研究了一種AHV進(jìn)氣道的邊界層轉(zhuǎn)捩特性。由于試驗耗時長且成本高，相比之下工程算法則非常高效，因此也有眾多學(xué)者采用各種工程算法進(jìn)行研究。其中，Skujins等[6]基于激波膨脹波法建立了一種鴨式布局AHV的氣動模型，重點研究了鴨翼對尾翼的氣動干擾現(xiàn)象，并和數(shù)值方法作了對比，其計算迎角為-1°～5°。曾開春等[7-9]使用激波膨脹波法結(jié)合當(dāng)?shù)亓骰钊碚摻⒘艘环NAHV在迎角為±5°以內(nèi)的氣動模型。

雖然工程算法非常高效，卻以犧牲精度為代價，而數(shù)值方法不僅計算精度高，且比試驗方法節(jié)省了大量時間和費用，尤其是近年來隨著計算機(jī)技術(shù)的快速發(fā)展，采用計算流體力學(xué)(CFD)方法進(jìn)行研究存在很大的優(yōu)勢，國內(nèi)外開展了一些關(guān)于高超聲速飛行器的CFD模擬研究[10-11]。目前關(guān)于高超聲速飛行器的CFD研究主要是針對RHV，而針對AHV的則相對較少。其中，Cui等[12]設(shè)計了一種具有雙進(jìn)氣道的AHV機(jī)身前體，指出機(jī)身前體是影響整機(jī)升阻力特性的重要部分，計算迎角為-2°～10°。Gollan和Smart[13]在一種錐形機(jī)身上設(shè)計了“矩形-橢圓形”過渡的進(jìn)氣道，表明進(jìn)氣道與機(jī)身前體之間存在氣動干擾。Taguchi等[14]計算了一種AHV的氣動力，但計算時未將發(fā)動機(jī)包括在內(nèi)，計算迎角為0°～10°，結(jié)果顯示升力系數(shù)與迎角幾乎呈線性關(guān)系。Mirmirani等[15]用機(jī)身和發(fā)動機(jī)的二維模型計算了迎角±5°以內(nèi)的氣動力，用于驗證理論計算的準(zhǔn)確性。

到目前為止，采用數(shù)值方法對AHV進(jìn)行的研究多以機(jī)身或進(jìn)氣道性能為重點，極少數(shù)關(guān)于整機(jī)的氣動分析也是處在較小迎角范圍內(nèi)，尚未見到針對大迎角的整機(jī)氣動特性研究。事實上，高速飛行過程中，在大氣紊流等外部干擾的作用下，飛行器很可能會出現(xiàn)大迎角工作狀態(tài)，由于超燃沖壓發(fā)動機(jī)的正常工作對飛行器姿態(tài)具有苛刻的限制，同時大迎角下為了恢復(fù)正常姿態(tài)還需要操縱面大角度偏轉(zhuǎn)，因此很有必要對包括發(fā)動機(jī)、操縱面在內(nèi)的整機(jī)大迎角氣動問題進(jìn)行研究。本文對一種典型AHV進(jìn)行了數(shù)值仿真模擬，重點分析大迎角氣動特性。

1 吸氣式高超聲速飛行器模型

AHV的特點是采用機(jī)體/發(fā)動機(jī)一體化設(shè)計，為避免波阻增大，推進(jìn)系統(tǒng)必須位于機(jī)身前體產(chǎn)生的激波內(nèi)，因此推進(jìn)系統(tǒng)通常置于非?？亢蟮奈恢?，從而使機(jī)體前部呈尖楔形。飛行器靠前體下表面產(chǎn)生附體激波，提供大部分的升力，并為發(fā)動機(jī)提供高壓空氣。這種機(jī)身/發(fā)動機(jī)一體化的特性，使得飛行器的氣動-推進(jìn)系統(tǒng)之間存在強(qiáng)烈的耦合作用，造成其氣動性能的復(fù)雜性。

根據(jù)典型AHV驗證機(jī)得出的標(biāo)準(zhǔn)模型如圖1所示。本文旨在研究縱向特性，操縱面僅考慮水平全動尾翼。定義坐標(biāo)系為：坐標(biāo)原點取機(jī)身縱向?qū)ΨQ面內(nèi)前緣處，x軸沿機(jī)身軸線指向機(jī)尾為正，z軸在機(jī)身縱向?qū)ΨQ面上垂直于x軸指向上為正，y軸垂直于x軸和z軸，按右手定則給定正方向。圖1(b)中：τ1、τ2和τ3分別為上表面傾角、下表面傾角和后楔角；V∞為來流速度；α為迎角。

圖1 AHV模型Fig.1 AHV model

2 數(shù)值方法與驗證

2.1 數(shù)值方法

由于高超聲速流動涉及激波、邊界層、湍流、化學(xué)反應(yīng)等復(fù)雜的物理化學(xué)現(xiàn)象，要在數(shù)值計算中模擬所有現(xiàn)象難度很大。本文忽略了化學(xué)反應(yīng)、高溫氣體效應(yīng)等因素，以雷諾平均Navier-Stokes(RANS)方程為控制方程，采用標(biāo)準(zhǔn)k-ε湍流模型計算雷諾應(yīng)力項，在近壁區(qū)域采用壁面函數(shù)進(jìn)行修正。

標(biāo)準(zhǔn)k-ε湍流模型方程為

Gk+Gb-ρ ε-YM+Sk

(1)

(2)

式中：ρ為密度；t為時間；ui為時均速度；xi和xj為速度分量；μ為黏度，μt為湍動黏度；Gk為由平均速度梯度引起的湍動能k的產(chǎn)生項；Gb為由浮力引起的湍動能k的產(chǎn)生項；YM表示可壓湍流中脈動擴(kuò)張的貢獻(xiàn)；C1ε、C2ε和C3ε為經(jīng)驗常數(shù)；σk和σε分別為與湍動能k和耗散率ε對應(yīng)的Prandtl數(shù)；Sk和Sε為用戶定義的源項。

對模型表面及流場的結(jié)構(gòu)網(wǎng)格劃分如圖 2所示。由于無橫側(cè)向運動，取半模計算，在氣動影響較為嚴(yán)重的地方對網(wǎng)格進(jìn)行加密，網(wǎng)格總數(shù)為238萬。在全動尾翼與機(jī)身連接處采用滑移網(wǎng)格處理，便于偏轉(zhuǎn)尾翼，同時為下一步的動態(tài)研究作準(zhǔn)備。

圖2 計算網(wǎng)格Fig.2 Computational grid

該飛行器在設(shè)計巡航狀態(tài)下，馬赫數(shù)Ma=8，高度為30 km，壓強(qiáng)、溫度等氣體參數(shù)由標(biāo)準(zhǔn)大氣表查得。由于前體下表面必須產(chǎn)生附體斜激波才能為發(fā)動機(jī)提供高壓空氣，因此低于-6.2°時前體下表面不再產(chǎn)生附體激波，不具備研究價值。同時考慮到迎角大于15°時，飛行器很可能已經(jīng)承受不住氣動加熱而燒毀，也沒有研究意義。所以，本文選取計算迎角在-6°～15°范圍內(nèi)。

2.2 算例驗證

由于數(shù)值計算結(jié)果與湍流模型、計算網(wǎng)格、算法等一系列因素相關(guān)，為了驗證所用數(shù)值方法的可靠性，選取一個典型再入飛行器機(jī)體前部的橢球體模型為驗證模型，將數(shù)值計算與實驗結(jié)果[16]進(jìn)行對比，重點考察大迎角下的吻合程度。計算條件為Ma=8，雷諾數(shù)Re=1.98×107，總溫T0=720 K，總壓P0=8.5 MPa。采用2.1節(jié)中的數(shù)值方法，對比結(jié)果如圖3所示，圖中顯示了迎角α=-5°,10°，25°時雙橢球體對稱面上下兩側(cè)的壓力系數(shù)Cp分布，其中橫坐標(biāo)表示沿軸向無量綱位置，機(jī)頭處為0。

圖3 數(shù)值計算驗證結(jié)果Fig.3 Numerical simulation validation results

從圖3可以看出，在對稱面上下側(cè)，數(shù)值計算結(jié)果與實驗結(jié)果幾乎完全一致，在25°大迎角下吻合程度依然非常高。由此可以認(rèn)為，本文所采用的數(shù)值計算方法完全可以用于研究高超聲速飛行器大迎角下的氣動特性。

3 結(jié)果與分析

3.1 整機(jī)氣動特性

本節(jié)主要研究飛行器整機(jī)氣動特性，圖4給出了整機(jī)及各部分升力系數(shù)CL隨迎角的變化曲線?？梢钥闯鲭S著迎角的變化，整機(jī)升力系數(shù)逐漸增加，在迎角為-6°～5°之間基本呈線性增長，迎角大于5°時，升力系數(shù)表現(xiàn)出非線性增加。機(jī)身是產(chǎn)生升力的主要部分，在迎角大于5°時升力系數(shù)也出現(xiàn)非線性增加。另外，隨著迎角的增大，尾翼產(chǎn)生的升力也逐漸增大，但是數(shù)值較小。而發(fā)動機(jī)部分產(chǎn)生的升力為負(fù)值，且迎角越大負(fù)升力越大，存在氣動/推進(jìn)耦合現(xiàn)象。

圖4 升力系數(shù)曲線Fig.4 Lift coefficient curves

圖5 機(jī)身表面附近流線分布Fig.5 Streamlines distribution near surface of fuselage

為進(jìn)一步理解非線性形成的原因，觀察圖5所示的機(jī)身表面附近流線可以發(fā)現(xiàn)，0°迎角時，機(jī)身上表面流線附著在物面，而15°迎角時，機(jī)身上表面處于背風(fēng)膨脹波區(qū)，氣流分離后形成剪切層，剪切層卷積形成旋渦，從而對表面附近的流動產(chǎn)生誘導(dǎo)作用，與未形成旋渦的“死水區(qū)”相比，減小了表面壓強(qiáng)，造成了非線性升力。

圖6給出了飛行器整機(jī)及各部分阻力系數(shù)CD隨迎角的變化曲線。整機(jī)阻力系數(shù)以0°迎角為對稱呈拋物線趨勢增長。其中產(chǎn)生阻力的主要部分仍然是機(jī)身，尾翼產(chǎn)生的阻力隨迎角增大而增大。發(fā)動機(jī)在迎角大于3°以后產(chǎn)生的阻力變?yōu)樨?fù)值，表現(xiàn)為“推力”。這是因為氣流在發(fā)動機(jī)腔內(nèi)來回折射，對發(fā)動機(jī)內(nèi)下表面產(chǎn)生壓力，壓力沿逆流方向上的分量即表現(xiàn)出的“推力”。

圖6 阻力系數(shù)曲線Fig.6 Drag coefficient curves

圖7 俯仰力矩系數(shù)曲線Fig.7 Pitching moment coefficient curves

圖7所示為整機(jī)及各部分俯仰力矩系數(shù)Cm隨迎角的變化曲線。整機(jī)和機(jī)身俯仰力矩系數(shù)隨迎角增大而增大，表現(xiàn)出靜不穩(wěn)定。這是由于整機(jī)升力主要由機(jī)身前體下表面產(chǎn)生，因此整機(jī)氣動中心位于重心之前，從而導(dǎo)致縱向靜不穩(wěn)定。整機(jī)和機(jī)身俯仰力矩系數(shù)同樣在5°迎角以后出現(xiàn)非線性增長。發(fā)動機(jī)部分的俯仰力矩系數(shù)幾乎為0，且基本保持不變。尾翼部分的力矩系數(shù)隨迎角增大而減小，且幅度較小，這主要是由于尾翼面積相比機(jī)身而言較小，因此操縱效率并不高。

圖8 升阻比和氣動中心變化曲線Fig.8 Lift-to-drag ratio and aerodynamic center change curves

圖9所示為飛行器表面壓力系數(shù)分布和截面周圍的馬赫數(shù)分布，左下角為對稱面上的壓力云圖?？梢钥闯?，飛行器的前體下表面為升力的主要產(chǎn)生面。0°迎角時，機(jī)體周圍產(chǎn)生完整的激波錐，而在15°迎角下，前體下表面附體激波開始脫離，氣流繞過邊緣在機(jī)身上側(cè)產(chǎn)生膨脹波，機(jī)身上表面壓力降低，從而使得升力增加，機(jī)身上側(cè)的膨脹波在邊緣內(nèi)側(cè)被嵌入的橫流激波終止，形成了低壓區(qū)，進(jìn)一步增大了升力。同時從左下角的壓

圖9 表面壓力系數(shù)和截面馬赫數(shù)分布Fig.9 Distribution of surface pressure coefficient and cross-section Mach number

力云圖可以看出，0°迎角時，氣流經(jīng)過前體下表面的壓縮作用直接進(jìn)入發(fā)動機(jī)，而15°迎角時，激波更加靠近前體下表面，并和進(jìn)氣道外部激波系發(fā)生干擾，使得激波向外折射，在引起波阻增量的同時，也會對發(fā)動機(jī)性能產(chǎn)生不利影響。

3.2 進(jìn)氣道性能

超燃沖壓發(fā)動機(jī)對于工作環(huán)境有著嚴(yán)苛的要求，機(jī)身前體下表面作為發(fā)動機(jī)的預(yù)壓縮面，其作用相當(dāng)于進(jìn)氣道。由于推進(jìn)系統(tǒng)即使在試驗中也難以集成到地面模擬模型中，因此在本文采用的數(shù)值模擬中，發(fā)動機(jī)始終處于冷噴狀態(tài)。

圖10給出了進(jìn)氣道增壓比p/p∞和總壓恢復(fù)系數(shù)σ隨迎角的變化趨勢?？梢钥闯觯窃酱?，增壓比越大，這是由于迎角越大，形成的激波越強(qiáng)，激波過后的壓強(qiáng)也就越大。-6°迎角時增壓比只有1.1，進(jìn)氣質(zhì)量遠(yuǎn)不能滿足發(fā)動機(jī)的要求?？倝夯謴?fù)系數(shù)在-3°迎角時達(dá)到最大值0.994，大迎角下，總壓恢復(fù)系數(shù)急劇降低，到15°迎角時只有0.24，總壓損失嚴(yán)重，會大大減小發(fā)動機(jī)產(chǎn)生的推力。

圖10 增壓比和總壓恢復(fù)系數(shù)變化曲線Fig.10 Pressure ratio and total pressure recovery coefficient change curves

圖11 捕獲質(zhì)量流量變化曲線Fig.11 Mass flow capture change curve

3.3 發(fā)動機(jī)性能

根據(jù)文獻(xiàn)[17]中的超燃沖壓發(fā)動機(jī)模型，可將發(fā)動機(jī)簡化為由擴(kuò)壓段、燃燒室和內(nèi)噴管組成，如圖12(a)所示。假設(shè)發(fā)動機(jī)內(nèi)部流動為一維流動，由數(shù)值模擬結(jié)果得到入流參數(shù)，在考慮燃油的基礎(chǔ)上計算出內(nèi)流參數(shù)，再根據(jù)發(fā)動機(jī)入口及出口處參數(shù)，結(jié)合沖量定理，計算出發(fā)動機(jī)推力。具體計算方法參見文獻(xiàn)[17]，計算公式為

(3)

圖12 發(fā)動機(jī)截面形狀和發(fā)動機(jī)推力曲線Fig.12 Cross-section and thrust curve of engine

從圖12可以看出，發(fā)動機(jī)所產(chǎn)生的推力在迎角大于5°時呈急劇下降趨勢，此時對應(yīng)的總壓恢復(fù)系數(shù)和質(zhì)量流量也均減小。充分說明大迎角下雖然氣流經(jīng)過更強(qiáng)烈的減速增壓，但也會帶來發(fā)動機(jī)進(jìn)氣質(zhì)量流量下降和總壓恢復(fù)系數(shù)降低的不利影響，從而直接導(dǎo)致發(fā)動機(jī)推力下降，并且可能出現(xiàn)不起動現(xiàn)象，總體來說不利于發(fā)動機(jī)的正常工作。

假設(shè)大迎角下發(fā)動機(jī)仍處于工作狀態(tài)，由于發(fā)動機(jī)推力沿發(fā)動機(jī)軸線指向機(jī)頭，因此產(chǎn)生抬頭力矩。圖13所示為考慮發(fā)動機(jī)推力前后整機(jī)俯仰力矩系數(shù)Cm曲線的對比?？紤]發(fā)動機(jī)推力產(chǎn)生的力矩之后，在迎角小于5°時，整機(jī)俯仰力矩系數(shù)曲線斜率Cmα有所增大，然而迎角大于5°時，由于發(fā)動機(jī)推力急劇減小使得Cmα減小，從而使整機(jī)不穩(wěn)定程度減小。

圖13 俯仰力矩系數(shù)曲線對比Fig.13 Comparison of pitching moment coefficient curves

3.4 全動尾翼氣動性能

考慮到如果尾翼偏角過大，可能發(fā)生熱燒蝕破壞，且產(chǎn)生的配平阻力會影響鉸鏈力矩(對作動器的性能要求)，本文假設(shè)尾翼偏角δt的范圍為±20°，上偏為正。圖14顯示迎角分別為0°和15°時，尾翼的偏轉(zhuǎn)對升力系數(shù)的影響?？梢钥闯觯ο禂?shù)與尾翼偏轉(zhuǎn)角基本呈線性增加趨勢，整機(jī)升力系數(shù)的增量絕大部分是由尾翼偏轉(zhuǎn)引起的。但尾翼偏轉(zhuǎn)同時也影響了機(jī)身的升力系數(shù)，造成機(jī)身部分的小幅變化，而對發(fā)動機(jī)部分基本沒有影響。這是因為尾翼與機(jī)身之間存在氣流干擾，而與發(fā)動機(jī)相距較遠(yuǎn)，不會造成干擾。另外，15°迎角時升力系數(shù)曲線的斜率與0°迎角相比較小。

圖14 尾翼偏轉(zhuǎn)對升力系數(shù)的影響Fig.14 Influences of tail deflection on lift coefficient

圖15顯示迎角分別為0°和15°時，尾翼偏轉(zhuǎn)對俯仰力矩系數(shù)的影響，與對升力系數(shù)的影響相似。另外可以看出，尾翼正偏時對機(jī)身俯仰力矩系數(shù)的影響要大于負(fù)偏時。這是因為，尾翼正偏時，前緣超出機(jī)身上表面，產(chǎn)生的激波對機(jī)身后段上表面流動造成影響，提高了機(jī)身后段上表面的壓強(qiáng)，從而增加了機(jī)身的抬頭力矩；而負(fù)偏時，尾緣超出機(jī)身，位置較前者靠后很多，產(chǎn)生的激波對機(jī)身后段的影響沒有前者大。同樣的，可以看出和0°迎角相比，15°迎角時尾翼的操縱效率有所降低。

圖15 尾翼偏轉(zhuǎn)對俯仰力矩系數(shù)的影響Fig.15 Influences of tail deflection on pitching moment coefficient

圖 16所示為尾翼表面及附近壓力系數(shù)分布圖，由于尾翼偏轉(zhuǎn)時，上表面處于背風(fēng)區(qū)，表面上的壓力幾乎可以忽略不計，因此圖中顯示的是尾翼下表面上的壓力系數(shù)分布?？梢钥吹剑?°迎角時，尾翼前緣是升力產(chǎn)生的主要位置，來流經(jīng)過機(jī)身之前受到尾翼干擾，通過尾翼與機(jī)身連接處改變了機(jī)身上表面的壓力分布，使得靠近尾翼部分的壓力增大；而15°迎角時，尾翼中后緣是升力產(chǎn)生的主要位置，處在與機(jī)身連接處之后，對機(jī)身的影響不如前者大。因此，15°時升力系數(shù)和俯仰力矩系數(shù)隨尾翼偏角的變化斜率絕對值均小于0°迎角時。

圖16 尾翼表面及附近壓力系數(shù)分布Fig.16 Distribution of surface pressure coefficient on and near tail

因此，迎角為0°時，尾翼偏轉(zhuǎn)8°左右就可以配平，而15°迎角下，尾翼偏轉(zhuǎn)20°依然遠(yuǎn)遠(yuǎn)無法達(dá)到配平要求，配平難度明顯增大。雖然在模擬中并沒有考慮發(fā)動機(jī)的工作，但是依據(jù)前面的推斷，發(fā)動機(jī)此時的工作狀態(tài)不佳，如果發(fā)生不起動現(xiàn)象，那么此時僅靠偏轉(zhuǎn)尾翼使飛行器恢復(fù)正常飛行姿態(tài)具有相當(dāng)大的難度，必須采取其他措施。

4 結(jié) 論

1) 大迎角下，飛行器的氣動參數(shù)表現(xiàn)出非線性；迎角超過8°之后升阻比開始減??；整機(jī)縱向表現(xiàn)為靜不穩(wěn)定，且不穩(wěn)定性隨迎角增大而增大。

2) 隨著迎角增大，氣流經(jīng)過進(jìn)氣道進(jìn)入發(fā)動機(jī)時的增壓比增大，進(jìn)氣道壓縮性增強(qiáng)，但是總壓恢復(fù)系數(shù)降低，15°迎角時只有0.24，總壓損失嚴(yán)重，且進(jìn)氣道捕獲質(zhì)量流量在5°以后急劇降低，從而直接減小了發(fā)動機(jī)產(chǎn)生的推力，甚至可能使發(fā)動機(jī)出現(xiàn)不起動現(xiàn)象，不利于發(fā)動機(jī)正常工作，但卻降低了整機(jī)縱向靜不穩(wěn)定的程度。

3) 大迎角下，全動尾翼的操縱效率會略有降低，配平難度增大，考慮到大迎角下發(fā)動機(jī)的工作狀態(tài)，此時僅靠偏轉(zhuǎn)尾翼使飛行器恢復(fù)正常飛行姿態(tài)具有相當(dāng)大的難度，必須采取其他措施。

[1] Yentsch R J, Gaitonde D V, Kimmel R. Performance of turbulence modeling in simulation of the HIFiRE-1 flight test[J]. Journal of Spacecraft and Rockets, 2014, 51(1): 117-127.

[2] Murphy K J, Nowak R J, Thompson R A, et al. X-33 hypersonic aerodynamic characteristics[J].Journal of Spacecraft and Rockets, 2001, 38(5): 670-683.

[3] Engelund W C, Holland S D, Cockrell C E, et al. Aerodynamic database development for the Hyper-X airframe-integrated scramjet propulsion experiments[J]. Journal of Spacecraft and Rockets, 2001, 38(6): 803-810.

[4] Holland S D, Woods W C, Engelund W C. Hyper-X research vehicle experimental aerodynamics test program overview[J]. Journal of Spacecraft and Rockets, 2001, 38(6): 828-835.

[5] Fan X Q, Li H, Jia D, et al. Experimental investigation on forced boundary-layer transition of axisymmetric inlet and numerical verification[J]. Journal of Aircraft, 2006, 43(5): 1544-1549.

[6] Skujins T, Cesnik C E S, Oppenheimer M W, et al. Canard-elevon interactions on a hypersonic vehicle[J]. Journal of Spacecraft and Rockets, 2010, 47(1): 90-100.

[7] Zeng K C, Xiang J W. Uncertainty analysis of flight dynamic characteristics for hypersonic vehicles[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(4): 798-808(in Chinese). 曾開春, 向錦武. 高超聲速飛行器飛行動力學(xué)特性不確定分析[J]. 航空學(xué)報, 2013, 34(4): 798-808.

[8] Zeng K C, Xiang J W, Li D C. Aeroservoelastic modeling and analysis of a canard-configured air-breathing hypersonic vehicles[J]. Chinese Journal of Aeronautics, 2013, 26(4): 831-840.

[9] Xiang J W, Zeng K C, Nie L. Elastic effects on flight mechanics of waverider[J]. Journal of Beijing University of Aeronautics and Astronautics, 2012, 38(10): 1306-1310 (in Chinese). 向錦武, 曾開春, 聶璐. 考慮彈性影響的乘波體飛行動力學(xué)特性[J]. 北京航空航天大學(xué)學(xué)報, 2012, 38(10): 1306-1310.

[10] Gupta K K, Voelker L S. Aeroelastic simulation of hy-personic flight vehicles[J]. AIAA Journal, 2012, 50(3): 717-723.

[11] Chen Q, Si F F, Chen J Q, et al. Study of protuberance in supersonic flow with RANS/LES method[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(7): 1531-1537 (in Chinese). 陳琦， 司芳芳， 陳堅強(qiáng)， 等. RANS/LES 在超聲速突起物繞流中的應(yīng)用研究[J]. 航空學(xué)報, 2013, 34(7): 1531-1537.

[12] Cui K, Hu S, Li G L, et al. Conceptual design and aerodynamic evaluation of hypersonic airplane with double flanking air inlets[J]. Science China Technological Sciences, 2013, 56(8): 1980-1988.

[13] Gollan R J, Smart M K. Design of modular shape-transition inlets for a conical hypersonic vehicle[J]. Journal of Propulsion and Power, 2013, 29(4): 832-838.

[14] Taguchi H, Kobayashi H, Kojima T, et al. Research on hypersonic aircraft using pre-cooled turbojet engines[J]. Acta Astronautica, 2012, 73(1): 164-172.

[15] Mirmirani M, Wu C, Clark A, et al. Modeling for con-trol of a generic airbreathing hypersonic vehicle[C]∥AIAA Guidance, Navigation, and Control Conference and Exhibit. San Francisco: AIAA, 2005: 6256.

[16] Li S X. Typical hypersonic flow characteristics[M]. Beijing: National Defence Industry Press, 2007: 18-39(in Chinese). 李素循. 典型外形高超聲速流動特征[M]. 北京: 國防工業(yè)出版社, 2007: 18-39.

[17] Bolender M A, Doman D B. Nonlinear longitudinal dynamical model of an air-breathing hypersonic vehicle[J]. Journal of Spacecraft and Rockets, 2007, 44(2): 374-387.

Tel: 010-82338786

E-mail: buaaluowenli@163.com

向錦武 男， 博士， 教授， 博士生導(dǎo)師。主要研究方向：飛行器設(shè)計、結(jié)構(gòu)動力學(xué)和氣動彈性等。

Tel: 010-82338786

E-mail: xiangjw@buaa.edu.cn

*Corresponding author. Tel.： 010-82338786 E-mail: xiangjw@buaa.edu.cn

Aerodynamic characteristics analysis of air-breathing hypersonic vehicles at high angle of attack

LUO Wenli, LI Daochun, XIANG Jinwu*

SchoolofAeronauticScienceandEngineering,BeihangUniversity,Beijing100191,China

During the flight of air-breathing hypersonic vehicles, the atmospheric turbulence or other external disturbances may result in high angle of attack conditions. In order to study the aerodynamic performances influenced by high angle of attack, numerical simulations are performed on the flow field of a typical air-breathing hypersonic vehicle. The flow field and aerodynamic characteristics are obtained by solving the Reynolds-averaged Navier-Stokes (RANS) equations using standardk-εturbulence models. The performances of vehicle, inlet and the all-moving tails are analyzed and explained based on the flow field characteristics especially at high angle of attack. The investigation indicates that there is coupling between aerodynamics and propulsive system. At high angle of attack, aerodynamic parameters exhibit nonlinear characteristics, meanwhile, the lift-to-drag ratio begins to reduce and the longitudinal instability enhances. The inlet has poor performance at high angle of attack, resulting in the decrease of thrust provided by the engine, which is not conducive to the normal operation of the engine. However, on the other hand, the decrease of thrust will reduce the longitudinal instability of the vehicle. In addition, the control efficiency of the all-moving tails reduces and consequently makes trimming difficult.

air-breathing; hypersonic; numerical simulation; high angle of attack; inlet; tail

2014-06-04； Revised： 2014-09-10； Accepted： 2014-09-22； Published online： 2014-10-19 13:14

s： National Natural Science Foundation of China (91216102, 11402014)

2014-06-04; 退修日期： 2014-09-10; 錄用日期： 2014-09-22; 網(wǎng)絡(luò)出版時間： 2014-10-09 13：14

www.cnki.net/kcms/detail/10.7527/S1000-6893.2014.0266.html

國家自然科學(xué)基金 (91216102， 11402014)

Luo W L, Li D C, Xiang J W. Aerodynamic characteristics analysis of air-breathing hypersonic vehicles at high angle of attack[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(1): 223-231. 羅文莉， 李道春， 向錦武. 吸氣式高超聲速飛行器大迎角氣動特性分析[J]. 航空學(xué)報, 2015, 36(1): 223-231.

http://hkxb.buaa.edu.cn hkxb@buaa.edu.cn

10.7527/S1000-6893.2014.0266

V211.4

A

1000-6893(2015)01-0223-09

羅文莉 女， 碩士研究生。主要研究方向：飛行器氣動彈性。

*通訊作者.Tel.： 010-82338786 E-mail: xiangjw@buaa.edu.cn

URL： www.cnki.net/kcms/detail/10.7527/S1000-6893.2014.0266.html