Data-driven Radiative Magnetohydrodynamics Simulations with the MURaM Code: the Emergence of Active Region 11158 and the X2.2 Flare

We present the application of the data-driven branch of the MURaM code to the extensively studied flare-productive active region 11158. We refine the hybrid model strategy, which was described in the earlier paper of this series, to model the emergence of the active region during 4 solar days starting shortly before 2011 February 11 and the eruption of an X2.2 flare on February 15. After 4 days of evolution, a major eruption of a magnetic flux rope occurs in the simulation at approximately 3 hours (3% difference) before the real flare. The eruption leads to magnetic reconnection that contributes to bulk heating in the chromosphere and corona. The deposition of flare energy in the chromosphere causes strong condensations and evaporations, which fill hot post-flare loops and bright flare ribbons that exhibit separation and extension similar to the observed ribbon evolution. The synthesized soft X-ray flux corresponds to X class, which is close to the real event. The upward eruption of the flux rope leads to a piston-driven shock and horizontal expansion that exerts a strong downward impact on the lower atmosphere and generate an apparently fast-propagating chromospheric Moreton wave. We conclude that the data-driven radiative simulation of this active region can reproduce the key observational results of the real flare and demonstrate the great potential of this method for studying solar eruptions in a realistic corona environment.

💡 Research Summary

This paper presents a comprehensive data‑driven radiative magnetohydrodynamic (RMHD) simulation of solar active region (AR) 11158 and its X2.2 flare on 2011 February 15, using the MURaM code. Building on earlier work that introduced a data‑driven bottom boundary for MURaM, the authors refine a three‑stage hybrid modeling strategy to capture the full evolution from flux emergence to eruption.

In the first stage, a zero‑β MHD run (256 × 256 × 1152 grid) is driven by observed photospheric vertical magnetic field (SDO/HMI) with a time‑acceleration factor of 12. The horizontal electric field needed for the induction equation is derived from the observed ∂Bz/∂t and a rotation parameter Ω, which is estimated from the shear motion of the two main sunspot pairs (N1–P2). An Ω of 1.0 × 10⁻⁵ s⁻¹ is found to be sufficient to inject free magnetic energy; control runs with Ω = 0 or half this value produce no eruption.

The second stage is a full radiative MHD simulation (the “evo. run”) that uses the actual physical size of AR 11158 (512 × 512 × 1920 grid, 540 km × 64 km spacing). The initial magnetic field is constructed by adding a potential field derived from the observed Bz, the non‑potential component extracted from the zero‑β snapshot, and a quiet‑Sun background field. Plasma density, internal energy, and velocities are taken from a relaxed quiet‑Sun simulation and scaled to match the larger domain. No time‑acceleration is applied in this stage, allowing realistic coronal heating and plasma response.

The third stage (“flare run”) extends the vertical domain to 245 Mm (512 × 512 × 3840 grid) to accommodate the upward motion of the erupting flux rope. It is initialized from the evo. run snapshot at 300 000 iterations (≈ Feb 14 23:05 UT). The upper half of the domain is filled with a potential field and hydrostatic plasma stratification, while the lower half inherits the fully evolved state from the evo. run.

Results show that the total magnetic energy grows to ≈ 4 × 10³³ erg, with free magnetic energy reaching ≈ 6.5 × 10³² erg. Around 3 hours before the observed flare (a 3 % timing error), the free energy drops sharply, releasing ≈ 1.94 × 10³² erg over ~1000 s, marking the simulated eruption. The flux rope undergoes a slow rise followed by a rapid acceleration; synthetic GOES‑15 soft X‑ray flux peaks at 1.97 × 10⁻⁴ W m⁻², corresponding to an X2 class flare, closely matching the observed X2.2 event.

The eruption drives magnetic reconnection that heats the chromosphere and corona, producing strong condensations and evaporations. These fill post‑flare loops and generate bright flare ribbons whose separation and extension reproduce the observed EUV ribbon evolution. The upward motion of the flux rope creates a piston‑driven shock that expands horizontally and exerts a strong downward impact on the lower atmosphere, generating a fast‑propagating chromospheric Moreton wave.

The authors discuss the strengths of the approach: (i) realistic coupling of observed photospheric magnetic evolution to a full RMHD corona, (ii) ability to capture both large‑scale CME dynamics and fine‑scale chromospheric responses, and (iii) quantitative agreement with key observables (timing, energy release, X‑ray flux, ribbon motion, Moreton wave). Limitations include the absence of non‑thermal particle physics, reliance on an empirically chosen Ω, and the use of a time‑acceleration factor only in the zero‑β stage.

In conclusion, the data‑driven radiative MHD simulation successfully reproduces the essential observational signatures of the AR 11158 X2.2 flare, demonstrating that such models can serve as powerful tools for investigating solar eruption mechanisms and for future predictive capabilities.