[論文レビュー] Some numerical aspects of the conservative PSM scheme in a 4D drift-kinetic code

Jean‐Philippe Braeunig, Nicolas Crouseilles+4 more|arXiv (Cornell University)|Mar 9, 2013

Magnetic confinement fusion research参考文献 14被引用数 7

ひとこと要約

本稿は、核融合プラズマシミュレーションの4次元ドリフト運動的コードにおける保存的ピecewise多項式再構成スキーム（PSM）の数値的安定性および保存性特性を調査する。有限体積形式のPSMは、特に長時間にわたる乱流状態下で、方向分割型と比較して最大原理を保ち、物理的に意味のない振動をより効果的に回避することを示している。スロープ制限付きスプライン（SLS）リミッターは、制御された拡散を加えることでさらに安定性を向上させる。

ABSTRACT

The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, kinetic Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy simulations because a 6D dimensional problem has to be solved, even if reduced to a 5D in so called gyrokinetic models. Accurate schemes, parallel algorithms need to be designed to bear these simulations. This paper describes the numerical studies to improve robustness of the conservative PSM scheme in the context of its development in the GYSELA code. In this paper, we only consider the 4D drift-kinetic model which is the backbone of the 5D gyrokinetic models and relevant to build a robust and accurate numerical method.

研究の動機と目的

• GYSELA 5次元ジェオキンティックコードにおける保存的PSMスキームの数値的安定性を向上させるために、低次元化された4次元ドリフト運動的モデルにおける挙動を調査すること。
• 離散的レベルで最大原理が保たれることを保証し、シミュレーションのクラッシュを引き起こす可能性のある物理的に意味のない過剰応答や不足応答を回避すること。
• 方向分割型と有限体積型のPSMバリエーションを比較し、長時間積分下での保存性特性および安定性を評価すること。
• SLSリミッターが保存則を破らないまま、不適切な振動を効果的に抑制する効果を評価すること。
• 質量およびエネルギー保存を維持しながら、乱流構造を正確に捉えるスキームの性能を検証すること。

提案手法

• イオン乱流の完全な5次元ジェオキンティック問題を簡略化した、円筒座標系(r, θ, z, v∥)における4次元ドリフト運動的モデルを用いる。
• 保存的PSMスキームを方向分割（DS）形式および有限体積（FV）形式で実装し、後者ではフラックスの発散なしを保証する。
• グリッドスケール未満での振動を制御するため、PSMにスロープ制限付き再構成（SLS）を適用する。
• 初期ガウスパルスを用いたベンチマークで、分布関数の極値の時間変化を追跡し、最大原理の保持を評価する。
• 非線形飽和の初期段階（t = 1800）と完全に発達した乱流状態（t = 4400）をそれぞれシミュレーションする。
• PSM-DS、PSM-FV、SLS-DS、SLS-FV、および参照用のBSLスキームとの間で、精度、振動レベル、極値の保持性を比較する。

実験結果

リサーチクエスチョン

• RQ1 長時間積分下で、PSMスキームの有限体積形式は方向分割型と比較して、最大原理をどの程度効果的に保っているか？
• RQ2 SLSリミッターは、保存則を破らないまま、PSMスキームにおける不適切な振動をどの程度効果的に抑制しているか？
• RQ3 乱流状態下で、保存的PSMスキームは長時間積分の間、分布関数の有界な極値を維持できるか？
• RQ4 有限体積形式の発散なし性質が、不安定性に対して解を安定化させる役割を果たすか？
• RQ5 異なる数値的実装（PSM-DS、PSM-FV、SLS-FV）は、物理的に意味のないアーチファクトを生じさせず、小スケール乱流構造をどの程度正確に解像できるか？

主な発見

• t = 1800におけるPSMの有限体積（FV）形式は、初期の極値（最小 = 0.331、最大 = 0.4187）を保持しているが、PSM方向分割型ではすでにずれが生じており（最小 = 0.3086、最大 = 0.4427）、その差が顕著である。
• t = 4400において、PSM-DSスキームは顕著な物理的に意味のない振動を生じ、シミュレーションのクラッシュリスクを伴うが、SLS-FVスキームは許容可能な範囲内に極値を維持している（最小 = 0.3263、最大 = 0.4221）。
• SLSリミッターは制御された拡散を加えることで振動を効果的に抑制し、全体の解構造に影響を与えることなく安定性を向上させる。
• 有限体積形式は、特に勾配が大きく複雑な流れ構造を示す領域で、方向分割型よりも優れた数値的安定性を示す。
• 参照用のBSLスキームも再構成ステップのため最大原理を破っており、極値のずれが生じている（t = 1800で最小 = 0.3124、最大 = 0.4430）が、PSM-DSに比べて振動制御性能が優れている。
• 有限体積スキームの発散なし性質は、特にr方向の中間半径付近で不安定性を防ぐ上で極めて重要であり、SLS-DSでは物理的に意味のないリング状の振動が観測されるが、SLS-FVではそのような振動は見られない。