Department of Mathematics
 Search | Help | Login | pdf version | printable version

Math @ Duke



Publications [#243720] of Anita T. Layton

Papers Published

  1. Layton, AT; Spotz, WF, A semi-Lagrangian double Fourier method for the shallow water equations on the sphere, Journal of Computational Physics, vol. 189 no. 1 (2003), pp. 180-196 [doi]
    (last updated on 2017/11/24)

    We describe a numerical method, based on the semi-Lagrangian semi-implicit approach, for solving the shallow water equations (SWEs) in spherical coordinates. The most popular spatial discretization method used in global atmospheric models is currently the spectral transform method, which generates high-order numerical solutions and provides an elegant solution to the pole problems induced by a spherical coordinate system. However, the standard spherical harmonic spectral transform method requires associated Legendre transforms, which for problems with resolutions of current interest, have a computational complexity of O(N3), where N is the number of spatial subintervals in one dimension. Thus, the double Fourier spectral method, which uses trigonometric series, may be a viable alternative. The advantage of the double Fourier method is that fast Fourier transforms, which have a computational complexity of O(N2 log N), can be used in both the longitudinal and latitudinal directions. In this implementation, the SWEs are discretized in time by means of the three-time-level semi-Lagrangian semi-implicit method, which integrates along fluid trajectories and allows large time steps while maintaining stability. Numerical results for the standard SWEs test suite are presented to demonstrate the stability and accuracy of the method. © 2003 Elsevier Science B.V. All rights reserved.
ph: 919.660.2800
fax: 919.660.2821

Mathematics Department
Duke University, Box 90320
Durham, NC 27708-0320