黃喜元，李小艷，楊 勇，陳 智，苗文博

（1.中國(guó)運(yùn)載火箭技術(shù)研究院，北京，100076；2.中國(guó)航天空氣動(dòng)力技術(shù)研究院，北京，100074）

0 引 言

熱防護(hù)一直是高超聲速飛行器攻關(guān)研究中的一項(xiàng)核心關(guān)鍵技術(shù)問(wèn)題[1,2]。傳統(tǒng)高超聲速飛行器熱防護(hù)系統(tǒng)設(shè)計(jì)主要采用被動(dòng)熱防護(hù)設(shè)計(jì)，或者依靠防熱結(jié)構(gòu)和材料本身將熱量吸收或輻射出去，如剛性陶瓷防熱瓦、柔性氈式防熱、高導(dǎo)熱碳復(fù)合材料防熱、蓋板式防熱等，達(dá)到保護(hù)飛行器內(nèi)部結(jié)構(gòu)的目的；或者通過(guò)自身材料質(zhì)量的損失，如燒蝕熱防護(hù)，吸收并帶走大量的熱，減少熱量向內(nèi)傳遞，從而保護(hù)飛行器內(nèi)部結(jié)構(gòu)能夠正常工作。受材料體系限制，依靠研制新型的耐高溫材料來(lái)解決高超聲速飛行器的熱防護(hù)問(wèn)題不僅技術(shù)難度越來(lái)越大，研制周期長(zhǎng)、研制成本高也成為制約飛行器性能進(jìn)一步提升的瓶頸問(wèn)題。

主動(dòng)熱防護(hù)技術(shù)由于冷卻能力強(qiáng)、可以適應(yīng)長(zhǎng)時(shí)間高熱流環(huán)境、設(shè)計(jì)靈活，可滿足極端苛刻、嚴(yán)酷的氣動(dòng)加熱環(huán)境，近年來(lái)受到了極大的重視。作為主動(dòng)熱防護(hù)技術(shù)之一的逆向噴流方法，由于可以顯著降低飛行器的氣動(dòng)加熱，成為高超聲速飛行器主動(dòng)熱防護(hù)系統(tǒng)設(shè)計(jì)可選擇的有效途徑之一。逆向噴流降熱技術(shù)的主要原理是在飛行器頭部設(shè)置一個(gè)噴口，噴口噴射出來(lái)的氣體將飛行器頭部弓形激波更遠(yuǎn)地推離飛行器壁面，從而降低氣動(dòng)加熱熱量向飛行器內(nèi)部傳遞，達(dá)到主動(dòng)改善飛行器氣動(dòng)加熱熱環(huán)境的目的。

本文研究高超聲速飛行器逆向噴流降熱技術(shù)，針對(duì)某高超聲速飛行器頭部，開展了單噴管逆向噴流降熱特性和多噴管逆向噴流降熱技術(shù)的研究。重點(diǎn)呈現(xiàn)了使用單噴管逆向噴流和多噴管逆向噴流實(shí)現(xiàn)高超聲速飛行器頭部降熱的一些重要特性和結(jié)論，反映了使用逆向噴流降低高超聲速飛行器頭部熱流的技術(shù)可行性，單噴管逆向噴流降熱的局限性和多噴管逆向噴流降熱的魯棒性和工程適應(yīng)性等。

1 逆向噴流降熱技術(shù)的主要研究成果

1.1 逆向噴流降熱機(jī)理

日本九州大學(xué)Shigeru Aso 教授等[3～5]在超聲速流場(chǎng)中逆向噴流的減阻降熱試驗(yàn)中揭示了逆向噴流降熱的作用機(jī)理，如圖1 所示。

圖1 逆向噴流在超聲速主流中的流場(chǎng)結(jié)構(gòu)[1]Fig.1 Flowfield Structure of Counterflowing Jet to Supersonic Mainstream

逆向噴流從噴口噴出后，先形成一個(gè)馬赫盤，以平衡噴流壓力與弓形激波后的來(lái)流壓力。噴流與來(lái)流相遇形成接觸面，逆向噴流受到自由來(lái)流阻擋而產(chǎn)生回流，重新附于物面，從而形成回流區(qū)，而自由來(lái)流由于噴流的阻擋而向周邊流去，在回流再附點(diǎn)附近形成可壓縮激波。由于噴流將脫體激波更遠(yuǎn)地推離壁面，被噴流覆蓋的區(qū)域氣動(dòng)加熱環(huán)境得到明顯的改善，如圖2 所示。

