李郝 楊鑫2)3) 張正平2)3)?

1) (貴州大學(xué)大數(shù)據(jù)與信息工程學(xué)院, 貴陽(yáng) 550025)

2) (半導(dǎo)體功率器件可靠性教育部工程研究中心, 貴陽(yáng) 550025)

3) (貴州省微納電子與軟件技術(shù)重點(diǎn)實(shí)驗(yàn)室, 貴陽(yáng) 550025)

為利用太赫茲波解決飛行器再入過(guò)程遇到的“黑障”問(wèn)題, 以散射矩陣方法為基礎(chǔ), 分別以非均勻磁化等離子體的磁化方向、電子密度、外加磁場(chǎng)強(qiáng)度和碰撞頻率為變量, 研究了垂直入射情形下它們對(duì)太赫茲波傳輸行為的影響.結(jié)果表明: 這些參數(shù)對(duì)太赫茲波傳輸性能影響明顯, 例如按某一方向改變磁化角度對(duì)左極化和右極化太赫茲波的傳輸功率有相反的影響; 降低磁化強(qiáng)度能一定程度地避開(kāi)等離子體對(duì)右極化波的吸收; 而降低碰撞頻率能縮小等離子體對(duì)右極化波的吸收頻帶.通過(guò)調(diào)整這些參數(shù), 有望在一定程度上緩解黑障現(xiàn)象.

1 引 言

高超音速飛行器再入大氣過(guò)程中會(huì)與空氣發(fā)生摩擦, 在產(chǎn)生大量熱量的同時(shí)飛行器表面的空氣會(huì)發(fā)生電離, 生成一層薄等離子體鞘套[1].飛行器的天線窗口被等離子體覆蓋后, 其發(fā)射和接受信號(hào)的能力大大下降, 嚴(yán)重的時(shí)候通信會(huì)完全中斷, 產(chǎn)生“黑障”現(xiàn)象[2].

為了解決“黑障”問(wèn)題, 大量的科研人員對(duì)等離子體特性[3]和電磁波在等離子體中的傳播行為進(jìn)行了研究.研究人員發(fā)現(xiàn): 電磁波頻率低于等離子體頻率是導(dǎo)致“黑障”問(wèn)題的重要原因之一, 因而就連高頻段的毫米波也無(wú)法完全穿透飛行器表面形成的高電子密度等離子體鞘套.近年來(lái)THz 波技術(shù)不斷進(jìn)步, 其優(yōu)秀的穿透性和抗干擾性使得“黑障”問(wèn)題有了解決的可能, 因此, 研究THz 波在等離子體鞘套中的傳播有著重要意義.另外, 研究人員開(kāi)發(fā)了多種算法來(lái)分析電磁波與等離子體的相互作用, 如時(shí)域有限差分法(FDTD)[4?7]、解析法[8]、溫澤爾-克萊默-布里淵法(WKB)[9,10]以及散射矩陣法(SMM)[11,12]等.本文采用SMM 方法進(jìn)行研究分析, SMM 方法較好地運(yùn)用了分層思想來(lái)實(shí)現(xiàn)對(duì)等離子體不均勻性的體現(xiàn), 能直接高效地計(jì)算出總的反射、透射和吸收率, 相較于FDTD 方法也不需要進(jìn)行大量的網(wǎng)格化操作, 所以對(duì)于飛行器表面等離子體鞘套的分析, SMM 方法具有較高的效率和較低的算法復(fù)雜度.正因如此, Chen 等[13]用SMM方法研究了時(shí)變等離子體鞘層對(duì)雷達(dá)回波信號(hào)的影響, Guo 和Guo[14]用SMM 方法分析了THz 波在非均勻磁化熱等離子體中的傳輸特性, Zhang 等[15]用SMM 方法研究了THz 波在右極化并且外加磁場(chǎng)角度與入射角度平行時(shí), 在非均勻磁化等離子體中的傳輸特性.我們發(fā)現(xiàn), 這些研究大多只是基于固定的外加磁場(chǎng)角度對(duì)傳輸特性進(jìn)行研究, 沒(méi)有考慮外加磁場(chǎng)方向的變化對(duì)THz 波傳輸?shù)挠绊?一些論文也只是基于單一極化方向的分析, 然而不同的極化方向也是影響THz 波傳輸性能的因素.所以, 為了更全面地探究磁化角度的改變對(duì)左極化和右極化THz 波在非均勻磁化等離子體中傳輸?shù)挠绊? 本文加入磁化角度作為分析傳輸特性的變量之一, 使用SMM 方法分別研究了左極化和右極化THz 波與等離子體之間不同的作用特點(diǎn).

