RAPPORT Spatial Discretization of the Shallow Water Equations in Spherical Geometry using Osher ’ s Scheme

The shallow water equations in spherical geometry provide a first prototype for developing and testing numerical algorithms for atmospheric circulation models. Since the seventies these models are often solved with spectral methods. Increasing demands on grid resolution combined with massive parallelism and local grid refinement seem to offer significantly better perspectives for gridpoint methods. In this paper we study the use of Osher’s finite-volume scheme for the spatial discretization of the shallow water equations on the rotating sphere. This high-order finite volume scheme of upwind type is well suited to solve a hyperbolic system of equations. Special attention is paid to the pole problem. To that end Osher’s scheme is applied on the common (reduced) latitude-longitude grid and on a stereographic grid. The latter is most appropriate in the polar region as in stereographic coordinates the pole singularity does not exist. The latitude-longitude grid is preferred on lower latitudes. Therefore, across the sphere we apply Osher’s scheme on a combined grid connecting the two grids at high latitude. We will show that this provides an attractive spatial discretization for explicit integration methods, as it can greatly reduce the time step limitation incurred by the pole singularity when using a latitude-longitude grid only. When time step limitation plays no significant role, the standard (reduced) latitude-longitude grid is advocated provided that the grid is kept sufficiently fine in the polar region to resolve flow over the poles. 1991 Mathematics Subject Classification: 65M06, 86-08, 86A10. 1991 ACM Computing Classification System: G.1.8, J.2.

[1]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme III. Upstream-centered finite-difference schemes for ideal compressible flow , 1977 .

[2]  The Shallow Water Test Cases for a Global Model with Documentation of the Results , .

[3]  G. L. Browning,et al.  A comparison of three numerical methods for solving differential equations on the sphere , 1989 .

[4]  G. Starius,et al.  Composite mesh difference methods for elliptic boundary value problems , 1977 .

[5]  S. Orszag.  Transform method for the calculation of vector-coupled sums: Application to the spectral form of the vorticity equation , 1970 .

[6]  S. Osher,et al.  Stable and entropy satisfying approximations for transonic flow calculations , 1980 .

[7]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[8]  S. Spekreijse,et al.  Multigrid Solution of the Steady Euler Equations. , 1989 .

[9]  David L. Williamson,et al.  Review of Numerical Approaches for Modeling Global Transport , 1992 .

[10]  N. A. Phillips,et al.  A Map Projection System Suitable for Large-scale Numerical Weather Prediction , 1957 .

[11]  Willem Hundsdorfer,et al.  Vectorization aspects of a spherical advection scheme on a reduced grid , 1994 .

[12]  S. Osher,et al.  One-sided difference approximations for nonlinear conservation laws , 1981 .

[13]  J. Lions,et al.  New formulations of the primitive equations of atmosphere and applications , 1992 .

[14]  J. Holton.  An introduction to dynamic meteorology , 2004 .

[15]  Stanley Osher,et al.  Upwind schemes and boundary conditions with applications to Euler equations in general geometries , 1983 .

[16]  A. Arakawa,et al.  Numerical methods used in atmospheric models , 1976 .

[17]  S. P. Spekreijse,et al.  Multiple grid and Osher''s scheme for the efficient solution of the steady Euler equations Applied N , 1986 .

[18]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[19]  David L. Williamson,et al.  Comparison of Grids and Difference Approximations for Numerical Weather Prediction Over a Sphere , 1973 .

[20]  G. J. Haltiner.  Numerical Prediction and Dynamic Meteorology , 1980 .