Elliptic coordinates

☆ Save On Wikipedia ↗
Elliptic coordinate system

In geometry, the elliptic coordinate system is a two-dimensional orthogonal coordinate system in which the coordinate lines are confocal ellipses and hyperbolae. The two foci F 1 {\displaystyle F_{1}} {\displaystyle F_{1}} and F 2 {\displaystyle F_{2}} {\displaystyle F_{2}} are generally taken to be fixed at − a {\displaystyle -a} {\displaystyle -a} and + a {\displaystyle +a} {\displaystyle +a}, respectively, on the x {\displaystyle x} {\displaystyle x}-axis of the Cartesian coordinate system.

Basic definition

The most common definition of elliptic coordinates ( μ , ν ) {\displaystyle (\mu ,\nu )} {\displaystyle (\mu ,\nu )} is

x = a   cosh ⁡ μ   cos ⁡ ν y = a   sinh ⁡ μ   sin ⁡ ν {\displaystyle {\begin{aligned}x&=a\ \cosh \mu \ \cos \nu \\y&=a\ \sinh \mu \ \sin \nu \end{aligned}}} {\displaystyle {\begin{aligned}x&=a\ \cosh \mu \ \cos \nu \\y&=a\ \sinh \mu \ \sin \nu \end{aligned}}}

where μ {\displaystyle \mu } {\displaystyle \mu } is a nonnegative real number and ν ∈ [ 0 , 2 π ] . {\displaystyle \nu \in [0,2\pi ].} {\displaystyle \nu \in [0,2\pi ].}

On the complex plane, an equivalent relationship is

x + i y = a   cosh ⁡ ( μ + i ν ) {\displaystyle x+iy=a\ \cosh(\mu +i\nu )} {\displaystyle x+iy=a\ \cosh(\mu +i\nu )}

These definitions correspond to ellipses and hyperbolae. The trigonometric identity

x 2 a 2 cosh 2 ⁡ μ + y 2 a 2 sinh 2 ⁡ μ = cos 2 ⁡ ν + sin 2 ⁡ ν = 1 {\displaystyle {\frac {x^{2}}{a^{2}\cosh ^{2}\mu }}+{\frac {y^{2}}{a^{2}\sinh ^{2}\mu }}=\cos ^{2}\nu +\sin ^{2}\nu =1} {\displaystyle {\frac {x^{2}}{a^{2}\cosh ^{2}\mu }}+{\frac {y^{2}}{a^{2}\sinh ^{2}\mu }}=\cos ^{2}\nu +\sin ^{2}\nu =1}

shows that curves of constant μ {\displaystyle \mu } {\displaystyle \mu } form ellipses, whereas the hyperbolic trigonometric identity

x 2 a 2 cos 2 ⁡ ν − y 2 a 2 sin 2 ⁡ ν = cosh 2 ⁡ μ − sinh 2 ⁡ μ = 1 {\displaystyle {\frac {x^{2}}{a^{2}\cos ^{2}\nu }}-{\frac {y^{2}}{a^{2}\sin ^{2}\nu }}=\cosh ^{2}\mu -\sinh ^{2}\mu =1} {\displaystyle {\frac {x^{2}}{a^{2}\cos ^{2}\nu }}-{\frac {y^{2}}{a^{2}\sin ^{2}\nu }}=\cosh ^{2}\mu -\sinh ^{2}\mu =1}

shows that curves of constant ν {\displaystyle \nu } {\displaystyle \nu } form hyperbolae.

Scale factors

In an orthogonal coordinate system the lengths of the basis vectors are known as scale factors. The scale factors for the elliptic coordinates ( μ , ν ) {\displaystyle (\mu ,\nu )} {\displaystyle (\mu ,\nu )} are equal to

h μ = h ν = a sinh 2 ⁡ μ + sin 2 ⁡ ν = a cosh 2 ⁡ μ − cos 2 ⁡ ν . {\displaystyle h_{\mu }=h_{\nu }=a{\sqrt {\sinh ^{2}\mu +\sin ^{2}\nu }}=a{\sqrt {\cosh ^{2}\mu -\cos ^{2}\nu }}.} {\displaystyle h_{\mu }=h_{\nu }=a{\sqrt {\sinh ^{2}\mu +\sin ^{2}\nu }}=a{\sqrt {\cosh ^{2}\mu -\cos ^{2}\nu }}.}

Using the double argument identities for hyperbolic functions and trigonometric functions, the scale factors can be equivalently expressed as

