Time-evolving coronal modelling of solar maximum around the May 2024 storm by COCONUT

Coronal simulations of the solar maximum struggle with poor numerical stability and low computational efficiency since the magnetic field is more complex and stronger and coronal structures evolve more rapidly. This paper aims to enhance the numerical stability of the time-evolving COCONUT coronal model to mitigate these issues, to evaluate differences between the time-evolving and quasi-steady-state coronal simulation results, and to assess the impact of spatial resolution on global MHD coronal modelling of solar maximum.After enhancing the positivity-preserving property of the time-evolving COCONUT, we employ it to simulate the evolution of coronal structures from the solar surface to 0.1 AU over two CRs around the May 2024 solar storm event. These simulations are performed on unstructured meshes containing 6.06, 1.52, and 0.38 M cells to assess the impact of grid resolution. We also conduct a quasi-steady-state coronal simulation, treating the solar surface as a rigidly rotating spherical shell, to demonstrate the impact of magnetic flux emergence and cancellation in global coronal simulations. Comparison with observations further validates the reliability of this model.This paper demonstrates that incorporating magnetic field evolution in inner-boundary conditions can significantly improve the fidelity of global MHD coronal simulations around solar maximum. The simulated magnetic field strength using a refined mesh with 6.06 M cells can be more than 40% stronger than that in the coarser mesh with 0.38 M cells. A time step of 5 minutes and the mesh containing 1.5 M cells can effectively capture the evolution of large-scale coronal structures and small-sized dipoles. Thus, this model shows promise for accurately conducting real-time global coronal simulations of solar maximum, making it suitable for practical applications.