非平行分層等離子體鞘套電波傳播的復(fù)射線方法

張粲宇 葉紅霞 付海洋 何鴻飛 梁子長(zhǎng)

(1.復(fù)旦大學(xué) 電磁波信息教育部重點(diǎn)實(shí)驗(yàn)室,上海 200433;2.電磁散射重點(diǎn)實(shí)驗(yàn)室,上海 200438)

?

非平行分層等離子體鞘套電波傳播的復(fù)射線方法

張粲宇1葉紅霞1付海洋1何鴻飛2梁子長(zhǎng)2

(1.復(fù)旦大學(xué) 電磁波信息教育部重點(diǎn)實(shí)驗(yàn)室,上海 200433;2.電磁散射重點(diǎn)實(shí)驗(yàn)室,上海 200438)

基于損耗媒質(zhì)的復(fù)Snell定律提出采用復(fù)射線方法研究電磁波在非均勻且非平行分層的等離子體媒質(zhì)中的傳播特性. 該方法考慮非均勻平面電磁波的復(fù)射線(包括等幅度面的傳播射線和等相位面的傳播射線),追蹤電磁波復(fù)射線在每層媒質(zhì)中的傳播路徑以及它們透過(guò)媒質(zhì)分界面時(shí)的折射射線,同時(shí)根據(jù)復(fù)射線的傳播方向計(jì)算電磁波在每層媒質(zhì)中的傳播衰減,數(shù)值累計(jì)整個(gè)傳播過(guò)程中的傳播衰減即可獲得電磁波穿過(guò)等離子體鞘套的總衰減.由鈍錐體仿真流場(chǎng)數(shù)據(jù)簡(jiǎn)化出非平行分層幾何模型并采用復(fù)射線方法進(jìn)行計(jì)算分析. 計(jì)算結(jié)果表明:飛行器頭部至尾部的傳播衰減相差巨大且呈現(xiàn)迅速減小趨勢(shì),非平行分層非均勻等離子體媒質(zhì)存在某特定入射角,能使傳播衰減達(dá)到最小值.

復(fù)矢量法；射線追蹤；傳播衰減；等離子體鞘套

DOI 10.13443/j.cjors.2015090801

引 言

臨近空間飛行器(飛船返回艙、彈頭等)以超高聲速進(jìn)入大氣層時(shí),由于粘性和激波的作用,飛行器周?chē)鷷?huì)形成一個(gè)溫度高達(dá)幾千攝氏度的高溫區(qū),高溫區(qū)內(nèi)分子被分解和電離形成一個(gè)包覆飛行器的“等離子體鞘套”. 當(dāng)電磁波在等離子體中傳播時(shí),等離子體不僅會(huì)通過(guò)反射和折射改變電磁波的傳播方向,還會(huì)因?yàn)閹щ娏Ｗ娱g的碰撞吸收一部分能量,從而減弱雷達(dá)波強(qiáng)度并對(duì)雷達(dá)波形和極化方向產(chǎn)生影響,對(duì)地面測(cè)控產(chǎn)生極大干擾. 因此,研究電磁波在等離子體鞘套內(nèi)的傳播特性,對(duì)雷達(dá)目標(biāo)跟蹤、黑障通信等問(wèn)題具有重要的理論意義和應(yīng)用價(jià)值.

等離子體問(wèn)題的研究可以追溯到20世紀(jì)50年代,Glick首次提出簡(jiǎn)化的等離子體微觀物理模型,推導(dǎo)了非磁化等離子體的介電常數(shù)計(jì)算式[1];趙漢章等人在非均勻等離子體密度的雙指數(shù)分布假設(shè)基礎(chǔ)上,導(dǎo)出了電磁波穿過(guò)等離子體鞘套時(shí)的反射及透射系數(shù)的解析表達(dá)式[2];Gregolre等人從理論上推導(dǎo)了非磁化等離子體中電磁波傳輸?shù)幕竟?分析了電磁波在等離子體中的吸收因子和透射因子,研究非磁化等離子體包覆物的雷達(dá)散射特性[3];Laroussi等用WKB(Wentzel. Kramers. Brillouin)方法計(jì)算以任意角入射的電磁波在不均勻等離子體中的衰減,并分析了電磁波在等離子體中的衰減與等離子體碰撞頻率、電磁波頻率、等離子體最大密度、電磁波入射角的關(guān)系[4];Hu等人采用散射矩陣方法(Scattering Matrix Method,SMM)對(duì)拋物線和指數(shù)型分布的磁化等離子體的反射、透射、吸收進(jìn)行分析研究[5];Raymond等人提出用時(shí)域有限差分(Finite Difference Time Domain,FDTD)法研究等離子體的電波傳播問(wèn)題[6],以及各種對(duì)FDTD方法的改進(jìn)方法,以期提高計(jì)算效率[7-10]. Qian等人采用時(shí)域積分方程(Time-Domain Integral Equation,TDIE)結(jié)合Born近似思想計(jì)算了疊加等離子體鞘套后的雷達(dá)散射面積(Radar-Cross Section,RCS),可估計(jì)等離子體鞘套對(duì)RCS的影響[11].

