The well-known method of images relates the solution of the heat equation on R n (typically n ¼ 2 or n ¼ 3) to the solution of the heat equation on certain spatial subdomains X of R n : By reformulating the method of images in terms of a convolution kernel, two novel extensions are obtained in this paper. First, the method of images is extended from thermal problems to thermoelastic problems, that is, it is demonstrated how the heat-induced deformations on R n can be related to the heat-induced deformations on certain subdomains X of R n : Secondly, an explicit expression for the convolution kernel for the disk is obtained. This enables the application of the method of images to circular domains to which it could not be applied before. The two obtained extensions lead to a computationally efficient simulation method for repetitive heat loads on a disk. In a representative simulation example of wafer heating, the proposed method is more than ten times faster than a conventional Finite Element approach.