高鐵鎖,江濤,傅楊?yuàn)W驍,丁明松,劉慶宗,董維中,許勇,李鵬

(中國(guó)空氣動(dòng)力研究與發(fā)展中心 計(jì)算空氣動(dòng)力研究所,四川 綿陽 621000)

0 引言

高超聲速技術(shù)是未來先進(jìn)飛行器發(fā)展的核心技術(shù)之一。在對(duì)高超聲速飛行器優(yōu)化設(shè)計(jì)時(shí),除了分析評(píng)估飛行器的氣動(dòng)特性和飛行性能,還需要評(píng)估飛行器周圍等離子體流動(dòng)對(duì)目標(biāo)特性和通信性能的影響[1-3]。高超聲速飛行器在大氣中飛行時(shí),與周圍氣體發(fā)生強(qiáng)烈作用產(chǎn)生強(qiáng)激波效應(yīng),使得波后氣體溫度急劇升高,氣體分子發(fā)生振動(dòng)激發(fā)、分解和電離等復(fù)雜氣動(dòng)物理化學(xué)現(xiàn)象,在飛行器周圍形成高溫電離氣體層俗稱等離子體鞘套。電磁波通過等離子鞘套傳播時(shí),被等離子體反射、折射和吸收,發(fā)生強(qiáng)度衰減、傳播方向偏折和相位畸變等效應(yīng),導(dǎo)致電磁波作用距離縮短和信噪比下降等電磁性能退化現(xiàn)象,情況嚴(yán)重時(shí),等離子體使電磁波傳輸中斷,出現(xiàn)“黑障”現(xiàn)象,這些現(xiàn)象一般統(tǒng)稱為氣動(dòng)電磁波傳輸效應(yīng)[4]。氣動(dòng)電磁波傳輸效應(yīng)特別是“黑障”現(xiàn)象嚴(yán)重影響飛行器與地面之間的微波通信性能,對(duì)飛行器實(shí)時(shí)控制和飛行安全構(gòu)成嚴(yán)重威脅[5-6]。因此,研究氣動(dòng)電磁波傳輸效應(yīng)對(duì)“黑障”現(xiàn)象的分析評(píng)估及高超聲速飛行器通信系統(tǒng)的設(shè)計(jì)都具有重要意義。

隨著臨近空間飛行器飛行速度的不斷提高,氣動(dòng)電磁波傳輸效應(yīng)對(duì)電磁通信的影響問題愈加突出。考慮到飛行試驗(yàn)?zāi)M的高成本及地面試驗(yàn)對(duì)實(shí)際飛行條件模擬能力限制等因素,對(duì)氣動(dòng)電磁波傳輸效應(yīng)問題的研究目前仍以數(shù)值模擬為主要手段。1995年,Lin等[7]總結(jié)分析了飛行器等離子體鞘套的計(jì)算和試驗(yàn)研究狀況,給出了均勻等離子體中各種電磁波傳輸效應(yīng)的解析表述;1998年,Nusca等[8]采用7組分化學(xué)非平衡N-S方程和波動(dòng)方程,計(jì)算分析了平面電磁波在無線電衰減測(cè)量C(Radio Attenuation Measurement C,RAM-C)再入等離子體鞘套中衰減效應(yīng);2002年,Funaki等[9]結(jié)合固體火箭發(fā)動(dòng)機(jī)地面試驗(yàn),預(yù)估分析了發(fā)動(dòng)機(jī)羽流等離子體對(duì)不同波段電磁波的衰減特性;2003年,Starkey等[10]采用Park化學(xué)模型和解析方法,分析了航天飛機(jī)等多種飛行器等離子體鞘套對(duì)電磁波通信的影響;2006年,White等[11]基于高階數(shù)值方法,分析了RAM-C再入等離子體鞘套的電磁波傳輸效應(yīng);2009年,Kim等[12]采用非平衡等離子體數(shù)值模擬方法及傳輸效應(yīng)解析方法,研究了高超聲速飛行器軌道再進(jìn)入試驗(yàn)(Obiter Re-entry Experiment,OREX)繞流等離子體對(duì)電磁波傳輸?shù)挠绊?探索外加磁場(chǎng)降低通信窗口等離子體傳輸效應(yīng)的方法;2018年,龔旻等[1]對(duì)臨近空間飛行器“黑障”現(xiàn)象的數(shù)值和地面模擬方法進(jìn)行了回顧總結(jié);2019年,左光等[13]采用數(shù)值及解析方法,分析了大鈍頭返回艙和類X-37B升力體飛行器再入等離子體鞘套對(duì)電磁波傳輸?shù)挠绊懱卣?。上述研究工作的共性特點(diǎn)是把流體力學(xué)和電磁學(xué)結(jié)合起來,分析高超聲速飛行器等離子體鞘套中電磁波傳輸效應(yīng),但對(duì)不同特征幾何尺度高超聲速飛行器等離子體鞘套及電磁波傳輸效應(yīng)的產(chǎn)生機(jī)制和規(guī)律分析較少。

