  We present novel numerical schemes for ideal magnetohydrodynamic (MHD)
simulations aimed at enhancing stability against numerical shock instability
and improving the accuracy of low-speed flows in multidimensions. Stringent
benchmark tests confirm that our scheme is more robust against numerical shock
instability and is more accurate for low-speed, nearly incompressible flows
than conventional shock-capturing schemes. Our scheme is a promising tool for
tackling MHD systems, including both high and low Mach number flows.
