段焰輝， 范召林， 吳文華

中國空氣動力研究與發(fā)展中心 計算空氣動力研究所， 綿陽 621000

定后掠角密切錐乘波體的生成和設(shè)計方法

段焰輝， 范召林， 吳文華*

中國空氣動力研究與發(fā)展中心 計算空氣動力研究所， 綿陽 621000

對定后掠角密切錐乘波體(OCWRCAS)的生成方法和考慮黏性的設(shè)計方法進(jìn)行了研究。定后掠角乘波體的前緣具有特定的后掠角，能夠在上表面產(chǎn)生穩(wěn)定分離渦從而改善乘波體的氣動性能。本文首先在傳統(tǒng)密切錐乘波體生成方法的基礎(chǔ)上給出了定后掠角密切錐乘波體的生成方法；從前緣后掠的幾何特征中提取了后掠角、激波角和前緣曲線程度等設(shè)計變量，并研究了設(shè)計變量的取值范圍；以遍歷設(shè)計空間的思路對兩類定后掠角密切錐乘波體進(jìn)行了設(shè)計分析，研究了升阻比、體積效率隨設(shè)計變量的變化規(guī)律，然后在設(shè)計空間內(nèi)進(jìn)行了多目標(biāo)尋優(yōu)；最后使用計算流體力學(xué)方法對定后掠角乘波體的乘波特性和渦升力特性進(jìn)行了驗證。結(jié)果表明，由本文生成方法得到的定后掠角密切錐乘波體具有明顯的乘波特性并且能夠在較高的升阻比時保證一定的體積效率；定后掠角前緣能夠在一定的迎角下在上表面產(chǎn)生穩(wěn)定的分離渦，產(chǎn)生渦升力。

乘波體； 密切錐； 黏性； 后掠角； 渦升力

乘波體經(jīng)過幾十年的發(fā)展，從早期的單一構(gòu)型逐漸發(fā)展為具有更多特點的復(fù)雜構(gòu)型，處理手段也愈加靈活。Nonweiler[1]在1959年提出了由二維流場構(gòu)造高超聲速飛行器的理論，并據(jù)此生成了“Λ”型乘波體，隨后Jones等[2]提出了基于軸對稱流場的乘波體生成理論，國內(nèi)也對錐導(dǎo)乘波體進(jìn)行了詳細(xì)研究[3-4]，將軸對稱流場擴(kuò)展至更一般的三維流場，就出現(xiàn)了源于帶迎角錐和橢圓錐體流場等非軸對稱流場的乘波體[5-6]。發(fā)動機(jī)一體化設(shè)計需要下表面有較為均勻的流場，源于楔型-錐型混合流場的乘波體[7]以及由反設(shè)計思路得到的密切錐或密切流場乘波體[8-10]具有“平坦”的下表面，符合上述設(shè)計條件。在諸多生成方法中，密切錐(或密切流場)乘波體生成方法靈活，可以生成具有更多特性的乘波體，國外已將這種方法發(fā)展為成熟軟件，能夠生產(chǎn)具有復(fù)雜平面形狀的乘波體[11]，國內(nèi)由賀旭照等[12-13]提出來的密切曲面錐乘波體也是基于這種反設(shè)計思路得到的。

根據(jù)乘波體的原理，其氣動性能主要來源于下表面，如果能有效利用乘波體的上表面將有可能改善乘波體的整體性能。傳統(tǒng)做法是將上表面整個或部分區(qū)域設(shè)計為膨脹面[14]，若只有膨脹而無壓縮則必然降低乘波體的體積，這對體積效率本就很低的乘波體是非常不利的(若采用先壓縮后膨脹表面，雖然有可能提高體積，但壓縮區(qū)域必然增大乘波體阻力)。定后掠角乘波體(Waverider with Constant Angle of Sweepback， WRCAS)，能夠通過具有固定后掠角的前緣在上表面產(chǎn)生穩(wěn)定分離渦從而提供更高的升力，而無需將上表面設(shè)計為犧牲體積的膨脹面；更重要的是這種乘波體的前緣后掠角能夠作為設(shè)計參數(shù)在設(shè)計階段予以控制。

定后掠角乘波體由來已久，但在早期只是其他型乘波體的“附屬品”。由Nonweiler提出來的“Λ”型乘波體就是一種定后掠角乘波體，但是由于其升阻比不高，體積利用率低等問題，難以實用。Jones在對錐導(dǎo)乘波體進(jìn)行介紹時，提到了一種流動捕獲管(Flow Capture Tube， FCT)過圓錐頂點的乘波體，這種乘波體也屬于定后掠角乘波體，也因為實用價值不高未受重視。真正以定后掠角前緣為研究重點的是洛克希德·馬丁公司的Rodi[15-17]，其在文獻(xiàn)中對定后掠角密切錐乘波體(Osculating Cone Waverider with Constant Angle of Sweepback，OCWRCAS)和定后掠角密切流場乘波體(Osculating Flow Waverider with Constant Angle of Sweepback， OFWRCAS)進(jìn)行了研究，其研究重點在于此類乘波體的幾何特征和渦升力特性，卻并未說明怎樣生成定后掠角前緣對應(yīng)的乘波體下表面，而這也正是這種乘波體的重點之一，因為只有合理的下表面才能確保乘波體的乘波特性。

