A spectral element semi-Lagrangian (SESL) method for the spherical shallow water equations
Perot, J. B.
Fischer, P. F.
MetadataShow full item record
A spectral element semi-Lagrangian (SESL) method for the shallow water equations on the sphere is presented. The sphere is discretized using a hexahedral grid although any grid imaginable can be used as long as it is comprised of quadrilaterals. The equations are written in Cartesian coordinates to eliminate the pole singularity which plagues the equations in spherical coordinates. In a previous paper [Int. J. Numer. Methods Fluids 35 (2001) 869] we showed how to construct an explicit Eulerian spectral element (SE) model on the sphere; we now extend this work to a semi-Lagrangian formulation. The novelty of the Lagrangian formulation presented is that the high order SE basis functions are used as the interpolation functions for evaluating the values at the Lagrangian departure points. This makes the method not only high order accurate but quite general and thus applicable to unstructured grids and portable to distributed memory computers. The equations are discretized fully implicitly in time in order to avoid having to interpolate derivatives at departure points. By incorporating the Coriolis terms into the Lagrangian derivative, the block LU decomposition of the equations results in a symmetric positive-definite pseudo-Helmholtz operator which we solve using the generalized minimum residual method (GMRES) with a fast projection method [Comput. Methods Appl. Mech. Eng. 163 (1998) 193]. Results for eight test cases are presented to confirm the accuracy and stability of the method. These results show that SESL yields the same accuracy as an Eulerian spectral element semi-implicit (SESI) while allowing for time-steps 10 times as large and being up to 70% more efficient.
The article of record as published may be located at http://dx.doi.org/10.1016/S0021-9991(03)00300-0
Showing items related by title, author, creator and subject.
Giraldo, F.X. (2001);The spectral element method for the two-dimensional shallow water equations on the sphere is presented. The equations are written in conservation form and the domains are discretized using quadrilateral elements obtained ...
Giraldo, F.X. (2000);The weak Lagrange-Galerkin finite element method fur the 2D shallow water equations on the sphere is presented. This method offers stable and accurate solutions because the equations are integrated along thr characteristics. ...
Giraldo, F.X.; Rosmond, T. E. (2004);A new dynamical core for numerical weather prediction (NWP) based on the spectral element method is presented. This paper represents a departure from previously published work on solving the atmospheric primitive equations ...