[論文レビュー] Gaussian-Moment Relaxation Closures for Verifiable Numerical Simulation of Fast Magnetic Reconnection in Plasma

Evan Johnson|arXiv (Cornell University)|Sep 24, 2014

Gas Dynamics and Kinetic Theory参考文献 32被引用数 15

ひとこと要約

本稿では、速度分布関数のガウス近似に基づくモーメント緩和閉じ込め法を二流体磁気力学（MHD）に導入し、運動論的シミュレーションと一致する検証可能で収束する高速磁気リコネクションのシミュレーションを可能にする。圧力テンソルの等方化と熱フラックスに対するガウス-BGK閉じ込めを組み込むことで、特異性を回避し、拡散領域における主要な物理を捉える。数値結果は、ピークリコネクションにおいて良好な一致を示すが、明示的な衝突項がなければ不安定になる。

ABSTRACT

The motivating question for this dissertation was to identify the minimal requirements for fluid models of plasma to allow converged simulations that agree well with converged kinetic simulations of fast magnetic reconnection. We show that truncation closure for the deviatoric pressure or for the heat flux results in singularities. Due to the strong pressure anisotropies that arise, we study magnetic reconnection with a Gaussian-moment two-fluid MHD with isotropization of the pressure tensor. For the GEM magnetic reconnection challenge problem, our deviatoric pressure tensor agrees well with published kinetic simulations at the time of peak reconnection, but sometime thereafter the numerical solution becomes unpredictable and develops near-singularities that crash the simulation unless positivity limiters are applied. To explain these difficulties, we show that steady reconnection requires heat flux. Specifically, for two-dimensional problems invariant under 180-degree rotation about the X-point, entropy production in the vicinity of the X-point is necessary for reconnection and causes solutions to become singular for conservative models that lack heat flux. This prompts the need for a 10-moment gyrotropic heat flux closure. Using a Chapman-Enskog expansion with a Gaussian-BGK collision operator yields a heat flux closure for a magnetized 10-moment charged gas which generalizes the closure of McDonald and Groth.

研究の動機と目的

• 高速磁気リコネクションの収束性・検証可能性のあるシミュレーションを実現するための最小限の流体モデル要件を特定すること。
• 磁気リコネクションにおける強い圧力異方性によって生じる、切断された流体モデルにおける特異性の問題を解消すること。
• 断熱的五モーメントまたは十モーメントモデルが、明示的な衝突演算子なしでは収束性・特異性のない定常状態のリコネクションを支持できないことを確立すること。
• 拡散領域における主要な運動論的物理を捉える最小限で物理的に整合性のあるモデルを構築すること。
• GEMリコネクションチャレンジ問題に対するモデルの妥当性を検証し、数値的安定性および収束性を分析すること。

提案手法

• 特異性を防ぐために圧力テンソルの等方化を組み込んだガウスモーメント二流体MHDモデルを導入する。
• ガウス-BGK衝突演算子を適用し、マクドナルドとグロスの閉じ込めを一般化した十モーメント熱フラックス閉じ込めを導出する。
• 磁化プラズマにおけるBGK衝突モデルから、チャップマン＝エンスコーグ展開を用いて熱フラックス閉じ込めを導出する。
• GEM磁気リコネクションチャレンジのシミュレーションに、等方化を施した断熱的二流体Maxwell系を用いてモデルを実装する。
• 収束性を検討するために、調整可能な等方化時間τ_sを有する種内衝突項（R_sおよびQ_s）を組み込む。
• 熱フラックスを無視した場合のシミュレーション安定化のため、正値制限子を適用し、数値的不安定性を明らかにする。

実験結果

リサーチクエスチョン

• RQ1 収束性・特異性のない運動論的シミュレーションと一致する高速磁気リコネクションのための最小限の流体モデル構造は何か？
• RQ2 断熱的閉じ込めを用いた切断された流体モデル（例：五モーメントまたは十モーメント）が、収束性のある定常状態リコネクションを支持できない理由は何か？
• RQ3 熱フラックスは定常磁気リコネクションを維持するために果たす役割は何か？なぜ数値的安定性に不可欠なのか？
• RQ4 ガウス-BGKに基づく熱フラックス閉じ込めは、エントロピーを尊重する代替手法よりも物理的に整合性があり、数値的に安定な代替手段を提供できるか？
• RQ5 明示的な衝突項（等方化）の導入が、二流体MHDシミュレーションの収束性および精度に与える影響は何か？

主な発見

• 歪みテンソルまたは熱フラックステンソルの切断閉じ込めは、定常状態の磁気リコネクションシミュレーションで特異性を引き起こす。
• 断熱的五モーメントまたは十モーメントモデルは、明示的な衝突演算子なしでは収束性・特異性のない定常状態リコネクションを支持できない。
• 圧力の等方化とガウス-BGK熱フラックス閉じ込めを組み込んだ提案されたガウスモーメント二流体MHDモデルは、ピークリコネクションにおいて運動論的シミュレーションと良好に一致する歪み圧力を生成する。
• 明示的な熱フラックス項がないシミュレーションは不安定になり、近似的に特異的になるため、クラッシュを防ぐために正値制限子が必要となる。
• リコネクションにおける持続的熱フラックスの必要性は、動的リコネクションが定常駆動過程に近いことを示唆し、十モーメントのギロトロープ的閉じ込めを必要とする。
• 導出されたガウス-BGK熱フラックス閉じ込めは、マクドナルドとグロスの閉じ込めを一般化しており、この文脈においてエントロピーを尊重する代替手法よりも物理的に整合性が高いかが示された。