本文首先根據(jù)OCWRCAS的幾何特征及密切錐乘波體的生成方法，給出了定后掠角前緣對應(yīng)的乘波體下表面生成方法；之后提取了兩種典型定后掠角密切錐乘波體的設(shè)計變量，以遍歷設(shè)計空間的思路進(jìn)行了多目標(biāo)優(yōu)化設(shè)計；最后使用CFD方法對此類乘波體的乘波特性和渦升力特性進(jìn)行了驗證。

1 OCWRCAS的生成及氣動性能計算方法

1.1 OCWRCAS的生成方法

OCWRCAS是密切錐乘波體的一種，生成方法與其類似，因此本文從密切錐乘波體的生成方法入手，引出OCWRCAS的生成方法。密切錐乘波體的生成方法是一種反設(shè)計方法，這種方法由給定的激波形狀來確定乘波外形。所以在設(shè)計時除了要給定流動捕獲管FCT，還需要給出激波形狀，即進(jìn)氣捕獲曲線(Inlet Capture Curve， ICC)，該曲線與進(jìn)氣道進(jìn)口截面形狀有關(guān)。這兩條曲線必須一階導(dǎo)數(shù)連續(xù)，ICC還要保證二階導(dǎo)數(shù)非負(fù)。對于任意的ICC，其對應(yīng)的流場一般都是三維的，直接計算這種流場計算量較大不利于快速生成乘波體。Sobieczky等[8]提出的相切錐理論大大簡化了由激波形狀反向計算流場的過程,首先定義過ICC上任意一點并且垂直于這個點切線方向的平面為密切平面，然后使用一系列密切平面內(nèi)的錐型流近似這種三維流動，這樣就大大提高了計算效率。密切錐的頂點可由密切平面內(nèi)對應(yīng)點的曲率半徑和激波角確定，當(dāng)曲率半徑無窮大時，當(dāng)做二維流場處理。密切錐確定之后，可根據(jù)錐型流理論求得解析解，再由FCT和激波形狀確定該密切面內(nèi)的前緣點，并以該點作為起始點在密切平面內(nèi)做流線追蹤，得到一條流線。所有密切平面上的流線組成的流面即為密切錐乘波體的下表面。生成密切錐乘波體的兩個充分條件是：① 每個密切面內(nèi)流場的激波必須與激波形狀(ICC)一致；② 相鄰密切面內(nèi)的橫向流動必須足夠小。要滿足這兩個條件，一般都限定每個密切面內(nèi)激波角相等。

在滿足上述條件的前提下，要得到具有一定后掠角的直線前緣，必須保證前緣對應(yīng)的激波面和流動捕獲管都是平面。鑒于此，定后掠角密切錐乘波體與一般密切錐乘波的不同之處在于：① ICC上使用一條直線段來生成帶有特定后掠角的前緣；② 對于ICC上曲率半徑無限大的情況要進(jìn)行判斷，如果是對應(yīng)發(fā)動機(jī)進(jìn)口的直線段，使用二維流場計算，如果對應(yīng)后掠前緣則使用錐型流計算；③ 直線ICC的每個密切面上的密切錐頂點需要滿足一定的分布規(guī)律以使激波面為平面，本文取與直線激波平行的直線；④ FCT上與直線ICC對應(yīng)的部分應(yīng)為直線，為了使直線前緣處于水平面內(nèi)，取水平直線段。

圖1 兩種乘波體的密切錐示意圖Fig.1 Schematic diagram of osculating cone of two types of waverider

圖1給出了密切錐乘波體(OCWR)和OCWRCAS的密切錐示意圖，實線分別表示FCT和ICC，虛線表示密切平面，“△”表示ICC上密切平面內(nèi)的點，“□”表示對應(yīng)密切錐的頂點。OCWR的ICC比較任意，密切錐頂點由曲線當(dāng)?shù)氐那拾霃酱_定，由圖1(a)可以看出，其密切錐頂點分布與FCT無關(guān)，靠近對稱面附近的水平直線曲率半徑為無限大，計算時使用二維方法。OCWRCAS的FCT和ICC如圖1(b)所示，OCWRCAS的ICC由圓弧(紅色部分)和直線段(藍(lán)色部分)組成，圓弧對應(yīng)曲線型頭部，直線對應(yīng)可控后掠前緣，其對應(yīng)密切錐頂點沿與其平行的直線分布。

1.2 圓錐流場及黏性力計算方法

