Abstract:The numerical approximation of the shallow-water equations, which support geophysical wave propagation, namely the fast external Poincaré and Kelvin waves and the slow large-scale planetary Rossby waves, is a delicate and difficult problem. Indeed, the coupling between the momentum and the continuity equations may lead to the presence of erratic or spurious solutions for these waves, e.g. the spurious pressure and inertial modes. In addition, the presence of spurious branches and the emergence of spectral gaps in the dispersion relation at specific wavenumbers may lead to anomalous dissipation/dispersion in the representation of both fast and slow waves. The aim of the present study is to propose a class of possible discretization schemes, via the discontinuous Galerkin method, that is not affected by the above-mentioned problems. A Fourier/stability analysis of the 2D shallow-water model is performed by analysing the finite volume method and the linear discontinuous P1DG and non conforming P1NC Galerkin discretizations. These are stabilized by means of a family of numerical fluxes via the Polynomial Viscosity Matrix approach. A long time stability result is proven for all schemes and fluxes examined here. Further, a super-convergent result is demonstrated for the P1NC discrete frequencies compared to the P1DG ones, except for the slow mode in the Roe flux case. Indeed, we show that the Roe flux yields spurious frequencies and sub-optimal rates of convergence for the slow mode, for both the P1DG and P1NC methods. Finally, numerical solutions of linear and non-linear test problems confirm the theoretical results and reveal the computational cost efficiency of the P1NC approach.

Frozen Gaussian approximation for high frequency wave propagation in periodic media

Convergence of frozen Gaussian approximation for high frequency wave propagation

Frozen Gaussian approximation for the fractional Schrödinger equation

Frozen Gaussian Sampling for Scalar Wave Equations

Efficient Calculation of Large Finite Periodic Structures Based on Surface Wave Analysis

Frozen Gaussian Series Representation of the Imaginary Time Propagator Theory and Numerical Tests

On the Bloch decomposition based spectral method for wave propagation in periodic media

Seismic modeling using the frozen Gaussian approximation

Gaussian beam approach for the boundary value problem of high frequency Helmholtz equation

Wavelet-Galerkin Method for Computations of Electromagnetic Fields

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Global Geometrical Optics Method

An Error Bound for the Time-Sliced Thawed Gaussian Propagation Method

A convergent low-wavenumber, high-frequency homogenization of the wave equation in periodic media with a source term

Global Geometrical Optics Approximation to the High Frequency Helmholtz Equation with Discontinuous Media

An FFT-based algorithm for efficient computation of Green's functions for the Helmholtz and Maxwell's equations in periodic domains

A super-convergent quasi-discontinuous Galerkin method for approximating inertia-gravity and Rossby waves in geophysical flows

Computing High Frequency Solutions of Symmetric Hyperbolic Systems with Polarized Waves

Phase-Informed Discontinuous Galerkin Method for Extremely High-Frequency Wave Modeling

Finite-Temperature Correlation Functions Obtained from Combined Real- and Imaginary-Time Propagation of Variational Thawed Gaussian Wavepackets

A hybrid numerical method for elastic wave propagation in discontinuous media with complex geometry