本文基于求解非平衡流場(chǎng)N-S控制方程和求解電磁場(chǎng)波動(dòng)方程的模擬方法,采用自主研發(fā)的氣動(dòng)物理流場(chǎng)數(shù)值模擬軟件和再入黑障預(yù)測(cè)分析軟件,研究高超聲速再入體周圍等離子體的產(chǎn)生機(jī)制及氣動(dòng)電磁波傳輸效應(yīng),重點(diǎn)分析不同頻率電磁波在再入等離子體鞘套中傳播的衰減效應(yīng),研究不同球頭半徑尺寸再入體的周圍等離子體的分布規(guī)律及其對(duì)電磁波衰減的影響特征。

1 高超聲速等離子體流場(chǎng)數(shù)值模擬研究

1.1 控制方程

高超聲速飛行器從外層空間再入大氣層過程中,由于熱化學(xué)過程與流動(dòng)過程的特征時(shí)間尺度效應(yīng),飛行器繞流等離子體一般要經(jīng)歷熱力學(xué)和化學(xué)非平衡過程。文獻(xiàn)[14]的研究發(fā)現(xiàn),熱力學(xué)非平衡效應(yīng)對(duì)高超聲速飛行器非平衡繞流中等離子體分布的影響很小,此時(shí)可以不單獨(dú)考慮氣體分子的振動(dòng)能量模態(tài),而是基于單溫度氣體模型,通過求解含化學(xué)反應(yīng)源項(xiàng)的化學(xué)非平衡流動(dòng)N-S方程,數(shù)值模擬高超聲速行器周圍等離子體鞘套,三維化學(xué)非平衡流動(dòng)的N-S方程的無量綱化形式[14]如下：

(1)

式中：Q為守恒變量,Q=(ρi,ρ,ρu,ρv,ρw,ρE)T,ρi和ρ是組分的分密度和混合氣體的總密度,u、v、w為直角坐標(biāo)下3個(gè)方向的速度,E為混合氣體的總能;Re是來流雷諾數(shù);F、G、H和FV、GV、HV分別對(duì)應(yīng)3個(gè)方向的無黏和黏性通量項(xiàng);W為化學(xué)非平衡源項(xiàng),W=(wi,0,0,0,0,0)T,wi為i組分生成源項(xiàng),

(2)

(3)

1.2 化學(xué)模型

高超聲速飛行器與來流空氣發(fā)生作用形成脫體激波,波后壓縮空氣溫度急劇增高,使得來流壓縮空氣不僅發(fā)生離解反應(yīng)和置換反應(yīng),而且發(fā)生各種電離反應(yīng),在不計(jì)空氣中微量組分的情況下,認(rèn)為來流空氣由O2和N2組成,此時(shí)高溫空氣中的電離反應(yīng)主要包括締合電離和碰撞電離反應(yīng)機(jī)制[14]：

2 數(shù)值計(jì)算方法

2.1 等離子體流場(chǎng)數(shù)值模擬方法

