[論文レビュー] A Unified Gas Kinetic Scheme for Multi-scale Plasma Transport

Chang Liu, Kun Xu|arXiv (Cornell University)|Sep 17, 2016

Ionosphere and magnetosphere dynamics被引用数 3

ひとこと要約

本稿では、電子およびイオンのVlasov-BGK方程式とマクスウェル方程式を統合することで、多スケールプラズマ輸送のための統一ガス動力学的スキーム（UGKS）を提案する。この手法により、運動論的から流体的プラズマ状態への滑らかな遷移が可能となり、粒子輸送、衝突、電磁場の結合を分割しない時間発展型でスケール依存のフラックス再構築によって、衝突なしから高衝突性プラズマ動態に至る包括的取り扱いが実現される。

ABSTRACT

A unified gas kinetic scheme (UGKS) for multi-scale and multi-component plasma transport is constructed. The current scheme is a direct modeling method, where the time evolution solutions from the Vlasov-BGK equations for both electron and ion, and the Maxwell equations are used to construct the scale-dependent plasma simulation. As a result, based on the modeling scales of mesh size and time step, the discretized governing equations for the whole plasma regimes are obtained. The UGKS takes into account the electron inertia, full electromagnetic field equations, and separate electron and ion evolution. The physics recovered in UGKS ranges from the kinetic Vlasov equation to the hydrodynamic magnetohydrodynamic (MHD) equations, with a unified treatment in all scales from the collisionless particle transport to the hydrodynamic wave interactions. The UGKS presents a plasma description which is more general than the Vlasov equation in the kinetic scale and all kinds of MHD equations in the hydrodynamic scale, such as Hall, Resistive, and Ideal Magnetohydrodynamics (MHD). All above single scale equations become the subsets of the UGKS. The key dynamics in UGKS is the non-splitting treatment of particle collision, acceleration, and transport in the construction of the numerical flux at cell interface, and this flux is evaluated from the scale and time evolving integral solution of the Vlasov-BGK model. The evolution of electromagnetic field and the plasma flow are coupled in a semi-implicit manner. The UGKS is able to give a physically reliable solution for plasma evolution from the collisionless limit to the highly collisional one.

研究の動機と目的

• 運動論的から流体的状態まで、あらゆるスケールのプラズマ挙動を捉える単一の計算フレームワークの開発を目的とする。
• Vlasov法で運動論的、MHD法で流体的挙動を扱う必要がある既存手法の限界を克服することを目的とする。
• 慣性および電磁場の時間発展を含む電子およびイオンの動態を、一貫した数値スキームで統合的に取り扱うことを目的とする。
• 別々のソルバーを必要としないように、同一のアルゴリズムフレームワーク内で理想MHD、ホールMHD、抵抗性MHDなどの既知の式が特殊ケースとして還元されることを目的とする。
• 粒子輸送、衝突、電磁場発展を分割しないフラックス定式化によって、スケールをまたいで物理的一致性を確保することを目的とする。

提案手法

• 本スキームは、電子およびイオンのVlasov-BGK方程式の直接的モデリングに基づき、時間発展解を用いてセル界面における数値フラックスを構築する。
• フラックスは、スケールおよび時間に依存するVlasov-BGKモデルの積分解から計算され、粒子輸送、衝突、加速を包括的にかつ正確に取り扱う。
• マクスウェル方程式は、フィールド-粒子の一貫性を保つために、半陰的結合戦略でプラズマ方程式と同時に解く。
• メッシュサイズおよび時間ステップの調整により、スケールに応じた解の有効スケールを決定することで、衝突なしから衝突性への自然な遷移を内蔵的に捉える。
• 粒子衝突、加速、輸送をフラックス評価において同時に取り扱うことで、オペレータ分割を回避し、スケールをまたいで物理的正確性を維持する。
• 系が流体的極限に達した場合に、理想MHD、抵抗性MHD、ホールMHDなどの既知の式に自然に還元されることで、既存のモデルと一貫性を保つ。

実験結果

リサーチクエスチョン

• RQ1 モデル切り替えを要せず、運動論的から流体的状態まで、あらゆるスケールのプラズマ輸送を正確にシミュレートできる単一の数値スキームが可能か？
• RQ2 複数スケールにわたり、粒子衝突、電磁場、輸送を分割しない方法で一貫して結合できるか？
• RQ3 UGKSが、理想MHD、ホールMHD、抵抗性MHDなどの既知のプラズマモデルを、極限ケースとしてどれほど正確に回復できるか？
• RQ4 同一のアルゴリズムフレームワークで、衝突なしと高衝突性プラズマ状態の両方において物理的正確性を維持できるか？
• RQ5 Vlasov-BGKモデルの時間発展型積分解が、スケール適応型プラズマシミュレーションを可能にする役割は何か？

主な発見

• UGKSは、運動論的状態におけるVlasov方程式と、流体的状態におけるすべての標準的MHDモデル（理想MHD、抵抗性MHD、ホールMHD）を包含する統一的プラズマ動態記述を提供する。
• 局所的なメッシュおよび時間ステップスケールに応じて適応することで、衝突なしから衝突性プラズマ挙動への遷移を的確に捉える。
• 粒子輸送、衝突、電磁場発展の分割しない取り扱いにより、あらゆるスケールで物理的一致性が保証される。
• 運動論的極限および流体的極限の両方において、正確性と安定性を維持し、別々のソルバーやモデル切り替えの必要がない。
• 適切な漸近的極限において、理想MHDや抵抗性MHDなどの既知の式に自然に還元されることで、既存のプラズマモデルと一貫性が裏付けられる。
• 電磁場とプラズマ流れの半陰的結合により、広範なプラズマ状態において安定かつ物理的に信頼できるシミュレーションが可能になる。