Calculating the potential of 2 spherical perfect conductors with the image method

I am searching for a way to calculate the potential on the surface of two perfect conductors that are spheres. I am not sure my method is correct.

Here is a diagram of what I am studying : They both have a radius (respectively $R_1$ and $R_2$) and are separated by a distance $D$.

I suppose they have respectively charges $Q1$ and $Q2$.

Using the method of images we have the image charges inside the spheres:

• $Q_1$ has an image $q'_1$ located at $O_2 - ( \frac{R_2^2}{D} )$ with $q'_1 = Q1 \frac{-R_2}{D}$
• $Q_2$ has an image $q'_2$ located at $O_1 + ( \frac{R_1^2}{D} )$ with $q'_2 = Q1 \frac{-R_1}{D}$

These images have also images charges in the other sphere :

• $q'_1$ has an image $q''_1$ located at $O_1 + ( \frac{R_1^2}{distance(q'_1, O_1)} )$ with $q''_1 = q'_1 \frac{-R_1}{distance(q'_1, O_1)}$
• $q'_2$ has an image $q''_2$ located at $O_2 - ( \frac{R_2^2}{distance(q'_2, O_2)} )$ with $q''_2 = q'_2 \frac{-R_2}{distance(q'_2, O_2)}$

We can keep going on until it converges

In order to conserve the charges on the surface of the spheres $Q_1$ and $Q_2$ we must place a charge at the centre of the spheres equal to :

• $Q_1-q'_2 - q''_1 - ...$ at the center of $O_1$

• $Q_2-q'_1 - q''_2 - ...$ at the center of $O_2$

The potential is this scenario at any point is :

$V(M) =\frac{1}{4\pi\epsilon_0} (\frac{Q_1-q'_2 - q''_1}{d(M,O_1)} +\frac{Q_2-q'_1 - q''_2}{d(M,O_2)}+ \frac{q'_2}{d(M,q'_2)} + \frac{q'_1}{d(M,q'_1)} + \frac{q''_2}{d(M,q''_2)} + \frac{q''_1}{d(M,q''_1)}) + ...$

This can be done for any number of images using a computer script (python)

Is the correct way to calculate the potential on the spheres to take the above formula and apply it to a point on the surface of the spheres ?

Comparing analytic and numerical results (Edit)

I've set up the converging series in a Python script, I then compared the results for different radius of the electrodes to a numerical calculation using Comsol Myltiphysics.

Here are the results : The distance in between the electrodes is 2 m. The charge $Q_1 = 1 C$ and $Q_2 = 0 C$. I then calculated the potential on the surface of the first electrode.

For small radius's the two methods give very similar results but when the radius gets to big the results diverge. An there is the little drop at the end that also bothers me. This seems to indicate that I am forgetting something but i can't see what.

Convergence (Second edit) As we can see it takes longer for it to converge with a greater radius than a smaller one

Thanks very much for your help

• Do you evaluate your series with fixed number of terms or to a certain precision? – Azad May 7 '15 at 8:44
• With a fixed number, i've tried 100 and 10000. but it seems to converge quite quickly. – Anthony Lethuillier May 7 '15 at 9:03
• That's true. But I expect slower convergence for bigger spheres. Was the result for 2.0 m case the same for 100 and 10000? – Azad May 7 '15 at 9:09
• It converges slower for bigger spheres (see new edit of question) – Anthony Lethuillier May 11 '15 at 13:16
• I think you should also check COMSOL documentation to find out how much is the default precision or number of terms evaluated and whether you can increase it – Azad May 11 '15 at 13:26