Effect of plasma on transmission characteristics of high-frequency microwave
-
摘要: 研究高频微波在等离子体中的传输特性能有效地分析评估在微波通信和雷达技术中信息的传递过程。通过使用数值仿真的方法分析了等离子体电子密度、厚度及入射波频率对微波反射、吸收和透射的影响。结果显示,等离子体厚度和电子密度增加会导致吸收增强、透射减弱;且反射会随厚度降低和电子密度升高而微弱升高;高频微波更易于穿透等离子体,透射随频率提高而增强。此外,研究结果表明电子密度不仅能影响能量的传输,还会影响电磁波波形,使其展宽。高密度等离子体会明显导致微波波形时空上延展增宽,非弹性碰撞使得增宽现象明显。波形的改变规律能为雷达回波和微波通信所携带信息的复原工作提供一定的理论支撑。Abstract: Studying the transmission characteristics of high-frequency microwaves in plasma can effectively analyze and assess the information transfer process in microwave communication and radar technology. Numerical simulations were employed to analyze the effects of plasma electron density, thickness, and incident wave frequency on microwave reflection, absorption, and transmission. The results indicate that increased plasma thickness and electron density lead to enhanced absorption and reduced transmission; reflection increases slightly with decreased thickness and increased electron density. High-frequency microwaves are more easily transmitted through plasma, with transmission enhancing as frequency increases. Furthermore, the research shows that electron density not only affects energy transmission but also alters the electromagnetic wave shape, causing its broadening. High-density plasma significantly broadens microwave waveforms both spatially and temporally, with inelastic collisions contributing prominently to this broadening. The patterns of waveform changes can provide theoretical support for the restoration of information carried by radar echoes and microwave communications.
-
在微波雷达探测与高频微波通信等现代电子技术领域,高频微波的传输特性对系统性能具有至关重要的影响。随着科学技术的发展,通信频率越来越高,电磁波能量越来越大。一些研究中发现高能、高频的电磁波会激发周围空气电离产生等离子体[1-4]。其独特的电磁性质逐渐引起了科研人员的广泛关注。等离子体,作为由大量带电粒子(如电子、离子)组成的非束缚态系统,其内部复杂的电磁相互作用对高频微波的传输会产生显著影响[5-7]。高频微波与等离子体的相互作用是一个复杂而有趣的研究课题[8-10]。国内外已有研究表明,等离子体的电子密度、温度、碰撞频率等参数会显著影响微波在其中的传输特性[11-13],包括反射、吸收、透射以及色散效应等[14-16]。特别是当微波频率与等离子体中的某些特征频率相近时,会发生共振现象,导致微波能量的高效转换或衰减[17],从而出现一些电磁波增强、增宽的现象[18-20]。使用实际实验的方法往往需要大量的资源,为此数值仿真应运而生,数值仿真目前主流方法有:微分方程模拟、PIC模拟等[21-23]。PIC模拟虽能极精确仿真粒子间的相互作用,但这也导致PIC模拟需要极高的计算能力,因此PIC只能模拟少量的粒子。要模拟等离子体,最理想的方法是用微分方程[24-27]。通过对比微分方程模拟结果和实验数据,研究发现数值仿真的分析方法是较为可靠的[28-30]。整体而言,国内外研究已有一定成果,重点分析能量的传输特性[31-33],对于电磁波波形的分析较少[34]。波形的改变会严重影响微波携带的信息,目前国内外研究结果一定程度上能说明波形的改变和能量传递密不可分[35]。一些研究中提到等离子体鞘套与微波的相互作用随微波电场幅值的变化呈现非线性[36]。本文旨在探讨等离子体对高频微波传输特性的影响,通过数值仿真分析了传播过程中波形的变化规律。结合理论对传播过程进行了解析。相比于已有研究,本工作不仅仅分析能量的传播,更重要的是分析了波形的改变。在不同参数(等离子体电子密度、厚度、入射波频率)下进行仿真,分析高频微波反射、吸收和透射的传播规律。同时,研究中发现等离子体会使微波时空增宽,通过将等离子体看作很多层,揭示了等离子体对微波传输特性的影响。这一工作一定程度上满足了微波通信和微波雷达对微波波形分析的需求。根据本文推导的波形改变规律,可以将接收到的微波波形进行适当变换,得到更为准确的目标信息。也可以通过优化等离子体的参数配置,实现对高频微波传输的调控,从而提升相关电子系统的性能与稳定性。本研究不仅有助于深化对等离子体电磁特性的理解,还为微波通信、雷达探测等领域的技术创新提供了一定理论数据支持。
1. 微波与等离子体相互作用的仿真模型
数值仿真使用的理论模型为麦克斯韦方程组。
∇×H=J+∂D∂t (1a) ∇×E=∂B∂t (1b) ∇⋅B=0 (1c) ∇⋅D=ρf (1d) {B=μrμ0HD=εrε0E (2) 方程(2)为电磁场的本构关系,麦克斯韦方程组的推导要求提供
B 和H 、D 和E 之间的关系。再使用公式{J=σE∇×A=B (3) 将所有方程代入方程(1a)中得到
∇×∇×Aμr+σμ0∂A∂t+μ0∂∂tε0εr∂A∂t=0 (4) 上述方程中:
H 为磁场强度,J 为电流密度,D 为电位移矢量,E 为电场强度,B 为磁感应强度,ρf 为自由电荷密度,μr 为相对磁导率,μ0 为真空磁导率,εr 为相对介电常数,ε0 为真空介电常数,σ 为复电导率,A 为电磁场矢势。目前,方程(4)找不到解析解,因此采取数值方法求解。对于等离子体而言。相对磁导率取1,复电导率可由相对介电常数求出,公式为σ=iε0ω(εr−1) 。ω 表示入射电磁场角频率。则求解方程(4)仅需要一个参数:相对介电常数εr ,即可求得电磁场矢势A 。此时要求解电场强度仅需假设电磁场规范(A,φ )中φ 在零点为0即可。本研究中采用COMSOL软件进行仿真,数值求解方程(4)采用时域有限差分法。在分析等离子体的数值方法中较为常见的是将等离子体看作流体,求解纳维-斯托克斯方程和各类守恒方程。看作流体时仿真能反映等离子体流动和密度变化等特性。本研究使用的方程(4)相较于纳维-斯托克斯方程而言,能更加准确反映电磁场的瞬时变化特性,且求解速度更快,对计算机的要求更低。综合而言,更适用于分析等离子体和电磁波的相互作用过程。几何模型及网格设置如图1所示。入射电磁波使用电场定义。模型中间为等离子体区域,两边为空气。等离子体密度采取二次分布,取边缘为中心密度的1%。根据CFL条件(the Courant-Friedrichs-Lewy condition),网格最大单元格取波长的1/6,时间步长取周期的1/6。所有计算前设置空白对照,即等离子体区域为空气。频率范围选取100 GHz~1 THz。结果如图2所示,从图2中可看出,刚开始时入射电场脉冲紧凑,随着时间推进,脉冲到达仿真区域中间,波形前端和后端未出现明显变形和展宽。最后入射电场几乎无改变地穿过。入射电场振幅前后皆为1 V/m。结果表明电场无变化,能量反射率为0.0%,透射率为100.0%。其余计算与空白对照做比较以排除模型本身的影响。
2. 数值模型运算结果与分析
2.1 等离子体厚度对传播的影响
改变等离子体厚度分析微波传播受到的影响。入射波长取0.3 mm,等离子体密度取1021 m−3,等离子体厚度取20 mm和26 mm,结果如表1所示。已有结果表明,等离子体越厚吸收越强,透射越低。能量衰减的主要原因是微波对电子做功,电子与中性粒子碰撞吸收能量而内能升高,特别是反射率降低。设置多组实验,每组实验仅改变厚度,不同实验组之间改变入射微波频率、等离子体特征频率等,结果表明均为上述趋势。分析可知模型中理想空气对电磁场无吸收。当等离子体变厚时,等离子体内粒子增多,系统自由度升高,电子与中性粒子对入射微波的吸收升高,从而能被反射的波变少。
表 1 不同等离子体厚度时反射、吸收、透射Table 1. Reflection, absorption and transmission in different plasma thicknessplasma thickness/mm reflectivity/% absorptivity/% transmissivity/% 20 1.2 74.8 24.0 26 0.9 91.1 8.10 2.2 等离子体电子密度对传播的影响
改变电子密度分析微波传播受到的影响。等离子体密度从1019 m−3开始取至1021 m−3,等离子体厚度设置为26 mm,入射波波长为0.3 mm,结果如图3所示。结果表明:在其他条件相同时,随着等离子体电子密度的增加,反射率在微弱地升高,透射率从95.88%逐步降低至8.10%,将等离子体看作很多按照密度划分的层。假设单层密度满足一定值才会反射,这是合理的。从粒子的角度分析,粒子过少则不能与入射光子产生相互作用,更不会发生反射。当电子密度整体升高时,电子与中性粒子变多使得总体吸收加强。同时电子密度整体各层都升高,使能产生反射的电子层变多。从而整体显现出反射率微弱升高的现象。
2.3 入射波频率对传播的影响
升高入射波频率分析对结果造成的影响。在等离子体密度取1019 m−3、厚度为1 mm时,入射波选取0.1~5 THz,结果如图4所示。结果表明:随着入射波频率提高,反射降低,透射变强,吸收降低。当入射波频率高达1 THz时,等离子体特征频率为0.178 THz。此时入射波频率相比于等离子体频率较高,入射波透射率可达99.5%以上,几乎全部透射。当入射波频率高到一定程度后,反射和透射都变得稳定。从整体上看,入射电磁场频率很高时,整体电场来不及通过做功将能量传递给粒子,从而粒子与光子相互作用变弱,反射和吸收都变弱,从而透射增强。结果表明高频微波很容易穿透等离子体传播。
2.4 等离子体对微波波形的影响
设置入射微波频率为1 THz,则不同电子密度时入射波和透射波的场强如图5所示。由图5可知,当等离子体电子密度为单一定值1019 m−3 时,此时密度较低,等离子体区对电磁场的吸收很弱,约为1%,所以电磁场穿过等离子体后场强随时间的变化图和电磁场穿过空气后场强随时间的变化图几乎一致;而当电子密度提高到1020 m−3时,等离子体区域吸收增强,当提高至1021 m−3时,等离子体区吸收接近90%。从图5(d)中可以看到等离子体内电磁场的增宽现象。这可以使用经典理论和量子理论解释。
基于电磁辐射理论,假设强度降低为原来的1/e所花费时间为
τ ,其表示波的寿命。当电子密度升高、碰撞频率升高时,有些波列会被打断,使得寿命降低。假设其关系式为τ=Cf(n,ν0) (5) 式中:
C 为常数,n为等离子体电子密度,ν0 为碰撞频率。上述假设关系表明函数f(n,ν0) 应为变量n、ν0 的单调递减函数。即当n、ν0 变大时,寿命τ 会变小。结合图5和经典理论则电磁波波形可表示为E(t)={0,t<0E0e−tτ×ei2πf0t,t⩾ (6) 式中:
{f}_{0} 表示入射波频率。式(6)描述了电场随时间的变化关系。选取计算时刻为起点,随着时间的流逝,电场振幅包络以e指数衰减,同时电场存在频率为{f}_{0} 的振荡。对电磁波波形进行傅里叶变换得到频域表达式E\left(f\right)={\displaystyle\int }_{-\infty }^{+\infty }E\left(t\right){\mathrm{e}}^{-\mathrm{i} 2\pi ft}{\mathrm{d}}t=\dfrac{{E}_{0}}{\mathrm{i} 2\pi \left(f-{f}_{0}\right)+\dfrac{1}{\tau }} (7) 其对应的的宽度为
{{ \Delta}} f=\dfrac{1}{2\pi \tau }=\dfrac{{C}f\left(n,{\nu }_{0}\right)}{2\pi } (8) 理论上可以得知,随着密度和频率的升高,电磁波的波形会增宽。结合量子理论不确定关系
{{ \Delta}} E \tau \geqslant \dfrac{\mathrm{\hslash }}{2} (9) 可以得知当密度和碰撞频率升高时,每个能级的宽度
{{ \Delta}} E 会增宽。{{ \Delta}} f=\dfrac{{\Delta E}_{1}+{\Delta E}_{2}}{h} (10) 由式(10)可知大量的跃迁过程就会展现总体上的增宽现象。
整体上看,随着电子密度升高,越来越多的粒子与入射波发生作用,不同密度层的粒子作用时间不一,整体电场空间上延展,再者等离子体内大量相互作用为非弹性碰撞,从而出现碰撞增宽的现象。从图5看,最终结果中透射电场在时空上有增宽。
3. 结 论
本文通过数值仿真探究了高能微波在等离子体中的传输特性,得出了如下结论:
(1)将等离子体看作很多按照密度划分的层。假设单层密度满足一定条件才会反射。站在粒子的角度看,当粒子过少时,入射光子几乎无法与电子相互作用,即不会发生吸收和反射。当粒子密度提高到一定值时,入射光子与电子相遇发生作用。站在电磁波的角度看,随着入射电磁场频率的升高,电子的运动会跟不上电场的频率。从而解释了等离子体厚度、电子密度、入射波频率改变时,微波与等离子体相互作用的能量变化。
(2)等离子体的电子密度会大幅度改变微波的形状。当电子密度较高时,微波的波形会在时空上延展增宽。这是因为密度的升高会使得更多的电子能与入射波发生作用。不同密度层作用时间不一,从而空间上延展。再者等离子体中碰撞多为非弹性碰撞,这会使得碰撞增宽现象较为明显。整体上看,波形变得混乱且时空上增宽。
本文结果也揭示着在未来研究中,随着高能微波技术的不断发展,会有越来越多的研究发现等离子体对微波的影响。随着能量的不断提高,研究二者之间的相互作用过程对微波的进一步应用是非常重要的。
-
表 1 不同等离子体厚度时反射、吸收、透射
Table 1. Reflection, absorption and transmission in different plasma thickness
plasma thickness/mm reflectivity/% absorptivity/% transmissivity/% 20 1.2 74.8 24.0 26 0.9 91.1 8.10 -
[1] Poorreza E, Gargari N D. Modeling and simulation of a microwave-assisted plasma with different input power for plasma-based applications[J]. Russian Journal of Physical Chemistry B, 2023, 17(3): 719-724. doi: 10.1134/S1990793123030235 [2] Ikeda Y, Soriano J K, Kawahara N, et al. Spatially and temporally resolved plasma formation on alumina target in microwave-enhanced laser-induced breakdown spectroscopy[J]. Spectrochimica Acta Part B: Atomic Spectroscopy, 2022, 197: 106533. doi: 10.1016/j.sab.2022.106533 [3] Zhang Chaoyang, Fu Wenjie, Hu Shijie, et al. Investigation on continuous and modulated microwave plasma filaments at atmospheric pressure[J]. IEEE Access, 2021, 9: 154318-154323. doi: 10.1109/ACCESS.2021.3128610 [4] Barnes B K, Ouro-Koura H, Derickson J, et al. Plasma generation by household microwave oven for surface modification and other emerging applications[J]. American Journal of Physics, 2021, 89(4): 372-282. doi: 10.1119/10.0002706 [5] Zhao Chengwei, Li Xiaoping, Liu Yanming, et al. Research on plasma electron density distribution based on microwave diffraction[J]. Plasma Sources Science and Technology, 2022, 31: 015007. doi: 10.1088/1361-6595/ac39ad [6] Ye Xin, Wang Yongge, Yao Jingfeng, et al. Plasma-enabled microwave modulation for continuous beam scanning[J]. Journal of Physics D: Applied Physics, 2022, 55: 435202. doi: 10.1088/1361-6463/ac8da1 [7] Zhang Jianwei, Luo Wei, Jiang Ming, et al. Plasma propagation in the microwave window breakdown at the air/dielectric interface[J]. Plasma Sources Science and Technology, 2020, 29: 025013. doi: 10.1088/1361-6595/ab6e59 [8] Mu Xiangchao, Dong Guoxiang, Li Xiaoping, et al. Analysis of the electromagnetic wave transmission characteristics in inhomogeneous plasma based on an equivalent circuit model[J]. IEEE Transactions on Plasma Science, 2024, 52(3): 698-706. doi: 10.1109/TPS.2024.3371100 [9] Noori E. Investigation of near cut-off properties of electromagnetic wave propagation in homogeneous, collisional plasma slab[J]. Contributions to Plasma Physics, 2022, 62: e202200016. doi: 10.1002/ctpp.202200016 [10] Zhang Chaoyang, Chen Chi, Fu Wenjie, et al. Investigation on the microwave excited plasma filament at atmospheric pressure[J]. IEEE Transactions on Plasma Science, 2021, 49(6): 1877-1881. doi: 10.1109/TPS.2021.3075574 [11] Zhao Yuexing, Su Ruiming, Ma Jinping, et al. SO-FDTD simulation on the transmission characteristics of terahertz waves in inhomogeneous magnetized dusty plasma[J]. The European Physical Journal D, 2024, 78: 32. doi: 10.1140/epjd/s10053-024-00824-8 [12] Wan Xiaohuan, Zhou Zhikun, Zhang Juan, et al. Propagation characteristics of obliquely incident terahertz waves in high-temperature magnetized plasma[J]. IEEE Transactions on Plasma Science, 2022, 50(2): 241-249. doi: 10.1109/TPS.2021.3139373 [13] Zhou Zhikun, Wan Xiaohuan, Zhang Juan, et al. Influence of temperature on terahertz waves propagating in magnetized plasma[J]. Physica Scripta, 2021, 96: 075607. doi: 10.1088/1402-4896/abfcf1 [14] Shen Fangfang, Zhang Zhongdao, Bai Bowen, et al. Research on the reflection characteristics of the broadband electromagnetic wave in nonuniform plasma[J]. IEEE Transactions on Plasma Science, 2024, 52(3): 657-665. doi: 10.1109/TPS.2024.3373645 [15] Zhang Jie, Li Miao, Han Bing. Analysis of electromagnetic waves reflected by re-entry plasma sheath based on CSO-FDTD[J]. Physica Scripta, 2023, 98: 095610. doi: 10.1088/1402-4896/acf0f9 [16] Zhang Youyi, Xu Guanjun, Zheng Zhengqi. Terahertz waves propagation in an inhomogeneous plasma layer using the improved scattering-matrix method[J]. Waves in Random and Complex Media, 2021, 31(6): 2466-2480. doi: 10.1080/17455030.2020.1757177 [17] Lyu Xingbao, Yuan Chengxun, Avtaeva S, et al. Attenuation of microwave radiation by post-anode plasma in a composite grid electrode structure[J]. IEEE Access, 2022, 10: 7675-7683. doi: 10.1109/ACCESS.2022.3143582 [18] Vhanmore B D, Rajmane S P, Sadale S B, et al. Dominance of polarization modes and absorption on self-focusing of laser beams in collisionless magnetized plasma[J]. Journal of Nonlinear Optical Physics & Materials, 2024. [19] Chen Peiqi, Nie Qiuyue, Zhang Zhonglin, et al. Integrative implementation of scattering reduction and radiation enhancement for an electrically small antenna by subwavelength plasmas[J]. Physics of Plasmas, 2024, 31: 073503. doi: 10.1063/5.0211317 [20] Metelskii I I, Kovalev V F, Bychenkov V Y. Nonlinear laser radiation absorption due to relativistic plasma resonance in an inhomogeneous plasma[J]. Journal of Experimental and Theoretical Physics, 2021, 133(2): 236-252. doi: 10.1134/S1063776121080069 [21] Sun Jinhai, Zhao Yarui, Yin Hongcheng, et al. An open simulation model for terahertz wave transmission in plasma[J]. Journal of Applied Physics, 2024, 135: 153104. doi: 10.1063/5.0199442 [22] Desai M, Ghosh P, Kumar A, et al. Deep-learning architecture-based approach for 2-D-simulation of microwave plasma interaction[J]. IEEE Transactions on Microwave Theory and Techniques, 2022, 70(12): 5359-5368. doi: 10.1109/TMTT.2022.3217138 [23] Janicek A, Thornton E, Garrett T, et al. Length dependence on broadband microwave emission from laser-generated plasmas[J]. IEEE Transactions on Plasma Science, 2020, 48(6): 1979-1983. doi: 10.1109/TPS.2020.2988762 [24] Chang Qi, Ma Yunpeng, Liu Ji, et al. Simulation of electromagnetic waves in plasma by subdomain-level nonconformal DGTD method[J]. IEEE Transactions on Plasma Science, 2022, 50(11): 4775-4781. doi: 10.1109/TPS.2022.3214215 [25] Zhou Zhikun, Wan Xiaohuan, Li Xiaolin, et al. SO-FDTD analysis on transmission characteristics of terahertz wave in plasma[J]. Physics of Plasmas, 2021, 28: 072105. doi: 10.1063/5.0053611 [26] Wang Ge, Pan Hui, Lai Shimiao, et al. Dynamic measurement of relative complex permittivity of microwave plasma at atmospheric pressure[J]. Processes, 2021, 9: 1812. doi: 10.3390/pr9101812 [27] Yao Jingfeng, Yu Zhi, Yuan Chengxun, et al. The influence of plasma distribution on microwave reflection in a plasma-metal model[J]. IEEE Transactions on Plasma Science, 2020, 48(2): 359-363. doi: 10.1109/TPS.2019.2943519 [28] Liu Yanming, Zhang Xi, Bai Bowen, et al. A calculation method of electromagnetic wave reflection in plasma sheath environment[J]. IEEE Transactions on Plasma Science, 2022, 50(7): 2030-2038. doi: 10.1109/TPS.2022.3181220 [29] Zhang Chaoyang, Lu Dun, Hu Shijie, et al. An economic real-time microwave plasma impedance measurement method[J]. IEEE Transactions on Plasma Science, 2021, 49(11): 3503-3508. doi: 10.1109/TPS.2021.3115846 [30] Pak I H, Kim Y H, Oh C H, et al. Transmission characteristics of electromagnetic waves in a semicircular plasma filament layer generated by a femtosecond laser[J]. Contributions to Plasma Physics, 2021, 61: e202100032. doi: 10.1002/ctpp.202100032 [31] 王彦同, 帕提曼·阿不都玛洪, 石雁祥, 等. 弱电离尘埃等离子体微波衰减理论的实验研究[J]. 电波科学学报, 2020, 35(6):967-973Wang Yantong, Abudoumahong P, Shi Yanxiang, et al. Experimental study on microwave attenuation theory of weakly ionized dusty plasma[J]. Chinese Journal of Radio Science, 2020, 35(6): 967-973 [32] Moshkov A V, Pozhidaev V N. Numerical simulation of very-low-frequency waves passing through the magnetoactive plane-layered plasma of Earth's lower ionosphere[J]. Journal of Communications Technology and Electronics, 2020, 65(5): 472-479. doi: 10.1134/S1064226920050101 [33] Ma Zhu, Wei Min, Li Meng, et al. Study on electromagnetic characteristics of plasma model-based on the symplectic multiresolution time-domain scheme[J]. Modern Physics Letters B, 2020, 34: 2050046. [34] Bao Yu, He Xiang, Su Wei, et al. Study on the generation of terahertz waves in collision plasma[J]. Physics of Plasmas, 2024, 31: 093302. doi: 10.1063/5.0219947 [35] 潘惠, 王舸, 杨阳. 大气压微波等离子体射流装置放电特性研究[J]. 强激光与粒子束, 2022, 34:049001 doi: 10.11884/HPLPB202234.210277Pan Hui, Wang Ge, Yang Yang. Design and study of atmospheric pressure microwave plasma jet[J]. High Power Laser and Particle Beams, 2022, 34: 049001 doi: 10.11884/HPLPB202234.210277 [36] 陈煜青, 王蕾, 赵立山, 等. 等离子体鞘套低频通信电磁波透射率与辐照微波场强关系仿真研究[J]. 强激光与粒子束, 2023, 35:089001 doi: 10.11884/HPLPB202335.220361Chen Yuqing, Wang Lei, Zhao Lishan, et al. Simulation study of the relationship between low-frequency communication EM wave transmissivity of plasma sheaths and irradiation microwave E-field strength[J]. High Power Laser and Particle Beams, 2023, 35: 089001 doi: 10.11884/HPLPB202335.220361 -