然而,現(xiàn)有對(duì)等離子體鞘套內(nèi)電波傳播特性的研究方法中,解析近似的WKB方法和SMM方法往往局限于簡(jiǎn)化的平行分層模型,如均勻平面或球面分層結(jié)構(gòu). 對(duì)實(shí)際分布情況更復(fù)雜的等離子體鞘套模型難以適用. 數(shù)值方法如FDTD方法和TDIE等,雖然能適用于任意分布的等離子體媒質(zhì),但是龐大的存儲(chǔ)空間和計(jì)算時(shí)間需求導(dǎo)致其無(wú)法應(yīng)用于實(shí)際電大尺寸問(wèn)題. 本文提出采用復(fù)射線追蹤方法研究這種既非平行分層又電大尺寸的等離子體鞘套的電波傳播特性.

1 等離子體鞘套的流場(chǎng)分析

高超聲速飛行器在再入過(guò)程中,周?chē)髿獍l(fā)生強(qiáng)烈的氣動(dòng)加熱和材料燒蝕現(xiàn)象,把這種現(xiàn)象稱(chēng)為真實(shí)氣體效應(yīng). 當(dāng)繞流流場(chǎng)的電離度達(dá)到一定程度時(shí),在飛行器周?chē)鷷?huì)形成含有大量自由電子的等離子體包覆流場(chǎng),即等離子體鞘套．為了數(shù)值模擬高速飛行器周?chē)牡入x子體流場(chǎng)參數(shù),本文采用七組元(O、N、O2、N2、NO+、NO、e-)的Navier-Stokes氣動(dòng)力學(xué)方程和K-ε湍流模型模擬計(jì)算獲得等離子體參數(shù)分布[12-13]. 考慮圖1中所示的鈍錐體模型(圖中白色區(qū)域),其尺寸參數(shù)如下:鈍錐體全長(zhǎng)L=0.5 m,底面半徑R=0.15 m,側(cè)面傾角φ=7°,鈍頭為半徑r=0.1 m的半球型,攻角β=0°,即飛行方向沿軸向.

由于鈍椎體模型幾何對(duì)稱(chēng),仿真飛行方向沿軸向,仿真計(jì)算的流場(chǎng)分布具有圓周對(duì)稱(chēng)性. 本文后面關(guān)于鈍椎體周?chē)邷氐入x子體媒質(zhì)的流場(chǎng)參數(shù)和電磁參數(shù)的分布圖中,都是只顯示了過(guò)中心軸線的剖面的一半.

圖1給出不同飛行速度和高度下模擬計(jì)算的電子密度和碰撞頻率分布情況.由圖可見(jiàn):鈍錐體周?chē)植加写蠓秶母邷仉婋x媒質(zhì),即等離子體媒質(zhì). 這些媒質(zhì)呈非均勻、非平行分層分布,飛行器頭部的氣體電離和碰撞更明顯.飛行速度提高將同時(shí)增大電子密度以及碰撞頻率,而飛行高度則主要影響碰撞頻率. 此外,低空區(qū)域由于空氣密度大,空氣各成分間的碰撞頻率更高.

電磁波在非磁化等離子體內(nèi)傳播時(shí),其反射、透射和吸收等物理作用主要由等離子體的相對(duì)介電常數(shù)決定,而相對(duì)介電常數(shù)εr與等離子體密度Ne、等離子體頻率ωpe、等離子體碰撞頻率ν和入射電磁波頻率ω0有關(guān),表示為[1]

