Astronomy:Double Fourier sphere method

In mathematics, the double Fourier sphere (DFS) method is a simple technique that transforms a function defined on the surface of the sphere to a function defined on a rectangular domain while preserving periodicity in both the longitude and latitude directions.

Introduction

First, a function [math]\displaystyle{ f(x, y, z) }[/math] on the sphere is written as [math]\displaystyle{ f(\lambda,\theta) }[/math] using spherical coordinates, i.e.,

[math]\displaystyle{ f(\lambda,\theta) = f(\cos\lambda\sin\theta,\sin\lambda\sin\theta, \cos\theta), (\lambda,\theta)\in[-\pi,\pi]\times[0,\pi]. }[/math]

The function [math]\displaystyle{ f(\lambda, \theta) }[/math] is [math]\displaystyle{ 2\pi }[/math]-periodic in [math]\displaystyle{ \lambda }[/math], but not periodic in [math]\displaystyle{ \theta }[/math]. The periodicity in the latitude direction has been lost. To recover it, the function is "doubled up” and a related function on [math]\displaystyle{ [-\pi,\pi]\times[-\pi,\pi] }[/math] is defined as

[math]\displaystyle{ \tilde{f}(\lambda,\theta) = \begin{cases} g(\lambda + \pi, \theta), & (\lambda, \theta) \in [-\pi, 0] \times [0, \pi],\\ h(\lambda, \theta), &(\lambda, \theta) \in [0, \pi] \times [0, \pi],\\ g(\lambda, -\theta), &(\lambda, \theta) \in [0, \pi] \times [-\pi, 0],\\ h(\lambda + \pi, -\theta), & (\lambda, \theta) \in [-\pi, 0] \times [-\pi, 0],\\ \end{cases} }[/math]

where [math]\displaystyle{ g(\lambda, \theta) = f(\lambda- \pi, \theta) }[/math] and [math]\displaystyle{ h(\lambda, \theta) = f(\lambda, \theta) }[/math] for [math]\displaystyle{ (\lambda, \theta) \in[0, \pi] \times [0, \pi] }[/math]. The new function [math]\displaystyle{ \tilde{f} }[/math] is [math]\displaystyle{ 2\pi }[/math]-periodic in [math]\displaystyle{ \lambda }[/math] and [math]\displaystyle{ \theta }[/math], and is constant along the lines [math]\displaystyle{ \theta = 0 }[/math] and [math]\displaystyle{ \theta = \pm\pi }[/math], corresponding to the poles.

The function [math]\displaystyle{ \tilde{f} }[/math] can be expanded into a double Fourier series

[math]\displaystyle{ \tilde{f} \approx \sum_{j=-n}^n \sum_{k=-n}^n a_{jk} e^{ij\theta}e^{ik\lambda} }[/math]

History

The DFS method was proposed by Merilees and developed further by Steven Orszag. The DFS method has been the subject of relatively few investigations since (a notable exception is Fornberg's work), perhaps due to the dominance of spherical harmonics expansions. Over the last fifteen years it has begun to be used for the computation of gravitational fields near black holes and to novel space-time spectral analysis.