由定后掠角密切錐乘波體的生成方法可知，生成該型乘波體需要錐型流場。本文使用Taylor-Maccoll方程[18]計算錐型流場，該方程的無量綱形式可以表示為

(1)

(2)

(3)

定義無量綱速度為

(4)

式中：Vmax為氣體膨脹至絕對零度時的理想最大速度。該方程為常微分方程，可采用四階Runge-Kutta求解。

乘波體的黏性阻力與波阻處于同一量級，因此升阻比分析必須考慮黏性阻力。本文采用文獻(xiàn)[19]中的參考溫度方法，該方法比邊界層積分方法的精度稍低但計算速度卻要快得多，因此具有很大的實用價值。本文假設(shè)乘波體表面全部為湍流，并對不可壓平板的黏性力計算方法進(jìn)行了可壓縮修正，提高了計算精度。首先由式(5)計算參考溫度：

(5)

式中：T*為參考溫度；Tδ為邊界層外緣溫度；Maδ為邊界層外緣馬赫數(shù)；TW為壁面溫度，可在計算時指定或者由氣動熱計算方法計算。參考雷諾數(shù)為

(6)

式中：Vδ為邊界層外緣速度；μ*=fT*為參考黏性系數(shù)，是參考溫度的函數(shù)，可由撒特蘭公式求得；ρ*=gp*,T*為參考密度，可根據(jù)氣體狀態(tài)方程求得。附面層內(nèi)壓力梯度較小，所以取p*=pδ，pδ為邊界層外緣壓力。不可壓平板湍流黏性力系數(shù)為

(7)

可壓縮黏性力系數(shù)的修正因子為

Fc=ρδ/ρ*

(8)

式中：ρδ為邊界層外緣密度。

最終求得當(dāng)?shù)仞ば粤ο禂?shù)為

(9)

2 OCWRCAS的設(shè)計方法

OCWRCAS的生成方法屬于密切錐乘波體的一種，為便于研究，本文將FCT設(shè)定為水平直線段，其外形就主要由ICC控制，因此生成此類乘波體的關(guān)鍵就是ICC曲線。一般密切錐乘波體的ICC只要滿足一階導(dǎo)數(shù)連續(xù)，二階導(dǎo)數(shù)非負(fù)的條件即可，而OCWRCAS則必須有一條位于遠(yuǎn)離對稱面一側(cè)的直線段來生成具有一定后掠角的前緣。

本文對其中兩種典型外形進(jìn)行研究，分別用OCWRCAS I和OCWRCAS II表示。OCWRCAS I頭部為尖頂點，即整個前緣為后掠角可控的直線段前緣；OCWRCAS II頭部為曲線，除頭部外的部分前緣為后掠角可控的直線段前緣。而且ICC中的曲線都選擇圓弧，這樣減少了設(shè)計變量的個數(shù)，便于將研究重點放在后掠角部分。通過簡化，這兩種乘波體的設(shè)計變量都很少，因此文中采用遍歷的思路進(jìn)行優(yōu)化設(shè)計。

乘波體的設(shè)計目標(biāo)一般都為升阻比最大，同時要求具有一定的體積效率。體積效率的定義方法有多種[20]，本文使用

τ=V2/3/Sp

(10)

式中：Sp為乘波體的平面面積。本文乘波體的設(shè)計狀態(tài)均為：30 km高度，來流馬赫數(shù)Ma∞=6，乘波體長度取20 m。

2.1 OCWRCAS I的設(shè)計

OCWRCAS I的ICC由兩部分組成：圓弧段和直線段，并且圓弧的圓心O都位于FCT上。這類乘波體的俯視圖和后視圖(側(cè)視圖基本相同，此處不予畫出)如圖2所示，整個前緣為一條直線段VE，具有相同的后掠角，頭部為一個尖點。本文從其幾何特征入手，分析影響乘波體性能的設(shè)計變量，并提出了具有針對性的設(shè)計方法。

圖2 OCWRCAS I的幾何特征Fig.2 Geometry characteristic of OCWRCAS I

2.1.1 幾何特征分析

如圖2所示，OCWRCAS I的ICC由兩部分組成：圓心在O點的圓弧RA和直線段AE，且兩者相切于A點。圓弧部分對應(yīng)頭部的頂點，直線部分對應(yīng)具有某一后掠角的前緣。這樣的曲線由兩個變量控制，一個是圓弧的半徑r，一個是直線與FCT的夾角γ。定義機(jī)身長度為l，機(jī)身寬度為s，激波角為β，后掠角為λ。

r可由機(jī)身長度和激波角確定，即

r=OA=OR=l·tanβ

γ可根據(jù)相切關(guān)系求得，即

sinγ=OA/OE

式中：OE為機(jī)身寬度，可由機(jī)身長度和后掠角確定，即

s=OE=l·tanλ

最終可得[16]

sinγ=tanλ·tanβ

