Cylindrical coordinate system

Coordinate system conversions

The cylindrical coordinate system is one of many three-dimensional coordinate systems. The following formulae may be used to convert between them.

Cartesian coordinates

For the conversion between cylindrical and Cartesian coordinates, it is convenient to assume that the reference plane of the former is the Cartesian xy-plane (with equation z = 0), and the cylindrical axis is the Cartesian z-axis. Then the z-coordinate is the same in both systems, and the correspondence between cylindrical (ρ, φ, z) and Cartesian (x, y, z) are the same as for polar coordinates, namely $${\displaystyle \begin{array}{cc}x& =\rho \cos \phi \\ y& =\rho \sin \phi \\ z& =z\end{array}}$$ in one direction, and $${\displaystyle \begin{array}{cc}\rho & =\sqrt{x^2+y^2}\\ \phi & =\{\begin{array}{ll}\text{indeterminate}& \text{if }x=0\text{ and }y=0\\ \arcsin \left(\frac{y}{\rho }\right)& \text{if }x\ge 0\\ -\arcsin \left(\frac{y}{\rho }\right)+\pi & \text{if }x<0\text{ and }y\ge 0\\ -\arcsin \left(\frac{y}{\rho }\right)-\pi & \text{if }x<0\text{ and }y<0\end{array}\end{array}}$$ in the other. The arcsine function is the inverse of the sine function, and is assumed to return an angle in the range [−⁠π/2⁠, +⁠π/2⁠] = [−90°, +90°]. These formulas yield an azimuth φ in the range [−180°, +180°].

By using the arctangent function that returns also an angle in the range [−⁠π/2⁠, +⁠π/2⁠] = [−90°, +90°], one may also compute ${\displaystyle \phi }$ without computing ${\displaystyle \rho }$ first $${\displaystyle \begin{array}{cc}\phi & =\{\begin{array}{ll}\text{indeterminate}& \text{if }x=0\text{ and }y=0\\ \frac{\pi }{2}\frac{y}{|y|}& \text{if }x=0\text{ and }y\ne 0\\ \arctan \left(\frac{y}{x}\right)& \text{if }x>0\\ \arctan \left(\frac{y}{x}\right)+\pi & \text{if }x<0\text{ and }y\ge 0\\ \arctan \left(\frac{y}{x}\right)-\pi & \text{if }x<0\text{ and }y<0\end{array}\end{array}}$$ For other formulas, see the article Polar coordinate system.

Many modern programming languages provide a function that will compute the correct azimuth φ, in the range (−π, π), given x and y, without the need to perform a case analysis as above. For example, this function is called by atan2(y, x) in the C programming language, and (atan y x) in Common Lisp.

Spherical coordinates

Spherical coordinates (radius r, elevation or inclination θ, azimuth φ), may be converted to or from cylindrical coordinates, depending on whether θ represents elevation or inclination, by the following:

Line and volume elements

In many problems involving cylindrical polar coordinates, it is useful to know the line and volume elements; these are used in integration to solve problems involving paths and volumes.

The line element is $${\displaystyle \mathrm{d}\mathit{r}=\mathrm{d}\rho \hat{\mathit{\rho }}+\rho \mathrm{d}\phi \hat{\mathit{\phi }}+\mathrm{d}z\hat{\mathit{z}}.}$$

The volume element is $${\displaystyle \mathrm{d}V=\rho \mathrm{d}\rho \mathrm{d}\phi \mathrm{d}z.}$$

The surface element in a surface of constant radius ρ (a vertical cylinder) is $${\displaystyle \mathrm{d}{S}_{\rho }=\rho \mathrm{d}\phi \mathrm{d}z.}$$

The surface element in a surface of constant azimuth φ (a vertical half-plane) is $${\displaystyle \mathrm{d}{S}_{\phi }=\mathrm{d}\rho \mathrm{d}z.}$$

The surface element in a surface of constant height z (a horizontal plane) is $${\displaystyle \mathrm{d}{S}_z=\rho \mathrm{d}\rho \mathrm{d}\phi .}$$

The del operator in this system leads to the following expressions for gradient, divergence, curl and Laplacian: $${\displaystyle \begin{array}{cc}\mathrm{\nabla }f& =\frac{\mathrm{\partial }f}{\mathrm{\partial }\rho }\hat{\mathit{\rho }}+\frac{1}{\rho }\frac{\mathrm{\partial }f}{\mathrm{\partial }\phi }\hat{\mathit{\phi }}+\frac{\mathrm{\partial }f}{\mathrm{\partial }z}\hat{\mathit{z}}\\ \mathrm{\nabla }\cdot \mathit{A}& =\frac{1}{\rho }\frac{\mathrm{\partial }}{\mathrm{\partial }\rho }\left(\rho A_{\rho }\right)+\frac{1}{\rho }\frac{\mathrm{\partial }{A}_{\phi }}{\mathrm{\partial }\phi }+\frac{\mathrm{\partial }{A}_z}{\mathrm{\partial }z}\\ \mathrm{\nabla }\times \mathit{A}& =\left(\frac{1}{\rho }\frac{\mathrm{\partial }{A}_z}{\mathrm{\partial }\phi }-\frac{\mathrm{\partial }{A}_{\phi }}{\mathrm{\partial }z}\right)\hat{\mathit{\rho }}+\left(\frac{\mathrm{\partial }{A}_{\rho }}{\mathrm{\partial }z}-\frac{\mathrm{\partial }{A}_z}{\mathrm{\partial }\rho }\right)\hat{\mathit{\phi }}+\frac{1}{\rho }\left(\frac{\mathrm{\partial }}{\mathrm{\partial }\rho }\left(\rho A_{\phi }\right)-\frac{\mathrm{\partial }{A}_{\rho }}{\mathrm{\partial }\phi }\right)\hat{\mathit{z}}\\ {\mathrm{\nabla }}^{2}f& =\frac{1}{\rho }\frac{\mathrm{\partial }}{\mathrm{\partial }\rho }\left(\rho \frac{\mathrm{\partial }f}{\mathrm{\partial }\rho }\right)+\frac{1}{{\rho }^2}\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{\phi }^2}+\frac{{\mathrm{\partial }}^{2}f}{\mathrm{\partial }{z}^2}\end{array}}$$

Cylindrical harmonics

The solutions to the Laplace equation in a system with cylindrical symmetry are called cylindrical harmonics.

Kinematics

In a cylindrical coordinate system, the position of a particle can be written as $${\displaystyle \mathit{r}=\rho \hat{\mathit{\rho }}+z\hat{\mathit{z}}.}$$ The velocity of the particle is the time derivative of its position, $${\displaystyle \mathit{v}=\frac{\mathrm{d}\mathit{r}}{\mathrm{d}t}=\dot{\rho }\hat{\mathit{\rho }}+\rho \dot{\phi }\hat{\mathit{\phi }}+\dot{z}\hat{\mathit{z}},}$$ where the term ${\displaystyle \rho \dot{\phi }\hat{\phi }}$ comes from the Poisson formula ${\displaystyle \frac{\mathrm{d}\hat{\rho }}{\mathrm{d}t}=\dot{\phi }\hat{z}\times \hat{\rho }}$. Its acceleration is $${\displaystyle \mathit{a}=\frac{\mathrm{d}\mathit{v}}{\mathrm{d}t}=\left(\ddot{\rho }-\rho {\dot{\phi }}^2\right)\hat{\mathit{\rho }}+\left(2\dot{\rho }\dot{\phi }+\rho \ddot{\phi }\right)\hat{\mathit{\phi }}+\ddot{z}\hat{\mathit{z}}}$$