Jet impingement has been widely used due to its high heat transfer rate, but numerical simulation of this flow is highly challenging because of complex flow phenomena. In this paper, a turbulence model based on the physics of jet impingement is developed using the shear stress transport (SST) four-equation transition model and the SST with cross-diffusion correction (SSTCD) model. The proposed method is evaluated for plane jet impingement under various nozzle-plate spacings (2.6, 4 and 6) and different Reynolds numbers ranging from 10,200 to 20,000 in terms of heat transfer and flow fields. The results show that the developed model has the ability to capture the second peak of heat transfer and skin friction at low nozzle-plate spacing. Even at high nozzle-plate spacing, this approach also gives accurate trends for the above two features. However, without the cross-diffusion correction, the transition model provides too high energy and low velocity peaks. With the developed model in hand, the effects of pressure gradient on heat transfer and skin friction coefficient are further investigated. It is found that when the adverse pressure gradient becomes strong, the position of the secondary peak of skin friction occurs earlier than that of the secondary peak of heat transfer rate. Conversely, if the adverse pressure gradient disappears, the positions of these features coincide.