Simulating the photospheric to coronal plasma using magnetohydrodynamic characteristics III: validation including gravity, flux emergence, and an eruption

TL;DR

This paper validates a data-driven MHD simulation approach for modeling solar active region emergence and eruptions, demonstrating high accuracy and stability over the entire active region lifecycle, aligning well with observational cadences.

Contribution

The study introduces a validated, stable MHD simulation method driven by observational data that accurately reproduces solar magnetic flux emergence and eruptions.

Findings

Simulation accurately reproduces active region emergence and eruptions.

Total emerged energy matches ground truth within 1%.

Method remains stable over entire active region lifetime.

Abstract

Solar eruptions arise from instabilities or loss of equilibria in the solar atmosphere, but routinely inferring the precise magnetic and plasma properties that lead to eruptions is not currently practical using synoptic solar observations. Data driven simulations offer an appealing alternative. We test our boundary data-driven magnetohydrodynamic (MHD) approach, based on the method of characteristics, on a simulation that includes full MHD, a stratified atmosphere, and the emergence of a model solar magnetic active region, from the photosphere upwards. The driven simulation is tested against a larger, ab initio ``Ground Truth'' simulation that extends downward into the convection zone. Our driven simulation accurately reproduces the dynamic emergence of the active region above the photosphere, the formation of key topological features throughout the corona, and the subsequent eruption of mass and magnetic field. The total emerged energy matches to better than one percent, the ratio of emerged to eruptive energy is $\approx2\%$, and the actual values of each energy term agree to within $10\%$ between the two cases. Crucially, the data injection cadence, when properly scaled, matches the cadence of synoptic observations of the Sun's surface magnetic field, and is three to four orders of magnitude longer than the inherent CFL time step of the simulations. The stability of the code and fidelity of the results over an entire active region lifetime, from emergence to eruption, strongly suggests that our method will produce reliable results when driven using solar synoptic observations from existing and anticipated ground and spaced based observatories.