楊 鑫， 魏 兵， 尹偉科， 譚 暢

(1. 西安電子科技大學(xué) 物理與光電工程學(xué)院，陜西 西安 710071； 2. 西安電子科技大學(xué) 信息感知技術(shù)協(xié)同創(chuàng)新中心，陜西 西安 710071; 3. 陜西省等離子體物理與應(yīng)用技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710100)

當(dāng)飛行器在大氣層中以超高音速飛行時，其表面與大氣激烈的相互作用將使周圍氣體電離，形成一層包裹飛行器的等離子體薄層，即等離子體鞘套[1-2]．等離子體鞘套在一定條件下會強(qiáng)烈干擾電波信號，甚至中斷通信而威脅飛行安全，產(chǎn)生所謂的“黑障”現(xiàn)象．從20世紀(jì)60年代起，人們對鞘套中電波的傳播特性進(jìn)行了大量的研究[3-13]．但因等離子體鞘套電波傳輸特性的實(shí)驗(yàn)條件苛刻，前期往往采用理論分析和數(shù)值計算相結(jié)合的方法進(jìn)行研究．?dāng)?shù)值方法在低頻時常用時域有限差分(Finite Difference Time Domain，F(xiàn)DTD)方法，高頻時可以采用射線跟蹤算法．近年來，采用高頻率的電磁波穿透鞘套實(shí)現(xiàn)無障礙通信，已成為克服黑障效應(yīng)的可能技術(shù)之一．然而，鞘套中的電波傳輸特性與其電子濃度分布有關(guān)，而電子濃度分布又與飛行速度、高度、飛行器外形等均有關(guān)系．在現(xiàn)有文獻(xiàn)對在不同飛行器外形條件下，鞘套中電波傳輸特性的研究并不多．鑒于此，筆者采用非均勻介質(zhì)中的射線跟蹤理論，研究了典型飛行目標(biāo)(鈍頭和尖頭)外部等離子體鞘套對高頻電波的干擾特征．結(jié)果表明，尖頭外形飛行器對應(yīng)的外部鞘套流場的電子濃度量級較低，更有利于電波的傳播．

1 典型飛行器外形及其流場分布特性

圖1給出了鈍頭和尖頭飛行器在高度65 km，速度為6 800 m/s時，流場中電子濃度的分布．為方便顯示，圖1(a)和圖1(c)為取對數(shù)后的結(jié)果．由圖可見，鈍頭等離子體鞘套中的電子濃度明顯地高于尖頭; 如圖圖1(b)和圖1(d)所示，鈍頭鞘套的電子濃度在飛行器表面最大，隨徑向距離的增大而減小，隨后上升到一個極大值，然后逐步減小為零．對于尖頭鞘套，電子濃度的這種徑向分布特性并不明顯．

圖1 典型飛行器外形及其流場中電子濃度分布

2 計算原理及數(shù)值方法

當(dāng)飛行器以高頻電波與外界通訊時，可采用幾何光學(xué)方法對其傳播特性進(jìn)行近似分析．此時，電波在等離子體中的傳播過程可由如下射線方程描述[14]:

(1)

其中，X=1-n2(n為媒質(zhì)的折射率);p′表示電磁射線的積分路徑;c為真空中電波速度;r、θ和φ為球坐標(biāo)系的空間坐標(biāo)值;kr、kθ和kφ為波矢量在球坐標(biāo)系中分量．在不同方向上求解式(1)，就可獲得電波在等離子中的傳播軌跡，進(jìn)一步便可得到電波的相關(guān)傳播特性．

為了從式(1)中求解電波的傳播軌跡，還需要知道等離子體的電磁特性．根據(jù)電波在等離子體中的傳播理論[15]，傳播參數(shù)可以下列方程計算:

(2)

其中，α為衰減系數(shù)，ε′和ε″為介電常數(shù)的實(shí)部和虛部，ne、ωp、v、ω、e、me和ε0分別為電子濃度、等離子體頻率、等離子體碰撞頻率、入射波頻率、電子電量、電子質(zhì)量和真空介電常數(shù)．結(jié)合式(1)和式(2)，式(3)給出了電波能量衰減量Att的計算方法[16]:

(3)

其中，E0為電波的初始振幅，E為觀察點(diǎn)處的電場值．

