[论文解读] Gyrokinetic simulations of fusion plasmas using a spectral velocity space representation

Joseph Thomas Parker|arXiv (Cornell University)|Mar 15, 2016

Magnetic confinement fusion research参考文献 155被引用 3

一句话总结

本文提出SpectroGK，一种高性能的磁流体力学模拟代码，采用基于汉克尔-埃米特多项式的谱速度空间表示，可精确守恒离散自由能，并仅用十个速度变量高效模拟朗道阻尼。该方法实现完全本地化的电磁场计算，消除了处理器间通信，加速了线性不稳定性搜索；模拟结果表明，非线性级联会抑制湍流等离子体中的线性相混叠，导致流体样行为。

ABSTRACT

Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operation, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the gyrokinetic equations, a nonlinear integro-differential system for evolving the particle distribution functions in a five-dimensional position and velocity space, and the consequent electromagnetic field. Due to the high dimensionality, simulations of small reactor sections require hundreds of thousands of CPU hours on High Performance Computing platforms. We develop a Hankel-Hermite spectral representation for velocity space that exploits structural features of the gyrokinetic system. The representation exactly conserves discrete free energy in the absence of explicit dissipation, while our Hermite hypercollision operator captures Landau damping with few variables. Calculation of the electromagnetic fields becomes purely local, eliminating inter-processor communication in, and vastly accelerating, searches for linear instabilities. We implement these ideas in SpectroGK, an efficient parallel code. Turbulent fusion plasmas may dissipate free energy through linear phase mixing to fine scales in velocity space, as in Landau damping, or through a nonlinear cascade to fine scales in physical space, as in hydrodynamic turbulence. Using SpectroGK to study saturated electrostatic drift-kinetic turbulence, we find that the nonlinear cascade suppresses linear phase mixing at energetically-dominant scales, so the turbulence is fluid-like. We use this observation to derive Fourier-Hermite spectra for the electrostatic potential and distribution function, and confirm these spectra with simulations.

研究动机与目标

• 开发一种基于汉克尔-埃米特多项式的谱速度空间表示，以在磁流体力学模拟中精确守恒离散自由能。
• 通过最少数量的速度空间变量，实现高效且精确的朗道阻尼建模。
• 通过完全本地化电磁场计算，消除处理器间通信，加速线性不稳定性搜索。
• 研究非线性级联在湍流聚变等离子体中抑制线性相混叠的作用。
• 推导并验证漂移-动力学湍流中静电势和分布函数的傅里叶-埃米特谱。

提出的方法

• 在速度空间中采用汉克尔-埃米特谱表示，以利用磁流体力学系统中平流、平均化和碰撞项的结构特性。
• 采用埃米特超碰撞算子，仅用十个速度变量即可捕捉朗道阻尼，实现高精度与自由能精确守恒。
• 通过谱结构实现电磁场的完全本地化计算，消除并行模拟中的全局通信。
• 开发SpectroGK，一种并行代码，结合时间积分、非线性项去卷积和速度空间耗散的递推关系以保证稳定性。
• 应用离散傅里叶-埃米特谱表示，分析静电势和分布函数的相空间谱。
• 基于汉克尔变换的的速度空间网格用于表示垂直速度空间，实现磁流体力学项的高效计算。

实验结果

研究问题

• RQ1 基于汉克尔-埃米特多项式的谱速度空间表示能否在磁流体力学模拟中实现精确的自由能守恒？
• RQ2 在谱表示下，仅用十个速度空间变量在多大程度上能准确捕捉朗道阻尼？
• RQ3 物理空间中的非线性级联是否抑制速度空间中的线性相混叠，导致湍流聚变等离子体中出现流体样行为？
• RQ4 静电势和分布函数的傅里叶-埃米特谱在饱和漂移-动力学湍流中的标度规律是什么？
• RQ5 基于速度空间结构推导的流体模型能否在显著降低分辨率的情况下准确复现动力学湍流？

主要发现

• 在物理空间中的非线性级联完全抑制了能量主导尺度上的线性相混叠，导致饱和静电漂移-动力学湍流中出现流体样行为。
• 埃米特谱表示仅用十个速度变量即可捕捉朗道阻尼，且在无显式耗散时实现高精度与精确的自由能守恒。
• 由于谱结构，电磁场计算变为完全本地化，消除了处理器间通信，加速了线性不稳定性搜索。
• 静电势和分布函数的傅里叶-埃米特谱遵循与理论预测一致的幂律标度，经SpectroGK模拟验证。
• 自由能耗散的超碰撞平台在不同分辨率和超碰撞参数下均保持稳健，且在中等分辨率下实现收敛。
• 净自由能向精细平行速度空间尺度的转移被抑制，表明仅使用少数埃米特模态即可构建降阶流体模型，实现约O(10^4)的分辨率降低。