[論文レビュー] Trinity: A Unified Treatment of Turbulence, Transport, and Heating in Magnetized Plasmas

M. Barnes|ArXiv.org|Jan 17, 2009

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

ひとこと要約

この論文では、複数のギャロアクティブラインシミュレーションを自己無撞撃的輸送方程式と結合することで、磁化プラズマにおける乱流、輸送、加熱を統合する新しい計算フレームワーク、Trinityを紹介する。適応的衝突性、速度空間診断、保存則を満たす衝突演算子を用いることで、従来のグローバルシミュレーションのコストの一部にとどまる第一原理的シミュレーションが可能となり、ITER関連プラズマの効率的かつ高精度なシミュレーションが実現される。

ABSTRACT

To faithfully simulate ITER and other modern fusion devices, one must resolve electron and ion fluctuation scales in a five-dimensional phase space and time. Simultaneously, one must account for the interaction of this turbulence with the slow evolution of the large-scale plasma profiles. Because of the enormous range of scales involved and the high dimensionality of the problem, resolved first-principles global simulations are very challenging using conventional (brute force) techniques. In this thesis, the problem of resolving turbulence is addressed by developing velocity space resolution diagnostics and an adaptive collisionality that allow for the confident simulation of velocity space dynamics using the approximate minimal necessary dissipation. With regard to the wide range of scales, a new approach has been developed in which turbulence calculations from multiple gyrokinetic flux tube simulations are coupled together using transport equations to obtain self-consistent, steady-state background profiles and corresponding turbulent fluxes and heating. This approach is embodied in a new code, Trinity, which is capable of evolving equilibrium profiles for multiple species, including electromagnetic effects and realistic magnetic geometry, at a fraction of the cost of conventional global simulations. Furthermore, an advanced model physical collision operator for gyrokinetics has been derived and implemented, allowing for the study of collisional turbulent heating, which has not been extensively studied. To demonstrate the utility of the coupled flux tube approach, preliminary results from Trinity simulations of the core of an ITER plasma are presented.

研究の動機と目的

• 磁化プラズマにおけるマルチスケールで五次元的な乱流をシミュレートする計算的課題、特にITERやその他の核融合装置を対象とする。
• 乱流フラックスと加熱を考慮しながら、大規模なプラズマプロファイルを自己無撞撃的に進化させる手法を開発すること。
• 特に乱流加熱が未十分に探求されてきたが、乱流プラズマにおける衝突効果の高精度なシミュレーションを可能とすること。
• ブルートフォースグローバルシミュレーションに代えて、結合されたフラックスチューブ－輸送アプローチを用いることで、計算コストを低減すること。
• 適応的散乱と分解能診断を通じて、速度空間ダイナミクスにおける数値的忠実性を保証すること。

提案手法

• ギャロアクティブラインシミュレーションにおける最小限の必要十分な散乱を保証するため、積分誤差推定とスペクトル法を用いた速度空間分解能診断を開発した。
• 局所的なプラズマ状態に応じて衝突周波数を適応的に調整する手法を実装し、過剰な数値的拡散を回避しながら精度を維持した。
• 局所モーメントを保存し、H定理を満たす新しい保存則を満たす線形化Fokker-Planck衝突演算子を考案した。
• 複数のギャロアクティブラインフラックスチューブシミュレーションを輸送方程式で結合し、複数のイオンおよび電子種の自己無撞撃的安定状態プロファイルを進化させた。
• 時間的に陰解法で、Sherman-Morrison加速を用いたコンact有限差分法を用いて、結合された輸送方程式を効率的に解いた。
• 電磁効果と現実的な磁場幾何学的構造をTrinityコードフレームワークに統合し、現実的なコアプラズマモデリングを可能にした。

実験結果

リサーチクエスチョン

• RQ1 磁化プラズマにおいて、複数の空間的・時間的スケールにわたって乱流、輸送、加熱を自己無撞撃的にシミュレートする方法は何か？
• RQ2 衝突は乱流加熱にどのような役割を果たすのか？ そして、ギャロアクティブラインシミュレーションにおいてこの現象を高精度にモデル化するにはどうすればよいか？
• RQ3 結合されたフラックスチューブ－輸送アプローチは、グローバルシミュレーションのコストの一部にとどめる範囲で、高精度な安定状態のプラズマプロファイルを達成できるか？
• RQ4 安定性と保存則の性質を維持しつつ、最小限の数値的散乱で速度空間ダイナミクスをどのように解像できるか？
• RQ5 ギャロアクティブラインコードにおける衝突演算子の最適な数値的処理は何か？ 物理的一致性と正確性を保証するには？

主な発見

• Trinityフレームワークは、結合されたフラックスチューブシミュレーションを用いて、ITERに類似したプラズマのコア領域におけるプラズマプロファイル（密度、温度）の自己無撞撃的進化を成功裏にシミュレートした。
• 適応的衝突周波数と速度空間診断により、人工的散乱を最小限に抑えつつ、キネティック効果の高精度な解像が可能となった。
• 新規のモデル衝突演算子は局所モーメントを保存し、H定理を満たしており、衝突的乱流加熱の信頼できる研究を可能にした。
• 結合アプローチにより、ブルートフォースグローバルシミュレーションと比較して計算コストを顕著に低減したが、物理的忠実性は維持された。
• ITERコアの初期Trinityシミュレーションは、安定した自己無撞撃的プロファイルを示し、現実的な乱流フラックスと熱負荷を再現した。
• コンパクト有限差分法にSherman-Morrison加速を組み合わせた実装により、輸送方程式の効率的時間的陰解法が可能となった。