(1)

(2)

由數(shù)值模擬獲得的繞流場(chǎng)參數(shù)數(shù)據(jù),根據(jù)式(1)計(jì)算鈍錐體周?chē)入x子體媒質(zhì)的電磁參數(shù). 圖2給出了相同飛行條件下(H=50 km、v=15 Ma)三種不同入射波頻率對(duì)等離子體電磁參數(shù)的影響. 可見(jiàn):入射波頻率較低時(shí),鞘套層較厚且密,介電常數(shù)呈現(xiàn)很強(qiáng)的等離子體特性;隨著入射波頻率的提高,鞘套變得更薄,介電常數(shù)也相對(duì)接近空氣. 此外還發(fā)現(xiàn):飛行器頭部介電常數(shù)的變化往往更加劇烈,且虛部值相比尾部要大得多.

(a) 不同速度下電子密度和碰撞頻率的分布(相同高度H=70 km)

(b) 不同高度下電子密度和碰撞頻率的分布(相同速度v=10 Ma)圖1 不同飛行條件下流場(chǎng)參數(shù)的分布

圖2 不同入射波頻率下介電常數(shù)的分布

2 復(fù)射線追蹤理論

復(fù)射線追蹤理論[15]從復(fù)矢量波動(dòng)方程理論出發(fā),將損耗媒質(zhì)中非均勻電磁波波動(dòng)方程的幅度項(xiàng)與相位項(xiàng)區(qū)分考慮,研究復(fù)射線(等幅度面射線和等相位面射線)在每層媒質(zhì)中的傳播特性. 在不同媒質(zhì)的分界面處,則根據(jù)復(fù)Snell定律(Complex Snell’s Law)確定復(fù)透射射線. 下面詳細(xì)描述復(fù)射線追蹤方法,以及它在等離子體鞘套媒質(zhì)層中的應(yīng)用.

2.1 復(fù)射線追蹤理論

當(dāng)均勻平面電磁波由空氣進(jìn)入損耗媒質(zhì)時(shí),電磁波的等相面與等幅面的傳播方向不能重合,變成非均勻平面波,電磁波的傳播矢量用復(fù)矢量來(lái)描述,如電場(chǎng)的復(fù)振幅矢量可描述為

E(r)=E0exp(ik·r).

(3)

(4)

圖3 兩種損耗媒質(zhì)分界面的反射與透射

N1sin θi=N2sin θt;

(5)

K1sin ξi=K2sin ξt.

(6)

結(jié)合上面式(1)和(4)得到關(guān)于N2和K2的計(jì)算式為

N1K1sin θisin φi)2=0;

(7)

(χre2χim2-N1K1sin θisin ξi)2=0.

(8)

由式(5)和(6)得到等相面和等幅面的透射角分別為:

聚類(lèi)分析又稱(chēng)集群分析，它是按“物以類(lèi)聚”原則研究事物分類(lèi)的一種多元統(tǒng)計(jì)分析方法，它根據(jù)樣本的多指標(biāo)、多個(gè)觀察樣品數(shù)據(jù)，定量地確定樣品、指標(biāo)之間存在的相似性或親疏關(guān)系，并據(jù)此聯(lián)結(jié)這些樣品或指標(biāo)，歸成大小類(lèi)群，構(gòu)成分類(lèi)樹(shù)結(jié)構(gòu)圖，是分區(qū)最常用的方法之一。本研究采用SPSS統(tǒng)計(jì)軟件進(jìn)行分類(lèi)，在一級(jí)分區(qū)的基礎(chǔ)上，完成二級(jí)分區(qū)。

(9)

(10)

因此,等相面和等幅面?zhèn)鞑サ姆较蚴噶繛?

(11)

(12)

(13)

與式(3)比較可見(jiàn),沿該路徑傳播后的傳播損耗為

(14)

