Spectral methods for the simulation of incompressible flows in spherical shells

SUMMARY A spatial discretization of the incompressible Navier‐Stokes equation is presented in which the velocity is decomposed using poloidal and toroidal scalars whose spatial dependence is given in terms of spherical harmonics and Chebychev polynomials. The radial resolution needs to be large enough at any given angular resolution in order to avoid instability in the simulation of rotating flows. Several semi-implicit time steps are discussed. The most accurate scheme is an integrating factor technique. Copyright © 1999 John Wiley & Sons, Ltd.