苗文博，程曉麗，艾邦成

(中國航天空氣動力技術(shù)研究院，北京 100074)

0 引言

未來高超聲速飛行器普遍具有高空以及高飛行馬赫數(shù)的特點，此時流動的熱非平衡特性也非常突出。非平衡效應會對流動產(chǎn)生較大的影響，它會改變激波后和壁面附近的溫度以及組分分布。高馬赫數(shù)狀態(tài)時，激波后溫度可達6000K～10000K，甚至10000K以上，氧分子、氮分子開始離解并伴隨置換反應和電離反應進行，流場中原子和離子在到達壁面時會發(fā)生表面催化復合反應釋放熱量，這部分熱量是由組分濃度梯度決定的，一般也稱之為熱流的組分擴散項。當流場中原子在壁面處的復合比例比較高時，壁面的催化復合放熱不可忽略，稱之為壁面催化效應。Sebo D.E.［1］發(fā)現(xiàn)多級再入軌道器要比普通再入受到更高的熱流加熱，認為壁面催化條件是產(chǎn)生這種熱流差異的重要因素。Papadopoulos P.［2］，N.E.Afonina［3］通過數(shù)值模擬火星再入飛行發(fā)現(xiàn)不同壁面催化條件熱流相差1倍還多。

壁面催化的概念很早就為人們所關(guān)注，Goulard R.J.［4］首先對比了幾種典型金屬和氧化物材料的壁面催化復合系數(shù)并通過試驗給出了這幾種材料對于氧原子的復合催化系數(shù)的大致數(shù)值。Inger G.R.［5－6］發(fā)現(xiàn)除了受材料的表面物理化學特性影響外，壁面溫度是影響壁面催化能力的一個重要因素，并給出了關(guān)于某碳基材料的壁面催化復合系數(shù)隨溫度的擬合關(guān)系式。一般認為原子在壁面上主要發(fā)生以下兩種類型的復合過程［7］:1)自由原子撞上吸附原子復合成氣體分子(Eley-Rideal復合模型);2)壁面上吸附原子相遇復合成氣體分子(Langmuir-Hinshelwood復合模型)。

在數(shù)值模擬研究中，通常建立以下三種壁面催化條件來模擬原子在壁面處的復合，分別是完全催化壁、完全非催化壁和有限催化壁，并使用壁面分子復合系數(shù)來表征催化能力的強弱。這一參數(shù)只能通過材料表面復合試驗來獲得，存在較高的實現(xiàn)成本和難度。因此，簡化的壁面催化模型——完全催化壁和非催化壁是最常用的壁面催化條件。一般的存在以下共識，壁面完全催化壁所得熱流要大于有限催化壁所得熱流，有限催化壁所得熱流大于非催化壁所得熱流。在實際應用中，這種不同壁面催化模型之間的差異受諸多因素影響，需要選取適當?shù)膮?shù)來表征這種差異。基于這一認識，本文通過分析典型來流狀態(tài)下完全催化壁與非催化壁的熱流差異來研究壁面催化條件的影響效應。并通過對比不同高度和馬赫數(shù)的飛行狀態(tài)，研究了來流條件與熱流組分擴散項的相關(guān)性。

化學動力學模型和多溫度模型的建立是非平衡流研究的一個重要內(nèi)容，Park［8－9］、Dunn ＆ Kang［10］基于理論和試驗研究提出了各自的化學動力學模型，能夠準確地描述高溫條件下的化學反應過程。7組分模型較多應用于近空間飛行條件，11組分模型或者23組分模型主要應用于再入飛行器研究［11］。本文選取了Park的7組分反應模型，利用精度較高的 AUSM+-up格式［12］求解了簡單球錐外形的非平衡流場。來流狀態(tài)選取RAM-CII［13］飛行試驗狀態(tài)，外形為頭部半徑為0.1524m，半錐角9°的球錐外形。

1 控制方程

以三維情形為例，在所研究的控制體積內(nèi)，將含化學反應源項的基本方程以積分形式表示如下:

此處Q是基本守恒量，六面體控制體△V的體積為Ω。假設(shè)控制體中心標為(I，J，K)，左面中心標為(i，J，K)，右面中心標為(i+1，J，K)，以此類推。其中大寫字母I表征單元體中心信息，小寫字母i表征控制體面元中心信息。積分方程離散如下:

