搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

螺旋波等离子体中螺旋波与Trivelpiece-Gould波模式耦合及线性能量沉积特性参量分析

李文秋 赵斌 王刚 相东

引用本文:
Citation:

螺旋波等离子体中螺旋波与Trivelpiece-Gould波模式耦合及线性能量沉积特性参量分析

李文秋, 赵斌, 王刚, 相东

Parametric analysis of mode coupling and liner energy deposition properties of helicon and Trivelpiece-Gould waves in helicon plasma

Li Wen-Qiu, Zhao Bin, Wang Gang, Xiang Dong
PDF
HTML
导出引用
  • 采用有限温度等离子体介电张量模型, 在考虑粒子热效应情形下, 通过求解传导边界条件下等离子体柱中本征模的色散关系, 分析了螺旋波等离子体中典型参量条件下螺旋波与Trivelpiece-Gould (TG) 波的耦合特性及线性能量沉积特性. 在ω/(2π) = 13.56 MHz和TeV,i = 0.1TeV,e参量条件下计算结果表明: 对于螺旋波, 存在截止静磁场B0,H,cutoff与截止等离子体密度n0,H,cutoff, 在B0 > B0,H,cutoffn0 < n0,H,cutoff条件下, 螺旋波变为消逝波; 在ω/ωce ∈ (0.01, 0.10)范围内, 对于m = 0 角向模, TG波Landau阻尼致使的能量沉积占主导地位, 而对于m = 1角向模, 螺旋波Landau阻尼或TG波Landau阻尼致使的能量沉积哪个占据主导地位则取决于B0的大小; 在ωpe/ωce ∈ (3, 100)范围内, TG波Landau阻尼致使的能量沉积占主导地位; 在整体能量沉积过程中, 对于m = 0模和m = 1模, Landau阻尼致使的能量沉积均占据主导地位.
    Based on the finite temperature plasma dielectric tensor model which contains the particle thermal effect, by numerically solving the eigenmode dispersion relation of electromagnetic waves propagating in radially uniform and magnetized warm plasma column which is surrounded by conducting boundary, the mode coupling characteristic and liner damping mechanism induced wave power deposition properties of helicon and Trivelpiece-Gould (TG) waves are parametrically analyzed. The detailed investigations show as follows. Under typical helicon plasma parameter conditions, i.e. wave frequency ω/(2π) = 13.56 MHz, ion temperature is much smaller than electron temperature, for the helicon wave, there exist a cut-off magnetic field B0,H,cutoff and a cut-off plasma density n0,H,cutoff, for which under the conditions of B0 > B0,H,cutoff or n0 < n0,H,cutoff, the helicon wave becomes an evanescent wave. When the magnetic field intensity changes from 48.4 to 484 G, i.e., ω/ωce ranges from 0.01 to 0.1, for the power deposition intensity, Landau damping of TG wave dominates for the m = 0 mode, meanwhile, for the m = 1 mode, which wave, i.e. helicon wave or TG wave, plays a major role in power deposition mainly depends on the magnitude of the magnetic field. On the other hand, for a given magnetic field B0 = 100 G, when ωpe/ωce changes from 3 to 100, for both the m = 0 mode and the m = 1 mode, the power deposition induced by Landau damping of TG wave plays a major role, further, one may notice that the power deposition of TG wave decreases while the power deposition of the helicon wave increases as plasma density increases. Finally, for both the m = 0 mode and the m = 1 mode, the power deposition due to the Landau damping plays a dominant role. All these conclusions provide us with some useful clues to better understanding the high ionization mechanism of helicon wave discharges.
      通信作者: 李文秋, beiste@163.com
    • 基金项目: 国家级-国家“万人计划”科技创新领军人才(Y8BF130272)
      Corresponding author: Li Wen-Qiu, beiste@163.com
    [1]

    Chen F F 2015 Plasma Sources Sci. Technol. 24 014001Google Scholar

    [2]

    Isayama S, Shinohara S, Hada T 2018 Plasma Fusion Res. 13 1101014Google Scholar

    [3]

    Shinohara S 2018 Adv. Phys.: X 3 1420424Google Scholar

    [4]

    Aigrain P 1960 Proceedings of the International Conference on Semiconductor Physics Prague, Czech Republic, August 4–8, 1960 p224

    [5]

    Chen F F 1991 Plasma Phys. Controlled Fusion 33 339Google Scholar

    [6]

    Shamrai K P, Taranov V B 1996 Plasma Sources Sci. Technol. 5 474Google Scholar

    [7]

    Shamrai K P, Taranov V B 1994 Plasma Phys. Controlled Fusion 36 1719Google Scholar

    [8]

    Shamrai K P 1998 Plasma Sources Sci. Technol. 7 499Google Scholar

    [9]

    Chen F F, Blackwell D D 1999 Phys. Rev. Lett. 82 2677Google Scholar

    [10]

    Blackwell D D, Chen F F 2001 Plasma Sources Sci. Technol. 10 226Google Scholar

    [11]

    Kline J, Scime E 2003 Phys. Plasmas 10 135Google Scholar

    [12]

    Kim S, Hwang Y 2008 Plasma Phys. Controlled Fusion 50 035007Google Scholar

    [13]

    Isayama S, Hada T, Shinohara S, et al. 2016 Phys. Plasmas 23 063513Google Scholar

    [14]

    成玉国, 程谋森, 王墨戈等 2014 物理学报 63 035203Google Scholar

    Cheng Y G, Cheng M S, Wang M G, et al. 2014 Acta Phys. Sin. 63 035203Google Scholar

    [15]

    平兰兰, 张新军, 杨桦等 2019 物理学报 68 205201Google Scholar

    Ping L L, Zhang X J, Yang H, et al. 2019 Acta Phys.Sin. 68 205201Google Scholar

    [16]

    Chen F F, Arnush D 1997 Phys. Plasmas 4 3411Google Scholar

    [17]

    Arnush D, Chen F F 1998 Phys. Plasmas 5 1239Google Scholar

    [18]

    Arnush D 2000 Phys. Plasmas 7 3042Google Scholar

    [19]

    Sakawa Y, Kunimatsu H, Kikuchi H, et al. 2003 Phys. Rev. Lett. 90 105001Google Scholar

    [20]

    Vey B 1984 PFC/RR-84-12 (Cambridge: Massachusetts Institute of Technology) pp30–31

    [21]

    Niemi K, Krämer M 2008 Phys. Plasmas 15 073503Google Scholar

    [22]

    Shamrai K P, Shinohara S 2001 Phys. Plasmas 8 4659Google Scholar

    [23]

    Huba J D 2016 NRL Plasma Formulary (Washington: Naval Research Laboratory) p34

    [24]

    Fried B D, Conte S D 2015 The Plasma Dispersion Function: the Hilbert Transform of the Gaussian (New York: Academic Press) pp1–3

    [25]

    Mouzouris Y, Scharer J E 1998 Phys. Plasmas 5 4253Google Scholar

    [26]

    Shoji T, Sakawa Y, Nakazawa S, et al. 1993 Plasma Sources Sci. Technol. 2 5Google Scholar

    [27]

    Swanson D G 1989 Plasma Waves (New York: Academic Press) pp375–376

    [28]

    Lafleur T, Charles C, Boswell R 2010 Phys. Plasmas 17 073508Google Scholar

  • 图 1  被传导边界包裹的等离子体柱横向截面示意图

    Fig. 1.  Cross section of plasma column surround by conducting boundary

    图 2  粒子温度对螺旋波与TG波耦合色散关系的影响 (a)电子温度的影响; (b)离子温度的影响

    Fig. 2.  Influence of particle temperature on dispersion relation between helicon and TG waves: (a) Electron temperature effect; (b) ion temperature effect.

    图 3  螺旋波与TG波横向波数对轴向静磁场的依赖关系

    Fig. 3.  The perpendicular wave number of helicon and TG waves given as functions of axial static magnetic field.

    图 4  螺旋波与TG波横向波数对等离子体密度的依赖关系

    Fig. 4.  The perpendicular wave number of helicon and TG waves given as functions of plasma density.

    图 5  螺旋波轴向波数随参量变化情况 (a)轴向波数随轴向静磁场变化; (b)轴向静磁场随等离子体密度变化

    Fig. 5.  The axial wave number of the right hand polarized wave is given as a function of (a) axial static magnetic field and (b) plasma density.

    图 6  螺旋波与TG波径向功率沉积分布 (a) m = 0 角向对称模; (b) m = 1 角向对称模

    Fig. 6.  Radial power deposition profiles of the helicon and TG waves for: (a) m = 0 mode; (b) m = 1 mode.

    图 7  螺旋波与TG波功率沉积随电子温度的变化关系 (a) m = 0 模; (b) m = 1 模

    Fig. 7.  Power deposition profiles of helicon and TG waves are given as functions of electron temperature for: (a) m = 0 mode; (b) m = 1 mode.

    图 8  螺旋波与TG波功率沉积随离子温度的变化关系 (a) m = 0 模; (b) m = 1 模

    Fig. 8.  Power deposition profiles of helicon and TG waves are given as functions of ion temperature for: (a) m = 0 mode; (b) m = 1 mode.

    图 9  螺旋波与TG波功率沉积随轴向静磁场的变化关系 (a) m = 0 模; (b) m = 1 模

    Fig. 9.  Power deposition profiles of helicon and TG waves are given as functions of axial static magnetic field for: (a) m = 0 mode; (b) m = 1 mode.

    图 10  螺旋波与TG波功率沉积随等离子体密度的变化关系 (a) m = 0 模; (b) m = 1 模

    Fig. 10.  Power deposition profiles of helicon and TG waves are given as functions of plasma density for: (a) m = 0 mode; (b) m = 1 mode.

  • [1]

    Chen F F 2015 Plasma Sources Sci. Technol. 24 014001Google Scholar

    [2]

    Isayama S, Shinohara S, Hada T 2018 Plasma Fusion Res. 13 1101014Google Scholar

    [3]

    Shinohara S 2018 Adv. Phys.: X 3 1420424Google Scholar

    [4]

    Aigrain P 1960 Proceedings of the International Conference on Semiconductor Physics Prague, Czech Republic, August 4–8, 1960 p224

    [5]

    Chen F F 1991 Plasma Phys. Controlled Fusion 33 339Google Scholar

    [6]

    Shamrai K P, Taranov V B 1996 Plasma Sources Sci. Technol. 5 474Google Scholar

    [7]

    Shamrai K P, Taranov V B 1994 Plasma Phys. Controlled Fusion 36 1719Google Scholar

    [8]

    Shamrai K P 1998 Plasma Sources Sci. Technol. 7 499Google Scholar

    [9]

    Chen F F, Blackwell D D 1999 Phys. Rev. Lett. 82 2677Google Scholar

    [10]

    Blackwell D D, Chen F F 2001 Plasma Sources Sci. Technol. 10 226Google Scholar

    [11]

    Kline J, Scime E 2003 Phys. Plasmas 10 135Google Scholar

    [12]

    Kim S, Hwang Y 2008 Plasma Phys. Controlled Fusion 50 035007Google Scholar

    [13]

    Isayama S, Hada T, Shinohara S, et al. 2016 Phys. Plasmas 23 063513Google Scholar

    [14]

    成玉国, 程谋森, 王墨戈等 2014 物理学报 63 035203Google Scholar

    Cheng Y G, Cheng M S, Wang M G, et al. 2014 Acta Phys. Sin. 63 035203Google Scholar

    [15]

    平兰兰, 张新军, 杨桦等 2019 物理学报 68 205201Google Scholar

    Ping L L, Zhang X J, Yang H, et al. 2019 Acta Phys.Sin. 68 205201Google Scholar

    [16]

    Chen F F, Arnush D 1997 Phys. Plasmas 4 3411Google Scholar

    [17]

    Arnush D, Chen F F 1998 Phys. Plasmas 5 1239Google Scholar

    [18]

    Arnush D 2000 Phys. Plasmas 7 3042Google Scholar

    [19]

    Sakawa Y, Kunimatsu H, Kikuchi H, et al. 2003 Phys. Rev. Lett. 90 105001Google Scholar

    [20]

    Vey B 1984 PFC/RR-84-12 (Cambridge: Massachusetts Institute of Technology) pp30–31

    [21]

    Niemi K, Krämer M 2008 Phys. Plasmas 15 073503Google Scholar

    [22]

    Shamrai K P, Shinohara S 2001 Phys. Plasmas 8 4659Google Scholar

    [23]

    Huba J D 2016 NRL Plasma Formulary (Washington: Naval Research Laboratory) p34

    [24]

    Fried B D, Conte S D 2015 The Plasma Dispersion Function: the Hilbert Transform of the Gaussian (New York: Academic Press) pp1–3

    [25]

    Mouzouris Y, Scharer J E 1998 Phys. Plasmas 5 4253Google Scholar

    [26]

    Shoji T, Sakawa Y, Nakazawa S, et al. 1993 Plasma Sources Sci. Technol. 2 5Google Scholar

    [27]

    Swanson D G 1989 Plasma Waves (New York: Academic Press) pp375–376

    [28]

    Lafleur T, Charles C, Boswell R 2010 Phys. Plasmas 17 073508Google Scholar

  • [1] 杨温渊, 董烨, 孙会芳, 杨郁林, 董志伟. 超宽带等离子体相对论微波噪声放大器的物理分析和数值模拟. 物理学报, 2023, 72(5): 058401. doi: 10.7498/aps.72.20222061
    [2] 李文秋, 唐彦娜, 刘雅琳, 马维聪, 王刚. 各向同性等离子体覆盖金属天线辐射增强现象分析. 物理学报, 2023, 0(0): . doi: 10.7498/aps.72.20230101
    [3] 李文秋, 唐彦娜, 刘雅琳, 王刚. 电子温度各向异性对螺旋波等离子体中电磁模式的传播及功率沉积特性的影响. 物理学报, 2023, 72(5): 055202. doi: 10.7498/aps.72.20222048
    [4] 苏瑞霞, 黄霞, 郑志刚. 耦合Frenkel-Kontorova双链的格波解及其色散关系. 物理学报, 2022, 71(15): 154401. doi: 10.7498/aps.71.20212362
    [5] 季佩宇, 黄天源, 陈佳丽, 诸葛兰剑, 吴雪梅. 螺旋波等离子体制备多种碳基薄膜原位诊断研究. 物理学报, 2021, 70(9): 097201. doi: 10.7498/aps.70.20201809
    [6] 杨建荣, 毛杰键, 吴奇成, 刘萍, 黄立. 强碰撞磁化尘埃等离子体中的漂移波. 物理学报, 2020, 69(17): 175201. doi: 10.7498/aps.69.20200468
    [7] 李文秋, 赵斌, 王刚. 电子温度对螺旋波等离子体中电磁模式能量沉积特性的影响. 物理学报, 2020, 69(21): 215201. doi: 10.7498/aps.69.20201018
    [8] 平兰兰, 张新军, 杨桦, 徐国盛, 苌磊, 吴东升, 吕虹, 郑长勇, 彭金花, 金海红, 何超, 甘桂华. 螺旋波等离子体原型实验装置中天线的优化设计与功率沉积. 物理学报, 2019, 68(20): 205201. doi: 10.7498/aps.68.20182107
    [9] 张凯, 杜春光, 高健存. 长程表面等离子体的增强效应. 物理学报, 2017, 66(22): 227302. doi: 10.7498/aps.66.227302
    [10] 李文秋, 王刚, 苏小保. 非磁化冷等离子体柱中的模式辐射特性分析. 物理学报, 2017, 66(5): 055201. doi: 10.7498/aps.66.055201
    [11] 李小泽, 滕雁, 王建国, 宋志敏, 张黎军, 张余川, 叶虎. 过模结构表面波振荡器模式选择. 物理学报, 2013, 62(8): 084103. doi: 10.7498/aps.62.084103
    [12] 刘三秋, 国洪梅. 极端相对论快电子分布等离子体中横振荡色散关系. 物理学报, 2011, 60(5): 055203. doi: 10.7498/aps.60.055203
    [13] 刘炳灿, 逯志欣, 于丽. 金属和Kerr非线性介质界面上表面等离子体激元的色散关系. 物理学报, 2010, 59(2): 1180-1184. doi: 10.7498/aps.59.1180
    [14] 季沛勇, 鲁楠, 祝俊. 量子等离子体中波的色散关系以及朗道阻尼. 物理学报, 2009, 58(11): 7473-7478. doi: 10.7498/aps.58.7473
    [15] 王 亮, 曹金祥, 王 艳, 牛田野, 王 舸, 朱 颖. 电磁脉冲在实验室等离子体中传播时间的实验研究. 物理学报, 2007, 56(3): 1429-1433. doi: 10.7498/aps.56.1429
    [16] 赵国伟, 徐跃民, 陈 诚. 等离子体天线色散关系和辐射场数值计算. 物理学报, 2007, 56(9): 5298-5303. doi: 10.7498/aps.56.5298
    [17] 周国成, 曹晋滨, 王德驹, 蔡春林. 无碰撞等离子体电流片中的低频波. 物理学报, 2004, 53(8): 2644-2653. doi: 10.7498/aps.53.2644
    [18] 谢鸿全, 刘濮鲲, 李承跃, 鄢 扬, 刘盛纲. 等离子体填充波纹波导中低频模式特性分析. 物理学报, 2004, 53(9): 3114-3118. doi: 10.7498/aps.53.3114
    [19] 于 威, 刘丽辉, 侯海虹, 丁学成, 韩 理, 傅广生. 螺旋波等离子体增强化学气相沉积氮化硅薄膜. 物理学报, 2003, 52(3): 687-691. doi: 10.7498/aps.52.687
    [20] 喻 胜, 李宏福, 谢仲怜, 罗 勇. 渐变复合腔回旋管高次谐波注-波互作用非线性模拟. 物理学报, 2000, 49(12): 2455-2459. doi: 10.7498/aps.49.2455
