莫 凡, 王鎖柱, 高振勛

(1. 北京航空航天大學(xué)航空科學(xué)與工程學(xué)院, 北京 100191;2. 北京航天長征飛行器研究所, 北京 100076)

引 言

高超聲速飛行器在大氣中進(jìn)行高M(jìn)ach數(shù)飛行時, 來流經(jīng)過頭部脫體激波的劇烈壓縮后往往能達(dá)到很高的溫度[1], 導(dǎo)致高溫空氣中發(fā)生復(fù)雜的化學(xué)反應(yīng). 一般情況下, 在2 500 K以上空氣中的氧氣開始離解, 4 000 K以上氮氣開始離解, 而到9 000 K以上氧原子和氮原子將進(jìn)一步電離.

當(dāng)離解之后的高溫空氣到達(dá)飛行器表面時, 表面材料會對空氣中的原子和離子產(chǎn)生顯著的催化復(fù)合作用, 這種催化復(fù)合作用將改變飛行器壁面附近各組分的質(zhì)量分?jǐn)?shù)梯度分布, 引起流場中輸運(yùn)特性的改變, 進(jìn)而顯著改變壁面熱流密度的分布. 已有研究表明, 在高M(jìn)ach數(shù)情況下關(guān)鍵部位防熱材料表面的催化復(fù)合效應(yīng)最高將帶來近50%以上的氣動熱載荷[2]. 因此在CFD模擬中, 催化邊界條件的確定往往非常重要, 完全非催化(non-catalytic wall, NCW)和完全催化(full-catalytic wall, FCW)[3]是高超聲速飛行器氣動熱環(huán)境數(shù)值模擬中最常用的兩種催化邊界條件. 完全非催化假設(shè)壁面催化復(fù)合速率為0, 即壁面組分質(zhì)量分?jǐn)?shù)的梯度為0. 而完全催化一般假設(shè)壁面來流原子的離子在壁面完全復(fù)合, 此時壁面組分質(zhì)量分?jǐn)?shù)的梯度也將最大.

然而, 以上兩種組分邊界條件均不能反映實際發(fā)生的壁面催化過程, 實際飛行中均為有限催化情況, 因此必須使用有限催化[3](finite-rate catalysis, FRC)邊界條件才能反映真實的壁面催化效應(yīng), 該邊界條件的關(guān)鍵在于催化復(fù)合系數(shù)的取值. 粟斯堯等[4]選取了催化復(fù)合系數(shù)為常數(shù)的有限催化邊界條件, 針對類聯(lián)盟號飛船返回艙外形研究了氣動熱環(huán)境隨壁面催化效率的變化規(guī)律, 并對壁面有限催化影響氣動加熱的物理機(jī)制進(jìn)行了探討. 丁明松等[5]開展了在有限催化條件下, 表面材料催化特性差異性對熱環(huán)境的影響, 得出壁面催化復(fù)合系數(shù)的差異性會帶來局部熱流的跳變, 使局部區(qū)域熱流明顯高于全表面催化復(fù)合系數(shù)為0.1的情況, 且壁面復(fù)合系數(shù)差異越大, 熱流跳變的程度越大. 然而對于實際的高超聲速飛行器表面催化效應(yīng), 其表面的復(fù)合系數(shù)取常數(shù)是不恰當(dāng)?shù)? Inger[6]的實驗表明催化復(fù)合系數(shù)往往與壁面材料和壁面溫度有關(guān), 并給出了某碳基材料催化復(fù)合系數(shù)隨表面溫度變化的擬合公式. Scott[7]開展了典型材料的催化復(fù)合系數(shù)研究, 并給出了2 000 K以下大部分碳基和硅基熱防護(hù)材料催化復(fù)合系數(shù)小于0.1的結(jié)論. Scott[8]和Zoby等[9]分別給出了實驗和飛行數(shù)據(jù)修正后氧原子和氮原子催化復(fù)合系數(shù)隨溫度擬合的公式. 國內(nèi)關(guān)于采用溫度擬合催化復(fù)合系數(shù)研究有限催化的工作尚少, 董維中等[10]采用源于實驗數(shù)據(jù)擬合得到的與壁溫相關(guān)的O原子的催化復(fù)合系數(shù), 開展了輻射平衡表面溫度分布與氣動熱的耦合計算, 但由于未考慮氮原子的催化復(fù)合效應(yīng), 因此模型并不完整.

