This paper describes the efficient and accurate solution of the twodimensional anelastic equations by a Fourier-Chebyshev spectral method. A fourth-order Runge-Kutta method is used for th6 time integration, with the diffusion terms treated implicitly and all other terms (including the pressure gradi
A multi-layer spectral model and the semi-implicit method
โ Scribed by B. J. Hoskins; A. J. Simmons
- Publisher
- John Wiley and Sons
- Year
- 1975
- Tongue
- English
- Weight
- 956 KB
- Volume
- 101
- Category
- Article
- ISSN
- 0035-9009
No coin nor oath required. For personal study only.
โฆ Synopsis
Abstract
The formulaton of a multiโlayer primitive equation model on the sphere is described. The horizontal representation is by means of spherical harmonics, truncated either in the triangular or rhomboidal manner. The time integration is performed using the semiโimplicit method in which the linearized gravity wave terms are time averaged and thus the fast moving waves of this type are slowed. For a 5โlayer hemispheric model with triangular truncation at wavenumber 21, storage of 38K words is needed and with the time scheme allowing a timeโstep of 90 minutes, one day's simulation requires 11 seconds of CDC 7600 time. The growth of a baroclinic wave on a simple basic state of differential solid body rotation is exhibited. The errors involved in this case in utilizing the large timeโstep allowed by the semiโimplicit scheme are thoroughly examined by comparing wave amplitudes and phases, conservation properties and gravity wave treatment for different timeโsteps. These errors are found to be negligible. The conservation properties of the model are in fact extremely good. The vertical finite differencing scheme of Arakawa is studied in the same baroclinic instability simulation. The growth is similar though the conservation of angular momentum is greatly improved. The transform method used in all these integrations allows some aliasing, but this is shown to be negligible.
๐ SIMILAR VOLUMES
Classical semi-implicit backward Euler/Adams-Bashforth time discretizations of the Navier-Stokes equations induce, for high-Reynolds number flows, severe restrictions on the time step. Such restrictions can be relaxed by using semi-Lagrangian schemes essentially based on splitting the full problem i