计量
  • 文章访问数:  3567
  • PDF下载量:  57
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-09
  • 修回日期:  2020-03-31
  • 刊出日期:  2020-06-05

螺旋波等离子体中螺旋波与Trivelpiece-Gould波模式耦合及线性能量沉积特性参量分析

  • 1. 中国科学院空天信息创新研究院, 北京 100094
  • 2. 普林斯顿大学, 普林斯顿等离子体物理实验室, 新泽西 08543, 美国
  • 3. 中国科学院大学电子电气与通信工程学院, 北京 100049
  • 4. 北京微电子技术研究所, 北京 100094
  • 通信作者: 李文秋, beiste@163.com
    基金项目: 国家级-国家“万人计划”科技创新领军人才(Y8BF130272)

摘要: 采用有限温度等离子体介电张量模型, 在考虑粒子热效应情形下, 通过求解传导边界条件下等离子体柱中本征模的色散关系, 分析了螺旋波等离子体中典型参量条件下螺旋波与Trivelpiece-Gould (TG) 波的耦合特性及线性能量沉积特性. 在ω/(2π) = 13.56 MHz和TeV,i = 0.1TeV,e参量条件下计算结果表明: 对于螺旋波, 存在截止静磁场B0,H,cutoff与截止等离子体密度n0,H,cutoff, 在B0 > B0,H,cutoffn0 < n0,H,cutoff条件下, 螺旋波变为消逝波; 在ω/ωce ∈ (0.01, 0.10)范围内, 对于m = 0 角向模, TG波Landau阻尼致使的能量沉积占主导地位, 而对于m = 1角向模, 螺旋波Landau阻尼或TG波Landau阻尼致使的能量沉积哪个占据主导地位则取决于B0的大小; 在ωpe/ωce ∈ (3, 100)范围内, TG波Landau阻尼致使的能量沉积占主导地位; 在整体能量沉积过程中, 对于m = 0模和m = 1模, Landau阻尼致使的能量沉积均占据主导地位.