復(fù)射線追蹤過(guò)程中,等相面的追蹤路徑即電磁波經(jīng)歷不同分層媒質(zhì)的傳播路徑,即相位傳播的方向,等幅面的方向即電磁波在每層媒質(zhì)中等幅衰減的方向. 由于該方法數(shù)值追蹤電磁波在每個(gè)媒質(zhì)分界面以及每層媒質(zhì)中的電磁波復(fù)射線,能適用于任意分層形狀的多分層介質(zhì)模型,克服了傳統(tǒng)WKB方法對(duì)媒質(zhì)平行分層的限制,具有更廣泛的應(yīng)用價(jià)值,尤其適合高速飛行器周?chē)姆蔷鶆蚍瞧叫蟹謱拥牡入x子鞘套模型的電磁分析.

2.2 算法驗(yàn)證

復(fù)射線追蹤算法和WKB方法計(jì)算不同入射波頻率和不同碰撞頻率下電磁波經(jīng)過(guò)該等離子體層的雙程衰減.

復(fù)射線追蹤算法中,該等離子體媒質(zhì)模型被平均分成500個(gè)離散均勻媒質(zhì)層進(jìn)行分析,相鄰層之間的分界面處考慮電磁波的折射,實(shí)現(xiàn)電磁波的復(fù)射線追蹤和傳播衰減的累計(jì)計(jì)算. 垂直入射情況下電磁波逐層垂直進(jìn)入等離子體層,遇到目標(biāo)后再按原路徑垂直返回,電磁波兩次經(jīng)過(guò)該等離子體媒質(zhì)層的總衰減為

(15)

WKB方法只適用于平行分層的等離子模型,等離子體媒質(zhì)對(duì)電磁波傳播的雙程衰減可由如下的解析積分計(jì)算[19]:

(16)

(17)

采用以上兩種方法對(duì)斜入射(θi=60°)情況下的雙程衰減計(jì)算結(jié)果如圖4所示,圖中直線為WKB方法的計(jì)算結(jié)果,圓點(diǎn)為射線追蹤方法的計(jì)算結(jié)果,不同線條對(duì)應(yīng)不同碰撞頻率或入射波頻率的情況. 結(jié)果表明:復(fù)射線追蹤方法對(duì)非均勻平行分布等離子體的雙程衰減計(jì)算結(jié)果與WKB方法的計(jì)算結(jié)果完全一致,證明了復(fù)射線方法用于非均勻等離子體傳播計(jì)算的可效性. 由于WKB方法只能用于幾何平行分層的二維不均勻性等離子體模型,針對(duì)非平行分層的等離子體鞘套模型將采用復(fù)射線追蹤方法進(jìn)行分析和計(jì)算.

(a) 雙程衰減與入射波頻率f0的關(guān)系

(b) 雙程衰減與碰撞頻率ν的關(guān)系圖4 雙程衰減與入射波頻率f0 、碰撞頻率ν的關(guān)系

3 等離子體鞘套的傳輸特性分析

取高度H=50 km、速度v=15 Ma、入射波頻率f0=28 GHz的典型再入?yún)?shù)為例[20],其周?chē)牡入x子體鞘套折射率虛部χim的分布如圖5所示. 考慮圖中沿X方向等間隔分布的十個(gè)參考點(diǎn)位置,計(jì)算電磁波從目標(biāo)表面位置沿垂直出射方向經(jīng)過(guò)等離子體層的傳播衰減. 忽略垂直入射情況下的等離子體分層的輕微不平行分布,將不同位置的流場(chǎng)參數(shù)由飛行器表面沿表面法向向外導(dǎo)出,代入平行分層等離子體模型進(jìn)行計(jì)算,其單程傳播結(jié)果如圖6所示.

圖5 鞘套側(cè)面折射率虛部χim分布及參考點(diǎn)位置

圖6的結(jié)果表明:隨著參考點(diǎn)位置向后,傳播衰減迅速減小.在該算例下,自點(diǎn)6(X=0.25 m)起,可以得到較為理想的低傳播衰減結(jié)果.

為了考慮等離子體鞘套流場(chǎng)實(shí)際存在的不平行分布特點(diǎn)以及電磁波入射角度對(duì)傳播衰減的影響,本文隨后將等離子體鞘套離散分層.由圖1(b)可以看出:從X=0.1 m起,鞘套層等離子體的電子密度分布可近似為斜劈分層模型. 一般而言,高速鈍錐體飛行器周?chē)那侍锥碱?lèi)似于此種形狀. 按圖5中相應(yīng)參考點(diǎn)位置在流場(chǎng)分層處作切線,可將等離子體鞘套模型近似為一頂角Υ=11°的斜劈模型,如圖7所示.

