Induction contour hardening of gear wheels belongs to effective heat treatment technologies especially recommended for high-tech applications in machinery, automotive and aerospace industries. In comparison with long term, energy consuming conventional heat treatment (carburizing and consequent quenching), its main positive features are characterized by high total efficiency, short duration and relatively low energy consumption. However, modeling of the process is relatively complicated. The numerical model should contain both multi-physic and non-linear formulation of the problem. The paper concentrates on the modeling of rapid induction heating being the first stage of the contour induction hardening process which is the time consuming part of the computations. It is taken into consideration that critical temperatures and consequently the hardening temperature are dependent on the velocity of the induction heating. Numerical modeling of coupled non-linear electromagnetic and temperature fields are shortly presented. Investigations are provided for gear wheels made of a special quality steel AISI 300M. In order to evaluate the accuracy of the proposed approach, exemplary computations of the full induction contour hardening process are provided. The exemplary results are compared with the measurements and a satisfactory accordance between them is achieved.