Exploring various flux vector splittings for the magnetohydrodynamic system

Our FVS successfully captures isolated, stationary contact discontinuities in MHD. However, we explain why any FVS for MHD is not adept at capturing isolated, stationary Alfvenic discontinuities. Several stringent one-dimensional Riemann problems are presented to show that the method works successfully and can effectively capture the full panoply of wave structures that arise in MHD. This includes compound waves and switch-on and switch-off shocks that arise because of the non-convex nature of the MHD system.