圖6 鞘套不同位置的傳播衰減

圖7 簡(jiǎn)化幾何模型及復(fù)射線追蹤

本文采用復(fù)射線追蹤法對(duì)靠尾部的媒質(zhì)區(qū)域進(jìn)行非平行分層媒質(zhì)的射線追蹤和衰減計(jì)算. 斜劈等角度間隔地分為10層媒質(zhì),每層媒質(zhì)的電磁參數(shù)根據(jù)沿鈍錐體表面法向向外導(dǎo)出的流場(chǎng)參數(shù)計(jì)算,入射電磁波從鈍椎導(dǎo)體表面沿外法向出發(fā)進(jìn)入等離子體媒質(zhì)層.之后采用復(fù)射線追蹤方法追蹤電磁波在該非平行分層媒質(zhì)中的傳播路徑,并累積計(jì)算其總傳播衰減.

圖8所示的是尾部三個(gè)參考點(diǎn)(分別距頂部0.3 m,0.4 m,0.45 m)的傳播衰減隨電磁波入射角的變化.結(jié)果表明:隨著入射角度的變化,傳播衰減呈現(xiàn)先減小后增大趨勢(shì),即存在一個(gè)特定的角度使得傳播衰減達(dá)到最小值. 在該算例中損耗最小角度為85°,進(jìn)一步研究發(fā)現(xiàn)不同入射角度下的各層損耗系數(shù)與計(jì)入損耗的有效距離都會(huì)發(fā)生變化,但在85°入射時(shí)各層累加的損耗系數(shù)與損耗有效距離的乘積是最小的. 針對(duì)不同幾何模型以及介電常數(shù)分布,都存在相應(yīng)的傳播損耗最小角,采用復(fù)射線追蹤方法可以計(jì)算各類(lèi)等離子體分布下傳播衰減最小的角度,為工程上確定發(fā)射天線主瓣方向提供參考價(jià)值.

圖8 不同入射角度的傳播衰減

4 結(jié)論與討論

本文針對(duì)等離子體鞘套非均勻且非平行分層的復(fù)雜分布特點(diǎn),基于復(fù)Snell定律提出用復(fù)射線追蹤方法分析每層媒質(zhì)中電磁波復(fù)射線的傳播過(guò)程,并根據(jù)等幅面射線和等相面射線的傳播路徑數(shù)值累計(jì)計(jì)算電磁波經(jīng)過(guò)等離子體鞘套的傳播衰減.該算法可以克服傳統(tǒng)的解析WKB方法對(duì)平行分層的限制以及FDTD計(jì)算復(fù)雜度高的缺陷.

本文針對(duì)旋轉(zhuǎn)對(duì)稱(chēng)的鈍椎體再入模型,數(shù)值仿真獲取再入體周?chē)牧鲌?chǎng)數(shù)據(jù),然后結(jié)合流場(chǎng)數(shù)據(jù)的旋轉(zhuǎn)對(duì)稱(chēng)性和大尺度曲面分層特點(diǎn),利用切平面近似代替大尺度曲面,獲得了一個(gè)近似的斜劈分層模型,并且每層媒質(zhì)的電磁參數(shù)由仿真的流場(chǎng)數(shù)據(jù)直接計(jì)算.其后本文采用復(fù)射線追蹤方法對(duì)該斜面分層模型進(jìn)行電磁波傳播路徑追蹤和數(shù)值電磁衰減計(jì)算. 計(jì)算結(jié)果表明:鈍椎體尾部區(qū)域的電波傳播衰減相對(duì)于頭部明顯減小;在相同的位置,不同角度的電磁波在等離子體媒質(zhì)中會(huì)因折射經(jīng)歷不同的傳播路徑,從而引起不同的傳播衰減.

[1]GLICK H S. Interaction of electromagnetic waves with plasmas of hypersonic flows[J]. ARS journal, 1962, 32: 1359-1360.

[2] 趙漢章, 吳是靜, 董乃涵. 不均勻等離子體鞘套中電磁波的傳播[J]. 地球物理學(xué)報(bào), 1983, 26(1): 9-16.

