Math @ Duke

Publications [#243720] of Anita T. Layton
Papers Published
 Layton, AT; Spotz, WF, A semiLagrangian double Fourier method for the shallow water equations on the sphere,
Journal of Computational Physics, vol. 189 no. 1
(2003),
pp. 180196 [doi]
(last updated on 2018/08/20)
Abstract: We describe a numerical method, based on the semiLagrangian semiimplicit 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 highorder 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 threetimelevel semiLagrangian semiimplicit 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.


dept@math.duke.edu
ph: 919.660.2800
fax: 919.660.2821
 
Mathematics Department
Duke University, Box 90320
Durham, NC 277080320

