Data-constrained Magnetohydrodynamic Simulation of a Filament Eruption in a Decaying Active Region 13079 on a Global Scale

Abstract

Filaments are special plasma phenomena embedded in the solar atmosphere, characterized by unique thermodynamic properties and magnetic structures. Magnetohydrodynamic (MHD) simulations are useful to investigate the eruption mechanisms of filaments. We conduct a data-constrained zero-$\beta$ MHD simulation in spherical coordinates to investigate a C3.5 class flare triggered by an eruptive filament on 2022 August 15 in a decaying weak active region 13079. We reconstruct the three-dimensional coronal magnetic field using vector magnetograms and synoptic maps from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI). We transform vector magnetic field into Stonyhurst heliographic spherical coordinates combined with a synoptic map and constructed a potential field source surface (PFSS) model with a magnetic flux rope (MFR) embedded using the Regularized Biot--Savart Laws (RBSL). Subsequently, we conduct a spherical zero-$\beta$ MHD simulation using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC) and replicated the entire dynamic process of the filament eruption consistent with observations. With the calculation of time-distance profile, Qusai-Separatrix Layers (QSL), and synthetic radiation from simulated current density, we find a good agreement between our simulation and observations in terms of dynamics and magnetic topology. Technically, we provide a useful method of advanced data-constrained simulation of weak active regions in spherical coordinates. Scientifically, the model allows us to quantitatively describe and diagnose the entire process of filament eruption.