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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:371

Article

Eldred, Christopher¹  Le Roux, Daniel Y.² 

[1] Univ Grenoble Alpes, Inria, CNRS, Grenoble INP,LJK, F-38000 Grenoble, France

[2] Univ Lyon 1, Univ Lyon, CNRS, Inst Camille Jordan, 43,Blvd 11 Novembre 1918, F-69622 Villeurbanne, France

关键词: Dispersion relationship;    Shallow water equations;    Geophysical fluid dynamics;    Mixed finite elements;    Finite element exterior calculus;    Mimetic Galerkin differences;   

DOI  :  10.1016/j.jcp.2018.06.007

来源: Elsevier

【 摘 要 】

In this work, we study the dispersion properties of two compatible Galerkin schemes for the 1D linearized shallow water equations: the P-n(C) - P-n-1(DG) and the GD(n) - DGD(n-1) element pairs. P-n is the order n Lagrange space, P-n-1(DG) is the order n - 1 discontinuous Lagrange space, GD(n) is the order n Galerkin difference space, and DGD(n-1) is the order n - 1 discontinuous Galerkin difference space. 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-2(C) - P-1(DG) pair has indeed indicated the presence of spectral gaps, and it is extended in this paper to the study of the P-n(C) - P-n-1(DG) pair for arbitrary n. Additionally, an alternative element pair is introduced, the GD(n) - DGD(n-1) 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 diagonalizes the velocity mass matrix) for the P-n(C) - P-n-1(DG) pair and reduced quadrature for the GD(n) - DGD(n-1) 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-3(C) - P-2(DG) 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. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free