-
Non-inductive current drive plays a crucial role in tokamak, especially for its steady state operations. Recently, the helicon wave (HW) has been regarded as a promising tool for driving off-axis plasma current in reactor-grade machine. The lower-hybrid wave (LHW) is the most effective radio-frequency current drive method, however, it has the drawback, which is limited by the conditions of wave accessibility in the high parameter tokamak, making the wave power usually damped at the plasma edge. HW can spiral towards the plasma centre directly under a high electron density. To obtain a long pulse steady state operation of reactor tokamak, the complementarity of HW and LHW in the aspect of driven current distribution in the high parameter tokamak is considered. The synergy current drive of the HW and the LHW is studied numerically in the steady-state scenario of HL-2M. According to the fast wave dispersion relation of plasma, the HW parameters, including its wave frequency and launched parallel refractive index, are obtained firstly. Results of GENRAY code simulation show that a single pass wave power absorption of the HW can be obtained generally through the electron Landau damping and transit time magnetic pumping effects. On the other hand, the LHW parameters are adopted from the equipped system on the machine. Results of single pass wave absorption are also obtained in the case of LHW. And then, the synergy effects of HW and LHW are studied numerically based on the GENRAY/CQL3D models. The cooperation of these two waves results in a broad plasma current distribution along the radial direction (
$\rho = $ 0.2-0.9) in the machine. Taking the electron distribution functions of these waves into account, it is clear that the electrons are accelerated by the HW in the parallel magnetic field direction, resulting in more electrons entering the region of LHW resonance area. As the consequence, a net plasma current appears. Furthermore, a fine-grained parametric scan is performed by changing the launched parallel refractive index of HW, and the results indicate that positive synergy effects can be generally observed once the related wave current drive profiles are overlapped. Finally, the synergy factor is shown to be proportional to this overlap and reaches its maximum value of 1.18 in HL-2M.-
Keywords:
- tokamaks /
- helicon wave /
- lower hybrid wave /
- synergy current drive
1. 引 言
射频波电流驱动是维持托卡马克装置稳态、高性能运行的重要手段. 离子回旋频率内的快波(FW)、低混杂波(LHW)和电子回旋波(ECW)作用下的电流驱动在EAST, HL-2A等装置上得到了充分的理论验证[1–3]及实验应用[4–6]. 近年来, 这些装置的等离子体放电参数显著提升. 最近, EAST装置成功实现可重复性的403 s高约束模(H模)等离子体放电, 芯部电子温度高达8 keV. 新一代HL-2M装置在Greenwald密度比例${f_{\text{G}}}$约为 0.5的条件下, 可以实现等离子体电流${I_{\text{p}}} = $1.0—1.4 MA的高性能运行. 因此, 在高温、高密(高比压)条件下, 如何有效地加热等离子体(驱动等离子体电流)成为了本领域研究的重点.
低混杂波电流驱动(LHCD)在当前托卡马克射频波电流驱动中效率最高, 被广泛应用于国内外大型托卡马克装置中. 然而, 受限于LHW的可近性条件, 高的等离子体密度(${n_{\text{e}}}$)阻碍了波向芯部等离子体传播; 同时, 随着等离子体温度的提高, 强的电子朗道阻尼导致大量波功率沉积在边缘等离子体区域[7]. 值得一提的是, 在高电子比压(${\beta _{\text{e}}}$)等离子体参数条件下, 螺旋波(HW)可以高效离轴驱动等离子体电流[8]. HW是高谐快波, 主要通过电子朗道阻尼和渡越期磁泵效应作用于电子, 其朗道共振原理与LHCD相同, 不同于LHW的是, HW即使在高密度等离子体中也能够直接传播到芯部等离子体区域; 并且理论和数值的模拟结果表明, 在DⅢ-D装置的放电参数条件下, HW的电流驱动效率是ECW或者中性束电流驱动的2—4倍[8]. 在ITER放电参数中的研究表明, HW电流驱动效率仍然高达60 kA/MW[9]. 2021年度的DⅢ-D实验显示, 将300 kW的HW功率成功注入等离子体后, 电子温度显著提升[10]. 此外, 也针对HL-2M装置中的HW电流驱动展开了深入研究, 结果表明, 在强阻尼条件下, HW能够显著地离轴驱动等离子体电流, 其电流驱动效率高达90 kA/MW[11]. 因此, 基于HL-2M装置现有的LHW系统, 理论上加入HW来研究它们之间是否存在正协同效应具有极其重要的意义.
托卡马克装置中的双波协同电流驱动早在20世纪80年代就受到了人们的关注. 1984年, Fidone等[12]理论分析了LHW和ECW的联合电流驱动, 发现ECW产生的高能电子能够有效加强LHW的电流驱动, 从而总的驱动电流大于这两支波单独注入时产生的驱动电流之和, 即发生正协同效应. 该现象当时在很多托卡马克实验中被证实[13–15]. Harvey等[16]基于三维Fokker-Planck方程研究了ECW分别和FW, LHW的协同电流驱动, 在DⅢ-D典型的放电参数下LHW和ECW普遍存在显著的正协同效应, 但是, 在ITER堆级参数下FW和ECW的协同效应可忽略. Yang等[17]基于一维Fuchs模型和CQL3D程序数值分析了EAST装置中两支频率分别为2.45 GHz和4.6 GHz的LHW间的协同效应, 结果表明, 协同因子依赖于两波共振区域的相对距离; 此外, Yang等[18]在探究不同频率LHW间协同的基础上加入了ECW, 更加深入讨论了三波协同的可能, 结果表明, ECW在垂直方向加速电子, 可以使得更多的电子进入到LHW的共振区间, 从而增大了LHW的驱动电流. 最近, Yin等[19]采用GENRAY/CQL3D程序研究了EAST装置中LHW和高谐快波(HHFW)间的协同效应. 由于选择的HHFW频率(1.0 GHz)与其模拟参数下对应的LHW共振频率[20](${f_{{\text{LH}}}} \approx \sqrt {{f_{{\text{ci}}}}{f_{{\text{ce}}}}} \approx $1.1 GHz)接近, 并且等离子体参数较低(${ \beta _{\text{e}}} \sim$1.3%), 使得这两支波在等离子体中都呈现出多次往返(波能量吸收较弱)的现象, 所以计算获得的协同因子随机性较大. 因此, 本文基于HL-2M装置稳态运行模式, 选择高${\beta _{\text{e}}}$等离子体参数条件下波射线能量以单次吸收为主的HW和LHW为研究对象, 利用GENRAY/CQL3D程序研究该装置中HW与LHW的协同效应. 研究结果将为HL-2M装置相关的工程应用提供理论数值参考.
2. 物理模型及参数设置
考虑存在射频波加热和库仑碰撞时, HW和LHW双波协同作用下的电子分布函数$ {f_{\text{e}}} $的演化满足:
$$ \frac{\partial f_{\rm e}}{\partial t} ({\boldsymbol u},\rho,t) = \langle \langle Q (f) \rangle \rangle +\langle \langle C(f)\rangle \rangle, $$ (1) 其中$\left\langle {\left\langle \cdot \right\rangle } \right\rangle $代表对其中的物理量作反弹平均; $C(f)$是库仑碰撞项; $Q(f)$是描述波与粒子相互作用的准线性项; ${\boldsymbol{u}} = {\boldsymbol{p}}/m$是单位静止质量的电子的动量, $m$为电子的静止质量; $\rho $是径向位置坐标. 双波作用下的$Q(f)$可表示为$\dfrac{\partial }{{\partial {\boldsymbol{u}}}} \cdot \Big({{\boldsymbol{D}}_{{\text{HW}}}} \cdot \dfrac{{\partial f}}{{\partial {\boldsymbol{u}}}}\Big) + \dfrac{\partial }{{\partial {\boldsymbol{u}}}} \cdot \Big({{\boldsymbol{D}}_{{\text{LH}}}} \cdot \dfrac{{\partial f}}{{\partial {\boldsymbol{u}}}}\Big)$, 其中, ${{\boldsymbol{D}}_{{\text{HW}}}}$和${{\boldsymbol{D}}_{{\text{LH}}}}$可由相对论下的Kennel和Engelmann表达式得到[21]
$$ \begin{split} {\boldsymbol{D}} =\,& {D_{/ /}}{\hat {\boldsymbol{u}}_{/ /}}{\hat {\boldsymbol{u}}_{/ /}} + {D_{ \bot / /}}{\hat {\boldsymbol{u}}_ {\bot} }{\hat {\boldsymbol{u}}_{/ /}} + {D_{/ / \bot }}{\hat {\boldsymbol{u}}_{/ /}}{\hat {\boldsymbol{u}}_ {\bot} } \\ & + {D_ {\bot} }{\hat {\boldsymbol{u}}_ {\bot} } . \end{split} $$ (2) 这里, ${\hat {\boldsymbol{u}}_{/ /}}$和${\hat {\boldsymbol{u}}_ \bot }$分别是${{\boldsymbol{u}}_{/ /}}$和${{\boldsymbol{u}}_ \bot }$方向上的单位矢量, ${{{u}}_{/ /}}$和${u_ \bot }$分别为单位静止质量的电子动量在平行和垂直于磁场方向上的分量, 且
$$ \begin{split} {D_{/ /}} =\,& \sum\limits_n {\frac{{{q^2}}}{{2{m^2}}}\pi \gamma \delta \left(\gamma \omega - {k_{/ /}}{u_{/ /}} - n\omega _{\rm c}^ \pm \right)}\\ \,&\times {\left| {\vartheta _{n,{\boldsymbol k}}^ \pm } \right|^2}{\left(\frac{{{k_{/ /}}{u_ \bot }}}{{\omega \gamma }}\right)^2}, \\ {D_{/ / \bot }} =\,& {D_{ \bot / /}} = \sum\limits_n {\frac{{{q^2}}}{{2{m^2}}}\pi \gamma \delta \left(\gamma \omega - {k_{/ /}}{u_{/ /}} - n\omega _{\rm c}^ \pm \right)} \\ \,&\times {\left| {\vartheta _{n,{\boldsymbol k}}^ \pm } \right|^2}\left(\frac{{{k_{/ /}}{u_ \bot }}}{{\omega \gamma }}\right)\left(\frac{{n\omega _{\rm c}^ \pm }}{{\omega \gamma }}\right), \\ {D_ \bot } = \,&\sum\limits_n {\frac{{{q^2}}}{{2{m^2}}}\pi \gamma \delta \left(\gamma \omega - {k_{/ /}}{u_{/ /}} - n\omega _{\rm c}^ \pm \right)} \\ \,&\times {\left| {\vartheta _{n,{\boldsymbol k}}^ \pm } \right|^2}{\left(\frac{{n\omega _{\rm c}^ \pm }}{{\omega \gamma }}\right)^2},\\[-1pt] \end{split} $$ (3) 并有
$$ \vartheta _{n,{\boldsymbol k}}^ \pm = {{\text{J}}_{n + 1}}\frac{{{E_x} - {\mathrm{i}}{E_y}}}{2} + {{\text{J}}_{n - 1}}\frac{{{E_x} - {\mathrm{i}}{E_y}}}{2} + \frac{{{u_{/ /}}}}{{{u_ \bot }}}{{\text{J}}_n}{E_{/ /}}, $$ (4) 式中, $q$为电子电量; $\gamma = \sqrt {1 + {{{u^2}}/{{c^2}}}} $是相对论因子, 其中$c$是真空中光速; $\omega $是入射波频率; $\omega _{\rm c}^ \pm $是局域回旋频率(上标“±”代表不同带电粒子的局域回旋频率, 其中“+”代表离子, “–”代表电子); ${k_{/ /}}$和${k_ \bot }$分别是波矢$\boldsymbol{k} $平行、垂直于磁场的分量; $n$是射频波的谐波数, 对于LHW和HW, $n = $0; ${{\text{J}}_n}$为$n$阶贝塞尔函数; ${E_x}$, ${E_y}$和${E_{/ /}}$是电场的各个分量.
在HL-2M装置稳态运行模式下, 基于其较高的等离子体密度和温度, LHW和HW都将是以波射线能量的单次强吸收为主[8,22]; 并且由于HW和LHW的径向波长(${\lambda _{{\text{ rHW}}}} \sim$1 cm, ${\lambda _{{\text{ rLH}}}}\sim$0.2 cm)远小于局域空间尺度(HL-2M装置小半径$a = $0.65 m), 可以采用几何光学近似来计算波射线的传播轨迹[23]. 此外, 我们前期的研究表明, 该运行模式下GENRAY计算获得的波射线轨迹与全波模型(AORSA) 所计算的电场轮廓吻合较好[12]. 因此, 本文采用GENRAY线性模型来计算波的射线轨迹, 再将GENRAY得到的结果耦合到CQL3D准线性模型来计算分析双波协同电流驱动.
HL-2M装置大半径$R = $1.78 m, 小半径$a = $0.65 m. 等离子体平衡由METIS程序[24]集成模拟给出, 环向磁场${B_{\text{T}}} = $2.0 T, 等离子体电流${{I_{\text{p}}} = }$1.2 MA. 等离子体温度、密度剖面, 有效电荷$ {Z_{{\text{eff}}}} $及安全因子$q(\rho )$剖面如图1所示.
在归一化小半径${\rho = }$0.5处, 容易获得${\beta _{\text{e}}} = {n_{\text{e}}}k{T_{\text{e}}}/\left( {{B^2}/2{\mu _0}} \right) \approx $2.0%, 这是一种典型的高${\beta _{\text{e}}}$运行模式. 通过求解HW强阻尼条件$ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a \geqslant $0.5[1](其中$ \overline {k_{ \bot {\text{I}}}^{(\rm e)}} $为电子垂直波矢虚部值的平均值, 可以反映波射线能量在电子中的沉积), 可以获得该模式下HW强吸收对应的波参数范围. 由图2(a)可知(等高线$ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a = $0.5处), 满足强阻尼条件的最小HW频率${f_{{\text{min}}}} = $0.43 GHz, 且HL-2M装置的环向磁场${B_{\text{T}}} = $2.0 T, 对应的LHW共振频率${f_{{\text{LH}}}} \approx $0.92 GHz[20], 又因HW频率${f_{{\text{HW}}}}$需低于LHW共振频率${f_{{\text{LH}}}}$, 则HW的频率范围为0.43—0.92 GHz. 从前期的研究结果分析, 选取${f_{{\text{HW}}}} = $0.6 GHz比较合适[11]. 图2(b)为${f_{{\text{HW}}}} = $0.6 GHz时, 强阻尼条件下${\beta _{\text{e}}}$及${\xi _{\text{e}}}$的取值要求, 一方面, ${\beta _{\text{e}}}$需大于等于1.8%; 另一方面, 在${\beta _{\text{e}}}{=}$2.0%处, 由${n_{/ /}} = c /(\xi_{\rm e} V_{{T_{\text{e}}}})$(${V_{{T_{\text{e}}}}}$表示电子热速度, ${\xi _{\text{e}}}$是波的平行相速度与电子热速度的比值)可以确定HW的平行折射率${n_{{{/ / \rm HW}}}}$[25], 即3.3$ \leqslant {n_{{{/ / \rm HW}}}} \leqslant $5.6. 考虑到HW电流驱动效率${n_{{{/ /}}}}$随增大而降低[1], 选取3.3$ \leqslant {n_{{{/ / \rm HW}}}} \leqslant $4.1为优化的${n_{{{/ /\rm HW}}}}$范围. LHW参数由HL-2M装置现有的LHW系统给定, 即${f_{{\text{LH}}}} = $3.7 GHz, ${n_{{{/ /\rm LH}}}} \sim$2.6.
3. 双波协同电流驱动
为研究高${\beta _{\text{e}}}$等离子体参数条件下HW和LHW的协同效应, 首先使用GENRAY程序计算两支波在HL-2M装置稳态运行模式下的传播轨迹. 波射线从中平面发射, 平行折射率展宽$\Delta {n_{/ /}} = \pm $0.2, 发射功率$P = $1 MW. 图3为LHW和HW在等离子体中传播的射线轨迹. 在图3中, 两支LHW频率均为3.7 GHz, 当${n_{{{/ /\rm LH}}}} = $2.2时, 较高的等离子体密度导致了LHW在边缘等离子体区域形成了多次反射(多次吸收); 而当${n_{{{/ /\rm LH}}}} = $2.6时, LHW在等离子体中的可近性变好, 值得注意的是, 由于强的电子朗道阻尼效应, 波射线反射的幅度不明显, 可以近似认为是单次吸收. 然而, 受限于密度极限问题, 该运行模式下LHW都不能有效地传播到芯部等离子体区域. 另一方面, 计算得到的HW (${f_{{\text{HW}}}} = $0.6 GHz, ${n_{/ /{\rm {HW}}}} = $3.7)轨迹表明, HW可直接传播到芯部等离子体区域且波射线能量被单次吸收. 图3中波轨迹的分布区域表明, HW和LHW在HL-2M装置等离子体电流驱动中有明显的互补性.
为进一步研究双波协同下的驱动电流分布, 由GENRAY和CQL3D程序共同给出波作用下的驱动电流密度剖面. 在GENRAY的计算中, ${P_{{\text{HW}}}} = {P_{{\text{LH}}}} = $4 MW, ${f_{{\text{LH}}}} = $3.7 GHz, ${f_{{\text{HW}}}} = $0.6 GHz, 波射线从3个点发射, 每个点采用20条波射线模拟发射功率谱, 即$\Delta {n_{/ /}} = \pm $0.2的展宽中每条波射线的功率按照$P = {{{{\sin }^2}(n_{/ /}^{{\text{peak}}} - {n_{/ /}})} /{{{(n_{/ /}^{{\text{peak}}} - {n_{/ /}})}^2}}}$分配, 其中LHW的发射谱中包含20%的负方向功率, 负方向的平行折射率峰值$n_{{{/ /\rm LH}}}^{{\text{peak}}}$= –6.5, 主平行折射率峰值仍为2.6, 而HW的负方向功率谱可忽略不计[8], $n_{{{/ / \rm HW}}}^{{\text{peak}}} = $3.3.
图4为由GENRAY程序模拟的波射线轨迹耦合到CQL3D程序中计算得到的驱动电流密度剖面. 其中, $ {I_{{\text{HW}}}} $, $ {I_{{\text{LH}}}} $和$ {I_{{\text{HW + LH}}}} $分别为HW单独驱动时、LHW单独驱动时以及双波协同驱动时的等离子体电流. 由图4可知, HW和LHW的电流驱动在径向分布上具有较好的互补性: HW驱动的等离子体电流主要分布在$\rho = $0.2—0.6, 而LHW驱动的等离子体电流主要分布在$\rho = $0.5—0.9, 这使得双波协同驱动的等离子体电流在较大的径向范围($\rho = $0.2—0.9)内均有分布. 在部分靠近HW驱动电流密度峰值的LHW电流驱动范围($\rho = $0.5—0.7), 双波明显呈正协同效应, 加强了等离子体中间区域的电流驱动; 且双波协同减小了LHW在边缘等离子体($\rho = $0.8—0.9)的能量沉积; 而在HW主导的电流驱动区域($\rho = $0.2—0.45), 双波协同驱动的电流密度剖面与HW单独驱动的电流密度剖面基本重合. 双波协同还使得波驱动等离子体电流的效率明显提高, 双波协同驱动的等离子体电流$ {I_{{\text{HW + LH}}}} $=1228.0 kA明显大于这两支波单独注入时产生的驱动电流之和$ {I_{{\text{HW}}}} $+$ {I_{{\text{LH}}}} $= 440.6 kA + 669.0 kA = 1109.6 kA.
双波协同驱动等离子体电流时, HW和LHW对总电流的贡献不同. 由${n_{/ /}} = {{c{k_{/ /}}}/\omega }$和(3)式中HW/LHW相对论性共振条件$\gamma \omega = {k_{/ /}}{u_{/ /}}$可知: 较高的${n_{/ /}}$对应较低的共振速度${u_{/ /}}$, 对于${n_{{{/ / \rm HW}}}} = $3.3, ${n_{{{/ /\rm LH}}}} = $2.6的情况, HW对应的${u_{/ /}}$较低, 即HW在低${u_{/ /}}$区与电子发生共振相互作用, 同时将这些电子推入相邻的LHW共振区域, 从而增强LHW的电流驱动. 这种协同效应还与共振区间内平行分布函数的平台高度密切相关[26].
图5为正协同效应最强处($\rho = $0.58)准线性扩散和碰撞弛豫共同作用下的电子平行分布, ${D_{{\text{Maxwell}}}}$代表无射频波注入情况下的麦克斯韦分布, ${D_{{\text{HW}}}}$, ${D_{{\text{LHW}}}}$和${D_{{\text{HW + LHW}}}}$分别代表注入HW, LHW, HW+LHW后的电子平行分布. 未进入波共振区的低速电子区域, 碰撞弛豫占主导地位, 电子平行分布均呈现麦克斯韦分布. 波共振区域的准线性扩散作用最为明显, 倾向于拉平分布, 建立平台. 由图5可知双波协同作用下的电子平行分布在$u/{{u_{{\text{norm}}}}} \lt 0.2$(${u_{{\text{norm}}}}$为电子动量的归一化因子)的范围内与HW单独作用下的电子平行分布基本重合, 而当$0.2 \lt {u /{{u_{{\text{norm}}}}}} \lt 0.23$时双波协同作用下的电子平行分布明显高于LHW单独注入下的电子平行分布, 且低于HW单独注入下的电子平行分布, 这证明HW在平行方向上加速电子, 使这些电子有更大概率进入相邻的LHW共振区间, 从而使得电子平行分布在HW和LHW共振区域的交叠区域($0.2 \lt {u /{{u_{{\text{norm}}}}}} \lt 0.23$)重新分配, 使电子的平行速度更加集中在$0.23 \lt {u /{{u_{{\text{norm}}}}}} \lt 0.30$的范围内. 由LHW的驱动电流正比于平行分布函数在LHW共振区内的平台高度[26]可知, 双波协同促进了LHW电流驱动效率从而有效地增大了双波的总驱动电流, 因此定义协同因子${F_{{\text{syn}}}} = ({I_{{\text{LH + HW}}}} - {I_{{\text{HW}}}}) /{{I_{{\text{LH}}}}}$作为后续描述协同效果的特征参量.
图6展示了在$ \rho = $0.58的通量面处, 波电场加速下的电子通量及准线性扩散下的电子分布. 此通量面由归一化速度空间(${{{u_{/ /}}} / {{u_{{\text{norm}}}}}}$, ${{{u_ \bot }} \mathord{\left/ {\vphantom {{{u_ \bot }} {{u_{{\text{norm}}}}}}} \right. } {{u_{{\text{norm}}}}}}$)表示, 并划分了捕获/通行粒子区域(图中的“V”字形折线为捕获/通行粒子边界). 图6(a)—(c)分别展示了HW (${P_{{\text{HW}}}} = $4 MW, ${f_{{\text{HW}}}} = $600 Hz, ${n_{{{/ / \rm HW}}}} = $3.3)单独作用下、LHW (${P_{{\text{LH}}}} = $4 MW, ${f_{{\text{LH}}}} = $3.7 GHz, ${n_{{{/ /\rm LH}}}} = $2.6)单独作用下和双波协同作用下的电子通量. 其中矢量长度正比于电子通量的大小, 方向代表通量输运方向, 矢量正负号的选取与电子通量大小无关, 而与电子运动方向是否与磁场方向同向相对应. 由图6(a)可知HW直接与速度处于–0.22${{ \lt {u_{/ /}}} /{{u_{{\text{norm}}}}}} \lt $–0.15范围内的电子相互作用, 并将净的平行动量传递给电子; 而LHW可直接加速–0.30${{ \lt {u_{/ /}}}/ {{u_{{\text{norm}}}}}} \lt $–0.20范围内的电子, 且其相互作用更为强烈, 如图6(b)所示. 对比图6(a)和图6(b)可知, HW和LHW的电场作用范围在–0.23${{ \lt {u_{/ /}}}/{{u_{{\text{norm}}}}}} \lt $–0.20内交叠, 且LHW共振区处于更高的电子相速度区间, 这反映HW在低平行速度区与共振电子相互作用, 可将共振电子推入相邻的LHW共振区, 产生正协同效应. 由图6(c)可知, 协同后射频电场作用的区域($-0.30 \lt { {u_{/ /}}} / {{u_{{\text{norm}}}}} \lt -0.15$)更为广泛, 使波在更宽的径向范围内驱动等离子体电流, 且在这3种情况下射频波电场都并不能直接加速捕获电子.
图6(d)—(f)分别展示了HW和碰撞作用下、LHW和碰撞作用下及双波协同和碰撞作用下电子总通量的对数, 其中矢量长度正比于电子总通量的对数的大小. 图6(d)—(f)清晰地展现了由射频场和碰撞弛豫共同作用下引起的旋涡, 因为在稳态下, 电子总通量必须是无散度的, 这种电子通量的演化以旋涡的形式展现, 产生旋涡的区域与产生电子通量的区域相对应. 可以看出当等离子碰撞弛豫和射频电场共同作用时, 捕获电子可间接受射频电场影响, 进而影响通行/捕获电子间的转换.
图6(g)—(i)分别为HW和碰撞作用下、LHW和碰撞作用下及双波协同和碰撞作用下的电子分布. 注入HW/LHW后, 电子分布相对于麦克斯韦分布发生偏离, 形成非对称分布, 这会造成等离子体的“非对称电阻性”, 从而驱动电流. 这种分布上的畸变越大, 驱动的等离子体电流会越大. 图6(g)中的电子分布相对于麦克斯韦分布的偏离较小, 这是由于HW的电子朗道阻尼效应较弱, 使得此偏离并不明显. 图6(h)中的电子分布偏离麦克斯韦分布非常明显, 与LHW的强电子朗道阻尼作用相对应. 对比图6(h)和图6(i)可知, 协同后处于–0.30$\lt {{ {u_{/ /}}} \mathord{\left/ {\vphantom {{ \lt {u_{/ /}}} {{u_{{\text{norm}}}}}}} \right. } {{u_{{\text{norm}}}}}} \lt $–0.20范围的电子分布相较于单独注入LHW时的电子分布偏离得更为明显, 表明在此区间内正协同效应显著, 并增强了总的驱动电流.
在强阻尼条件下, 当${f_{{\text{HW}}}} = $0.6 GHz时, ${n_{{{/ / \rm HW}}}} = $3.3—4.1. 为探究不同${n_{{{/ / \rm HW}}}}$值下的协同效果, 表1列出了固定LHW参数(${P_{{\text{LH}}}} = $4 MW, ${f_{{\text{LH}}}} = $3.7 GHz, ${n_{{{/ /\rm LH}}}} = $2.6)时, 不同${n_{{{/ / \rm HW}}}}$值对协同效果的影响(其余HW参数为${P_{{\text{HW}}}} = $4 MW, ${f_{{\text{HW}}}} = $0.6 GHz). 由表1可知, 当${n_{{{/ / \rm HW}}}} = $3.3时, 协同效果最好, 协同因子${F_{{\text{syn}}}} = $1.18, 这是由于此时的HW驱动电流效率最高, 可以加速更多电子使其进入LHW共振区域, 从而有效提高双波联合驱动电流的效率. 此外, 随着${n_{{{/ / \rm HW}}}}$的增大, ${F_{{\text{syn}}}}$的值逐渐减小, 主要因为${n_{{{/ / \rm HW}}}}$的增大导致HW的波阻尼增强, 波与等离子体相互作用的时间变短, 所以波射线在传播的过程中很快被吸收(传播距离近), 从而使得HW的驱动电流减小, 这样会减少加速至LHW共振速度的电子数量, 导致${F_{{\text{syn}}}}$减小.
${n_{{{/ / \rm HW}}}}$ ${I_{{\text{HW}} + {\text{LH}}}}/{\text{kA}}$ ${I_{{\text{HW}}}}/{\text{kA}}$ ${I_{{\text{LH}}}}/{\text{kA}}$ ${F_{{\text{syn}}}}$ 3.3 1228.0 440.6 669.0 1.18 3.5 1211.1 433.3 669.0 1.16 3.7 1180.6 416.0 669.0 1.14 3.9 1151.6 396.5 669.0 1.13 4.1 1120.5 377.0 669.0 1.11 4. 讨论与结论
基于GENRAY/CQL3 D程序, 依托HL-2M装置研究了高$ {\beta _{\text{e}}} $参数条件下以单次吸收为主的HW和LHW间的协同效应. 在典型的HL-2M稳态放电模式中, ${\beta _{\text{e}}}\sim $2.0%, 据此本文结合工程应用, 理论给出了HW强吸收的波参数范围, 即${f_{{\text{HW}}}} = $0.43—0.90 GHz, ${n_{{{/ / \rm HW}}}} = $3.3—4.1. 通过研究准线性扩散和等离子体碰撞弛豫共同作用下的电子平行分布及相空间内的电子通量分布发现, 当两支波的共振区域有交叠时, HW通过在平行方向上加速电子, 为LHW提供更多共振电子产生双波正协同效应; 另一方面, 双波协同驱动的电流密度剖面具有较好的互补性, 可以在较大的径向范围($ \rho = $0.2—0.9)内有效驱动等离子体电流. 对于HW和LHW, 在注入功率均为4 MW条件下, 双波协同总驱动电流甚至高达1.2 MA, 这将非常有利于HL-2M装置的高性能运行. 另外在强阻尼条件下, 随着HW平行折射率的增大, 协同因子的值会逐渐减小, 主要因为平行折射率的增大导致HW的波阻尼增强, 波与等离子体相互作用的时间变短, 从而减小HW的驱动电流, 这样会使加速至LHW共振速度的电子数量变少, 导致协同因子减小. 在高比压等离子体参数条件下, HW和LHW功率沉积的径向位置分布重叠较少, 因此计算得到的协同因子较小. 但是, 前期的数值模拟表明, HW的电流驱动效率是ECW的2—4倍[8], 这有利于提高等离子体电流驱动效率, 实现托卡马克长脉冲稳态运行. 本文的研究为堆级托卡马克装置中的多波协同电流驱动提供积极理论参考.
[1] Chiu S C, Chan V S, Harvey R W, Porkolab M 1989 Nucl. Fusion 29 2175Google Scholar
[2] Peysson Y, Decker J, Nilsson E, Artaud J F, Ekedahl A, Goniche M, Hillairet J, Ding B, Li M, Bonoli P T, Shiraiwa S, Madi M 2016 Plasma Phys. Controlled Fusion 58 044008Google Scholar
[3] Wang Z T, Long Y X, Dong J Q, He Z X 2013 Chin. Phys. B 22 095201Google Scholar
[4] Phillips C K, Bell R E, Berry L A, Bonoli P T, Harvey R W, Hosea J C, Jaeger E F, LeBlanc B P, Ryan P M, Taylor G, Valeo E J, Wilgen J B, Wilson J R, Wright J C, Yuh H, NSTX Team 2009 Nucl. Fusion 49 075015Google Scholar
[5] Huang C B, Gao X, Liu Z X, Han X, Zhang T, Wang Y M, Zang S B, Kong D F, EAST Team 2016 Plasma Phys. Controlled Fusion 58 075005Google Scholar
[6] Duan X R, Ding X T, Dong J Q, Yan L W, Liu Y, Huang Y, Song X M, Zou X L, Xu M, Yang Q W, Liu D Q, Rao J, Xuan W M, Chen L Y, Mao W C, Wang Q M, Cao J Y, Lei G J, Zhang J H, Li X D, Chen W, Zhao K J, Xiao W W 2013 Nucl. Fusion 53 104009Google Scholar
[7] Li X X, Xiang N, Gan C Y 2015 Chin. Phys. Lett. 32 035202Google Scholar
[8] Prater R, Moeller C P, Pinsker R I, Porkolab M, Meneghini O, Vdovin V L 2014 Nucl. Fusion 54 083024Google Scholar
[9] Vdovin V L 2013 Plasma Phys. Rep. 39 95Google Scholar
[10] Van Compernolle B, Brookman M, Pinsker R, Moeller C, Squire J, Garofalo A M, Nagy A, Torrezan A, Ponce D, Pawley C, Chowdury S, Crocker N, Degrandchamp G, Hinson E, Lohr J, Marinoni A, Martin E, Petty C, Porkolab M, Rost C, Schmitz O, Thome K, Wang Q H, Watkins J, Zeller K 2021 63rd Annual Meeting of the APS Division of Plasma Physics Pittsburgh, November 8–12, 2021 UO07
[11] Liu H B, Liu G N, Sun A P, Xiao Z Y, Li X X 2022 J. Korean Phys. Soc. 81 397Google Scholar
[12] Fidone I, Giruzzi G, Granata G, Meyer R L 1984 Phys. Fluids. 27 2468Google Scholar
[13] Kawashima H, Yamamoto T, Hoshino K, Uesugi Y, Mori M, Suzuki N 1991 Nucl. Fusion 31 495Google Scholar
[14] Maekawa T, Maehara T, Minami T, Kishigami Y, Kishino T, Makino K, Hanada K, Nakamura M, Terumichi Y, Tanaka S 1993 Phys. Rev. Lett. 70 2561Google Scholar
[15] Maehara T, Yoshimura S, Minami T, Hanada K, Nakamura M, Maekawa T, Terumichi Y 1998 Nucl. Fusion 38 39Google Scholar
[16] Harvey R W, Chiu S C, McCoy M G, Kerbel G D, Smith G R, Mau T K 1991 Proc. of IAEA TCM on Fast Wave Current Drive in Reactor Scale Tokamaks Arles, France, September 23–25, 1991 p135
[17] Yang Y L, Xiang N, Hu Y M 2017 Phys. Plasmas 24 032502Google Scholar
[18] Yang Y L, Xiang N, Hu Y M 2017 Phys. Plasmas 24 082503Google Scholar
[19] Yin L, Zheng P W, Gong X Y, Yang C, Yin X H, Song C Y, Huang Q H, Chen Y, Zhong Y J 2022 Nucl. Fusion 62 066023Google Scholar
[20] Pinsker R I, Porkolab M, Petty C C, Prater R, Moeller C P 2015 AIP Conference Proceedings 1689 080012Google Scholar
[21] Kennel C F, Engelmann F 1966 Phys. Fluids 9 2377Google Scholar
[22] 刘祖光 2020 硕士学位论文 (衡阳: 南华大学)
Liu Z G 2020 M. S. Thesis (Hengyang: University of South China
[23] The GENRAY Ray Tracing Code, Smirnov A P, Harvey R W https://compxco.com/Genray_manual.pdf [2003-03-17
[24] Artaud J F, Imbeaux F, Garcia J, Giruzzi G, Aniel T, Basiuk V, Bécoulet A, Bourdelle C, Buravand Y, Decker J, Dumont R, Eriksson L G, Garbet X, Guirlet R, Hoang G T, Huynh P, Joffrin E, Litaudon X, Maget P, Moreau D, Nouailletas R, Pégourié B, Peysson Y, Schneider M, Urban J 2018 Nucl. Fusion 58 105001Google Scholar
[25] 李新霞, 李国壮, 刘洪波 2020 物理学报 69 145201Google Scholar
Li X X, Li G Z, Liu H B 2020 Acta Phys. Sin. 69 145201Google Scholar
[26] 杨友磊, 胡业民, 项农 2017 物理学报 66 245202Google Scholar
Yang Y L, Hu Y M, Xiang N 2017 Acta Phys. Sin. 66 245202Google Scholar
-
图 2 在HL-2M装置放电条件下通过HW强阻尼条件求得的$ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a $值的等高线图 (a) ${\beta _{\text{e}}}\sim$2.0%时, $ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a $关于${\xi _{\text{e}}}$和$f$的等高线图; (b) $ {f_{{\text{HW}}}} = $0.6 GHz时, $ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a $关于${\xi _{\text{e}}}$和${\beta _{\text{e}}}$的等高线图
Figure 2. Contours of $ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a $ as a function of (a) ${\xi _{\text{e}}}$ and $f$ with ${\beta _{\text{e}}}\sim$2.0% for the strong damping condition of the HW of HL-2M; contours of $ 2\overline {k_{ \bot {\text{I}}}^{(\rm e)}} a $ as a function of (b) ${\xi _{\text{e}}}$ and ${\beta _{\text{e}}}$ with ${f_{{\text{HW}}}} = $0.6 GHz for the strong damping condition of the HW of HL-2M.
图 3 HL-2M装置LHW/HW波射线传播轨迹, 其中, ${f_{{\text{LH}}}} = $3.7 GHz, ${n_{{{/ /\rm LH}}}}$分别取2.2和2.6; ${f_{{\text{HW}}}} = $0.6 GHz, ${n_{{{/ / \rm HW}}}} = $3.7
Figure 3. Ray trajectories of the HW with ${f_{{\text{HW}}}} = $0.6 GHz and ${n_{{{/ / \rm HW}}}} = $3.7 in HL-2M, as well as the LHW with ${f_{{\text{LH}}}} = $3.7 GHz and ${n_{{{/ /\rm LH}}}}$ of 2.2 and 2.6 respectively.
图 5 (a)麦克斯韦(${D_{{\text{Maxwell}}}}$), HW(${D_{{\text{HW}}}}$), LHW(${D_{{\text{LH}}}}$)单独作用下和双波协同作用(${D_{{\text{HW}} + {\text{LH}}}}$)下的电子平行分布; (b) 图5(a)中黑色矩形框的放大区域
Figure 5. (a) Parallel electron distributions of the Maxwell, the HW, the LHW, and HW+LHW; (b) the enlarged area of the black rectangular box in Fig. 5(a).
图 6 波电场加速下的电子通量、波电场及碰撞作用下电子总通量的对数及射频波准线性扩散和碰撞作用下的电子分布 (a)—(c)分别为HW, LHW以及双波协同下的电子通量; (d)—(f) 分别为HW, LHW以及双波协同下的电子总通量的对数; (g)—(i)分别为HW, LHW以及双波协同下的电子分布
Figure 6. (a)–(c) Electron flux of the HW, the LHW, and the HW+LHW respectively; (d)–(f) the logarithms of the total electron fluxes for the case of Fig. 6(a)-(c); (g)–(i) contours of the electron distribution function for the case of Fig. 6(a)-(c).
表 1 不同HW平行折射率下的协同效果
Table 1. Synergistic effect in different HW parallel refractive indexes.
${n_{{{/ / \rm HW}}}}$ ${I_{{\text{HW}} + {\text{LH}}}}/{\text{kA}}$ ${I_{{\text{HW}}}}/{\text{kA}}$ ${I_{{\text{LH}}}}/{\text{kA}}$ ${F_{{\text{syn}}}}$ 3.3 1228.0 440.6 669.0 1.18 3.5 1211.1 433.3 669.0 1.16 3.7 1180.6 416.0 669.0 1.14 3.9 1151.6 396.5 669.0 1.13 4.1 1120.5 377.0 669.0 1.11 -
[1] Chiu S C, Chan V S, Harvey R W, Porkolab M 1989 Nucl. Fusion 29 2175Google Scholar
[2] Peysson Y, Decker J, Nilsson E, Artaud J F, Ekedahl A, Goniche M, Hillairet J, Ding B, Li M, Bonoli P T, Shiraiwa S, Madi M 2016 Plasma Phys. Controlled Fusion 58 044008Google Scholar
[3] Wang Z T, Long Y X, Dong J Q, He Z X 2013 Chin. Phys. B 22 095201Google Scholar
[4] Phillips C K, Bell R E, Berry L A, Bonoli P T, Harvey R W, Hosea J C, Jaeger E F, LeBlanc B P, Ryan P M, Taylor G, Valeo E J, Wilgen J B, Wilson J R, Wright J C, Yuh H, NSTX Team 2009 Nucl. Fusion 49 075015Google Scholar
[5] Huang C B, Gao X, Liu Z X, Han X, Zhang T, Wang Y M, Zang S B, Kong D F, EAST Team 2016 Plasma Phys. Controlled Fusion 58 075005Google Scholar
[6] Duan X R, Ding X T, Dong J Q, Yan L W, Liu Y, Huang Y, Song X M, Zou X L, Xu M, Yang Q W, Liu D Q, Rao J, Xuan W M, Chen L Y, Mao W C, Wang Q M, Cao J Y, Lei G J, Zhang J H, Li X D, Chen W, Zhao K J, Xiao W W 2013 Nucl. Fusion 53 104009Google Scholar
[7] Li X X, Xiang N, Gan C Y 2015 Chin. Phys. Lett. 32 035202Google Scholar
[8] Prater R, Moeller C P, Pinsker R I, Porkolab M, Meneghini O, Vdovin V L 2014 Nucl. Fusion 54 083024Google Scholar
[9] Vdovin V L 2013 Plasma Phys. Rep. 39 95Google Scholar
[10] Van Compernolle B, Brookman M, Pinsker R, Moeller C, Squire J, Garofalo A M, Nagy A, Torrezan A, Ponce D, Pawley C, Chowdury S, Crocker N, Degrandchamp G, Hinson E, Lohr J, Marinoni A, Martin E, Petty C, Porkolab M, Rost C, Schmitz O, Thome K, Wang Q H, Watkins J, Zeller K 2021 63rd Annual Meeting of the APS Division of Plasma Physics Pittsburgh, November 8–12, 2021 UO07
[11] Liu H B, Liu G N, Sun A P, Xiao Z Y, Li X X 2022 J. Korean Phys. Soc. 81 397Google Scholar
[12] Fidone I, Giruzzi G, Granata G, Meyer R L 1984 Phys. Fluids. 27 2468Google Scholar
[13] Kawashima H, Yamamoto T, Hoshino K, Uesugi Y, Mori M, Suzuki N 1991 Nucl. Fusion 31 495Google Scholar
[14] Maekawa T, Maehara T, Minami T, Kishigami Y, Kishino T, Makino K, Hanada K, Nakamura M, Terumichi Y, Tanaka S 1993 Phys. Rev. Lett. 70 2561Google Scholar
[15] Maehara T, Yoshimura S, Minami T, Hanada K, Nakamura M, Maekawa T, Terumichi Y 1998 Nucl. Fusion 38 39Google Scholar
[16] Harvey R W, Chiu S C, McCoy M G, Kerbel G D, Smith G R, Mau T K 1991 Proc. of IAEA TCM on Fast Wave Current Drive in Reactor Scale Tokamaks Arles, France, September 23–25, 1991 p135
[17] Yang Y L, Xiang N, Hu Y M 2017 Phys. Plasmas 24 032502Google Scholar
[18] Yang Y L, Xiang N, Hu Y M 2017 Phys. Plasmas 24 082503Google Scholar
[19] Yin L, Zheng P W, Gong X Y, Yang C, Yin X H, Song C Y, Huang Q H, Chen Y, Zhong Y J 2022 Nucl. Fusion 62 066023Google Scholar
[20] Pinsker R I, Porkolab M, Petty C C, Prater R, Moeller C P 2015 AIP Conference Proceedings 1689 080012Google Scholar
[21] Kennel C F, Engelmann F 1966 Phys. Fluids 9 2377Google Scholar
[22] 刘祖光 2020 硕士学位论文 (衡阳: 南华大学)
Liu Z G 2020 M. S. Thesis (Hengyang: University of South China
[23] The GENRAY Ray Tracing Code, Smirnov A P, Harvey R W https://compxco.com/Genray_manual.pdf [2003-03-17
[24] Artaud J F, Imbeaux F, Garcia J, Giruzzi G, Aniel T, Basiuk V, Bécoulet A, Bourdelle C, Buravand Y, Decker J, Dumont R, Eriksson L G, Garbet X, Guirlet R, Hoang G T, Huynh P, Joffrin E, Litaudon X, Maget P, Moreau D, Nouailletas R, Pégourié B, Peysson Y, Schneider M, Urban J 2018 Nucl. Fusion 58 105001Google Scholar
[25] 李新霞, 李国壮, 刘洪波 2020 物理学报 69 145201Google Scholar
Li X X, Li G Z, Liu H B 2020 Acta Phys. Sin. 69 145201Google Scholar
[26] 杨友磊, 胡业民, 项农 2017 物理学报 66 245202Google Scholar
Yang Y L, Hu Y M, Xiang N 2017 Acta Phys. Sin. 66 245202Google Scholar
期刊类型引用(1)
1. 樊浩,陈少永,牟茂淋,刘泰齐,张业民,唐昌建. 低杂波注入对剥离气球模的作用. 物理学报. 2024(09): 210-221 . 百度学术 其他类型引用(0)
Catalog
Metrics
- Abstract views: 2417
- PDF Downloads: 78
- Cited By: 1