圖2 逆向噴流對(duì)稱面流線、馬赫數(shù)云圖和溫度云圖Fig.2 Contours of Mach Number and Temperature with Counterflowing Jet

逆向噴流在降低飛行器頭部熱流的同時(shí)，也會(huì)帶來(lái)流場(chǎng)的振蕩和不穩(wěn)定問(wèn)題。Venkatachari 等[6～9]針對(duì)鈍頭體駐點(diǎn)附近噴管逆向噴流的流場(chǎng)結(jié)構(gòu)進(jìn)行了深入的研究，鈍頭體弓形激波與逆向噴流兩者間作用的流場(chǎng)結(jié)構(gòu)模態(tài)與逆向噴流質(zhì)量流率密切相關(guān)，隨著噴流質(zhì)量流率的增大，相互作用的流場(chǎng)逐漸變得不穩(wěn)定，流場(chǎng)結(jié)構(gòu)呈現(xiàn)出長(zhǎng)穿透模態(tài)（Long Penetration Mode，LPM），當(dāng)噴流質(zhì)量流率超過(guò)了某個(gè)臨界值，相互作用的流場(chǎng)又退化為接近穩(wěn)定狀態(tài)，流場(chǎng)結(jié)構(gòu)由長(zhǎng)穿透模態(tài)形態(tài)轉(zhuǎn)換為短穿透模態(tài)（Short Penetration Mode，SPM）[6]。

1.2 逆向噴流降熱主要研究成果

綜合國(guó)內(nèi)外研究成果，在逆向噴流降熱方面，大體有以下研究結(jié)論：

a）逆向噴流的降熱效果由激波脫體距離決定[10]。

b）由文獻(xiàn)[11]的研究發(fā)現(xiàn)，噴流氣體的物理屬性是影響降熱效果的重要因素之一，氦氣由于具有高擴(kuò)散系數(shù)，采用氦氣噴流可形成更大的回流區(qū)，降熱效果更好。文獻(xiàn)[12]的研究則顯示，氮?dú)庀啾群饨禑嵝Ч谩?/p>

c）文獻(xiàn)[13]的研究顯示，將單噴管噴流設(shè)計(jì)成微間距陣列噴管噴流，在相同的總質(zhì)量流量條件下，由于多噴管噴流可以更容易地?cái)U(kuò)散到邊界層，能夠?qū)崿F(xiàn)更好的降熱效果。

2 單噴管逆向噴流降熱特性研究

選取某高超聲速飛行器頭部開展逆向噴流降熱特性分析研究，飛行器頭部球頭半徑為56 mm，噴管布置于飛行器頭部幾何駐點(diǎn)處，研究了Φ3 mm、Φ6 mm、Φ12 mm 3 種噴口直徑的降熱效果，如圖3 所示。計(jì)算工況選為高度45 km、馬赫數(shù)10、攻角0°、側(cè)滑角0°。此工況下波后總壓為19.26 kPa，考慮噴流參數(shù)的影響，選取了20 kPa、40 kPa、60 kPa 3 個(gè)噴流總壓進(jìn)行降熱特性分析，噴流總溫取為300 K，噴流工質(zhì)為空氣。

圖3 單噴管逆向噴流降熱研究模型Fig.3 Models of Single Counterflowing Jet for a Nose Cap

氣動(dòng)熱特性計(jì)算采用數(shù)值求解可壓縮雷諾平均Navier-Stokes 方程的CFD 方法；空間離散采用Roe 的FDS 格式，時(shí)間離散采用有限體積半離散格式，無(wú)粘通量采用隱式格式，粘性通量采用顯式方法。研究發(fā)現(xiàn)逆向噴流的流場(chǎng)結(jié)構(gòu)呈現(xiàn)出明顯非定常效應(yīng)，理論上對(duì)于非定常流場(chǎng)計(jì)算應(yīng)采用時(shí)間精確的顯式迭代方法或者雙時(shí)間步進(jìn)行模擬，考慮到巨大的計(jì)算量，本文采用定常計(jì)算方法獲得基本穩(wěn)定解后對(duì)若干迭代步計(jì)算結(jié)果進(jìn)行平均的方法獲得流場(chǎng)結(jié)構(gòu)以及面熱流、壓力。

