Bipolar coordinates

Bipolar coordinates are a two-dimensional orthogonal coordinate system based on the Apollonian circles.[1]. Confusingly, the same term is also sometimes used for two-center bipolar coordinates. There is also a third system, based on two poles (biangular coordinates).

See also: two-center bipolar coordinates

Bipolar coordinate system

The term "bipolar" is further used on occasion to describe other curves having two singular points (foci), such as ellipses, hyperbolas, and Cassini ovals. However, the term bipolar coordinates is reserved for the coordinates described here, and never used for systems associated with those other curves, such as elliptic coordinates.

Geometric interpretation of the bipolar coordinates. The angle σ is formed by the two foci and the point P, whereas τ is the logarithm of the ratio of distances to the foci. The corresponding circles of constant σ and τ are shown in red and blue, respectively, and meet at right angles (magenta box); they are orthogonal.

Definition

The system is based on two foci F₁ and F₂. Referring to the figure at right, the σ-coordinate of a point P equals the angle F₁ P F₂, and the τ-coordinate equals the natural logarithm of the ratio of the distances d₁ and d₂:

${\displaystyle \tau =\ln {\frac {d_{1}}{d_{2}}}.}$

If, in the Cartesian system, the foci are taken to lie at (−a, 0) and (a, 0), the coordinates of the point P are

${\displaystyle x=a\ {\frac {\sinh \tau }{\cosh \tau -\cos \sigma }},\qquad y=a\ {\frac {\sin \sigma }{\cosh \tau -\cos \sigma }}.}$

The coordinate τ ranges from ${\displaystyle -\infty }$ (for points close to F₁) to ${\displaystyle \infty }$ (for points close to F₂). The coordinate σ is only defined modulo 2π, and is best taken to range from -π to π, by taking it as the negative of the acute angle F₁ P F₂ if P is in the lower half plane.

Proof that coordinate system is orthogonal

The equations for x and y can be combined to give

${\displaystyle x+iy=ai\cot \left({\frac {\sigma +i\tau }{2}}\right).}$[2][3]

This equation shows that σ and τ are the real and imaginary parts of an analytic function of x+iy (with logarithmic branch points at the foci), which in turn proves (by appeal to the general theory of conformal mapping) that the curves of σ and τ intersect at right angles, i.e., that the coordinate system is orthogonal.

Curves of constant σ and τ

The curves of constant σ correspond to non-concentric circles

${\displaystyle x^{2}+\left(y-a\cot \sigma \right)^{2}={\frac {a^{2}}{\sin ^{2}\sigma }}}$

that intersect at the two foci. The centers of the constant-σ circles lie on the y-axis. Circles of positive σ are centered above the x-axis, whereas those of negative σ lie below the axis. As the magnitude |σ| increases, the radius of the circles decreases and the center approaches the origin (0, 0), which is reached when |σ| = π.

The curves of constant ${\displaystyle \tau }$ are non-intersecting circles of different radii

${\displaystyle y^{2}+\left(x-a\coth \tau \right)^{2}={\frac {a^{2}}{\sinh ^{2}\tau }}}$

that surround the foci but again are not concentric. The centers of the constant-τ circles lie on the x-axis. The circles of positive τ lie in the right-hand side of the plane (x > 0), whereas the circles of negative τ lie in the left-hand side of the plane (x < 0). The τ = 0 curve corresponds to the y-axis (x = 0). As the magnitude of τ increases, the radius of the circles decreases and their centers approach the foci.

Reciprocal relations

The passage from the Cartesian coordinates towards the bipolar coordinates can be done via the following formulas:

${\displaystyle \tau ={\frac {1}{2}}\ln {\frac {(x+a)^{2}+y^{2}}{(x-a)^{2}+y^{2}}}}$

and

${\displaystyle \pi -\sigma =2\arctan {\frac {2ay}{a^{2}-x^{2}-y^{2}+{\sqrt {(a^{2}-x^{2}-y^{2})^{2}+4a^{2}y^{2}}}}}.}$

The coordinates also have the identities:

${\displaystyle \tanh \tau ={\frac {2ax}{x^{2}+y^{2}+a^{2}}}}$

and

