Fast Legendre spectral method for computing the perturbation of a gradient temperature field in an unbounded region due to the presence of two spheres

The temperature distribution around two spheres is considered when the main field has a constant gradient at infinity. Bispherical coordinates are used, together with a transformation of the dependent variable that leads to separation of variables. Then the solution can be sought in Legendre series with respect to one of the bispherical coordinates. An important element of the proposed work is the effective way to reduce an essentially 3D problem to a set of three 2D problems. The Legendre spectral method is shown to have an exponential convergence which is confirmed by the computations. The efficiency is so high that even for the hard cases of two closely situated spheres, an accuracy of 10^?10 is achieved with as few as 20 terms in the expansion. Solutions with both longitudinal and transverse gradients at infinity are obtained, and the contour lines of the temperature field are presented graphically. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010