Hard Electrostat Prob - Find E-Field at Point h

[Swapnil]

Homework Statement

A hemisphere of radius a centered at the origin carriers a total charge Q uniformly distributed over its curved surface. Find the e-field at a point h on the z-axis.

Homework Equations

d\vec{E} = \frac{1}{4\pi\epsilon} \frac{dq}{R^2}\hat{R}

The Attempt at a Solution

Can someone please tell me if this answer is correct (my professor and I both had a real hard time figuring this one out):

\vec{E}(0,0,h) = \frac{Q}{4\pi\epsilon}\Big(\frac{a}{h^2\sqrt{h^2+a^2}} + \frac{1}{h^2}\Big)\hat{z}

 

[siddharth]

Can someone please tell me if this answer is correct (my professor and I both had a real hard time figuring this one out):

\vec{E}(0,0,h) = \frac{Q}{4\pi\epsilon}\Big(\frac{a}{h^2\sqrt{h^2+a^2}} + \frac{1}{h^2}\Big)\hat{z}

I don't think that's right.

One way of doing this, is to first calculate V(z)[/tex], and then get \vec{E} from that.

 

[lalbatros]

Swapnil,

Tha Gauss law should be the simplest way.
The equation you mention for calculating the electric field is a consequence of the Gauss law.
You can use you "relevant equation" too. But in this case, try to simplify the calculations by using the symmetries of your problem.

Michel

 

[siddharth]

Tha Gauss law should be the simplest way.

The Gauss law won't be useful to obtain the field in this question.

Since the object is a hemisphere, you won't be able to construct a gaussian surface which will give you the Electric field. There's no symmetry.

 

[Swapnil]

Here is another answer I got using the principle of superposition:

\vec{E}(0,0,h) = \frac{Q}{4\pi\epsilon h^2}\Big(\frac{1}{\sqrt{a^2+h^2}} - \frac{1}{\sqrt{a^2-2ha+h^2}} + \frac{h}{a\sqrt{a^2-2ha+h^2}}\Big) \hat{z}

Is this closer to the right answer?

 

[Kurdt]

As siddharth mentioned it is best to obtain the potential at a point first. Do you know the general equation for obtaining the potential of a charge distribution?

 

[Swapnil]

Using the potential method, I got the following equation for the electric field:

\vec{E}(0,0,h) = \frac{Q}{4\pi\epsilon}\Big(\frac{h^3-a^3}{a^2 h^2\sqrt{h^2+a^2}}\Big)\hat{z}

How is this?

 

[lalbatros]

siddharth,

Sorry, I read you question too fast and did not notice the "hemi"-sphere.
I calculated the potential for the hemisphere.
For a charge density of 1, and droping the permitivity, this gives:

V = (r^2 + z^2)^(1/2)/(2*z) - (r^2 - 2*r*z + z^2)^(1/2)/(2*z)

and the electric field is a longer expression:

E = -(r - z)/(2*((r - z)^2)^(1/2)*z) - ((-r + z)^2)^(1/2)/(2*z^2)
- 1/(2*(r^2 + z^2)^(1/2)) + (r^2 + z^2)^(1/2)/(2*z^2)

In the above expressions, you have to be careful if taking the square roots.
I did the same exercice for the full-sphere, and I got also similar square roots. Taking the correct sign is essential to get the right solution inside or outside the sphere. For the hemisphere this should be also taken care for.

michel

 

[Swapnil]

BTW, would you get a different answer if you assume that h < a?

 

[Swapnil]

Here is another solution I got by redoing the surface integral for the E-field very carefully. This time I am pretty confident on my answer. I have checked and double checked everything.

\vec{E}(0,0,h) = \frac{Q}{4\pi\epsilon h^2}\Big(\frac{h-a}{|h-a|} + \frac{a}{\sqrt{a^2+h^2}} \Big)\hat{z}

Anyone else agrees?

 

[Swapnil]

I don't think that's right.

One way of doing this, is to first calculate V(z)[/tex], and then get \vec{E} from that.

<br /> <br /> Actually, that is the right answer (but only for h &gt; a). I got the same answer both by doing the surface integral of the E-field and also by calculating the potential and then finding the E-field.

 