基于有限差分方法對(duì)流動(dòng)控制方程式(1)進(jìn)行數(shù)值求解：對(duì)于方程中的無黏通量項(xiàng),采用上下對(duì)稱的高斯-賽德爾(Lower-Upper Symmetric Gauss Seidel,LU-SGS)隱式處理分法[18],以解決無黏通量雅可比矩陣直接求逆帶來的運(yùn)算量大的問題;為了解決流動(dòng)和化學(xué)過程特征時(shí)間尺度效應(yīng)導(dǎo)致的剛性問題,采用全隱式耦合方法對(duì)方程中化學(xué)反應(yīng)源項(xiàng)進(jìn)行處理[19],即用同一時(shí)間尺度同時(shí)求解流動(dòng)方程和化學(xué)反應(yīng)方程以保證計(jì)算的穩(wěn)定性和收斂性;為了準(zhǔn)確模擬流動(dòng)中激波結(jié)構(gòu)和邊界層特性,同時(shí)保證計(jì)算的魯棒性,采用高分辨壓力權(quán)函數(shù)修正的迎風(fēng)型矢通量分裂(Advection Upstream Splitting Method by Pressure-based Weight functions ,AUSMPW+)格式[20]離散無黏通量,采用中心差分格式離散項(xiàng)黏性通量項(xiàng);為了模擬高空稀薄效應(yīng)的影響,引入壁面參數(shù)滑移修正模型[21];為了模擬多組分混合氣體輸運(yùn)特性,混合氣體的黏性系數(shù)和熱傳導(dǎo)系數(shù)用Wilke半經(jīng)驗(yàn)公式計(jì)算,各組分的輸運(yùn)系數(shù)基于Blotter曲線擬合公式和Eucken關(guān)系式計(jì)算,擴(kuò)散系數(shù)采用等效二元擴(kuò)散模型計(jì)算,具體計(jì)算方法詳見文獻(xiàn)[18-19]。

2.2 等離子體中電磁波傳輸特性預(yù)測(cè)方法

考慮一維情況,即等離子體參數(shù)僅在z軸方向非均勻分布,且平面電磁波沿z軸正向傳播,此時(shí)電場(chǎng)只是隨著z軸發(fā)生變化,設(shè)電場(chǎng)平行于y軸,此時(shí)波動(dòng)方程[22-24]為

(4)

式中：Ey為電場(chǎng)強(qiáng)度;k為波數(shù)。對(duì)于緩變介質(zhì),即介質(zhì)的電磁參數(shù)在z軸方向變化較小,該方程的Wentzel-Kramers-Brillouin(WKB)解[24]為

(5)

式中：Ey0為z=0 m處的電場(chǎng)強(qiáng)度。設(shè)電磁波從z=0 m處垂直入射到等離子體內(nèi)部,并在z=d界面處透射出來,電磁波在此處的能量衰減為

(6)

在各向同性非磁化等離子體介質(zhì)中波數(shù)k為復(fù)數(shù),其表達(dá)式為

(7)

式中：ω與ωp分別表示電磁波與等離子體的角頻率;c為光速;r為等離子體介質(zhì)的相對(duì)介電常數(shù);ν為等離子體的碰撞頻率,碰撞頻率基于工程經(jīng)驗(yàn)關(guān)系進(jìn)行計(jì)算[4]。

3 數(shù)值計(jì)算結(jié)果

3.1 不同天線位置等離子體分布及其電磁傳輸特性

研究狀態(tài)為RAM-C再入飛行狀態(tài)[25]：再入體為球錐對(duì)稱體外形,再入高度H=71 km,飛行速度V=7.65 km/s,壁面催化條件為非催化壁面(Non Catalytic Wall,NCW),壁面溫度Tw=1 500 K,飛行攻角近似為0°。在零攻角飛行狀態(tài)下等離子體流動(dòng)為軸對(duì)稱流動(dòng),等離子體分布沿軸向變化是主要的,因此通信天線的軸向安裝位置將嚴(yán)重影響電磁波的傳輸效果。

圖1 不同軸向位置法向峰值與臨界等離子體分布Fig.1 Peak and critical plasma parameters distribution at different axial positions

