A Nodal Triangle-based Spectral Element Method for the Shallow Water Equations on the Sphere
Abstract
A nodal triangle-based spectral element method for the shallow water equations on the sphere is presented. The original SE method uses quadrilateral elements and high-order nodal Lagrange polynomials, constructed from a tensor-product of the Legendre-Gauss-Lobatto points. In this work we construct the high-order Lagrange polynomials directly on the triangle using nodal sets obtained from the electrostatics principle and Fekete points. These points have good approximation properties and far better Lebesgue constants than any other nodal set derived for the triangle. By employing triangular elements as the basic building blocks of the SE method and the Cartesian coordinate form of the equations, we can use any grid imaginable including adaptive unstructured grids. Results for six test cases are presented to confirm the accuracy and stability of the method. The results show that the triangle-based SE method yields the expected exponential convergence and that it can be more accurate than the quadrilateral-based SE method even while using 30% to 60% fewer grid points especially when adaptive grids are used to align the grid with the flow direction. However, at the moment, the quadrilateral-based SE method is twice as fast as the triangle-based SE method because the latter does not yield a diagonal mass matrix.
Document Details
- Document Type
- Technical Report
- Publication Date
- Jan 01, 2005
- Accession Number
- ADA511198
Entities
People
- F. X. Giraldo
- T. Warburton
Organizations
- United States Naval Research Laboratory