Failed ejection and oscillations of a current-carrying filament balanced by gravity

问题描述
太阳丝状结构（日冕物质抛射）是悬浮在太阳日冕中的高密度、低温等离子体结构，通常与太阳爆发和日冕物质抛射（CME）相关。尽管许多模型描述了丝状结构的稳定性与爆发，但现有数值模拟的一个显著局限在于假设等离子体贝塔值（plasma beta）为零，从而忽略了重力的影响。真实的日冕环境涉及有限的等离子体压力和重力分层，这对于理解丝状结构的稳定性、质量负载以及热结构至关重要。此外，关于丝状结构爆发的观测数据经常显示出“失败”的事件，即抛射过程被中止，随后丝状结构回落，这一现象需要能够充分捕捉底层物理机制的模型来考虑重力效应和真实的电流配置。

研究方法
作者利用 Solov'ev (2010) 提出的重力平衡模型研究了丝状结构在失稳后的演化过程。该研究采用了双重方法：

1. 解析研究：作者基于 Solov'ev 模型，推导出了电流密度、等离子体密度及其空间分布的解析表达式。在该模型中，携带电流的高密度丝状结构通过其亚光球镜像电流与自身电流之间的排斥洛伦兹力来抵消重力。
2. 数值模拟：使用 Lare2d 数值代码构建了一个二维磁流体力学（MHD）模型。该模型在笛卡尔坐标系下求解电阻性、可压缩 MHD 方程，并明确包含了重力和电阻率。模拟区域（$x \in [-10, 10]$ Mm，$y \in [0, 40]$ Mm）由 $2000 \times 4000$ 的网格划分。初始平衡态的特征为磁场 $B_0 = 10^{-3}$ T，质量密度 $\rho_0 \sim 1.3 \times 10^{-9}$ kg m$^{-3}$，温度 $T \sim 13000$ K。研究采用了一种修正的温度剖面来模拟光球层与日冕层之间的过渡区。

系统通过两种不同的情景进行失稳：(a) 增加丝状结构中的电流（通过增加磁场强度 $B_0$），以及 (b) 降低丝状结构的质量密度。通过追踪系统的演化过程，以观察抛射、停止和回落的动力学过程。

主要贡献

• 解析表达式的推导：论文提供了电流密度、等离子体密度及其作为模型参数（包括总电流和丝状结构重量）函数的空间分布的解析公式。
• 重力平衡 MHD 模型：本研究实现了一个 2D MHD 模型，该模型明确平衡了排斥性的磁力与重力，解决了忽略重力效应的零贝塔模型的局限性。
• 失败抛射的特征化：作者系统地模拟并表征了“失败抛射”现象，即失稳的丝状结构被抛射上升，在特定高度停止，随后回落并进入振荡循环。
• 定量依赖关系：最大抛射高度被确定为关于失稳电流量级和丝状结构等离子体稀释程度的函数。

结果

• 振荡行为：在两种失稳机制（电流增加和密度降低）下，丝状结构都会向上抛射，在一定高度停止并回落。该循环以约 600 秒的周期重复。
• 抛射速度：最大抛射速度根据失稳机制的不同而不同：电流增加情景下约为 80 km s$^{-1}$，密度降低情景下约为 40 km s$^{-1}$。
• 电流片形成：在被抛射的丝状结构下方，形成了一个发生磁重联的电流片。磁重联破坏了连接上方和下方电流的磁力线。
• 参数依赖性：
  • 增加电流（磁场强度）会导致更高的抛射高度。
  • 降低丝状体密度也会导致更高的抛射高度，较低的密度比（$\rho/\rho_0$）会导致更高的上升高度。
  • 与电流增加情景相比，密度降低情景下的丝状结构能更快达到最大高度并停留更长时间（约 200 秒），尽管振荡周期相似。
• 不对称性：在电流增加的情况下，丝状结构下方的密度结构最终会产生不对称性，导致丝状结构偏离垂直轴。在密度降低的情况下，丝状结构沿垂直轴振荡直至达到稳定状态。
• 验证：通过在扩展垂直域中的对比模拟确认，原始数值域的大小不会影响磁场配置或达到的最大高度。

意义
本文证明了将重力和有限等离子体贝塔值纳入模型，对于理解爆发前的平衡状态以及丝状结构失去稳定性的具体条件是必不可少的。通过重现观测到的“失败爆发”现象及随后的振荡行为，本研究为解释有关丝状结构形成与失稳的观测数据提供了理论框架。研究结果强调，质量负载和热结构是影响不稳定性阈值及爆发过程非线性动力学的关键因素，为理想真空模型提供了一个更具现实意义的替代方案。