A discontinuous Galerkin method for the shallow water equations in spherical triangular coordinates


Contact
Matthias.Laeuter [ at ] awi.de

Abstract

A global model of the atmosphere is presented governed by the shallow water equations and discretized by a Runge-Kutta discontinuous Galerkin method on an unstructured triangular grid. The shallow water equations on the sphere, a two-dimensional surface in R3, are locally represented in terms of spherical triangular coordinates, the appropriate local coordinate mappings on triangles. On every triangular grid element, this leads to a two-dimensional representation of tangential momentum and therefore only two discrete momentum equations. The discontinuous Galerkin method consists of an integral formulation which requires both area (elements) and line (element faces) integrals. Here, we use a Rusanov numerical flux to resolve the discontinuous fluxes at the element faces. A strong stability-preserving third-order Runge-Kutta method is applied for the time discretization. The polynomial space of order k on each curved triangle of the grid is characterized by a Lagrange basis and requires high-order quadature rules for the integration over elements and element faces. For the presented method no mass matrix inversion is necessary, except in a preprocessing step. The validation of the atmospheric model has been done considering standard tests from Williamson et al. [D.L. Williamson, J.B. Drake, J.J. Hack, R. Jakob, P.N. Swarztrauber, A standard test set for numerical approximations to the shallow water equations in spherical geometry, J. Comput. Phys. 102 (1992) 211-224], unsteady analytical solutions of the nonlinear shallow water equations and a barotropic instability caused by an initial perturbation of a jet stream. A convergence rate of O (Δ xk + 1) was observed in the model experiments. Furthermore, a numerical experiment is presented, for which the third-order time-integration method limits the model error. Thus, the time step Δ t is restricted by both the CFL-condition and accuracy demands. Conservation of mass was shown up to machine precision and energy conservation converges for both increasing grid resolution and increasing polynomial order k. © 2008 Elsevier Inc.



Item Type
Article
Authors
Divisions
Programs
Publication Status
Published
Eprint ID
17879
DOI https://www.doi.org/10.1016/j.jcp.2008.08.019

Cite as
Läuter, M. , Giraldo, F. X. , Handorf, D. and Dethloff, K. (2008): A discontinuous Galerkin method for the shallow water equations in spherical triangular coordinates , Journal of Computational Physics, 227 (24), pp. 10226-10242 . doi: https://www.doi.org/10.1016/j.jcp.2008.08.019


Download
[thumbnail of Fulltext]
Preview
PDF (Fulltext)
Lut2007b.pdf

Download (860kB) | Preview
Cite this document as:

Share


Citation

Research Platforms
N/A

Campaigns
N/A


Actions
Edit Item Edit Item