  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.
