Skip to main content

An Intuitive Way to the Spherical Gradient and Laplacian

It's that time of year again when physics students everywhere are deriving the spherical and cylindrical del, nabla, gradient, or Laplacian operators.  Every derivation I saw prior to this week involved lots of algebra and the chain rule... even mine.  Fortunately for me, a comment on my derivation, and a homework assignment from Rutgers [pdf] led me to a far simpler and more intuitive way of doing things.  You just start from the differential displacement in a given coordinate system and go from there.

The differential displacement in spherical coordinates is:



The element in the r direction is easy to understand.  A small displacement along the r direction is represented as dr.  The theta and phi displacements might not be as obvious.  The graphic to the left illustrates what's going on.  With small displacement along the theta direction you're moving along a circle with radius r.  The distance you've moved is equal to the length of the arc which is equal to the radius of the circle, r times the small angle, hence, .  The element in the in the phi direction is similar to the element in the theta direction in that you're moving along a circle.  However, this time you have to keep in mind that the radius of the circle is not r.  The r vector is propped up by angle of theta and so, the projection of r on the r/phi plane is , and the distance traveled along the phi displacement is .
Once we have our differential displacement element, the rest is kind of easy.  First, we get the gradient operator.  The gradient is just the derivative of a function with respect to the differential displacement, so we get:


The next step is to derive the divergence in spherical coordinates.  First, we need the differential surface area element in all three dimensions.  I don't have a great picture for this yet, so you'll just have to visualize a small square on the surface of a sphere.  The surface area element in the r direction is a rectangle on the surface of a sphere with a length of  on one side and a length of  on the other.  So, the element of area is the product of the sides or  .  In the theta direction, the sides are a small displacement in the r direction and a small displacement in the phi direction giving us .  Finally, in the phi direction the sides of the rectangle are a small displacement along r and a small displacement along theta giving us , (.

You might be wondering what all the surface elements had to do with divergence.  This is where the tip I mentioned at the start of the post comes in.  Divergence is defined to be the flux of a vector field, (the vector field dot producted with the surface area), per unit volume.  Now that we have a differential surface area vector, we can dot product our vector of interest, (say A), with it to get the flux through a differential surface area.  We wind up with:

Now, we need the differential element of volume for spherical coordinates.  Once again, we can use purely geometric considerations to get where we want to go.  Visualize making a small cube out of the small rectangles we defined above for the surface area differentials.  If you took the volume of that cube, you'd just multiply the three sides:

Now to get the divergence, we just divide each term of the flux by the differential volume element.  For the r term we get:

the theta and phi differentials cancel out.  You can bring the sin theta term outside of the derivative and cancel it since it doesn't depend on r.

For the theta term we get:

The r and phi differentials cancel.  The r can be taken outside the derivative since it is constant with respect to the remaining theta differential.

Finally, for the phi term we get:

The r and theta differentials cancel and we can bring the r outside the remaining phi derivative.

All together, we now have the three terms of the divergence in spherical coordinates:

We didn't use the chain rule once, just simple geometry!

The final task is to derive the Laplacian.  The Laplacian is the divergence of the gradient.  We simply substitute our gradient result in place of the vector A in our divergence experession to get the following three terms:
r:

theta:

phi:

 

Comments

  1. Dear Blogger, first of thanks for this post, it got me started in the right direction.

    But, unfortunately your way of deriving had two connected oversimplifications which I think to have solved:

    First of, you need to account for the flux going through the box. To do so you dot the surface of ONE side with the vector. But here you do not account for the other side. A constant vector in the dot product with both sides would yield no net flux...

    Solution: You actually skip a step where you should say that the difference in the vector A in the direction of the surface on opposite sides of the box is the infinitesimal difference dA between A(r) and A(r + dr) for instance.

    The other problem was to me a manifestation of the previous simplification:

    In calculating the end result you do the following (1/dr) x r^2 x Ar = (d/dr) x r^2 x Ar. This is a big oversimplication/ error to make.

    Solution: Use the dA instead from before this way you`ll get 1/dr x r^2 x dA = d(r^2 x A)/dr.

    I would really appreciate if you responded.

    ReplyDelete
  2. Willem,
    Thanks very much for all the info. I think you're correct and it clears up a bit of confusion I had regarding some of the steps. Thanks very much! I'm going to write more when I have time to review what you've said more completely. Thanks again!

    ReplyDelete
  3. Sorry >__< when i try to use the equation for Flux/volume

    i didnt get what you get in the equation

    for example: (Ar r^2 sin(theta) d(theta) dphi)/(r^2 sin(pheta) dr d(theta) d(phi))

    =Ar/(dr) ~~ what happen????(not the above one of d(Ar r^2)/dr * (1/r^2))

    thank you so much

    ReplyDelete

Post a Comment

Please leave your comments on this topic:

Popular posts from this blog

More Cowbell! Record Production using Google Forms and Charts

First, the what : This article shows how to embed a new Google Form into any web page. To demonstrate ths, a chart and form that allow blog readers to control the recording levels of each instrument in Blue Oyster Cult's "(Don't Fear) The Reaper" is used. HTML code from the Google version of the form included on this page is shown and the parts that need to be modified are highlighted. Next, the why : Google recently released an e-mail form feature that allows users of Google Documents to create an e-mail a form that automatically places each user's input into an associated spreadsheet. As it turns out, with a little bit of work, the forms that are created by Google Docs can be embedded into any web page. Now, The Goods: Click on the instrument you want turned up, click the submit button and then refresh the page. Through the magic of Google Forms as soon as you click on submit and refresh this web page, the data chart will update immediately. Turn up the:

Cool Math Tricks: Deriving the Divergence, (Del or Nabla) into New (Cylindrical) Coordinate Systems

Now available as a Kindle ebook for 99 cents ! Get a spiffy ebook, and fund more physics The following is a pretty lengthy procedure, but converting the divergence, (nabla, del) operator between coordinate systems comes up pretty often. While there are tables for converting between common coordinate systems , there seem to be fewer explanations of the procedure for deriving the conversion, so here goes! What do we actually want? To convert the Cartesian nabla to the nabla for another coordinate system, say… cylindrical coordinates. What we’ll need: 1. The Cartesian Nabla: 2. A set of equations relating the Cartesian coordinates to cylindrical coordinates: 3. A set of equations relating the Cartesian basis vectors to the basis vectors of the new coordinate system: How to do it: Use the chain rule for differentiation to convert the derivatives with respect to the Cartesian variables to derivatives with respect to the cylindrical variables. The chain

The Valentine's Day Magnetic Monopole

There's an assymetry to the form of the two Maxwell's equations shown in picture 1.  While the divergence of the electric field is proportional to the electric charge density at a given point, the divergence of the magnetic field is equal to zero.  This is typically explained in the following way.  While we know that electrons, the fundamental electric charge carriers exist, evidence seems to indicate that magnetic monopoles, the particles that would carry magnetic 'charge', either don't exist, or, the energies required to create them are so high that they are exceedingly rare.  That doesn't stop us from looking for them though! Keeping with the theme of Fairbank[1] and his academic progeny over the semester break, today's post is about the discovery of a magnetic monopole candidate event by one of the Fairbank's graduate students, Blas Cabrera[2].  Cabrera was utilizing a loop type of magnetic monopole detector.  Its operation is in concept very sim