In this study, heat transfer characteristics in the single slot jet impinging cooling process of constant heat flux surface are numerically investigated. It is assumed that the flow is turbulent, two-dimensional and in steady state. Governing equations are solved by using Galerkin finite element method by employing five two-equation turbulence models based on Reynolds-averaged Navier-Stokes (RANS) approach. Although the most satisfactory results are obtained with nonlinear algebraic stress model of Shih-Zhu-Lumley in stagnation region, overall performance of RNG and standard k- models are better in comparison with other models by considering entire region. Subsequent computations are performed with RNG and standard k- models for nozzle to plate spacing and Reynolds numbers in the ranges of 4 <= zlD(h)<= 10 and 4000 <= Re <= 12000, respectively. Also, inlet turbulence intensity and heat flux boundary conditions effects on heat transfer are investigated. Property variation and buoyancy effects are considered to decrease possible discrepancy with experimental results and capture the turbulence intensity effects more accurately. Acceptable agreement with the measured values in published literature are obtained and discussed.