2 SMM 理論

對(duì)于分布不均勻的等離子體可以設(shè)定為圖1所示的分層模型, 將等離子體平板平均分成n層等離子體薄片, 每一薄片中的電子密度都視為均勻分布, 總的等離子體層電子密度按給定的密度函數(shù)分布.電磁波垂直入射等離子板表面, 入射區(qū)(0)和透射區(qū)(p)都為自由空間.

圖1 太赫茲波在分層等離子體中傳播的模型圖Fig.1.Model diagram of THz waves propagation in the layered plasma.

在此模型中每一層薄片等離子體中都有其傳播常數(shù)k(m):

其中, 正負(fù)號(hào)代表左極化或者右極化,θ是波傳播方向與磁化方向的夾角,X=1?iv/ω,v是碰撞角頻率,Y=ωce2/ω2,ωce是電子回旋角頻率,Z=是 第m 層的等離子體角頻率并由(3)式表示[20]:

其中Ne為自由電子密度,e和me分別為電子電荷量和電子質(zhì)量,ε0為真空中的介電常數(shù).

將圖1 中入射區(qū)域(0)和薄層(1)的電磁場(chǎng),用入射(透射)電場(chǎng)、反射電場(chǎng)、入射(透射)磁場(chǎng)和反射磁場(chǎng)表示, 并分別寫(xiě)為以下形式:

其中A是射入到第一個(gè)等離子層表面的反射系數(shù),B1和C1為電磁波從區(qū)域(0)射入等離子體薄層(1)的透射系數(shù)和薄層(2)反射回薄層(1)的反射系數(shù),k(0)是自由空間的波數(shù).這樣可以確定入射區(qū)域(0)和薄層(1)的總電磁場(chǎng)關(guān)系, 通過(guò)匹配x=0處分界面的邊界條件能得到矩陣關(guān)系:

其中

和上面相同, 繼續(xù)匹配第2 到n個(gè)分界面能得到矩陣關(guān)系:

其中

在最后一個(gè)分界面的左右是薄層(n)和自由空間(p), 在區(qū)域(p)中只有透射電磁波, 沒(méi)有反射電磁波, 所以匹配x=d處的邊界條件能得到矩陣關(guān)系:

其中

通過(guò)上面的(6)式、(8)式和(10)式可以遞歸變換得到:

其中

接著將Sg表示為Sg=(Sg1,Sg2) , 其中Sg1和Sg2表示Sg的第一列和第二列向量, 所以(12)式可以變換為

最后解得系數(shù)A和D的矩陣:

在接下來(lái)的分析中將電磁波功率歸一化處理,歸一化反射功率表示為|A|2, 歸一化透射功率表示為|D|2, 那么歸一化吸收功率為 1?|A|2?|D|2.

3 結(jié)果與分析

根據(jù)NASA 鈍頭RAM 飛行器的公開(kāi)數(shù)據(jù)[16],在不同高度下有不同的電子密度分布, 并且通常來(lái)說(shuō)可以用指數(shù)分布、高斯分布、拋物分布等模型來(lái)模擬等離子體中的電子密度分布[21].但單一的分布模型不能模擬實(shí)際分布的非對(duì)稱(chēng)性, 所以本文采用雙高斯分布模型來(lái)模擬計(jì)算以更好地貼近實(shí)際分布情況, 分布函數(shù)表示如(16)式, 函數(shù)圖像如圖2 所示.

圖2 電子密度分布與等離子體厚度的關(guān)系Fig.2.Electron density profile versus plasma thickness.

3.1 外加磁場(chǎng)方向?qū)ψ髽O化和右極化THz波傳輸功率的影響

根據(jù)現(xiàn)有飛行器表面形成的等離子體鞘套的公開(kāi)數(shù)據(jù)和流場(chǎng)計(jì)算[22], 在此設(shè)置等離子體平板的厚度為0.1 m, 為保證計(jì)算收斂性將其分成50 層薄片, 碰撞頻率fv=0.1 THz , 最大電子密度N0=1×1018m?3, 外加磁場(chǎng)強(qiáng)度B=5 T.根據(jù)(2)式, 可以對(duì)外加磁場(chǎng)的角度進(jìn)行分析.

