[论文解读] Evolution of Relative Magnetic Helicity: New Boundary Conditions for the Vector Potential

Shangbin Yang, Joerg Buechner+2 more|MPG.PuRe (Max Planck Society)|Dec 26, 2017

Solar and Space Plasma Dynamics参考文献 34被引用 4

一句话总结

该论文提出了一种新的边界条件，用于在有限体积内计算矢量势，无需在所有边界上保持通量平衡，从而实现开放系统中相对磁螺旋度的精确计算。通过采用库仑规范和基于正弦函数的边势插值方法，该方法在验证中误差低于1.2%，在AR 11429的模拟中仅出现8%的螺旋度耗散和2%的逃逸，与Woltjer-Taylor猜想及磁云螺旋度输运的观测约束一致。

ABSTRACT

We recently proposed a method to calculate the relative magnetic helicity in a finite volume for a given magnetic field which however required the flux to be balanced separately on all the sides of the considered volume. In order to allow finite magnetic fluxes through the boundaries, a Coulomb gauge is constructed that allows for global magnetic flux balance. We tested and verified our method in a theoretical fore-free magnetic field model. We apply the new method to the former calculation data and found a difference of less than 1.2\%. We also applied our method to the magnetic field above active region NOAA 11429 obtained by a new photospheric-data-driven MHD model code GOEMHD3. We analyzed the magnetic helicity evolution in the solar corona using our new method. It was found that the normalized magnetic helicityis equal to -0.038 when fast magnetic reconnection is triggered. This value is comparable to the previous value (-0.029) in the MHD simulations when magnetic reconnection happened and the observed normalized magnetic helicity (-0.036) from the eruption of newly emerging active regions. We found that only 8\% of the accumulated magnetic helicity is dissipated after it is injected through the bottom boundary. This is in accordance with the Woltjer conjecture. Only 2\% of magnetic helicity injected from the bottom boundary escapes through the corona. This is consistent with the observation of magnetic clouds, which could take away magnetic helicity into the interplanetary space, in the case considered here, several halo CMEs and two X-class solar flares origin from this active region.

研究动机与目标

• 消除在有限体积中计算相对磁螺旋度时对所有边界通量平衡的严格要求。
• 开发一种规范不变的矢量势重构方法，允许边界存在有限磁通量。
• 利用真实的MHD模拟研究太阳日冕活动区中的螺旋度耗散与逃逸。
• 在太阳日冕动力学背景下，检验Woltjer-Taylor关于螺旋度守恒的猜想。

提出的方法

• 构建库仑规范，以确保尽管边界存在有限磁通量，仍能保持全局磁通量平衡。
• 使用正弦函数对位于活动区上方的长方体体积的12条边上的矢量势值进行插值。
• 利用理论力自由磁场模型（Low & Lou, 1990）和Yang et al. (2013) 的先前模拟数据对方法进行验证。
• 将该方法应用于GOEMHD3 MHD模型的磁场数据，该模型利用光球层数据驱动日冕磁场演化。
• 使用Finn-Antonsen公式结合新矢量势计算相对磁螺旋度，确保规范不变性。
• 分析螺旋度演化，重点关注磁重联和CME事件期间的注入率、耗散率和逃逸率。

实验结果

研究问题

• RQ1 当磁通量在所有边界上未平衡时，是否仍能准确计算有限体积内的相对磁螺旋度？
• RQ2 在太阳爆发期间，注入的磁螺旋度中有多大比例发生耗散或逃逸至行星际空间？
• RQ3 在真实日冕模型中，快速磁重联期间归一化的磁螺旋度（H/Φ²）如何演化？
• RQ4 模拟结果在多大程度上支持Woltjer-Taylor关于准理想MHD系统中螺旋度守恒的猜想？
• RQ5 计算得到的螺旋度值与磁云和CME的观测估计值相比如何？

主要发现

• 与先前的通量平衡方法相比，该新方法的相对误差低于1.2%，验证了其准确性。
• 在太阳黑子群NOAA 11429的模拟中，仅有8%的注入磁螺旋度发生耗散，与Woltjer-Taylor猜想一致。
• 通过光球层边界注入的磁螺旋度中仅有2%逃逸至日冕，表明螺旋度输运效率较低。
• 磁重联被触发后，归一化磁螺旋度（H/Φ²）达到-0.038，与爆发活动区观测值-0.036相当。
• 结果支持磁螺旋度在很大程度上被守恒且未被高效损失的观点，即使在多次 halo CME 和 X级耀斑期间亦如此。
• -0.038 的值比 kink/磁扭不稳定性理想MHD模拟结果小一个数量级，表明异常电阻率可能在增强耗散中起作用。