Dynamics/Kinematics/Coordinate Systems/Cylindrical

Introduction
Content from Cylindrical coordinate system



A cylindrical coordinate system is a three-dimensional coordinate system that specifies point positions by the distance from a chosen reference axis, the direction from the axis relative to a chosen reference direction, and the distance from a chosen reference plane perpendicular to the axis. The latter distance is given as a positive or negative number depending on which side of the reference plane faces the point.

The origin of the system is the point where all three coordinates can be given as zero. This is the intersection between the reference plane and the axis. The axis is variously called the cylindrical or longitudinal axis, to differentiate it from the polar axis, which is the ray that lies in the reference plane, starting at the origin and pointing in the reference direction. Other directions perpendicular to the longitudinal axis are called radial lines.

The distance from the axis may be called the radial distance or radius, while the angular coordinate is sometimes referred to as the angular position or as the azimuth. The radius and the azimuth are together called the polar coordinates, as they correspond to a two-dimensional polar coordinate system in the plane through the point, parallel to the reference plane. The third coordinate may be called the height or altitude (if the reference plane is considered horizontal), longitudinal position, or axial position.

They are sometimes called "cylindrical polar coordinates" and "polar cylindrical coordinates", and are sometimes used to specify the position of stars in a galaxy ("galactocentric cylindrical polar coordinates").

Definition
The three coordinates ($O$, $A$, $L$) of a point $r$ are defined as:
 * The axial distance or radial distance $θ$ is the Euclidean distance from the $z$-axis to the point $ρ$.
 * The azimuth $φ$ is the angle between the reference direction on the chosen plane and the line from the origin to the projection of $z$ on the plane.
 * The axial coordinate or height $P$ is the signed distance from the chosen plane to the point $ρ$.

Unique cylindrical coordinates
As in polar coordinates, the same point with cylindrical coordinates $ρ = 4$ has infinitely many equivalent coordinates, namely $φ = 130°$ and $z = 4$ where $z$ is any integer. Moreover, if the radius $P$ is zero, the azimuth is arbitrary.

In situations where someone wants a unique set of coordinates for each point, one may restrict the radius to be non-negative ($(r,θ,z)$) and the azimuth $φ$ to lie in a specific interval spanning 360°, such as $Z = v_{bz}t$ or $(ρ, φ, z)$.

Conventions
The notation for cylindrical coordinates is not uniform. The ISO standard 31-11 recommends $(ρ, φ ± n×360°, z)$, where $P$ is the radial coordinate, $z$ the azimuth, and $P$ the height. However, the radius is also often denoted $n$ or $ρ$, the azimuth by $φ$ or $ρ$, and the third coordinate by $φ$ or (if the cylindrical axis is considered horizontal) $z$, or any context-specific letter.





In concrete situations, and in many mathematical illustrations, a positive angular coordinate is measured counterclockwise as seen from any point with positive height.

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 $r$-plane (with equation $(−ρ, φ ± (2n + 1)×180°, z),$), and the cylindrical axis is the Cartesian $s$-axis. Then the $θ$-coordinate is the same in both systems, and the correspondence between cylindrical $ρ ≥ 0$ and Cartesian $[−180°,+180°]$ are the same as for polar coordinates, namely

\begin{align} x &= \rho \cos \varphi \\ y &= \rho \sin \varphi \\ z &= z \end{align} $$ in one direction, and
 * $$\begin{align} \rho &= \sqrt{x^2+y^2} \\

\varphi &= \begin{cases} 0 & \mbox{if } x = 0 \mbox{ and } y = 0\\ \arcsin\left(\frac{y}{\rho}\right) & \mbox{if } x \geq 0 \\ \arctan\left(\frac{y}{x}\right) & \mbox{if } x > 0 \\ -\arcsin\left(\frac{y}{\rho}\right) + \pi & \mbox{if } x < 0 \end{cases} \end{align}$$ in the other. The arcsin function is the inverse of the sine function, and is assumed to return an angle in the range $[0,360°]$ = $(ρ, φ, z)$. These formulas yield an azimuth $t$ in the range $(ρ, φ, z)$. For other formulas, see the polar coordinate article.

