Based on the Reynolds-Averaged Navier-Stokes (RANS) method, this work has developed a combinational turbulence model for the numerical simulation of turbulent slot jet impingement heat transfer. The combinational model is constructed by taking the original Shear Stress Transport (SST) turbulence model as the parent model and coupling with the Kato-Launder model, the intermittency transition model and the crossflow transition model, which has been abbreviated as SST-KIC model. Detailed study has been conducted for the effect of different component models in terms of the turbulent kinetic energy, the local Nusselt number as well as the flow structures including the mean velocity, skin friction and the law of wall. By comparing with both experimental data and available numerical results of other researchers, it has been found that the proposed model possesses good performance in terms of different parameters for both two nozzle-plate spacing of 4 and 9.2. Finally, the effect of Reynolds number on the flow field characteristics was studied and it is noted that the positive gradient of turbulent intensity has similar variation tendency as the positive gradient of the Nusselt number. (C) 2018 Elsevier Ltd. All rights reserved.