圖4 給出了Φ3 mm 噴管無(wú)噴流時(shí)對(duì)稱面馬赫數(shù)及表面熱流，無(wú)噴流時(shí)頭部熱流峰值約為1200 kW/m2。圖5 給出了Φ3 mm 噴管各噴流總壓下表面熱流分布。由圖5 可知，噴流總壓為20 kPa 時(shí)，對(duì)稱面激波形狀觀察不到明顯變化，頭部熱流較無(wú)噴流工況略有減小，隨著噴流總壓的增大，駐點(diǎn)區(qū)域激波脫體距離逐漸增大，駐點(diǎn)熱流逐漸減小，但無(wú)干擾區(qū)域熱流變化不大。

圖4 Φ3mm 無(wú)噴流時(shí)的降熱特性Fig.4 Heat Flux Reduction Characteristics without Jet (Φ3mm)

圖5 Φ3mm 各噴流總壓對(duì)稱面表面熱流Fig.5 Surface Heat Flux with Differen Jet Total Pressure(Φ3mm)

為進(jìn)一步量化分析噴流總壓、噴管直徑對(duì)駐點(diǎn)區(qū)域降熱影響，圖6～8 分別給出了Φ3 mm、Φ6 mm、Φ12 mm 噴管各噴流總壓對(duì)稱面中心線熱流及壓力分布。由圖6 至圖8 可知，相同噴流總壓條件下隨著噴管直徑增大，駐點(diǎn)區(qū)域熱流顯著降低，非駐點(diǎn)區(qū)域峰值熱流也降幅顯著；相同噴管直徑條件下隨著噴流總壓提高，駐點(diǎn)區(qū)域熱流降幅顯著，但局部非駐點(diǎn)局域呈現(xiàn)不同規(guī)律：較小噴管直徑（Φ3 mm、Φ6 mm）條件下非干擾區(qū)峰值熱流隨距離理論駐點(diǎn)位置的距離增大先降低再增高；較大噴管直徑（Φ12 mm）條件下非干擾區(qū)峰值熱流呈單調(diào)降低趨勢(shì)。從壓力分布上也可以觀察到，噴管周邊存在明顯高壓再附帶，而且隨著噴流總壓提高再附愈發(fā)明顯；Φ12 mm 噴管直徑由于噴流質(zhì)量流率增大，再附影響稍微減弱。

圖6 Φ3mm 對(duì)稱面中心線熱流及壓力分布Fig.6 Distributions of Surface Heat Flux and Pressure Along Middle Line (Φ3mm)

圖7 Φ6mm 對(duì)稱面中心線熱流及壓力分布Fig.7 Distributions of Surface Heat Flux and Pressure Along Middleline (Φ6mm)

圖8 Φ12mm 對(duì)稱面中心線熱流及壓力分布Fig.8 Distributions of Surface Heat Flux and Pressure Along Middleline (Φ12mm)

綜上所述，對(duì)于完全逆向噴流工況，大噴口直徑、大噴流總壓噴流可以顯著降低駐點(diǎn)表面熱流。為了進(jìn)一步研究單噴管逆向噴流降熱的適應(yīng)性，進(jìn)一步針對(duì)Φ6 mm、Φ12 mm 噴口直徑開展了高度45 km、馬赫數(shù)10、攻角15°、側(cè)滑角0°工況的計(jì)算分析研究。

圖9、圖10 分別給出了15°攻角下Φ6 mm、Φ12 mm噴口直徑不同噴流總壓下對(duì)稱面中心線熱流分布?？梢钥吹?，在噴流不完全正對(duì)來(lái)流的工況下，由于幾何駐點(diǎn)與物理駐點(diǎn)不重合，物理駐點(diǎn)處的熱流與無(wú)噴流工況相當(dāng)，而且噴流孔周圍還出現(xiàn)了高于駐點(diǎn)熱流高熱流帶，且在相同的噴流總壓下，隨著噴口直徑的增加，高熱流帶的范圍也進(jìn)一步增大。

圖9 Φ6mm 對(duì)稱面中心線熱流分布（攻角影響）Fig.9 Distributions of Surface Heat Flux Along the Middleline(Φ6mm,AOA’s Effects)

圖10 Φ12mm 對(duì)稱面中心線熱流分布（攻角影響）Fig.10 Distributions of Surface Heat Flux Along the Middleline(Φ12mm,AOA’s Effects)

分析發(fā)現(xiàn)，當(dāng)噴流不正對(duì)來(lái)流時(shí)，噴流甚至?xí)鸶哂跓o(wú)噴流工況下的干擾熱流，而且增大噴流孔直徑會(huì)引起干擾熱流帶面積的增大。因此，將逆向噴流降熱技術(shù)應(yīng)用在高超聲速飛行器上時(shí)，采用單噴管逆向噴流方案的適應(yīng)性是非常有限的。

