Exploring the effects of numerical methods and slope limiters in heliospheric modeling
Tinatin Baratashvili,Christine Verbeke,Rony Keppens,Stefaan Poedts
2023-05-24
Abstract:Coronal mass ejections (CMEs) are large eruptions close to the solar surface, where plasma is ejected outwards into space at large speeds. When directed towards Earth, they interfere with Earth's magnetic fields and cause strong geo-effective storms. In order to mitigate the potential damage, forecasting tools are implemented. Recently, a novel heliospheric modelling tool, Icarus, has been implemented, which exploits the open-source framework MPI-AMRVAC as its core MHD solver. This new model efficiently performs 3D MHD simulations of the solar wind and the evolution of interplanetary CMEs with the help of advanced techniques, such as adaptive mesh refinement and gradual radial grid stretching. The numerical methods applied in the simulations can have significant effects on the simulation results and on the efficiency of the model. In this study, the effect of different combinations of numerical schemes and slope limiters, for reconstructing edge-based variabes used in fluxes, is considered. We explore frequently exploited combinations from the available numerical schemes in MPI-AMRVAC: TVDLF, HLL and HLLC along with the slope limiters 'woodward', 'minmod', 'vanleer', and 'koren'. For analysis purposes, we selected one particular solar wind configuration and studied the influence on variables at 1 AU in the equatorial plane. The goal is to find the optimal combination to produce accurate results fast and in a robust way so that the model can be reliable for day-to-day use by space weather scientists. As a conclusion, the best result assessed with these two criteria is the combination of the TVDLF scheme with the 'woodward' limiter.
Solar and Stellar Astrophysics,Computational Physics,Plasma Physics,Space Physics