Elliptic Cylinder Coordinates

Background

Building on our general introduction to Direction Cosines in the context of orthogonal curvilinear coordinate systems, here we detail the properties of Elliptic Cylinder Coordinates. First, we will present this coordinate system in the manner described by [MF53]; second, we will provide an alternate presentation, obtained from Wikipedia; then, third, we will investigate whether or not a related coordinate system based on concentric (rather than confocal) elliptic surfaces can be satisfactorily described.

It is useful to keep in mind various properties of a set of confocal ellipses in which the location of the pair of foci is fixed at, <math>~(x, y) = (\pm~ c, 0)</math>, and the semi-major axis, <math>~a</math>, is the parameter. The relevant prescriptive relation is,

The semi-minor axis length, <math>~b</math>, and the eccentricity, <math>~e</math>, of the ellipse are, respectively,

and,

<math>~e\equiv \biggl[1 - \frac{b^2}{a^2} \biggr]^{1 / 2}</math>

The length, <math>~\ell_1</math>, of the chord that connects one focus to a point, <math>~P(x,y)</math>, on the ellipse is,

<math>~a + \biggl(\frac{c}{a}\biggr)x \, ;</math>

and the length, <math>~\ell_2</math>, of the chord that connects the second focus to that same point on the ellipse is,

<math>~a - \biggl(\frac{c}{a}\biggr)x \, .</math>

It is easy to see that, for any point on the ellipse, the sum of these two lengths is, <math>~2a</math>. It is worth noting as well that the associated <math>~y</math> coordinate of the relevant point can be obtained from the relation,

<math>~ \ell_1^2</math>	<math>~=</math>	<math>~y^2 + (c+x)^2</math>
<math>~\Rightarrow~~~ (ay)^2</math>	<math>~=</math>	<math>~(a \ell_1)^2 - (ac+ ax)^2</math>
	<math>~=</math>	<math>~(a^2 + cx )^2 - (ac+ ax)^2</math>
	<math>~=</math>	<math>~(a^4 + 2a^2 cx + c^2x^2) - (a^2c^2 + 2a^2 cx + a^2x^2)</math>
	<math>~=</math>	<math>~(a^4 + c^2x^2) - (a^2c^2 + a^2x^2)</math>
	<math>~=</math>	<math>~(a^2-x^2)(a^2 - c^2) </math>
<math>~\Rightarrow ~~~ y</math>	<math>~=</math>	<math>~\pm~\frac{1}{a}\biggl[ (a^2-x^2)(a^2 - c^2) \biggr]^{1 / 2} \, .</math>

MF53

Definition

From MF53's Table of Separable Coordinates in Three Dimensions (see their Chapter 5, beginning on p. 655), we find the following description of Elliptic Cylinder Coordinates (p. 657).

Elliptic Cylindrical Coordinates
(MF53 Primary Definition)

<math>~x</math>	<math>~=</math>	<math>~\xi_1 \xi_2 </math>
<math>~y</math>	<math>~=</math>	<math>~\biggl[ (\xi_1^2 - d^2)(1 - \xi_2^2) \biggr]^{1 / 2} </math>
<math>~z</math>	<math>~=</math>	<math>~\xi_3 </math>

Alternate Definition

Making the substitutions, <math>~\xi_3 \rightarrow z</math>, <math>~\xi_2 \rightarrow \cos\nu</math>, and <math>~\xi_1 \rightarrow d\cosh\mu</math>, we equally well obtain:

<math>~x</math>	<math>~=</math>	<math>~d\cosh\mu \cdot \cos\nu </math>
<math>~y</math>	<math>~=</math>	<math>~d \sinh\mu \cdot \sin\nu </math>
<math>~z</math>	<math>~=</math>	<math>~z </math>

Notice that,

Hence, as is pointed out in a related Wikipedia discussion, "… this shows that curves of constant <math>~\mu</math> form ellipses." For a given choice of <math>~\mu</math> — say, <math>~\mu_0</math> — let's see how the shape of the resulting ellipse relates to the background ellipses discussed above.

