A flux-limited numerical method for the MHD equations to simulate propulsive plasma flows