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

作者：Eldred Christopher^*; Le Roux Daniel Y

来源：Journal of Computational Physics, 2018, 371: 779-800.

DOI：10.1016/j.jcp.2018.06.007

• Dispersion relationship
• Shallow water equations
• Geophysical fluid dynamics
• Mixed finite elements
• Finite element exterior calculus
• Mimetic Galerkin differences

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 th

• 出版日期2018-10-15
• 单位INRIA