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}} and F 2 {\displaystyle F_{2}} are generally taken to be fixed at − a {\displaystyle -a} and + a {\displaystyle +a} , respectively, on the x {\displaystyle x} -axis of the Cartesian coordinate system.
Basic definition
The most common definition of elliptic coordinates ( μ , ν ) {\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}}}where μ {\displaystyle \mu } is a nonnegative real number and ν ∈ [ 0 , 2 π ] . {\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 )}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}shows that curves of constant μ {\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}shows that curves of constant ν {\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 )} 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 }}.}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 )}}.}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}}}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}}}Other differential operators such as ∇ ⋅ F {\displaystyle \nabla \cdot \mathbf {F} } and ∇ × F {\displaystyle \nabla \times \mathbf {F} } can be expressed in the coordinates ( μ , ν ) {\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 )} are sometimes used, where σ = cosh μ {\displaystyle \sigma =\cosh \mu } and τ = cos ν {\displaystyle \tau =\cos \nu } . Hence, the curves of constant σ {\displaystyle \sigma } are ellipses, whereas the curves of constant τ {\displaystyle \tau } are hyperbolae. The coordinate τ {\displaystyle \tau } must belong to the interval [-1, 1], whereas the σ {\displaystyle \sigma } coordinate must be greater than or equal to one.
The coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} have a simple relation to the distances to the foci F 1 {\displaystyle F_{1}} and F 2 {\displaystyle F_{2}} . For any point in the plane, the sum d 1 + d 2 {\displaystyle d_{1}+d_{2}} of its distances to the foci equals 2 a σ {\displaystyle 2a\sigma } , whereas their difference d 1 − d 2 {\displaystyle d_{1}-d_{2}} equals 2 a τ {\displaystyle 2a\tau } . Thus, the distance to F 1 {\displaystyle F_{1}} is a ( σ + τ ) {\displaystyle a(\sigma +\tau )} , whereas the distance to F 2 {\displaystyle F_{2}} is a ( σ − τ ) {\displaystyle a(\sigma -\tau )} . (Recall that F 1 {\displaystyle F_{1}} and F 2 {\displaystyle F_{2}} are located at x = − a {\displaystyle x=-a} and 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 )} , so the conversion to Cartesian coordinates is not a function, but a multifunction.
x = a σ τ {\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).}Alternative scale factors
The scale factors for the alternative elliptic coordinates ( σ , τ ) {\displaystyle (\sigma ,\tau )} are
h σ = a σ 2 − τ 2 σ 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}}}}.}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 }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].}Other differential operators such as ∇ ⋅ F {\displaystyle \nabla \cdot \mathbf {F} } and ∇ × F {\displaystyle \nabla \times \mathbf {F} } can be expressed in the coordinates ( σ , τ ) {\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:
- The elliptic cylindrical coordinates are produced by projecting in the z {\displaystyle z} -direction.
- The prolate spheroidal coordinates are produced by rotating the elliptic coordinates about the 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} -axis, i.e., the axis separating the foci.
- 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} } and q {\displaystyle \mathbf {q} } that sum to a fixed vector r = p + q {\displaystyle \mathbf {r} =\mathbf {p} +\mathbf {q} } , where the integrand was a function of the vector lengths | p | {\displaystyle \left|\mathbf {p} \right|} and | q | {\displaystyle \left|\mathbf {q} \right|} . (In such a case, one would position r {\displaystyle \mathbf {r} } between the two foci and aligned with the x {\displaystyle x} -axis, i.e., r = 2 a x ^ {\displaystyle \mathbf {r} =2a\mathbf {\hat {x}} } .) For concreteness, r {\displaystyle \mathbf {r} } , p {\displaystyle \mathbf {p} } and 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
- "Elliptic coordinates", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
- Korn GA and Korn TM. (1961) Mathematical Handbook for Scientists and Engineers, McGraw-Hill.
- Weisstein, Eric W. "Elliptic Cylindrical Coordinates." From MathWorld — A Wolfram Web Resource. http://mathworld.wolfram.com/EllipticCylindricalCoordinates.html