李郝 杨鑫2)3) 张正平2)3)†
1) (贵州大学大数据与信息工程学院, 贵阳 550025)
2) (半导体功率器件可靠性教育部工程研究中心, 贵阳 550025)
3) (贵州省微纳电子与软件技术重点实验室, 贵阳 550025)
为利用太赫兹波解决飞行器再入过程遇到的“黑障”问题, 以散射矩阵方法为基础, 分别以非均匀磁化等离子体的磁化方向、电子密度、外加磁场强度和碰撞频率为变量, 研究了垂直入射情形下它们对太赫兹波传输行为的影响.结果表明: 这些参数对太赫兹波传输性能影响明显, 例如按某一方向改变磁化角度对左极化和右极化太赫兹波的传输功率有相反的影响; 降低磁化强度能一定程度地避开等离子体对右极化波的吸收; 而降低碰撞频率能缩小等离子体对右极化波的吸收频带.通过调整这些参数, 有望在一定程度上缓解黑障现象.
高超音速飞行器再入大气过程中会与空气发生摩擦, 在产生大量热量的同时飞行器表面的空气会发生电离, 生成一层薄等离子体鞘套[1].飞行器的天线窗口被等离子体覆盖后, 其发射和接受信号的能力大大下降, 严重的时候通信会完全中断, 产生“黑障”现象[2].
为了解决“黑障”问题, 大量的科研人员对等离子体特性[3]和电磁波在等离子体中的传播行为进行了研究.研究人员发现: 电磁波频率低于等离子体频率是导致“黑障”问题的重要原因之一, 因而就连高频段的毫米波也无法完全穿透飞行器表面形成的高电子密度等离子体鞘套.近年来THz 波技术不断进步, 其优秀的穿透性和抗干扰性使得“黑障”问题有了解决的可能, 因此, 研究THz 波在等离子体鞘套中的传播有着重要意义.另外, 研究人员开发了多种算法来分析电磁波与等离子体的相互作用, 如时域有限差分法(FDTD)[4−7]、解析法[8]、温泽尔-克莱默-布里渊法(WKB)[9,10]以及散射矩阵法(SMM)[11,12]等.本文采用SMM 方法进行研究分析, SMM 方法较好地运用了分层思想来实现对等离子体不均匀性的体现, 能直接高效地计算出总的反射、透射和吸收率, 相较于FDTD 方法也不需要进行大量的网格化操作, 所以对于飞行器表面等离子体鞘套的分析, SMM 方法具有较高的效率和较低的算法复杂度.正因如此, Chen 等[13]用SMM方法研究了时变等离子体鞘层对雷达回波信号的影响, Guo 和Guo[14]用SMM 方法分析了THz 波在非均匀磁化热等离子体中的传输特性, Zhang 等[15]用SMM 方法研究了THz 波在右极化并且外加磁场角度与入射角度平行时, 在非均匀磁化等离子体中的传输特性.我们发现, 这些研究大多只是基于固定的外加磁场角度对传输特性进行研究, 没有考虑外加磁场方向的变化对THz 波传输的影响.一些论文也只是基于单一极化方向的分析, 然而不同的极化方向也是影响THz 波传输性能的因素.所以, 为了更全面地探究磁化角度的改变对左极化和右极化THz 波在非均匀磁化等离子体中传输的影响, 本文加入磁化角度作为分析传输特性的变量之一, 使用SMM 方法分别研究了左极化和右极化THz 波与等离子体之间不同的作用特点.
对于分布不均匀的等离子体可以设定为图1所示的分层模型, 将等离子体平板平均分成n层等离子体薄片, 每一薄片中的电子密度都视为均匀分布, 总的等离子体层电子密度按给定的密度函数分布.电磁波垂直入射等离子板表面, 入射区(0)和透射区(p)都为自由空间.
图1 太赫兹波在分层等离子体中传播的模型图Fig.1.Model diagram of THz waves propagation in the layered plasma.
在此模型中每一层薄片等离子体中都有其传播常数k(m):
其中, 正负号代表左极化或者右极化,θ是波传播方向与磁化方向的夹角,X=1−iv/ω,v是碰撞角频率,Y=ωce2/ω2,ωce是电子回旋角频率,Z=是 第m 层的等离子体角频率并由(3)式表示[20]:
其中Ne为自由电子密度,e和me分别为电子电荷量和电子质量,ε0为真空中的介电常数.
将图1 中入射区域(0)和薄层(1)的电磁场,用入射(透射)电场、反射电场、入射(透射)磁场和反射磁场表示, 并分别写为以下形式:
其中A是射入到第一个等离子层表面的反射系数,B1和C1为电磁波从区域(0)射入等离子体薄层(1)的透射系数和薄层(2)反射回薄层(1)的反射系数,k(0)是自由空间的波数.这样可以确定入射区域(0)和薄层(1)的总电磁场关系, 通过匹配x=0处分界面的边界条件能得到矩阵关系:
其中
和上面相同, 继续匹配第2 到n个分界面能得到矩阵关系:
其中
在最后一个分界面的左右是薄层(n)和自由空间(p), 在区域(p)中只有透射电磁波, 没有反射电磁波, 所以匹配x=d处的边界条件能得到矩阵关系:
其中
通过上面的(6)式、(8)式和(10)式可以递归变换得到:
其中
接着将Sg表示为Sg=(Sg1,Sg2) , 其中Sg1和Sg2表示Sg的第一列和第二列向量, 所以(12)式可以变换为
最后解得系数A和D的矩阵:
在接下来的分析中将电磁波功率归一化处理,归一化反射功率表示为|A|2, 归一化透射功率表示为|D|2, 那么归一化吸收功率为 1−|A|2−|D|2.
根据NASA 钝头RAM 飞行器的公开数据[16],在不同高度下有不同的电子密度分布, 并且通常来说可以用指数分布、高斯分布、抛物分布等模型来模拟等离子体中的电子密度分布[21].但单一的分布模型不能模拟实际分布的非对称性, 所以本文采用双高斯分布模型来模拟计算以更好地贴近实际分布情况, 分布函数表示如(16)式, 函数图像如图2 所示.
图2 电子密度分布与等离子体厚度的关系Fig.2.Electron density profile versus plasma thickness.
根据现有飞行器表面形成的等离子体鞘套的公开数据和流场计算[22], 在此设置等离子体平板的厚度为0.1 m, 为保证计算收敛性将其分成50 层薄片, 碰撞频率fv=0.1 THz , 最大电子密度N0=1×1018m−3, 外加磁场强度B=5 T.根据(2)式, 可以对外加磁场的角度进行分析.
图3 THz 波归一化传输功率与外加磁场角度的关系 (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 可以看出, 随着磁场角度由0°增加到90°时, 右极化THz 波传输功率基本逐渐升高, 同时左极化THz 波的传输功率基本逐渐降低.由此可以看出, 当磁场方向与波的传播方向垂直时, 最有利于右极化太赫兹波信号的接收; 当磁场方向与波的传播方向平行时, 最有利于左极化太赫兹波信号的接收.另外, 当电磁波在外加静态磁场的等离子体中传播时, 等离子体中的电子将受到磁场作用做回旋运动, 而右极化波会与回旋运动的电子同步旋转形成回旋共振, 导致能量吸收.所以在相同磁化角度情况下, 左极化波和右极化波的传输功率曲线产生了不同的变化趋势.
那么以磁化角度分别为90°和0°来分析右极化太赫兹波和左极化太赫兹波在非均匀磁化等离子体中的传输特性.当磁化角度θ=90°且为右极化时相对介电常数如(17)式; 当磁化角度θ=0°且为左极化时相对介电常数如(18)式.
在这两种条件下就可以进一步分析在不同电子密度、碰撞频率和外加磁场强度条件下右极化和左极化THz 波在等离子体中的传输特性.
图4(a)—图4(c)分别给出了右极化THz 波在不同最大电子密度下的反射、透射和吸收功率.其中等离子体平板的总厚度d=10 cm , 碰撞频率fv=0.1 THz , 外加磁场强度为B=5 T , 最大电子密度的值分别为N0=1×1016m−3,N0=1×1017m−3,N0=1×1018m−3和N0=1×1019m−3.可以看出, 较高电子密度的等离子体对太赫兹波的反射能力非常弱.随着电子密度的增加, THz 波传输功率相应降低, 且吸收功率增加, 在电子密度为N0=1×1019m−3时出现了传输阻带, 阻带内几乎所有的电磁波功率都被等离子体层吸收.另外, 当电子密度N0<1×1019m−3时, 右极化电磁波的吸收功率在0— 0.6 THz 内先增加后减少, 在0.1—0.2 THz间出现一个突起峰值, 此时外加磁场条件下的电子回旋频率fce=0.14 THz , 正好在峰值区域, 所以可以推断由于右旋波的偏振方向与电子回旋方向相同, 引起了电子回旋共振, 导致能量吸收增强, 并形成了突起峰值.而电子密度的增加也引起了吸收峰值的升高, 这是因为电子密度的增加意味着等离子体中电子数量的增加, 从而导致电子从电磁波中吸收的能量增加.图4(d)—图4(f) 分别给出了左极化THz 波在不同最大电子密度下的反射、透射和吸收功率, 参数与右极化THz 波保持一致.可以看到的是, 由于左旋波的偏振方向与电子回旋方向不同, 左极化THz 波相对于右极化THz波最大的差异在于没有了电子回旋共振的影响, 并且在相同的电子密度下左极化THz 波对等离子体有更好的穿透性能.
图4 THz 波归一化反射、传输、吸收功率与电子密度最大值的关系 (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 波在等离子体中的反射功率非常低,传输功率和吸收功率的变化规律相反, 但传输功率在一定程度上能更好地体现信号传输的能力, 所以接下来主要对传输功率进行分析.
图5 THz 波归一化传输功率与外加磁场强度的关系 (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 波归一化传输功率与碰撞频率的关系 (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)分别给出了外加磁场强度对右极化和左极化THz 波传输功率的影响, 此时等离子体平板的总厚度d=10 cm , 碰撞频率fv=0.1 THz , 最大电子密度的值N0=1×1018m−3.从图5(a)和图5(c)可以看出, 左极化波传输功率在整个频带上随外加磁场强度的增加而增加.而右极化波传输功率随着外加磁场强度的增加在低频端和高频端有着不同的变化趋势: 在低频端右极化波传输功率随外加磁场强度的增加而增加, 但由于磁场强度的增加引起回旋共振的频带向高频端偏移, 进而在高频端产生更多的能量吸收, 所以在高频端右极化波传输功率会随外加磁场强度的增加而降低.
图6(a)和图6(b)以及图6(c)和图6(d)分别给出了碰撞频率对右极化和左极化THz 波传输功率的影响, 此时等离子体平板的总厚度d=10 cm ,外加磁场强度为B=5T , 最大电子密度的值N0=1×1018m−3.可以看出, 右极化波传输功率的最小值随碰撞频率的增加而减小, 但吸收频带增大.左极化波传输功率随着碰撞频率的增大而减少.总体来说左极化波更能穿透磁化冷等离子体, 有助于更好地改善飞行器再入“黑障”问题.
本文以NASA 钝头RAM 实验数据为基础,研究了左极化和右极化太赫兹波在非均匀磁化冷等离子体中的传输特性.结果表明, 外加磁场的角度对左右极化太赫兹波的影响有所不同, 在磁场方向与波传输方向垂直时最有利于右极化太赫兹波的传输; 在磁场方向与波传播方向平行时最有利于左极化太赫兹波的传输.另外, 随着等离子体最大电子密度的增加, 等离子体对太赫兹波的能量吸收增强, 电磁波在右极化时会发生回旋共振, 在回旋频率附近产生吸收峰, 并且外加磁场的强度会改变吸收峰在频谱上的左右位置, 此外碰撞频率大小也会影响太赫兹波的传输.总的来说左极化太赫兹波在非均匀磁化等离子体中的传输性能要优于右极化太赫兹波.这些结果可为改善和解决再入飞行“黑障”问题提供理论参考.