Summary: In this paper the Kirchhoff transformation for shape sensitivity analysis of bodies with steady-state nonlinear heat conduction is considered. The approach reduces the nonlinear thermal problem to the standard Laplace problem. The variation of a general integral functional is described in terms of adjoint quantities. The design sensitivities are calculated using the material derivative concept. A thermal diffuser problem is considered to illustrate the method proposed.