計算鞘套中電磁射線軌跡時，需知任意計算點(diǎn)處介質(zhì)參數(shù)．因測試(或仿真)數(shù)據(jù)通常以三維離散方式給出，為獲得任意點(diǎn)的介質(zhì)參數(shù)值，就需要通過已知參數(shù)點(diǎn)來計算．這一過程，文中采用高維Lagrange插值方法[17]來實(shí)現(xiàn)．電磁射線方程(1)的求解采用四階Runge-Kutta法[17]．

3 結(jié)果及分析

3.1 鈍頭飛行器鞘套的傳輸特性

計算中，天線均以貼片式放置于飛行器前部 (如圖1(a)); 工作頻率取為 10 GHz、20 GHz 和 50 GHz; 飛行器的飛行高度為 65 km，速度分別為 6 120 m/s 、6 800 m/s、7 480 m/s 和 8 160 m/s; 滿足下列條件之一時，射線跟蹤終止:電波穿透鞘套;電場幅值E衰減至初始幅值E0的 1/10-5;電波頻率小于等離子體頻率．

圖2表示天線頻率為10 GHz時，不同發(fā)射角β(射線與飛行器軸線的夾角)情形對應(yīng)的射線軌跡和衰減情況(衰減強(qiáng)度是觀察點(diǎn)電場幅值比上初始電場幅值的百分比)．由圖可見，當(dāng)飛行器以較低速度飛行時，電波穿透能力強(qiáng)，能量傳輸率高，可達(dá)80%以上．但能量傳輸顯示出一定方向特性．β越接近90°，能量傳輸率越高;越靠近飛行器表面，能量衰減越強(qiáng)．由于在 65 km 高處大氣已變得稀薄，在相對低速的運(yùn)動狀態(tài)下，飛行器與大氣的相互作用較弱，鞘套的電離程度較低以至于還無法對電波傳播產(chǎn)生強(qiáng)烈的影響，所以能量傳輸效率較高；當(dāng)射線發(fā)射角β接近90°時，由于穿透鞘套所需的路徑較短，且軌跡的彎折效應(yīng)較弱，所以能量傳輸率相對較高．

圖2 不同速度和發(fā)射角度下的電磁軌跡及能量衰減

隨著速度增加，圖2顯示電波軌跡出現(xiàn)明顯的彎折．在靠近飛行器表面處甚至顯現(xiàn)出波導(dǎo)的特性，形成電波在鞘套中傳播而無法穿透鞘套的現(xiàn)象．這一現(xiàn)象可由圖1(b)中鞘套的分布特性加以解釋，圖中以點(diǎn)A和B為分界點(diǎn)分為3個區(qū)域: 飛行器表面至電子濃度極小值點(diǎn)A為內(nèi)層，點(diǎn)A至B為中間層，點(diǎn)B往外劃為外層．射線傳播過程中，若在到達(dá)分界點(diǎn)B之前就已折向鞘套內(nèi)部，如圖2(c)中20°發(fā)射角的射線，則形成波導(dǎo)現(xiàn)象; 若穿過點(diǎn)B，則其在外層的傳播過程相當(dāng)于由光疏介質(zhì)向光密介質(zhì)傳播，因而將呈現(xiàn)往垂直于飛行器軸線方向偏折的趨勢，如圖2(c)以30°發(fā)射的射線．此時，電磁射線將穿透鞘套．

圖3給出了頻率為20 GHz時典型射線軌跡，以及20 GHz和50 GHz時的電波衰減特性．由圖可見，電波傳播特征與 10 GHz 時有相似之處的，例如能量傳播的方向性及射線彎折等．但在提高頻率的情形下，射線的直線性顯著地增強(qiáng)，能量傳輸效率得到了有效的提升，比如當(dāng)頻率升高至 50 GHz 時，幾乎50%以上的能量可穿透鞘套．

圖3 不同飛行速度下電波傳播及電磁參數(shù)特性

