Register to reply 
Calculating equipotential surfaces 
Share this thread: 
#1
Jun2211, 02:35 AM

PF Gold
P: 7,120

So this may be a bit silly, but one thing I've never really learned in all my years is how one actually goes about calculating equipotential surfaces for arbitrary potentials? Let's say I have a potential that goes like
[tex]\Phi = A_0e^{\left({{r}\over{r_0}}\right)^2}\cos^2 (\phi) \sin^2(\theta)[/tex] where [itex] r, \phi[/itex] are your typical spherical coordinates and [itex]A_0[/itex] is a constant and [itex]r_0[/itex] is there for dimensional purposes and I want to determine even a single equipotential surface, how does one go about doing this? You typically see very simple examples in physics and you're typically finding something trivial like the shells you see around point particles, but this relies on plain ol intuition. How does one determine the surface with arbitrary (but wellbehaved) functions such as my example? Edit: I changed my function a few times so that it's Laplacian has no weird features for what I'm doing. 


#2
Jun2211, 04:35 AM

HW Helper
P: 3,464

I'd suggest doing the gradient in spherical coordinates on the potential. This would give the normal to the equipotential surface at all points in space. Then just specify a particular value for the potential to get one particular surface.



#3
Jun2211, 04:55 AM

PF Gold
P: 7,120

Hmm but that doesn't actually give me the surface as far as I can tell. I assume however one can define the surface, it'll be in terms of some parameter or relationship between the other coordinates.
My initial hunch would be you would set [itex]\Phi = \Phi_0[/itex] for generality, then determine the relationship [itex]\phi = \phi(r,\theta)[/itex] and my surface would be defined running through the [itex]0 \le \phi < 2\pi[/itex] and [itex]0 \le \theta < \pi[/itex]. I eventually have to figure out a surface integral which makes me wonder how do I even define a differential area element. 


#4
Jun2211, 05:26 AM

Sci Advisor
Thanks
P: 4,160

Calculating equipotential surfaces
The total differential of Φ is dΦ = ∂Φ/∂x dx + ∂Φ/∂y dy + ∂Φ/∂z dz. To get an equipotential surface, set this equal to zero and solve for dz:
dz =  (∂Φ/∂x dx + ∂Φ/∂y dy)/(∂Φ/∂z) This gives you a pair of PDEs to solve for the surface expressed as z(x,y): ∂z/∂x =  (∂Φ/∂x)/(∂Φ/∂z) ∂z/∂y =  (∂Φ/∂y)/(∂Φ/∂z) 


#5
Jun2211, 06:12 AM

HW Helper
P: 3,464




#6
Jun2211, 06:13 AM

PF Gold
P: 7,120

So what I really need to figure out is given a potential, how can I define equipotential surfaces for it so that I can later integrate the gradient of the potential over those equipotential surfaces. And to be clear, I would be integrating [itex]\nabla \Phi \cdot d\vec{A}[/itex] 


#7
Jun2211, 06:47 AM

HW Helper
P: 3,464

[tex] \nabla \Phi \cdot d \vec{A} = \nabla^2 \Phi dV [/tex]
But the Laplacian is zero for your problem , so luckily for you the integral is zero. If [itex] \Phi [/itex] was not Laplacian, then Bill_K gives the method for finding the change in z when you change x or y (while still keeping on the equipotential surface). So this effectively defines an equipotential surface when you start from a particular value. To use this, you'd need to redefine your problem in cartesian coordinates. Alternatively, you could do the total differential in spherical coordinates, which is more timeconsuming, but you could get r as dependent on polar angle and azimuthal angle, which is more convenient for your problem. This gives the 'parameter relationship' you were looking for. Then to integrate over this surface, you'd have to define the surface as being equal to the equipotential surface which you have just defined. I'm not sure how you'd do this. Maybe by creating a coordinate system where one variable is constant over the equipotential surface i.e. you'd need to define one of your variables as equal to [itex] \Phi [/itex], and then integrate over the other two variables? 


#8
Jun2211, 06:58 AM

Sci Advisor
PF Gold
P: 1,777

