Authors: Nian Liu (NJIT), Satoshi Inoue (NJIT), Ju Jing (NJIT, BBSO), Haimin Wang (NJIT, BBSO)

We conducted both data-constrained and data-driven magnetohydrodynamic (MHD) simulations of solar active region 13663, which produced two X1.3-class flares and one M7.5-class flare on May 5th, 2024. The simulations employed the zero-beta MHD approximation, with initial conditions derived from the nonlinear force-free field (NLFFF) extrapolated from photospheric magnetograms taken one hour before each flare, specifically for the data-constrained simulations. The first NLFFF at 04:48 UT also served as the initial condition for the data-driven simulation. Anomalous resistivity was applied in the data-constrained simulations to locally enhance magnetic reconnection. In data-driven simulation, the bottom magnetic fields are driven by the electric field that derived from 12 minutes cadence hmi.sharp_cea magnetograms. The resultant 3D magnetic field structures were then compared with AIA EUV observations. Our findings indicate that all three data-constrained simulations consistently reproduced the magnetic flux rope eruption.