[論文レビュー] Enhanced conservation properties of Vlasov codes through coupling with conservative fluid models

Thomas Alexander Trost, S. Lautenbach+1 more|arXiv (Cornell University)|Feb 1, 2017

Gas Dynamics and Kinetic Theory参考文献 21被引用数 3

ひとこと要約

本稿では、運動論的コードをモーメントベースのフィードバックによって補正する保守的流体モデルを用いて、Vlasovシミュレーションにおけるエネルギー保存を向上させる新規結合手法を提案する。一般化された10モーメントフィッティング手順を適用することで、数値的加熱が顕著に低減されつつ、運動論的特徴が維持される。10モーメントフィッティングは安定的かつ有効であるが、等方的または5モーメントアプローチとは異なり、実用的でない。

ABSTRACT

Many phenomena in collisionless plasma physics require a kinetic description. The evolution of the phase space density can be modeled by means of the Vlasov equation, which has to be solved numerically in most of the relevant cases. One of the problems that often arise in such simulations is the violation of important physical conservation laws. Numerical diffusion in phase space translates into unphysical heating, which can increase the overall energy significantly, depending on the time scale and the plasma regime. In this paper, a general and straightforward way of improving conservation properties of Vlasov schemes is presented that can potentially be applied to a variety of different codes. The basic idea is to use fluid models with good conservation properties for correcting kinetic models. The higher moments that are missing in the fluid models are provided by the kinetic codes, so that both kinetic and fluid codes compensate the weaknesses of each other in a closed feedback loop.

研究の動機と目的

• 位相空間離散化と人工的加熱によって引き起こされるVlasovシミュレーションにおける数値的エネルギー非保存性を是正すること。
• 運動論的プラズマコードにおける保存性特性を向上させるための一般化された、スキームに依存しない手法を開発すること。
• 流体モデルの保存性の強みと運動論的コードの高次モーメント忠実性を、相互フィードバックによって活用すること。
• 線形および非線形プラズマ領域におけるモーメントベース補正の有効性を検証すること。
• 保守的モーメントフィッティングを用いて、流体・運動論的コード間の強固な結合を、マルチリージョンシミュレーションで実現すること。

提案手法

• Vlasovソルバと保守的流体モデルを閉ループフィードバックで結合し、流体モーメントが運動論的位相空間分布を補正する。
• 速度空間における線形変換を用いて、運動論的分布を所定の流体モーメントにフィットさせる。これは、従来の等方的フィッティングを非等方的ケースに一般化したものである。
• フィッティング行列は勾配降下法による正確な計算または近似計算のいずれかで算出され、両手法の安定性と精度が検証された。
• 空間微分にはCWENO法、時間積分には時間分割法を用いた有限差分Vlasovソルバに実装された。
• 10モーメントフィッティングは密度、平均速度、エネルギー密度テンソルを一致させ、5モーメントフィッティングは比較のためのベースラインとして使用された。
• GEMチャレンジベンチマークを用いて、再結合率、エネルギー変化、分布関数の形状を比較することでシミュレーションの妥当性を検証した。

実験結果

リサーチクエスチョン

• RQ1 保守的流体モデルは、不自然な効果を引き起こさずに、Vlasovシミュレーションにおける数値的加熱を効果的に補正できるか？
• RQ2 非等方的10モーメントフィッティングは、等方的または5モーメントフィッティングと比較して、運動論的特徴と保存則の維持にどの程度優れているか？
• RQ3 フィッティング行列の正確な計算と近似計算の違いが、エネルギー保存性とシミュレーション安定性に与える影響は何か？
• RQ4 この結合手法は、マルチリージョンプラズマシミュレーションにおける流体-運動論的界面問題に一般化可能か？
• RQ5 フィードバックループは、非線形領域においても人工的加熱を顕著に低減し、物理的正確性を維持できるか？

主な発見

• 近似行列計算を用いた10モーメントフィッティングは、補正なしのVlasovシミュレーションと比較して、全エネルギー増加を最大50％まで低減し、エネルギー保存性が顕著に向上した。
• 等方的フィッティング手順は数値的不安定性を引き起こし、一般手法の特殊ケースではあるが、実用的ではないことが判明した。
• 5モーメントフィッティングは、重要な運動論的特徴を破壊し、不自然な不安定性を誘発するため、シミュレーションは物理的に意味を持たなくなる。
• GEMチャレンジのシミュレーションでは、補正済みVlasovコードが、高時間分解能でも参照解に近い再結合フラックスおよび電流プロファイルを維持した。
• 本手法は、さまざまな解像度および時間刻みに対して安定しており、CWENO法における正則化パラメータεの選択にほとんど感度を示さなかった。
• フィードバックループは、位相空間の拡散を補正することでシミュレーションを安定化させつつ、プラズマの運動論的記述を損なわずに保った。