ZHAO H Z, WU S J, DONG N H. The propagation of electromagnetic waves in an inhomogeneous plasma sheath[J]. Chinese journal geophysics, 1983, 26(1): 9-16. (in Chinese)

[3]GREGOIRE D J, SANTORU J, SCHUMACHER R W. Electromagnetic wave propagationin unmagnetized plasmas[R]. Washington D C: Air Force Office of Scientific Research, 1992.

[4]LAROUSSI M, ROTH J R. Numerical calculation of the reflection, absorption, and transmission of microwaves by a nonuniform plasma slab[J]. IEEE transactions on plasma science, 1993, 21(4): 366-372.

[5]HU B, WEI G, LAI S L. SMM analysis of reflection, absorption, and transmission from nonuniform magnetized plasma slab[J]. IEEE transactions on plasma science, 1999, 27(4): 1131-1135.

[6]RAYMOND J L, FORREST H, KARL S K. A frequency-dependent finite-difference time-domain formulation for transient propagation in plasma[J]. IEEE transactions on antennas and propagation, 1991, 39(1): 29-34.

[7]莫錦軍, 劉少斌, 袁乃昌.非均勻等離子體覆蓋目標(biāo)隱身研究[J]. 電波科學(xué)學(xué)報(bào), 2002, 17(1): 69-73.

MO J J, LIU S B, YUAN N C. On the stealth effect of non-uniform plasma covered radar targets[J].Chinese journal of radio science, 2002, 17(1): 69-73.(in Chinese)

[8]劉少斌, 張光甫, 袁乃昌. 等離子體覆蓋立方散射體目標(biāo)雷達(dá)散射截面的時(shí)域有限差分法分析[J]. 物理學(xué)報(bào), 2004, 53(8): 2633-2637.

LIU S B, ZHANG G F, YUAN N C. Finite-difference time-domain analysis on radar cross section of conducting cube scatter covered with plasmas[J]. Acta physica sinica, 2004, 53(8): 2633-2637. (in Chinese)

[9]LIU J F, XI X L, WAN G B, et al. Simulation of electromagnetic wave propagation through plasma sheath using the moving-window finite-difference time-domain method[J]. IEEE transactions on plasma science, 2011, 39(3): 852-855.

[10]LI J, GUO L X, JIAO Y C, et al. Composite scattering of a plasma-coated target above dispersive sea surface by the ADE-FDTD method[J]. IEEE geoscience and remote sensing letters, 2013, 10(1): 4-8.

[11] QIAN J W, XIA M Y. Simulation of scattering by a rotating hypersonic object with plasma sheath[C]//IEEE International Conference on Computational Electromagnetics(ICCEM). Hong Kong, February 2-5, 2015:258-260.

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

LI J T, GUO L X, JIN S S, et al. EM wave propagation characteristic in plasma sheath[J]. Chinese journal of radio science, 2011, 26(3):494-499. (in Chinese)

[13]常雨, 陳偉芳, 羅寧, 等. 基于物理光學(xué)法的再入等離子體包覆體空間散射特性分析[J]. 微波學(xué)報(bào), 2008, 24(2): 1-6.

CHANG Y, CHEN W F, LUO N, et al. Analysis of the spatial scattering characteristic for the reentry target cloaked by plasma based on the physical optics method[J]. Journal of microwaves, 2008, 24(2): 1-6. (in Chinese)

[14]常雨. 超聲速/高超聲速等離子體流場(chǎng)數(shù)值模擬及其電磁特性研究[D]. 長(zhǎng)沙: 國(guó)防科技大學(xué), 2009.

CHANG Y. The study of the supersonic/hypersonic plasma flow numerical simulation and its electromagnetic characteristics[D]. Changsha: National University of Defense Technology, 2009. (in Chinese)

[15]CHANG C Y, PETER J G, HOPCRAFT K I. Ray tracing in absorbing media[J]. Journal of quantitative spectroscopy and radiative transfer, 2005, 96: 327-341.

[16]DUPERTUIS M A, PROCTOR M, ACKLIN B. Generalization of complex Snell-Descartes and Fresnel laws[J]. Journal of the optical society America, 1994, 11(3): 1159-1166.

