Using a (2+1)-d viscous hydrodynamical model, we study the dependence of flow observables on the collision energy ranging from √ s = 7.7 A GeV at the Relativistic Heavy Ion Collider (RHIC) to √ s = 2760 A GeV at the Large Hadron Collider (LHC). With a realistic equation of state, Glauber model initial conditions and a small specific shear viscosity η/s = 0.08, the differential charged hadron elliptic flow v ch 2 (pT , √ s) is found to exhibit a very broad maximum as a function of √ s around top RHIC energy, rendering it almost independent of collision energy for 39 ≤ √ s ≤ 2760 A GeV. Compared to ideal fluid dynamical simulations, this "saturation" of elliptic flow is shifted to higher collision energies by shear viscous effects. For color-glass motivated MC-KLN initial conditions, which require a larger shear viscosity η/s = 0.2 to reproduce the measured elliptic flow, a similar "saturation" is not observed up to LHC energies, except for very low pT . We emphasize that this "saturation" of the elliptic flow is not associated with the QCD phase transition, but arises from the interplay between radial and elliptic flow which shifts with √ s depending on the fluid's viscosity and leads to a subtle cancellation between increasing contributions from light and decreasing contributions from heavy particles to v2 in the √ s range where v ch 2 (pT , √ s) at fixed pT is maximal. By generalizing the definition of spatial eccentricity ǫx to isothermal hyper-surfaces, we calculate ǫx on the kinetic freeze-out surface at different collision energies. Up to top RHIC energy, √ s = 200 A GeV, the fireball is still out-of-plane deformed at freeze out, while at LHC energy the final spatial eccentricity is predicted to approach zero.