[論文レビュー] On the Impact of the Numerical Method on Magnetic Reconnection and Particle Acceleration -- I. The MHD case

Eleonora Puzzoni, A. Mignone+1 more|arXiv (Cornell University)|Sep 27, 2021

Solar and Space Plasma Dynamics参考文献 85被引用数 14

ひとこと要約

本研究は、2次元MHDシミュレーションにおける剥離不安定な電流シートの磁気再結合と粒子加速に、数値解法が与える影響を調査する。PLUTOコードを用い、リーマン解法、再構成法、グリッド解像度を変化させた結果、再結合率は十分な解像度および有限のランジエール数の下でのみ収束し、粒子加速は強固で、指数が ≈1.7 のパワーロー spectrum を示し、数値的詳細に依存しない。これは非線形・高速再結合領域に入った後である。

ABSTRACT

We present 2D MHD numerical simulations of tearing-unstable current sheets coupled to a population of non-thermal test-particles, in order to address the problem of numerical convergence with respect to grid resolution, numerical method and physical resistivity. Numerical simulations are performed with the PLUTO code for astrophysical fluid dynamics through different combinations of Riemann solvers, reconstruction methods, grid resolutions at various Lundquist numbers. The constrained transport method is employed to control the divergence-free condition of magnetic field. Our results indicate that the reconnection rate of the background tearing-unstable plasma converges only for finite values of the Lundquist number and for sufficiently large grid resolutions. In general, it is found that (for a 2nd-order scheme) the minimum threshold for numerical convergence during the linear phases requires the number of computational zones covering the initial current sheet width to scale roughly as $\sim \sqrt{\bar{S}}$, where $\bar{S}$ is the Lundquist number defined on the current sheet width. On the other hand, the process of particle acceleration is found to be nearly independent of the underlying numerical details inasmuch as the system becomes tearing-unstable and enters in its nonlinear stages. In the limit of large $\bar{S}$, the ensuing power-law index quickly converge to $p \approx 1.7$, consistently with the fast reconnection regime.

研究の動機と目的

• グリッド解像度、数値スキーム、物理的抵抗率に対する磁気再結合率の収束を評価すること。
• 異なる数値手法が、剥離不安定な電流シートのダイナミクスおよびそれに続く粒子加速に与える影響を調査すること。
• 再結合の線形段階および非線形段階における数値収束に必要な最小グリッド解像度を特定すること。
• さまざまな数値手法および抵抗率レベルにおいて、粒子エネルギースペクトルの頑健性を評価すること。
• 特に高速再結合領域において、粒子加速が数値的詳細に依存しなくなる条件を確立すること。

提案手法

• 発散自由な磁場を維持するための制約付き伝搬法を用いた2次元MHD電流シートの数値シミュレーション。
• リーマン解法（HLL、HLLD）、再構成法（線形、5次精度WENO-Z）、および異なるランジエール数におけるグリッド解像度を体系的に変化。
• 高解像度シミュレーションにおける精度と安定性を向上させるために、UCT-HLLD電磁場平均化スキームを採用。
• 時間変化するMHD背景場における粒子エネルギー増幅を追跡するために、テスト粒子アプローチを採用。
• 複数のシミュレーション設定において、再結合率、線形成長段階、粒子エネルギースペクトルを分析。
• 収束性および数値拡散効果を評価するために、異なる数値スキーム間の結果を比較。

実験結果

リサーチクエスチョン

• RQ1 剥離不安定な電流シートにおける再結合率の数値収束に必要な最小グリッド解像度は何か？
• RQ2 異なるリーマン解法および再構成スキームは、剥離不安定性の発生および進化にどのように影響するか？
• RQ3 粒子エネルギースペクトルは、数値手法および抵抗率の選択にどの程度依存するか？
• RQ4 非線形・高速再結合段階において、粒子エネルギー分布のパワーロー指数が普遍的値に収束するか？
• RQ5 数値的抵抗率のみが、理想MHDシミュレーションにおいてプラズモイドの形成と粒子加速を引き起こすか？

主な発見

• 再結合率は、有限のランジエール数および十分なグリッド解像度がある場合にのみ収束し、収束には初期電流シート幅に対して約 √S̄ 個の計算ゾーンが必要である。
• 2次精度スキームでは、線形段階の収束が、電流シート幅あたりのゾーン数が ∼√S̄/10⁴ に比例するとき達成される。S̄ は電流シート幅上で定義されるランジエール数である。
• HLLD リーマン解法と5次精度WENO-Z再構成、およびUCT-HLLD EMF平均化を組み合わせた手法は、S̄ = 10⁴ において a/Δx ≈ 10 で収束を達成し、精度が低いスキームを上回る性能を示す。
• 非線形段階において、空間的に平均化された横磁場はランジエール数に依存せず同じ値に達する。これは乱流ダイナミクスによって支配される普遍的挙動を示している。
• S̄ ≳ 10⁴ の場合、粒子エネルギースペクトルは速やかに指数 p ≈ 1.7 のパワーローに収束し、高速再結合領域と整合的である。
• システムが非線形・高速再結合領域に入った後、粒子加速は数値的詳細および抵抗率レベルにほとんど依存せず、最もエネルギーの高い粒子はプラズモイドコアの周囲の磁化リングに集中している。