The main objective is to develop an efficient BEM scheme for the numerical solution of two‐dimensional heat problems. Our scheme will be of the re‐initialization type, in which the domain integrals are computed by a recursion relation which depends only on the boundary temperature and flux at previous time step. To obtain the re‐initialization approach, we will use in the integral representation formula a Green function corresponding to zero temperature in a box containing the original domain, instead of using the classical free space fundamental solution. This Green function is given in terms of the original fundamental solution plus a regular solution of the heat equation inside the domain under consideration. It can therefore be used in the integral representation formula of the heat equation (direct formulation) to obtain the solution of a heat problem in such a domain. The Green function mentioned can be obtained by the images method, and the resulting source series can also be rewritten in terms of a double Fourier series, that we will use in the domain integral of the integral representation formula to transform such integral into equivalent surface integrals.