實(shí)際上,再入體周圍等離子體鞘套處于空間(定常情況下)非均勻分布狀態(tài),繞流等離子體厚度和大小分布均在變化。為了定量分析再入體周圍等離子體鞘套對(duì)電磁波衰減的影響,以再入等離子體參數(shù)為基礎(chǔ),認(rèn)為電磁波沿物面法向傳播,采用上述WKB方法預(yù)測(cè)電磁波通過等離子體的衰減量,如圖2所示。從圖2(a)可以看出,對(duì)于C頻段和X頻段高頻電磁波,其在等離子體中的衰減沿軸向總體上變小,而低頻VHF電磁波的衰減量在球頭附近區(qū)域出現(xiàn)極小值和極大值,使得電磁波衰減沿軸向呈現(xiàn)非線性變化。盡管越靠近頭部峰值電子數(shù)密度越高,但頭部等離子體鞘套的厚度較小,且等離子體參數(shù)沿物面法向和流向非均勻分布,電磁波衰減隨等離子體參數(shù)(電子數(shù)密度和碰撞頻率)非線性變化,導(dǎo)致不同頻率電磁波衰減沿軸向變化出現(xiàn)不同變化規(guī)律。總之,不管是低頻(VHF頻段)還是高頻(C、X頻段)電磁波,其通過飛行器后部等離子體鞘套的衰減量均小于通過頭部等離子體的衰減量,而且在同一軸向位置處,電磁波頻率越高其衰減量越小,因此把天線位置安裝在靠近飛行器后身部,或適當(dāng)選擇頻率較高的電磁頻段,有利于減少飛行器等離子體鞘套對(duì)電磁通信的影響。

圖2 不同軸向位置的電磁波能量衰減和峰值碰撞頻率Fig.2 Electromagnetic energy attenuation and peak collision frequency at different axial positions

為了分析等離子體碰撞頻率的影響,圖2(a)中還給出了理想無碰撞等離子體(假定碰撞頻率ν=0 Hz)的預(yù)測(cè)結(jié)果。從圖2(a)可見,碰撞頻率對(duì) C頻段和X頻段電磁波的衰減的幾乎沒有影響,而對(duì)低頻VHF電磁波在球頭附近區(qū)域等離子體中的衰減有較大影響。圖2(b)給出了物面法向等離子體剖面的峰值碰撞頻率隨軸向的變化情況。由圖2(b)可知：C頻段和X頻段的電磁波頻率顯著高于等離子體的碰撞頻率,而VHF頻段電磁波頻率低于再入體球頭附近等離子體的碰撞頻率;碰撞頻率對(duì)電磁波衰減的影響主要與碰撞頻率與電磁波頻率的相對(duì)大小有關(guān)[22],即當(dāng)電磁波頻率接近或小于等離子體碰撞頻率時(shí),碰撞頻率對(duì)電磁波衰減才會(huì)表現(xiàn)出來,二者差異越大,等離子體碰撞頻率對(duì)電磁波衰減的影響越大。

3.2 不同尺寸飛行器等離子體分布及電磁傳輸效應(yīng)

以球頭半徑Rn=15.24 cm的RAM-C再入等離子體鞘套[25]為基礎(chǔ),改變球頭半徑尺寸,保持長(zhǎng)度不變(L=1.295 m),NCW條件,壁面溫度Tw=1 500 K,飛行攻角近似為0°。分析不同球頭半徑(Rn為15.24 cm、7.5 cm、2.5 cm)對(duì)再入體周圍等離子體分布的影響。