本文在現(xiàn)有考慮離子組分的有限催化邊界條件的基礎(chǔ)上, 引入了由溫度擬合得到的氮氧原子催化復(fù)合系數(shù), 研究了在輻射平衡壁溫下有限催化模型對典型再入高超聲速飛行器氣動熱環(huán)境的影響規(guī)律. 文中首先給出了流動控制方程及有限催化邊界條件的數(shù)學(xué)形式, 在此基礎(chǔ)上建立了高超聲速飛行器化學(xué)非平衡流場有限催化氣動熱環(huán)境數(shù)值計算方法及程序, 并通過與風(fēng)洞實驗氣動熱結(jié)果對比驗證了數(shù)值模擬結(jié)果的可靠性. 在此基礎(chǔ)上考慮到有限催化中催化復(fù)合系數(shù)與壁面溫度相關(guān), 采用了Scott和Zoby給出的氧原子和氮原子催化復(fù)合系數(shù)隨溫度擬合的公式, 并研究鈍錐外形在輻射平衡壁溫下有限催化對壁面熱環(huán)境的影響規(guī)律.

1 計算方法

1.1 化學(xué)非平衡流動的瞬時N-S方程

在Descartes坐標(biāo)系下, 含化學(xué)反應(yīng)流動的瞬時N-S方程可寫成如下守恒形式

(1)

式中,ρ,ui,p,T分別為氣體的密度、 速度、 壓強(qiáng)與靜溫;E為單位(體積)質(zhì)量總能量

(2)

(3)

其中，μ為混合氣體的分子黏性系數(shù). 各組分分子黏性系數(shù)、 熱傳導(dǎo)系數(shù)、 質(zhì)量擴(kuò)散系數(shù)等輸運(yùn)系數(shù)和比熱、 比焓均采用Gupta的擬合公式得到[11], 而空氣混合物整體的輸運(yùn)系數(shù)則根據(jù)Wilke的公式計算得到[12].

1.2 數(shù)值求解方法

在對方程(1)無量綱化處理后, 采用有限差分方法進(jìn)行離散, 其中對流項采用Roe格式離散, 并采用Harten-Yee型熵修正, 且熵修正系數(shù)取0.25, 黏性項采用2階中心差分離散. 時間推進(jìn)上為提高CFL步長以及求解的穩(wěn)定性, 采用LU-SGS隱式方法, 并采用當(dāng)?shù)貢r間步長加速收斂. 對于三維復(fù)雜多塊結(jié)構(gòu)網(wǎng)格, 采用MPI分區(qū)并行求解技術(shù)以提高計算效率.

1.3 有限催化邊界條件

在高超聲速化學(xué)非平衡流動的CFD模擬中, 壁面催化效應(yīng)通過壁面邊界條件影響對應(yīng)的控制方程的解. 下面以7組分電離空氣為例, 給出了包含離子在內(nèi)的組分質(zhì)量分?jǐn)?shù)邊界條件.

對于包含離子組分的高溫氣體, 在壁面上發(fā)生的催化復(fù)合反應(yīng)主要有兩類, 即原子的復(fù)合和離子的復(fù)合. 現(xiàn)僅考慮以下3個催化復(fù)合反應(yīng)

N+N→N2; O+O→O2; NO++e-→NO

(4)

有限催化壁邊界條件一般通過求解壁面附近質(zhì)量守恒方程的方法構(gòu)建. 不考慮質(zhì)量引射時, 氣體組分由催化反應(yīng)生成(或消耗)的質(zhì)量與擴(kuò)散出(或入)壁面微元的質(zhì)量相等.

由催化反應(yīng)引起的單位時間, 單位面積壁面上各組分質(zhì)量消耗(生成)率為

(5)

其中，kws為反應(yīng)速率常數(shù). 假設(shè)壁面原子服從Maxwell-Boltzmann分布, 則kws與壁面催化復(fù)合系數(shù)γs具有如下關(guān)系[13]

(6)

式中,R0為普適氣體常數(shù);Tw為壁面溫度;Ms為組分s的摩爾質(zhì)量. 壁面催化復(fù)合系數(shù)γs表示壁面發(fā)生催化復(fù)合的原子數(shù)與入射到壁面總原子數(shù)之比. 其取值范圍在0到1之間, 取0時表示完全非催化(NCW), 取1時表示完全催化(FCW), 取0到1之間則為有限催化(FRC).

