On divergence
I think it's a bit hard to come from a cartesian understanding of divergence (and Laplacian operator for that matter) and then understand the same in terms of spherical polar coordinates intuitively?
I mean, say that a vector field is v = x i where i is the x-direction unit normal vector. This vector field is very easy to visualize. It's divergence is 1 because it's x value is increasing in the x direction.
Thus if v is a velocity field, the divergence value of 1 can mean one of two things:
Case 1: The fluid density is constant. The velocity vectors of the fluid entering the left-half of any closed surface are smaller in magnitude than the velocity vectors in the right-half, i.e. more fluid is entering any closed surface than is leaving it, which means that there must be a source of fluid everywhere.
Case 2: There are no fluid sources or sinks. Again the velocity vectors of the fluid entering in the left-half of any closed surface are smaller in magnitude less than the velocity vectors in the right-half, i.e. the fluid is diverging in the x-direction, it's speeding up making the fluid density less and less in the x-direction.
Now let's look at the vector field v = k*1/r^2 where k is just a constant and r is the distance from the origin. The formula for divergence in spherical polar coordinates gives a ZERO divergence everywhere except the origin, at which it is undefined. I'm sure that if you calculate the divergence of the field v = k*1/sqrt(x^2+y^2+z^2) in cartesian coordinates, you'll get zero as well, but if you look at any closed surface with the fluid entering as before, the velocity vectors entering the surface closest to the origin are larger in magnitude than the ones leaving the surface, therefore you should think that the density should be increasing / there should be a source everywhere?
I think this divergence intuition inherited from the cartesian coordinate system is confusing me when trying to appy it to spherically symmetrical situations and spherical polar coordinates.
On the Laplacian
Tac-Tics recently wrote: "The idea that it's just the sum of the second derivatives in each orthogonal direction just didn't sink in until I came back to it."
I kind of got the hang of thinking about it like this as well, but again when I apply it to the eletrostatic potential, phi = k*1/r, I just don't find it that intuitive anymore that the Laplacian of phi should give zero on the basis on my intuition of Laplacian from the cartesian world.
Take an arbitrary point along the x-axis (x=x0, y=0, z=0) in the electrostatic potental scalarfield of a positive point charge in the origin, then d^2(phi)/dx^2 is negative, but both d^2(phi)/dy^2 and d^2(phi)/dz^2 are also negative since all three increases distance to origin, all three which makes the value of the scalar potential lower. How can this give zero?
Again, you can probably show it's zero by during the proper calculation, but it shouldn't be too hard to point out the the fault in my logic and make it to make sense intuitively?
Thanks in advance to anyone who can help me sort this out :)