[17]鄭浩, 葉紅霞, 徐豐. 基于雙向射線追蹤的損耗介質(zhì)目標(biāo)散射技術(shù)[J]. 電波科學(xué)學(xué)報(bào), 2015, 30(5):896-902.

ZHENG H, YE H X, XU F. Extended bi-directional analytical ray tracing algorithm for lossy dielectric object[J]. Chinese journal of radio science, 2015, 30(5):896-902.(in Chinese)

[18]孫愛(ài)萍, 李麗瓊, 邱孝明, 等. 電磁波與非磁化等離子體的相互作用[J].核聚變與等離子體物理， 2002, 22(3): 135-138.

SUN A P, LI L Q, QIU X M, et al. Interaction of the electromagnetic waves and non-magnetized plasmas[J]. Nuclear fusion and plasma physics, 2002, 22(3): 135-138. (in Chinese)

[19]劉少斌, 袁乃昌, 莫錦軍. 等離子體隱身技術(shù)[M]. 北京: 科學(xué)出版社, 2005.

[20]劉嘉興. 飛行器測(cè)控通信工程[M]. 北京: 國(guó)防工業(yè)出版社, 2010.

[21]劉少斌, 莫錦軍, 袁乃昌. 非磁化等離子體密度與目標(biāo)雷達(dá)隱身的關(guān)系[J]. 電波科學(xué)學(xué)報(bào), 2003, 18(1): 57-61.

LIU S B, MO J J, YUAN N C. Research on the relation between the unmagnetized plasma density and the stealth of target[J]. Chinese journal of radio science,2003,18(1):57-61. (in Chinese)

張粲宇 (1992-),男,江西人,復(fù)旦大學(xué)電磁波信息科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室研究生,主要研究方向?yàn)閃KB近似和等離子體鞘套傳播等.

葉紅霞 (1976-),女,江蘇人,復(fù)旦大學(xué)信息科學(xué)與工程學(xué)院副教授.已發(fā)表論文40 余篇,出版專(zhuān)著一本.曾獲教育部自然科學(xué)一等獎(jiǎng).研究方向?yàn)橛?jì)算電磁學(xué)、粗糙面散射與遙感、散射與逆散射等.

EM wave propagation in non-parallel plasma sheath using complex ray tracing algorithm

ZHANG Canyu1YE Hongxia1FU Haiyang1HE Hongfei2LIANG Zichang2

(1.KeyLaboratoryforInformationScienceofElectromagneticWaves,Shanghai200433,China; 2.ScienceandTechnologyonElectromagneticScatteringLaboratory,Shanghai200438,China)

This paper introduces a complex ray tracing algorithm based on the complex Snell law to investigate the electromagnetic wave propagation in Reentry Plasma Sheath with non-uniform and non-parallel stratification. The algorithm takes the non-uniform electromagnetic complex rays (including propagation ray of equal-amplitude plane and equal-phase plane) into account, traces the path of complex rays in every layer and the refract rays in each interface. The attenuation is calculated along the direction of equal-amplitude plane in each layer and the total attenuation through the entire propagation path can be obtained by accumulate the attenuation in each layer. This paper further analyzes non-parallel stratified plasma sheath model with the complex ray tracing algorithm. Numerical results show that the attenuation sharp decreases from the top of aircraft to the tail. There exists a specific incident angle with minimum attenuation in non-parallel stratification layer.

complex vector; ray tracing; propagation attenuation; plasma sheath

10.13443/j.cjors.2015090801

2015-09-08

國(guó)家自然科學(xué)基金(No.61001007, 41404122)； 教育部科技支撐基金

TN011；P354

A

1005-0388(2016)04-0625-07

張粲宇, 葉紅霞, 付海洋, 等. 非平行分層等離子體鞘套電波傳播的復(fù)射線方法[J]. 電波科學(xué)學(xué)報(bào),2016,31(4):625-631.

ZHANG C Y, YE H X, FU H Y, et al. EM wave propagation in non-parallel plasma sheath using complex ray tracing algorithm[J]. Chinese journal of radio science,2016,31(4):625-631. (in Chinese).DOI: 10.13443/j.cjors.2015090801

聯(lián)系人： 葉紅霞 E-mail：yehongxia@fudan.edu.cn