觀察圖3中能量衰減與飛行速度的關(guān)系可見，當(dāng)速度超過閾值vt后，比如 7 480 m/s 與 8 160 m/s，能量傳輸率與飛行速度的關(guān)系將與低速時的相反．這一現(xiàn)象可由圖3(d)解釋，它給出鞘套的電子密度與飛行速度的關(guān)系，并根據(jù)曲線的特征分為2個區(qū)域．在區(qū)域1內(nèi)，隨飛行器速度的增加，飛行器表面與大氣的相互作用增強(qiáng)，電離效應(yīng)隨之升高，粒子含量亦升高，使得該區(qū)域衰減特性與飛行速度呈正相關(guān)．但當(dāng)飛行速度超過臨界速度vt后，電離效應(yīng)與飛行速度則呈現(xiàn)與低速時相反的變化關(guān)系，所以電磁衰減特性將負(fù)相關(guān)于飛行速度．在區(qū)域2內(nèi)，由于飛行器的速度已接近甚至超過第一宇宙速度，這時飛行器表面與大氣存在極強(qiáng)的相互作用．在這種極端的環(huán)境下，氣動加熱的增強(qiáng)將使周圍氣體的密度急劇降低[18-19]，從而間接導(dǎo)致離子密度的下降．同時，其他研究顯示出了類似的現(xiàn)象，例如文獻(xiàn)[20]指出，當(dāng)鈍頭飛行器的飛行速度由低速逐漸接近 7 480 m/s 時，鞘套中的電子濃度將趨于飽和．

圖4 不同飛行速度下尖頭鞘套中的電磁軌跡

3.2 尖頭飛行器鞘套的傳輸特性

由圖1(c)可見，在相同飛行條件下，由于尖頭鞘套的電子密度遠(yuǎn)低于鈍頭情形，并且兩者電子密度的空間分布形式也有顯著的差別．因此，與鈍頭鞘套中的電波傳播特性相比，尖頭鞘套將更有利于電波的傳播．為顯示這一特性，圖4給出了尖頭鞘套中電磁射線的傳播情況，計算參數(shù)設(shè)置為: 飛行高度為 65 km，電波頻率為 10 GHz．由計算結(jié)果可見，尖頭鞘套對電磁射線的影響較弱，直線特性良好．而通過對能量傳輸特性的分析，發(fā)現(xiàn)98%以上的能量均可穿透鞘層．可見，與相同飛行條件下的鈍頭飛行器相比，尖頭外形更有利于電波的傳播．

4 結(jié) 束 語

文中以鈍頭和尖頭超高音速飛行器為目標(biāo)，研究了飛行器外形對電波傳播的影響特征．研究表明，鈍頭飛行器產(chǎn)生的鞘套對電波傳播有著強(qiáng)烈的影響，顯示出豐富的傳播特性，比如強(qiáng)烈的能量衰減，嚴(yán)重時甚至出現(xiàn)“波導(dǎo)”現(xiàn)象．隨著電波頻率的升高，這些現(xiàn)象將逐漸減弱．而對于尖頭飛行器，由于其外形利于流場的通過，鞘套中的等離子體分布與鈍頭飛行器有著明顯的區(qū)別．因而，在相同的計算條件下，這些現(xiàn)象并未在尖頭鞘套中顯現(xiàn)．

綜上所述，由于飛行器的氣動外形強(qiáng)烈地影響著流場的特性及分布，所以對電磁信號的傳輸亦影響明顯．其中，尖頭氣動外形更有利于電磁波的傳輸，而增加電磁波的頻率，將有效地降低鞘套對波傳播的干擾，提升其整體的傳輸效率．這些結(jié)果可進(jìn)一步促進(jìn)人們對高超音速鞘套電磁特性的理解，并為緩解飛行器在高超音速條件下的通信困難提供理論參考．

[1] 李江挺, 郭立新, 金莎莎, 等. 等離子體鞘套中的電波傳播特性研究[J]. 電波科學(xué)學(xué)報, 2011, 26(3): 494-500.

LI Jiangting, GUO Lixin, JIN Shasha, et al. EM Wave Propagation Characteristic in Plasma Sheath[J]. Chinese Journal of Radio Science, 2011, 26(3): 494-500.

[2] KEIDAR M, KIM M, BOYD I D. Electromagnetic Reduction of Plasma Density During Atmospheric Reentry and Hypersonic Flights[J]. Journal of Spacecraft and Rockets, 2008, 45(3): 445-453.

[3] SCHARFMAN W E. The Use of Langmuir Probes to Determine the Electron Density Surrounding Re-entry Vehicles: NASA-CR-66108 [R]. Washington: NASA, 1965.

