Practical evaluation of five partly discontinuous finite element pairs for the non‐conservative shallow water equations

This paper provides a comparison of five finite element pairs for the shallow water equations. We consider continuous, discontinuous and partially discontinuous finite element formulations that are supposed to provide second‐order spatial accuracy. All of them rely on the same weak formulation, using Riemann solver to evaluate interface integrals. We define several asymptotic limit cases of the shallow water equations within their space of parameters. The idea is to develop a comparison of these numerical schemes in several relevant regimes of the subcritical shallow water flow. Finally, a new pair, using non‐conforming linear elements for both velocities and elevation (P?P), is presented, giving optimal rates of convergence in all test cases. P?P₁ and P?P₁ mixed formulations lack convergence for inviscid flows. P?P₂ pair is more expensive but provides accurate results for all benchmarks. P?P provides an efficient option, except for inviscid Coriolis‐dominated flows, where a small lack of convergence is observed. Copyright © 2009 John Wiley & Sons, Ltd.     

A triangular spectral/hp discontinuous Galerkin method for modelling 2D shallow water equations

We present a spectral/hp element discontinuous Galerkin model for simulating shallow water flows on unstructured triangular meshes. The model uses an orthogonal modal expansion basis of arbitrary order for the spatial discretization and a third‐order Runge–Kutta scheme to advance in time. The local elements are coupled together by numerical fluxes, evaluated using the HLLC Riemann solver. We apply the model to test cases involving smooth flows and demonstrate the exponentially fast convergence with regard to polynomial order. We also illustrate that even for results of ‘engineering accuracy’ the computational efficiency increases with increasing order of the model and time of integration. The model is found to be robust in the presence of shocks where Gibbs oscillations can be suppressed by slope limiting. Copyright 2004 John Wiley & Sons, Ltd.     

Edge‐based finite element method for shallow water equations

This paper describes an edge‐based implementation of the generalized residual minimum (GMRES) solver for the fully coupled solution of non‐linear systems arising from finite element discretization of shallow water equations (SWEs). The gain in terms of memory, floating point operations and indirect addressing is quantified for semi‐discrete and space–time analyses. Stabilized formulations, including Petrov–Galerkin models and discontinuity‐capturing operators, are also discussed for both types of discretization. Results illustrating the quality of the stabilized solutions and the advantages of using the edge‐based approach are presented at the end of the paper. Copyright © 2001 John Wiley & Sons, Ltd.     

High‐order semi‐implicit time‐integrators for a triangular discontinuous Galerkin oceanic shallow water model

We extend the explicit in time high‐order triangular discontinuous Galerkin (DG) method to semi‐implicit (SI) and then apply the algorithm to the two‐dimensional oceanic shallow water equations; we implement high‐order SI time‐integrators using the backward difference formulas from orders one to six. The reason for changing the time‐integration method from explicit to SI is that explicit methods require a very small time step in order to maintain stability, especially for high‐order DG methods. Changing the time‐integration method to SI allows one to circumvent the stability criterion due to the gravity waves, which for most shallow water applications are the fastest waves in the system (the exception being supercritical flow where the Froude number is greater than one). The challenge of constructing a SI method for a DG model is that the DG machinery requires not only the standard finite element‐type area integrals, but also the finite volume‐type boundary integrals as well. These boundary integrals pose the biggest challenge in a SI discretization because they require the construction of a Riemann solver that is the true linear representation of the nonlinear Riemann problem; if this condition is not satisfied then the resulting numerical method will not be consistent with the continuous equations. In this paper we couple the SI time‐integrators with the DG method while maintaining most of the usual attributes associated with DG methods such as: high‐order accuracy (in both space and time), parallel efficiency, excellent stability, and conservation. The only property lost is that of a compact communication stencil typical of time‐explicit DG methods; implicit methods will always require a much larger communication stencil. We apply the new high‐order SI DG method to the shallow water equations and show results for many standard test cases of oceanic interest such as: standing, Kelvin and Rossby soliton waves, and the Stommel problem. The results show that the new high‐order SI DG model, that has already been shown to yield exponentially convergent solutions in space for smooth problems, results in a more efficient model than its explicit counterpart. Furthermore, for those problems where the spatial resolution is sufficiently high compared with the length scales of the flow, the capacity to use high‐order (HO) time‐integrators is a necessary complement to the employment of HO space discretizations, since the total numerical error would be otherwise dominated by the time discretization error. In fact, in the limit of increasing spatial resolution, it makes little sense to use HO spatial discretizations coupled with low‐order time discretizations. Published in 2009 by John Wiley & Sons, Ltd.     