English Abstract

    • 螺旋波等离子体源(helicon plasma sources, HPS)因其高电离率产生的高等离子体密度特性, 在纳米电路刻蚀、深空电磁推进系统、卫星姿态调控、激光等离子体尾波场加速及引力波探测等领域正受到越来越多的关注[1-3]. 自1960年 Aigrain[4]首次使用螺旋波(helicon wave, H-mode)的概念描述高电导率介质中传播的低频电磁波, 1991年Chen[5]首次利用螺旋波的Landau阻尼机制解释HPS中存在的高电离率现象以来, 关于HPS中的高电离率物理机制至今仍作为热点问题得到许多物理学家的研究. 1994年—1998年, Shamrai 等[6-8]发现在螺旋波放电中存在一种静电性质的 Trivelpiece-Gould (TG)波, 其通过在等离子体柱表面附近与螺旋波进行模式耦合而将能量沉积入等离子体中, 但其无法解释在等离子体柱中心处出现的峰值密度现象. 1999年—2001年, Chen 和Blackwell[9,10]通过测量螺旋波等离子体中电子能量分布函数, 发现在螺旋波放电中并不存在强Landau阻尼引起的高能电子分布, 从而推翻了1991年提出的螺旋波强Landau阻尼机制. 2003年, Kline和Scime[11]认为螺旋波(泵浦波)可以通过参数不稳定性模式转换为离子声波和TG波, 而离子声湍流可在局部通过散射加热电子. 2008年, Kim和Hwang[12]指出碰撞效应在体模式转换过程中起到重要作用, 并进一步指出模式耦合层(mode coupling surface, MCS)附近螺旋波的碰撞阻尼率足以解释等离子体柱中心密度峰值现象. 2016年, Shinohara[3]和Isayama等[2,13]指出, 当碰撞效应非常小时, 螺旋波与TG波之间的体模式转换将在波能量沉积过程中起到重要作用. 国内方面, 成玉国等[14]采用线性扰动波假设, 分析了轴向静磁场对螺旋波和TG波能量沉积的影响; 平兰兰等[15]利用Helic code[16-18], 分析了轴向静磁场和等离子体密度对功率沉积及其分布的影响. 至今, 虽然关于螺旋波等离子体中能量沉积的全部物理机理并未完全明朗, 但越来越多的细节正在被不断揭晓, 这为我们最终彻底理解螺旋波等离子体的高电离机制提供了有力支撑.

      由于考虑粒子热效应会将求解本征模色散关系过程变得异常复杂, 故至今有关求解螺旋波和TG波横向波数的研究均采用冷等离子体模型, 为研究精确、完整起见, 实际情形中粒子的热效应应当予以考虑. 故本文利用包含粒子热效应的有限温度等离子体介电张量模型, 在不同等离子体色散函数(plasma dispersion function, PDF)近似条件下, 通过理论求解螺旋波和TG波本征模色散方程, 分析了螺旋波等离子体典型参量情形下螺旋波和TG波的耦合色散特性, 并参数分析了碰撞阻尼、Landau阻尼、多普勒频移回旋阻尼及异常多普勒阻尼四种线性能量沉积机制致使的波能量沉积特性. 本文所得分析结果对揭示螺旋波等离子体高电离率物理机理具有一定的理论意义.

    • 图1所示为典型螺旋波等离子体源中放电管横向界面示意图, 半径为$a$的有限温度等离子体柱浸没在沿$z$轴的均匀静磁场${B_0}{e_z}$中, 且被传导金属边界包裹. 等离子体柱中径向等离子体密度均匀分布, 且只存在电子及携带单位正电荷的氩离子, 其中电子温度远高于粒子温度.

      图  1  被传导边界包裹的等离子体柱横向截面示意图

      Figure 1.  Cross section of plasma column surround by conducting boundary

    • 假定扰动电磁场的传播因子为${{\rm{e}}^{{\rm{j}}\left( {m\theta + {k_{zm}}z - \omega t} \right)}}$, 利用${\rm{Maxwell}}$方程组:

      $\nabla \times {{E}} = - \frac{{\partial {{B}}}}{{\partial t}},$

      $\nabla \times {{B}} = - \frac{{{\rm{j}}\omega }}{{{c^2}}}\epsilon \cdot {{E}},$

      其中${k_{zm}}$为第$m$个角向模的轴向波数($m$为角向模数); $\epsilon$为温等离子体介电张量. 借助等价符号$\partial /\partial z = {\rm{j}}{k_z}$, $\partial /\partial t = - {\rm{j}}\omega $, $\partial /\partial \theta = {\rm{j}}m$, 求得如下关于横向电磁场分量${E_z}$${B_z}$满足的波动方程:

      $\left( {\nabla _ \bot ^2 + k_{ \bot m,H}^2} \right)\left( {\nabla _ \bot ^2 + k_{ \bot m,TG}^2} \right)\left( {\begin{aligned} {{E_{zm}}}\\ {{B_{zm}}} \end{aligned}} \right)=0,\;\;r \leqslant a,$

      其中$\nabla _ \bot ^2 = \dfrac{1}{r}\frac{\partial }{{\partial r}}\left( {r\dfrac{\partial }{{\partial r}}} \right) + \dfrac{1}{{{r^2}}}\dfrac{{{\partial ^2}}}{{\partial {\theta ^2}}}$为横向拉普拉斯算子; ${k_{ \bot m, H}}$${k_{ \bot m, TG}}$分别为螺旋波、TG横向波数, 其精确解可由以下等离子体中电磁波的一般色散关系得到[19]:

      $\begin{split} &{\epsilon_ \bot }k_ \bot ^4 + k_ \bot ^2\left[ {k_z^2\left( {{\epsilon_ \bot } + {\epsilon_\parallel }} \right) - k_0^2\left( {{\epsilon_ \bot }{\epsilon_\parallel } + RL} \right)} \right] \\ &+ {\epsilon_\parallel }\left( {k_z^2 - k_0^2R} \right)\left( {k_z^2 - k_0^2L} \right) = 0,\end{split}$

      其中$R = {\epsilon_ \bot } + {\epsilon_† }$, $L = {\epsilon_ \bot } - {\epsilon_† }$. ${\epsilon_ \bot }$, ${\epsilon_† }$, ${\epsilon_\parallel }$为温度各向同性、无漂移、零阶有限拉莫尔半径效应近似下的温等离子体介电张量元素[20]:

      $ { \epsilon} = \left( {\begin{array}{*{20}{c}} {{\epsilon_ \bot }}&{ - {\rm{j}}{\epsilon_† }}&0\\ {{\rm{j}}{\epsilon_† }}&{{\epsilon_ \bot }}&0\\ 0&0&{{\epsilon_\parallel }} \end{array}} \right), $

      其中张量各元素可表示为:

      ${\epsilon_ \bot } = 1 + \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{2\omega {k_z}{v_{{\rm{th}}\ell }}}}} \left[ {Z\left( {{\xi _{ - 1,\ell }}} \right) + Z\left( {{\xi _{1,\ell }}} \right)} \right],\tag{6a}$

      ${\epsilon_† } = \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{2\omega {k_z}{v_{{\rm{th}}\ell }}}}} \left[ {Z\left( {{\xi _{ - 1,\ell }}} \right) - Z\left( {{\xi _{1,\ell }}} \right)} \right],\tag{6b}$

      ${\epsilon_\parallel } = 1 - \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{k_z^2v_{{\rm{th}}\ell }^2 + {\rm{j}}{\nu _\ell }{k_z}{v_{{\rm{th}}\ell }}Z\left( {{\xi _{0\ell }}} \right)}}} Z'\left( {{\xi _{0\ell }}} \right),\tag{6c}$

      其中${k_z}$为轴向波数; ${\omega _{{\rm{p}}\ell }} = \sqrt {{n_0}{e^2}/({\epsilon_0}{m_\ell })} $$\ell $粒子等离子体频率; ${\omega _{{\rm{c}}\ell }} = e{B_0}/{m_\ell }$$\ell $粒子回旋频率;${\xi _{m\ell }} = (\omega + {\rm{j}}{\nu _\ell } + m{\omega _{c\ell }})/{k_z}{v_{th\ell }}$, ${v_{{\rm{th}}\ell }} = \sqrt {2 e{T_{{\rm{eV}}, \ell }}/{m_\ell }} $$\ell $粒子热速度, ${\nu _\ell }$$\ell $粒子碰撞频率[5,21,22]:

      $ \begin{split} & {\nu _{\rm{e}}} = {\nu _{{\rm{ei}}}} + {\nu _{{\rm{en}}}} \\ =&\; 2.9 \times {10^{ - 12}}{\lambda _{{\rm{ei}}}}\frac{{{n_0}}}{{T_{{\rm{eV,e}}}^{3/2}}}\; + 1.3 \times {10^6}{p_{{\rm{Ar}}}}{T_{{\rm{eV,e}}}}, \end{split} $

      ${\nu _{\rm{i}}} = 4.8 \times {10^{ - 14}}{\lambda _{{\rm{ei}}}}\frac{{{n_0}}}{{\sqrt \varLambda T_{{\rm{eV,i}}}^{3/2}}},$

      其中等离子体密度${n_0}$的单位为${{\rm{m}}^{{\rm{ - 3}}}}$, 氩气气压${p_{{\rm{Ar}}}}$的单位为mTorr (1 mTorr≈0.133 Pa), 电子温度${T_{{\rm{eV, e}}}}$与离子温度${T_{{\rm{eV, i}}}}$的单位为eV, ${\lambda _{{\rm{ei}}}} = 23 - 0.5 {\rm{ln}}({{10}^{ - 6}}{n_0}/ T_{{\rm{eV, e}}}^3)$为电子-离子库仑对数[23], $\varLambda = 40$为氩元素原子量, ${m_\ell }$$\ell $粒子质量, $e = 1.602 \times {10^{ - 19}}\;{\rm{C}}$为电子电量, $\omega $为波频率, ${k_0} = \omega /c$为自由空间波数, $c$为光速; $Z(\xi )$为等离子体色散函数[24]:

      $Z(\xi ) = \frac{1}{{\sqrt {\text{π}} }}\int_{ - \infty }^{ + \infty } {\frac{{{{\rm{e}}^{ - {\tau ^2}}}}}{{\tau - \xi }}} \;{\rm{d}}\tau, \;{\rm{Im}} \left\{ \xi \right\} > 0,$

      $Z'(\xi ) = - 2\left[ {1 + \xi \cdot Z(\xi )} \right].$

      在螺旋波等离子体频率范围内, ${\omega _{{\rm{ci}}}} \ll \omega < {\omega _{{\rm{ce}}}} \ll {\omega _{{\rm{pe}}}}$, (4)式可化简为[25]

      $\frac{{{c^2}{k^2}}}{{{\omega ^2}}} = \frac{{{\omega _{{\rm{pe}}}}^2}}{{\omega \left[ {{\omega _{{\rm{ce}}}}\cos \theta - \omega \left( {1 + {\rm{j}}{\gamma _{\rm{e}}}} \right)} \right]}},$

      其中${\gamma _{\rm{e}}} = {\nu _{\rm{e}}}/\omega $, $k$为总波数, $\cos \theta = {k_z}/k$. 求解(11)式, 得到关于总波数$k$的两个解:

      ${k_{{\rm{H,TG}}}} = \frac{{{k_z}}}{{2{\zeta _{\rm{e}}}}}\left[ {1 \mp \sqrt {1 - 4{{\left( {\frac{{{\zeta _{\rm{e}}}}}{{{k_z}{\delta _0}}}} \right)}^2}} }\; \right],$

      其中${\zeta _{\rm{e}}} = \left( {\omega + {\rm{j}}{\nu _{\rm{e}}}} \right)/{\omega _{{\rm{ce}}}}$, ${\delta _0} = c/{\omega _{{\rm{pe}}}}$为等离子体趋肤深度. 至此, 可根据(12)式、波数间关系得到横向波数${k_{ \bot m, {\rm{H}}}}$${k_{ \bot m, {\rm{TG}}}}$的具体数值.

      求解波动方程(3)式, 得到等离子体区域磁场分量分布[15]:

      $\begin{split}{B_{zm}} =\;& {B_{zm,{\rm{H}}}} + {B_{zm,{\rm{TG}}}} \\ =\;& - 2{\rm{j}}[{f_m}{k_{ \bot m,{\rm{H}}}}{{\rm{J}}_m}\left( {{k_{ \bot m,{\rm{H}}}}r} \right) \\ & + {g_m}{k_{ \bot m,{\rm{TG}}}}{{\rm{J}}_m}\left( {{k_{ \bot m,{\rm{TG}}}}r} \right)],\qquad \quad \end{split}$

      $\begin{split} {B_{\theta m}} =\; & {B_{\theta m,{\rm{H}}}} + {B_{\theta m,{\rm{TG}}}} \\=\; & {\rm{j}}{f_m}[\left( {{k_{\rm{H}}} + {k_{zm}}} \right){{\rm{J}}_{m - 1}}\left( {{k_{ \bot m,{\rm{H}}}}r} \right) \\ &- \left( {{k_{\rm{H}}} - {k_{zm}}} \right){{\rm{J}}_{m + 1}}\left( {{k_{ \bot m,{\rm{H}}}}r} \right)] \\ & + {\rm{j}}{g_m}[\left( {{k_{{\rm{TG}}}} + {k_{zm}}} \right){{\rm{J}}_{m - 1}}\left( {{k_{ \bot m,{\rm{TG}}}}r} \right) \\ &- \left( {{k_{{\rm{TG}}}} - {k_{zm}}} \right){{\rm{J}}_{m + 1}}\left( {{k_{ \bot m,{\rm{TG}}}}r} \right)], \end{split} $

      $\begin{split} {B_{rm}} =\;& {B_{rm,{\rm{H}}}} + {B_{rm,{\rm{TG}}}} \\ =\; & \;{f_m}[\left( {{k_{\rm{H}}} + {k_{zm}}} \right){{\rm{J}}_{m - 1}}\left( {{k_{ \bot m,{\rm{H}}}}r} \right) \\ &+ \left( {{k_{\rm{H}}} - {k_{zm}}} \right){{\rm{J}}_{m + 1}}\left( {{k_{ \bot m,{\rm{H}}}}r} \right)] \\ & + {g_m}[\left( {{k_{{\rm{TG}}}} + {k_{zm}}} \right){{\rm{J}}_{m - 1}}\left( {{k_{ \bot m,{\rm{TG}}}}r} \right) \\ &+\left( {{k_{{\rm{TG}}}} - {k_{zm}}} \right){{\rm{J}}_{m + 1}}\left( {{k_{ \bot m,{\rm{TG}}}}r} \right)]. \end{split} $

      其中, Jmm阶 Bessel函数. 利用${\rm{Maxwell}}$方程, 得到电场分量分布:

      ${E_{zm}} = - {\rm{j}}\omega \delta _{\rm{0}}^{\rm{2}}\left[ {{k_{\rm{H}}}{B_{zm,{\rm{H}}}} + {k_{{\rm{TG}}}}{B_{zm,{\rm{TG}}}}} \right],$

      ${E_{rm}} = \frac{\omega }{{{k_{zm}}}}{B_{\theta m}} - {\rm{j}}\frac{1}{{{k_{zm}}}}\frac{{\partial {E_{zm}}}}{{\partial r}},$

      ${E_{\theta m}} = \frac{m}{{{k_{zm}}r}}{E_{zm}} - \frac{\omega }{{{k_{zm}}}}{B_{rm}},$

      其中${f_m}$${g_m}$为场幅值系数.

      最后, 利用边界条件: ${E_z}$${E_\theta }$在边界$r = a$处为零, 即

      ${E_{zm}}\left( {r = a} \right) = 0,$

      ${E_{\theta m}}\left( {r = a} \right) = 0.$

      将场分布代入(19)式和(20)式, 经过整理, 得到角向模数$m$对应的径向模式色散关系如下所示:

      $\begin{split} &\frac{{{k_{\rm{H}}}{k_{ \bot m,{\rm{H}}}}}}{{{k_{{\rm{TG}}}}{k_{ \bot m,{\rm{TG}}}}}} \cdot \frac{{{{\rm{J}}_m}\left( {{k_{ \bot m,{\rm{H}}}}a} \right)}}{{{{\rm{J}}_m}\left( {{k_{ \bot m,{\rm{TG}}}}a} \right)}} \\ =\; & \big[\left( {{k_{\rm{H}}} + {k_{zm}}} \right){{\rm{J}}_{m - 1}}\left( {{k_{ \bot m,{\rm{H}}}}a} \right) \\ &+ \left( {{k_{\rm{H}}} - {k_{zm}}} \right){{\rm{J}}_{m + 1}}\left( {{k_{ \bot m,{\rm{H}}}}a} \right)\big]\\ &\times\big[\left( {{k_{{\rm{TG}}}} + {k_{zm}}} \right){{\rm{J}}_{m - 1}}\left( {{k_{ \bot m,{\rm{TG}}}}a} \right) \\ &+ \left( {{k_{{\rm{TG}}}} - {k_{zm}}} \right){{\rm{J}}_{m + 1}}\left( {{k_{ \bot m,{\rm{TG}}}}a} \right)\big]^{-1}.\end{split}$

      可以看到, (21)式是边界约束条件下波数${k_{zm}}$, ${k_{ \bot m}}$$k$($k = \sqrt {k_{zm}^2 + k_{ \bot m}^2} $)之间关于贝塞尔函数的隐函数方程, 联立(12)式和(21)式, 可得到精确场分布.

    • 在螺旋波等离子体中, 存在由碰撞阻尼(collision damping, CD)、朗道阻尼(Landau damping, LD)、多普勒频移回旋阻尼(Doppler shifted cyclotron damping, DSCD)及异常多普勒阻尼(anomalous Doppler damping, ADD)致使的螺旋波与TG波线性能量沉积机制. 四种阻尼效应通过波-粒互作用机制将电磁波的能量传递给等离子体, 从而实现电子和离子的加热. 利用复数形式坡印廷定理, 得到螺旋波与TG模第$m$个角向模的功率沉积${P_{{\rm{abs}}}}$

      $\begin{split} & {P_{{\rm{abs}}}}(m,r) \\ =\;& {\mathop{\rm Im}\nolimits} \bigg\{ {1 + \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{\omega {k_z}{v_{{\rm{th}}\ell }}}}Z\left( {{\xi _{ - 1,\ell }}} \right)} } \bigg\}{\epsilon_0}{\left| {{E_{ - m}}(r)} \right|^2} \\ & + {\mathop{\rm Im}\nolimits} \bigg\{ {1 + \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{\omega {k_z}{v_{{\rm{th}}\ell }}}}} Z\left( {{\xi _{1,\ell }}} \right)} \bigg\}{\epsilon_0}{\left| {{E_{ + m}}(r)} \right|^2} \\ & + {\mathop{\rm Im}\nolimits} \bigg\{ {1 - \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{k_z^2v_{th\ell }^2 + {\rm{j}}{\nu _\ell }{k_z}{v_{{\rm{th}}\ell }}Z\left( {{\xi _{0\ell }}} \right)}}} Z'\left( {{\xi _{0\ell }}} \right)} \bigg\}\\ & \times {\epsilon_0}{\left| {{E_{zm}}(r)} \right|^2},\\[-12pt] \end{split}$

      其中, ${E_{\_m}} = ({E_{rm}} - {\rm{j}}{E_{\theta m}})/\sqrt 2 $为右旋极化电场分量, ${E_{ + m}} = ({E_{rm}} + {\rm{j}}{E_{\theta m}})/\sqrt 2 $为左旋极化电场分量. (22)式等号右侧第一项为多普勒频移回旋阻尼, 第二项为异常多普勒阻尼, 第三项为联合碰撞 - Landau阻尼. 为区分碰撞阻尼与Landau阻尼, 我们将纵向介电张量元素${\epsilon_\parallel }$重写为如下形式[26]:

      $ \left\{ {\begin{aligned} &{{\epsilon_{{\rm{\parallel,CD}}}} \to 1 - \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{\omega \left( {\omega + {\rm{j}}{\nu _\ell }} \right)}}\;,\;\;{\text{碰撞阻尼}},} } \\ & {{\epsilon_{{\rm{\parallel,LD}}}} \to 1 \!-\! \sum\limits_{\ell = {\rm{i,e}}} {\frac{{\omega _{{\rm{p}}\ell }^2}}{{k_z^2v_{{\rm{th}}\ell }^2}}} Z'\left( {{\xi _{0\ell }}} \right)\;,\;\;{\text{朗道阻尼}}.} \end{aligned}} \right. $

      本文计算中, 中性工质气体为氩气; 等离子体柱半径$a = 2.5\;{\rm{cm}}$; 射频波频率$\omega = 2{\text{π}} \times 13.56\;{\rm{MHz}}$; 离子温度${T_{{\rm{eV}}, {\rm{i}}}} = {T_{{\rm{eV, e}}}}/10$; 根据宗量$\xi $的大小, 等离子体色散函数采取以下两种渐进形式[27]:

      $\begin{split} Z(\xi ) =\; & {\rm{j}}\sqrt {\text{π}} \sum\limits_{n = 0}^\infty {\frac{{{{\left( {{\rm{j}}\xi } \right)}^n}}}{{\Gamma \left( {\dfrac{n}{2} + 1} \right)}}} \\ & = {\rm{j}}\sqrt {\text{π}} {{\rm{e}}^{ - {\xi ^2}}} - 2\xi \bigg[1 - \frac{2}{3}{\xi ^2} + \frac{4}{{15}}{\xi ^4} \bigg. \\ & \bigg.+ {{\left( { - 2} \right)}^n}\frac{{{\xi ^{2n}}}}{{\left( {2n + 1} \right)\left( {2n - 1} \right) \cdots 3 \cdot 1}} \bigg], \\ & \qquad \left| \xi \right| \ll 1, \end{split}$

      $\begin{split} Z(\xi ) = \; &- \frac{1}{\xi }\sum\limits_{n = 0}^\varPi {\frac{{\left( {2n - 1} \right)!!}}{{{{\left( {2{\xi ^2}} \right)}^n}}}} = {\rm{j}}\sigma \sqrt {\text{π}} {{\rm{e}}^{ - {\xi ^2}}} \\ &- \frac{1}{\xi }\left[ {1 + \frac{1}{{2{\xi ^2}}} + \frac{3}{{4{\xi ^4}}} + \cdot \cdot \cdot \frac{{\left( {2n - 1} \right)!!}}{{{{\left( {2{\xi ^2}} \right)}^n}}}} \right], \\ & \qquad \left| \xi \right| \gg 1, \\ &\sigma = \left\{ {\begin{aligned} & {0\;,\;\;\operatorname{Im} (\xi ) > 0,} \\ & {1\;,\;\;\operatorname{Im} (\xi ) = 0,} \\ & {2\;,\;\;\operatorname{Im} (\xi ) < 0,} \end{aligned}} \right.\\[-25pt] \end{split} $

      其中, Π 为某一子项的值不再单调递减时对应的该子项的项序.

    • 在螺旋波等离子体典型参量条件下, ${n_0} = 1 \times {10^{12}}\;{\rm{c}}{{\rm{m}}^{{\rm{ - 3}}}}$, ${B_0} = 50\;{\rm{G}}$(1 G = 104 T), ${p_{{\rm{Ar}}}} = 3\;{\rm{mTorr}}$, 利用(4)式, 精确计算得到不同${T_{{\rm{eV, e}}}}$值与不同${T_{{\rm{eV, i}}}}$值条件下螺旋波与TG波的耦合色散曲线. 由图2(a)可知, 在${T_{{\rm{eV, e}}}}$约为几电子伏条件下, 等离子体柱中同时存在螺旋波与TG波; 每个${k_z}$值对应两个不同数值的横向波数, ${k_{ \bot, {\rm{H}}}}$${k_{ \bot, {\rm{TG}}}}$(其中较大者${k_{ \bot, {\rm{TG}}}}$为TG波的横向波数, 较小者${k_{ \bot, {\rm{H}}}}$为螺旋波的横向波数模); 螺旋波与TG模在MCS处发生模式耦合: ${k_z} = {k_{z, {\rm{MCS}}}} = 0.38/{\rm{c}}{{\rm{m}}}$, ${\rm{Real}}({k_ \bot }) \!=\! {\rm{Real}}({k_{ \bot, {\rm{MCS}}}}){|_{{\rm{Real}}({k_{ \bot, {\rm{H}}}}) = {\rm{Real}}({k_{ \bot, {\rm{TG}}}})}} $=1.8/cm; 曲线上, MCS上侧为TG-mode branch, 下侧为H-mode branch; 此外, MCS的位置对电子温度的依赖较为明显: 随着${T_{{\rm{eV, e}}}}$的增大, ${k_{z, {\rm{MCS}}}}$逐渐增大, 而${\rm{Real}}\;({k_{ \bot, {\rm{MCS}}}})$变化较为微弱; 电子温度的变化对螺旋波与TG波横向波数的影响亦有不同: 在远离MCS处, 电子温度的变化对螺旋波横向波数的影响可以忽略, 而对TG波横向波数的影响甚为显著; 在MCS附近, 电子温度的变化对螺旋波与TG波两者横向波数的影响均较为显著. 另一方面, 当电子温度给定, ${T_{{\rm{eV, e}}}} = 3\;{\rm{eV}}$, 图2(b)表明离子温度对螺旋波与TG波耦合色散关系的影响极其微小, 完全可以忽略, 这是因为离子质量远大于电子质量, 导致其热速度远远小于电子热速度, 从而使得离子热速度与波相速的相互作用远不如电子热速度与波相速的相互作用那样强烈.

      图  2  粒子温度对螺旋波与TG波耦合色散关系的影响 (a)电子温度的影响; (b)离子温度的影响

      Figure 2.  Influence of particle temperature on dispersion relation between helicon and TG waves: (a) Electron temperature effect; (b) ion temperature effect.

      作为影响螺旋波与TG模耦合关系的核心参量之一, 轴向静磁场${T_{{\rm{eV, e}}}} = 3\;{\rm{eV}}$对螺旋波与TG波横向波数的决定关系如图3所示. 在参量TeV, e = 3 eV, kz = 0.25 cm–1, n0 = 1 × 1011 cm–3, pAr = 0.1 mTorr条件下, 图3计算结果表明: 存在一个临界轴向静磁场${B_{0, {\rm{MCS}}}}$, 螺旋波与TG波在${B_0} = {B_{0, {\rm{MCS}}}}$处发生模式耦合; 对于螺旋波, 存在一个截止轴向静磁场${B_{{\rm{0, H, cutoff}}}}$, 当${B_0} > {B_{0, {\rm{H, cutoff}}}}$时螺旋波成为消逝模; 碰撞效应的存在使得${\rm{Real}}\;({k_{ \bot, {\rm{H}}}}) = {\rm{Real}}\;({k_{ \bot, {\rm{TG}}}})$的位置转移至修正模式耦合层(modified mode coupling surface, mMCS), 在B0 = B0, mMCS处, $\operatorname{Im} {\rm{ag}}\;({k_{ \bot, {\rm{H}}}}) \ne 0$, 暗示B0 = B0, mMCS处并非螺旋波与TG波发生模式耦合的位置; TG模可在${B_0} > {B_{0, {\rm{mMCS}}}}$区域内传播, 而螺旋波只能在${B_{0, {\rm{mMCS}}}} < {B_0} < {B_{0, {\rm{H, cutoff}}}}$区域内传播. 作为影响螺旋波与TG模耦合关系的另一个核心参量, 图4所示为等离子体密度参量${n_0}$对螺旋波与TG模横向波数的决定关系, 在参量条件TeV, e = 3 eV, kz = 0.25 cm–1, B0 = 60 G, ${p_{{\rm{Ar}}}} = 0.1\;{\rm{mTorr}}$下, 计算结果表明: 存在一个临界${n_{0, {\rm{MCS}}}}$, 螺旋波与TG波在${n_0} = {n_{0, {\rm{MCS}}}}$处发生模式耦合; 对于螺旋波, 存在一个截止${n_{{\rm{0, H, cutoff}}}}$, 当${n_0} < {n_{0{\rm{, H, cutoff}}}}$时螺旋波成为消逝模; 碰撞效应的存在使得${\rm{Real}}\;({k_{ \bot, {\rm{H}}}}) = {\rm{Real}}\;({k_{ \bot, {\rm{TG}}}})$的位置转移至mMCS, 在${n_0} = {n_{{\rm{0, mMCS}}}}$处, $\operatorname{Im} {\rm{ag}}\;({k_{ \bot, {\rm{H}}}}) \ne 0$; TG模可在${n_0} < {n_{0, {\rm{mMCS}}}}$区域内传播, 而螺旋波只能在${n_{{\rm{0, H, cutoff}}}} < {n_0} < {n_{{\rm{0, mMCS}}}}$区域内传播.

      图  3  螺旋波与TG波横向波数对轴向静磁场的依赖关系

      Figure 3.  The perpendicular wave number of helicon and TG waves given as functions of axial static magnetic field.

      图  4  螺旋波与TG波横向波数对等离子体密度的依赖关系

      Figure 4.  The perpendicular wave number of helicon and TG waves given as functions of plasma density.

      图5描述了在${T_{{\rm{eV, e}}}} = 3\;{\rm{eV}}$, ${p_{{\rm{Ar}}}} = 3\;{\rm{mTorr}}$参数条件下, 螺旋波轴向波数的实部与虚部随轴向静磁场/等离子体密度的变化关系. 图5(a)暗示, 在${n_0} = 1 \times {10^{12}}\;{\rm{c}}{{\rm{m}}^{{\rm{ - 3}}}}$等离子体密度条件下, 在$\omega /{\omega _{{\rm{ce}}}} \approx 0.2$(${B_0} = 30\;{\rm{G}}$)处螺旋波开始出现回旋阻尼, 且随着轴向静磁场的减小, 回旋阻尼强度显著增大; 当波频率$\omega /(2{\text{π}}) = 1\;{\rm{GHz}}$时, 在$\omega /{\omega _{{\rm{ce}}}} \approx $0.8—0.9范围内开始出现回旋阻尼[28]. 图5(b)暗示, 在${B_0} = 30\;{\rm{G}}$条件下, 螺旋波在${\omega _{{\rm{pe}}}}/{\omega _{{\rm{ce}}}} \approx 100$(${n_0} = 1 \times {10^{12}}\;{\rm{c}}{{\rm{m}}^{{\rm{ - 3}}}}$)处开始出现回旋阻尼, 且随着等离子体密度的增大, 回旋阻尼强度逐渐增大.

      图  5  螺旋波轴向波数随参量变化情况 (a)轴向波数随轴向静磁场变化; (b)轴向静磁场随等离子体密度变化

      Figure 5.  The axial wave number of the right hand polarized wave is given as a function of (a) axial static magnetic field and (b) plasma density.

      图6描述了${T_{{\rm{eV, e}}}} = 3\;{\rm{eV}}$, ${n_0} = 1 \times {10^{12}}\;{\rm{c}}{{\rm{m}}^{{\rm{ - 3}}}}$, ${B_0} = 48.4\;{\rm{G}}$($\omega /{\omega _{{\rm{ce}}}} = 0.1$), ${p_{{\rm{Ar}}}} = 3\;{\rm{mTorr}}$参数条件下, 螺旋波与TG波的径向功率沉积分布. 图6(a)显示, 对于$m = 0$模, 螺旋波与TG波碰撞阻尼与Landau阻尼致使的功率沉积均在中心处取得峰值, 且TG波Landau阻尼致使的功率沉积占据主导地位; 图6(b)显示, 对于$m = 1$模, 回旋阻尼与异常多普勒阻尼致使的功率沉积在中心处取得峰值, 而碰撞阻尼与Landau阻尼致使的功率沉积在偏离中心处取得峰值且占据主导地位.

      图  6  螺旋波与TG波径向功率沉积分布 (a) m = 0 角向对称模; (b) m = 1 角向对称模

      Figure 6.  Radial power deposition profiles of the helicon and TG waves for: (a) m = 0 mode; (b) m = 1 mode.

      图7图8描述了${n_0} = 1 \times {10^{12}}\;{\rm{c}}{{\rm{m}}^{{\rm{ - 3}}}}$, B0 = 48.4 G ($\omega /{\omega _{{\rm{ce}}}} = 0.1$), ${p_{{\rm{Ar}}}} = 0\;{\rm{mTorr}}$参数条件下, 螺旋波与TG波的功率沉积随电子温度/离子温度的变化关系. 在螺旋波等离子体典型电子温度范围内, 对于$m = 0$模, 图7(a)显示在${T_{{\rm{eV, e}}}} < 0.5\;{\rm{eV}}$范围内, 螺旋波和TG波碰撞阻尼致使的功率沉积占据主要地位, 而在${T_{{\rm{eV, e}}}} > 0.5\;{\rm{eV}}$范围内, 螺旋波和TG波Landau阻尼致使的功率沉积在整个功率沉积机制中占据主导地位; 对于$m = 1$模, 图7(b)显示在${T_{{\rm{eV, e}}}} < 0.5\;{\rm{eV}}$范围内, TG波碰撞阻尼致使的功率沉积占据主要地位, 而在${T_{{\rm{eV, e}}}} > 0.5\;{\rm{eV}}$范围内, TG波Landau阻尼致使的功率沉积在整个功率沉积机制中占据主导地位. 另一方面, 在给定电子温度(${T_{{\rm{eV, e}}}} = 1\;{\rm{eV}}$)条件下, 图8描述了螺旋波与TG波的功率沉积随离子温度的变化关系, 由图可知, 对于$m = 0$模与$m = 1$模, 离子温度的变化对螺旋波与TG波各类阻尼致使的功率沉积的影响完全可以忽略不计, 与功率沉积随电子温度变化不同的是, 在${T_{{\rm{eV, i}}}}/{T_{{\rm{eV, e}}}} \in (0.1\;, \;10)$范围内, TG波Landau阻尼致使的功率沉积始终在整个功率沉积机制中占据主导地位. 此外, 一个重要的结论是粒子热效应的引入显著地改变了波功率沉积特性: 与仅包含碰撞效应的冷等离子体模型计算结果不同的是, 热效应的计入导致的朗道阻尼、回旋阻尼及异常多普勒阻尼为我们提供了更加清晰的波能量沉积细节特性, 即, 对于$m = 0$模, 图7(a)揭示了螺旋波和TG波Landau阻尼在功率沉积中的主导地位; 而对于$m = 1$模, 图7(b)则揭示了TG波Landau阻尼在功率沉积中的主导地位; 回旋阻尼与异常多普勒阻尼亦对功率沉积有所贡献, 但在当前参量条件下其占比很小.

      图  7  螺旋波与TG波功率沉积随电子温度的变化关系 (a) m = 0 模; (b) m = 1 模

      Figure 7.  Power deposition profiles of helicon and TG waves are given as functions of electron temperature for: (a) m = 0 mode; (b) m = 1 mode.

      图  8  螺旋波与TG波功率沉积随离子温度的变化关系 (a) m = 0 模; (b) m = 1 模

      Figure 8.  Power deposition profiles of helicon and TG waves are given as functions of ion temperature for: (a) m = 0 mode; (b) m = 1 mode.

      图9描述了${T_{{\rm{eV, e}}}} = 3\;{\rm{eV}}$, ${n_0} = 1 \times {10^{12}}\;{\rm{c}}{{\rm{m}}^{{\rm{ - 3}}}}$, ${p_{{\rm{Ar}}}} = 3\;{\rm{mTorr}}$参数条件下, 螺旋波与TG波的功率沉积在$\omega /{\omega _{{\rm{ce}}}} \in (0.01\;, \;0.1)$范围内的变化情况. 由图9(a)可知, 对于$m = 0$模, TG波Landau阻尼致使的功率沉积在整个功率沉积中占据主导地位, 且随着$\omega /{\omega _{{\rm{ce}}}}$的增大这种主导特性逐渐增强; 对于$m = 1$模, 图9(b)表明, 在$\omega /{\omega _{{\rm{ce}}}} \in (0.01\;, \;0.05)$范围内, 螺旋波的Landau阻尼致使功率沉积占据主导地位, 而在$\omega /{\omega _{{\rm{ce}}}} \in (0.05\;, \;0.1)$范围内, TG波的Landau阻尼致使功率沉积占据主导地位. 这些结论表明: 对于不同角向模数, 轴向静磁场对波能量沉积影响不同; 此外我们应注意到, 相比碰撞阻尼和Landau阻尼, 回旋阻尼与异常多普勒阻尼致使的功率沉积始终很小.

      图  9  螺旋波与TG波功率沉积随轴向静磁场的变化关系 (a) m = 0 模; (b) m = 1 模

      Figure 9.  Power deposition profiles of helicon and TG waves are given as functions of axial static magnetic field for: (a) m = 0 mode; (b) m = 1 mode.

      图10描述了${T_{{\rm{eV, e}}}} = 3\;{\rm{eV}}$, ${B_0} = 100\;{\rm{G}}$, ${p_{{\rm{Ar}}}} = $3 mTorr参数条件下, 螺旋波与TG波的功率沉积在${\omega _{{\rm{pe}}}}/{\omega _{{\rm{ce}}}} \in (3\;, \;100)$范围内的变化情况. 对于螺旋波, 其功率沉积随等离子体密度的增大总体呈现上升趋势; 而对于TG波, 其功率沉积随等离子体密度的增大总体呈现下降趋势. 对于$m = 0$模和$m = 1$模, 图10(a)图10(b)表明TG波的功率沉积在整个功率沉积机制中占据主导地位, 更精确地说, 是TG波Landau阻尼致使的能量沉积占据主导作用; 在两个角向模式中, $m = 1$模在TG波Landau阻尼致使的能量沉积过程占据主导地位.

      图  10  螺旋波与TG波功率沉积随等离子体密度的变化关系 (a) m = 0 模; (b) m = 1 模

      Figure 10.  Power deposition profiles of helicon and TG waves are given as functions of plasma density for: (a) m = 0 mode; (b) m = 1 mode.

    • 基于有限温度等离子体均匀填充圆柱传导边界物理模型, 在考虑粒子热效应条件下, 通过理论分析螺旋波等离子体典型参量条件下螺旋波与TG波的耦合色散特性及线性阻尼致使的能量沉积特性, 得到了一些有用的结论: 1)电子温度对螺旋波与TG波模式耦合层的位置有显著影响; 相较于TG波可在更广参量范围内传播, 螺旋波仅能在一定参量范围内传播; 2)螺旋波与TG波角向对称模($m = 0$模)的碰撞阻尼与Landau阻尼致使的能量沉积均在等离子体柱中心轴处取得峰值, 而角向非对称模($m = 1$角)的碰撞阻尼与Landau阻尼致使的能量沉积均在偏离中心轴处取得峰值; 3)在螺旋波等离子体典型电子温度范围${T_{{\rm{eV, e}}}} \in $(3 eV, 5 eV)内, 对于$m = 0$模, 螺旋波与TG波Landau阻尼致使的功率沉积占据主导地位, 而对于$m = 1$角, TG波Landau阻尼致使的功率沉积占据主导地位; 4)在$\omega /{\omega _{{\rm{ce}}}} \in (0.01\;, \;0.10)$范围内, 对于$m = 0$模, TG波Landau阻尼致使的功率沉积占据主导地位, 而对于$m = 1$模, 螺旋波Landau阻尼或TG波Landau阻尼致使的功率沉积何者占据主导地位则取决于轴向静磁场大小; 5)在${\omega _{{\rm{pe}}}}/{\omega _{{\rm{ce}}}} \in (3\;, \;100)$范围内, 随着等离子体密度增大, 螺旋波功率沉积总体呈现上升趋势, 而TG波功率沉积总体呈现下降趋势; 对于$m = 0$模和$m = 1$模, TG波Landau阻尼致使的功率沉积在整个功率沉积中占据主导地位. 综上分析, 在低磁场中等密度螺旋波等离子体中, TG模Landau阻尼致使的功率沉积占据主导地位. 这些结论为我们揭示螺旋波等离子体高电离率物理机理提供了重要的线索.

参考文献 (28)

目录

    /

    返回文章
    返回