(11)

與一般的密切錐乘波體相同，當(dāng)給定來流馬赫數(shù)和設(shè)計激波角后，乘波體下表面只和ICC有關(guān)。根據(jù)上述分析，對于OCWRCAS I，控制ICC的變量為l、β和λ。所以影響OCWRCAS I性能的變量有4個：Ma∞、l、β和λ。Ma∞和l在設(shè)計前給定，所以O(shè)CWRCAS I最終的設(shè)計變量為激波角和后掠角，即β和λ。接下來分析β和λ的變化范圍。

當(dāng)給定后掠角后，由式(10)可確定激波角上限為

tanλ·tanβ<1

得

β<90°-λ

對于來流馬赫數(shù)Ma∞，馬赫波可以看做是該馬赫數(shù)下最弱的激波，因此可認(rèn)為馬赫角是最小的激波角，所以激波角的范圍為

arcsin1/Ma∞<β<90°-λ

(12)

同時可得后掠角的上限為

λ<90°-arcsin(1/Ma∞)

(13)

當(dāng)Ma∞=6時，arcsin1/Ma∞=9.59°，此時λ<80.41°，取λ=70° 的乘波體為例，激波角的變化范圍為：9.59°<β<20°；若后掠角λ=65°，則激波角的變化范圍為：9.59°<β<35°。高超聲速飛行器的后掠角不能太小，因此本文取后掠角的下限為50°。

通過上述分析，本文分別在后掠角為75°、70°、65°、60°、55°、50° 的情況下，對不同激波角下的乘波體進(jìn)行了分析，具體分配如表1所示，每個后掠角下激波角的變化間隔都為0.5°，共對210個乘波體進(jìn)行了分析。

表1 OCWRCAS I設(shè)計變量分布Table 1 Distribution of design variables of OCWRCAS I (°)

2.1.2 設(shè)計結(jié)果分析

圖3 不同后掠角下升阻比和體積效率隨激波角變化曲線Fig.3 Lift to drag ratio and volumetric efficiency vs shock angle at different angles of sweepback

不同后掠角下氣動性能和體積效率隨激波角變化如圖3所示。圖3(a)給出了不同后掠角時升阻比隨激波角的變化曲線，可以看出，不同后掠角時升阻比CL/CD隨激波角的變化趨勢類似，都呈現(xiàn)先增大后減小的趨勢。升阻比的最大值都出現(xiàn)在11° 激波角時，而且最大值隨著后掠角的增大而減小，但總體變化幅度不大。圖3(b)給出了不同后掠角時體積效率τ隨激波角的變化曲線，不同后掠角時體積效率都隨激波角的增大而增大，若給定激波角，體積效率隨后掠角的增大而增大。

圖4給出了OCWRCAS I隨后掠角和激波角的變化規(guī)律，但并不利于從中選擇“最優(yōu)結(jié)果”。本文對圖3中的結(jié)果按升阻比隨體積效率的變化輸出，如圖4所示，這樣就產(chǎn)生了類似Pareto前沿線的多目標(biāo)優(yōu)化結(jié)果，該前沿線上的結(jié)果包含了升阻比最大、體積效率最大以及升阻比和體積效率適度的所有乘波體。由圖可知，后掠角為75°，激波角大于10.5°，所有乘波體都處于前沿線上，并且屬于升阻比較高部分。根據(jù)文獻(xiàn)[20]，一般高超聲速飛行器的體積效率不應(yīng)低于0.2，本文以此為體積效率的下限，選擇升阻比最大的乘波體，圖中箭頭標(biāo)出了“最優(yōu)”的乘波體外形(Selected)。表2給出了升阻比最大、體積效率最大和“最優(yōu)”乘波體的性能參數(shù)，此3種外形的三視圖如圖5所示?！白顑?yōu)”乘波體升阻比為6.58，體積效率為0.222 1，后掠角為75°，激波角為 12.0°；升阻比最大和體積效率最大的乘波體都為50° 后掠角乘波體，最大升阻比為7.53，最大體積效率為1.14。

圖4 OCWRCAS I升阻比隨體積效率的變化規(guī)律 Fig.4 Lift to drag ratio vs volumetric efficiency for OCWRCAS I

表2 OCWRCAS I優(yōu)化設(shè)計結(jié)果Table 2 Optimal design results of OCWRCAS I

TypeCL/CDV2/3/Spλ/(°)β/(°)MaxCL/CD7．530．11355011．0MaxV2/3/Sp1．140．58615039．0Selected6．580．22217512．0

圖5 OCWRCAS I 3種外形的三視圖Fig.5 Three views of three types of OCWRCAS I

2.2 OCWRCAS II的設(shè)計

OCWRCAS II是指頭部為光滑曲線的定后掠角密切錐乘波體。將圖2中圓心O的位置上移，其他方法不變，便可得到生成OCWRCAS II的ICC，如圖6所示。