Two‐dimensional dispersion analyses of finite element approximations to the shallow water equations

Dispersion analysis of discrete solutions to the shallow water equations has been extensively used as a tool to define the relationships between frequency and wave number and to determine if an algorithm leads to a dual wave number response and near 2Δx oscillations. In this paper, we explore the application of two‐dimensional dispersion analysis to cluster based and Galerkin finite element‐based discretizations of the primitive shallow water equations and the generalized wave continuity equation (GWCE) reformulation of the harmonic shallow water equations on a number of grid configurations. It is demonstrated that for various algorithms and grid configurations, contradictions exist between the results of one‐dimensional and two‐dimensional dispersion analysis as a result of subtle changes in the mass matrix. Numerical experiments indicate that the two‐dimensional dispersion analysis correctly predicts the existence and onset of near 2Δx noise in the solution. Copyright © 2004 John Wiley & Sons, Ltd.     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION （Ⅰ）——THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations, silt transport equation, the equation of bottom topography change, and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element (MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived. The error estimates are optimal.     

An adaptive discretization of shallow‐water equations based on discontinuous Galerkin methods

In this paper, we present a discontinuous Galerkin formulation of the shallow‐water equations. An orthogonal basis is used for the spatial discretization and an explicit Runge–Kutta scheme is used for time discretization. Some results of second‐order anisotropic adaptive calculations are presented for dam breaking problems. The adaptive procedure uses an error indicator that concentrates the computational effort near discontinuities like hydraulic jumps. Copyright © 2006 John Wiley & Sons, Ltd.     

An eigenvector‐based linear reconstruction approach for time stepping in discontinuous Galerkin scheme used to solve shallow water equations

Discontinuous Galerkin (DG) methods have shown promising results for solving the two‐dimensional shallow water equations. In this paper, the classical Runge–Kutta (RK) time discretisation is replaced by the eigenvector‐based reconstruction (EVR) that allows the second‐order time accuracy to be achieved within a single time‐stepping procedure. Moreover, the EVRDG approach yields stable solutions near drying and wetting fronts, whereas the classical RKDG approach yields instabilities. The proposed EVRDG technique is compared with the original RKDG approach on various test cases with analytical solutions. The EVRDG solutions are shown to be as accurate as those obtained with the RKDG scheme. Besides, the EVRDG scheme is 1.6 times faster than the RKDG method. Simulating dambreaks involving dry beds confirms that EVRDG scheme gives correct solutions, whereas the RKDG method yields instabilities. Copyright © 2012 John Wiley & Sons, Ltd.     

The shallow flow equations solved on adaptive quadtree grids

This paper describes an adaptive quadtree grid‐based solver of the depth‐averaged shallow water equations. The model is designed to approximate flows in complicated large‐scale shallow domains while focusing on important smaller‐scale localized flow features. Quadtree grids are created automatically by recursive subdivision of a rectangle about discretized boundary, bathymetric or flow‐related seeding points. It can be fitted in a fractal‐like sense by local grid refinement to any boundary, however distorted, provided absolute convergence to the boundary is not required and a low level of stepped boundary can be tolerated. Grid information is stored as a tree data structure, with a novel indexing system used to link information on the quadtree to a finite volume discretization of the governing equations. As the flow field develops, the grids may be adapted using a parameter based on vorticity and grid cell size. The numerical model is validated using standard benchmark tests, including seiches, Coriolis‐induced set‐up, jet‐forced flow in a circular reservoir, and wetting and drying. Wind‐induced flow in the Nichupté Lagoon, México, provides an illustrative example of an application to flow in extremely complicated multi‐connected regions. Copyright © 2001 John Wiley & Sons, Ltd.     

Discontinuous boundary implementation for the shallow water equations

Quasi‐bubble finite element approximations to the shallow water equations are investigated focusing on implementations of the surface elevation boundary condition. We first demonstrate by numerical results that the conventional implementation of the boundary condition degrades the accuracy of the velocity solution. It is also shown that the degraded velocity leads to a critical instability if the advection term is present in the momentum equation. Then we propose an alternative implementation for the boundary condition. We refer to this alternative implementation as a discontinuous boundary (DB) implementation because it introduces at each boundary node two independent mass–flux values that result in a discontinuity at the boundary. Numerical results show that the proposed DB implementation is consistent, stabilizes the quasi‐bubble scheme, and leads to second‐order accuracy at the surface elevation specified boundary. Copyright © 2004 John Wiley & Sons, Ltd.     

A well‐balanced explicit/semi‐implicit finite element scheme for shallow water equations in drying–wetting areas

Numerical solutions of the shallow water equations can be used to reproduce flow hydrodynamics occurring in a wide range of regions. In hydraulic engineering, the objectives include the prediction of dam break wave propagation, fluvial floods and other catastrophic flooding phenomena, the modeling of estuarine and coastal circulations, and the design and optimization of hydraulic structures. In this paper, a well‐balanced explicit and semi‐implicit finite element scheme for shallow water equations over complex domains involving wetting and drying is proposed. The governing equations are discretized by a fractional finite element method using a two‐step Taylor–Galerkin scheme. First, the intermediate increment of conserved variable is obtained explicitly neglecting the pressure gradient term. This is then corrected for the effects of pressure once the pressure increment has been obtained from the Poisson equation. In order to maintain the ‘well‐balanced’ property, the pressure gradient term and bed slope terms are incorporated into the Poisson equation. Moreover, a local bed slope modification technique is employed in drying–wetting interface treatments. The proposed model is well validated against several theoretical benchmark tests. Copyright © 2014 John Wiley & Sons, Ltd.     

An hp‐adaptive discontinuous Galerkin method for shallow water flows

An adaptive spectral/hp discontinuous Galerkin method for the two‐dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non‐conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p?1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h‐type refinement, the parent element is subdivided into four similar sibling elements. The time‐stepping is performed using a third‐order Runge–Kutta scheme. The performance of the hp‐adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p‐adaptivity is more efficient than h‐adaptivity with respect to degrees of freedom and computational time. Copyright © 2010 John Wiley & Sons, Ltd.     

Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes

A 2D, depth-integrated, free surface flow solver for the shallow water equations is developed and tested. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov-type second-order upwind finite volume formulation, whereby the inviscid fluxes of the system of equations are obtained using Roe's flux function. The eigensystem of the 2D shallow water equations is derived and is used for the construction of Roe's matrix on an unstructured mesh. The viscous terms of the shallow water equations are computed using a finite volume formulation which is second-order-accurate. Verification of the solution technique for the inviscid form of the governing equations as well as for the full system of equations is carried out by comparing the model output with documented published results and very good agreement is obtained. A numerical experiment is also conducted in order to evaluate the performance of the solution technique as applied to linear convection problems. The presented results show that the solution technique is robust. © 1997 John Wiley & Sons, Ltd.     

Adaptive finite element method for analysis of pollutant dispersion in shallow water

A finite element method for analysis of pollutant dispersion in shallow water is presented. The analysis is divided into two parts ： （ 1 ） computation of the velocity flow field and water surface elevation, and （2） computation of the pollutant concentration field from the dispersion model. The method was combined with an adaptive meshing technique to increase the solution accuracy, as well as to reduce the computational time and computer memory. The finite element formulation and the computer programs were validated by several examples that have known solutions. In addition, the capability of the combined method was demonstrated by analyzing pollutant dispersion in Chao Phraya River near the gulf of Thailand.     

Stability/dispersion analysis of the discontinuous Galerkin linearized shallow‐water system

The frequency or dispersion relation for the discontinuous Galerkin mixed formulation of the 1‐D linearized shallow‐water equations is analysed, using several basic DG mixed schemes. The dispersion properties are compared analytically and graphically with those of the mixed continuous Galerkin formulation for piecewise‐linear bases on co‐located grids. Unlike the Galerkin case, the DG scheme does not exhibit spurious stationary pressure modes. However, spurious propagating modes have been identified in all the present discontinuous Galerkin formulations. Numerical solutions of a test problem to simulate fast gravity modes illustrate the theoretical results and confirm the presence of spurious propagating modes in the DG schemes. Copyright © 2005 John Wiley & Sons, Ltd.     

Compatibility between finite volumes and finite elements using solutions of shallow water equations for substance transport

This paper formulates a finite volume analogue of a finite element schematization of three‐dimensional shallow water equations. The resulting finite volume schematization, when applied to the continuity equation, exactly reproduces the set of matrix equations that is obtained by the application of the corresponding finite element schematization to the continuity equation. The procedure allows the consistent and mass conserving coupling of the finite element Telemac model for three‐dimensional flow with the finite volume Delft3D‐WAQ model for water quality. The work has been carried out as part of a joint development by LNHE and WL∣Delft Hydraulics to explore the mutual interaction of their software. Copyright © 2006 John Wiley & Sons, Ltd.     

Mesh adaptation strategies for shallow water flow

This paper shows how the mesh adaptation technique can be exploited for the numerical simulation of shallow water flow. The shallow water equations are numerically approximated by the Galerkin finite element method, using linear elements for the elevation field and quadratic elements for the unit width discharge field; the time advancing scheme is of a fractional step type. The standard mesh refinement technique is coupled with the numerical solver; movement and elimination of nodes of the initial triangulation is not allowed. Two error indicators are discussed and applied in the numerical examples. The conclusion focuses the relevant advantages obtained by applying this adaptive approach by considering specific test cases of steady and unsteady flows. Copyright © 1999 John Wiley & Sons, Ltd.     

A flux‐splitting solver for shallow water equations with source terms

This paper introduces a stable flux‐splitting solver for one‐dimensional (1D) shallow water equations. This solver is specifically designed to satisfy a strengthened consistency condition for stationary solutions that ensures the stability and accuracy of the scheme. It applies to channels with variable depth and width, including terms modelling friction at bottom and vertical walls. Some numerical tests by comparison to both analytical solutions and experimental measurements show the good performances of the scheme. Copyright © 2003 John Wiley & Sons, Ltd.     

Stabilized finite element formulation for incompressible flow on distorted meshes

Flow computations frequently require unfavourably meshes, as for example highly stretched elements in regions of boundary layers or distorted elements in deforming arbitrary Lagrangian Eulerian meshes. Thus, the performance of a flow solver on such meshes is of great interest. The behaviour of finite elements with residual‐based stabilization for incompressible Newtonian flow on distorted meshes is considered here. We investigate the influence of the stabilization terms on the results obtained on distorted meshes by a number of numerical studies. The effect of different element length definitions within the elemental stabilization parameter is considered. Further, different variants of residual‐based stabilization are compared indicating that dropping the second derivatives from the stabilization operator, i.e. using a streamline upwind Petrov–Galerkin type of formulation yields better results in a variety of cases. A comparison of the performance of linear and quadratic elements reveals further that the inconsistency of linear elements equipped with residual‐based stabilization introduces significant errors on distorted meshes, while quadratic elements are almost unaffected by moderate mesh distortion. Copyright © 2008 John Wiley & Sons, Ltd.     

Improvements in mass conservation using alternative boundary implementations for a quasi‐bubble finite element shallow water model

Finite element approaches generally do not guarantee exact satisfaction of conservation laws especially when Dirichlet‐type boundary conditions are imposed. This article discusses improvement of the global mass conservation property of quasi‐bubble finite element solutions for the shallow water equations, focusing on implementations of the surface‐elevation boundary conditions. We propose two alternative implementations, which are shown by numerical verification to be effective in improving the smoothness of solutions near the boundary and in reducing the mass conservation error. The improvement of the mass conservation property contributes to augmenting the reliability and robustness of long‐term time integrations. Copyright © 2006 John Wiley & Sons, Ltd.