Fully implicit higher-order schemes applied to polymer flooding