低温等离子体用于高功率微波防护研究*
doi: 10.11887/j.cn.202306012
李志刚 , 邱志楠 , 汪家春 , 刘丽萍 , 王俊儒 , 陈宗胜
国防科技大学 脉冲功率激光技术国家重点实验室, 安徽 合肥 230037
基金项目: 安徽省自然科学基金面上资助项目(1908085MF205)
Study of high-power microwave protection technology based on low-temperature plasma
LI Zhigang , QIU Zhinan , WANG Jiachun , LIU Liping , WANG Junru , CHEN Zongsheng
State Key Laboratory of Pulsed Power Laser Technology, National University of Defense Technology, Hefei 230037 , China
摘要
等离子体对于高功率微波的攻击具有独特的防护效果。基于等离子体流体近似方法,利用COMSOL软件研究了高功率微波与柱状等离子体阵列相互作用过程中入射电场随时间的演变过程,分析了等离子体防护高功率微波的物理过程和作用机理。研究结果表明,入射的高功率微波会使等离子体参数发生剧烈变化,特别是其电子密度将急剧增加,从而使等离子体对入射的高功率微波表现出类似金属的电磁特性,最终实现对入射高功率微波的有效防护。此外,利用高频辉光放电产生柱状等离子体阵列,通过实验验证了等离子体对高功率微波的防护作用。最后,总结了基于等离子体的高功率微波防护技术需解决的主要问题。
Abstract
Plasma has a unique protective effect against high-power microwave attack. Based on the plasma fluid approximation method, the time dependent evolution of the incident electric field during the interaction between high-power microwaves and columnar plasma arrays was studied by using COMSOL software, and the physical process and mechanism of plasma protection against high-power microwaves were analyzed. The results show that the incident high-power microwave will change the plasma parameters drastically, especially the electron density will increase sharply, so that the plasma will show metal-like electromagnetic properties to the incident high-power microwave, and finally realize the effective protection against the incident high-power microwave. In addition, the plasma protection against high-power microwaves was experimentally verified by using columnar plasma arrays generated by high-frequency glow discharge. Finally, the main problems to be solved for plasma-based high-power microwave protection technology were summarized.
高功率微波攻击方式的出现给电子信息系统安全带来了巨大的威胁。在未来信息化条件下,电子信息系统一旦遭受破坏,再强大的武器装备也会变成“瞎子”“聋子”,甚至变为一堆废铁。可见,高功率微波攻击或将改变未来战争的作战样式,对未来战争形态产生重大影响[1]
正是认识到高功率微波攻击的强大作战效能,美、俄、英、日、法等国都非常重视,并已有相关装备报道。以美国为例,已投入使用的装备有主动拒止系统、警惕鹰系统、携带高功率微波弹头的AGM-86巡航导弹以及用于空域封锁的高功率微波炮等[2]。目前,美军在这一领域仍在继续加大投入。因此,提升装备对高功率微波的防护能力具有重要的军事效益和现实意义。
目前基于固态半导体器件的电路级防护,其能防护的最大脉冲功率一般在数千瓦以内(具体数值与微波脉冲的宽度有关)。对于更高功率的微波脉冲,可能的防护方法主要有频率选择表面、能量选择表面和等离子体[3-5]。其中等离子体对高功率微波的能量具有敏感性,能够随着入射高功率微波的强弱改变自身的状态,进而影响入射波的传输,兼具了频率选择表面和能量选择表面的优点,具有很好的研究价值和应用前景[6-8]
本文针对等离子体对高功率微波的防护效果开展了仿真与实验分析,研究了高功率微波与柱状等离子体阵列相互作用过程中入射电场随时间的演变过程,分析了等离子体防护高功率微波的物理过程和作用机理;在此基础上,通过实验验证了等离子体对高功率微波的防护效果;最后结合相关研究工作,对基于等离子体的高功率微波防护技术需解决的主要问题进行了总结。
1 基本理论
为了研究在高功率微波作用下等离子体产生的电磁屏蔽作用,本文采用流体近似方法[9-10]进行仿真分析,分别采用等离子体中的波动方程、电子传递方程和重物质传递方程,对高功率微波在等离子体中的传播过程和等离子体内部电子及其他带电粒子受入射电磁能量的影响进行表征。
1.1 波动方程
当电磁波入射到密度均匀的等离子体上时,电磁波在等离子体中的麦克斯韦方程[11]为:
×H=jωε0εrE
(1)
×E=-jωμ0H
(2)
可得波动方程:
××E-k02εrE=0
(3)
其中:H为磁场强度;E为电场强度;ω为入射波的角频率;ε0μ0k0分别为真空中的介电常数、磁导率和波数;εr为等离子体的相对介电常数,可通过式(4)求解。
εr=1-ωpω(ω-iv)
(4)
式中,ωp为等离子频率,υ为电子碰撞频率,i为复数标号。
1.2 电子传递方程
在外加电磁场的作用下,等离子体内部电子密度的变化可通过电子传递方程来进行分析,方程可表示为:
tne+Γe=Re
(5)
Γe=-μeEne-Dene
(6)
其中:ne为电子数密度;μeDe为电子迁移率和电子扩散率;Re为电子源项,表征了内部碰撞反应导致的电子的产生与消失,可通过式(7)求得。
Re=j=1M xjkjNnne
(7)
式中,xjj反应中碰撞粒子的摩尔质量分数,kjj反应的反应速率,Nn为等离子体中总的粒子数密度。
1.3 重物质传递方程
外加电磁场还会引发等离子体内部其他带电粒子的分布产生变化,可通过重物质传递方程来进行分析,方程为:
ρtωk+ρ(u)ωk=jk+Rk
(8)
式中:ωk为第k种粒子的摩尔分数;ρ为气体密度;u为平均流体速率;jk为第k种粒子的扩散通量,可通过式(9)表示。
jk=ρωkVk
(9)
Vk=Dk,mωkωk+Dk,mMnMn+DkTρωkTT-zkμk,mE
(10)
式中,Dkm为重物质平均扩散系数,Mn为重物质平均摩尔质量,DTkzkμk,m分别为第k种粒子的热扩散系数、电荷数以及平均迁移率,T为气体温度。
2 仿真计算模型
采用COMSOL软件进行计算,具体模型如图1所示。模型模拟的物理过程为:电磁波从上边界进入模型,向下传播,经空气介质和等离子体介质后,传输至设置的完美匹配层(避免波反射产生的二次激励)被完全吸收。等离子体层由多根前后延伸的柱状等离子体单元紧密排列组成,假设每根等离子体单元的参数完全一致。
1高功率微波与柱状等离子体阵列相互作用的仿真计算模型示意图
Fig.1Simulation schematic diagram of the interactions between high-power microwave and columnar plasma arrays
假设入射电磁波场强为E0,电磁波频率为f,选择常用的惰性气体Ar作为工作气体,所涉及的粒子种类和碰撞反应在表1给出[12],气体压强为P
1氩等离子体内部碰撞反应方程及类型
Table1Collision reaction functions and types inside Ar plasma
3 仿真计算结果分析
假设初始时刻气体温度为300 K,压强为400 Pa,电子初始为均匀分布,密度为1.0×1016 m-3,等离子体直径为25 mm,玻璃管壁厚度为0.3 mm,初始电子能为3 eV,碰撞频率为5×109 Hz,高功率微波频率为6 GHz,电场强度为2×106 V/m,极化方向与等离子体管轴向方向一致,计算时间设置为0~0.001 s。
图2为不同时刻计算模型中的电场分布,图中采用统一的颜色图例(E:0~4×106 V/m)。高功率微波脉冲从上边界入射,经空气层和等离子体层后被完美匹配层吸收。
图2可以看出,在t=2.15×10-10 s时刻,等离子体区域的电场值与上下空气区域基本一致,说明在初始时刻入射波几乎可以无损透过等离子体区域。从t=3.98×10-10 s到t=5.99×10-10 s,可以明显看到等离子体下方空气区域的电场值越来越小,上方空气电场峰值不断增大。说明这一时间内,等离子体与高功率微波发生了剧烈的相互作用。从t=9.03×10-10 s到t=0.001 s可以看到,等离子体上方区域存在明显的驻波场,下方区域仅有很少的能量透过,电场值接近于零,等离子体类似于金属,对入射电磁波表现出较强的屏蔽特性。
为了更好地解释上述结果的产生原因,图3给出了计算模型中轴线电子密度和电场分布的演变曲线,图中灰色区域代表玻璃放电管的物理尺寸。
图3中表征的物理过程可描述为:①初始阶段,等离子体内部电子受到高功率微波的加速作用,由低能电子转化为高能电子,并通过弹性碰撞反应,将吸收的电磁能量传递给中性粒子,提高等离子体的内能,在这一阶段非弹性碰撞反应并不显著,电子密度增加缓慢,等离子体本身的电磁参数基本保持不变,入射波几乎无损透过等离子体区域;②随着高能电子不断积累和等离子体内能不断增加,等离子体中非弹性碰撞反应速率显著增大,大量的高能电子被消耗,平均电子能急剧下降,电子雪崩效应产生,大量新生电子参与到等离子体振荡中,等离子体电导率急剧增大,对入射波的衰减急剧增大,等离子体电磁特性发生较大改变,其对入射波的屏蔽作用开始显现,透过等离子体区域的电磁能量不断减小;③当高能电子被大量消耗,直接电离反应(见表1中反应4)开始减弱,间接电离反应(见表1中反应5)增强,激发态氩原子被大量消耗用于继续产生新的带电粒子,等离子体电磁特性进一步发生改变,电导率继续增大,入射波在等离子体中趋肤深度不断减小,等离子体对入射波的屏蔽性能更为显著,透过等离子体区域的电磁能量进一步减小;④随着等离子体区域电子密度不断增加,入射的电磁能量除小部分被吸收外,大部分被反射,等离子体内部碰撞反应趋于缓和,电子密度增加速率变缓,空间分布开始变得均匀,等离子体内部逐渐达到新的平衡。
2计算模型中电场强度空间分布随激励时间的演变过程
Fig.2Time dependence of the electric field intensity in simulation calculation model
3模型中轴线电子密度和电场分布随激励时间的演变过程
Fig.3Time dependence of the axial electron density and electric field distribution in simulation calculation model
从粒子平衡角度对这一过程进行解释:等离子体与高功率微波之间的相互作用是通过状态改变来维持粒子平衡的。等离子体作为特殊的介质,内部存在一定的粒子平衡,并且具有强烈维持这一平衡的意愿。高功率微波的入射破坏了等离子体中的粒子平衡,等离子体内部随之产生强烈的振荡和碰撞反应,建立新的状态来平衡外加电磁场产生的影响。
从上面的物理过程分析可知,入射的高功率微波会使等离子体的参数发生剧烈变化,特别是其电子密度将急剧增加,这会引起等离子体电磁参数发生同步改变,从而使等离子体对入射的高功率微波表现出类似金属的电磁特性,最终实现对入射高功率微波的有效防护。
4 实验验证
4.1 实验方案
图4为实验设备连接示意图。实验中,高功率微波信号由L波段高功率微波源产生,经定向耦合器后分成两路信号:一路信号经同轴线缆接入高精度示波器,作为参考信号;另一路信号通过同轴线缆传输至发射喇叭以产生高功率微波辐射信号,高功率辐射波传输通过等离子体阵列后被接收喇叭接收,而后经衰减器、检波器、同轴线缆接入高精度示波器,作为接收信号。柱状等离子体阵列与收发喇叭口面平行,距离发射喇叭3 m(满足远场条件),紧贴接收喇叭并使等离子体单元完全覆盖整个接收喇叭口面,从而能够保证到达等离子体阵列处的高功率微波功率密度与到达接收口面的功率密度相一致。等离子体采用高频辉光放电方式产生,其初始平均电子密度约为1.2×1016 m-3
4实验测试示意图
Fig.4Schematic diagram of the experimental test
4.2 实验结果
实验采用1.3 GHz高功率微波源,脉冲宽度为1 000 ns。图5为脉冲源发射功率为170 kW时(此时辐照到等离子体上的电场强度约为8 100 V/m),等离子体未开启情况下的测试结果,其中黄色线为发射端输入的参考信号测试结果,绿色线为接收信号测试结果。图6为脉冲源发射功率约为144 kW时(此时辐照到等离子体上的电场强度约为7 450 V/m),等离子体开启情况下的测试结果。对比图5图6可以看出,等离子体对1.3 GHz高功率微波产生了明显的防护作用(透射衰减大于20 dB)。
5发射功率为170 kW、等离子体未开启情况下的测试结果
Fig.5Test results in condition of emission power 170 kW, plasma turned off
实验过程中还对同样条件下的双层等离子体进行了防护验证,高功率微波透射衰减值超过26 dB,显然,等离子体厚度越大,防护效果越好。另外,由于是初步探索实验,在设置上缺少了对极化方式、微波频率、等离子体状态等多因素的考虑,在后续的研究工作中将进一步完善。
6发射功率为144 kW、等离子体开启情况下的测试结果
Fig.6Test results in condition of emission power 144 kW, plasma turned on
5 结论
本文基于等离子体流体近似方法,建立了高功率微波与等离子体相互作用仿真计算模型,研究了高功率微波与柱状等离子体阵列相互作用过程中入射电场随时间的演变过程,分析了等离子体防护高功率微波的机理和物理过程。仿真结果表明,高功率微波通过等离子体时,会被后者吸收造成电子雪崩,产生更高电子密度的等离子体,使得等离子体对高功率微波的反射效应不断增强。这一过程不断发生,直至仅能有少量的微波能量进入等离子体,用于后产生电子的维持。此后,当进入等离子体中的微波能量继续增大时,电子雪崩效应再次发生,电子密度也会相应继续增大,直至达到新的动态平衡;而当高功率微波消失时,由于缺乏外界能量维持,电子和离子迅速复合,等离子体迅速恢复到与高功率微波相互作用前的状态。同时,利用本文建立的模型,还可以对高功率微波与等离子体相互作用过程中的影响因素(如高功率微波频率、脉宽、等离子体初始电子密度、初始电子能量、放电气体种类和气压等)进行分析。
此外,利用高频辉光放电产生柱状等离子体阵列,通过实验证实了等离子体对高功率微波具有很好的防护作用。目前正在开展基于等离子体的高功率微波小型化防护器件研究,主要解决的问题有:①对于被保护电子系统工作频带内的小功率电磁信号,防护器件不能影响其正常传输;②在高功率微波作用下,等离子体发生非线性效应的时间尽可能短,以免造成漏过微波功率过大问题;③等离子体产生非线性效应的响应频率要尽可能宽;④等离子体产生装置自身能经得起高功率微波的攻击,同时其体积、质量和功耗满足使用要求。
1高功率微波与柱状等离子体阵列相互作用的仿真计算模型示意图
Fig.1Simulation schematic diagram of the interactions between high-power microwave and columnar plasma arrays
2计算模型中电场强度空间分布随激励时间的演变过程
Fig.2Time dependence of the electric field intensity in simulation calculation model
3模型中轴线电子密度和电场分布随激励时间的演变过程
Fig.3Time dependence of the axial electron density and electric field distribution in simulation calculation model
4实验测试示意图
Fig.4Schematic diagram of the experimental test
5发射功率为170 kW、等离子体未开启情况下的测试结果
Fig.5Test results in condition of emission power 170 kW, plasma turned off
6发射功率为144 kW、等离子体开启情况下的测试结果
Fig.6Test results in condition of emission power 144 kW, plasma turned on
1氩等离子体内部碰撞反应方程及类型
Table1Collision reaction functions and types inside Ar plasma
刘振林, 杨光, 段难, 等. 高功率微波导弹对战场环境的影响及对抗技术研究[J]. 微波学报,2020,36(增刊1):358-361. LIU Z L, YANG G, DUAN N,et al. Research on the impact of CHAMP on the battlefield environment and the countermeasure technology[J]. Journal of Microwaves,2020,36(Suppl 1):358-361.(in Chinese)
冯奇, 傅镇波. 高功率微波武器典型场景应用分析[J]. 中国电子科学研究院学报,2021,16(9):916-920. FENG Q, FU Z B. Application analysis of HPM weapon in typical scenarios[J]. Journal of China Academy of Electronics and Information Technology,2021,16(9):916-920.(in Chinese)
刘洋, 程立. 电磁脉冲防护技术研究现状[J]. 材料导报,2016,30(增刊2):272-275. LIU Y, CHENG L. Research status of electromagnetic pulse weapon and its protection technology[J]. Materials Review,2016,30(Suppl 2):272-275.(in Chinese)
宋玮, 邵浩, 张治强, 等. 射频击穿等离子体对高功率微波传输特性的影响[J]. 物理学报,2014,63(6):158-162. SONG W, SHAO H, ZHANG Z Q,et al. High power microwave propagation properties in radio frequency breakdown plasma[J]. Acta Physica Sinica,2014,63(6):158-162.(in Chinese)
KOURTZANIDIS K, BOEUF J P, ROGIER F. Three dimensional simulations of pattern formation during high-pressure,freely localized microwave breakdown in air[J]. Physics of Plasmas,2014,21(12):123513.
赵朋程, 郭立新, 李慧敏.110 GHz高功率微波在大气击穿等离子体中的传输、反射和吸收[J]. 电波科学学报,2016,31(3):512-515. ZHAO P C, GUO L X, LI H M. Transmission,reflection and absorption of 110 GHz high-power microwave in air breakdown plasma[J]. Chinese Journal of Radio Science,2016,31(3):512-515.(in Chinese)
PAYNE K, XU K, CHOI J H,et al. Multiphysics analysis of plasma-based tunable absorber for high-power microwave applications[J]. IEEE Transactions on Antennas and Propagation,2021,69(11):7624-7636.
WANG H Y, HU F, HU B,et al. Characteristics of microwave breakdown in cavity filter under high power microwave environment[C]//Proceedings of 2020 IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization(NEMO),2020.
郑灵, 赵青, 罗先刚, 等. 等离子体中电磁波传输特性理论与实验研究[J]. 物理学报,2012,61(15):343-349. ZHENG L, ZHAO Q, LUO X G,et al. Theoretical and experimental studies of electromagnetic wave transmission in plasma[J]. Acta Physica Sinica,2012,61(15):343-349.(in Chinese)
李志刚, 陈宗胜. 等离子体与高功率微波相互作用中电子分布特性[J]. 国防科技大学学报,2020,42(1):10-17. LI Z G, CHEN Z S. Distributions of the electron in the interactions between high power microwave and plasma[J]. Journal of National University of Defense Technology,2020,42(1):10-17.(in Chinese)
袁忠才, 时家明. 高功率微波与等离子体相互作用理论和数值研究[J]. 物理学报,2014,63(9):255-264. YUAN Z C, SHI J M. Theoretical and numerical studies on interactions between high-power microwave and plasma[J]. Acta Physica Sinica,2014,63(9):255-264.(in Chinese)
HE W, LIU X H, XIAN R C,et al. Kinetics characteristics and bremsstrahlung of argon DC discharge under atmospheric pressure[J]. Plasma Science and Technology,2013,15(4):335-342.