[4] WEAVER W L, BOWEN J T. Entry Trajectory, Entry Environment, and Analysis of Spacecraft Motion for the RAM C-3 Flight Experiment: NASA-TM-X-2562 [R]. Washington: NASA, 1972.

[5] PETRIN A B. On the Transmission of Microwaves through Plasma Layer[J]. IEEE Trasactions on Plasma Science, 2000, 28(3): 1000-1008.

[6] THOMA C, ROSE D V, MILLER C L, et al. Electromagnetic Wave Propagation through an Overdense Magnetized Collisional Plasma Layer[J]. Journal of Applied Physics, 2009, 106(4): 043301.

[7] KUNDRAPU M, LOVERICH J, BECKWITH K, et al. Modeling Radio Communication Blackout and Blackout Mitigation in Hypersonic Vehicles[J]. Journal of Spacecraft and Rockets, 2015, 52(3): 853-862.

[8] LIU Z W, BAO W M, LI X P, et al. Influence of Plasma Pressure Fluctuation on RF Wave Propagation[J]. Plasma Science and Technology, 2016, 18(2): 131-137.

[9] SONG W J, ZHANG H. Analysis of Electromagnetic Wave Propagation and Scattering Characteristics of Plasma Shealth via High Order ADE-ADI FDTD[J]. Journal of Electromagnetic Waves and Applications, 2016, 30(10): 1321-1333.

[10] MOUSAVI A, ESFANDIARI-KALEJAHI A, AKBARI-MOGHANJOUGHI M. Optical Properties of Nonextensive Inhomogeneous Plasma Sheath[J]. Physics of Plasmas, 2016, 23(7): 073511.

[11] STARKEY R P. Hypersonic Vehicle Telemetry Blackout Analysis[J]. Journal of Spacecraft and Rockets, 2015, 52(2): 426-438.

[12] 田媛, 韓一平, 牛化恒, 等. 弱電離高超聲速流場對太赫茲波傳播影響[J]. 紅外與激光工程, 2015, 44(12): 3690-3695.

TIAN Yuan, HAN Yiping, NIU Huaheng, et al. Effects of Weakly Ionized Hypersonic Flow on Propagation of Terahertz Wave[J]. Infrared and Laser Engineering, 2015, 44(12): 3690-3695.

[13] PODDAR S, SHARMA D. Blackout Mitigation During Space Vehicle Re-entry[J]. Optik, 2015, 126(24): 5899-5902.

[14] JONES R M, STEPHENSON J J. A Versatile Three-dimensional Ray Tracing Computer Program for Radio Waves in the Ionosphere: OT Report 75-76 [R]. Washington: US Department of Commerce, 1975.

[15] 王蕊, 郭立新, 楊閣. 等離子體涂層橢球目標(biāo)的電磁散射[J]. 西安電子科技大學(xué)學(xué)報, 2008, 35(5): 903-909.

WANG Rui, GUO Lixin, YANG Ge. Investigation on the Electromagnetic Scattering from the Spheroid Covered with the Plasma[J]. Journal of Xidian University, 2008, 35(5): 903-909.

[16] 袁忠才, 時家明. 高功率微波與等離子體相互作用理論和數(shù)值研究[J]. 物理學(xué)報, 2014, 63(9): 095202.

YUAN Zhongcai, SHI Jiaming. Theoretical and Numerical Studies on Interactions Between High-power Microwave and Plasma[J]. Acta Physica Sinica, 2014, 63(9): 095202.

[17] KINCAID D R, CHENEY E W. Numerical Analysis: Mathematics of Scientific Computing[M]. 3. Pacific Grove, California: Brooks Cole, 2001.

[18] RYBAK J P. Causes, Effects and Diagnostic Measurements of the Reentry Plasma Sheath: AD-0718428[R]. Fort Collins: Department of Electrical Engineering, Colorado State University, 1970.

[19] KUNDRAPU M, LOVERICH J, BECKWITH K, et al. Modeling Radio Communication Blackout and Blackout Mitigation in Hypersonic Vehicles[J]. Journal of Spacecraft and Rockets, 2015, 52(3): 853-862.

[20] LI J T, GUO L X. Research on Electromagnetic Scattering Characteristics of Reentry Vehicles and Blackout Forecast Model[J]. Journal of Electromagnetic Waves and Applications, 2012, 26(13): 1767-1778.