E，F(xiàn)，G 為原始通量表達式，σi+1，J，K為所在面面積，(nx，ny，nz)為該面法向量。

2 數(shù)值格式

AUSM格式的主要思想是認為流場在傳播中存在對流影響與聲波影響，為了分別考慮兩個過程，將無粘項分為對流項和壓力項進行處理，基于馬赫數(shù)對兩者分別進行特征分裂，是適用于化學反應流求解的一類高精度格式。AUSM+-up是該類格式的最新改進型，它在構(gòu)造界面馬赫數(shù)時引入壓力耗散機制，在構(gòu)造界面壓力時引入速度耗散機制，在低速流動區(qū)域抑制數(shù)值震蕩。下面給出其通量表達式。僅以項為例說明。

其中，ψ =(ci，u，v，w，H，eV)，g=(ci，nx，ny，nz，0，0)，界面馬赫數(shù)以及壓力定義見文獻［12］。

3 定解條件

本文使用了Park的有限速率7組元6反應模型，各組分為 N，O，N2，O2，NO，NO+，e。選用Park的雙溫度模型，使用平動溫度T表征平動以及轉(zhuǎn)動能，使用振動溫度Tv表征分子振動能以及粒子電子激發(fā)能?；瘜W反應特征溫度按照C.Park的設(shè)置給定［14］。

本文對三種壁面催化條件進行了數(shù)值模擬。根據(jù)Fick定律，有限催化壁條件滿足公式(kiρCi)w=［ρDi是各組分催化復合速率，α是各組分的壁面復合系數(shù)，定義為壁面復合原子數(shù)與到達壁面原子數(shù)之比。對于有限催化壁，本文僅考慮了O和N在壁面處的復合，并且認為O和N具有相同的壁面復合系數(shù)α=0.03。

4 算例驗證

為驗證算法可靠性，選取J.Muylaret(1992)文獻［15］引用的飛行試驗數(shù)據(jù)為驗證對象。幾何外形為頭部半徑175mm，半錐角4.6°長約2m的尖錐。來流馬赫數(shù)13，來流溫度T=265K，來流壓力P=53Pa，壁面溫度取飛行試驗駐點溫度Tw=343K。進行了完全催化壁、非催化壁的熱流與飛行試驗熱流的對比。圖1給出了本文完全催化壁、非催化壁熱流結(jié)果與飛行試驗數(shù)據(jù)的對比，在駐點高熱流區(qū)完全催化壁熱流比非催化壁熱流大30%左右，后部低熱流區(qū)完全催化壁熱流比非催化壁熱流大30%－40%，飛行試驗數(shù)據(jù)基本介于兩者之間。證明本程序具有較好的非平衡數(shù)值模擬能力，可以滿足熱環(huán)境模擬要求。

圖1 熱流對比Fig.1 Heat-flux along the axial direction

5 結(jié)果及分析

本文選擇RAM-CII的外形為研究對象，研究了不同壁面催化條件時的熱流預測差異，分析了不同壁面催化條件對熱環(huán)境預測的影響，研究了來流參數(shù)對熱流組分擴散項的影響。計算模型為球頭半徑r=0.1524m，半錐角9°，后部長度0.3m的球錐。

首先針對 M=25.9，H=71km，Tw=1500K 的典型飛行狀態(tài)對比分析了三種壁面催化條件對熱流預測的影響。此時的自由流速度約為7650km/s，流動的非平衡特性非常強。圖2給出了三種壁面催化條件下的表面熱流對比以及熱流的組分擴散項對比和熱傳導項對比，由圖可以看出不同壁面催化條件之間的差異主要體現(xiàn)在熱流的組分擴散項部分，不同壁面催化條件時熱流的熱傳導項貢獻相差不大，在2%之內(nèi)。由此可知，在固定壁溫的前提下，原子的壁面催化復合過程對流場的影響非常微弱，主要影響壁面熱流分布。

圖2 不同壁面催化條件下熱流對比Fig.2 Heat-flux of different catalytic conditions