圖6 OCWRCAS II的幾何特征Fig.6 Geometry characteristic of OCWRCAS II

2.2.1 幾何特征分析

如圖6所示，ICC的直線段AE對應(yīng)乘波體的后掠前緣ME，這部分與OCWRCAS I的前緣相同，所以式(10)同樣適用，不同之處在于OCWRCAS I的O點與O1是重合的，而OCWRCAS II的O點與O1點存在一定距離，定義為Δr，有限的Δr確保了乘波體靠近頭部的前緣為曲線形狀，該值越大前緣曲線部分所占比例就越大。由Δr與其他已知變量可得圓弧半徑r，聯(lián)立：

最終得

(14)

保持機(jī)身長度和機(jī)身厚度不變，Δr越大，前緣曲線部分所占的寬度就越大，機(jī)身寬度也越大，這在一定程度上增大了乘波體的體積。但是Δr不能無限制地增大，圖7給出了Δr的變化范圍示意圖，保持激波角和后掠角不變，只增大Δr。圖7中Δr2對應(yīng)的ICC為RA2E2曲線，相比較于Δr1對應(yīng)的RA1E1曲線機(jī)身寬度增加，可以預(yù)見體積也有明顯增大。當(dāng)達(dá)到Δrmax時，ICC為RE3整個一條圓弧，相應(yīng)乘波體為錐導(dǎo)乘波體的一部分，但是其前緣離圓錐軸線非常遠(yuǎn)，體積效率很低。Δr作為控制變量，便于說明問題，但是不宜作為設(shè)計變量，因為不同激波角的Δr變化范圍不同，難以進(jìn)行統(tǒng)一量化。本文使用前緣曲線程度變量Δw/s來替換Δr，其中Δw為圓弧對應(yīng)FCT上的寬度，很明顯該變量表征了前緣曲線部分所占比例的大小，Δw/s與Δr和其他已知量之間的關(guān)系為

(15)

通過幾何特征分析，OCWRCAS II相對于OCWRCAS I多了一個設(shè)計變量。對其設(shè)計時，后掠角和激波角仍然沿用OCWRCAS I的設(shè)定，對于Δw/s，為0時生成的乘波體為OCWRCAS I，為1時生成的乘波體為錐導(dǎo)乘波體，因此本文選定其變化范圍為0.04～0.96，變化間隔為0.04，如表3所示，共對5 040個乘波體進(jìn)行了分析。

圖7 Δr變化范圍示意圖Fig.7 Schematic diagram of variation scope of Δr

表3 OCWRCAS II設(shè)計變量分布Table 3 Distribution of design variables of OCWRCAS II

λ/(°)βmin/(°)βmax/(°)Δβ/(°)(Δw/s)min(Δw/s)maxΔ(Δw/s)5010．039．50．50．040．960．045510．034．50．50．040．960．046010．029．50．50．040．960．046510．024．50．50．040．960．047010．019．50．50．040．960．047510．014．50．50．040．960．04

2.2.2 設(shè)計結(jié)果分析

圖8 不同激波角下升阻比和體積效率隨Δw/s的變化曲線Fig.8 Lift to drag ratio and volumetric efficiency vs Δw/s at different shock angles

在每個特定Δw/s下，OCWRCAS II隨激波角的變化趨勢與圖3(a)一致，并且都是在11° 激波角時出現(xiàn)升阻比最大的情況，本節(jié)不再給出OCWRCAS II隨激波角變化曲線，重點分析Δw/s大小對乘波體性能的影響。以75° 后掠角的設(shè)計結(jié)果為代表進(jìn)行研究，圖8給出了該型乘波體的升阻比和體積效率隨Δw/s的變化曲線，圖中曲線根據(jù)激波角的大小分為3類：小于11°，升阻比隨激波角的增大而增大，圖中用虛線畫出；大于11°，升阻比隨激波角的增大而減小，圖中用點線畫出；等于11°，圖中用實線表示。圖8(a)也給出了升阻比隨Δw/s的變化曲線，小于11° 的曲線升阻比隨著Δw/s的增大先增后減，減小幅度隨激波角的增大逐漸降低，11° 時的升阻比也為先增大后減小，但減小趨勢非常低，到大于11° 時就變?yōu)殡S著Δw/s的增大而增大。圖8(b)給出了體積效率隨Δw/s的變化曲線，其變化規(guī)律較升阻比相比要簡單，都是隨著Δw/s的增大而減小。圖8中的結(jié)果說明Δw/s的增大對升阻比而言有有利情況也有不利情況，對體積效率則完全是不利的，因此要根據(jù)上述結(jié)論在設(shè)計空間中選擇“最優(yōu)”的結(jié)果很困難，還需要給出升阻比隨體積效率的變化趨勢。

