The main objective of the present work was to investigate computationally the flow and the dispersion of atmospheric pollutants over three dimensional complex topographies in laboratory scale. The investigations included the numerical simulation on the neutral and stably stratified flows over hills and flat terrain. The mathematical model was based on the solution of the general conservation equations and included the Reynolds stress model for turbulence and a two layer zonal model for the flow treatment in the near wall region. The commercial code Fluent (Version 6.0.12), which is based on the finite volume method, was employed in the computational simulations. The numerical results were compared to data obtained in wind tunnel experiments, available in the literature. Comparisons were also made with results obtained by employing the standard (k less épsilon) model for turbulence. The comparisons between the experimental data and the numerical results showed that the combined use of the Reynolds stress model and the two layer treatment provided the best results for the flow representation. This modeling approach was particularly superior in representing the flow recirculation on the leeside of the hill. The predicted concentrations results were reasonably good at regions far away from the emission source. In the near source regions, the ground level concentrations were overestimated by the numerical modeling. These discrepancies were attributed to the employment of an isotropic turbulent diffusivity model in the turbulent dispersion calculations. Nevertheless, the calculated concentration fields represented well important qualitative features of PUC Rio.