Dispersion analysis of compatible Galerkin schemes for the 1D shallow water model

Christopher Eldred 1 Daniel Le Roux 2
1 AIRSEA - Mathematics and computing applied to oceanic and atmospheric flows
Grenoble INP - Institut polytechnique de Grenoble - Grenoble Institute of Technology, UGA - Université Grenoble Alpes, LJK - Laboratoire Jean Kuntzmann, Inria Grenoble - Rhône-Alpes
2 MMCS - Modélisation mathématique, calcul scientifique
ICJ - Institut Camille Jordan [Villeurbanne]
Abstract : In this work, we study the dispersion properties of two compatible Galerkin schemes for the 1D linearized shallow water equations: the P C n −P DG n−1 and the GD n −DGD n element pairs. Compatible Galerkin methods have many desirable properties, including energy conservation, steady geostrophic modes and the absence of spurious stationary modes, such as pressure modes. However, this does not guarantee good wave dispersion properties. Previous work on the P C 2 − P DG 1 pair has indeed indicated the presence of spectral gaps, and it is extended in this paper to the study of the P C n − P DG n−1 pair for arbitrary n. Additionally, an alternative element pair is introduced, the GD n − DGD n pair, that is free of spectral gaps while benefiting from the desirable properties of compatible elements. Asymptotic convergence rates are established for both element pairs, including the use of inexact quadrature (which di-agonalizes the velocity mass matrix) for the P C n − P DG n−1 pair and reduced quadrature for the GD n − DGD n pair. Plots of the dispersion relationship and group velocities for a wide range of n and Rossby radii are shown. A brief investigation into the utility of mass lumping to remove the spectral gaps for the P C 3 − P DG 2 pair is performed. Finally, a pair of numerical simulations are run to investigate the consequences of the spectral gaps and highlight the main differences between the two elements.
Document type :
Journal articles
Complete list of metadatas

Cited literature [27 references]  Display  Hide  Download

Contributor : Christopher Eldred <>
Submitted on : Wednesday, June 6, 2018 - 10:41:52 AM
Last modification on : Thursday, November 21, 2019 - 2:01:31 AM


Files produced by the author(s)



Christopher Eldred, Daniel Le Roux. Dispersion analysis of compatible Galerkin schemes for the 1D shallow water model. Journal of Computational Physics, Elsevier, 2018, 371, pp.779-800. ⟨10.1016/j.jcp.2018.06.007⟩. ⟨hal-01669048v2⟩



Record views


Files downloads