Div an curl in different coordinate systems

[Sunshine]

To calculate the divergence of a vectorfield in cartesian coordinates, you can think of it as a dot product, and to calculate the curl, you can think of it as a cross product. But how can you calculate the div and curl when you have spherical or cylindrical coordinates, without explicitely converting them to the cartesian system?

For example the radial unity vector e_r = (sin a * cos b, sin a * sin b, cos a)
(standard: a = Theta, b = phi) should affect the curl and div by more than the scaling factor (h)

P.S Is it possible to use tex when writing a post?

 

[dextercioby]

Yes,we have LaTex enabled.

http://mathworld.wolfram.com/Curl.html

http://mathworld.wolfram.com/Divergence.html

Daniel.

 

[quasar987]

I think a good starting point is to define the nabla operator as the operator \nabla such that

df=\nabla f\cdot d\vec{r}

This allows us to recover easily the form of \nabla in any coordinate system.

 

[quasar987]

For exemple, in cylindrical coordinate: Consider a function f = f(r,\theta,z).Then

df = \frac{\partial f}{\partial r}dr + \frac{\partial f}{\partial \theta}d\theta + \frac{\partial f}{\partial z }dz

and

d\vec{r} = \hat{r}dr + \hat{\theta}rd\theta + \hat{z}dz

So after messing around a bit you find that

df=\nabla f\cdot d\vec{r} \Leftrightarrow \nabla = \hat{r}\frac{\partial}{\partial r} + \frac{\hat{\theta}}{r}\frac{\partial}{\partial \theta}+ \hat{z}\frac{\partial}{\partial z }

 

[Sunshine]

d\vec{r} = \hat{r}dr + \hat{\theta}rd\theta + \hat{z}dz

So after messing around a bit you find that

df=\nabla f\cdot d\vec{r} \Leftrightarrow \nabla = \hat{r}\frac{\partial}{\partial r} + \frac{\hat{\theta}}{r}\frac{\partial}{\partial \theta}+ \hat{z}\frac{\partial}{\partial z }

Hm, shouldn't it be
d\vec{r} = \hat{r}dr + \mathbf{\dfrac{\hat{\theta}}{r}}d\theta + \hat{z}dz
if you divide with the scalefactor?

Otherwise the second equation would become
df=\nabla f\cdot d\vec{r} \Leftrightarrow \nabla = \hat{r}\frac{\partial}{\partial r} + \hat{\theta}r\frac{\partial}{\partial \theta}+ \hat{z}\frac{\partial}{\partial z }

I'm still not sure if I've got it right, but I'll think over this, and return with more questions if they arise. Thanks for the help :)

 

[quasar987]

Hm, shouldn't it be
d\vec{r} = \hat{r}dr + \mathbf{\dfrac{\hat{\theta}}{r}}d\theta + \hat{z}dz
if you divide with the scalefactor?

Mmh.. the position vector in cylindrical coordinate is just \vec{r}=r\hat{r}+z\hat{z}. If you differentiate that with respect to an hypothetical variable t, you get

\frac{d\vec{r}}{dt} = \frac{dr}{dt}\hat{r}+r\frac{d\hat{r}}{dt}+\frac{dz}{dt}\hat{z}+\z\frac{d\hat{z}}{dt} = \frac{dr}{dt}\hat{r}+r\frac{d\hat{r}}{d\theta}\frac{d\theta}{dt}+\frac{dz}{dt}\hat{z}+0 = \frac{dr}{dt}\hat{r}+r\frac{d\theta}{dt}\hat{\theta}+\frac{dz}{dt}\hat{z}

And multiplying by dt gives d\vec{r}.

Otherwise the second equation would become
df=\nabla f\cdot d\vec{r} \Leftrightarrow \nabla = \hat{r}\frac{\partial}{\partial r} + \hat{\theta}r\frac{\partial}{\partial \theta}+ \hat{z}\frac{\partial}{\partial z }

Better double check your calculations.

 

[Sunshine]

You're right, my mistake

However my biggest problem remains.

Let's look at the position vector
\vec{r} = (r cos \theta, r sin \theta, z)

\frac{\partial \vec{r}}{\partial r} = (cos \theta, sin \theta, 0)

h_r = \sqrt{cos^2 \theta + sin^2 \theta} = 1

and according to the definition

\hat{r} = \frac{1}{h_r} \frac{\partial \vec{r}}{\partial r} = (cos \theta, sin \theta, 0)

the question is, which of the equations below are correct

\nabla \cdot \hat{r} = \nabla \cdot (1,0,0)
OR
\nabla \cdot \hat{r} = \nabla \cdot (cos \theta, sin \theta, 0)

...if I use the \nabla that you gave for cylindrical coordinates.
If it's #1, would #2 apply if the \nabla was defined on another way, i.e not as

\nabla = \hat{r}\frac{\partial}{\partial r} + \hat{\theta}r\frac{\partial}{\partial \theta}+ \hat{z}\frac{\partial}{\partial z }

 

[quasar987]

Both

\nabla \cdot \hat{r} = \nabla \cdot (1,0,0)
AND
\nabla \cdot \hat{r} = \nabla \cdot (cos \theta, sin \theta, 0)

are correct, if you notice the distinction that in #1, (1,0,0) means (r=1,\theta=0,z=0) (i.e. (1,0,0) represents the vector \hat{r} = 1\hat{r}+0\hat{\theta}+0\hat{z}), while in #2, (cos\theta,sin\theta,0) means (x=cos\theta,y=sin\theta,z=0) (i.e. represents the vector \hat{r}=cos\theta \hat{x} + sin\theta \hat{y} + 0\hat{z}

The base of \mathbb{R}^3 used in #1 is different than that used in #2. Therfor if you're going to take the dot product of either of them with nabla, you got to express nabla in the same base.

 

[Sunshine]

so you mean that
(\frac{\partial}{\partial x}, \frac{\partial}{\partial y},\frac{\partial}{\partial z }) \cdot (cos \theta, sin \theta, 0) = (\frac{\partial}{\partial r}, \frac{1}{r}\frac{\partial}{\partial \theta},\frac{\partial}{\partial z }) \cdot (1,0,0)

Well, I see that it's true in this case where both sides are zero but is it so for an arbitrary vector field?

And is curl similar?

 

[dextercioby]

Of course.Remember that the divergence of a vector field is a scalar wrt general coordinates transformations,so eveluating it in cartesian is no difference than eveluating it in cylindrical or spherical...

As for the curl,well,it gives a pseudotensor of a higher rank,so it's not invariant under general coordinate transformations...

Daniel.

 

[quasar987]

Look at equations (7), (8) at http://mathworld.wolfram.com/Divergence.html

and (7), (8), (9), at http://mathworld.wolfram.com/Curl.html

 

[Sunshine]

I've turned those formulas upside down by now, and still very confused...

The thing is that I'm supposed to calculate \nabla \cdot \vec{e_x} and \nabla \times \vec{e_x} where \vec{e_x} are the unit vectors with x = r, \rho, \theta, \phi (cylindrical and spherical coordinates).
It means a total of 8 "questions" and all of them (!) becomes 0 if I use any of the equations in my last post (you always get 0 when differentiating a scalar)... that's why I've been suspecting that I'm supposed to use
\vec{e_r} = (cos \theta, sin \theta, 0)

\vec{e_{\theta}} = (-sin \theta, cos \theta, 0)

etc. and differentiate with respect to r, \theta, \phi ...

But that isn't the case...?