The steady motion of a rotating sphere is analysed through two contrasting viscoelastic models, a constant viscosity (FENE-CR) model and a shear-thinning (LPTT) model. Giesekus (Rheol. Acta 9:30-38, 1970) presented an intriguing rotating viscoelastic flow, which to date has not been completely explained. In order to investigate this flow, sets of parameters have been explored to analyse the significant differences introduced with the proposed models, while the momentum-continuity-stress equations are solved through a hybrid finite-element/finite volume numerical scheme. Solutions are discussed for first, sphere angular velocity increase (Ω), and second, through material velocity-scale increase (α). Numerical predictions for different solvent-ratios (β) show significant differences as the sphere angular velocity increases. It is demonstrated that an emerging equatorial anticlockwise vortex emerges in a specific range of Ω. As such, this solution matches closely with the Giesekus experimental findings. Additionally, inside the emerging inertial vortex, a contrasting positive second normal stress-difference (N 2 (γ ) = τ rr − τ θθ ) region is found compared against the negative N 2 -enveloping layer.