h μ = h ν = a 1 2 ( cosh ⁡ 2 μ − cos ⁡ 2 ν ) . {\displaystyle h_{\mu }=h_{\nu }=a{\sqrt {{\frac {1}{2}}(\cosh 2\mu -\cos 2\nu )}}.} {\displaystyle h_{\mu }=h_{\nu }=a{\sqrt {{\frac {1}{2}}(\cosh 2\mu -\cos 2\nu )}}.}

Consequently, an infinitesimal element of area equals

d A = h μ h ν d μ d ν = a 2 ( sinh 2 ⁡ μ + sin 2 ⁡ ν ) d μ d ν = a 2 ( cosh 2 ⁡ μ − cos 2 ⁡ ν ) d μ d ν = a 2 2 ( cosh ⁡ 2 μ − cos ⁡ 2 ν ) d μ d ν {\displaystyle {\begin{aligned}dA&=h_{\mu }h_{\nu }d\mu d\nu \\&=a^{2}\left(\sinh ^{2}\mu +\sin ^{2}\nu \right)d\mu d\nu \\&=a^{2}\left(\cosh ^{2}\mu -\cos ^{2}\nu \right)d\mu d\nu \\&={\frac {a^{2}}{2}}\left(\cosh 2\mu -\cos 2\nu \right)d\mu d\nu \end{aligned}}} {\displaystyle {\begin{aligned}dA&=h_{\mu }h_{\nu }d\mu d\nu \\&=a^{2}\left(\sinh ^{2}\mu +\sin ^{2}\nu \right)d\mu d\nu \\&=a^{2}\left(\cosh ^{2}\mu -\cos ^{2}\nu \right)d\mu d\nu \\&={\frac {a^{2}}{2}}\left(\cosh 2\mu -\cos 2\nu \right)d\mu d\nu \end{aligned}}}

and the Laplacian reads

∇ 2 Φ = 1 a 2 ( sinh 2 ⁡ μ + sin 2 ⁡ ν ) ( ∂ 2 Φ ∂ μ 2 + ∂ 2 Φ ∂ ν 2 ) = 1 a 2 ( cosh 2 ⁡ μ − cos 2 ⁡ ν ) ( ∂ 2 Φ ∂ μ 2 + ∂ 2 Φ ∂ ν 2 ) = 2 a 2 ( cosh ⁡ 2 μ − cos ⁡ 2 ν ) ( ∂ 2 Φ ∂ μ 2 + ∂ 2 Φ ∂ ν 2 ) {\displaystyle {\begin{aligned}\nabla ^{2}\Phi &={\frac {1}{a^{2}\left(\sinh ^{2}\mu +\sin ^{2}\nu \right)}}\left({\frac {\partial ^{2}\Phi }{\partial \mu ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \nu ^{2}}}\right)\\&={\frac {1}{a^{2}\left(\cosh ^{2}\mu -\cos ^{2}\nu \right)}}\left({\frac {\partial ^{2}\Phi }{\partial \mu ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \nu ^{2}}}\right)\\&={\frac {2}{a^{2}\left(\cosh 2\mu -\cos 2\nu \right)}}\left({\frac {\partial ^{2}\Phi }{\partial \mu ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \nu ^{2}}}\right)\end{aligned}}} {\displaystyle {\begin{aligned}\nabla ^{2}\Phi &={\frac {1}{a^{2}\left(\sinh ^{2}\mu +\sin ^{2}\nu \right)}}\left({\frac {\partial ^{2}\Phi }{\partial \mu ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \nu ^{2}}}\right)\\&={\frac {1}{a^{2}\left(\cosh ^{2}\mu -\cos ^{2}\nu \right)}}\left({\frac {\partial ^{2}\Phi }{\partial \mu ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \nu ^{2}}}\right)\\&={\frac {2}{a^{2}\left(\cosh 2\mu -\cos 2\nu \right)}}\left({\frac {\partial ^{2}\Phi }{\partial \mu ^{2}}}+{\frac {\partial ^{2}\Phi }{\partial \nu ^{2}}}\right)\end{aligned}}}

Other differential operators such as ∇ ⋅ F {\displaystyle \nabla \cdot \mathbf {F} } {\displaystyle \nabla \cdot \mathbf {F} } and ∇ × F {\displaystyle \nabla \times \mathbf {F} } {\displaystyle \nabla \times \mathbf {F} } can be expressed in the coordinates ( μ , ν ) {\displaystyle (\mu ,\nu )} {\displaystyle (\mu ,\nu )} by substituting the scale factors into the general formulae found in orthogonal coordinates.

Alternative definition

