In this paper, a numerical modeling of the impact, spreading, and eventually rebound of a viscoelastic droplet is reported. The numerical model is based on the volume of fluid (VOF) method coupled with the FENE-CR constitutive equations, and accounts for both the surface tension and the substrate wettability. The FENE-CR constitutive equations are used to model the polymer solution, while taking advantage of its rheological characterization. The comparison is performed between droplets of Newtonian solvent and a monodisperse polymer solution. The droplet impact on both hydrophilic and superhydrophobic substrate is analyzed through a detailed analysis of the spreading diameter evolution. It is found that while the droplet kinematic phase seems independent of the substrate and fluids properties, the recoiling phase is highly related to all of them. In addition, the model infers a critical polymer concentration above which the droplet rebound from a superhydrophobic substrate is suppressed. The simulation is of particular interest to ink-jet processing, and demonstrates the capability of the model to handle complex non-Newtonian droplet dynamics.