圖3給出不同球頭半徑飛行器等離子體分布,圖中每條等值線上的數(shù)字表示對(duì)應(yīng)的電子數(shù)密度。可見隨著球頭半徑減小,等離子體鞘套層的厚度明顯減小,鞘套層的外邊界越貼近壁面。圖4進(jìn)一步給出不同尺寸再入體頭部駐點(diǎn)線的流場(chǎng)溫度和電子數(shù)密度分布情況,圖中R為壁面法向距離。隨著球頭半徑增大,沿駐點(diǎn)線激波層的厚度也隨著增大,激波電離層電子數(shù)密度的峰值也隨之增大,如圖4(b)所示。大尺寸球頭激波層氣體電離度的增強(qiáng)意味著較強(qiáng)電離效應(yīng),而氣體離解和電離等吸熱反應(yīng)使得大尺寸球頭來流氣體過激波后沿駐點(diǎn)線溫度下降得更快,如圖4(a)所示。實(shí)際上,在同一飛行條件下,飛行器頭部尺寸越大,脫體激波面的曲率半徑越大。曲面激波可看作是多個(gè)微元段平面斜激波段的組合,曲率半徑越大,對(duì)應(yīng)每個(gè)微元面的激波角越大,頭部接近正激波面元段就越多,激波對(duì)來流氣體的壓縮效應(yīng)就越強(qiáng),波后壓縮氣體溫升使得化學(xué)電離效應(yīng)就越強(qiáng)。此時(shí)高速來流氣體通過激波滯止過程中動(dòng)能減少,一部分動(dòng)能轉(zhuǎn)化為氣體粒子的熱運(yùn)動(dòng)能,另一部則轉(zhuǎn)化為離解和電離產(chǎn)物的零點(diǎn)能。

圖3 不同球頭半徑再入體周圍等離子體鞘套特征Fig.3 Plasma sheaths around the reentry body with different sphere radii

圖4 不同球頭半徑再入體壁面法向等離子體分布(x=0 m)Fig.4 Wall-normal plasma distribution (x=0 m)

圖5給出再入體身部1.234 m處(通信天線附近)壁面法向剖面的NO+質(zhì)量分?jǐn)?shù)和電子數(shù)密度分布,比較圖4(b)和圖5(b)不難看出,從頭部駐點(diǎn)區(qū)域到尾部天線附近1.234 m處,對(duì)于不同尺寸的再入體,等離子體鞘套剖面內(nèi)平均電子數(shù)密度降低均超過兩個(gè)量級(jí),而且頭部半徑越小,等離子體中電子數(shù)密度衰減越快,這是因?yàn)檠卦偃塍w長(zhǎng)度方向等離子體特征電子數(shù)密度與球頭半徑的冪次方呈正比關(guān)系[6,22],球頭半徑越大,從頭部開始等離子體層厚度相對(duì)于上游增加的速度越快,等離子體電子數(shù)密度相對(duì)于上游衰減越慢。從圖5還可以發(fā)現(xiàn),在天線附近1.234 m處,等離子體剖面內(nèi)平均電子數(shù)密度和NO+質(zhì)量分?jǐn)?shù)隨頭部半徑的變化規(guī)律基本一致,二者均隨著頭部半徑的減小而減小,表明此時(shí)NO的電離效應(yīng)對(duì)等離子體電子數(shù)密度分布起主導(dǎo)作用,這可以從圖6中各種離子組分的分布特征得到進(jìn)一步印證。從圖6可以看出,對(duì)于不同頭部半徑尺寸的再入體,在此尾部天線附近位置處,沿物面法向的所有離子組分分布中,NO+質(zhì)量分?jǐn)?shù)最大,是最主要的電離組分,說明有關(guān)NO電離反應(yīng)機(jī)制對(duì)等離子體中電子數(shù)密度分布的貢獻(xiàn)最大。

圖5 再入體壁面法向等離子體分布(x=1.234 m)Fig.5 Wall-normal plasma distribution (x=1.234 m)

圖6 再入體壁面法向離子組分質(zhì)量分?jǐn)?shù)分布(x=1.234 m)Fig.6 Ion mass fraction along wall normal (x=1.234 m)

圖7給出球頭半徑對(duì)等離子體中電磁波傳輸特性影響,可見,隨著球頭半徑增加,由于峰值電子數(shù)密度隨之增加,等離子體對(duì)電磁波傳輸過程中產(chǎn)生的衰減量也隨之增加。實(shí)際上,隨著球頭半徑增加,等離子體的厚度也隨之增加,即球頭半徑對(duì)等離子體參數(shù)大小與厚度分布均產(chǎn)生影響,從而影響電磁波在其中的傳輸特性。從圖7中還可以看出,再入體的球頭半徑改變時(shí),等離子體對(duì)不同頻率電磁波衰減的影響的程度有所不同,但電磁波衰減量沿軸向的變化規(guī)律基本一致。

