30.12 Generalized and Coulomb Spheroidal Functions30.14 Wave Equation in Oblate Spheroidal Coordinates

§30.13 Wave Equation in Prolate Spheroidal Coordinates

Contents

§30.13(i) Prolate Spheroidal Coordinates

Prolate spheroidal coordinates \xi,\eta,\phi are related to Cartesian coordinates x,y,z by

30.13.1
x=c\sqrt{(\xi^{2}-1)(1-\eta^{2})}\mathop{\cos\/}\nolimits\phi,
y=c\sqrt{(\xi^{2}-1)(1-\eta^{2})}\mathop{\sin\/}\nolimits\phi,
z=c\xi\eta,

where c is a positive constant. The (x,y,z)-space without the z-axis corresponds to

30.13.2
1<\xi<\infty,
-1<\eta<1,
0\leq\phi<2\pi.

The coordinate surfaces \xi=\mbox{const}. are prolate ellipsoids of revolution with foci at x=y=0, z=\pm c. The coordinate surfaces \eta=\mbox{const}. are sheets of two-sheeted hyperboloids of revolution with the same foci. The focal line is given by \xi=1, -1\leq\eta\leq 1, and the rays \pm z\geq c, x=y=0 are given by \eta=\pm 1, \xi\geq 1.

§30.13(iv) Separation of Variables

The wave equation

30.13.7 \nabla^{2}w+\kappa^{2}w=0,

transformed to prolate spheroidal coordinates (\xi,\eta,\phi), admits solutions

30.13.8 w(\xi,\eta,\phi)=w_{1}(\xi)w_{2}(\eta)w_{3}(\phi),

where w_{1}, w_{2}, w_{3} satisfy the differential equations

30.13.9 \frac{d}{d\xi}\left((1-\xi^{2})\frac{dw_{1}}{d\xi}\right)+\left(\lambda+\gamma^{2}(1-\xi^{2})-\frac{\mu^{2}}{1-\xi^{2}}\right)w_{1}=0,
30.13.10 \frac{d}{d\eta}\left((1-\eta^{2})\frac{dw_{2}}{d\eta}\right)+\left(\lambda+\gamma^{2}(1-\eta^{2})-\frac{\mu^{2}}{1-\eta^{2}}\right)w_{2}=0,
30.13.11 \frac{{d}^{2}w_{3}}{{d\phi}^{2}}+\mu^{2}w_{3}=0,

with \gamma^{2}=\kappa^{2}c^{2}\geq 0 and separation constants \lambda and \mu^{2}. Equations (30.13.9) and (30.13.10) agree with (30.2.1).

In most applications the solution w has to be a single-valued function of (x,y,z), which requires \mu=m (a nonnegative integer) and

30.13.12 w_{3}(\phi)=a_{3}\mathop{\cos\/}\nolimits\!\left(m\phi\right)+b_{3}\mathop{\sin\/}\nolimits\!\left(m\phi\right).

Moreover, w has to be bounded along the z-axis away from the focal line: this requires w_{2}(\eta) to be bounded when -1<\eta<1. Then \lambda=\mathop{\lambda^{{m}}_{{n}}\/}\nolimits\!\left(\gamma^{2}\right) for some n=m,m+1,m+2,\dots, and the general solution of (30.13.10) is

30.13.13 w_{2}(\eta)=a_{2}\mathop{\mathsf{Ps}^{{m}}_{{n}}\/}\nolimits\!\left(\eta,\gamma^{2}\right)+b_{2}\mathop{\mathsf{Qs}^{{m}}_{{n}}\/}\nolimits\!\left(\eta,\gamma^{2}\right).

The solution of (30.13.9) with \mu=m is

30.13.14 w_{1}(\xi)=a_{1}\mathop{S^{{m(1)}}_{{n}}\/}\nolimits\!\left(\xi,\gamma\right)+b_{1}\mathop{S^{{m(2)}}_{{n}}\/}\nolimits\!\left(\xi,\gamma\right).

If b_{1}=b_{2}=0, then the function (30.13.8) is a twice-continuously differentiable solution of (30.13.7) in the entire (x,y,z)-space. If b_{2}=0, then this property holds outside the focal line.

§30.13(v) The Interior Dirichlet Problem for Prolate Ellipsoids

Equation (30.13.7) for \xi\leq\xi _{0}, and subject to the boundary condition w=0 on the ellipsoid given by \xi=\xi _{0}, poses an eigenvalue problem with \kappa^{2} as spectral parameter. The eigenvalues are given by c^{2}\kappa^{2}=\gamma^{2}, where \gamma is determined from the condition

30.13.15 \mathop{S^{{m(1)}}_{{n}}\/}\nolimits\!\left(\xi _{0},\gamma\right)=0.

The corresponding eigenfunctions are given by (30.13.8), (30.13.14), (30.13.13), (30.13.12), with b_{1}=b_{2}=0. For the Dirichlet boundary-value problem of the region \xi _{1}\leq\xi\leq\xi _{2} between two ellipsoids, the eigenvalues are determined from

30.13.16 w_{1}(\xi _{1})=w_{1}(\xi _{2})=0,

with w_{1} as in (30.13.14). The corresponding eigenfunctions are given as before with b_{2}=0.

For further applications see Meixner and Schäfke (1954), Meixner et al. (1980) and the references cited therein; also Ong (1986), Müller et al. (1994), and Xiao et al. (2001).