另一方面, 由Fick定律得壁面處由擴(kuò)散作用產(chǎn)生的質(zhì)量通量為

(7)

其中，n為壁面法向單位矢量, 方向由壁面指向流場內(nèi)部. 根據(jù)質(zhì)量守恒定律, 各組分單位時間單位面積由催化反應(yīng)產(chǎn)生(消耗)的質(zhì)量與擴(kuò)散出(到)壁面的質(zhì)量應(yīng)相等, 即各組分在壁面處的凈質(zhì)量通量為零. 故下述關(guān)系式成立

(8)

式(7)中存在導(dǎo)數(shù)項, 在CFD計算時需對其進(jìn)行恰當(dāng)?shù)碾x散, 本文中消耗組分(s=O,N,NO+)采用以下離散方式

(9)

整理得

(10)

(11)

整理得

(12)

1.4 輻射平衡壁溫邊界條件

輻射平衡模型依據(jù)的主要原理是在壁面處能量守恒. 飛行器壁面一方面要接收流場的熱流(壁面熱流密度)Qflux, 另一方面要向內(nèi)部結(jié)構(gòu)傳熱Qinside, 另外還要通過電磁輻射向外發(fā)射能量(灰體輻射能量)Qradiation.當(dāng)流場處于平衡狀態(tài)時, 上述3種能量的傳輸也應(yīng)滿足以下的平衡方程

Qflux=Qinside+Qradiation

(13)

現(xiàn)對結(jié)構(gòu)的內(nèi)部傳熱忽略不計, 平衡方程則可以進(jìn)一步化簡為

Qflux=Qradiation

(14)

即流場向壁面的傳熱量應(yīng)與壁面出射的輻射能量達(dá)到平衡. 當(dāng)流場氣體考慮組分時, 壁面熱流密度可以表示為對流項熱流密度Qconv和組分?jǐn)U散項熱流密度Qdiff之和,可由式(15)進(jìn)行計算

(15)

其中,n代表壁面的法向方向,ρ為壁面流場密度,ns為氣體組分個數(shù),Ds代表組分s的擴(kuò)散系數(shù),hs是組分s的絕對焓,Ys代表組分s的質(zhì)量分?jǐn)?shù).

若將壁面視為發(fā)射率等于ε的灰體, 則壁面出射的輻射能根據(jù)Stefan-Boltzmann定律, 可以寫作

(16)

其中,σ為Stefan-Boltzmann常數(shù), 取值為σ=5.67×10-8W/(m2·K4);Tw為壁面溫度. 聯(lián)立式(13), (15), 并將Qflux離散化, 則最終可得平衡方程

(17)

其中,T1和Tw分別為第1層網(wǎng)格和壁面處的溫度;Ys1和Ysw分別為第1層網(wǎng)格和壁面處組分s的質(zhì)量分?jǐn)?shù); Δy為第1層網(wǎng)格與壁面間的距離. 求解式(17)采用了Newton迭代法, 迭代關(guān)系如式(18)所示

(18)

2 等溫壁有限催化計算

卡爾斯本大學(xué)巴法羅研究中心(CUBRC)利用LENS系列激波風(fēng)洞對標(biāo)模返回艙進(jìn)行了一系列的研究[14]. 其中幾何尺寸如圖1所示, 實驗工況如表1所示. 根據(jù)實驗來流靜溫和Mach數(shù)可計算得到來流總溫超過8 000 K, 因此采用了Gupta 5組分6反應(yīng)化學(xué)反應(yīng)動力學(xué)模型, 壁面上僅考慮O原子和N原子的催化復(fù)合反應(yīng). 計算網(wǎng)格如圖2所示. 為了更好地捕捉頭部脫體激波, 在返回艙頭部前體附近對網(wǎng)格加密至0.2 mm.

圖1 返回艙模型幾何尺寸(單位： in[mm])Fig. 1 Geometric dimensions of re-entry capsule model(unit: in[mm])

表1 返回艙實驗工況

