Selective Catalytic Reduction (SCR) of NOx through injection of Urea-Water-Solution (UWS) into the hot exhaust gas stream is an effective and extensively used strategy in internal combustion engines. Even though actual SCR systems have 95-96% de-NOx efficiency over test cycles, real driving emissions of NOx are much higher, hence proving that there is room for improvement. The efficiency of the NOx conversion is highly dependent on the size of UWS droplets and their spatial distribution. These factors are, in turn, mainly determined by the spray characteristics and its interaction with the exhaust gas flow. The main purpose of this study is to numerically investigate the sensitivity to the modelling framework of the evaporation and mixing of the spray upstream of the catalyst. The dynamics of discrete droplets is handled through the Lagrangian Particle Tracking framework, with models that account for droplet breakup and coalescence, turbulence effects, and water evaporation. All simulations have been run in the commercial CFD code Ansys Fluent 16.0. Experimental validation of droplet size distribution is carried out through PDPA measurements. Through the numerical study we have identified the most suitable modelling setup that provides accurate results with a competitive computational cost. Results also show the importance of accounting for the effects of evaporation and turbulent fluctuations in the droplet phase.