[论文解读] 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

一句话总结

本研究探讨了数值方法对二维MHD模拟中撕裂不稳定性电流片磁重联与粒子加速的影响。利用PLUTO代码，通过改变Riemann求解器、重构方法和网格分辨率，发现重联率仅在足够高分辨率和有限Lundquist数条件下实现收敛；而粒子加速则具有鲁棒性，产生幂律能谱，指数≈1.7，一旦系统进入非线性快速重联阶段，该结果与数值细节无关。

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.

研究动机与目标

• 评估磁重联率在网格分辨率、数值格式和物理电阻率方面的收敛性。
• 研究不同数值方法对撕裂不稳定性电流片动力学及相应粒子加速的影响。
• 确定线性与非线性重联阶段实现数值收敛所需的最小网格分辨率。
• 评估不同数值方法与电阻率水平下粒子能量谱的鲁棒性。
• 确定粒子加速在何种条件下可独立于数值细节，特别是在快速重联阶段。

提出的方法

• 使用PLUTO代码对二维MHD电流片进行数值模拟，采用约束传输法以保持磁场无散度。
• 系统性地改变Riemann求解器（HLL、HLLD）、重构方法（线性、五阶WENO-Z）和不同Lundquist数下的网格分辨率。
• 采用UCT-HLLD电磁场平均方案以提升高分辨率模拟中的精度与稳定性。
• 采用测试粒子方法，追踪在演化MHD背景场中的粒子能量增益。
• 分析多种模拟配置下的重联率、线性增长阶段及粒子能量谱。
• 对比不同数值格式的结果，以评估收敛性与数值扩散效应。

实验结果

研究问题

• RQ1 撕裂不稳定性电流片中，重联率实现数值收敛所需的最小网格分辨率是多少？
• RQ2 不同Riemann求解器与重构方法如何影响撕裂不稳定的触发与演化？
• RQ3 粒子能量谱在多大程度上依赖于数值方法与电阻率的选择？
• RQ4 在非线性快速重联阶段，粒子能量分布的幂律指数是否收敛至一个普适值？
• RQ5 在理想MHD模拟中，仅靠数值电阻率是否足以触发磁岛形成与粒子加速？

主要发现

• 重联率仅在有限Lundquist数和足够高分辨率下实现收敛，收敛要求初始电流片宽度范围内约√S̄个计算单元。
• 对于二阶格式，线性阶段收敛要求每个电流片宽度的单元数与∼√S̄/10⁴成比例，其中S̄为基于电流片宽度定义的Lundquist数。
• 采用HLLD Riemann求解器、五阶WENO-Z重构与UCT-HLLD电磁场平均方案时，在a/Δx ≈ 10（S̄ = 10⁴）条件下实现收敛，优于精度较低的方案。
• 在非线性阶段，平均横向磁场达到与Lundquist数无关的相同值，表明其行为具有普适性，受湍流动力学主导。
• 当S̄ ≳ 10⁴时，粒子能量谱迅速收敛为幂律谱，幂律指数p ≈ 1.7，与快速重联阶段一致。
• 一旦系统进入非线性快速重联阶段，粒子加速几乎与数值方法和电阻率水平无关，能量最高的粒子集中于磁岛核心周围的磁化环状区域。