[論文レビュー] Verification of Gyrokinetic codes: theoretical background and applications

N. Tronko, A. Bottino+4 more|arXiv (Cornell University)|Mar 22, 2017

Magnetic confinement fusion research参考文献 28被引用数 20

ひとこと要約

本稿では、ラグランジュ形式の変分法とコード間ベンチマークを組み合わせることで、ギロキンティックコードの検証を統一的かつ理論的・数値的に実現するフレームワークを提示する。ORB5（PIC）とGENE（Eulerian）コードがフルFLR補正を含む場合、線形モデルとして同等であることが示され、電磁的βスキャンおよびモード構造ベンチマークにおいて優れた一致が得られ、数値的整合性が検証されるとともに、微細構造に対する正確なポアソンソルバーの重要性が強調された。

ABSTRACT

In fusion plasmas the strong magnetic field allows the fast gyro-motion to be systematically removed from the description of the dynamics, resulting in a considerable model simplification and gain of computational time. Nowadays, the gyrokinetic (GK) codes play a major role in the understanding of the development and the saturation of turbulence and in the prediction of the subsequent transport. Naturally, these codes require thorough verification and validation. Here we present a new and generic theoretical framework and specific numerical applications to test the faithfulness of the implemented models to theory and to verify the domain of applicability of existing GK codes. For a sound verification process, the underlying theoretical GK model and the numerical scheme must be considered at the same time, which has rarely been done and therefore makes this approach pioneering. At the analytical level, the main novelty consists in using advanced mathematical tools such as variational formulation of dynamics for systematization of basic GK code's equations to access the limits of their applicability. The verification of numerical scheme is proposed via the benchmark effort. In this work, specific examples of code verification are presented for two GK codes: the multi-species electromagnetic ORB5 (PIC) and the radially global version of GENE (Eulerian). The proposed methodology can be applied to any existing GK code. We establish a hierarchy of reduced GK Vlasov-Maxwell equations implemented in the ORB5 and GENE codes using the Lagrangian variational formulation. At the computational level, detailed verifications of global electromagnetic test cases developed from the CYCLONE Base Case are considered, including a parametric $\beta$-scan covering the transition from ITG to KBM and the spectral properties at the nominal $\beta$ value.

研究の動機と目的

• ギロキンティックコードの厳密な二重検証フレームワークを確立し、理論的モデル検証と数値スキーム検証を統合する。
• 異なるコードで用いられる順序付けと近似を明確にすることで、理論的ギロキンティックモデルとコード実装との乖離を解消する。
• 電磁的テストケースにおける詳細なコード間ベンチマークを通じて、ORB5（PIC）およびGENE（Eulerian）コードの数値的整合性を検証する。
• 特に、ORB5における長波長ポアソンソルバーという数値的近似が、静電ポテンシャルの微細構造の解像度に与える影響を特定する。
• 微小不安定性遷移の正確なシミュレーションには、径方向グローバルモデルが不可欠であることを示す。特に有限βおよびρ∗効果が存在する場合に顕著である。

提案手法

• ラグランジュ変分法を用いてギロキンティックVlasov-Maxwell方程式を導出し、その背後にある理論的モデルとその適用限界を体系的に特定する。
• 同一の変分フレームワークを用いて、ORB5（PIC）とGENE（Eulerian）コードの簡略化された粒子運動およびMaxwell-Vlasovモデルを解析的に比較する。
• 2つのテストケースにおけるコード間ベンチマークを実施：固定トロイダルモード数n = 19のグローバル電磁的βスキャン、および固定βにおける径方向および磁層モード構造の詳細な比較。
• βスキャンでは線形成長率と実周波数を観測量とし、モード構造比較では静電ポテンシャルおよび平行磁場ポテンシャルの全磁層および径方向プロファイルを用いる。
• 両コードに同一の低β順序付けを適用して一貫性を確保し、GENEのフルFLRソルバーとORB5の長波長近似を対比する。
• グローバルおよびローカル（フラックスチューブ）GENE結果の比較を通じて有限サイズ（ρ∗）効果を定量化し、径方向構造効果によりモード遷移閾値が異なることが示された。

実験結果

リサーチクエスチョン

• RQ1 ORB5およびGENEコードに実装された理論的ギロキンティックモデルは、同一の変分フレームワークから導出された場合、どの程度一致するか？
• RQ2 数値的近似、特にポアソンソルバーが、静電ポテンシャルの微細構造の解像度にどのように影響を与えるか？
• RQ3 有限径方向サイズ（ρ∗）は、グローバル電磁的シミュレーションにおけるイオン温度勾配（ITG）モードとキネティックバルーンモード（KBM）の遷移にどのような影響を及ぼすか？
• RQ4 グローバル電磁的βスキャンにおいて、β値の範囲にわたるORB5およびGENEの線形成長率と周波数はどの程度一致するか？
• RQ5 2つのコード間における静電ポテンシャルおよび平行磁場ポテンシャルのモード構造にどのような差が生じるか、その原因は何か？

主な発見

• 同一のラグランジュ変分法を用いて導出した場合、ORB5およびGENEの理論的モデルは線形領域で同一であることが判明し、物理的整合性が確認された。
• n = 19のグローバル電磁的βスキャンにおいて、ORB5とGENEの両コードがβ範囲にわたり類似した線形成長率と実周波数を示し、優れた一致が観測された。
• グローバルとローカル（フラックスチューブ）シミュレーションの間で、ITGモードからKBMモードへの遷移が異なるβ値で発生することが判明し、有限サイズ（ρ∗ ≈ 1/182）効果が無視できないこと、かつコード依存的であることが示された。
• ORB5は長波長ポアソンソルバーを用いているが、GENEのフルFLRソルバーが捉える微細な径方向構造（特に有理数表面付近）を解像できない。
• ORB5ではこれらの欠落構造は低nでは成長率や周波数に顕著な影響を及ぼさないが、高トロイダルモード数（n ≳ 40）では著しくずれが生じ、フルFLRソルバーの必要性が強調された。
• ベンチマーク結果は、一貫した近似を用いる場合、ギロキンティックモデルの数値的実装が理論と忠実に一致することを確認した。また、核融合プラズマシミュレーションにおける信頼性の高い予測には、コード間検証が不可欠であることが示された。