${\displaystyle \tan \sigma ={\frac {2ay}{x^{2}+y^{2}-a^{2}}}.}$

Scale factors

To obtain the scale factors for bipolar coordinates, we take the differential of the equation for ${\displaystyle x+iy}$, which gives

${\displaystyle dx+i\,dy={\frac {-ia}{\sin ^{2}{\bigl (}{\tfrac {1}{2}}(\sigma +i\tau ){\bigr )}}}(d\sigma +i\,d\tau ).}$

Multiplying this equation with its complex conjugate yields

${\displaystyle (dx)^{2}+(dy)^{2}={\frac {a^{2}}{{\bigl [}2\sin {\tfrac {1}{2}}{\bigl (}\sigma +i\tau {\bigr )}\sin {\tfrac {1}{2}}{\bigl (}\sigma -i\tau {\bigr )}{\bigr ]}^{2}}}{\bigl (}(d\sigma )^{2}+(d\tau )^{2}{\bigr )}.}$

Employing the trigonometric identities for products of sines and cosines, we obtain

${\displaystyle 2\sin {\tfrac {1}{2}}{\bigl (}\sigma +i\tau {\bigr )}\sin {\tfrac {1}{2}}{\bigl (}\sigma -i\tau {\bigr )}=\cos \sigma -\cosh \tau ,}$

from which it follows that

${\displaystyle (dx)^{2}+(dy)^{2}={\frac {a^{2}}{(\cosh \tau -\cos \sigma )^{2}}}{\bigl (}(d\sigma )^{2}+(d\tau )^{2}{\bigr )}.}$

Hence the scale factors for σ and τ are equal, and given by

${\displaystyle h_{\sigma }=h_{\tau }={\frac {a}{\cosh \tau -\cos \sigma }}.}$

Many results now follow in quick succession from the general formulae for orthogonal coordinates. Thus, the infinitesimal area element equals

${\displaystyle dA={\frac {a^{2}}{\left(\cosh \tau -\cos \sigma \right)^{2}}}\,d\sigma \,d\tau ,}$

and the Laplacian is given by

${\displaystyle \nabla ^{2}\Phi ={\frac {1}{a^{2}}}\left(\cosh \tau -\cos \sigma \right)^{2}\left({\frac {\partial ^{2}\Phi }{\partial \sigma ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \tau ^{2}}}\right).}$

Expressions for ${\displaystyle \nabla f}$, ${\displaystyle \nabla \cdot \mathbf {F} }$, and ${\displaystyle \nabla \times \mathbf {F} }$ can be expressed obtained by substituting the scale factors into the general formulae found in orthogonal coordinates.

Applications

The classic applications of bipolar coordinates are in solving partial differential equations, e.g., Laplace's equation or the Helmholtz equation, for which bipolar coordinates allow a separation of variables. An example is the electric field surrounding two parallel cylindrical conductors with unequal diameters.

Polar plotters use bipolar coordinates to describe the drawing paths required to draw a target image.

Extension to 3-dimensions

Bipolar coordinates form the basis for several sets of three-dimensional orthogonal coordinates.

• The bipolar cylindrical coordinates are produced by translating the bipolar coordinates along the z-axis, i.e., the out-of-plane axis.

• The bispherical coordinates are produced by rotating the bipolar coordinates about the x-axis, i.e., the axis connecting the foci.

• The toroidal coordinates are produced by rotating the bipolar coordinates about the y-axis, i.e., the axis separating the foci.

References

1. Eric W. Weisstein, Concise Encyclopedia of Mathematics CD-ROM, Bipolar Coordinates, CD-ROM edition 1.0, May 20, 1999 Archived December 12, 2007, at the Wayback Machine
2. Polyanin, Andrei Dmitrievich (2002). Handbook of linear partial differential equations for engineers and scientists. CRC Press. p. 476. ISBN 1-58488-299-9.
3. Happel, John; Brenner, Howard (1983). Low Reynolds number hydrodynamics: with special applications to particulate media. Mechanics of fluids and transport processes. 1. Springer. p. 497. ISBN 978-90-247-2877-0.

• "Bipolar coordinates", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
• Korn GA and Korn TM. (1961) Mathematical Handbook for Scientists and Engineers, McGraw-Hill.