Numerical investigation of multistability in the unstable flow of a polymer solution through porous media