LENS激波風(fēng)洞實驗給出了標(biāo)模返回艙對稱面壁面上測量點的壓力與熱流測量值. 針對Test-1的CFD計算中組分邊界條件采用完全催化(FCW)、 完全非催化(NCW)以及γ分別取0.015和0.1的有限催化(FRC), 在圖3中與實驗熱流進(jìn)行了對比. 由于催化效應(yīng)幾乎沒有改變壁面的壓強(qiáng)分布, 因此壓強(qiáng)曲線在圖2中僅給出一條并與實驗結(jié)果對比.

針對 Test-2的CFD 計算中組分邊界條件采用完全催化(FCW)、 完全非催化(NCW)以及γ取0.005 的有限催化(FRC), 在圖4中與實驗壁面熱流及壓強(qiáng)進(jìn)行了對比. 圖5給出了壁面的壓強(qiáng)分布云圖. 圖6～8分別給出了完全非催化, 有限催化和完全催化邊界條件下的熱流密度分布云圖.

從Test-1與Test-2的CFD計算結(jié)果可以得到壁面熱流密度分布QFCW>QFRC>QNCW.這與以往經(jīng)驗和其他文獻(xiàn)結(jié)論相符.

圖2 返回艙模型計算網(wǎng)格Fig. 2 Three-dimensional structural grids for capsule geometry

圖3 Test-1返回艙壁面熱流與壓強(qiáng)分布Fig. 3 Heat flux and pressure distribution on the capsule wall in Test-1

圖4 Test-2返回艙壁面熱流與壓強(qiáng)分布Fig. 4 Heat flux and pressure distribution on the capsule wall in Test-2

圖5 Test-2返回艙壁面壓強(qiáng)分布云圖Fig. 5 Contour of pressure on the capsule wall in Test-2

圖6 Test-2返回艙壁面完全非催化熱流密度分布云圖Fig. 6 Contour of heat flux at the NCW in Test-2

圖7 Test-2返回艙壁面有限催化熱流密度分布云圖Fig. 7 Contour of heat flux at the FRC in Test-2

圖8 Test-2返回艙壁面完全催化熱流密度分布云圖Fig. 8 Contour of heat flux at the FCW in Test-2

而在使用有限催化模型時, 當(dāng)催化復(fù)合系數(shù)γ取合適值時可以得到與實驗結(jié)果更加吻合的結(jié)果. 這從側(cè)面反應(yīng)了高超聲速化學(xué)非平衡流動的壁面催化反應(yīng)是有限催化, 與以往人們的觀點相符. 此外, 當(dāng)催化復(fù)合系數(shù)γ取0時, 組分質(zhì)量分?jǐn)?shù)邊界條件自動退化為完全非催化壁條件.

由圖9分析Test-1壁面熱流組成可得, 壁面熱流由熱傳導(dǎo)引起的熱流和各組分質(zhì)量擴(kuò)散引起的熱流組成, 且在本算例中熱傳導(dǎo)引起的熱流大于各組分質(zhì)量擴(kuò)散引起的熱流. 在質(zhì)量擴(kuò)散引起的熱流中, 氧原子質(zhì)量擴(kuò)散引起的熱流密度占主要成分, 其他組分(如氮原子)僅僅占很小一部分, 這與各組分質(zhì)量分?jǐn)?shù)在壁面的梯度有關(guān). 以氧原子和氮原子為例, 如圖10所示, 由于激波后靜溫不到7 000 K, 氮氣分子離解程度遠(yuǎn)小于氧氣分子, 因此在壁面處氮原子質(zhì)量分?jǐn)?shù)梯度遠(yuǎn)小于氧原子質(zhì)量分?jǐn)?shù)梯度, 進(jìn)而由氮原子質(zhì)量擴(kuò)散引起的熱流密度小于氧原子質(zhì)量擴(kuò)散引起的熱流密度.

圖9 Test-1返回艙壁面熱流組成Fig. 9 Composition of the heat flux in Test-1

圖10 Test-1有限催化駐點線氧原子和氮原子質(zhì)量分?jǐn)?shù)分布圖Fig. 10 Mass fraction distribution of oxygen and nitrogrn along the stagnation line at the FRC in Test-1