圖7 球頭半徑對(duì)電磁波衰減特性的影響Fig.7 Effect of sphere radius on energy attenuation of electromagnetic wave

3.3 OREX通信中斷預(yù)測(cè)分析

在對(duì)飛行器OREX通信中斷預(yù)測(cè)分析之前,先對(duì)本文計(jì)算模型方法及軟件的可靠性進(jìn)行驗(yàn)證分析。沿飛行彈道的高超聲速飛行器周圍等離子體及通信中斷數(shù)據(jù)可以通過飛行測(cè)量獲得,但目前公開發(fā)表飛行測(cè)量數(shù)據(jù)很少,而地面試驗(yàn)還不能完全模擬實(shí)際飛行條件,且缺乏天地?fù)Q算相似準(zhǔn)則。RAM-C飛行器再入試驗(yàn)提供了相對(duì)全面包括等離子體分布和通信中斷的測(cè)量數(shù)據(jù)[12,25],該數(shù)據(jù)是目前少有的能夠同時(shí)驗(yàn)證等離子體流場(chǎng)及其對(duì)電磁波通信影響預(yù)測(cè)方法的測(cè)量數(shù)據(jù),這里基于此測(cè)量數(shù)據(jù)及彈道條件開展驗(yàn)證分析。數(shù)值模擬等離子體流場(chǎng)采用7組分Park化學(xué)反應(yīng)模型[16],壁面采用完全催化壁面(Full Catalytic Wall,FCW),壁面溫度為1 500 K,通信天線位置在x=6.4Rn處,飛行器以零攻角再入。圖8(a)給出典型再入高度條件下剖面峰值電子數(shù)密度沿軸向分布,可見數(shù)值模擬結(jié)果與飛行測(cè)量結(jié)果符合很好。在沿彈道對(duì)RAM-C飛行器再入流場(chǎng)等離子體分布數(shù)值模擬的基礎(chǔ)上,計(jì)算了對(duì)應(yīng)飛行彈道條件下電磁波通過天線附近等離子體的能量衰減,計(jì)算中認(rèn)為電磁波沿天線壁面位置的法向傳播。一般認(rèn)為電磁波衰減30 dB以上時(shí)發(fā)生通信中斷[23],這樣就可以通過預(yù)測(cè)的衰減量判斷是否發(fā)生通信中斷,獲得沿彈道的通信中斷區(qū)間,如圖8(b)所示,其中測(cè)量數(shù)據(jù)包括VHF波段通信中斷區(qū)間,中斷時(shí)間持續(xù)約半分鐘[6],而飛行測(cè)量只獲得C與X波段通信中斷區(qū)間,因此圖中只比較了這兩個(gè)波段通信中斷的起始高度。從圖(8)可見,RAM C-II沿彈道再入過程中,VHF、C和X波段電磁波通信中斷預(yù)測(cè)與測(cè)量數(shù)據(jù)具有較好的一致性,預(yù)測(cè)的通信中斷起始或結(jié)束高度最大誤差在5 km以內(nèi)。

圖8 RAM-C飛行器再入預(yù)測(cè)結(jié)果與飛行測(cè)量對(duì)比Fig.8 Comparison between predictional and experimental data for RAM-C vehicle reentry