因此，可以通過評估完全催化壁條件下組分擴散項熱流部分所占比例來初步估計真實熱流值的大小。從而避免需要通過試驗獲得真實條件下壁面催化復合系數(shù)的困難。圖3給出了同樣高度條件H=71km，馬赫數(shù)分別為15、20、25.9時的熱流組分擴散項與總熱流的對比?？梢钥闯鲈诋斍案叨认拢S著馬赫數(shù)增大，熱流組分擴散項在總熱流中所占比例越來越高。圖4給出了馬赫數(shù)M=20，高度分別為71、61、51km時的熱流組分擴散項與總熱流的對比。可以看出當馬赫數(shù)不變，隨著高度增高熱流組分擴散項在總熱流中所占比例也越來越高。表1給出了馬赫數(shù)以及高度參數(shù)對于熱流組分擴散項比例的影響關(guān)系。對于高超聲速流動來說，馬赫數(shù)越高激波后溫度就越高，分子振動離解就越容易，氣體化學反應進行的越劇烈，同時隨著高度增高，空氣密度降低，這又減弱了氣體分子通過碰撞完成能量交換的可能性，化學反應相應變?nèi)酢?/p>

表1 馬赫數(shù)和高度對熱流組分擴散項的影響(單位:kW/m2)Table 1 Effects of height and Mach number on mass diffusion heat-flux(unit:kW/m2)

6 結(jié)論

針對高超聲速熱化學非平衡熱環(huán)境預測問題，以RAM-CII飛行試驗模型為參考，本文研究了壁面催化條件對熱環(huán)境預測的影響以及來流條件與熱流組分擴散項的關(guān)系。研究發(fā)現(xiàn):完全催化壁時組分擴散項熱流所占比率表征了壁面催化效應的強弱;在同一高度下，隨著馬赫數(shù)升高熱流組分擴散項越來越顯著，而馬赫數(shù)相同時，隨著高度增高熱流組分擴散項所占比例越來越小，主要原因在于當?shù)鼗瘜W反應進行程度不同。

［1］SEBO D E.Heating effects of multiple skip reentry trajectories［R］.The Aerospace Corporation，Report No.ATR －9 I(6822)－ l，El Segundo，California，May 1991.

［2］PAPADOPOULOS P.Aero-thermal heating simulations with surface catalysis for the Mars 2001 aerospace mission［R］.AIAA 97－0473.

［3］AFONINA N E，GROMOV V G ＆ KOVALEV V L.Catalysis modeling for thermal protection systems of vehicles entering into martian atmosphere［R］.AIAA01 －2832.

［4］GOULARD R J.On catalytic recombination rates in hypersonic stagnation on heat transfer［J］.Jet Propulsion，1958，28(11):737－745.

［5］INGER G R.Correlation of surface temperature effect on nonequilibrium heat transfer［J］.ARS Jour.1962，32:1743－44.

［6］INGER G R.Non-equilibrium hypersonic stagnation flow with arbitrary surface catalycity including low Reynolds number effects［J］.Int.Jour.of Heat and Mass Transfer，1966，9:755－772.

［7］KUROTAKI T.Construction of catalytic model on SiO2-based surface and application to real trajectory［A］.34th AIAA Thermophysics conference［C］.19 －22 June 2000/Denver CO，2000－2366.

［8］PARK，CHUH.Assessment of two-temperature kinetic model for ionizing air［R］.AIAA －87－1574，June 1987.

［9］PARK CHUL.Two-temperature interpretation of dissociation rate data for N2and O2［R］.AIAA －88－0458，Jan.1988.

［10］DUNN，MICHAEL G and KANG SANG-WOOK.Theoretical and experimental studies of reentry plasmas［R］.NASACR－2232，1973.

［11］ROOP N GUPTA，JERROLD M YOS.A review of reaction rates and thermodynamic and transport properties for the 11-species air model for chemical and thermal non-equilibrium calculations to 30000K［R］.NASA-TM-101528，F(xiàn)ebruary 1989.

［12］MENG-SING LIOU.A further development of the AUSM+scheme towards robust and accurate solutions for all speeds［R］.AIAA Paper 2003－4116.

［13］CANDLER G V，MACCORMACK R W.The computation of hypersonic ionized flows in chemical and thermal non-equilibrium［R］.AIAA 1988－0511.

［14］PETER A GNOFFO，ROOP N GUPTA ＆ JUDY L SHINN.Conservation equations and physical models for hypersonic air flows in thermal and chemical equilibrium［R］.NASA TP－2867.

［15］MUYLARET J，WALPOT L，HAUSER J.Standard model testing in the European high enthalpy facility F4 and extrapolation to flight［R］.AIAA 92 －3905.