圖11～12給出了Test-1有限催化氧分子沿駐點線質(zhì)量分?jǐn)?shù)分布, 可得壁面氧原子質(zhì)量分?jǐn)?shù)下降, 氧氣上升, 表明在催化邊界的作用下, 氧原子在壁面復(fù)合為氧氣. 完全催化的復(fù)合程度最大且壁面氧原子質(zhì)量分?jǐn)?shù)降到接近于0, 小于來流的氧原子質(zhì)量分?jǐn)?shù).

完全非催化條件下, 氧原子質(zhì)量分?jǐn)?shù)在壁面也有一定程度的降低. 這是由于壁面溫度僅有300 K, 激波后到壁面之間溫度逐漸下降, 在壁面附近溫度急劇下降, 而溫度影響了化學(xué)反應(yīng)速率進(jìn)而使當(dāng)?shù)亓鲃于呄蛴谛碌幕瘜W(xué)平衡, 因此由激波后高溫離解產(chǎn)生的原子在壁面附近自動復(fù)合.

由以上可以更清晰地得出, 壁面復(fù)合效應(yīng)有兩個組成部分, 第一是由壁面材料催化直接導(dǎo)致的復(fù)合, 第二是由壁面附近低溫導(dǎo)致的化學(xué)反應(yīng)速率變化而間接引起的復(fù)合.

圖11 Test-1有限催化氧原子沿駐點線質(zhì)量分?jǐn)?shù)Fig. 11 Mass fraction distribution of oxygen along the stagnation line at different catalytic conditions in Test-1

圖12 Test-1有限催化氧分子沿駐點線質(zhì)量分?jǐn)?shù)Fig. 12 Mass fraction distribution of oxygen molecule along the stagnation line at different catalytic conditions in Test-1

3 輻射平衡模型驗證

北卡羅萊納州立大學(xué)的 Keenan 等對高超聲速球繞流問題進(jìn)行了數(shù)值模擬研究[15], 且在文中考慮了壁面輻射的影響. 本節(jié)擬參考該數(shù)值模擬的結(jié)果, 對本文的輻射平衡模型進(jìn)行驗證. 計算的具體參數(shù)如下: 球的半徑R=1 m, 壁面的發(fā)射率ε=0.9, 海拔高度為H=65 km, 來流速度取 8 km/s.

圖13給出了高超聲速球頭繞流計算所采用的結(jié)構(gòu)網(wǎng)格, 并在激波附近對網(wǎng)格進(jìn)行了加密, 以便提高對激波的捕捉能力. 圖14給出了高超聲速球頭繞流壓強(qiáng)的分布云圖, 可以發(fā)現(xiàn)來流在經(jīng)過頭部脫體激波后壓強(qiáng)迅速增加, 其中在駐點附近由于激波強(qiáng)度最大, 因此壓強(qiáng)增加也最為顯著.

圖13 高超聲速球繞流網(wǎng)格Fig. 13 Structural grids for the hypersonic spherical body

圖14 高超聲速球繞流壓強(qiáng)分布云圖Fig. 14 Contour of pressure for the hypersonic spherical body

圖15, 16分別給出了球頭沿流向壁面上的溫度和熱流密度分布. 可以看出, 壁面溫度與文獻(xiàn)結(jié)果十分吻合, 熱流密度與文獻(xiàn)結(jié)果在中間部分有一定偏差, 但最大偏差不超過10%. 由此驗證了輻射平衡模型在計算壁面溫度分布時的正確性, 所以可以將此模型應(yīng)用于之后的飛行器外形溫度分布計算.

圖15 溫度分布, v=8 km/sFig. 15 Distribution of temperature along the wall， v=8 km/s

圖16 熱流密度分布, v=8 km/sFig. 16 Distribution of heat flux along the wall， v=8 km/s

4 輻射平衡壁溫有限催化計算

本節(jié)借助三維高超聲速鈍頭雙錐繞流實驗?zāi)Ｐ蚚16], 將該模型等比例放大, 放大之后的模型外形參數(shù)如下: 鈍錐體頭部的前緣半徑為 0.153 4 m, 前錐體的半錐角為12.84°, 后錐體的半錐角為7°, 兩錐的連接處與頭部的距離為2.782 m, 全錐總長4.889 6 m. 計算采用的結(jié)構(gòu)網(wǎng)格如圖17所示， 其中對雙錐頭部網(wǎng)格進(jìn)行了加密處理以便更精確地捕捉頭部脫體激波. 計算工況采用STS-2飛行實驗來流參數(shù), 如表2所示.

