[論文レビュー] Accuracy of unperturbed motion of particles in a gyrokinetic semi-Lagrangian code

Guillaume Latu, V. Grandgirard+10 more|arXiv (Cornell University)|Sep 3, 2012

Magnetic confinement fusion research参考文献 9被引用数 3

ひとこと要約

本稿は、ギャロキナetic半ラグランジュコードにおける摂動のない粒子運動の精度を評価し、前方半ラグランジュ（FSL）スキームが長期シミュレーションにおける質量およびエネルギー保存に顕著に寄与することを示している。Vlasov作用素を線形および非線形部分に分割し、両者にFSLを適用することで、特にトロイダル幾何における大きな位相空間的変位が引き起こす方向性の移流に起因する数値誤差を低減している。従来の後方半ラグランジュ（BSL）スキームでは、このような誤差が保存性を著しく損なうが、FSLによりその問題が緩和される。

ABSTRACT

Inaccurate description of the equilibrium can yield to spurious effects in gyrokinetic turbulence simulations. Also, the Vlasov solver and time integration schemes impact the conservation of physical quantities, especially in long-term simulations. Equilibrium and Vlasov solver have to be tuned in order to preserve constant states (equilibrium) and to provide good conservation property along time (mass to begin with). Several illustrative simple test cases are given to show typical spurious effects that one can observes for poor settings. We explain why Forward Semi-Lagrangian scheme bring us some benefits. Some toroidal and cylindrical GYSELA runs are shown that use FSL.

研究の動機と目的

• 不適切な平衡記述およびVlasovソルバー選択に起因するギャロキナeticシミュレーションにおける数値誤差の原因を調査すること。
• 時間積分法および移流スキームが長期シミュレーションにおける質量およびエネルギー保存に与える影響を特定すること。
• 線形／非線形分割および前方半ラグランジュ（FSL）スキームの使用が、数値安定性および保存性に与える影響を評価すること。
• 円柱およびトロイダルテストケースにおいて、FSLと後方半ラグランジュ（BSL）スキームの質量およびエネルギー保存性を比較すること。
• FSLが不適切なゼロナルフローの励起を低減し、自己無撞着場評価の精度を向上させることを実証すること。

提案手法

• 本研究では、磁場配置における粒子運動の移流項（E×Bおよび曲率移流を含む）を有するギャロキナetic Vlasov方程式を用いる。
• Vlasov方程式を線形（L）および非線形（N）作用素に分割し、大きな変位と非線形力学を分離する。
• 前方半ラグランジュ（FSL）スキームをLおよびN作用素の両方に適用し、粒子軌道追跡および場の評価を高精度に保証する。
• 時間積分にはストラング分割法を用い、方向ごとの移流を (v∥/2, φ/2, rθ, φ/2, v∥/2) の順序で独立にLおよびN作用素に適用する。
• μ=0 の4次元円柱およびトロイダル幾何でシミュレーションを実施し、同一の初期条件の下でFSLとBSLスキームを比較する。
• 質量およびエネルギー保存性を時間経過とともにモニタリングし、数値的安定性および発散の有無を評価するため、電場ポテンシャル場を可視化する。

実験結果

リサーチクエスチョン

• RQ1 ギャロキナetic平衡記述およびVlasovソルバー選択の不正確さが、長期シミュレーションにおける質量およびエネルギー保存性に与える影響は何か？
• RQ2 線形／非線形分解を施さない標準的なストラング分割法では、なぜ粒子移流に顕著な数値誤差が生じるのか？
• RQ3 前方半ラグランジュ（FSL）スキームが、後方半ラグランジュ（BSL）スキームと比較して保存性特性をどのように改善するのか？
• RQ4 Vlasov作用素の線形／非線形分割が、トロイダル幾何における方向性移流に起因する誤差をどのように低減するのか？
• RQ5 FSLは、自己無撞着シミュレーションにおける電場評価および場の進化の精度をどの程度向上させるのか？

主な発見

• FSLスキームでは、長期シミュレーションにおいて質量保存性が10−6の精度で維持される一方、BSLスキームでは特にトロイダル配置において顕著な質量劣化が観察された。
• トロイダル幾何において、BSLシミュレーションの主な質量損失源は、位相空間内での大きな変位を伴う線形作用素Lであった。
• FSLは不適切なゼロナルフローの振動および場の発散を低減し、電場ポテンシャル場が時間経過とともに安定を保つ一方、BSLではt=8000以降に場が発散する。
• 全エネルギーはFSLでより良好に保存され、理論的期待値（ゼロ）に近い状態を維持するが、BSLでは質量保存性の悪化に起因して偏差が増大する。
• 線形／非線形分割と、両作用素にFSLを適用する組み合わせにより、方向性移流誤差が排除され、各サブステップにおける場の評価が高精度に保たれる。
• 円柱およびトロイダル両テストケースにおいて、FSLはより正確で安定した電場ポテンシャル場を生成し、視覚的比較でBSL結果が長期にわたり顕著に発散していることが明確に示された。