An alternative and geometrically intuitive set of elliptic coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} {\displaystyle (\sigma ,\tau )} are sometimes used, where σ = cosh ⁡ μ {\displaystyle \sigma =\cosh \mu } {\displaystyle \sigma =\cosh \mu } and τ = cos ⁡ ν {\displaystyle \tau =\cos \nu } {\displaystyle \tau =\cos \nu }. Hence, the curves of constant σ {\displaystyle \sigma } {\displaystyle \sigma } are ellipses, whereas the curves of constant τ {\displaystyle \tau } {\displaystyle \tau } are hyperbolae. The coordinate τ {\displaystyle \tau } {\displaystyle \tau } must belong to the interval [-1, 1], whereas the σ {\displaystyle \sigma } {\displaystyle \sigma } coordinate must be greater than or equal to one.

The coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} {\displaystyle (\sigma ,\tau )} have a simple relation to the distances to the foci F 1 {\displaystyle F_{1}} {\displaystyle F_{1}} and F 2 {\displaystyle F_{2}} {\displaystyle F_{2}}. For any point in the plane, the sum d 1 + d 2 {\displaystyle d_{1}+d_{2}} {\displaystyle d_{1}+d_{2}} of its distances to the foci equals 2 a σ {\displaystyle 2a\sigma } {\displaystyle 2a\sigma }, whereas their difference d 1 − d 2 {\displaystyle d_{1}-d_{2}} {\displaystyle d_{1}-d_{2}} equals 2 a τ {\displaystyle 2a\tau } {\displaystyle 2a\tau }. Thus, the distance to F 1 {\displaystyle F_{1}} {\displaystyle F_{1}} is a ( σ + τ ) {\displaystyle a(\sigma +\tau )} {\displaystyle a(\sigma +\tau )}, whereas the distance to F 2 {\displaystyle F_{2}} {\displaystyle F_{2}} is a ( σ − τ ) {\displaystyle a(\sigma -\tau )} {\displaystyle a(\sigma -\tau )}. (Recall that F 1 {\displaystyle F_{1}} {\displaystyle F_{1}} and F 2 {\displaystyle F_{2}} {\displaystyle F_{2}} are located at x = − a {\displaystyle x=-a} {\displaystyle x=-a} and x = + a {\displaystyle x=+a} {\displaystyle x=+a}, respectively.)

A drawback of these coordinates is that the points with Cartesian coordinates (x,y) and (x,-y) have the same coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} {\displaystyle (\sigma ,\tau )}, so the conversion to Cartesian coordinates is not a function, but a multifunction.

x = a σ τ {\displaystyle x=a\left.\sigma \right.\tau } {\displaystyle x=a\left.\sigma \right.\tau }
y 2 = a 2 ( σ 2 − 1 ) ( 1 − τ 2 ) . {\displaystyle y^{2}=a^{2}\left(\sigma ^{2}-1\right)\left(1-\tau ^{2}\right).} {\displaystyle y^{2}=a^{2}\left(\sigma ^{2}-1\right)\left(1-\tau ^{2}\right).}

Alternative scale factors

The scale factors for the alternative elliptic coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} {\displaystyle (\sigma ,\tau )} are

h σ = a σ 2 − τ 2 σ 2 − 1 {\displaystyle h_{\sigma }=a{\sqrt {\frac {\sigma ^{2}-\tau ^{2}}{\sigma ^{2}-1}}}} {\displaystyle h_{\sigma }=a{\sqrt {\frac {\sigma ^{2}-\tau ^{2}}{\sigma ^{2}-1}}}}
h τ = a σ 2 − τ 2 1 − τ 2 . {\displaystyle h_{\tau }=a{\sqrt {\frac {\sigma ^{2}-\tau ^{2}}{1-\tau ^{2}}}}.} {\displaystyle h_{\tau }=a{\sqrt {\frac {\sigma ^{2}-\tau ^{2}}{1-\tau ^{2}}}}.}

Hence, the infinitesimal area element becomes

d A = a 2 σ 2 − τ 2 ( σ 2 − 1 ) ( 1 − τ 2 ) d σ d τ {\displaystyle dA=a^{2}{\frac {\sigma ^{2}-\tau ^{2}}{\sqrt {\left(\sigma ^{2}-1\right)\left(1-\tau ^{2}\right)}}}d\sigma d\tau } {\displaystyle dA=a^{2}{\frac {\sigma ^{2}-\tau ^{2}}{\sqrt {\left(\sigma ^{2}-1\right)\left(1-\tau ^{2}\right)}}}d\sigma d\tau }

and the Laplacian equals

