[论文解读] Magnetic Fluctuations in Gyrokinetic Simulations of Tokamak Scrape-Off Layer Turbulence

Noah Mandell|arXiv (Cornell University)|Mar 30, 2021

Magnetic confinement fusion research参考文献 211被引用 7

一句话总结

本论文开发了首个针对托卡马克偏滤器层（SOL）中开放磁力线的全f电磁磁流体力学模拟能力，采用不连续伽辽金（DG）格式，确保能量守恒并避免安培抵消问题。该方法成功实现了对电磁涨落的稳健模拟，并揭示了其对NSTX类SOL条件下团块动力学与输运的影响。

ABSTRACT

Understanding turbulent transport physics in the tokamak edge and scrape-off layer (SOL) is critical to developing a successful fusion reactor. The dynamics in these regions plays a key role in achieving high fusion performance by determining the edge pedestal that suppresses turbulence in H-mode. Additionally, the survivability of a reactor is set by the heat load to the vessel walls, making it important to understand turbulent spreading of heat as it flows along open magnetic field lines in the SOL. Large-amplitude fluctuations, magnetic X-point geometry, and plasma interactions with material walls make simulating turbulence in the edge/SOL more challenging than in the core region, necessitating specialized gyrokinetic codes. Further, the inclusion of electromagnetic effects in gyrokinetic simulations that can handle the unique challenges of the boundary plasma is critical to the understanding of phenomena such as the pedestal and ELMs, for which electromagnetic dynamics are expected to be important. In this thesis, we develop the first capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines. This is an important step towards comprehensive electromagnetic gyrokinetic simulations of the coupled edge/SOL system. By using a continuum full-f approach via an energy-conserving discontinuous Galerkin (DG) discretization scheme that avoids the Ampere cancellation problem, we show that electromagnetic fluctuations can be handled in a robust, stable, and efficient manner in the gyrokinetic module of the Gkeyll code. We then present results which roughly model the scrape-off layer of the National Spherical Torus Experiment (NSTX), and show that electromagnetic effects can affect blob dynamics and transport. We also formulate the gyrokinetic system in field-aligned coordinates for modeling realistic edge and scrape-off layer geometries in experiments.

研究动机与目标

• 开发一种针对托卡马克边缘和偏滤器层（SOL）的稳健、能量守恒的电磁磁流体力学模拟能力，该区域存在开放磁力线和强烈的磁涨落，使湍流建模变得复杂。
• 通过实现连续全f不连续伽辽金（DG）空间离散化方案，克服磁流体力学模拟中的安培抵消问题。
• 在真实SOL几何构型下，实现对电磁效应在湍流输运和团块动力学中作用的精确建模，这对理解堆垛形成和边缘局部模（ELMs）至关重要。
• 在磁场对齐坐标系中构建磁流体力学系统，以更好地表征实验中的边缘等离子体构型。
• 通过开发一种新型DG算法，确保在存在电磁效应时，分布函数保持正定性，同时保持质量、动量和能量守恒，从而保证数值鲁棒性。

提出的方法

• 采用连续全f方法，结合不连续伽辽金（DG）空间离散化，以处理开放磁力线上磁流体力学模拟中的电磁涨落。
• 实现一种能量守恒的DG格式，避免电磁磁流体力学代码中已知的安培抵消问题，即一种不稳定性。
• 使用磁场对齐坐标系，以模拟包括X点构型和等离子体-壁相互作用在内的真实边缘和SOL等离子体几何构型。
• 开发一种分布函数的正性保持算法，确保在存在电磁效应时仍能保持质量、动量和能量守恒。
• 将该方法应用于简化NSTX类偏滤器层构型的模拟，以研究电磁效应在团块动力学和输运中的作用。
• 采用简化的Dougherty模型集成碰撞算子，并探索未来向非线性福克-普朗克模型和GPU加速的扩展。

实验结果

研究问题

• RQ1 在托卡马克边缘和偏滤器层（SOL）中，如何稳健且稳定地模拟电磁效应，该区域存在开放磁力线和大尺度涨落？
• RQ2 在真实的NSTX类条件下，电磁涨落在SOL中对团块动力学和横场输运的影响程度如何？
• RQ3 具有能量守恒和正性保持特性的不连续伽辽金（DG）格式，能否在开放磁力线上实现稳定全f电磁磁流体力学模拟？
• RQ4 电磁效应在堆垛输运和边缘局部模（ELMs）中起什么作用？它们如何影响堆垛与SOL之间的耦合？
• RQ5 如何将数值框架扩展以包含真实碰撞算子和GPU加速，以支持高性能计算？

主要发现

• 所开发的基于DG格式的全f电磁磁流体力学框架成功避免了安培抵消问题，实现了在开放磁力线上的稳定、能量守恒模拟。
• 模拟结果表明，电磁效应显著影响NSTX类偏滤器层中的团块动力学和湍流输运，表明其在边缘等离子体行为中具有关键作用。
• 正性保持DG算法在强电磁涨落存在下仍能保持数值鲁棒性，同时实现质量、动量和能量的守恒。
• 磁场对齐坐标系的实现使得真实边缘等离子体几何构型（包括X点构型和等离子体-壁相互作用）得以精确建模。
• 该框架可扩展以集成更先进的碰撞模型和GPU加速，目前已在NVIDIA GPU上完成求解核的CUDA移植工作。
• 该能力实现了对耦合堆垛/SOL系统的自洽研究，这对理解聚变反应堆中高约束模（H-mode）性能和功率排散挑战至关重要。