Scale Factors

Primary

Appreciating that,

<math>~\frac{\partial y}{\partial \xi_1}</math>	<math>~=</math>	<math>~ +\biggl[ (\xi_1^2 - d^2)(1 - \xi_2^2) \biggr]^{- 1 / 2}\xi_1(1-\xi_2^2) \, , </math> and that,
<math>~\frac{\partial y}{\partial \xi_2}</math>	<math>~=</math>	<math>~ - \biggl[ (\xi_1^2 - d^2)(1 - \xi_2^2) \biggr]^{- 1 / 2}\xi_2(\xi_1^2 - d^2) \, , </math>

we find that the respective scale factors are given by the expressions,

<math>~ h_1^2</math>	<math>~=</math>	<math>~\biggl(\frac{\partial x}{\partial\xi_1} \biggr)^2 + \biggl(\frac{\partial y}{\partial\xi_1} \biggr)^2 + \biggl(\frac{\partial z}{\partial\xi_1} \biggr)^2 </math>
	<math>~=</math>	<math>~\xi_2^2 +\biggl[ (\xi_1^2 - d^2)(1 - \xi_2^2) \biggr]^{- 1 }\xi_1^2 (1-\xi_2^2)^2 </math>
	<math>~=</math>	<math>~ (\xi_1^2 - d^2)^{- 1 } [ (\xi_1^2 - d^2)\xi_2^2 +\xi_1^2 (1-\xi_2^2) ]</math>
	<math>~=</math>	<math>~ \biggl[ \frac{ \xi_1^2 - d^2 \xi_2^2 }{\xi_1^2 - d^2} \biggr] \, ;</math>
<math>~ h_2^2</math>	<math>~=</math>	<math>~\biggl(\frac{\partial x}{\partial\xi_2} \biggr)^2 + \biggl(\frac{\partial y}{\partial\xi_2} \biggr)^2 + \biggl(\frac{\partial z}{\partial\xi_2} \biggr)^2 </math>
	<math>~=</math>	<math>~\xi_1^2 + \biggl[ (\xi_1^2 - d^2)(1 - \xi_2^2) \biggr]^{- 1 }\xi_2^2(\xi_1^2 - d^2)^2 </math>
	<math>~=</math>	<math>~(1 - \xi_2^2)^{- 1 } [\xi_1^2(1 - \xi_2^2) + \xi_2^2(\xi_1^2 - d^2) ]</math>
	<math>~=</math>	<math>~\biggl[ \frac{ \xi_1^2 - d^2 \xi_2^2 }{1 - \xi_2^2} \biggr] \, ;</math>
<math>~ h_3^2</math>	<math>~=</math>	<math>~\biggl(\frac{\partial x}{\partial\xi_3} \biggr)^2 + \biggl(\frac{\partial y}{\partial\xi_3} \biggr)^2 + \biggl(\frac{\partial z}{\partial\xi_3} \biggr)^2 </math>
	<math>~=</math>	<math>~1 \, . </math>

These match the scale-factor expressions found in MF53.

Alternatively

Alternatively, the Wikipedia discussion gives,

<math>~h_\mu = h_\nu</math>	<math>~=</math>	<math>~d\sqrt{ \sinh^2\mu + \sin^2\nu}</math>
<math>~\nabla^2\Phi</math>	<math>~=</math>	<math>~ \frac{1}{d^2(\sinh^2\mu + \sin^2\nu)} \biggl[ \frac{\partial^2 \Phi}{\partial \mu^2} + \frac{\partial^2 \Phi}{\partial \nu^2} \biggr] + \frac{\partial^2 \Phi}{\partial z^2} \, . </math>

Inverting Coordinate Mapping

Inverting the original coordinate mappings, we find,