圖3 THz 波歸一化傳輸功率與外加磁場(chǎng)角度的關(guān)系 (a), (b)右極化; (c), (d)左極化Fig.3.Relationship between normalized transmission power of THz wave and external magnetic field angle: (a), (b) Right-handed polarization; (c), (d) left-handed polarization.

由圖3 可以看出, 隨著磁場(chǎng)角度由0°增加到90°時(shí), 右極化THz 波傳輸功率基本逐漸升高, 同時(shí)左極化THz 波的傳輸功率基本逐漸降低.由此可以看出, 當(dāng)磁場(chǎng)方向與波的傳播方向垂直時(shí), 最有利于右極化太赫茲波信號(hào)的接收; 當(dāng)磁場(chǎng)方向與波的傳播方向平行時(shí), 最有利于左極化太赫茲波信號(hào)的接收.另外, 當(dāng)電磁波在外加靜態(tài)磁場(chǎng)的等離子體中傳播時(shí), 等離子體中的電子將受到磁場(chǎng)作用做回旋運(yùn)動(dòng), 而右極化波會(huì)與回旋運(yùn)動(dòng)的電子同步旋轉(zhuǎn)形成回旋共振, 導(dǎo)致能量吸收.所以在相同磁化角度情況下, 左極化波和右極化波的傳輸功率曲線產(chǎn)生了不同的變化趨勢(shì).

那么以磁化角度分別為90°和0°來(lái)分析右極化太赫茲波和左極化太赫茲波在非均勻磁化等離子體中的傳輸特性.當(dāng)磁化角度θ=90°且為右極化時(shí)相對(duì)介電常數(shù)如(17)式; 當(dāng)磁化角度θ=0°且為左極化時(shí)相對(duì)介電常數(shù)如(18)式.

在這兩種條件下就可以進(jìn)一步分析在不同電子密度、碰撞頻率和外加磁場(chǎng)強(qiáng)度條件下右極化和左極化THz 波在等離子體中的傳輸特性.

3.2 電子密度對(duì)右極化和左極化THz波反射、傳輸、吸收功率的影響

圖4(a)—圖4(c)分別給出了右極化THz 波在不同最大電子密度下的反射、透射和吸收功率.其中等離子體平板的總厚度d=10 cm , 碰撞頻率fv=0.1 THz , 外加磁場(chǎng)強(qiáng)度為B=5 T , 最大電子密度的值分別為N0=1×1016m?3,N0=1×1017m?3,N0=1×1018m?3和N0=1×1019m?3.可以看出, 較高電子密度的等離子體對(duì)太赫茲波的反射能力非常弱.隨著電子密度的增加, THz 波傳輸功率相應(yīng)降低, 且吸收功率增加, 在電子密度為N0=1×1019m?3時(shí)出現(xiàn)了傳輸阻帶, 阻帶內(nèi)幾乎所有的電磁波功率都被等離子體層吸收.另外, 當(dāng)電子密度N0<1×1019m?3時(shí), 右極化電磁波的吸收功率在0— 0.6 THz 內(nèi)先增加后減少, 在0.1—0.2 THz間出現(xiàn)一個(gè)突起峰值, 此時(shí)外加磁場(chǎng)條件下的電子回旋頻率fce=0.14 THz , 正好在峰值區(qū)域, 所以可以推斷由于右旋波的偏振方向與電子回旋方向相同, 引起了電子回旋共振, 導(dǎo)致能量吸收增強(qiáng), 并形成了突起峰值.而電子密度的增加也引起了吸收峰值的升高, 這是因?yàn)殡娮用芏鹊脑黾右馕吨入x子體中電子數(shù)量的增加, 從而導(dǎo)致電子從電磁波中吸收的能量增加.圖4(d)—圖4(f) 分別給出了左極化THz 波在不同最大電子密度下的反射、透射和吸收功率, 參數(shù)與右極化THz 波保持一致.可以看到的是, 由于左旋波的偏振方向與電子回旋方向不同, 左極化THz 波相對(duì)于右極化THz波最大的差異在于沒(méi)有了電子回旋共振的影響, 并且在相同的電子密度下左極化THz 波對(duì)等離子體有更好的穿透性能.

