[论文解读] Sunspot simulations with MURaM -- I. Parameter study using potential field initial conditions

Markus Schmassmann, N. Bello González+2 more|arXiv (Cornell University)|Jan 30, 2026

Solar and Space Plasma Dynamics被引用 0

一句话总结

该论文使用 MURaM 代码在潜在场初始条件下进行日斑模拟的参数研究，改变初始场强、通量内容、盒子尺寸和分辨率，以评估日斑及晕轮的形成并与观测进行比较。

ABSTRACT

Context. Existing sunspot simulations fail to reproduce the observed magnetic field distribution due to an artificially increased $B_{hor}$ at the upper boundary. Aims. We explore alternative ways to better reproduce the magnetic and dynamic properties of observed sunspots. Methods. We use the radiative MHD code MURaM. As initial conditions, we placed a potential magnetic field into small-scale dynamo simulations and used potential-field extrapolation at the top. Results. We find that: (1) Simulations with increasing initial magnetic field strengths (20, 40, 80, and 160 kG) show larger spots, umbrae, and penumbrae. (2) The penumbral-to-spot sizes are smaller than those measured in observed sunspots. (3) In none of the runs are pure Evershed (radially outward) flows. Instead, bi-directional flows with inflows in the inner penumbra and outflows in the outer penumbra were measured, similar to early observations of penumbra formation for runs with $\ge80$ kG at 96/32 km resolution, whereas runs with 40 kG or less showed pure inflows. (4) Simulations with 160 kG at 32/16 km resolution contain filaments with bi-directional and Evershed flows. (5) Simulations with fluxes $>10^{22}$ Mx show unrealistically strong fields in the umbra. (6) The best runs with 160 kG and $10^{22}$ Mx give realistic radial profiles of $B_z$ and $B_r$, although stronger fields than observed. (7) Increasing the width of the box and reducing the overall flux by subtracting a uniform opposing vertical field have little influence on internal spot dynamics and fields, but change the mean vertical field outside the spot. Conclusions. Simulations of small ($10^{22}$) sunspots with an initial potential field and intensified bottom magnetic field strength best reproduce observations of the initial stages of sunspot formation. Numerical resolution may be critical for achieving fully developed penumbrae.

研究动机与目标

• 评估潜在场初始条件是否比以往设置更能产生更现实的日斑磁场与流场性质。
• 确定初始底部场强（B0）、总通量、盒子大小和分辨率如何影响日斑/晕轮的大小和动力学。
• 将模拟的磁场剖面与流场模式与观测基准在晕轮形成阶段进行比较。
• 明确在何种条件下模拟能够再现观测到的艾弗斯（Evershed）样外流和双向晕轮流动。

提出的方法

• 使用带有开放底边界和潜在场顶边界条件的 MURaM 放射性 MHD 代码。
• 在底部施加初始轴向磁场 B0，并用高斯或均匀对立通量 Bopp 以设定底部的初始 Bz(r)，并外推为潜在场。
• 在不同运行中改变参数：B0 = 20、40、80、160 kG；FGauss = 1e22 Mx（或带 Bopp 的 Ftot）；盒宽为 49.152 或 98.304 Mm；分辨率为 96/32 km 或 32/16 km；以及若干 Bopp 值（0–300 G）。
• 在 τ=1 时计算方位平均量（I、Bz、Br、γ），从强度等值线提取黑子/斑点边界，并分析径向流 vz、vr 及其符号混合行为。
• 对活动进行时间平均拟合，以获得代表性状态值，并与观测的日斑（如 NOAA AR 11591）进行比较。
• 在高分辨率运行中包含非灰朗 Radiative Transfer 以实现正向建模和仪器效应。

Figure 1: Quadrants of bolometric intensity maps for four selected simulation runs with different initial magnetic field strengths. From top-left clock-wise: $B_{0}=20$ kG, 40 kG, 80 kG and 160 kG. The magnetic flux is the same for all four runs: $F=10^{22}$ Mx.

实验结果

研究问题

• RQ1 从潜在场初始条件出发并加强底部磁场，日斑与晕轮属性是否比传统设置更接近观测？
• RQ2 初始场强、总通量、域大小和数值分辨率如何影响日斑尺寸、晕轮范围和流动模式（类似艾弗斯德外流与内流）在晕轮形成过程中的表现？
• RQ3 模拟斑点的磁场剖面（Bz、Br）及其径向依赖性为何，与观测基准有何比较？
• RQ4 在何种条件下模拟能产生双向晕轮流动和真实的晕轮丝状结构？
• RQ5 更高分辨率或更大域是否缓解或加剧晕轮发展和通量集中度的差异？

主要发现

• 随着初始底部磁场强度 B0 的增加，晕轮的范围增大；在 20 kG 时没有晕轮，在 160 kG 时出现长而细的丝状结构。
• FGauss > 1e22 Mx 在后期阶段导致日点垂直场过强，难以现实。
• 在 B0 ≥ 80 kG 的模拟中出现双向晕轮流动，内晕轮显示 inward 流、外晕轮显示 outward 流，且在高分辨率下表现更明显。
• 高分辨率运行（32/16 km）产生具有双向与艾弗斯德样流动的丝结构；而低分辨率运行则显示更强的内流和不太现实的晕轮模式。
• 超过 1e22 Mx 的通量会导致日冕中的 Bz 过强，提示在这些设定下存在现实性模拟的上限。
• 在 160 kG、1e22 Mx（及相关高分辨率）运行中获得与半径相关的 Bz、Br 分布的最佳匹配剖面，尽管总体场强仍可能强于太阳观测。
• 更宽的计算盒或减去均匀相对垂直场对内部斑点动力学影响有限，但会改变外部平均垂直场。

Figure 2: Azimuthally averaged intensity $I/I_{\textrm{qs}}$ (solid lines) as a function of radius in Mm and fraction of the spot radius. The colours indicate the simulations listed in the inset legend. For the older runs (with names ending with 'o', red), the bolometric intensity is shown, and for