圖17 雙錐結(jié)構(gòu)網(wǎng)格Fig. 17 Structural grids for bi-cone geometry

表2 雙錐計算工況

由于Mach數(shù)等于26.3, 靜溫202.72 K, 總溫達(dá)到了15 000 K以上, 因此需要考慮電離反應(yīng), 采用Gupta 7組分9反應(yīng)化學(xué)反應(yīng)動力學(xué)模型, 壁面催化效應(yīng)考慮式(4)中的3個反應(yīng). CFD計算分別采用完全催化、 有限催化和完全非催化邊界條件. 其中， 有限催化模型的氧原子和氮原子催化復(fù)合系數(shù)采用Scott實驗數(shù)據(jù)擬合結(jié)果和Zoby飛行數(shù)據(jù)修正結(jié)果. 考慮到催化復(fù)合系數(shù)隨溫度變化的連續(xù)性, 因此本節(jié)最終采用式(19)的連續(xù)形式, 且γNO+取常數(shù).

(19)

圖18給出了雙錐輻射平衡壁溫下有限催化流場溫度云圖. 圖19～21分別給出了完全催化、 有限催化和完全非催化條件下駐點線各組分質(zhì)量分?jǐn)?shù)分布, 其中NO+和e-由于電離程度較低, 組分質(zhì)量分?jǐn)?shù)含量極少, 最終這兩種組分在圖中幾乎與x軸重合. 對于氧氣、 氮氣質(zhì)量分?jǐn)?shù), 壁面處完全催化大于有限催化大于完全非催化. NO質(zhì)量分?jǐn)?shù), 有限催化大于完全催化大于完全非催化,這與壁面溫度及各組分質(zhì)量分?jǐn)?shù)有關(guān).

圖18 雙錐輻射平衡壁溫下有限催化溫度云圖Fig. 18 Contour of temperature at the radiation balanced wall under FRC

圖19 完全催化駐點線組分質(zhì)量分?jǐn)?shù)分布Fig. 19 Mass fraction distribution of components along the stagnation line under FCW

圖20 有限催化駐點線組分質(zhì)量分?jǐn)?shù)分布Fig. 20 Mass fraction distribution of components along the stagnation line under FRC

圖21 完全非催化駐點線組分質(zhì)量分?jǐn)?shù)分布Fig. 21 Mass fraction distribution of components along the stagnation line under NCW

圖22給出了在輻射平衡邊界條件下, 完全催化, 有限催化和完全非催化駐點線溫度分布曲線. 表明在高空高M(jìn)ach情況下, 來流空氣稀薄, 激波厚度較厚, 同時壁面催化復(fù)合效應(yīng)也將改變激波位置. 其中完全非催化激波位置離壁面更遠(yuǎn), 且駐點溫度更低.

圖23, 24分別給出了完全催化、 有限催化和完全非催化壁面的溫度分布和熱流密度分布. 圖23表明有限催化下壁面溫度的峰值約為2 265 K, 而在大面積區(qū)域壁面溫度在1 150 K左右, 同時通過式(19)的催化復(fù)合系數(shù)計算公式求得在駐點附近氧原子的催化復(fù)合系數(shù)約為0.17, 氮原子的催化復(fù)合系數(shù)約為0.026, 而在大面積區(qū)域, 氧原子的催化復(fù)合系數(shù)約為0.005 3, 氮原子的催化復(fù)合系數(shù)約為 0.01. 另一方面對比了圖24中3種催化邊界的壁面熱流密度, 結(jié)果表明完全催化的壁面熱流峰值比有限催化高約21%, 而完全非催化的壁面熱流峰值比有限催化低約29%, 對于大面積的熱流分布, 完全催化的壁面熱流峰值比有限催化高約10%, 而完全非催化的壁面熱流峰值比有限催化低約8.7%.

圖22 完全催化、 有限催化和非催化駐點線溫度對比Fig. 22 Comparison of temperature along the stagnation line between FCW, FRC and NCW

圖23 完全催化、 有限催化和非催化壁面溫度對比Fig. 23 Comparison of wall temperature between FCW, FRC and NCW

圖24 完全催化、 有限催化和非催化壁面熱流對比Fig. 24 Comparison of heat flux between FCW, FRC and NCW

