Induction hardening of gear wheels is modeled. The model consists of two nonlinear partial
differential equations describing the distributions of magnetic and temperature fields
in the system. All material parameters are supposed to be functions of temperature. The
model is then solved numerically in the 3D hard-coupled formulation using the professional
code FLUX3D supplemented with a number of own scripts and procedures. The
methodology is illustrated with a typical example whose results are discussed.