For a thermal isotropic body, the first and second variations of an arbitrary thermal functional corresponding to variation of material parameters is derived by using both direct and adjoint approaches. Both the steady-state and transient cases are considered, and the analysis is equally applicable to analytical or numerical solution techniques. A simple example is used to demonstrate the various approaches.