[siddharth]

Actually, that is the right answer (but only for h > a)

I don't think so.

I got the same answer both by doing the surface integral of the E-field and also by calculating the potential and then finding the E-field.

Can you post the actual steps instead of just the final answer? It'll be easier to spot the errors that way, if any.

V = (r^2 + z^2)^(1/2)/(2*z) - (r^2 - 2*r*z + z^2)^(1/2)/(2*z)

The form of this solution seems correct.

 

[lalbatros]

siddharth,

The potential is given by the following integral (Mathematica notations):

Integrate[(2 Pi r Sin[teta]) / (4 Pi Sqrt[r^2Sin[teta]^2 + (r cos[teta]-z)^2]),{teta,0,Pi/2}]​

The numerator (2 Pi r Sin[teta])*d(teta) is the elementary charge from a circle on the sphere.
(for a charge density of 1, and dropping the eps factor)

Going to the http://integrals.wolfram.com/index.jsp" you can check the following result:
(but it is not very difficult with paper and pencil)

Integrate[(2 Pi r Sin[x]) / (4 Pi Sqrt[r^2 Sin[x]^2 + (r Cos[x]-z)^2]),x] = ...​

From this you can get the potential for a sphere or for an hemisphere:

VSphere = ...​

VHemisphere = ...​

You get the end result by introducing the correct charge density and eps.

Michel

 

[siddharth]

Michel, my question was directed to Swapnil, and I apologize for not stating that clearly.

Your method is indeed correct. But, could you modify the post so as to not give away the answer? I'd prefer it if Swapnil could figure it out on his/her own.

 

[lalbatros]

Sorry siddharth,

I saw many answers and I did not realize this was the homework section.
(I see the PF through an RSS feed and didn't look further than the question asked)
I hope anyway that I motivated Swapnil and others.
I really think that using Mathematica and other software can boost students,
even if paper and pencil is absolutely needed to get all the details.

Michel

 

[Swapnil]

I don't think so.

Can you post the actual steps instead of just the final answer? It'll be easier to spot the errors that way, if any.

**V = (r^2 + z^2)^(1/2)/(2*z) - (r^2 - 2*r*z + z^2)^(1/2)/(2*z)**
The form of this solution seems correct.

I think that michel and I are on the same track(and so are you siddharth since you agree with michel's answer). My answer is just simplified further than michel's. Notice that

\sqrt{r^2 - 2rz + z^2} = \sqrt{(r - z)^2} = |r - z|.

If you realize this fact and then take the derivative of the above expression for the potential with respect to z, you will get an answer which has the same form as the answer I originally mentioned. I would have posted all the steps myself but with my typing speed, it'll take me hours.

 

[Mr.4]

I tried the question under the assumption that you're only looking for the z component of the field and here's my answer:

<equation removed>

I just plain integrated it dividing it into small half rings if infinitesmal width.

 

[Swapnil]

I tried the question under the assumption that you're only looking for the z component of the field and here's my answer:

{2Q\ln({2 { [\frac{a }{h } + 1] }^{2} - 1)}/{4\pi\epsilon({a}^{2} + 2ah )}

I just plain integrated it dividing it into small half rings if infinitesmal width.

This doesn't make sense. How can the e-field be proportional to \frac{1}{h}\ln\Big(\frac{1}{h}}\Big)? It should more or less follow an inverse square realtionship.

 

[Mr.4]

Oops sorry I made a 'slight' mistake in the calculation. I've some doubts:

In what direction is the hemisphere oriented? Besides the mistake I made in my previous answer I also assumed that the hemisphere had its circular face lying on the XZ plane. That is the say that its axis of symmetry was assumed to be coincident with the y axis. Is this the case? Or is the hemisphere having its axis as Z-axis?

 

[Swapnil]

The axis of symmetry is the z-axis in my OP.

 

[Mr.4]

In that case the answer I'm getting is exactly same as your first answer. So I think that is indeed correct.

PS: I solved using the potential method wherein I calculated potential at any point on the z-axis and then differentiated it and put -ve sign.