Space-time discontinuous Galerkin finite element method for shallow water flows

A space-time discontinuous Galerkin (DG) finite element method is presented for the shallow water equations over varying bottom topography. The method results in non-linear equations per element, which are solved locally by establishing the element communication with a numerical HLLC flux. To deal with spurious oscillations around discontinuities, we employ a dissipation operator only around discontinuities using Krivodonova's discontinuity detector. The numerical scheme is verified by comparing numerical and exact solutions, and validated against a laboratory experiment involving flow through a contraction. We conclude that the method is second order accurate in both space and time for linear polynomials.\u

A well-balanced Runge--Kutta Discontinuous Galerkin method for the Shallow-Water Equations with flooding and drying

We build and analyze a Runge--Kutta Discontinuous Galerkin method to approximate the one- and two-dimensional Shallow-Water Equations. We introduce a flux modification technique to derive a well-balanced scheme preserving steady-states at rest with variable bathymetry and a slope modification technique to deal satisfactorily with flooding and drying. Numerical results illustrating the performance of the proposed scheme are presented

A water wave model with horizontal circulation and accurate dispersion

We describe a new water wave model which is variational, and combines a depth-averaged vertical (component of) vorticity with depth-dependent potential flow. The model facilitates the further restriction of the vertical profile of the velocity potential to n-th order polynomials or a finite element profile with a small number of elements (say), leading to a framework for efficient modelling of the interaction of steepening and breaking waves near the shore with a large-scale horizontal flow. The equations are derived from a constrained variational formulation which leads to conservation laws for energy, mass, momentum and vertical vorticity (or circulation). We show that the potential flow water wave equations and the shallow-water equations are recovered in the relevant limits, and provide approximate shock relations for the model which can be used in numerical schemes to model breaking waves

Discontinuous Galerkin scheme for the spherical shallow water equations with applications to tsunami modeling and prediction

We present a novel high-order discontinuous Galerkin discretization for the spherical shallow water equations, able to handle wetting/drying and non-conforming, curved meshes in a well-balanced manner. This requires a well-balanced discretization, that cannot rely on exact quadrature, due to the curved mesh. Using the strong form of the discontinuous Galerkin discretization, we achieve a splitting of the well-balanced condition into individual problems for the flux and volume terms, which has significant advantages: It allows for the construction of non-conforming, well-balanced flux discretizations, i.e. we can perform non- conforming mesh refinement while preserving the well-balanced property of the scheme. More importantly, this approach enables the development of a new method for handling wet/dry transitions. In contrast to other wetting/drying methods, it is well-balanced and able to handle wetting/drying robustly at any polynomial order, without the introduction of physical model assumptions such as viscosity, artificial porosity or cancellation of gravity. We perform a series of one-dimensional tests and analyze the properties of our scheme. In order to validate our method for the simulation of large-scale tsunami events on the rotating sphere, we perform numerical simulations of the 2011 Tohoku tsunami and compare our results to real-world buoy data. The method is able to predict arrival times and wave amplitudes accurately even over long distances. This indicates that our method accurately captures all physical phenomena relevant to the long-term evolution of tsunami waves

Variational water-wave model with accurate dispersion and vertical vorticity

A new water-wave model has been derived which is based on variational techniques and combines a depth-averaged vertical (component of) vorticity with depth-dependent potential flow. The model facilitates the further restriction of the vertical profile of the velocity potential to n-th order polynomials or a finite-element profile with a small number of elements (say), leading to a framework for efficient modeling of the interaction of steepening and breaking waves near the shore with a large-scale horizontal flow. The equations are derived from a constrained variational formulation which leads to conservation laws for energy, mass, momentum and vertical vorticity. It is shown that the potential-flow water-wave equations and the shallow-water equations are recovered in the relevant limits. Approximate shock relations are provided, which can be used in numerical schemes to model breaking waves

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

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

An adaptive discontinuous Galerkin method for the simulation of hurricane storm surge

Numerical simulations based on solving the 2D shallow water equations using a discontinuous Galerkin (DG) discretisation have evolved to be a viable tool for many geophysical applications. In the context of flood modelling, however, they have not yet been methodologically studied to a large extent. Systematic model testing is non-trivial as no comprehensive collection of numerical test cases exists to ensure the correctness of the implementation. Hence, the first part of this manuscript aims at collecting test cases from the literature that are generally useful for storm surge modellers and can be used to benchmark codes. On geographic scale, hurricane storm surge can be interpreted as a localised phenomenon making it ideally suited for adaptive mesh refinement (AMR). Past studies employing dynamic AMR have exclusively focused on nested meshes. For that reason, we have developed a DG storm surge model on a triangular and dynamically adaptive mesh. In order to increase computational efficiency, the refinement is driven by physics-based refinement indicators capturing major model sensitivities. Using idealised numerical test cases, we demonstrate the model’s ability to correctly represent all source terms and reproduce known variability of coastal flooding with respect to hurricane characteristics such as size and approach speed. Finally, the adaptive mesh significantly reduces computing time with no effect on storm waves measured at discrete wave gauges just off the coast which shows the model’s potential for use as a robust simulation tool for real-time predictions

A well-balanced Runge--Kutta Discontinuous Galerkin method for the Shallow-Water Equations with flooding and drying

We build and analyze a Runge--Kutta Discontinuous Galerkin method to approximate the one- and two-dimensional Shallow-Water Equations. We introduce a flux modification technique to derive a well-balanced scheme preserving steady-states at rest with variable bathymetry and a slope modification technique to deal satisfactorily with flooding and drying. Numerical results illustrating the performance of the proposed scheme are presented

Simulation of a viscous fluid spreading by a bidimensional shallow water model

In this paper we propose a numerical method to solve the Cauchy problem based on the viscous shallow water equations in an horizontally moving domain. More precisely, we are interested in a flooding and drying model, used to modelize the overflow of a river or the intrusion of a tsunami on ground. We use a non conservative form of the two-dimensional shallow water equations, in eight velocity formulation and we build a numerical approximation, based on the Arbitrary Lagrangian Eulerian formulation, in order to compute the solution in the moving domain