圖4 THz 波歸一化反射、傳輸、吸收功率與電子密度最大值的關(guān)系 (a), (b), (c)右極化; (d), (e), (f)左極化Fig.4.Relationship between normalized reflection, transmission, and absorption power of the THz wave and the maximum value of electron density: (a), (b), (c) Right-handed polarization; (d), (e), (f) left-handed polarization.

由于THz 波在等離子體中的反射功率非常低,傳輸功率和吸收功率的變化規(guī)律相反, 但傳輸功率在一定程度上能更好地體現(xiàn)信號(hào)傳輸?shù)哪芰? 所以接下來(lái)主要對(duì)傳輸功率進(jìn)行分析.

3.3 外加磁場(chǎng)強(qiáng)度對(duì)右極化和左極化THz 波傳輸功率的影響

圖5 THz 波歸一化傳輸功率與外加磁場(chǎng)強(qiáng)度的關(guān)系 (a), (b)右極化; (c), (d)左極化Fig.5.Relationship between normalized transmission power of THz wave and applied magnetic field intensity: (a), (b) Righthanded polarization; (c), (d) left-handed polarization.

圖6 THz 波歸一化傳輸功率與碰撞頻率的關(guān)系 (a), (b)右極化; (c), (d)左極化Fig.6.Relationship between normalized transmission power of THz wave and collision frequency: (a), (b) Right-handed polarization; (c), (d) left-handed polarization.

圖5 (a)和圖5(b)以及圖5(c)和圖5(d)分別給出了外加磁場(chǎng)強(qiáng)度對(duì)右極化和左極化THz 波傳輸功率的影響, 此時(shí)等離子體平板的總厚度d=10 cm , 碰撞頻率fv=0.1 THz , 最大電子密度的值N0=1×1018m?3.從圖5(a)和圖5(c)可以看出, 左極化波傳輸功率在整個(gè)頻帶上隨外加磁場(chǎng)強(qiáng)度的增加而增加.而右極化波傳輸功率隨著外加磁場(chǎng)強(qiáng)度的增加在低頻端和高頻端有著不同的變化趨勢(shì): 在低頻端右極化波傳輸功率隨外加磁場(chǎng)強(qiáng)度的增加而增加, 但由于磁場(chǎng)強(qiáng)度的增加引起回旋共振的頻帶向高頻端偏移, 進(jìn)而在高頻端產(chǎn)生更多的能量吸收, 所以在高頻端右極化波傳輸功率會(huì)隨外加磁場(chǎng)強(qiáng)度的增加而降低.

3.4 碰撞頻率對(duì)右極化和左極化THz 波傳輸功率的影響

圖6(a)和圖6(b)以及圖6(c)和圖6(d)分別給出了碰撞頻率對(duì)右極化和左極化THz 波傳輸功率的影響, 此時(shí)等離子體平板的總厚度d=10 cm ,外加磁場(chǎng)強(qiáng)度為B=5T , 最大電子密度的值N0=1×1018m?3.可以看出, 右極化波傳輸功率的最小值隨碰撞頻率的增加而減小, 但吸收頻帶增大.左極化波傳輸功率隨著碰撞頻率的增大而減少.總體來(lái)說(shuō)左極化波更能穿透磁化冷等離子體, 有助于更好地改善飛行器再入“黑障”問(wèn)題.

4 結(jié) 論

本文以NASA 鈍頭RAM 實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ),研究了左極化和右極化太赫茲波在非均勻磁化冷等離子體中的傳輸特性.結(jié)果表明, 外加磁場(chǎng)的角度對(duì)左右極化太赫茲波的影響有所不同, 在磁場(chǎng)方向與波傳輸方向垂直時(shí)最有利于右極化太赫茲波的傳輸; 在磁場(chǎng)方向與波傳播方向平行時(shí)最有利于左極化太赫茲波的傳輸.另外, 隨著等離子體最大電子密度的增加, 等離子體對(duì)太赫茲波的能量吸收增強(qiáng), 電磁波在右極化時(shí)會(huì)發(fā)生回旋共振, 在回旋頻率附近產(chǎn)生吸收峰, 并且外加磁場(chǎng)的強(qiáng)度會(huì)改變吸收峰在頻譜上的左右位置, 此外碰撞頻率大小也會(huì)影響太赫茲波的傳輸.總的來(lái)說(shuō)左極化太赫茲波在非均勻磁化等離子體中的傳輸性能要優(yōu)于右極化太赫茲波.這些結(jié)果可為改善和解決再入飛行“黑障”問(wèn)題提供理論參考.