[論文レビュー] 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を紹介する。この手法により、離散的自由エネルギーの正確な保存が可能となり、わずか10個の並進変数で効率的にランダム減衰をモデル化できる。本手法により、完全に局所的な電磁場計算が可能となり、プロセッサ間通信が不要となり、線形不安定性探索が高速化される。また、シミュレーションでは非線形スピンオーバーが乱流プラズマにおける線形位相混合を抑制し、流体的挙動を引き起こすことが明らかになった。

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.

研究の動機と目的

• ギャロアクティブラジエーションシミュレーションにおける離散的自由エネルギーの正確な保存を実現するスペクトル的並進空間表現の開発。
• 最小限の並進空間変数を用いて、正確かつ効率的なランダム減衰のモデル化を可能にする。
• スペクトル構造を活用することで、電磁場計算を完全に局所化し、並列シミュレーションにおけるグローバル通信を排除し、線形不安定性探索を加速する。
• 乱流核融合プラズマにおける非線形スピンオーバーが線形位相混合をどのように抑制するかを調査する。
• ドリフト・キネティック乱流における静電ポテンシャルおよび分布関数のフーリエ–エルミートスペクトルの導出と検証。

提案手法

• ギャロアクティブラジエーション系におけるストリーミング、ジロ平均化、衝突項の構造的特徴を活用するため、並進空間にハンケル–エルミートスペクトル表現を採用する。
• ランダム減衰をわずか10個の並進変数で捉えることができるエルミートハイパーキャスティング作用素を採用し、高い精度と自由エネルギー保存を達成する。
• スペクトル構造を活用することで、電磁場計算を完全に局所化し、並列シミュレーションにおけるグローバル通信を排除する。
• 時間積分、非線形項のデコンボリューション、並進空間における散逸を組み合わせた再帰的処理を用いて、安定性を確保する並列コードSpectroGKを開発する。
• 静電ポテンシャルおよび分布関数の位相空間スペクトルを分析するため、離散的フーリエ–エルミートスペクトル表現を適用する。
• 垂直並進空間を表すためにハンケル変換に基づく並進空間グリッドを採用し、ギャロアクティブラジエーション項の効率的計算を可能にする。

実験結果

リサーチクエスチョン

• RQ1 ハンケル–エルミート多項式に基づくスペクトル的並進空間表現は、ギャロアクティブラジエーションシミュレーションにおいて、自由エネルギーの正確な保存を達成できるか？
• RQ2 スペクトル表現においてわずか10個の並進空間変数で、ランダム減衰をどの程度正確に捉えることができるか？
• RQ3 物理空間における非線形スピンオーバーは、並進空間における線形位相混合を抑制し、乱流核融合プラズマで流体的挙動を引き起こすか？
• RQ4 飽和したドリフト・キネティック乱流における静電ポテンシャルおよび分布関数のフーリエ–エルミートスペクトルのスケーリング則は何か？
• RQ5 並進空間構造から導かれる流体的モデルは、顕著に低解像度で運動論的乱流を正確に再現できるか？

主な発見

• 物理空間における非線形スピンオーバーは、エネルギー的に支配的スケールで線形位相混合を完全に抑制し、飽和した静電的ドリフト・キネティック乱流において流体的挙動を生じさせる。
• エルミートスペクトル表現は、明示的な散逸が存在しない状況でも、わずか10個の並進変数でランダム減衰を高精度に捉え、自由エネルギーの正確な保存を達成する。
• スペクトル構造のおかげで電磁場計算が完全に局所化され、プロセッサ間通信が不要となり、線形不安定性探索が高速化される。
• 静電ポテンシャルおよび分布関数のフーリエ–エルミートスペクトルは、理論的予測と整合的なべき乗則スケーリングを示し、SpectroGKシミュレーションによって確認された。
• 自由エネルギー散逸のハイパーキャスティングプラトーは、さまざまな解像度およびハイパーキャスティングパラメータに対して安定しており、中程度の解像度で収束が観察された。
• 微細な並進空間スケールへの自由エネルギー移動の抑制は、数個のエルミートモードのみを用いた低次元流体モデルの可能性を示唆しており、解像度を約O(10^4)まで削減可能である。