[论文解读] Verification and validation of gyrokinetic and kinetic-MHD simulations for internal kink instability in DIII-D tokamak

G. Brochard, Jian Bao+15 more|arXiv (Cornell University)|Sep 20, 2021

Magnetic confinement fusion research参考文献 59被引用 2

一句话总结

本研究验证并确认了在DIII-D托卡马克放电#141216中，对内脆性扭转向不稳定性（internal kink instability）的gyrokinetic（GTC）和kinetic-MHD（GAM-solver, M3D-C1-K, NOVA-K, XTOR-K）模拟结果。当抑制了动力学效应时，所有代码在增长率和模结构方面均表现出极好的一致性，且在将q=1面位置调整至实验不确定度范围内后，模拟得到的径向模结构与电子回旋辐射测量结果高度吻合。

ABSTRACT

Verification and validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from the same equilibrium reconstruction of the DIII-D shot #141216, these codes exhibit excellent agreement for the growth rate and mode structure of the internal kink mode when all kinetic effects are suppresed. The simulated radial mode structures agree quantitatively with the electron cyclotron emission measurement after adjusting, within the experimental uncertainty, the safety factor q=1 flux-surface location in the equilibrium reconstruction. Compressible magnetic perturbations strongly destabilize the kink, while poloidal variations of the equilibrium current density reduce the growth rate of the kink. Furthermore, kinetic effects of thermal ions are found to decrease the kink growth rate in kinetic-MHD simulations, but increase the kink growth rate in gyrokinetic simulations, due to the additional drive of the ion temperature gradient and parallel electric field. Kinetic thermal electrons are found to have negligible effects on the internal kink instability.

研究动机与目标

• 对真实托卡马克实验中电流驱动内脆性扭转向不稳定性进行综合的验证与确认（V&V），针对gyrokinetic与kinetic-MHD代码。
• 使用DIII-D放电#141216的相同平衡条件，对多个代码的线性模拟结果进行基准测试。
• 将模拟得到的增长率与模结构与电子回旋辐射（ECE）测量结果进行验证。
• 研究压缩性磁扰动、极向电流变化以及热离子与电子动力学效应对该扭转向不稳定性的影响。

提出的方法

• 所有代码均采用DIII-D放电#141216的相同平衡重构结果，确保磁几何结构与离子分布一致。
• 所有代码在相同初始条件与模型参数下执行线性全局模拟。
• 在相同条件下（包括抑制动力学效应）比较各代码的增长率与径向模结构。
• 通过调整平衡重构中q=1磁面的径向位置，使其与实验ECE测量结果匹配。
• 系统性地改变物理效应——包括压缩性磁扰动、极向电流变化以及热离子/电子动力学效应——以分离其对不稳定性增长率的影响。
• 采用多种代码：GTC（gyrokinetic）、GAM-solver、M3D-C1-K、NOVA-K与XTOR-K（kinetic-MHD），各自使用不同的数值方法与粒子模型。

实验结果

研究问题

• RQ1 当抑制动力学效应时，gyrokinetic与kinetic-MHD代码在内脆性扭转向模的增长率与模结构方面是否产生一致结果？
• RQ2 在将q=1面位置调整至实验不确定度范围内后，模拟得到的径向模结构在多大程度上与电子回旋辐射测量结果吻合？
• RQ3 压缩性磁扰动如何影响内脆性扭转向模的稳定性？
• RQ4 平衡电流密度的极向变化对扭转向模增长率有何影响？
• RQ5 热离子与电子的动力学效应在不同模拟框架中对内脆性扭转向模增长率的定量影响如何？

主要发现

• 当抑制动力学效应时，五种代码——GTC、GAM-solver、M3D-C1-K、NOVA-K与XTOR-K——在增长率与模结构方面表现出极佳的一致性，证实了代码间的一致性。
• 在将q=1磁面位置调整至实验不确定度范围内后，模拟得到的径向模结构与电子回旋辐射测量结果在定量上完全吻合。
• 压缩性磁扰动显著加剧内脆性扭转向模的不稳定性，与解析理论及以往gyrokinetic研究结果一致。
• 平衡电流密度的极向非均匀性会降低扭转向模的增长率，表明电流剖面非均匀性具有稳定作用。
• 在kinetic-MHD模拟中，热离子动力学效应会降低增长率；而在gyrokinetic模拟中，由于离子温度梯度与平行电场的额外驱动作用，增长率反而增加。
• 在两种模拟框架中，热捕获电子对内脆性扭转向不稳定性的影响均可忽略。