Many modern programming languages provide a function that will compute the correct azimuth $h$, in the range $ρ = 2$, 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 $x$, elevation or inclination $z$, azimuth $x$), may be converted into cylindrical coordinates by:

Cylindrical coordinates may be converted into spherical coordinates by:

Position
Material taken from Vector fields in cylindrical and spherical coordinates

Vectors are defined in cylindrical coordinates by (ρ, $$\phi$$, z), where
 * ρ is the length of the vector projected onto the xy-plane,
 * $$\phi$$ is the angle between the projection of the vector onto the xy-plane (i.e. ρ) and the positive x-axis (0 ≤ $$\phi$$ < 2π),
 * z is the regular z-coordinate.

(ρ, $$\phi$$, z) is given in cartesian coordinates by:


 * $$\begin{bmatrix} \rho \\ \phi \\ z \end{bmatrix} =

\begin{bmatrix} \sqrt{x^2 + y^2} \\ \operatorname{arctan}(y / x) \\ z \end{bmatrix},\ \ \ 0 \le \phi < 2\pi, $$

or inversely by:


 * $$\begin{bmatrix} x \\ y \\ z \end{bmatrix} =

\begin{bmatrix} \rho\cos\phi \\ \rho\sin\phi \\ z \end{bmatrix}.$$

Any vector field can be written in terms of the unit vectors as:
 * $$\mathbf A = A_x \mathbf{\hat x} + A_y \mathbf{\hat y} + A_z \mathbf{\hat z}

= A_\rho \mathbf{\hat \rho} + A_\phi \boldsymbol{\hat \phi} + A_z \mathbf{\hat z}$$ The cylindrical unit vectors are related to the Cartesian unit vectors by:
 * $$\begin{bmatrix}\boldsymbol{\hat \rho} \\ \boldsymbol{\hat\phi} \\ \mathbf{\hat z}\end{bmatrix}

= \begin{bmatrix} \cos\phi & \sin\phi & 0 \\ -\sin\phi & \cos\phi & 0 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} \mathbf{\hat x} \\ \mathbf{\hat y} \\ \mathbf{\hat z} \end{bmatrix}$$

Note: the matrix is an orthogonal matrix, that is, its inverse is simply its transpose.

An interesting but confusing point is the abstraction of $$\boldsymbol{\hat \phi}$$, which can be thought of as being an orthogonal vector to $$\boldsymbol \hat \rho$$:


 * $$\boldsymbol{\hat \phi}=-\sin(\phi)\mathsf{I}+\cos(\phi)\mathsf{J} $$

It is also possible to represent a position vector in Cartesian and cylindrical coordinates as follows:


 * $$\mathsf{r}_P=X_P\mathsf{I}+Y_P\mathsf{J}+Z_P\mathsf{K}=\rho \boldsymbol{\hat \rho}+Z_P\mathsf{K}$$

where we can leave out the term $$\phi \boldsymbol{\hat \phi}$$ and describe the position fully.

Velocity
Material taken from Vector fields in cylindrical and spherical coordinates

To find out how the vector field A changes in time we calculate the time derivatives. For this purpose we use Newton's notation for the time derivative ($$\dot{\mathbf{A}}$$). In Cartesian coordinates this is simply:
 * $$\dot{\mathbf{A}} = \dot{A}_x \hat{\mathbf{x}} + \dot{A}_y \hat{\mathbf{y}} + \dot{A}_z \hat{\mathbf{z}}$$

However, in cylindrical coordinates this becomes:
 * $$\dot{\mathbf{A}} = \dot{A}_\rho \hat{\boldsymbol{\rho}} + A_\rho \dot{\hat{\boldsymbol{\rho}}}

