[论文解读] High Lundquist Number Resistive MHD Simulations of Magnetic Reconnection: Searching for Secondary Island Formation

C. S. Ng, Srivathsan Ragunathan|arXiv (Cornell University)|Jun 2, 2011

Ionosphere and magnetosphere dynamics被引用 3

一句话总结

本研究对 Lundquist 数高达 $2 \times 10^5$ 的磁重联过程进行了高分辨率阻抗 MHD 模拟，以探究二次岛（等离子体团）的形成。在无外部噪声的情况下，重联遵循 Sweet-Parker 标度 ($S^{-1/2}$)，且未形成等离子体团；然而，当加入小尺度噪声后，等离子体团出现，且重联速率不再依赖于 $S$，表明噪声诱导的湍流重联机制可在阻抗 MHD 中实现快速重联。

ABSTRACT

Recently, secondary island formation due to the tearing instability of the Sweet-Parker current sheet was identified as a possible mechanism that can lead to fast reconnection (less sensitive dependence on Lundquist number $S$) both in numerical simulations using Particle-in-Cell (PIC) method [Daughton et al. 2009], as well as using resistive magnetohydrodynamics (MHD) [Lapenta 2008; Bhattacharjee et al. 2009]. This instability is thought to appear when $S$ is greater than a certain threshold. These recent results prompt us to perform more resistive MHD simulations of a basic reconnection configuration based on the island coalescence instability, using much higher resolutions and larger $S$. Our simulations are based on a fairly standard pseudo spectral code, which has been tested for accuracy, convergence, and compared well with codes using other methods [Ng et al. 2008]. In our simulations, formation of plasmoids were not found, except when insufficient resolution was used, or when a small amount of noise was added externally. The reconnection rate is found to follow the Sweet-Parker scaling when no noise is added, but increases to a level independent of $S$ with noise, when plasmoids form. Latest results with $S$ up to $2 imes 10^5$ will be presented.

研究动机与目标

• 探究在高 Lundquist 数阻抗 MHD 模拟中，二次岛（等离子体团）的形成是否能自然发生。
• 确定先前研究中观察到的快速重联速率是源于数值伪影还是物理解释。
• 评估数值分辨率和外部噪声在触发等离子体团形成和增强重联速率中的作用。
• 比较无噪声与注入噪声的模拟中的重联动力学，以隔离快速重联出现的条件。
• 在受控条件下评估 Sweet-Parker 标度在高分辨率阻抗 MHD 中的有效性。

提出的方法

• 模拟使用伪谱代码求解归一化的二维不可压缩阻抗 MHD 方程，形式为 $\partial_t \Omega + [\phi, \Omega] = [A, J] + \nu \nabla_\perp^2 \Omega$ 和 $\partial_t A + [\phi, A] = \eta \nabla_\perp^2 A$。
• 通过电阻率 $\eta$ 控制 Lundquist 数 $S$，通过 $\eta \sim S^{-1}$ 实现高达 $2 \times 10^5$ 的 $S$ 值。
• 数值分辨率达到 $8192^2$ 个网格点，以确保收敛性并最小化因分辨率不足导致的虚假等离子体团形成。
• 在谱空间中添加外部随机噪声，控制其振幅（能量比约 ~$10^{-6}$）和关联时间（约 0.25），以模拟湍流扰动。
• 通过 $dA_0/dt$ 量化重联速率，其中 $A_0$ 为磁通函数最大幅值，并在重联阶段进行平均。
• 通过在不同 $S$ 值和分辨率下对比有噪声与无噪声的模拟，隔离噪声对等离子体团形成和重联速率的影响。

实验结果

研究问题

• RQ1 在高 Lundquist 数（$S \sim 10^5$）且充分解析的阻抗 MHD 模拟中，Sweet-Parker 电流片是否会对撕裂模不稳定性产生响应？
• RQ2 在无外部扰动的高分辨率阻抗 MHD 模拟中，二次岛（等离子体团）的形成是否能自然发生？
• RQ3 数值噪声或外部扰动在触发等离子体团形成和实现快速重联中起什么作用？
• RQ4 在无噪声和有噪声条件下，重联速率如何随 $S$ 变化？
• RQ5 模拟结果在多大程度上依赖于分辨率以及小振幅噪声的引入？

主要发现

• 在无外部噪声的高分辨率模拟中，重联速率遵循 Sweet-Parker 标度 $S^{-1/2}$，即使在 $S = 2 \times 10^5$ 时也未观察到等离子体团形成。
• 仅当向系统中加入小振幅外部噪声时，才会观察到等离子体团的形成与喷射，表明其存在类似阈值的不稳定性触发机制。
• 加入噪声后，重联速率近似与 $S$ 无关，且在大 $S$ 条件下显著高于 Sweet-Parker 预测值。
• 噪声的引入导致重联速率波动但整体增强，且在 $S$ 值从 $1.25 \times 10^4$ 到 $2.5 \times 10^4$ 的范围内，平均 $dA_0/dt$ 保持恒定。
• 结果表明，阻抗 MHD 中的等离子体团形成并非高 $S$ 条件下 Sweet-Parker 电流片的内在特征，而是外部扰动或数值噪声的结果。
• 本研究强调了分辨率以及小尺度扰动的存在在决定模拟中是否出现快速重联机制中的关键作用。