Borehole heat exchangers (BHE) are the most common solution for ground coupled heat pumps (GCHP). The thermal interactions with ground and the GCHP performance are strictly related to a number of factors, including the ground properties, the BHE field geometry, the sequence and strength of heat loads required by the building. A typical assumption to solve this problem is assuming the ground volume as a purely conductive medium with constant properties: in such a way the transient Fourier equation can be solved either analytically or numerically with a reasonable computational effort provided that a simple BHE geometry is considered. According to this approach suitable temperature response factors can be calculated: they account for constant heat transfer rate conditions and they can be later applied to superposition techniques for simulating transient heat loads to/from the ground. In literature a number of reference geometries and related temperature response factors for BHE modelling are available. They include the infinite and finite line source and the infinite cylindrical source, at imposed heat transfer rate or even imposed temperature. The aim of this paper is to develop a reliable numerical model for generating temperature response factors for the single BHE according to different geometrical and boundary conditions. A detailed analysis is carried out to recast the Fourier equation and to efficiently solve it with suitable numerical parameters. The results are compared with literature analytical data and new insights are provided to the temperature response factor approach for simulating GCHP systems. (C) 2015 Elsevier Ltd. All rights reserved.