+ \dot{A}_\phi \hat{\boldsymbol{\phi}} + A_\phi \dot{\hat{\boldsymbol{\phi}}} + \dot{A}_z \hat{\boldsymbol{z}} + A_z \dot{\hat{\boldsymbol{z}}}$$

We need the time derivatives of the unit vectors. They are given by:
 * $$\begin{align}

\dot{\hat{\mathbf{\rho}}} & = \dot{\phi} \hat{\boldsymbol{\phi}} \\ \dot{\hat{\boldsymbol{\phi}}} & = - \dot\phi \hat{\mathbf{\rho}} \\ \dot{\hat{\mathbf{z}}} & = 0 \end{align}$$

So the time derivative simplifies to:
 * $$\dot{\mathbf{A}} = \hat{\boldsymbol{\rho}} (\dot{A}_\rho - A_\phi \dot{\phi})

+ \hat{\boldsymbol{\phi}} (\dot{A}_\phi + A_\rho \dot{\phi}) + \hat{\mathbf{z}} \dot{A}_z$$

It is also possible to represent position and velocity vectors in Cartesian and cylindrical coordinates as follows:


 * $$\mathsf{r}_P=X_P\mathsf{I}+Y_P\mathsf{J}+Z_P\mathsf{K}=\rho \boldsymbol{\hat \rho}+Z_P\mathsf{K}$$

and


 * $$\mathsf{v}_P=\dot X_P\mathsf{I}+\dot Y_P\mathsf{J}+\dot Z_P\mathsf{K}=\dot \rho \boldsymbol{\hat \rho}+\rho \dot \phi \boldsymbol{\hat \phi}+\dot Z_P\mathsf{K}$$

Acceleration
Material taken from Vector fields in cylindrical and spherical coordinates

The second time derivative is of interest in physics, as it is found in equations of motion for classical mechanical systems. The second time derivative of a vector field in cylindrical coordinates is given by:
 * $$\mathbf{\ddot A} = \mathbf{\hat \rho} (\ddot A_\rho - A_\phi \ddot\phi - 2 \dot A_\phi \dot\phi - A_\rho \dot\phi^2)

+ \boldsymbol{\hat\phi} (\ddot A_\phi + A_\rho \ddot\phi + 2 \dot A_\rho \dot\phi - A_\phi \dot\phi^2) + \mathbf{\hat z} \ddot A_z$$

To understand this expression, we substitute A = P, where p is the vector (\rho, θ, z).

This means that $$\mathbf{A} = \mathbf{P} = \rho \mathbf{\hat \rho} + z \mathbf{\hat z}$$.

After substituting we get:
 * $$\ddot\mathbf{P} = \mathbf{\hat \rho} (\ddot \rho - \rho \dot\phi^2)

+ \boldsymbol{\hat\phi} (\rho \ddot\phi + 2 \dot \rho \dot\phi) + \mathbf{\hat z} \ddot z$$

In mechanics, the terms of this expression are called:
 * $$\begin{align}

\ddot \rho \mathbf{\hat \rho}         &= \mbox{central outward acceleration} \\ -\rho \dot\phi^2 \mathbf{\hat \rho}   &= \mbox{centripetal acceleration} \\ \rho \ddot\phi \boldsymbol{\hat\phi}     &= \mbox{angular acceleration} \\ 2 \dot \rho \dot\phi \boldsymbol{\hat\phi} &= \mbox{Coriolis effect} \\ \ddot z \mathbf{\hat z}              &= \mbox{z-acceleration}

\end{align}$$

Another form for describing the acceleration of position is given as follows:


 * $$\mathsf{a}_P=\ddot X_P\mathsf{I}+\ddot Y_P\mathsf{J}+\ddot Z_P\mathsf{K}=(\ddot \rho -\rho \dot \phi ^2) \boldsymbol{\hat \rho}+ (\rho \ddot \phi + 2 \dot \rho \dot \phi) \boldsymbol{\hat \phi}+\ddot Z_P\mathsf{K}$$

See also Centripetal force, Angular acceleration, and Coriolis effect.