Space-time discontinuous Galerkin discretization of rotating shallow water equations on moving grids

A space-time discontinuous Galerkin (DG) discretization is presented for the (rotating) shallow water equations over varying topography. We formulate the space-time DG finite element discretization in an efficient and conservative discretization. The HLLC flux is used as numerical flux through the finite element boundaries. When discontinuities are present, we locally apply dissipation around these discontinuities with the help of Krivodonova's discontinuity indicator such that spurious oscillations are suppressed. The non-linear algebraic system resulting from the discretization is solved using a pseudo-time integration with a second-order five-stage Runge-Kutta method. A thorough verification of the space-time DG finite element method is undertaken by comparing numerical and exact solutions. We also carry out a discrete Fourier analysis of the one dimensional linear rotating shallow water equations to show that the method is unconditionally stable with minimal dispersion and dissipation error. The numerical scheme is validated in a novel way by considering various simulations of bore-vortex interactions in combination with a qualitative analysis of PV generation by non-uniform bores. Finally, the space-time DG method is particularly suited for problems where dynamic grid motion is required. To demonstrate this we simulate waves generated by a wave maker and verify these for low amplitude waves where linear theory is approximately valid

Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes

We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of spherical symmetry within a 3+1 splitting of spacetime. As a direct application, we have constructed a one-dimensional code, EDGES, which has been used to asses the viability of these methods via a series of tests involving highly relativistic flows in strong gravity. Our results show that discontinuous Galerkin methods are able not only to handle strong relativistic shock waves but, at the same time, to attain very high orders of accuracy and exponential convergence rates in smooth regions of the flow. Given these promising prospects and their affinity with a pseudospectral solution of the Einstein equations, discontinuous Galerkin methods could represent a new paradigm for the accurate numerical modelling in relativistic astrophysics.Comment: 24 pages, 19 figures. Small changes; matches version to appear in PR

Mixed finite elements for numerical weather prediction

We show how two-dimensional mixed finite element methods that satisfy the conditions of finite element exterior calculus can be used for the horizontal discretisation of dynamical cores for numerical weather prediction on pseudo-uniform grids. This family of mixed finite element methods can be thought of in the numerical weather prediction context as a generalisation of the popular polygonal C-grid finite difference methods. There are a few major advantages: the mixed finite element methods do not require an orthogonal grid, and they allow a degree of flexibility that can be exploited to ensure an appropriate ratio between the velocity and pressure degrees of freedom so as to avoid spurious mode branches in the numerical dispersion relation. These methods preserve several properties of the C-grid method when applied to linear barotropic wave propagation, namely: a) energy conservation, b) mass conservation, c) no spurious pressure modes, and d) steady geostrophic modes on the $f$-plane. We explain how these properties are preserved, and describe two examples that can be used on pseudo-uniform grids: the recently-developed modified RT0-Q0 element pair on quadrilaterals and the BDFM1-\pdg element pair on triangles. All of these mixed finite element methods have an exact 2:1 ratio of velocity degrees of freedom to pressure degrees of freedom. Finally we illustrate the properties with some numerical examples.Comment: Revision after referee comment

A limiter-based well-balanced discontinuous Galerkin method for shallow-water flows with wetting and drying: Triangular grids

A novel wetting and drying treatment for second-order Runge-Kutta discontinuous Galerkin (RKDG2) methods solving the non-linear shallow water equations is proposed. It is developed for general conforming two-dimensional triangular meshes and utilizes a slope limiting strategy to accurately model inundation. The method features a non-destructive limiter, which concurrently meets the requirements for linear stability and wetting and drying. It further combines existing approaches for positivity preservation and well-balancing with an innovative velocity-based limiting of the momentum. This limiting controls spurious velocities in the vicinity of the wet/dry interface. It leads to a computationally stable and robust scheme -- even on unstructured grids -- and allows for large time steps in combination with explicit time integrators. The scheme comprises only one free parameter, to which it is not sensitive in terms of stability. A number of numerical test cases, ranging from analytical tests to near-realistic laboratory benchmarks, demonstrate the performance of the method for inundation applications. In particular, super-linear convergence, mass-conservation, well-balancedness, and stability are verified