同樣給出OCWRCAS II的升阻比隨體積效率的變化規(guī)律，如圖9所示，得到類似Pareto最優(yōu)解的前沿線，占據(jù)前沿線最大比例的為后掠角為75° 的乘波體，這個結(jié)論與OCWRCAS I一致，說明對于定后掠角密切錐乘波體，后掠角最大的乘波體是升阻比和體積效率都較高的一類乘波體，具有最高的實用價值。

如圖10所示，本節(jié)使用與OCWRCAS I相同的選擇標(biāo)準(zhǔn)，對OCWRCAS II進(jìn)行分析。表4給出了升阻比最大、體積效率最大和本文選擇結(jié)果的性能參數(shù)，升阻比最大的乘波體生成激波角為12°，后掠角為50°，Δw/s為最大值0.96，升阻比可達(dá)9.05，圖10(a)給出了升阻比最大的OCWRCAS II三視圖，從三視圖可以看出該型乘波體非常接近錐導(dǎo)乘波體，已經(jīng)沒有明顯的后掠前緣，而且厚度非常小，所以有很大的升阻比，但體積效率僅有0.083 9；體積效率最大的乘波體后掠角同樣為50°，生成激波角為最大激波角 39.5°，Δw/s為0.04，圖10(b)給出了體積效率最大的OCWRCAS II三視圖，該外形曲線型頭部非常小，具有與圓錐接近的外形，所以具有很大的體積效率，但升阻比只有1.15；本節(jié)選擇的“最優(yōu)”結(jié)果的生成激波角為12.5°，后掠角為75°，Δw/s為0.36，該型乘波體的三視圖如圖10(c)所示，其升阻比為7.15，優(yōu)于OCWRCAS I選擇的外形，但體積效率較OCWRCAS I要低。

圖9 OCWRCAS II升阻比隨體積效率的變化規(guī)律Fig.9 Lift to drag ratio vs volumetric efficiency for OCWRCAS II

圖10 OCWRCAS II 3種外形的三視圖Fig.10 Three views of three types of OCWRCAS II

表4 OCWRCAS II優(yōu)化設(shè)計結(jié)果Table 4 Optimal design results of OCWRCAS II

3 CFD驗證

3.1 乘波特性

以O(shè)CWRCAS II中的“最優(yōu)”乘波體為例，使用CFD方法進(jìn)行黏性流場分析，驗證乘波體的乘波氣動特性。氣動力結(jié)果比較如表5所示，可見CFD計算結(jié)果與理論結(jié)果是相近的。

圖11給出乘波體下表面流場中的壓力等值線，說明該型乘波體激波整體附著在前緣，使流場在前緣處泄露非常少，具有很好的乘波特性。圖12 給出了乘波體底部處的壓力分布等值線以及設(shè)計該型乘波體時的ICC曲線，ICC即底部的理論激波形狀，由等值線反映出的激波形狀與ICC符合得很好，說明本文給出的設(shè)計方法是一種有效的密切錐乘波體設(shè)計方法。文獻(xiàn)[21]通過計算指出，雖然下表面的理論結(jié)果存在一定誤差，但由于數(shù)值差別較小而且面積不大，最終的升力系數(shù)和阻力系數(shù)仍然能與CFD計算結(jié)果保持一致，所以使用理論結(jié)果進(jìn)行設(shè)計是有效的。

表5“最優(yōu)”O(jiān)CWRCASII的氣動力結(jié)果比較

Table5Comparisonofaerodynamicforceresultsof“optimal”O(jiān)CWRCASII

MethodCL(CD)inviscous(CD)viscousCL/CDTheory0．02030．00160．00127．15CFD0．02090．00170．00146．74

圖11 設(shè)計狀態(tài)時下表面流場的壓力等值線 Fig.11 Pressure contour at lower surface with design status

圖12 理論激波形狀(”)與CFD計算結(jié)果(壓力等值線)比較Fig.12 Comparison of shock wave between theoretical shape (”) and CFD (pressure contour) computational results

3.2 渦升力

本節(jié)使用CFD方法分別對OCWRCAS II的“最優(yōu)”外形進(jìn)行了研究。計算狀態(tài)為：30 km，馬赫數(shù)6，全湍流。計算迎角α為0°、2°、4°、6°、8°、10°和12°。OCWRCAS II上表面的壓力分布系數(shù)Cp如圖13所示，可知該型乘波體隨著迎角增大，上表面壓力漸趨減小，由渦誘導(dǎo)的低壓區(qū)域逐漸增大。圖14給出了升力系數(shù)隨迎角變化曲線，在6° 迎角之后，升力曲線呈現(xiàn)明顯非線性增加趨勢。本節(jié)計算說明，定后掠角密切錐乘波體在大迎角時具有高升力特性。但是高超聲速時，飛行器上表面的分離情況是非常復(fù)雜的，其生成、發(fā)展機(jī)理以及影響因素有待進(jìn)一步研究。

