下载:
-
等离子体相对论微波噪声放大器(plasma relativistic microwave noise amplifier, PRNA)是一种新型等离子体相对论微波发生器 (plasma relativistic microwave generator, PRMG)[1,2], 其工作机理为强流相对论电子束与加载具有陡峭边界的等离子体束的波导中激发的等离子体慢波发生契伦柯夫共振相互作用, 使得电子束的能量向等离子体慢波转移, 从而实现波的放大[3,4]. 与为提高器件输出性能而填充等离子体的电真空相对论微波器件不同[5-7], PRMG中的等离子体束是实现器件运行的必要条件, 辐射频率的调节可以通过改变等离子体密度来实现. 通过束波互作用区结构的设计和工作参数的选取, PRMG能以普通放大器或振荡器形式工作, 如等离子体相对论微波放大器(plasma relativistic microwave amplifier, PRMA) [8]或振荡器(plasma relativistic microwave oscillator, PRMO)[9], 也能够以噪声放大器的形式工作, 即PRNA. PRNA可以产生超宽带微波辐射, 与此同时, PRNA还具有输出功率高和可重复频率运行等潜质. 因此, PRNA在脉冲雷达、通信、电子对抗、物体探测、地质以及医学等诸多领域均具有良好的应用前景[10,11].
与传统的利用快速开关直接把直流电能转换成微波输出的高功率超宽带电磁脉冲辐射源[12] (electromagnetic pulse sources, EMPS)相比, PRNA具有以下优势: 不同脉冲的辐射频率可以在较大范围内变化, 且脉冲宽度不随中心频率的增加而缩短; 单脉冲能量可以较大, 甚至可以达到10 J[13]. 当然, 超宽带EMPS无需相对论电子束, 也就不需高能强流电子束源和真空电子设备, 原理相对简单, 成本较低. 因此, 可以根据实际需求, 合理选取超宽带微波的产生方式. PRNA的理论和实验研究以及输出性能的提升将可能为某些重要应用提供更适合的宽带微波源.
与同属PRMG器件的其他器件相比, PRMA工作时需要微波种子源的激励, PRMO调频方式为阶跃型调频, 不能在频率区间内连续调频, 阶跃宽度主要由波束互作用区具体结构参数决定. PRNA则可以实现在连续频谱内的调频, 且无需微波种子源和反馈机制, 因此结构更加紧凑. 但相对于PRMA和PRMO, PRNA的提出相对较晚. 2013年和2014年, Ernyleva等[14,15]首次提出了PRNA, 数值模拟得到了频率范围为4—17 GHz、带宽为2 GHz、功率为150 MW和效率为15%的微波输出. 随后, 科学家对PRNA进行了进一步的理论研究, 模拟得到了频率范围2—12 GHz、功率20 MW、效率4%—9% [16]和频率范围3—9 GHz、功率40 MW、效率约8%[17]的微波输出. 2019年, Strelkov 等[13]首次进行了PRNA的实验研究. 实验结果证实了PRNA产生连续谱超宽带的可行性, 并获得了单脉冲10 J的微波能量, 功率百MW量级, 输出微波的中心频率可以通过改变等离子体束密度进行调节, 频率范围为2.0—3.5 GHz. 随后, 通过进一步增加等离子体束密度, 实验观察到微波中心频率可以从3 GHz提高到25 GHz[18-20].
尽管科学家对PRNA的研究获得了令人鼓舞的实验结果, 但总体来说, 目前, 器件输出效率低于10%, 输出功率和效率均有待进一步提高. 同时, 仍然有不少实验现象需要深入分析和物理解释, 实验也遇到了一些问题需要解决和克服. 因此需要对PRNA器件开展更加细致的研究, 以深入理解其物理机制并获取更多的运行规律, 提高其输出性能, 从而更好地满足实际应用需求.
本文拟利用全电磁粒子模拟方法对PRNA进行物理分析和数值模拟, 其中等离子体束和相对论电子束中所有的粒子均将在电磁场的作用下自洽运动. 因此, 可以更多地考虑等离子体中电子和离子的运动与慢等离子体波之间的耦合关系, 从而更加细致地描述波束互作用过程. 通过对PRNA中波束互作用区进行数值模拟, 可以加深对束互作用物理机制和规律的理解和认识, 同时还可以给出主要物理参数, 尤其是等离子体密度、厚度以及相对论电子束电流和电压等对波束色散关系和器件主要输出特性的影响规律, 为进一步优化设计器件提供参考.
第2节将进行物理建模, 并对构建的PRNA模型中波束互作用区中的色散关系进行模拟计算, 同时利用已有近似解析公式[4]对不同物理参数下的线性增长率进行计算, 为后续整体模拟时参数的选择提供依据. 第3节对构建的PRNA模型进行整体数值模拟, 给出波束互作用过程的典型物理图像, 重点考查等离子体束密度np和厚度Δrp、电子束电流和电压, 以及等离子体束和电子束径向间距Δrpb等主要物理参数对器件输出微波频率和带宽的影响规律. 第4节对全文进行简要总结.
-
图1为构建的将要进行模拟和计算的PRNA结构示意图, 主要由波束互作用区和输出区组成. 由于该结构为轴对称结构, 只给出了rz截面示意图. 其中波束互作用区包括: 圆柱形金属波导1、薄环形强流相对论电子束2、薄环形等离子体束3以及粒子收集极4; 输出区是以收集极为内导体、金属波导为外导体的同轴波导. 图1中
$ {r_{\text{b}}} $ ,$ {r_{\text{p}}} $ , R0和R分别相对论电子束中心半径, 等离子体束中心半径、收集极半径和波导管半径,$ {r_{{\text{b, p}}}} < R $ , L为互作用区的长度,$ \Delta {r_{{\text{bp}}}} $ 为等离子体束和电子束径向间距,$ \Delta {r_{{\text{bp}}}}{\text{ = }}{r_{\text{p}}} - {r_{\text{b}}} $ .
图 1 PRNA结构示意图, rz截面. 1, 金属波导; 2, 相对论电子束; 3, 等离子体束; 4, 粒子收集极
Figure 1. Schematic drawing of the PRNA in the rz cross section, here 1 is the metal waveguide, 2 is the relativistic electron beam, 3 is the plasma beam and 4 is the particle collector.
模拟过程中具体参数的选取参考了相关的实验和文献[14, 15], 除非另有说明, 选取的基本参数: 互作用区半径R = 1.8 cm, 收集极半径R0 = 1.2 cm, 环形等离子体束通过低能电子束电离气压约为10–4—10–3 Torr (1 Torr = 133 Pa)惰性气体产生[21], 等离子体中组分为电子和氙离子, 等离子体束半径、径向厚度和密度分别为
$ {r_{\text{p}}}={\text{1.05 cm}} $ ,$ \Delta {r_{\text{p}}} = 0{\text{.1 cm}} $ 和np = 1.4×1019/m3, 外加磁场Bz = 2.0 T. 电子束电压和电流分别为500 kV和2 kA, 电子束半径和径向厚度分别为$ {r_{\text{b}}}{\text{ = 0}}{\text{.75 cm}} $ 和$ \Delta {r_{\text{b}}}{\text{ = 0}}{\text{.1 cm}} $ , 电子束和等离子体束径向间距为$ \Delta {r_{{\text{bp}}}} = 0{\text{.3 cm}} $ . -
利用全电磁粒子模拟程序计算了互作用区中的工作模式, 即慢等离子体波的最低模式TM01(下称P-TM01模)的色散特性. 由于相对论电子束密度相对较小, 因此参与波束互作用的电子束中的慢空间电荷波(下称慢空间电荷波)的色散曲线由解析公式[4]计算得到. 主要计算结果见图2—图4.
图 2 互作用区不同等离子体密度对应的P-TM01模与慢空间电荷波的色散关系图
Figure 2. Dispersion relations between the P-TM01 modes at different plasma density and the slow space charge wave in the beam-wave interaction region.
图 4 互作用区P-TM01模与不同电压和电流的慢空间电荷波的色散关系图
Figure 4. Dispersion relations between the P-TM01 mode and slow space charge waves with different voltage and current.
由图2可以看到, 随着等离子体密度np的增加, 相同波矢kz对应的P-TM01模的频率增加, kz越大, 增加越明显. 其他参数不变的情况下, P-TM01模与慢空间电荷波的交点对应的共振频率随着np的增加而升高.
由图3可以看到, 等离子体束的厚度
$ \Delta {r_{\text{p}}} $ 对P-TM01模的色散特性影响也较大, 相同波矢kz对应的P-TM01模的频率随$ \Delta {r_{\text{p}}} $ 的增加而增加, kz越大, 增加越明显; 其他参数不变的情况下, P-TM01模与慢相对论电子束波的交点对应的共振频率随着$ \Delta {r_{\text{p}}} $ 的增加而显著升高.
图 3 不同等离子体厚度对应的P-TM01模与慢空间电荷波的色散关系图
Figure 3. Dispersion relations between the P-TM01 modes at different plasma radial thickness and the slow space charge wave.
由图4可以看到, 电子束电流强度越高, 或者电子束电压越低, 相同波矢kz对应的慢空间电荷波的频率越低. 其他参数不变的情况下, P-TM01模与慢空间电荷波的交点对应的共振频率随着电子束电压的增高而有所降低, 随着电子束电流的增加而有所提高. 对比np和Δrp, 一定范围内, 电子束电压和电流的变化引起的共振频率的变化幅度相对较小.
综合上述研究结果: 等离子体束和电子束参数的变化均会对波束共振点的频率产生或大或小的影响, 因此实际情况下可以根据需求选取合适的物理参数.
-
辐射场的线性增长率一般定义为
$ {G_{\text{L}}} = 2{\rm Im} ({k_z}) , $ 即单位长度的小信号增益, GL反映了线性增长区波束互作用的强度. 数值求解了文献[4]给出的色散方程, 得到了GL随等离子体束和电子束参数的变化情况. 需要说明的是, 由于色散方程中做了较多的近似和假设, 且未考虑等离子体束内部对电磁场俘获效应, 随着频率和纵向波数kz的增加, GL的计算误差将会增大. 图5—图8展示主要物理参数变化时增益带宽曲线的变化情况, 除了图中标明的自变参数外, 其他参数保持缺省值不变. 由图5—图8可以看到, PRNA具有可以产生宽带微波输出的特性.
图 5 不同等离子体密度对应的线性增长率随频率的变化曲线
Figure 5. Variations of the linear growth rate with frequency for different plasma beam density.
图 8 不同等离子体束和电子束径向距离对应的线性增长率随频率的变化曲线
Figure 8. Variations of the linear growth rate with frequency for different plasma-beam radial gap.
由图5还可以看到, 线性增长率的最大值的位置基本对应共振频率, 随着等离子体密度np的增大, 被激发的辐射场频谱范围增加, 线性增长率的最大值对应的频率相应提高, 与色散关系图2所得结论一致; 同时, 在np为1.0×1019—1.8×1019/m3的区间内, 束等离子体波相互作用均以单粒子效应(single-particle Cherenkov resonance)[21]为主, 因此线性增长率的最大值随着np的增大而增加, 产生的微波带宽下限几乎为0, 上限超过了共振频率. 随着np的进一步增大(如np = 3.0×1019/m3和4.8×1019/m3), 共振频率增加, 束等离子体波相互作用由单粒子效应为主转变为集体效应(the collective Cherenkov resonance)[21]为主, 因而产生的微波带宽收缩到共振频率附近. 线性增长率的最大值也逐渐减小.
由图6还可以看到, 随着等离子体厚度
$\Delta {r_{\text{p}}}$ 的增加, 线性增长率最大值对应的频率显著提高, 与色散关系图3所得结论一致. 一定范围内, 随着$\Delta {r_{\text{p}}}$ 的增加, 被激发的辐射场频谱范围显著增加; 但随着$\Delta {r_{\text{p}}}$ 进一步增加, 辐射场频谱范围区域发生较大变化. 这是由于共振频率随着$\Delta {r_{\text{p}}}$ 的增加而显著升高, 束等离子体波相互作用由单粒子效应为主转变为集体效应为主, 因而产生的微波带宽收缩到共振频率附近.
图 6 不同等离子体厚度对应的线性增长率随频率的变化曲线
Figure 6. Variations of the linear growth rate with frequency for different plasma radial thickness.
如果电子束电压增加到600 kV, 对于相同的径向厚度Δrp = 0.3 cm时, 由于共振频率的降低(参见图3电子束电压增加到600 kV时的色散曲线), 束等离子体波相互作用由集体效应为主转回到单粒子效应为主. 此时, 产生的微波带宽下限几乎为0, 上限超过共振频率(见图6电子束电压为600 kV时的增益带宽曲线).
上述计算结果表明, 为了获得超宽带输出, 在PRNA设计时需要合理地选取相对论电子束参数和等离子体束参数, 如果需要覆盖低频段, 则需要使器件工作在单粒子效应为主的参数区间.
由图7不难看出, 随着电子束电压增加, 被激发的辐射场频谱范围略有减小, 线性增长率的最大值降低, 对应频率有所降低. 随着电子束电流增加, 辐射场频谱范围增加, 线性增长率的最大值增加, 对应频率有所提高, 频率的变化情况与色散关系图3所得结论一致.
图 7 具有不同电压和电流的相对论电子束对应的线性增长率随频率的变化曲线
Figure 7. Variations of the linear growth rate with frequency for different electron beam voltage and current.
由于我们是通过改变电子束半径
$ {r_{\text{b}}} $ 来改变$ \Delta {r_{{\text{bp}}}} $ 的值, 等离子体束中心半径$ {r_{\text{p}}} $ 保持不变, 因此$ \Delta {r_{{\text{bp}}}} $ 的改变对等离子体慢波P-TM01的色散特性没有影响. 由图8可以看到, 当$\Delta {r_{{\text{bp}}}}$ 减小时, 被激发的辐射场频谱范围基本保持不变. 但随着$\Delta {r_{{\rm{bp}}}}$ 的减小, 电子束逐渐靠近等离子体束, 参与波束互作用的P-TM01的纵向电场增大, 波束相互作用增强, 线性增长率有所增加.以上计算和分析表明: 1) PRNA具有产生宽带微波输出的优势; 2) 通过改变等离子体束的密度和厚度、电子束电压和电流等相关物理量, 可以方便地调节束波互作用共振点对应的频率. 因此, PRNA同时具有良好的可调谐性.
-
为了证实PRNA在带宽和可调谐性方面的输出优势, 下面将利用全电磁粒子模拟程序对其进行整体数值模拟. 首先选取相对论电子束脉冲的形状为梯形脉冲, 电子束脉宽小于电磁波在互作用区往返传输一次所需的时间, 以确保互作用区端口可能的微波反射不参与波束互作用, 这里取电子束总脉宽为2 ns, 上升沿和下降沿各为0.5 ns, 束波互作用区长度L = 35 cm, 除非另有说明, 其他参数同2.1节给的缺省值.
下面给出主要的模拟结果, 包括粒子在实空间和相空间的相位图、不同时刻互作用区和输出区径向电场和轴向电场的等高图、输出功率时间演化图及其频谱图等. 同时将给出主要物理参数的变化对输出性能的影响情况, 并与前面分析结果进行比较.
-
图9和图10为不同时刻等离子体束中电子、离子以及电子束电子在实空间和相空间中的相位图. 可以看到, 在时间为1.689 ns时, 相对论电子束的动量已经出现明显的调制效果, 到2.295 ns时动量调制进一步加剧, 电子束达到了强群聚状态, 电子束密度出现明显群聚现象, 此时波束互作用程度加剧, 产生的辐射场强度将大幅增加. 等离子体中的电子动量也随着辐射场的增加出现振荡, 与相对论电子束动量振荡频率一致. 由于离子质量远大于电子质量, 因此离子状态几乎没有发生变化.
图 9 不同时刻等离子体电子(红色)和离子(蓝色)、电子束电子(橙色)在实空间中的分布图 (a) t = 1.689 ns; (b) t = 2.295 ns
Figure 9. Real space plot of the plasma electrons (red) and ions (blue), and beam electrons (orange) at different times: (a) t = 1.689 ns; (b) t = 2.295 ns.
图 10 不同时刻等离子体电子(红色)和离子(蓝色)、电子束电子(橙色)在相空间中的分布图 (a) t = 1.689 ns; (b) t = 2.295 ns
Figure 10. Phase space plot of the plasma electrons (red) and ions (blue), and beam electrons (orange) at different times: (a) t = 1.689 ns; (b) t = 2.295 ns.
图11—图13分别为不同时刻波束互作用区和输出区轴向和径向电场的等高图. 通过图11—图13, 可以看到辐射场产生、发展和耦合输出的过程. 其中, 由图11和图12可知, 互作用区电磁场分布与P-TM01模的本征场分布完全一致; 由图12和图13可以看到, 互作用区和输出区的径向电场基本自然过渡, 圆柱波导中的P-TM01模转化为同轴TEM模输出.
图 11 不同时刻互作用区纵向截面轴向电场的等高图 (a) t = 1.689 ns; (b) t = 1.916 ns
Figure 11. Contour plot of the electric field Ez in the interaction region at different times: (a) t = 1.689 ns; (b) t = 1.916 ns.
图 12 不同时刻互作用区纵向截面径向电场的等高图 (a) t = 1.689 ns; (b) t = 1.916 ns
Figure 12. Contour plot of the electric field Er in the interaction region at different times: (a) t = 1.689 ns; (b) t = 1.916 ns.
图 13 不同时刻输出区纵向截面径向电场的等高图 (a) t = 2.295 ns; (b) t = 2.975 ns
Figure 13. Contour plot of the electric field Er in the output region at different times: (a) t = 2.295 ns; (b) t = 2.975 ns.
图14和图15给出了瞬时输出功率和平均功率随时间的变化曲线及其FFT变化图, 数值证实了PRNA作为一种高功率宽带器件的可行性. 从图14和图15可知, 瞬时输出功率接近600 MW, 平均峰值功率达到200 MW, 功率效率达到了20%. 图16给出了输出区辐射场Er随时间的变化曲线及其FFT变化图, 可以明显观察到, 辐射场频率范围约从7.0 GHz到大于9.0 GHz, 带宽达到了2 GHz左右.
图 14 瞬时输出功率随时间的变化曲线(a)及其FFT变换图(b)
Figure 14. Time plots of (a) the instantaneous and (b) the corresponding Fourier transform.
图 15 周期平均输出功率随时间的变化曲线
Figure 15. Time plots of the periodic-average output power.
图 16 输出区辐射场Er (a) 随时间的变化曲线及其(b) FFT变换图
Figure 16. Variations of Er (a) with time and (b) the corresponding Fourier transform in the output region.
-
本节将通过改变等离子体密度和厚度、相对论电子束电压和电流以及等离子体束和电子束径向间距等主要物理参数的取值范围, 对PRNA进行了全电磁粒子模拟并对模拟结果进行总结和分析, 以获得上述参数的变化对PRNA输出性能的影响规律, 重点考虑对带宽和频率的影响.
图17和图18分别给出了输出微波频谱与效率随等离子体密度的变化图. 从图17和图18可以看到, 随着等离子体密度的增加, 输出微波频率呈上升趋势, 中心频率约从7.0 GHz提高到10 GHz, 变化趋势与前面无耦合时色散关系的模拟结果一致. 保持其他参数不变的情况下, 尤其是波束互作用区长度固定时, 在一定等离子体密度范围内, 微波输出效率保持大于15%.
图 17 输出微波频谱随等离子体密度的变化
Figure 17. Variations of the output frequency spectrum with different plasma density.
图 18 输出效率随等离子体密度的变化
Figure 18. Variations of the output efficiency with different plasma density.
图19给出了输出微波频谱随等离子体厚度的变化图. 从图19可以看到, 随着等离子体厚度的增加, 输出微波频率显著提高, 等离子体束径向厚度从0.1 cm增加到0.3 cm时, 中心频率约从8.3 GHz提高到17.5 GHz, 变化趋势与前面无耦合时色散关系的模拟结果一致. 由于波束互作用区长度L固定, 输出效率随等离子体束厚度的增加而明显减小. 但通过对L的优化, 输出效率仍然能够达到约20%.
图 19 微波频谱随等离子体束厚度的变化
Figure 19. Variations of the output frequency spectrum with different plasma radial width.
图20和图21给出了输出微波频谱随电子束电压和电流的变化图. 从图20和图21可以看到, 输出频率随着电子束电压的增加而降低, 随着电子束电流的增加而升高, 与前面的冷腔分析一致. 但与等离子体密度和厚度相比, 一定范围内, 电流和电压的变化对输出频率的影响比较小, 因此只能作为频率调节的辅助手段. 对于给定的电子束电流和电压, 由于波束互作用区长度L存在最优值, 因此, 当L固定时, 输出效率随电子束电压和电流的变化而减小. 但通过对L的优化, 一定变化范围内, 输出效率能够达到约20%.
图 20 微波频谱随电子束电压的变化
Figure 20. Variations of the output frequency spectrum with different beam voltage.
图 21 输出微波带宽随电子束电流的变化
Figure 21. Variations of the output frequency spectrum with different beam current.
模拟结果表明, 一定范围内径向间距
$\Delta {r_{{\text{bp}}}}$ 对输出微波频率的影响很小, 中心频率基本保持不变, 但对输出效率有一定影响. 图22给出了输出微波效率随$\Delta {r_{{\text{bp}}}}$ 的变化情况, 可以看到, 尽管没有随着$\Delta {r_{{\text{bp}}}}$ 的变化对L进行优化, 但输出效率均超过了15%.
图 22 输出效率随等离子体束和电子束径向间距的变化
Figure 22. Variations of the output efficiency with different plasma and electron beam gaps.
-
利用全电磁粒子模拟方法对PRNA进行了物理分析和数值模拟. 首先利用全电磁粒子模拟方法对无耦合时的波束色散关系进行了模拟分析. 结果表明PRNA具有良好的可调谐性, 通过改变等离子体束的密度和厚度、电子束电压或电流等相关物理量, 可以方便地调节束波互作用共振点对应的频率. 接着利用近似解析公式得到了微波线性增长率与带宽的变化规律. 在上述模拟和分析的基础上, 对PRNA进行了整体模拟, 验证了PRNA在带宽和可调谐性方面的输出优势, 得到了相关物理参数对器件主要输出特性的影响规律. 当等离子体密度为1.4×1019/m3, 外加磁场为2.0 T时, 电子束电压和电流分别为500 kV和2 kA, 模拟获得了功率约200 MW, 效率为20%的微波输出, 辐射场频谱范围约为7.0—9.0 GHz, 带宽达到了2 GHz左右, 输出模式为同轴TEM模. 模拟结果还表明, 等离子体束的密度np和厚度
$ \Delta {r_{\text{p}}} $ 对束波色散关系影响较大, 随着np和$ \Delta {r_{\text{p}}} $ 的增加, 输出微波频率呈明显上升趋势. 当等离子体束径向厚度从0.1 cm 增加到0.3 cm时, 中心频率约从8.3 GHz提高到17.5 GHz. 电子束电流和电压的变化对输出频率的影响相对较小, 因此改变电子束电流和电压只能作为频率调节的辅助手段. 等离子体束和电子束径向间距的变化则对输出频率基本没有影响. 以上研究结果可为器件的进一步优化设计提供参考.
-
利用全电磁粒子模拟方法对等离子体相对论微波噪声放大器(plasma relativistic microwave noise amplifier, PRNA)进行了物理分析和数值模拟. 首先对无耦合时的波束色散关系进行了模拟分析, 接着计算了微波线性增长率与带宽的变化规律, 为后续整体模拟时参数的选择提供了理论依据. 最后对PRNA进行了整体模拟, 验证了PRNA在带宽和可调谐性方面的输出优势. 当等离子体束密度为1.4×1019/m3, 电子束电压和电流分别为500 kV和2 kA, 外加磁场为2.0 T时, 模拟获得了功率约200 MW, 效率为20%的微波输出, 辐射场频谱范围约为7.0—9.0 GHz, 带宽达到了2 GHz, 输出模式为同轴TEM模. 模拟结果还表明: 等离子体束的密度np和厚度
$ \Delta {r_{\text{p}}} $ 对束波色散关系影响较大, 随着np和$ \Delta {r_{\text{p}}} $ 的增加, 输出微波频率呈明显上升趋势, 电子束电流和电压的变化对输出频率的影响相对较小, 等离子体束和电子束径向间距的变化则对输出频率基本没有影响. 研究结果可为器件的进一步的优化设计提供参考依据.-
关键词:
- 相对论微波噪声放大器 /
- 等离子体 /
- 色散关系 /
- 粒子模拟
The physical mechanism and output properties of the plasma relativistic microwave noise amplifier (PRNA) are studied numerically by using the all electromagnetic particle-in-cell (PIC) code. Firstly, the dispersion relation between the operating mode and the slow space charge wave of relativistic electron beam without coupling is simulated and analyzed. Simulation results show that both the plasma density np and radial thickness Δrp affect the dispersion characteristics markedly and their increasing can lead the frequency at the beam-wave resonant point to be enhanced. The beam voltage and current also affect the resonant frequency, but the effect is relatively slight. Secondly, variation of the linear growth rate and the bandwidth are then evaluated by using the linear theory. Calculations show that the PRNA has the virtue of wideband output. Its bandwidth can reach a GHz level. By adjusting the plasma parameters np and Δrp, the relativistic electron beam voltage and current, the operating frequency can be tuned over a wide frequency range. Therefore the PRNA also has virtue of fine frequency tunability. Based on the above calculation results, the whole PIC simulations of the PRNA are then carried out to verify the virtues of wideband microwave output and frequency tunability. The basic features of the field distributions of the operating in the evolution process and out coupling process are given. The bunching process and the energy release process of relativistic electron beam are also plotted. Simulations show that with a plasma density of 1.4×1019 /m3, beam voltage and current of 500 kV and 2 kA and applied magnetic field of 2.0 T, 200 MW output microwave with efficiency about 20% can be obtained. The frequency ranges from about 7.0 to 9.0 GHz, the band width reaches 2 GHz. And the output mode is the TEM mode of the coaxial waveguide. Both np and Δrp affect the dispersion relations markedly and the output frequency increases clearly with np and Δrp increasing. The influence of beam voltage and current on the output frequency are both relatively small and the gap distance between the plasma and electron beam has little effect on the output frequency. The research results will provide useful reference for further designing the PRNA.-
Keywords:
- relativistic microwave noise amplifier /
- plasma /
- dispersion relation /
- numerical simulation
[1] Kuzelev M V, Mukhametzyanov F Kh, Rabinovich M S, Rukhadze A A, Strelkov P S, Shkvarunets A G 1982 Sov. Phys. JETP 56 780
[2] Kuzelev M V, Loza O T, Ponomarev A V, Rukhadze A A, Strelkov P S, Shkvarunets A G, Ulyanov D K 1996 Sov. Phys. JETP 82 1102
[3] Bogdankevich L S, Kuzelev M V, Rukhadze A A 1981 Sov. Phys. Usp. 24 1
[4] Kuzelev M V, Rukhadze A A 2000 Plasma Phys. Rep. 26 231
Google Scholar
[5] Carmel Y, Lou W R, Antonsen Jr T M, Rodgers J, Levush B, Destler W W, Granatstein V L 1992 Phys. Fluids B 4 2286
Google Scholar
[6] Shkvarunets A G, Kobayashi S, Weaver J, Carmel Y, Rodgers J, Antonsen Jr T M, Granatstein V L, Destler W W, Ogura K, Minami K 1996 Phys. Rev. E 53 2045
Google Scholar
[7] 王宇, 陈再高, 雷奕安 2013 物理学报 62 125204
Google Scholar
Wang Y, Chen Z G, Lei Y A 2013 Acta Phys. Sin. 62 125204
Google Scholar
[8] Ponomarev A V, Strelkov P S, Shkvarunets A G 1998 Plasma Phys. Rep. 24 48
[9] Strelkov P S, Ul’yanov D K 2000 Plasma Phys. Rep. 26 303
Google Scholar
[10] Barker R J, Schamiloglu E 2001 High-Power Microwaves Sources and Technologies (New York: Institute of Electrical and Electronic Engineer, Inc.) pp25–27
[11] Prather W D, Baum C E, Torres R J, Sabath F, Nitsch D 2004 IEEE Trans. Electromagn. Compat. 46 335
Google Scholar
[12] 黄裕年, 任国光 2002 微波学报 18 90
Google Scholar
Huan Y N, Ren G G 2002 J. Microwaves 18 90
Google Scholar
[13] Strelkov P S, Tarakanov V P, Mikh D E D, Ivanov I E, Shumeiko D V 2019 Plasma Phys. Rep. 45 345
Google Scholar
[14] Ernyleva S E, Litvin V O, Loza O T, Bogdankevich I L 2013 Prob. At. Sci. Technol. 2013 3
[15] Ernyleva S E, Litvin V O, Loza O T, Bogdankevich I L 2014 Tech. Phys. 59 1228
Google Scholar
[16] Ernyleva S E, Loza O T 2017 Phys. Wave Phenom. 25 56
Google Scholar
[17] Ulyanov D K, Bogdankevich I L, Ernyleva S E, Andreev S E 2019 Plasma Phys. Rep. 45 980
Google Scholar
[18] Kartashov I N, Kuzelev M V 2021 Plasma Phys. Rep. 47 548
Google Scholar
[19] Buleyko A B, Ponomarev A V, Loza O T, Ulyanov D K, Andreev S E 2021 Phys. Plasmas 28 023303
Google Scholar
[20] Buleyko A B, Ponomarev A V, Loza O T, Ulyanov D K, Sharypov K A, Shunailov S A, Yalandin M I 2021 Phys. Plasmas 28 023304
Google Scholar
[21] Kuzelev M V, Loza O T, Rukhadze A A, Strelkov P S, Shkvarunets A G 2001 Plasma Phys. Rep. 27 669
Google Scholar
-
图 1 PRNA结构示意图, rz截面. 1, 金属波导; 2, 相对论电子束; 3, 等离子体束; 4, 粒子收集极
Fig. 1. Schematic drawing of the PRNA in the rz cross section, here 1 is the metal waveguide, 2 is the relativistic electron beam, 3 is the plasma beam and 4 is the particle collector.
图 2 互作用区不同等离子体密度对应的P-TM01模与慢空间电荷波的色散关系图
Fig. 2. Dispersion relations between the P-TM01 modes at different plasma density and the slow space charge wave in the beam-wave interaction region.
图 4 互作用区P-TM01模与不同电压和电流的慢空间电荷波的色散关系图
Fig. 4. Dispersion relations between the P-TM01 mode and slow space charge waves with different voltage and current.
图 3 不同等离子体厚度对应的P-TM01模与慢空间电荷波的色散关系图
Fig. 3. Dispersion relations between the P-TM01 modes at different plasma radial thickness and the slow space charge wave.
图 5 不同等离子体密度对应的线性增长率随频率的变化曲线
Fig. 5. Variations of the linear growth rate with frequency for different plasma beam density.
图 8 不同等离子体束和电子束径向距离对应的线性增长率随频率的变化曲线
Fig. 8. Variations of the linear growth rate with frequency for different plasma-beam radial gap.
图 6 不同等离子体厚度对应的线性增长率随频率的变化曲线
Fig. 6. Variations of the linear growth rate with frequency for different plasma radial thickness.
图 7 具有不同电压和电流的相对论电子束对应的线性增长率随频率的变化曲线
Fig. 7. Variations of the linear growth rate with frequency for different electron beam voltage and current.
图 9 不同时刻等离子体电子(红色)和离子(蓝色)、电子束电子(橙色)在实空间中的分布图 (a) t = 1.689 ns; (b) t = 2.295 ns
Fig. 9. Real space plot of the plasma electrons (red) and ions (blue), and beam electrons (orange) at different times: (a) t = 1.689 ns; (b) t = 2.295 ns.
图 10 不同时刻等离子体电子(红色)和离子(蓝色)、电子束电子(橙色)在相空间中的分布图 (a) t = 1.689 ns; (b) t = 2.295 ns
Fig. 10. Phase space plot of the plasma electrons (red) and ions (blue), and beam electrons (orange) at different times: (a) t = 1.689 ns; (b) t = 2.295 ns.
图 11 不同时刻互作用区纵向截面轴向电场的等高图 (a) t = 1.689 ns; (b) t = 1.916 ns
Fig. 11. Contour plot of the electric field Ez in the interaction region at different times: (a) t = 1.689 ns; (b) t = 1.916 ns.
图 12 不同时刻互作用区纵向截面径向电场的等高图 (a) t = 1.689 ns; (b) t = 1.916 ns
Fig. 12. Contour plot of the electric field Er in the interaction region at different times: (a) t = 1.689 ns; (b) t = 1.916 ns.
图 13 不同时刻输出区纵向截面径向电场的等高图 (a) t = 2.295 ns; (b) t = 2.975 ns
Fig. 13. Contour plot of the electric field Er in the output region at different times: (a) t = 2.295 ns; (b) t = 2.975 ns.
图 14 瞬时输出功率随时间的变化曲线(a)及其FFT变换图(b)
Fig. 14. Time plots of (a) the instantaneous and (b) the corresponding Fourier transform.
图 16 输出区辐射场Er (a) 随时间的变化曲线及其(b) FFT变换图
Fig. 16. Variations of Er (a) with time and (b) the corresponding Fourier transform in the output region.
图 17 输出微波频谱随等离子体密度的变化
Fig. 17. Variations of the output frequency spectrum with different plasma density.
图 18 输出效率随等离子体密度的变化
Fig. 18. Variations of the output efficiency with different plasma density.
图 19 微波频谱随等离子体束厚度的变化
Fig. 19. Variations of the output frequency spectrum with different plasma radial width.
图 20 微波频谱随电子束电压的变化
Fig. 20. Variations of the output frequency spectrum with different beam voltage.
图 21 输出微波带宽随电子束电流的变化
Fig. 21. Variations of the output frequency spectrum with different beam current.
-
[1] Kuzelev M V, Mukhametzyanov F Kh, Rabinovich M S, Rukhadze A A, Strelkov P S, Shkvarunets A G 1982 Sov. Phys. JETP 56 780
[2] Kuzelev M V, Loza O T, Ponomarev A V, Rukhadze A A, Strelkov P S, Shkvarunets A G, Ulyanov D K 1996 Sov. Phys. JETP 82 1102
[3] Bogdankevich L S, Kuzelev M V, Rukhadze A A 1981 Sov. Phys. Usp. 24 1
[4] Kuzelev M V, Rukhadze A A 2000 Plasma Phys. Rep. 26 231
Google Scholar
[5] Carmel Y, Lou W R, Antonsen Jr T M, Rodgers J, Levush B, Destler W W, Granatstein V L 1992 Phys. Fluids B 4 2286
Google Scholar
[6] Shkvarunets A G, Kobayashi S, Weaver J, Carmel Y, Rodgers J, Antonsen Jr T M, Granatstein V L, Destler W W, Ogura K, Minami K 1996 Phys. Rev. E 53 2045
Google Scholar
[7] 王宇, 陈再高, 雷奕安 2013 物理学报 62 125204
Google Scholar
Wang Y, Chen Z G, Lei Y A 2013 Acta Phys. Sin. 62 125204
Google Scholar
[8] Ponomarev A V, Strelkov P S, Shkvarunets A G 1998 Plasma Phys. Rep. 24 48
[9] Strelkov P S, Ul’yanov D K 2000 Plasma Phys. Rep. 26 303
Google Scholar
[10] Barker R J, Schamiloglu E 2001 High-Power Microwaves Sources and Technologies (New York: Institute of Electrical and Electronic Engineer, Inc.) pp25–27
[11] Prather W D, Baum C E, Torres R J, Sabath F, Nitsch D 2004 IEEE Trans. Electromagn. Compat. 46 335
Google Scholar
[12] 黄裕年, 任国光 2002 微波学报 18 90
Google Scholar
Huan Y N, Ren G G 2002 J. Microwaves 18 90
Google Scholar
[13] Strelkov P S, Tarakanov V P, Mikh D E D, Ivanov I E, Shumeiko D V 2019 Plasma Phys. Rep. 45 345
Google Scholar
[14] Ernyleva S E, Litvin V O, Loza O T, Bogdankevich I L 2013 Prob. At. Sci. Technol. 2013 3
[15] Ernyleva S E, Litvin V O, Loza O T, Bogdankevich I L 2014 Tech. Phys. 59 1228
Google Scholar
[16] Ernyleva S E, Loza O T 2017 Phys. Wave Phenom. 25 56
Google Scholar
[17] Ulyanov D K, Bogdankevich I L, Ernyleva S E, Andreev S E 2019 Plasma Phys. Rep. 45 980
Google Scholar
[18] Kartashov I N, Kuzelev M V 2021 Plasma Phys. Rep. 47 548
Google Scholar
[19] Buleyko A B, Ponomarev A V, Loza O T, Ulyanov D K, Andreev S E 2021 Phys. Plasmas 28 023303
Google Scholar
[20] Buleyko A B, Ponomarev A V, Loza O T, Ulyanov D K, Sharypov K A, Shunailov S A, Yalandin M I 2021 Phys. Plasmas 28 023304
Google Scholar
[21] Kuzelev M V, Loza O T, Rukhadze A A, Strelkov P S, Shkvarunets A G 2001 Plasma Phys. Rep. 27 669
Google Scholar
计量
- 文章访问数: 434
- PDF下载量: 21
- 被引次数: 0
