### Div, Grad and Curl in Orthogonal Curvilinear Coordinates

The treatment here is standard, following that in Abraham and Becker, Classical Theory of Electricity and Magnetism.

Problems with a particular symmetry, such as cylindrical or spherical, are best attacked using coordinate systems that take full advantage of that symmetry.  For example, the Schrödinger equation for the hydrogen atom is best solved using spherical polar coordinates.  For this and other differential equation problems, then, we need to find the expressions for differential operators in terms of the appropriate coordinates.

We only look at orthogonal coordinate systems, so that locally the three axes (such as $r,\theta ,\phi$ ) are a mutually perpendicular set.

We denote the curvilinear coordinates by $\left({u}_{1},{u}_{2},{u}_{3}\right).$ The standard Cartesian coordinates for the same space are as usual $\left(x,y,z\right).$

Suppose now we take an infinitesimally small cube with edges parallel to the local curvilinear coordinate directions, and therefore with faces satisfying ${u}_{i}$ = constant, $i=1,2,3$ for the three pairs of faces.

The lengths of the edges are then ${h}_{1}d{u}_{1},{h}_{2}d{u}_{2}$ and ${h}_{3}d{u}_{3},$ where ${h}_{1},{h}_{2},{h}_{3}$ are in general functions of ${u}_{1},{u}_{2},{u}_{3}.$ That is to say, the distance across the cube from one corner to the opposite corner

$d{s}^{2}={h}_{1}^{2}d{u}_{1}^{2}+{h}_{2}^{2}d{u}_{2}^{2}+{h}_{3}^{2}d{u}_{3}^{2}=d{x}^{2}+d{y}^{2}+d{z}^{2}$

It is clear that the gradient of a function $\psi$ in the ${u}_{1}$ direction is

${\left(\nabla \psi \right)}_{1}=\underset{d{u}_{1}\to 0}{\mathrm{lim}}\frac{\psi \left(A\right)-\psi \left(0\right)}{{h}_{1}d{u}_{1}}=\frac{1}{{h}_{1}}\frac{\partial \psi }{\partial {u}_{1}}$.

The divergence of a vector field $\stackrel{\to }{V}$ in curvilinear coordinates is found using Gauss’ theorem, that the total vector flux through the six sides of the cube equals the divergence multiplied by the volume of the cube, in the limit of a small cube.

The area of the face bracketed by ${h}_{2}d{u}_{2}$ and ${h}_{3}d{u}_{3}$ is ${h}_{2}d{u}_{2}{h}_{3}d{u}_{3}.$ For that face, the component of the vector field contributing to the flow from the cube is $-{V}_{1},$ so the flow across the face is $-{V}_{1}{h}_{2}{h}_{3}d{u}_{2}d{u}_{3}.$ To find the flow across the opposite (parallel) face of the cube, corresponding to an increase in ${u}_{1}$ of $d{u}_{1},$ we must bear in mind that ${h}_{2},{h}_{3}$ and ${V}_{1}$ all vary with ${u}_{1},$ so the flow will be:

${V}_{1}{h}_{2}{h}_{3}d{u}_{2}d{u}_{3}+\frac{\partial }{\partial {u}_{1}}\left({h}_{2}{h}_{3}{V}_{1}\right)d{u}_{1}d{u}_{2}d{u}_{3}$

The first term here of course cancels the contribution from the other face. The remaining term, plus the terms with 123 replaced with 231 and 312 from the two other pairs of opposite faces, must, applying Gauss’ theorem, add to give .

This gives:

$\stackrel{\to }{\nabla }\cdot \stackrel{\to }{V}=\frac{1}{{h}_{1}{h}_{2}{h}_{3}}\left\{\frac{\partial }{\partial {u}_{1}}\left({h}_{2}{h}_{3}{V}_{1}\right)+\frac{\partial }{\partial {u}_{2}}\left({h}_{3}{h}_{1}{V}_{2}\right)+\frac{\partial }{\partial {u}_{3}}\left({h}_{1}{h}_{2}{V}_{3}\right)\right\}$.

