In this study, an analytical solution to the heat conduction equation in an annulus between eccentric spheres with isothermal boundaries and with heat generation is obtained using Green's function method. Deriving Green's function in bispherical coordinates for eccentric spherical annuli, an exact solution of heat equation for eccentric spheres with constant surface temperatures is expressed in terms of Green's function and source distribution. The solution is general and can easily be applicable to any space dependent heat source. Results are presented as temperature distributions and local Nusselt numbers for sources having practical and theoretical importance: uniform source, impulse source, shell source. Analytical results are compared with the results of the Computational Fluids Dynamics (CFD) solver Fluent and perfect agreement is observed. (C) 2013 Elsevier Masson SAS. All rights reserved.