Toroidal and poloidal coordinates

From HandWiki
Short description: Coordinate system relative to a torus
A diagram depicting the poloidal (θ) direction, represented by the red arrow, and the toroidal (ζ or φ) direction, represented by the blue arrow.

The terms toroidal and poloidal refer to directions relative to a torus of reference. They describe a three-dimensional coordinate system in which the poloidal direction follows a small circular ring around the surface, while the toroidal direction follows a large circular ring around the torus, encircling the central void.

The earliest use of these terms cited by the Oxford English Dictionary is by Walter M. Elsasser (1946) in the context of the generation of the Earth's magnetic field by currents in the core, with "toroidal" being parallel to lines of latitude and "poloidal" being in the direction of the magnetic field (i.e. towards the poles).

The OED also records the later usage of these terms in the context of toroidally confined plasmas, as encountered in magnetic confinement fusion. In the plasma context, the toroidal direction is the long way around the torus, the corresponding coordinate being denoted by z in the slab approximation or ζ or φ in magnetic coordinates; the poloidal direction is the short way around the torus, the corresponding coordinate being denoted by y in the slab approximation or θ in magnetic coordinates. (The third direction, normal to the magnetic surfaces, is often called the "radial direction", denoted by x in the slab approximation and variously ψ, χ, r, ρ, or s in magnetic coordinates.)

Example

As a simple example from the physics of magnetically confined plasmas, consider an axisymmetric system with circular, concentric magnetic flux surfaces of radius [math]\displaystyle{ r }[/math] (a crude approximation to the magnetic field geometry in an early tokamak but topologically equivalent to any toroidal magnetic confinement system with nested flux surfaces) and denote the toroidal angle by [math]\displaystyle{ \zeta }[/math] and the poloidal angle by [math]\displaystyle{ \theta }[/math]. Then the toroidal/poloidal coordinate system relates to standard Cartesian coordinates by these transformation rules:

[math]\displaystyle{ x = (R_0 +r \cos \theta) \cos\zeta }[/math]
[math]\displaystyle{ y = s_\zeta (R_0 + r \cos \theta) \sin\zeta }[/math]
[math]\displaystyle{ z = s_\theta r \sin \theta. }[/math]

where [math]\displaystyle{ s_\theta = \pm 1, s_\zeta = \pm 1 }[/math].

The natural choice geometrically is to take [math]\displaystyle{ s_\theta = s_\zeta = +1 }[/math], giving the toroidal and poloidal directions shown by the arrows in the figure above, but this makes [math]\displaystyle{ r,\theta,\zeta }[/math] a left-handed curvilinear coordinate system. As it is usually assumed in setting up flux coordinates for describing magnetically confined plasmas that the set [math]\displaystyle{ r,\theta,\zeta }[/math] forms a right-handed coordinate system, [math]\displaystyle{ \nabla r\cdot\nabla\theta\times\nabla\zeta \gt 0 }[/math], we must either reverse the poloidal direction by taking [math]\displaystyle{ s_\theta = -1, s_\zeta = +1 }[/math], or reverse the toroidal direction by taking [math]\displaystyle{ s_\theta = +1, s_\zeta = -1 }[/math]. Both choices are used in the literature.

Kinematics

To study single particle motion in toroidally confined plasma devices, velocity and acceleration vectors must be known. Considering the natural choice [math]\displaystyle{ s_\theta = s_\zeta = +1 }[/math], the unit vectors of toroidal and poloidal coordinates system [math]\displaystyle{ \left(r,\theta,\zeta\right) }[/math] can be expressed as:

[math]\displaystyle{ \mathbf{e}_r = \begin{pmatrix} \cos\theta \cos\zeta \\ \cos\theta \sin\zeta \\ \sin\theta \end{pmatrix} \quad \mathbf{e}_\theta = \begin{pmatrix} -\sin\theta \cos\zeta \\ -\sin\theta \sin\zeta \\ \cos\theta \end{pmatrix} \quad \mathbf{e}_\zeta = \begin{pmatrix} -\sin\zeta \\ \cos\zeta \\ 0 \end{pmatrix} }[/math]

according to Cartesian coordinates. The position vector is expressed as:

[math]\displaystyle{ \mathbf{r} = \left( r + R_0 \cos\theta \right) \mathbf{e}_r - R_0 \sin\theta \mathbf{e}_\theta }[/math]

The velocity vector is then given by:

[math]\displaystyle{ \mathbf{\dot{r}} = \dot{r} \mathbf{e}_r + r\dot{\theta} \mathbf{e}_\theta + \dot{\zeta} \left( R_0 + r \cos\theta \right) \mathbf{e}_\zeta }[/math]

and the acceleration vector is:

[math]\displaystyle{ \begin{align} \mathbf{\ddot{r}} = {} & \left( \ddot{r} - r \dot{\theta}^2 - r \dot{\zeta}^2 \cos^2\theta - R_0 \dot{\zeta}^2 \cos\theta \right) \mathbf{e}_r \\[5pt] & {} + \left( 2\dot{r}\dot{\theta} + r\ddot{\theta} + r\dot{\zeta}^2\cos\theta\sin\theta + R_0 \dot{\zeta}^2 \sin\theta \right) \mathbf{e}_\theta \\[5pt] & {} + \left( 2 \dot{r}\dot{\zeta}\cos\theta - 2 r \dot{\theta}\dot{\zeta} \sin\theta + \ddot{\zeta} \left( R_0 + r\cos\theta \right) \right) \mathbf{e}_\zeta \end{align} }[/math]

See also

References