下面針對(duì)OREX條件[12,26]開展計(jì)算分析,其外形尺寸、等離子體探針及通信天線布置情況見圖9。在OREX飛行器以零攻角再入[12]飛行條件下,飛行器繞流對(duì)稱,探針附近等離子體分布可以反映通信天線附近的等離子體特征。圖10給出了FCW和NCW條件下電子數(shù)密度的數(shù)值模擬結(jié)果與電子探針測(cè)量結(jié)果的比較,其中橫坐標(biāo)表示探針分布方向,離散點(diǎn)表示飛行器肩部分布探針的測(cè)量結(jié)果。由圖10可以看出,數(shù)值模擬結(jié)果在總體上和飛行測(cè)量結(jié)果具有較好的一致性,高度84.01 km的符合程度更好,高度79.9 km的計(jì)算結(jié)果相對(duì)測(cè)量結(jié)果偏高,但計(jì)算誤差在3倍左右,這與之前對(duì)RAM-C飛行器再入等離子體預(yù)測(cè)精度相當(dāng)[14],表明了數(shù)值模擬等離子體模型算法及軟件的可行性。至于高度79.9 km的數(shù)值模擬結(jié)果偏高的原因,很可能是高溫空氣反應(yīng)特別是電離反應(yīng)速率常數(shù)數(shù)據(jù)的不確定性所致,特別是高溫條件下,不同化學(xué)動(dòng)力學(xué)模型對(duì)等離子體分布的數(shù)值模擬結(jié)果可存在量級(jí)上差異[27],而高度79.9 km下流場(chǎng)中氣體電離效應(yīng)相對(duì)84.01 km更強(qiáng),有關(guān)電離反應(yīng)化學(xué)動(dòng)力學(xué)數(shù)據(jù)對(duì)數(shù)值結(jié)果的影響更大。圖11給出OREX飛行器在高度59.6 km條件下繞流等離子體及頭部駐點(diǎn)線上離子組分的分布情況。從圖11(a)可見,飛行器頭部激波層內(nèi)等離子體效應(yīng)很強(qiáng),從頭部到飛行器接近尾部的肩部位置,等離子體濃度沿流向雖然有所減弱,但衰減速率比較緩慢,這與RAM-C飛行器的等離子體沿流向的變化特征存在明顯不同(參見圖1和圖3),RAM-C飛行器周圍等離子體濃度沿流向的衰減效應(yīng)相比OREX飛行器更強(qiáng),這是由于二者在飛行條件接近的情況下,OREX飛行器屬于大鈍頭體外形,其球頭半徑和錐角都明顯大于RAM-C飛行器球頭半徑和錐角,從3.2節(jié)的分析可知,飛行器頭部球頭尺寸及錐角的大小造成了飛行器周圍等離子體分布及沿流向衰減程度的不同。從圖11(b)中電離組分的分布可見,頭部激波層等離子體流場(chǎng)中,NO的電離機(jī)制是最主要的,NO+離子對(duì)等離子體濃度的貢獻(xiàn)最大,其次是O的電離效應(yīng),其他組分電離效應(yīng)較弱。圖12給出了OREX飛行器天線附近不同高度和不同催化壁條件下的壁面法向等離子體的分布情況,圖中標(biāo)識(shí)H84.0V7416表示再入高度84.0 km與對(duì)應(yīng)再入速度7 416 m/s,其他類同。從圖12不難看出,FCW和NCW的計(jì)算結(jié)果只是在靠近壁面附近有所差別,在離開壁面的大部分流場(chǎng)區(qū)域二者非常一致,沿再入彈道發(fā)展過程中,等離子體強(qiáng)度在高度 67 km 左右達(dá)到最高。圖13是以此等離子體流場(chǎng)數(shù)據(jù)為基礎(chǔ)獲得的不同頻段電磁波能量衰減沿再入彈道高度的變化情況,計(jì)算時(shí)認(rèn)為電磁波的傳播方向沿壁面安裝天線位置的法線方向,不難發(fā)現(xiàn),沿再入彈道等離子體對(duì)電磁波衰減沿再入高度降低逐漸增大并在高度67 km左右達(dá)到峰值然后開始減小。這是因?yàn)閷?duì)于特定飛行器外形和飛行姿態(tài),來流大氣密度和再入速度對(duì)飛行器周圍等離子體的分布其主導(dǎo)作用[12,22],而等離子體強(qiáng)度的沿高度分布規(guī)律決定了通信電磁波的衰減規(guī)律,在電磁波衰減出現(xiàn)峰值以前,再入速度沿高度降低減速較慢,大氣密度沿高度降低增大,對(duì)等離子體分布和衰減起主導(dǎo)作用,而在電磁波衰減出現(xiàn)峰值以后,再入飛行速度沿高度快速降低成為等離子體衰減的主導(dǎo)因素。從圖13還可以看出,FCW與NCW的衰減計(jì)算值相差不大,在認(rèn)為電磁衰減30 dB以上發(fā)生通信中斷情況下[22-23],X波段電磁波在高度84 km左右開始出現(xiàn)中斷現(xiàn)象,而L和C波段出現(xiàn)中斷的高度更高,由于再入段末端飛行減速很快,3個(gè)波段電磁波在高度50 km左右時(shí)衰減迅速降低而恢復(fù)通信,可見隨著通信電磁波頻率提高,等離子體對(duì)電磁波衰減及通信中斷的影響減弱,發(fā)生通信中斷的高度區(qū)間縮小。

