In this paper the Kirchhoff transformation for shape sensitivity analysis of bodies with steady-state non-linear heat conduction is considered. The approach reduces the non-linear 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.