<math>~y^2</math>	<math>~=</math>	<math>~(\xi_1^2 - a^2)\biggl[ 1 - \biggl(\frac{x}{\xi_1}\biggr)^2 \biggr] </math>
<math>~\Rightarrow ~~~0</math>	<math>~=</math>	<math>~(\xi_1^2 - a^2) ( \xi_1^2 - x^2 ) - \xi_1^2 y^2</math>
	<math>~=</math>	<math>~(\xi_1^2 - a^2) \xi_1^2 - (\xi_1^2 - a^2) x^2 - \xi_1^2 y^2</math>
	<math>~=</math>	<math>~ \xi_1^4 - \xi_1^2 (a^2 + x^2 + y^2) + a^2 x^2 </math>
<math>~\Rightarrow~~~ \xi_1^2</math>	<math>~=</math>	<math>~ \frac{1}{2}\biggl\{ (a^2 + x^2 + y^2) \pm \biggl[ (a^2 + x^2 + y^2)^2 - 4a^2 x^2 \biggr]^{1 / 2} \biggr\} </math>

Only the superior — that is, only the positive — sign will ensure positive values of <math>~\xi_1^2</math>, so in summary we have,

Coordinate Transformation
<math>~\xi_1</math>	<math>~=</math>	<math>~ \frac{1}{\sqrt{2}}\biggl\{ \biggl[ (a^2 + x^2 + y^2)^2 - 4a^2 x^2\biggr]^{1 / 2} + (a^2 + x^2 + y^2) \biggr\}^{1 / 2} \, ; </math>
<math>~\xi_2</math>	<math>~=</math>	<math>~ \frac{x}{\xi_1} \, ; </math>
<math>~\xi_3</math>	<math>~=</math>	<math>~ z \, . </math>

Alternative Wikipedia Definition

This same MF53 coordinate system — with different variable notation — is referred to in a Wikipedia discussion as an "alternative and geometrically intuitive set of elliptic coordinates." The relevant mapping is, <math>~(d\sigma, \tau, z)_\mathrm{Wikipedia} = (\xi_1, \xi_2, \xi_3)_\mathrm{MF53}</math>. The identified mapping to Cartesian coordinates is,

<math>~x</math>	<math>~=</math>	<math>~(d\sigma)\tau </math>	<math>~=</math>	<math>~\xi_1 \xi_2 \, ;</math>
<math>~y</math>	<math>~=</math>	<math>~d \biggl[ (\sigma^2 - 1 )(1 - \tau^2) \biggr]^{1 / 2} </math>	<math>~=</math>	<math>~\biggl[ (\xi_1^2 - d^2)(1 - \xi_2^2) \biggr]^{1 / 2} \, ;</math>
<math>~z</math>	<math>~=</math>	<math>~z</math>	<math>~=</math>	<math>~\xi_3 \, .</math>

According to the Wikipedia discussion, the three scale factors are,

<math>~h_\sigma^2</math>

<math>~ d^2\biggl[\frac{\sigma^2 - \tau^2}{\sigma^2 - 1} \biggr] \, ; </math>

<math>~ d^2\biggl[\frac{\sigma^2 - \tau^2}{1 - \tau^2} \biggr] \, ; </math>

and,

Interestingly, the Wikipedia discussion also includes the following expression for the Laplacian in this elliptic cylindrical coordinate system:

<math>~\nabla^2\Phi</math>

<math>~ \frac{1}{d^2(\sigma^2 - \tau^2)} \biggl[ \sqrt{\sigma^2 - 1} \frac{\partial}{\partial\sigma}\biggl( \sqrt{\sigma^2 - 1} \frac{\partial\Phi}{\partial\sigma} \biggr) + \sqrt{1 - \tau^2 } \frac{\partial}{\partial\tau}\biggl( \sqrt{1 - \tau^2} \frac{\partial\Phi}{\partial\tau} \biggr) \biggr] + \frac{\partial^2\Phi}{\partial z^2} \, . </math>

User:Tohline/Appendix/Ramblings/EllipticCylinderCoordinates

Contents