[论文解读] Gyrokinetic Continuum Simulation of Turbulence in Open-Field-Line Plasmas

E. L. Shi|arXiv (Cornell University)|Aug 24, 2017

Magnetic confinement fusion research参考文献 184被引用 12

一句话总结

本论文提出了首个能够模拟托卡马克偏滤器层中开放磁力线上的湍流的 gyrokinetic 连续性模拟代码，采用能量守恒的不连续伽辽金方法，并引入了新型导体鞘层边界条件。主要贡献在于开发了一种稳定且保持矩的算法，实现了对边界区域等离子体湍流（包括电磁效应和鞘层物理）的精确模拟，且在直线与螺旋形磁力线几何结构上均通过了验证。

ABSTRACT

The properties of the boundary plasma in a tokamak are now recognized to play a key role in determining the achievable fusion power and the lifetimes of plasma-facing components. Accurate quantitative modeling and improved qualitative understanding of the boundary plasma ultimately require five-dimensional gyrokinetic turbulence simulations, which have been successful in predicting turbulence and transport in the core. The additional challenges of boundary-plasma simulation necessitate the development of new gyrokinetic codes or major modifications to existing core gyrokinetic codes. In this thesis, we develop the first gyrokinetic continuum code capable of simulating plasma turbulence on open magnetic field lines, which is a key feature of a tokamak scrape-off layer. In contrast to prior attempts at this problem, we use an energy-conserving discontinuous Galerkin discretization in space. To model the interaction between the plasma and the wall, we design conducting-sheath boundary conditions that permit local currents into and out of the wall. We start by designing spatially one-dimensional kinetic models of parallel SOL dynamics and solve these systems using novel continuum algorithms. By generalizing these algorithms to higher dimensions and adding a model for collisions, we present results from the first gyrokinetic continuum simulations of turbulence on two types of open-field-line systems. The first simulation features uniform and straight field lines, such as found in some linear plasma devices. The second simulation is of a hypothetical model we developed of the NSTX scrape-off layer featuring helical field lines. These developments comprise a major step towards a gyrokinetic continuum code for quantitative predictions of turbulence and transport in the boundary plasma of magnetic fusion devices.

研究动机与目标

• 开发一种 gyrokinetic 连续性代码，以模拟开放磁力线等离子体中的湍流，特别是托卡马克偏滤器层中的湍流。
• 解决边界等离子体模拟中大振幅涨落、电磁效应和复杂磁几何结构带来的挑战。
• 实施能量守恒和矩守恒的数值格式，以实现长时间模拟的高精度。
• 通过允许与壁面局部交换电流的导体鞘层边界条件，对等离子体-壁面相互作用进行建模。
• 在代表线性装置和 NSTX 类几何结构的直线与螺旋形磁力线构型上对代码进行验证。

提出的方法

• 采用不连续伽辽金有限元方法进行空间离散化，确保相空间中的能量守恒和高阶精度。
• 开发了允许向壁面注入和提取局部电流的导体鞘层边界条件，以模拟等离子体-壁面相互作用。
• 在速度空间中引入指数加权多项式，相比标准多项式基底，显著提升了分辨率效率。
• 提出了一种改进的不连续伽辽金算法，可保持分布函数的所有矩，纠正了标准实现中的缺陷。
• 将一维平行偏滤器层动力学的动能模型推广至高维，并引入碰撞模型以实现更真实的湍流模拟。
• 采用 gyrokinetic 框架，对离子和电子动力学采用线性化方程，通过改进的平行电流模型引入鞘层电流效应。

实验结果

研究问题

• RQ1 能否利用 gyrokinetic 连续性代码准确模拟托卡马克偏滤器层中开放磁力线上的湍流？
• RQ2 导体鞘层边界条件如何影响开放磁力线系统中湍流的稳定性和动力学特性？
• RQ3 与标准多项式相比，速度空间中的指数加权多项式在多大程度上提升了离散化的分辨率效率？
• RQ4 改进的不连续伽辽金算法是否能保持分布函数的矩？为何这对长时间模拟至关重要？
• RQ5 鞘层电流对开放磁力线构型中湍流模的增长率有何影响？

主要发现

• 首次成功实现了使用具有能量守恒特性的不连续伽辽金方法的 gyrokinetic 连续性模拟，用于开放磁力线上的湍流。
• 导体鞘层边界条件实现了对等离子体-壁面电流交换的合理建模，这对偏滤器层物理至关重要。
• 改进的不连续伽辽金算法成功保持了分布函数的所有矩，而标准实现则无法做到这一点。
• 速度空间中的指数加权多项式显著提升了分布函数的分辨率效率，尤其在梯度较大的区域表现更优。
• 在强鞘层电流极限（εₜ ≫ 1）下，湍流模的增长率降低了 √εₜ 分之一，表明存在显著的稳定化效应。
• 在直线与螺旋形磁力线几何结构上的模拟结果证实了该代码的鲁棒性及其在真实托卡马克边界构型中的适用性。