Abstract
Spectral elements combine the accuracy and exponential convergence of conventional spectral methods with the geometric flexibility of finite elements. Additionally, there are several apparent computational advantages to using spectral element methods on microprocessors. In particular, the computations are naturally cache-blocked and derivatives may be computed using nearest neighbor communications. Thus, an explicit spectral element atmospheric model has demonstrated close to linear scaling on a variety of distributed memory computers including the IBM SP and Linux Clusters. Explicit formulations of PDE's arising in geophysical fluid dynamics, such as the primitive equations on the sphere, are time-step limited by the phase speed of gravity waves. Semi-implicit time integration schemes remove the stability restriction but require the solution of an elliptic BVP. By employing a weak formulation of the governing equations, it is possible to obtain a symmetric Helmholtz operator that permits the solution of the implicit problem using conjugate gradients. We find that a block-Jacobi preconditioned conjugate gradient solver accelerates the simulation rate of the semi-implicit relative to the explicit formulation for practical climate resolutions by about a factor of three.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
REFERENCES
Bernardi, C., and Maday, Y. (1992). Approximations Spectrales de Problèmes aux Limites Elliptiques, Mathematiques et Applications, Vol. 10, Springer-Verlag, Paris, France, p. 242.
D'Azevedo, E., Eijkhout, V., and Romine, C. (1992). Reducing communication costs in the conjugate gradient algorithm on distributed memory multiprocessors. LAPACK working note 56, University of Tennessee.
Fischer, P. F. (1997). An overlapping Schwarz method for spectral element solution of the incompressible Navier-Stokes equations. J. Comp. Phys. 133, 84–101.
Held, I. H., and Suarez, M. J. (1994). A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models. Bull. Amer. Met. Soc. 75, 1825–1830.
Iskandarani, M., Haidvogel, D. B., and Boyd, J. P. (1995). A staggered spectral element model with application to the oceanic shallow water equations. Int. J. Numer. Meth. Fluids 20, 394–414.
Karniadakis, G. M., and Sherwin, S. J. (1999). Spectral/hp Element Methods for CFD, Oxford University Press, Oxford, England, p. 390.
Patera, A. T. (1984). A spectral element method for fluid dynamics: Laminar flow in a channel expansion. J. Comp. Phys. 54, 468.
Rancic, M., Purser, R. J., and Mesinger, F. (1996). A global shallow-water model using an expanded spherical cube: Gnomic versus conformal coordinates. Q. J. R. Meteorol. Soc. 122, 959–982.
Rivier, L., Polvani, L. M., and Loft, R. (2000). An efficient spectral general circulation model for distributed memory computers. To appear.
Robert, A. (1969). The integration of a spectral model of the atmosphere by the implicit method. In Proceedings of WMO/I UGG Symposium on NWP, VII, Tokyo, Japan, pp. 19–24.
Ronquist, E. M. (1988). Optimal Spectral Element Methods for the Unsteady Three Dimensional Navier Stokes Equations, Ph.D Thesis, Massachusetts Institute of Technology, p. 176.
Ronquist, E. M. (1991). A domain decomposition method for elliptic boundary value problems: Application to unsteady incompressible fluid flow. In Proceedings of the Fifth International Symposium on Domain Decomposition Methods for Partial Differential Equations, SIAM, Philadelphia, pp. 545–557.
Taylor, M., Tribbia, J., and Iskandarani, M. (1997a). The spectral element method for the shallow water equations on the sphere. J. Comp. Phys. 130, 92–108.
Taylor, M., Loft, R., and Tribbia, J. (1997b). Performance of a spectral element atmospheric model (SEAM) on the HP Exemplar SPP2000. NCAR Technical Note 439 + EDD.
Williamson, D. L., Drake, J. B., Hack, J. J., Jakob, R., and Swarztrauber, P. N. (1992). A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comp. Phys. 102, 211–224.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Thomas, S.J., Loft, R.D. Parallel Semi-Implicit Spectral Element Methods for Atmospheric General Circulation Models. Journal of Scientific Computing 15, 499–518 (2000). https://doi.org/10.1023/A:1011188832645
Issue Date:
DOI: https://doi.org/10.1023/A:1011188832645