[论文解读] Theoretical analysis of the saturation phase of the $1/1$ energetic-ion-driven resistive interchange mode

J. Varela, D. A. Spong+5 more|arXiv (Cornell University)|Jan 3, 2022

Magnetic confinement fusion research参考文献 40被引用 15

一句话总结

本研究采用FAR3d代码进行非线性磁流体模拟，模拟了大型螺旋装置（LHD）中1/1高能离子驱动的电阻性磁岛模（EIC）的饱和阶段。模拟再现了关键实验特征，包括随高能粒子β增加而增长的模振幅与宽度、阈值激发、9.4 kHz的突发事件及其泛音、由于二次模耦合引起的向内传播，以及与压强剖面恢复相关的稳定化现象。

ABSTRACT

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. A set of non linear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The non linear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP $\beta$ increases, the EP beta threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the non linear destabilization of the 3/4 and 2/3 EPMs, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

研究动机与目标

• 模拟大型螺旋装置（LHD）中1/1高能离子驱动的电阻性磁岛模（EIC）的非线性饱和阶段。
• 再现实验观测到的EIC多阶段演化过程，包括突发事件和模传播。
• 研究螺旋阱高能粒子（EPs）在EIC失稳与饱和中的作用。
• 分析模诱导的压强平坦化与随后RIC再激活之间的反馈机制。
• 验证模拟模型在实验频率、模结构和β阈值行为方面的准确性。

提出的方法

• 采用FAR3d磁流体代码求解与高能粒子（EPs）矩方程耦合的简化电阻性MHD方程。
• 对热 plasma 使用简化MHD模型，对EPs使用磁流体模型，并采用与解析TAE增长率匹配的闭合关系。
• 修改平均漂移速度算子，以包含螺旋阱EPs的旋转变迁漂移。
• 引入线性波-粒子共振效应（Landau阻尼/增益）和平行动量响应，以实现与几何声波的耦合。
• 使用VMEC代码生成的LHD放电116190固定边界平衡态，参数与实验剖面匹配。
• 在径向方向采用有限差分法，沿环向和极向角采用傅里叶展开，并采用半隐式、显式两步法以实现(Δt)³阶精度。

实验结果

研究问题

• RQ1 在存在螺旋阱高能粒子的情况下，1/1 EIC如何演化其饱和阶段？
• RQ2 LHD实验中观测到的突发事件（包括频率上移和复杂径向模结构）由何种机制驱动？
• RQ3 如2/3和3/4 EPM等二次模如何在突发阶段促进1/1 EIC的向内传播？
• RQ4 压强剖面平坦化与恢复在RIC稳定化及EIC突发后重新激活过程中的作用是什么？
• RQ5 FAR3d代码在多大程度上能再现实验测量的1/1 EIC频率、振幅和模结构？

主要发现

• 模拟再现了频率为9.4 kHz的1/1 EIC突发事件，与实验测量结果一致。
• 在突发期间，2/2和3/3泛音被激发，与在r/a = [0.6, 0.85]范围观测到的复杂径向模结构相符。
• 1/1 EIC的向内传播由与3/4和2/3 EPMs的非线性耦合驱动，导致突发期间共振重叠。
• 稳定化阶段由EIC引起的压强梯度平坦化消失后压强梯度的恢复触发，使RIC得以重新激活。
• EP β阈值对1/1 EIC激发的再现准确，且随着β升高，模振幅和宽度均增加。
• 该模型成功捕捉了完整的四阶段演化过程：RIC起始（阶段I）、EIC增长（阶段II）、突发伴随频率跳跃（阶段III），以及岛结构形成与RIC再激活（阶段IV）。