[论文解读] Electromagnetic gyrokinetic instabilities in STEP

Daniel Kennedy, M. Giacomin+5 more|arXiv (Cornell University)|Jul 4, 2023

Magnetic confinement fusion researchPhysics and Astronomy被引用 3

一句话总结

本研究针对STEP球形托卡马克设计中的电磁微不稳定性，开展了线性 gyrokinetic 分析，识别出由离子温度梯度和束缚电子驱动的混合型动能 ballooning 模式（KBM）在离子横向尺度上占主导地位。碰撞型微撕裂模（MTM）虽为次主导，但具有显著影响，且具有独特的径向尺度。通过多代码基准测试和分辨率收敛性验证，结果具备可靠性，为后续非线性研究（在附录论文中）提供了信心支持。

ABSTRACT

We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production, https://step.ukaea.uk). We examine a range of flux surfaces between the deep core and the pedestal top for two candidate flat-top operating points of the prototype device. Local linear gyrokinetic analysis is performed to determine the type of microinstabilities that arise under these reactor-relevant conditions. We find that the equilibria are dominated at ion binormal scales by a hybrid version of the Kinetic Ballooning Mode (KBM) instability that has significant linear drive contributions from the ion temperature gradient and from trapped electrons, while collisional Microtearing Modes (MTMs) are sub-dominantly also unstable at similar binormal scales. The hybrid-KBM and MTM exhibit very different radial scales. We study the sensitivity of these instabilities to physics parameters, and discuss potential mechanisms for mitigating them. The results of this investigation are compared to a small set of similar conceptual reactor designs in the literature. A detailed benchmark of the linear results is performed using three gyrokinetic codes; alongside extensive resolution testing and sensitivity to numerical parameters providing confidence in the results of our calculations, and paving the way for detailed nonlinear studies in a companion article.

研究动机与目标

• 评估高β、反应堆尺度、小长径比球形托卡马克（STEP）平衡态中电磁微不稳定的特性。
• 在从核心到边缘的磁面范围内，识别两种候选STEP运行点中主导与次主导不稳定性。
• 评估未来使用多种gyrokinetic代码进行非线性模拟时的数值收敛性与分辨率需求。
• 在三个gyrokinetic代码（GS2、CGYRO、GENE）之间进行基准测试，验证线性稳定性预测结果。
• 通过确定分辨率与参数敏感性，为后续非线性湍流模拟研究（论文II）奠定基础。

提出的方法

• 使用GS2代码对两种候选STEP平衡态（STEP-EC-HD与STEP-EB-HD）进行局部线性gyrokinetic分析。
• 在多个磁面（Ψₙ = 0.49 至 0.95）上执行模拟，覆盖从核心到边缘区域。
• 采用三种gyrokinetic代码——GS2、CGYRO与GENE——进行交叉验证与基准测试，比较线性增长率与模态频率。
• 在径向、平行方向（ballooning空间）与速度空间网格上执行严格的分辨率收敛性测试，以确保数值稳健性。
• 研究对碰撞模型（Lorentz、Sugama与完整线性化Fokker-Planck）及粒子种类数（2–5）的敏感性，尤其关注增长率稳定性。
• 通过详细参数扫描（nθ、nperiod、npass、nε）确定非线性模拟所需的最小分辨率要求。

实验结果

研究问题

• RQ1 在高β、反应堆相关、STEP平衡态中，离子横向尺度上的主导电磁微不稳定性是什么？
• RQ2 在磁面层级中，混合型动能Ballooning模（KBM）与碰撞型微撕裂模（MTM）在增长率与径向尺度上的对比如何？
• RQ3 为实现准确的线性稳定性预测以支持非线性模拟，需要多高的数值分辨率？
• RQ4 线性增长率与模态频率对碰撞算子选择及等离子体粒子种类数的敏感性如何？
• RQ5 多代码基准测试（GS2、CGYRO、GENE）在多大程度上证实了不稳定性预测的稳健性？

主要发现

• 在离子横向尺度（k_yρ_i ≲ 1）上，平衡态主要受一种混合型动能Ballooning模（KBM）主导，其线性驱动来自离子温度梯度与束缚电子的共同作用。
• 碰撞型微撕裂模（MTM）为次主导，但在相似横向尺度上仍不稳定，其径向尺度明显宽于混合-KBM。
• 当k_yρ_s ≤ 0.4时，nθ ≥ 32 即可实现分辨率收敛；但当k_yρ_s > 0.4时，需nθ ≥ 64，表明高波数模态需要更高分辨率。
• 速度空间分辨率对增长率与频率的依赖性极弱，在npass与nε的测试范围内均表现出收敛性。
• 线性谱对碰撞模型选择的敏感性较弱，但使用CGYRO中的Sugama算子时，增长率较Lorentz或完整Fokker-Planck模型低约20%。
• 在GS2、CGYRO与GENE之间进行的多代码基准测试，证实了增长率与模态频率的强一致性，验证了结果在非线性模拟中的可靠性。