圖9 OREX飛行器外形尺寸、靜電探針及天線位置Fig.9 Geometry with electrostatic probes and antenna of OREX vehicle

圖10 OREX等離子體數(shù)值模擬與飛行測(cè)量值的比較Fig.10 Comparison between computational and experimental data of OREX

圖11 OREX等離子體組分質(zhì)量分?jǐn)?shù)(H=59.6 km)Fig.11 Plasma and ion mass fraction for OREX (H=59.6 km)

圖12 天線附近壁面法向電子數(shù)密度分布Fig.12 Electron number density near antenna

圖13 OREX等離子體對(duì)電磁波衰減影響Fig.13 Electromagnetic wave attenuation due to OREX

4 結(jié)論

本文針對(duì)RAM-C和OREX鈍頭體高超聲速試驗(yàn)飛行器,開展不同特征尺度下飛行器周圍等離子體分布特性及氣動(dòng)電磁波傳輸效應(yīng)的數(shù)值模擬和分析。得出以下主要結(jié)論：

1)在同一再入條件下,對(duì)于鈍頭體飛行器,保持長(zhǎng)度不變而改變其球頭尺寸,隨著球頭尺寸的增加,球頭脫體激波內(nèi)壓縮氣體的化學(xué)效應(yīng)更強(qiáng),來流氣體動(dòng)能更多地轉(zhuǎn)化為波后氣體離解和電離的化學(xué)能,導(dǎo)致飛行器周圍等離子體鞘套厚度、等離子體的濃度以及對(duì)傳播電磁波的衰減均隨之增大,NO電離反應(yīng)對(duì)飛行器天線附近等離子體分布特性起主導(dǎo)作用。

2)再入飛行器等離子體鞘套空間分布決定了對(duì)電磁波傳輸效應(yīng)的不同特征：飛行器頭部等離子體對(duì)電磁波衰減效應(yīng)更強(qiáng),隨著流動(dòng)向下游發(fā)展等離子體效應(yīng)減弱;而在同樣等離子體分布條件下,隨著電磁波頻率的提高,等離子體對(duì)電磁波的衰減越弱;可以通過合理選擇天線安裝位置和適當(dāng)提高通信電磁波頻率的方法減緩或消除等離子體鞘套對(duì)微波通信的影響,天線安裝最好選擇飛行器的后身部位。

3)大鈍頭高超聲速飛行器(OREX飛行器)周圍等離子體鞘套厚度大且電離度強(qiáng),沿流動(dòng)方向衰減慢,在沿再入彈道飛行過程中,隨著飛行高度降低,等離子體鞘套對(duì)電磁波衰減逐漸增強(qiáng)出現(xiàn)峰值,然后隨之逐漸減弱;隨著通信電磁波工作頻率提高,等離子體鞘套對(duì)電磁波通信中斷的影響區(qū)間縮小。

4)針對(duì)自主發(fā)展的計(jì)算模型方法和軟件開展驗(yàn)證分析,典型飛行狀態(tài)下等離子體分布的數(shù)值模擬結(jié)果與飛行測(cè)量數(shù)據(jù)符合較好;典型飛行器沿彈道再入過程中,通信中斷預(yù)測(cè)與測(cè)量結(jié)果較為一致。該計(jì)算分析工具可為高超聲速飛行器的通信系統(tǒng)設(shè)計(jì)提供技術(shù)支持。