[论文解读] Magnetohydrodynamics dynamical relaxation of coronal magnetic fields. I. Parallel untwisted magnetic fields in 2D

Jorge Fuentes-Fernández, C. E. Parnell+1 more|St Andrews Research Repository (St Andrews Research Repository)|Oct 24, 2011

Solar and Space Plasma Dynamics参考文献 23被引用 10

一句话总结

本研究探讨了具有有限β值的二维日冕等离子体在非电阻磁流体动力学（MHD）下的弛豫过程，表明等离子体压强梯度显著改变了与传统无力场模型相比的平衡态。通过粘性阻尼磁声波，系统弛豫至一种非无力场磁流体静力平衡态，其中磁力与等离子体力达到平衡，即使在远离线性区的区域，解析预测结果也与数值结果高度一致。

ABSTRACT

Context. For the last thirty years, most of the studies on the relaxation of stressed magnetic fields in the solar environment have onlyconsidered the Lorentz force, neglecting plasma contributions, and therefore, limiting every equilibrium to that of a force-free field. Aims. Here we begin a study of the non-resistive evolution of finite beta plasmas and their relaxation to magnetohydrostatic states, where magnetic forces are balanced by plasma-pressure gradients, by using a simple 2D scenario involving a hydromagnetic disturbance to a uniform magnetic field. The final equilibrium state is predicted as a function of the initial disturbances, with aims to demonstrate what happens to the plasma during the relaxation process and to see what effects it has on the final equilibrium state. Methods. A set of numerical experiments are run using a full MHD code, with the relaxation driven by magnetoacoustic waves damped by viscous effects. The numerical results are compared with analytical calculations made within the linear regime, in which the whole process must remain adiabatic. Particular attention is paid to the thermodynamic behaviour of the plasma during the relaxation. Results. The analytical predictions for the final non force-free equilibrium depend only on the initial perturbations and the total pressure of the system. It is found that these predictions hold surprisingly well even for amplitudes of the perturbation far outside the linear regime. Conclusions. Including the effects of a finite plasma beta in relaxation experiments leads to significant differences from the force-free case.

研究动机与目标

• 探究有限等离子体β值在日冕磁场MHD弛豫过程中的作用，突破标准无力场近似。
• 确定等离子体压强梯度如何影响弛豫过程中最终的平衡态，特别是当磁力由等离子体压强与磁张力共同平衡时的影响。
• 评估基于线性解析预测在粘性阻尼MHD系统中非线性弛豫过程中的有效性。
• 考察等离子体在弛豫过程中的热力学演化，特别是能量从等离子体向磁场的转移过程。
• 确立最终平衡态仅取决于初始扰动和总压强，而与粘性系数或初始扰动的性质（密度、温度或内能变化）无关。

提出的方法

• 采用包含粘性耗散的全二维MHD数值模拟代码，驱动受扰动的磁场与等离子体状态的弛豫过程。
• 初始扰动通过局部改变等离子体压强、密度或内能引入，同时保持磁力线冻结条件。
• 磁声波（快波与慢波）的粘性阻尼作为唯一的耗散机制，使系统弛豫至最终的MHD平衡态。
• 基于绝热假设及沿磁力线与垂直于磁力线的波传播，推导出线性解析解，用于预测最终的压强、密度与磁场分布。
• 将基于总压强守恒与绝热演化的解析预测结果，与不同初始扰动幅度和等离子体β值下的数值结果进行对比。
• 采用二维高斯压强扰动测试非线性区行为，并量化与绝热预测的偏差，尤其关注等离子体密度的变化。

实验结果

研究问题

• RQ1 与无力场假设相比，引入有限等离子体β值如何改变弛豫后日冕磁场的最终平衡态？
• RQ2 基于线性绝热演化的解析预测在MHD弛豫的非线性区域内有多大的适用性？
• RQ3 粘性耗散在确定最终平衡态中起什么作用？它是否仅影响弛豫 timescale，而不影响最终状态？
• RQ4 在有限β值等离子体中，弛豫过程中能量如何在等离子体与磁场之间转移？
• RQ5 最终磁流体静力平衡是否可由初始扰动与总压强唯一确定，而与初始扰动的类型无关？

主要发现

• 最终平衡态为非无力场磁流体静力态，其中等离子体压强梯度与磁洛伦兹力平衡，无残余磁张力。
• 基于绝热演化与总压强守恒的解析预测与数值结果高度吻合，即使在初始扰动幅度较大、远超线性区的情况下亦成立。
• 最终等离子体密度与压强分布仅取决于初始扰动与总压强，而与初始扰动的性质（密度、温度或内能变化）无关。
• 粘性仅影响弛豫 timescale，不影响最终平衡态，只要粘性阻尼足以耗散波能。
• 非理想效应相对于初始扰动增长缓慢，因此即使在强非线性区域，线性预测仍保持有效性。
• 在弛豫过程中有显著量的能量从等离子体转移至磁场，表明日冕环境中存在显著的热力学耦合，尤其在磁零点附近。