# Charge Density of a Conducting Disc

1. Jul 7, 2006

### Blistering Peanut

Oh since this got moved: It's not for any course or class so I'm not following any book.

I'm not sure if anyone will be able to answer this but I'll ask anyway:

I've solved laplaces equation in cyclindrical coordinates for a disc of radius 'a' and constant potential V on the disk (disc in z=0 plane, centered on z axis) and got this as the potential everywhere

$$\phi(r,z) = \frac{2V}{\pi} \arcsin{\frac{2 a}{\sqrt{z^2 + (r + a)^2} + \sqrt{z^2 + (r-a)^2}}}$$

I want to now find the charge density on the disc. From Gauss and the fact that at the surface the electric field will be perpendicular to the surface we have $$E_{z} = 4 \pi \sigma$$ where sigma is the charge density.

So I differentiate Phi w.r.t z to get the z component of the electric field and take the limit as z->0. I get zero as my answer.

What am I doing wrong? I can only get this to work in oblate spheroidal coordinates. They are very convenient for this problem but not so when I try to generalize to more discs of different radii and potential.

Last edited: Jul 7, 2006
2. Jul 7, 2006

### lalbatros

Hi BP,

It would also be useful to plot a numerical example.

Michel

Last edited: Jul 7, 2006
3. Jul 7, 2006

### Blistering Peanut

Dear Michel,

Last edited: Jul 7, 2006
4. Jul 7, 2006

### nrqed

I am sure this is a really stupid question and I am missing something obvious, but if one sets z=0 in your solution, one gets a function of r. But you stated that the potential is constant on the disk. I am sure I am missing something obvious! (is your solution for r>a??)

5. Jul 7, 2006

### Blistering Peanut

Do a taylor expansion of what is left inside the arcsin for r<a and you will end up with 1. arcsin[1]=Pi/2

6. Jul 7, 2006

### nrqed

Oh, Ok I see. Actually, that was obvious even without a Taylor expansion since, for r<a, one sees right away that the r's cancel out. Sorry!

Would it be too much of a mess to show your result for the derivative?

7. Jul 7, 2006

### Blistering Peanut

Well I wasn't doing anything wrong with the derivative since looking at some graphs Ez does go to zero as z->0.

8. Jul 7, 2006

### Blistering Peanut

Here are my contour plots for the electric fields. You can imagine where the disc of radius 1 is.
Ez:
http://www.maths.tcd.ie/~dleen/ez.jpg
horizontal axis is r and vertical z

and Er:
http://www.maths.tcd.ie/~dleen/er.jpg
horizontal axis is r and vertical z

and while I'm at it:
http://www.maths.tcd.ie/~dleen/potential.jpg
z axis is phi, x and y are r and z

I was beginning to think that I would calculate Er...

Last edited: Jul 7, 2006
9. Jul 8, 2006

### lalbatros

BP,

I think you are using Mathematica or some similar program.
These programs must be handled with care when sigularities are involved.
In your problem, the electric field has a discontinuity across the plate.
There are many consequences like:
the limit for z->0 does not exist, it exists on the left and the right only,
the Taylor series does not exist​

In this situation, it is useful to consider special cases before handling the full question.
Here, the analysis on the axis of symetry (r=0) is easily done "manually" and you can spend a little time to understand all aspects of the problem. Eventually you can manage later to find out the charge distribution everywhere, but in a first step calculate it in the center.

In the attachment you will find the fields on the axis of symmetry. Clearly Ez is not zero. Proceeding carefully at other places will give you the charge distribution. Note that in general f(x) diagrams are easier to read than contour plots or 3D's .

Michel

Attachment URL:

.

Last edited: Nov 22, 2011
10. Jul 8, 2006

### Blistering Peanut

I have this so far

$$E_z = \frac{\frac{2 a z}{\sqrt{((r-a)^2 +z^2)((r+a)^2 +z^2)}}}{\sqrt{2(r^2 + z^2 - a^2 + \sqrt{((r-a)^2 +z^2)((r+a)^2 +z^2)})}}$$

I'm unsure how to proceed taking the limit. My first attempt gave

$$E_z = \frac{\frac{2 a z}{\sqrt{((r-a)^2 +z^2)((r+a)^2 +z^2)}}}{\sqrt{2 z^2}}} =\frac{\sqrt{2} a}{\sqrt{(a^2 -r^2)^2}}$$

Along the symetry axis I got 1/a as my answer

Your attachment is still pending approval so I can't see it.

Last edited: Jul 8, 2006
11. Jul 8, 2006

### lalbatros

BP,

There are some braket problems in your expression.

You need to analyse the limit for z->0 carefully.
First, factor out these few things that don't tend to zero.
Then observe the parts that produce the indetermination:
the 2az factor on the numerator
the big square root on the denominator​
Note indeed that the big square root goes to zero as z->0 .
Be careful to choose the right sign when taking the square root of a square!
Once you have observed the indetermination, you should be conviced that the limit is finite.
Patience, care, and a few Taylor series for the big denominator should tell you the answer.

Michel

Postscriptum

Here is the formula I got for Ez² (taking V=1 and a=1):

16*z²/(Pi²*u²*v²*(-2 + u + v)*(2 + u + v))​
with

u = Sqrt((1 + r)² + z²)
v = Sqrt((1 - r)² + z²)

You can check if this is consistent with your formula.

Last edited: Jul 8, 2006
12. Jul 8, 2006

### Blistering Peanut

Yes mine is the same as yours.

