Suvarna Garge (Editor)

Oblate spheroidal wave function

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit

In applied mathematics, oblate spheroidal wave functions are involved in the solution of the Helmholtz equation in oblate spheroidal coordinates. When solving this equation, Δ Φ + k 2 Φ = 0 , by the method of separation of variables, ( ξ , η , ϕ ) , with:

  x = ( d / 2 ) ξ η ,   y = ( d / 2 ) ( ξ 2 + 1 ) ( 1 η 2 ) cos ϕ ,   z = ( d / 2 ) ( ξ 2 + 1 ) ( 1 η 2 ) sin ϕ ,   ξ 0  and  | η | 1.

the solution Φ ( ξ , η , ϕ ) can be written as the product of a radial spheroidal wave function R m n ( i c , i ξ ) and an angular spheroidal wave function S m n ( i c , η ) by e i m ϕ . Here c = k d / 2 , with d being the interfocal length of the elliptical cross section of the oblate spheroid.

The radial wave function R m n ( i c , i ξ ) satisfies the linear ordinary differential equation:

  ( ξ 2 + 1 ) d 2 R m n ( i c , i ξ ) d ξ 2 + 2 ξ d R m n ( i c , i ξ ) d ξ ( λ m n ( c ) c 2 ξ 2 m 2 ξ 2 + 1 ) R m n ( i c , i ξ ) = 0 .

The angular wave function satisfies the differential equation:

  ( 1 η 2 ) d 2 S m n ( i c , η ) d η 2 2 η d S m n ( i c , η ) d η + ( λ m n ( c ) + c 2 η 2 m 2 1 η 2 ) S m n ( i c , η ) = 0 .

It is the same differential equation as in the case of the radial wave function. However, the range of the radial coordinate ξ is different from that of the angular coordinate η .

The eigenvalue λ m n ( i c ) of this Sturm-Liouville differential equation is fixed by the requirement that S m n ( i c , η ) must be finite for | η | = 1 .

For c = 0 these two differential equations reduce to the equations satisfied by the associated Legendre polynomials. For c 0 , the angular spheroidal wave functions can be expanded as a series of Legendre functions.

The differential equations given above for the oblate radial and angular wave functions can be obtained from the corresponding equations for the prolate spheroidal wave functions by the substitution of i c for c and i ξ for ξ . The notation for the oblate spheroidal functions reflects this relationship.

There are different normalization schemes for spheroidal functions. A table of the different schemes can be found in Abramowitz and Stegun. Abramowitz and Stegun (and the present article) follow the notation of Flammer.

Originally, the spheroidal wave functions were introduced by C. Niven, which lead to a Helmholtz equation in spheroidal coordinates. Monographs tying together many aspects of the theory of spheroidal wave functions were written by Strutt, Stratton et al., Meixner and Schafke, and Flammer.

Flammer provided a thorough discussion of the calculation of the eigenvalues, angular wavefunctions, and radial wavefunctions for both the oblate and the prolate case. Computer programs for this purpose have been developed by many, including Van Buren et al., King and Van Buren, Baier et al., Zhang and Jin, and Thompson. Van Buren has recently developed new methods for calculating oblate spheroidal wave functions that extend the ability to obtain numerical values to extremely wide parameter ranges. These results are based on earlier work on prolate spheroidal wave functions. Fortran source code that combines the new results with traditional methods is available at http://www.mathieuandspheroidalwavefunctions.com.

Tables of numerical values of oblate spheroidal wave functions are given in Flammer, Hanish et al., and Van Buren et al.

Asymptotic expansions of angular oblate spheroidal wave functions for large values of c have been derived by Müller.

The Digital Library of Mathematical Functions http://dlmf.nist.gov provided by NIST is an excellent resource for spheroidal wave functions.

References

Oblate spheroidal wave function Wikipedia