Take for example the function, [tex]\Phi(r,\theta,\phi) = r [/tex] This function has spherical shells of equipotential surfaces. Thus, [tex] \nabla\Phi = \hat{r}[/tex] [tex] \nabla\Phi\cdot d\mathbf{S} = \hat{r}\cdot\hat{r} r^2\sin\theta d\theta d\phi [/tex] Thus the value of the integral over the equipotential surface comes out to be 4 \pi r^2. 


#9
Jun2211, 07:49 AM

HW Helper
P: 3,464

Woops, sorry for some reason I thought that Pengwuino said in the OP that his potential [itex] \Phi [/itex] satisfied Laplace's equation. But he was just saying the Laplacian has no weird features. My mistake.



#10
Jun2211, 07:52 AM

Sci Advisor
PF Gold
P: 1,777




#11
Jun2211, 03:59 PM

PF Gold
P: 7,120




#12
Jun2311, 09:50 PM

HW Helper
P: 3,464

I think I agree with Pengwuino. Surely the Laplacian isn't necessarily zero anywhere for a general potential?
Back to the main problem: since the potential isn't spherical, I would rewrite it in terms of the cartesian coordinates. Also, the gradient of the potential will be perpendicular to the equipotential surface, so: [tex] \nabla \Phi \ \cdot \ d \vec{A} = \mid \nabla \Phi \mid \ dA [/tex] Now, to find what the equipotential surface is, we can use Bill K's method: [tex] \frac{ \partial z}{ \partial x} =  \frac{ \frac{ \partial \Phi }{ \partial x} }{ \frac{ \partial \Phi }{ \partial z } } [/tex] This gives us the partial derivative of z with respect to x. And we can use this to define the length integral which is perpendicular to y. ie [tex] ds = \sqrt{1 + { \frac{ \partial z}{ \partial x} }^2 } \ dz [/tex] So the integral becomes: [tex] \mid \nabla \Phi \mid \ \sqrt{1 + { \frac{ \partial z}{ \partial x} }^2 } \ dz \ dy [/tex] Someone please say if this is wrong, its been ages since I did this kind of maths, but I think its right Edit: After you've done these steps, I think you'll need to put x in terms of y and z before you do the integration, since x doesn't get integrated over? You can get the dependence of x on y and z by using Bill K's method. Basically, calculate the total differential of the potential along the surface, which we know is zero, and then calculate the total differential for x, then use these, along with boundary conditions, to get x(y,z). 


#13
Jun2311, 11:27 PM

Sci Advisor
PF Gold
P: 1,777




#14
Jun2411, 05:17 AM

HW Helper
P: 3,464




#15
Jun2711, 09:24 PM

PF Gold
P: 7,120

Ok so I've worked with a couple of potential functions and determined various surfaces. Now how does one form the surface element [itex]d\vec A[/itex]? Gah, I need to grab a good differential geometry text.
I have 2 things I need to do. 1) Determine the differential area element and eventually integrate this function over that surface and 2) determine the limits of integration for a corresponding volume integral. I wonder if Spivak would help in this? His intro text. 


#16
Jun2811, 12:07 AM

Sci Advisor
PF Gold
P: 1,777

For example, in cartesian coordinates then it is always something like [itex]\hat{z}dxdy[/itex]. If it is more complicated, then you should try to parameterize your equation for the surface. 


#17
Jun2811, 03:24 AM

PF Gold
P: 7,120

I've decided to use the potential [itex]A_0e^{{{r}\over{r_0}}^2}rcos(\theta)[/itex]
So far, I've determined the equipotential surfaces are given by [tex]x^2 + y^2 + {{z^2}\over{2}} = {{r_0^2ln(z)}\over{2}}+C[/tex] where the C is a constant that will denote the equipotential surface I believe. These surfaces are very complicated so I can't use any symmetry arguments to make things simpler. I wish I knew how to even look at the surface in Mathematica. 


Register to reply 
Related Discussions  
Equipotential surfaces  Introductory Physics Homework  0  
Equipotential surfaces  Introductory Physics Homework  6  
Equipotential Surfaces  Classical Physics  4  
Equipotential surfaces  Advanced Physics Homework  8 