What do you mean by a few taylor series for the big denominator? Should I expand it around z=0? I haven't had any experience with this so thanks for bearing with me.

13. Jul 8, 2006

### nrqed

But I thought your problem was that you were getting Ez=0 at z=0,no?

14. Jul 8, 2006

### lalbatros

Hi BP,

I used Mathematica v2.0 to find out Ez.

I had to control each step by myself because Mathematica v2.0 seems very unreliable for these kind of limits. Series expansions are not more reliable.

I worked with the expressions from my last postscriptum.
Essentialy I had to calculate the limit of z²/(-2+u+v) , other factors make no problem.
I expanded the divergent factor (-2+u+v) from the denominator in series of z.
I injected that back in Ez² and this removed the indeterminacy.
Going back to Ez and taking care for the correct sign gives a simple answer.

Michel

Postscriptum

I got the following result: Ez = -2/Pi * 1/sqrt(a²-r²) .
By the way, how did you find the formula for the potential?

Last edited: Jul 8, 2006
15. Jul 8, 2006

### nrqed

EDIT: I see how you got your answer. But it's incorrect. Th emistake is that you neglect z^e relative to a^2 and r^2 in part of the expression (under the square root) while retaining z^2 in another part of th eexpression. This is inconsistent. I work out the correct expnasion below.

I don't get sqrt(2 z^2) for the expansion of the
denominator. Using (I am assuming r<a)
$${\sqrt { (a-r)^2 + z^2 }}\approx a-r + {1 \over 2 } { z^2 \over a-r}$$
$${\sqrt { (a+r)^2 + z^2}} \approx a+r + {1 \over 2 } { z^2 \over a+r}$$
SO the product of the two, for small z, gives
$$a^2 - r^2 + {1\over 2} { a+r \over a-r} z^2 + {1 \over 2} { a-r \over a+r } z^2$$
Therefore, if I add r^2 -a^2 +z^2 to the above, I get
$$z^2 + {1\over 2} { a+r \over a-r} z^2 + {1 \over 2} { a-r \over a+r }z^2 = z^2 + z^2 { a^2 +r^2 \over a^2 -r^2}$$

All this must still be multiplied by two and put under a square root, giving
$${\sqrt {2 z^2}} {\sqrt { { 2 a^2 \over a^2 -r^2}}}$$

I will post and then check that there is no typo.
Edit: it looks right now.

Last edited: Jul 8, 2006
16. Jul 8, 2006

### Blistering Peanut

I solved the laplace equation in cyclindrical coordinates with the appropriate boundry condition. Weber's discontinuous integrals really helped.

Are you using the Series function in mathematica?

Why does the sign matter so much when taking the square root of the square in this problem? This seems like where mathematica is having the problem.

17. Jul 8, 2006

### Blistering Peanut

I think I've got it now. Thank you very much for your help. I have learned a lot today and hopefully I'll be able to apply this to the rest of the problem.

nrged sorry, I was still confused at that stage and I don't know why I posted that.

18. Jul 9, 2006

### pervect

Staff Emeritus
Buried deep in a 56 page postscript file that can be found at
http://www.ttc-cmc.net/~fme/captance.html

the above potential function and the associated charge distribution is given. I've attached this as a .png file.

#### Attached Files:

• ###### cap.png
File size:
2.4 KB
Views:
141
19. Jul 9, 2006

### lalbatros

BP,

Yes I do.

However, it fails if the function has a discontinuity or its derivative has a discontinuity and so on. All derivatives should be continuous, otherwise the series expansion needs further specifications that Mathematica does not handle (on the right or on the left).

Still the functionality was useful for you problem. I did not use it directly on the considered function (Ez), but on elements that can be developped in series, like one the square root. Assembling the expanded parts leads to the result.

Michel

Postscriptum

Here is how I did it with Mathematica.

V = Vo 2/Pi ArcSin[2 a / (Sqrt[z*z+(a+r)*(a+r)] + Sqrt[z*z+(a-r)*(a-r)])]
Ez = Simplify[D[V,z]]
Ezuv = Ez //.{Sqrt[(a+r)^2+z^2]->u,Sqrt[(a-r)^2+z^2]->v}
Ezuv = Ezuv//.{1/Sqrt[(a+r)^2+z^2]->1/u,1/Sqrt[(a-r)^2+z^2]->1/v}
Ezuv2 = Factor[Ezuv^2]
Ezuv2 = Ezuv2//.{u+v->upv}
upv = Sqrt[(a+r)^2+z^2] + Sqrt[(a-r)^2+z^2]
upv = Simplify[Normal[Series[upv,{z,0,2}]//.{Sqrt[(a-r)^2]->a-r,Sqrt[(a+r)^2]->a+r}]]
Ezuv2 = Simplify[(Ezuv2//.{1/u^2->1/((a+r)^2+z^2),1/v^2->1/((a-r)^2+z^2)})//.{z->0}]

Clearly I could not have the answer automaticaly.
It had to make sure that the term I called "upv" was properly handled.
I had also to guide Mathematica with the square root, to be sure it comes with positive sqrt.

I have a very old version of Mathematica.
I would like to know how you did it yourself and the tool you used (paper and pencil ?).

Last edited: Jul 9, 2006
20. Jul 9, 2006

### Blistering Peanut

Yes I have been mainly using pen and paper. I was getting sick of those long expressions (and they double in length for the capacitor) so I started using mathematica. However since I was not used to watching out for singularities etc I started running into problems. I've learned now though which is great.