Bipolar coordinates

Bipolar coordinate system

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).

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 F1 and F2. Referring to the figure at right, the σ-coordinate of a point P equals the angle F1 P F2, and the τ-coordinate equals the natural logarithm of the ratio of the distances d1 and d2:

${\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 F1) to ${\displaystyle \infty }$ (for points close to F2). The coordinate σ is only defined modulo , and is best taken to range from to π, by taking it as the negative of the acute angle F1 P F2 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 can be proven by first differentiating x and y with respect to sigma and tau and then reversing the logic of the section below for finding the scale factors.). 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) (the Cauchy-Riemann equations) that these particulse particular curves of σ and τ intersect at right angles, i.e., that the coordinate system is orthogonal. This can be proven by first differentiating x and y with respect to sigma and tau and thenreversing the logic of the section below for finding the scale factors.

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 |σ|- π/2 decreases, the radius of the circles decreases and the center approaches the origin (0, 0), which is reached when |σ| = π/2. (From elementary geometry, all triangles on a circle with 2 vertices on opposite ends of a diameter are right triangles.)

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}}}.}$

which is the limit one would get a x = 0 from the definition in the section above. And all the limits look quite ordinary at x =0.

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 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 copy". Archived from the original on December 12, 2007. Retrieved December 9, 2006.CS1 maint: archived copy as title (link)
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.