[論文レビュー] Mixed semi-Lagrangian/finite difference methods for plasma simulations

Francis Filbet, Chang Yang|arXiv (Cornell University)|Sep 30, 2014

Magnetic confinement fusion research参考文献 13被引用数 4

ひとこと要約

本稿では、任意の領域を有する直交メッシュ上で4次元ドリフト運動論的モデルを用いた長期間のプラズマシミュレーションのため、混合半ラグランジュ/有限差分法を提案する。線形期には効率性を発揮する半ラグランジュ法を、非線形期には安定性および保存性を確保する有限差分法を組み合わせることで、複雑な構造を有するイオン乱流のシミュレーションにおいて、高い精度、高いロバスト性、優れた保存性を実現する。

ABSTRACT

In this paper, we present an efficient algorithm for the long time behavior of plasma simulations. We will focus on 4D drift-kinetic model, where the plasma's motion occurs in the plane perpendicular to the magnetic field and can be governed by the 2D guiding-center model. Hermite WENO reconstructions, already proposed in \cite{YF15}, are applied for solving the Vlasov equation. Here we consider an arbitrary computational domain with an appropriate numerical method for the treatment of boundary conditions. Then we apply this algorithm for plasma turbulence simulations. We first solve the 2D guiding-center model in a D-shape domain and investigate the numerical stability of the steady state. Then, the 4D drift-kinetic model is studied with a mixed method, i.e. the semi-Lagrangian method in linear phase and finite difference method during the nonlinear phase. Numerical results show that the mixed method is efficient and accurate in linear phase and it is much stable during the nonlinear phase. Moreover, in practice it has better conservation properties.

研究の動機と目的

• 4次元ドリフト運動論的モデルに従うプラズマ乱流の長期間シミュレーションのための安定的かつ高精度な数値法の開発。
• 特にD字型などの非長方形幾何形状を含む任意の計算領域における境界条件の実装という課題に対処する。
• 半ラグランジュ法の欠点（質量保存の欠如）と有限差分法の制限（CFL制約）を補い、両者の長所を統合する。
• 特に非線形期における小スケール構造を伴う長期的シミュレーションにおいて、質量、エネルギー、エントロピーといった保存量が保持されることを保証する。
• 2次元ガイドセンターおよび4次元ドリフト運動論的モデルの両方で手法を検証し、現実のプラズマ乱流シナリオにおけるロバスト性を示す。

提案手法

• 半ラグランジュおよび有限差分フレームワークの両方において、高次元で振動を抑えるVlasov方程式の解法として、Hermite WENO再構成を適用する。
• 4次元ドリフト運動論的モデルおよび2次元ガイドセンター・モデルの離散化に直交メッシュを用い、複雑な幾何形状の柔軟な取り扱いを可能にする。
• Vlasov方程式に対して、D字型領域を含む任意の領域における境界条件を実装するため、逆Lax-Wendroff法を導入する。
• 電場ポテンシャルを計算するためのポアソン方程式の解法に補間技術を用い、境界処理と整合性を保つ。
• 小スケール構造が出現し、数値的安定性が重要になる非線形期の開始時に、半ラグランジュ法から有限差分法に切り替える。
• メッシュの細分化に伴う収束性を評価するため、精密なメッシュ（$64\times64\times32\times65$）からの参照解を用いる。

実験結果

リサーチクエスチョン

• RQ1 半ラグランジュ法と有限差分法を組み合わせた混合数値法は、複雑な幾何形状を有する長期間のプラズマシミュレーションにおいて、高精度および安定性を維持できるか？
• RQ2 プラズマ乱流の線形期から非線形期への遷移に伴い、混合法は保存量（質量、エネルギー、エントロピー）をどの程度正確に保持できるか？
• RQ3 Hermite WENO再構成の使用は、分布関数に小スケール構造が存在する状況で、数値的安定性および精度をどの程度向上させるか？
• RQ4 逆Lax-Wendroff法は、D字型などの非長方形領域における境界条件処理に、どの程度効果的か？
• RQ5 純粋な半ラグランジュ法や有限差分法と比較して、混合法は非線形期における保存性および微細構造の解像度において優れているか？

主な発見

• 線形期には解が滑らかであるため、半ラグランジュ法の効率性が最大限に発揮され、混合法は高い精度と安定性を達成する。
• 非線形期においては、有限差分部が不適切な振動を効果的に抑制し、保存性を維持する。これにより、長期的シミュレーションにおいて半ラグランジュ法を上回る性能を示す。
• $x$および$y$方向のメッシュ細分化は$L^2$ノルムおよびエントロピー保存性を顕著に向上させるが、$v$方向の細分化はわずかな改善にとどまる。
• エネルギー保存性は$x$および$y$方向の細分化によって最も向上し、配置空間における空間分解能が物理的不変量を維持するために極めて重要であることが確認された。
• 分布関数は線形期に指数関数的増大を示し、$t=3000$では5つの回転する渦とフィラメント構造に飽和する。$t=4000$ではさらに複雑な小スケール構造が形成される。
• $t=8000$では、系は相対的に安定状態に達し、微細構造が解像される。混合法は、純粋な半ラグランジュ法が失敗する状況でも保存性を的確に保持した。