∇ 2 Φ = 1 a 2 ( σ 2 − τ 2 ) [ σ 2 − 1 ∂ ∂ σ ( σ 2 − 1 ∂ Φ ∂ σ ) + 1 − τ 2 ∂ ∂ τ ( 1 − τ 2 ∂ Φ ∂ τ ) ] . {\displaystyle \nabla ^{2}\Phi ={\frac {1}{a^{2}\left(\sigma ^{2}-\tau ^{2}\right)}}\left[{\sqrt {\sigma ^{2}-1}}{\frac {\partial }{\partial \sigma }}\left({\sqrt {\sigma ^{2}-1}}{\frac {\partial \Phi }{\partial \sigma }}\right)+{\sqrt {1-\tau ^{2}}}{\frac {\partial }{\partial \tau }}\left({\sqrt {1-\tau ^{2}}}{\frac {\partial \Phi }{\partial \tau }}\right)\right].} {\displaystyle \nabla ^{2}\Phi ={\frac {1}{a^{2}\left(\sigma ^{2}-\tau ^{2}\right)}}\left[{\sqrt {\sigma ^{2}-1}}{\frac {\partial }{\partial \sigma }}\left({\sqrt {\sigma ^{2}-1}}{\frac {\partial \Phi }{\partial \sigma }}\right)+{\sqrt {1-\tau ^{2}}}{\frac {\partial }{\partial \tau }}\left({\sqrt {1-\tau ^{2}}}{\frac {\partial \Phi }{\partial \tau }}\right)\right].}

Other differential operators such as ∇ ⋅ F {\displaystyle \nabla \cdot \mathbf {F} } {\displaystyle \nabla \cdot \mathbf {F} } and ∇ × F {\displaystyle \nabla \times \mathbf {F} } {\displaystyle \nabla \times \mathbf {F} } can be expressed in the coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} {\displaystyle (\sigma ,\tau )} by substituting the scale factors into the general formulae found in orthogonal coordinates.

Extrapolation to higher dimensions

Elliptic coordinates form the basis for several sets of three-dimensional orthogonal coordinates:

  1. The elliptic cylindrical coordinates are produced by projecting in the z {\displaystyle z} {\displaystyle z}-direction.
  2. The prolate spheroidal coordinates are produced by rotating the elliptic coordinates about the x {\displaystyle x} {\displaystyle x}-axis, i.e., the axis connecting the foci, whereas the oblate spheroidal coordinates are produced by rotating the elliptic coordinates about the y {\displaystyle y} {\displaystyle y}-axis, i.e., the axis separating the foci.
  3. Ellipsoidal coordinates are a formal extension of elliptic coordinates into 3-dimensions, which is based on confocal ellipsoids, hyperboloids of one and two sheets.

Note that (ellipsoidal) Geographic coordinate system is a different concept from above.

Applications

The classic applications of elliptic coordinates are in solving partial differential equations, e.g., Laplace's equation or the Helmholtz equation, for which elliptic coordinates are a natural description of a system thus allowing a separation of variables in the partial differential equations. Some traditional examples are solving systems such as electrons orbiting a molecule or planetary orbits that have an elliptical shape.

The geometric properties of elliptic coordinates can also be useful. A typical example might involve an integration over all pairs of vectors p {\displaystyle \mathbf {p} } {\displaystyle \mathbf {p} } and q {\displaystyle \mathbf {q} } {\displaystyle \mathbf {q} } that sum to a fixed vector r = p + q {\displaystyle \mathbf {r} =\mathbf {p} +\mathbf {q} } {\displaystyle \mathbf {r} =\mathbf {p} +\mathbf {q} }, where the integrand was a function of the vector lengths | p | {\displaystyle \left|\mathbf {p} \right|} {\displaystyle \left|\mathbf {p} \right|} and | q | {\displaystyle \left|\mathbf {q} \right|} {\displaystyle \left|\mathbf {q} \right|}. (In such a case, one would position r {\displaystyle \mathbf {r} } {\displaystyle \mathbf {r} } between the two foci and aligned with the x {\displaystyle x} {\displaystyle x}-axis, i.e., r = 2 a x ^ {\displaystyle \mathbf {r} =2a\mathbf {\hat {x}} } {\displaystyle \mathbf {r} =2a\mathbf {\hat {x}} }.) For concreteness, r {\displaystyle \mathbf {r} } {\displaystyle \mathbf {r} }, p {\displaystyle \mathbf {p} } {\displaystyle \mathbf {p} } and q {\displaystyle \mathbf {q} } {\displaystyle \mathbf {q} } could represent the momenta of a particle and its decomposition products, respectively, and the integrand might involve the kinetic energies of the products (which are proportional to the squared lengths of the momenta).

See also

References