圖13 OCWRCAS II “最優(yōu)”外形不同迎角時上表面壓力云圖Fig.13 Pressure contour at upper surface of “Optimal” OCWRCAS II shape with different angles of attack

圖14 OCWRCAS II “最優(yōu)”外形升力系數(shù)隨迎角變化曲線 Fig.14 Lift coefficient of “Optimal” OCWRCAS II vs angle of attack

4 結(jié) 論

1) 本文給出的OCWRCAS定后掠角前緣對應(yīng)下表面的生成方法切實有效，得到的乘波體乘波特性明顯。

2) 后掠角最大的定后掠角密切錐乘波體，絕大部分都位于“Pareto前沿線”上，即這類乘波體處于升阻比和體積效率均“較優(yōu)”的位置。

3) 對于某一后掠角的定后掠角密切錐乘波體，升阻比隨激波角的增大呈現(xiàn)先增大后減小的趨勢，體積效率隨激波角的增大逐漸增大。

4) 對于OCWRCAS II，給定激波角和后掠角時，升阻比隨前緣曲線程度的變化趨勢與激波角的大小有關(guān)，體積效率則是隨前緣曲線程度的增大而增大。

5) 通過對多個迎角的CFD計算，說明定后掠角密切錐乘波體在一定迎角時，上表面會產(chǎn)生穩(wěn)定的分離，從而改善升力特性。

6) 本文在對定后掠角密切錐乘波體進(jìn)行設(shè)計時，組成ICC的曲線部分使用了圓弧，限制了設(shè)計空間，若采用一般曲線進(jìn)行優(yōu)化設(shè)計有望得到性能更優(yōu)的乘波體。

[1] NONWEILER T R F. Aerodynamic problems of manned space vehicles[J]. Journal of the Royal Aeronautical Society, 1959, 63(585): 521-528.

[2] JONES J G, MOORE K C, PIKE J, et al. A method for designing lifting configurations for high supersonic speeds using axisymmetric flow field[J]. Archive of Applied Mechanics, 1968, 37(1): 56-72.

[3] 耿永兵, 劉宏, 姚文秀, 等. 錐形流乘波體優(yōu)化設(shè)計研究[J]. 航空學(xué)報, 2006, 27(1): 23-28.

GENG Y B, LIU H, YAO W X, et al. Viscous optimized design of waverider derived from cone flow [J]. Acta Aeronautica et Astronautica Sinica, 2006, 27(1): 23-28 (in Chinese).

[4] 耿永兵, 劉宏, 雷麥芳, 等. 高升阻比乘波構(gòu)型優(yōu)化設(shè)計[J]. 力學(xué)學(xué)報, 2006, 38(4): 540-546.

GENG Y B, LIU H, LEI M F, et al. Optimized design of waverider with high lift over drag ratio[J]. Chinese Journal of Theoretical and Applied Mechanics, 2006, 38(4): 540-546 (in Chinese).

[5] RASMUSSEN M L. Waverider configurations derived from inclined circular and elliptic cones[J]. Journal of Spacecraft and Rockets, 1980, 17(6): 537-545.

[6] 樂貴高, 馬大為, 李自勇. 橢圓錐乘波體高超聲速流場數(shù)值計算[J]. 南京理工大學(xué)學(xué)報(自然科學(xué)版), 2006, 30(3): 257-260.

LE G H, MA D W, LI Z Y. Computation of hypersonic flowfields for elliptic-cone-derived waverider[J]. Journal of Nanjing University of Science and Technology (Natural Science), 2006, 30(3): 257-260 (in Chinese).

[7] TAKASHIMA N, LEWIS M J. Waverider configurations based on non-axisymmetric flow fields for engine-airframe integration: AIAA-1994-0380[R]. Reston: AIAA, 1994.

[8] SOBIECZKY H, DOUGHERTY F C, JONES K. Hypersonic waverider design from given shock wave[C]//First International Waverider Symposium. Maryland: University of Maryland, 1990.

[9] CENTER K, SOBIECZKY H, DOUGHERTY F C. Interactive design and analysis of hypersonic waverider geometries: AIAA-1991-1697[R]. Reston: AIAA, 1991.

[10] KONTOGIANNIS K, SOBESTER A, TAYLOR N J. On the conceptual design of waverider forebody feometries: AIAA-2015-1009[R]. Reston: AIAA, 2015.

[11] SZEMA K, LIU Z, MUNIPALLI R. An efficient GUI design tool for high-speed airbreathing propulsion integration: AIAA-2010-4362[R]. Reston: AIAA, 2010.

[12] 賀旭照, 倪鴻禮. 密切曲面錐乘波體-設(shè)計方法與性能分析[J]. 力學(xué)學(xué)報, 2011, 43(6): 1077-1082.