圖25給出了有限催化壁面總熱流、 熱傳導(dǎo)熱流與各組分熱流分布. 對比圖24， 25中有限催化和完全非催化的熱傳導(dǎo)熱流可以發(fā)現(xiàn), 有限催化條件下壁面的熱傳導(dǎo)熱流顯著增加, 這是由于在有限催化和輻射平衡耦合作用下, 使得激波位置相較于完全非催化發(fā)生了改變, 且更靠近壁面(如圖22所示), 從而使得壁面溫度梯度上升, 同時在這種耦合作用下壁面溫度也將上升, 進(jìn)而導(dǎo)致壁面熱傳導(dǎo)系數(shù)上升, 最終使得熱傳導(dǎo)熱流增強(qiáng). 圖26給出了有限催化壁面熱傳導(dǎo)熱流與各組分熱流占總熱流比值. 結(jié)合圖25分析可知, 在此來流條件下, 熱流密度中熱傳導(dǎo)引起的熱流占主要成分, 質(zhì)量擴(kuò)散引起的熱流中氧原子擴(kuò)散引起的熱流占主要成分, 且質(zhì)量擴(kuò)散引起的熱流在頭部駐點占比最大, 并隨著壁面離駐點的距離增大而下降.

圖25 有限催化壁面總熱流、 熱傳導(dǎo)熱流與各組分熱流對比Fig. 25 Comparison of total heat flux, thermal heat flux and heat flux of each component under FRC

圖26 有限催化壁面熱傳導(dǎo)熱流與各組分熱流占總熱流比值Fig. 26 Ratio of thermal heat flux and heat flux of each component to total heat flux under FRC

5 結(jié)論

本文發(fā)展了高超聲速飛行器輻射平衡壁溫下有限催化的數(shù)值方法, 基于返回艙外形的風(fēng)洞實驗數(shù)據(jù)進(jìn)行了數(shù)值模擬結(jié)果對比, 并進(jìn)一步針對典型高超聲速飛行器鈍雙錐研究了輻射平衡壁溫下有限催化對氣動熱環(huán)境的影響規(guī)律, 主要結(jié)論如下:

(1)針對返回艙外形的數(shù)值實驗結(jié)果表明, 完全催化與完全非催化邊界條件下壁面熱流密度均與風(fēng)洞實驗結(jié)果偏差較大, 不能表征真實壁面的催化復(fù)合效應(yīng), 而采用合適的有限催化模型, 則壁面熱流密度與風(fēng)洞試驗結(jié)果符合良好, 表明高超聲速飛行器壁面的催化復(fù)合效應(yīng)為有限催化.

(2)針對典型高超聲速飛行器鈍雙錐的研究表明, 在輻射平衡溫度邊界條件下駐點附近氧原子的催化復(fù)合系數(shù)約為0.17, 氮原子的催化復(fù)合系數(shù)約為0.026, 大面積區(qū)則分別降為0.005 3和0.01. 另一方面, 與完全催化和完全非催化相比, 有限催化對壁面溫度和熱流分布具有顯著影響. 駐點壁面溫度方面, 完全催化的壁面溫度峰值比有限催化高約5%, 而完全非催化的壁面溫度峰值比有限催化低約8%. 駐點熱流密度方面, 完全催化的壁面熱流峰值比有限催化高約21%, 而完全非催化的壁面熱流峰值比有限催化低約29%.

(3)高超聲速化學(xué)非平衡流場中飛行器壁面的組分復(fù)合效應(yīng)有兩個組成部分, 第一是由壁面材料催化直接導(dǎo)致的復(fù)合, 第二是由壁面附近低溫導(dǎo)致的化學(xué)反應(yīng)速率變化而間接引起的復(fù)合. 輻射平衡壁溫條件下, 與完全非催化相比, 有限催化效應(yīng)將間接改變壁面溫度梯度及壁面熱傳導(dǎo)系數(shù), 進(jìn)而影響熱傳導(dǎo)熱流.

本文建立的高超聲速飛行器化學(xué)非平衡流場輻射平衡與有限催化邊界耦合計算方法具有較大的拓展空間, 未來可結(jié)合分子動力學(xué)計算得到的壁面催化復(fù)合系數(shù), 實現(xiàn)跨尺度的CFD數(shù)值模擬.