3 多噴管逆向噴流降熱特性研究

仍然選取上述分析采用的某高超聲速飛行器頭部開展研究，在飛行器頭部球頭上對(duì)稱的布置多個(gè)噴管，單個(gè)噴管噴口直徑為Φ3 mm，計(jì)算工況選為高度50 km、馬赫數(shù)15、攻角15°、側(cè)滑角0°，該工況無(wú)噴流時(shí)頭部熱流峰值約為3100 kW/m2，噴流總溫取為300 K，噴流工質(zhì)為空氣。

19 個(gè)噴管均沿飛行器頭部球頭表面法向方向布局，其中噴管1 位于頭部幾何中心頂點(diǎn)，噴管1、10、11 中心線位于頭部縱對(duì)稱面內(nèi)，噴管1、2、7、12、17 位于頭部橫對(duì)稱面內(nèi)，噴管2、3、4、5、6、7、8、9 與噴管12、13、14、15、16、17、18、19 相對(duì)縱對(duì)稱面對(duì)稱。定義噴管軸線與球頭頂點(diǎn)法線夾角為θ，定義噴管軸線在頭部橫截面內(nèi)的投影與橫對(duì)稱面的夾角為φ，各噴管詳細(xì)的布局如表1 所示。

表1 噴管詳細(xì)布局Tab.1 Detailed Layout of Nozzles

圖11 給出了多噴管噴流總壓分別取為0.45 MPa和0.9 MPa 下多噴管逆向噴流下頭部對(duì)稱面馬赫數(shù)云圖，可以看出，噴流總壓0.45 MPa 下頭部表面最大熱流約為1600 kW/m2，噴流總壓0.9 MPa 下頭部表面最大熱流約為1100 kW/m2，2 種噴流總壓下，激波均被顯著的推離壁面，表面熱流均顯著減小。

圖11 多噴管逆向噴流下頭部對(duì)稱面馬赫數(shù)云圖Fig.11 Contours of Mach Number with Micro Counterflowing Jets Array

噴流總壓取為0.9 MPa，圖12 分別給出了高度50 km、馬赫數(shù)15、攻角30°、側(cè)滑角0°工況和高度50 km、馬赫數(shù)15、攻角15°、側(cè)滑角5°工況下的多噴管逆向噴流下頭部的表面熱流，2 種工況下頭部最大熱流密度分別不超過(guò)100 kW/m2和550 kW/m2，顯示該多噴管布局對(duì)于大攻角和有側(cè)滑的情況下均獲得了顯著的降熱效果。

圖12 多噴管逆向噴流下頭部的表面熱流Fig.12 Contours of Surface Heat Flux with Micro Counterflowing Jets Array

綜合上述的研究結(jié)果表明，相比單噴管逆向噴流降熱，多噴管逆向噴流降熱具有更好的魯棒性，更能滿足工程應(yīng)用的需求。

4 結(jié)束語(yǔ)

本文對(duì)高超聲速飛行器逆向噴流技術(shù)進(jìn)行研究，介紹了逆向噴流的降熱機(jī)理，簡(jiǎn)要總結(jié)了國(guó)內(nèi)外在逆向噴流降熱技術(shù)方面取得重要研究成果，并針對(duì)某高超聲速飛行器頭部，開展了單噴管逆向噴流降熱特性和多噴管逆向噴流降熱技術(shù)的研究，取得了以下成果：

a）揭示了完全逆向噴流、大噴口直徑、大噴流總壓（對(duì)應(yīng)大的噴流來(lái)流壓力比和大的噴流質(zhì)量流率）噴流可以顯著降低駐點(diǎn)表面熱流，但同時(shí)也可以看出，單噴管逆向噴流降熱對(duì)于高超聲速飛行器飛行攻角等有明確的限制，反映單噴管逆向噴流降熱工程適應(yīng)性較低。

b）通過(guò)對(duì)多噴管逆向噴流降熱特性的研究，顯示多噴管逆向噴流降熱是一種更為魯棒的降熱技術(shù)，具有較好的工程適應(yīng)性，但同時(shí)也可以看出，多噴管逆向噴流的流場(chǎng)結(jié)構(gòu)非常復(fù)雜，給多噴流逆向噴流降熱的應(yīng)用帶來(lái)了一定的困難，包括多噴管噴流的布局設(shè)計(jì)問(wèn)題等。