HE X Z, NI H L. Osculating curved cone(OCC) waverider: Design methods and performance analysis[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(6): 1077-1082 (in Chinese).

[13] 賀旭照, 周正, 毛鵬飛, 等. 密切曲面內(nèi)錐乘波前體進(jìn)氣道設(shè)計和試驗研究[J]. 實驗流體力學(xué), 2014, 28(3): 39-44.

HE X Z, ZHOU Z, MAO P F, et al. Design and experimental study of osculating inward turning cone waverider/inlet(OICWI)[J]. Journal of Experiments in Fluid Mechanics, 2014, 28(3): 39-44 (in Chinese).

[14] BOWCUTT K G. Optimization of hypersonic waveriders derived from cone flows-including viscous effects[D]. Maryland： University of Maryland, 1986.

[15] RODI P E. The osculating flowfield method of waverider geometry generation： AIAA-2005-0511[R]. Reston: AIAA, 2005.

[16] RODI P E. Geometrical relationships for osculating cones and osculating flowfield waverider: AIAA-2011-1188[R]. Reston: AIAA, 2011.

[17] RODI P E. Vortex lift waverider configurations: AIAA-2012-1238[R]. Reston: AIAA, 2012.

[18] ANDERSON J D. Modern compressible flow[M]. 2nd ed. New York: McGraw-Hill Publishing Company, 1999: 294-307.

[19] CORDA S, ANDERSON J. Viscous optimized hypersonic waveriders designed from axisymmetric flow fields: AIAA-1988-0369[R]. Reston: AIAA, 1988.

[20] DIETER J, GOTTFRIED S, SIEGFRIED W. Basic research and technologies for two-stage-to-orbit vehicles[M]. Weinheim: WILEY-VCH Verlag Gmbh & Co. KGaA, 2005: 215-260.

[21] CENTER K, SOBIECZKY H, DOUGHERTY F. Interactive design of hypersonic waverider geometries: AIAA-1991-1697[R].Reston: AIAA, 1991.

段焰輝男， 博士。主要研究方向： 計算流體力學(xué)， 氣動外形優(yōu)化設(shè)計。

E-mail: duanyanhui@foxmail.com

吳文華男， 研究員。主要研究方向： 空氣動力學(xué)， 氣動外形優(yōu)化設(shè)計。

Tel.: 0816-2463132

E-mail: 619677947@qq.com

URL：www.cnki.net/kcms/detail/11.1929.V.20160226.1344.002.html

Generationanddesignmethodsofosculatingconewaveriderwithconstantangleofsweepback

DUANYanhui,FANZhaolin,WUWenhua*

ComputationalAerodynamicsResearchInstitute，ChinaAerodynamicsResearchandDevelopmentCenter，Mianyang621000，China

Inthispaper,thegenerationanddesignmethodsofosculatingconewaveriderwithconstantangleofsweepback(OCWRCAS)arestudied,andviscousisconsideredduringdesign.Stablevortexcanbegeneratedbytheleadingedgewithconstantangleofsweepback,whichwillimprovetheaerodynamicabilityofthewaverider.Firstly,thegenerationmethodofOCWRCASispresentedbasedonthegenerationmethodoftraditionalosculatingconewaverider.Thedesignvariablesofsweepbackangle,shockangleandthecurveshapeofheadareextractedbyanalyzingthegeometrycharacteroftheOCWRCAS,andthevariationtrendoflifttodragratioandvolumetricefficiencywiththesevariablesisalsostudied.Themulti-objectiveoptimalsolutionsarefoundedfromtwoclassicaltypesofOCWRCASbysearchingthetotaldesignspace.Finally,themethodofcomputationalfluiddynamicsisusedtoverifythecharacterofwaveridingandvortexlift.TheresultsshowthatOCWRCASwithgoodabilityofwaveridingandhighlifttodragratiokeepsrelationalvolumetricefficiency;vortexliftcanbegeneratedbytheleadingedgewithconstantangleofsweepbackatcertainangleofattack.

waverider;osculatingcone;viscous;angleofsweepback;vortexlift

2015-12-14；Revised2015-12-30；Accepted2016-01-18；Publishedonline2016-02-261344

.Tel.：0816-2463132E-mail619677947@qq.com

2015-12-14；退修日期2015-12-30；錄用日期2016-01-18； < class="emphasis_bold">網(wǎng)絡(luò)出版時間

時間：2016-02-261344

www.cnki.net/kcms/detail/11.1929.V.20160226.1344.002.html

.Tel.：0816-2463132E-mail619677947@qq.com

段焰輝， 范召林， 吳文華． 定后掠角密切錐乘波體的生成和設(shè)計方法J． 航空學(xué)報,2016,37(10):3023-3034.DUANYH,FANZL,WUWH.GenerationanddesignmethodsofosculatingconewaveriderwithconstantangleofsweepbackJ.ActaAeronauticaetAstronauticaSinica,2016,37(10):3023-3034.

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

10.7527/S1000-6893.2016.0024

V211.5

A

1000-6893(2016)10-3023-12