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