Putting this together with the expression for the gradient gives immediately the expression for the Laplacian operator in curvilinear coordinates:

${\nabla }^{2}\psi =\frac{1}{{h}_{1}{h}_{2}{h}_{3}}\left\{\frac{\partial }{\partial {u}_{1}}\left(\frac{{h}_{2}{h}_{3}}{{h}_{1}}\frac{\partial \psi }{\partial {u}_{1}}\right)+\frac{\partial }{\partial {u}_{2}}\left(\frac{{h}_{3}{h}_{1}}{{h}_{2}}\frac{\partial \psi }{\partial {u}_{2}}\right)+\frac{\partial }{\partial {u}_{3}}\left(\frac{{h}_{1}{h}_{2}}{{h}_{3}}\frac{\partial \psi }{\partial {u}_{3}}\right)\right\}$

The curl of a vector field $\stackrel{\to }{A}$ is found by integrating around one of the square faces. Thus, the 1-component of $\stackrel{\to }{\nabla }×\stackrel{\to }{A}$ is given by integrating $\stackrel{\to }{A}\cdot \stackrel{\to }{ds}$ around the (23) square with two of its sides ${h}_{2}d{u}_{2}$ and ${h}_{3}d{u}_{3}.$ The integral must equal ${\left(\stackrel{\to }{\nabla }×\stackrel{\to }{A}\right)}_{1}$ multiplied by the area ${h}_{2}d{u}_{2}{h}_{3}d{u}_{3}.$ This gives

${\left(\stackrel{\to }{\nabla }×\stackrel{\to }{A}\right)}_{1}=\frac{1}{{h}_{2}{h}_{3}}\left\{\frac{\partial }{\partial {u}_{2}}\left({A}_{3}{h}_{3}\right)-\frac{\partial }{\partial {u}_{3}}\left({A}_{2}{h}_{2}\right)\right\}.$

Cylindrical Coordinates:

$\begin{array}{l}x=r\mathrm{cos}\phi \\ y=r\mathrm{sin}\phi \\ z=z\\ d{s}^{2}=d{r}^{2}+{r}^{2}d{\phi }^{2}+d{z}^{2}\end{array}$

Here $\left({u}_{1},{u}_{2},{u}_{3}\right)=\left(r,\phi ,z\right),$ and $\left({h}_{1},{h}_{2},{h}_{3}\right)=\left(1,r,1\right).$

Therefore, for example,

${\nabla }^{2}\psi =\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial \psi }{\partial r}\right)+\frac{1}{{r}^{2}}\frac{{\partial }^{2}\psi }{\partial {\phi }^{2}}+\frac{{\partial }^{2}\psi }{\partial {z}^{2}}.$

Spherical Polar Coordinates:

$\begin{array}{l}x=r\mathrm{sin}\theta \mathrm{cos}\phi \\ y=r\mathrm{sin}\theta \mathrm{sin}\phi \\ z=r\mathrm{cos}\theta \\ d{s}^{2}=d{r}^{2}+{r}^{2}d{\theta }^{2}+{r}^{2}{\mathrm{sin}}^{2}\theta d{\phi }^{2}\end{array}$

So $\left({u}_{1},{u}_{2},{u}_{3}\right)=\left(r,\theta ,\phi \right)$ and  $\left({h}_{1},{h}_{2},{h}_{3}\right)=\left(1,r,r\mathrm{sin}\theta \right).$

Here

${\nabla }^{2}\psi =\frac{1}{{r}^{2}}\frac{\partial }{\partial r}\left({r}^{2}\frac{\partial \psi }{\partial r}\right)+\frac{1}{{r}^{2}\mathrm{sin}\theta }\frac{\partial }{\partial \theta }\left(\mathrm{sin}\theta \frac{\partial \psi }{\partial \theta }\right)+\frac{1}{{r}^{2}{\mathrm{sin}}^{2}\theta }\frac{{\partial }^{2}\psi }{\partial {\phi }^{2}}.$