[論文レビュー] 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）におけるオープン磁場線を対象として、エネルギー保存を保証し、アームペルキャンセレーションを回避する不連続ガラーキン（DG）法を用いた、最初の完全f電磁気的ギャロアクイネティクスシミュレーション能力を開発した。電磁気的揺動の安定なシミュレーションを示し、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スキームを実装し、電磁気的ギャロアクイネティクスコードにおける既知の不安定性を解消する。
• 実際のエッジおよびSOLプラズマ幾何構造（Xポイント構成やプラズマ-壁相互作用を含む）をモデル化するために、場磁気座標系を用いる。
• 電磁気的効果が存在する中でも質量、運動量、エネルギーを保存することができる正値性を保つアルゴリズムを、分布関数に対して開発する。
• ブロブダイナミクスおよび輸送への電磁気的効果を調査するために、簡略化されたNSTXに類似したスクラップ・オフ層構造を用いてシミュレーションを実施する。
• 衝突演算子には簡略化されたドーキンジー・モデルを適用し、将来的には非線形フォッカープランクモデルやGPUアクセcelerationへの拡張を検討する。

実験結果

リサーチクエスチョン

• RQ1 オープン磁場線と大きな揺動が支配するトカマクエッジおよびスクラップ・オフ層（SOL）において、電磁気的効果を安定かつ強固にシミュレートする方法は何か？
• RQ2 現実的なNSTXに類似した条件下で、電磁気的揺動はSOLにおけるブロブダイナミクスおよび横方向輸送にどの程度の影響を及ぼすか？
• RQ3 エネルギー保存および正値性保持を備えた不連続ガラーキン（DG）スキームを用いることで、オープン磁場線における安定な完全f電磁気的ギャロアクイネティクスシミュレーションが可能になるか？
• RQ4 電磁気的効果はペディメント輸送およびエッジ局在モード（ELMs）において果たす役割は何か？また、ペディメントとSOLのカップリングにどのように影響を与えるか？
• RQ5 数値フレームワークを、現実的な衝突演算子およびGPUアクセcelerationを含めるようにどのように拡張できるか？

主な発見

• DG法を用いた完全f電磁気的ギャロアクイネティクスフレームワークは、アームペルキャンセレーション問題を効果的に回避し、オープン磁場線における安定でエネルギー保存型のシミュレーションを可能にした。
• 電磁気的効果が、シミュレートされたNSTXに類似したスクラップ・オフ層におけるブロブダイナミクスおよび乱流輸送に顕著な影響を与えることが示され、エッジプラズマ挙動においてその重要性が明らかになった。
• 正値性を保つDGアルゴリズムは、強い電磁気的揺動が存在する中でも、質量、運動量、エネルギーの保存を維持しながら、数値的安定性を確保した。
• 場磁気座標系の導入により、Xポイント構成やプラズマ-壁相互作用を含む現実的なエッジプラズマ幾何構造の正確なモデル化が可能になった。
• フレームワークは、より高度な衝突モデルやGPUアクセcelerationへの拡張が可能であり、NVIDIA GPU向けにソルバーコアをCUDAに移植する初期の進捗が得られた。
• この能力により、ペディメント/SOL系の自己一貫的かつ統合的な研究が可能となり、核融合炉における高閉じ込めモード（Hモード）性能およびパワー排出の課題を理解する上で不可欠となった。