A well‐balanced upstream flux‐splitting finite‐volume scheme for shallow‐water flow simulations with irregular bed topography

This study extends the upstream flux‐splitting finite‐volume (UFF) scheme to shallow water equations with source terms. Coupling the hydrostatic reconstruction method (HRM) with the UFF scheme achieves a resultant numerical scheme that adequately balances flux gradients and source terms. The proposed scheme is validated in three benchmark problems and applied to flood flows in the natural/irregular river with bridge pier obstructions. The results of the simulations are in satisfactory agreement with the available analytical solutions, experimental data and field measurements. Comparisons of the present results with those obtained by the surface gradient method (SGM) demonstrate the superior stability and higher accuracy of the HRM. The stability test results also show that the HRM requires less CPU time (up to 60%) than the SGM. The proposed well‐balanced UFF scheme is accurate, stable and efficient to solve flow problems involving irregular bed topography. Copyright © 2009 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.     

Semi‐Lagrangian semi‐implicit time‐splitting scheme for the shallow water equations

Time‐splitting technique applied in the context of the semi‐Lagrangian semi‐implicit method allows the use of extended time steps mainly based on physical considerations and reduces the number of numerical operations at each time step such that it is approximately proportional to the number of the points of spatial grid. To control time growth of the additional truncation errors, the standard stabilizing correction method is modified with no penalty for accuracy and efficiency of the algorithm. A linear analysis shows that constructed scheme is stable for time steps up to 2h. Numerical integrations with actual atmospheric fields of pressure and wind confirm computational efficiency, extended stability and accuracy of the proposed scheme. Copyright © 2006 John Wiley & Sons, Ltd.     

A dual‐weighted trust‐region adaptive POD 4‐D Var applied to a finite‐volume shallow water equations model on the sphere

In this paper we study solutions of an inverse problem for a global shallow water model controlling its initial conditions specified from the 40‐yr ECMWF Re‐analysis (ERA‐40) data sets, in the presence of full or incomplete observations being assimilated in a time interval (window of assimilation) with or without background error covariance terms. As an extension of the work by Chen et al. (Int. J. Numer. Meth. Fluids 2009), we attempt to obtain a reduced order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4D‐Var for a finite volume global shallow water equation model based on the Lin–Rood flux‐form semi‐Lagrangian semi‐implicit time integration scheme. Different approaches of POD implementation for the reduced inverse problem are compared, including a dual‐weighted method for snapshot selection coupled with a trust‐region POD adaptivity approach. Numerical results with various observational densities and background error covariance operator are also presented. The POD 4‐D Var model results combined with the trust‐region adaptivity exhibit similarity in terms of various error metrics to the full 4D Var results, but are obtained using a significantly lesser number of minimization iterations and require lesser CPU time. Based on our previous and current work, we conclude that POD 4‐D Var certainly warrants further studies, with promising potential of its extension to operational 3‐D numerical weather prediction models. Copyright © 2011 John Wiley & Sons, Ltd.     

Hybrid finite‐volume finite‐difference scheme for the solution of Boussinesq equations

A hybrid scheme composed of finite‐volume and finite‐difference methods is introduced for the solution of the Boussinesq equations. While the finite‐volume method with a Riemann solver is applied to the conservative part of the equations, the higher‐order Boussinesq terms are discretized using the finite‐difference scheme. Fourth‐order accuracy in space for the finite‐volume solution is achieved using the MUSCL‐TVD scheme. Within this, four limiters have been tested, of which van‐Leer limiter is found to be the most suitable. The Adams–Basforth third‐order predictor and Adams–Moulton fourth‐order corrector methods are used to obtain fourth‐order accuracy in time. A recently introduced surface gradient technique is employed for the treatment of the bottom slope. A new model ‘HYWAVE’, based on this hybrid solution, has been applied to a number of wave propagation examples, most of which are taken from previous studies. Examples include sinusoidal waves and bi‐chromatic wave propagation in deep water, sinusoidal wave propagation in shallow water and sinusoidal wave propagation from deep to shallow water demonstrating the linear shoaling properties of the model. Finally, sinusoidal wave propagation over a bar is simulated. The results are in good agreement with the theoretical expectations and published experimental results. Copyright © 2005 John Wiley & Sons, Ltd.     

Performance of finite volume solutions to the shallow water equations with shock‐capturing schemes

Numerical methods have become well established as tools for solving problems in hydraulic engineering. In recent years the finite volume method (FVM) with shock capturing capabilities has come to the fore because of its suitability for modelling a variety of types of flow; subcritical and supercritical; steady and unsteady; continuous and discontinuous and its ability to handle complex topography easily. This paper is an assessment and comparison of the performance of finite volume solutions to the shallow water equations with the Riemann solvers; the Osher, HLL, HLLC, flux difference splitting (Roe) and flux vector splitting. In this paper implementation of the FVM including the Riemann solvers, slope limiters and methods used for achieving second order accuracy are described explicitly step by step. The performance of the numerical methods has been investigated by applying them to a number of examples from the literature, providing both comparison of the schemes with each other and with published results. The assessment of each method is based on five criteria; ease of implementation, accuracy, applicability, numerical stability and simulation time. Finally, results, discussion, conclusions and recommendations for further work are presented. Copyright © 2002 John Wiley & Sons, Ltd.     

Adaptive Q‐tree Godunov‐type scheme for shallow water equations

This paper presents details of a second‐order accurate, Godunov‐type numerical model of the two‐dimensional shallow water equations (SWEs) written in matrix form and discretized using finite volumes. Roe's flux function is used for the convection terms and a non‐linear limiter is applied to prevent unwanted spurious oscillations. A new mathematical formulation is presented, which inherently balances flux gradient and source terms. It is, therefore, suitable for cases where the bathymetry is non‐uniform, unlike other formulations given in the literature based on Roe's approximate Riemann solver. The model is based on hierarchical quadtree (Q‐tree) grids, which adapt to inherent flow parameters, such as magnitude of the free surface gradient and depth‐averaged vorticity. Validation tests include wind‐induced circulation in a dish‐shaped basin, two‐dimensional frictionless rectangular and circular dam‐breaks, an oblique hydraulic jump, and jet‐forced flow in a circular reservoir. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite volume shock‐capturing solver of the fully coupled shallow water‐sediment equations

This paper describes a numerical solver of well‐balanced, 2D depth‐averaged shallow water‐sediment equations. The equations permit variable horizontal fluid density and are designed to model water‐sediment flow over a mobile bed. A Godunov‐type, Harten–Lax–van Leer contact (HLLC) finite volume scheme is used to solve the fully coupled system of hyperbolic conservation laws that describe flow hydrodynamics, suspended sediment transport, bedload transport and bed morphological change. Dependent variables are specially selected to handle the presence of the variable density property in the mathematical formulation. The model is verified against analytical and semi‐analytical solutions for bedload transport and suspended sediment transport, respectively. The well‐balanced property of the equations is verified for a variable‐density dam break flow over discontinuous bathymetry. Simulations of an idealised dam‐break flow over an erodible bed are in excellent agreement with previously published results, validating the ability of the model to capture the complex interaction between rapidly varying flow and an erodible bed and validating the eigenstructure of the system of variable‐density governing equations. Flow hydrodynamics and final bed topography of a laboratory‐based 2D partial dam breach over a mobile bed are satisfactorily reproduced by the numerical model. Comparison of the final bed topographies, computed for two distinct sediment transport methods, highlights the sensitivity of shallow water‐sediment models to the choice of closure relationships. Copyright © 2016 John Wiley & Sons, Ltd.     

Solution of the 2‐D shallow water equations with source terms in surface elevation splitting form

A vertex‐centred finite‐volume/finite‐element method (FV/FEM) is developed for solving 2‐D shallow water equations (SWEs) with source terms written in a surface elevation splitting form, which balances the flux gradients and source terms. The method is implemented on unstructured grids and the numerical scheme is based on a second‐order MUSCL‐like upwind Godunov FV discretization for inviscid fluxes and a classical Galerkin FE discretization for the viscous gradients and source terms. The main advantages are: (1) the discretization of SWE written in surface elevation splitting form satisfies the exact conservation property (??‐Property) naturally; (2) the simple centred‐type discretization can be used for the source terms; (3) the method is suitable for both steady and unsteady shallow water problems; and (4) complex topography can be handled based on unstructured grids. The accuracy of the method was verified for both steady and unsteady problems, including discontinuous cases. The results indicate that the new method is accurate, simple, and robust. Copyright © 2007 John Wiley & Sons, Ltd.     

An unconditionally stable semi‐Lagrangian method for the spherical atmospherical shallow water equations

A semi‐implicit, semi‐Lagrangian, mixed finite difference–finite volume model for the shallow water equations on a rotating sphere is introduced and discussed. Its main features are the vectorial treatment of the momentum equation and the finite volume approach for the continuity equation. Pressure and Coriolis terms in the momentum equation and velocity in the continuity equation are treated semi‐implicitly. Moreover, a splitting technique is introduced to preserve symmetry of the numerical scheme. An alternative asymmetric scheme (without splitting) is also introduced and the efficiency of both is discussed. The model is shown to be conservative in geopotential height and unconditionally stable for 0.5≤θ≤1. Numerical experiments on two standard test problems confirm the performance of the model. Copyright © 2000 John Wiley & Sons, Ltd.     

A mesh patching method for finite volume modelling of shallow water flow

A new mesh‐patching model is presented for shallow water flow described by the 2D non‐linear shallow water (NLSW) equations. The mesh‐patching model is based on AMAZON, a high‐resolution NLSW engine with an improved HLLC approximate Riemann solver. A new patching algorithm has been developed, which not only provides improved spatial resolution of flow features in particular parts of the mesh, but also simplifies and speeds up the (structured) grid generation process for an area with complicated geometry. The new patching technique is also compatible with increasingly popular parallel computing and adaptive grid techniques. The patching algorithm has been tested with moving bores, and results of test problems are presented and compared to previous work. Copyright © 2005 John Wiley & Sons, Ltd.     

A CIP/multi‐moment finite volume method for shallow water equations with source terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume‐integrated average (VIA) for each mesh cell, the surface‐integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi‐Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux‐based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non‐oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright © 2007 John Wiley & Sons, Ltd.     

A parametric finite‐difference method for shallow sea waves

This paper presents a parametric finite‐difference scheme concerning the numerical solution of the one‐dimensional Boussinesq‐type set of equations, as they were introduced by Peregrine (J. Fluid Mech. 1967; 27 (4)) in the case of waves relatively long with small amplitudes in water of varying depth. The proposed method, which can be considered as a generalization of the Crank‐Nickolson method, aims to investigate alternative approaches in order to improve the accuracy of analogous methods known from bibliography. The resulting linear finite‐difference scheme, which is analysed for stability using the Fourier method, has been applied successfully to a problem used by Beji and Battjes (Coastal Eng. 1994; 23 : 1–16), giving numerical results which are in good agreement with the corresponding results given by MIKE 21 BW (User Guide. In: MIKE 21, Wave Modelling, User Guide. 2002; 271–392) developed by DHI Software. Copyright © 2006 John Wiley & Sons, Ltd.     

A finite volume solver for 1D shallow‐water equations applied to an actual river

This paper describes the numerical solution of the 1D shallow‐water equations by a finite volume scheme based on the Roe solver. In the first part, the 1D shallow‐water equations are presented. These equations model the free‐surface flows in a river. This set of equations is widely used for applications: dam‐break waves, reservoir emptying, flooding, etc. The main feature of these equations is the presence of a non‐conservative term in the momentum equation in the case of an actual river. In order to apply schemes well adapted to conservative equations, this term is split in two terms: a conservative one which is kept on the left‐hand side of the equation of momentum and the non‐conservative part is introduced as a source term on the right‐hand side. In the second section, we describe the scheme based on a Roe Solver for the homogeneous problem. Next, the numerical treatment of the source term which is the essential point of the numerical modelisation is described. The source term is split in two components: one is upwinded and the other is treated according to a centred discretization. By using this method for the discretization of the source term, one gets the right behaviour for steady flow. Finally, in the last part, the problem of validation is tackled. Most of the numerical tests have been defined for a working group about dam‐break wave simulation. A real dam‐break wave simulation will be shown. Copyright © 2002 John Wiley & Sons, Ltd.     

A well‐balanced stable generalized Riemann problem scheme for shallow water equations using adaptive moving unstructured triangular meshes

We propose a well‐balanced stable generalized Riemann problem (GRP) scheme for the shallow water equations with irregular bottom topography based on moving, adaptive, unstructured, triangular meshes. In order to stabilize the computations near equilibria, we use the Rankine–Hugoniot condition to remove a singularity from the GRP solver. Moreover, we develop a remapping onto the new mesh (after grid movement) based on equilibrium variables. This, together with the already established techniques, guarantees the well‐balancing. Numerical tests show the accuracy, efficiency, and robustness of the GRP moving mesh method: lake at rest solutions are preserved even when the underlying mesh is moving (e.g., mesh points are moved to regions of steep gradients), and various comparisons with fixed coarse and fine meshes demonstrate high resolution at relatively low cost. Copyright © 2013 John Wiley & Sons, Ltd.     

A structured but non‐uniform Cartesian grid‐based model for the shallow water equations

In large‐scale shallow flow simulations, local high‐resolution predictions are often required in order to reduce the computational cost without losing the accuracy of the solution. This is normally achieved by solving the governing equations on grids refined only to those areas of interest. Grids with varying resolution can be generated by different approaches, e.g. nesting methods, patching algorithms and adaptive unstructured or quadtree gridding techniques. This work presents a new structured but non‐uniform Cartesian grid system as an alternative to the existing approaches to provide local high‐resolution mesh. On generating a structured but non‐uniform Cartesian grid, the whole computational domain is first discretized using a coarse background grid. Local refinement is then achieved by directly allocating a specific subdivision level to each background grid cell. The neighbour information is specified by simple mathematical relationships and no explicit storage is needed. Hence, the structured property of the uniform grid is maintained. After employing some simple interpolation formulae, the governing shallow water equations are solved using a second‐order finite volume Godunov‐type scheme in a similar way as that on a uniform grid. Copyright © 2010 John Wiley & Sons, Ltd.     

High resolution flux‐difference‐splitting scheme on adaptive grid for open‐channel flows

The effectiveness and usefulness of further enhancing the shock resolution of a second‐order accurate scheme for open‐channel flows by using an adaptive grid is investigated. The flux‐difference‐splitting (FDS) scheme based on the Lax–Wendroff numerical flux is implemented on a fixed as well as on a self‐adjusting grid for this purpose. The grid‐adjusting procedure, developed by Harten and Hyman, adjusts the grid by averaging the local characteristic velocities with respect to the signal amplitude in such a way that a shock always lies on a mesh point. This enables a scheme capable of perfectly resolving a stationary shock to capture a shock that moves from mesh point to mesh point. The Roe's approximate Jacobian is used for conservation and consistency, while theoretically sound treatment for satisfying entropy inequality conditions ensures physically realistic solutions. Details about inclusion of source terms, often left out of analyses for the homogeneous part of governing equations, are also explained. The numerical results for some exacting problems are compared with analytical as well as experimental results for examining improvements in resolution of discontinuities by the adaptive grid. Copyright © 2001 John Wiley & Sons, Ltd.     

Unstructured finite volume discretization of two‐dimensional depth‐averaged shallow water equations with porosity

This paper deals with the numerical discretization of two‐dimensional depth‐averaged models with porosity. The equations solved by these models are similar to the classic shallow water equations, but include additional terms to account for the effect of small‐scale impervious obstructions which are not resolved by the numerical mesh because their size is smaller or similar to the average mesh size. These small‐scale obstructions diminish the available storage volume on a given region, reduce the effective cross section for the water to flow, and increase the head losses due to additional drag forces and turbulence. In shallow water models with porosity these effects are modelled introducing an effective porosity parameter in the mass and momentum conservation equations, and including an additional drag source term in the momentum equations. This paper presents and compares two different numerical discretizations for the two‐dimensional shallow water equations with porosity, both of them are high‐order schemes. The numerical schemes proposed are well‐balanced, in the sense that they preserve naturally the exact hydrostatic solution without the need of high‐order corrections in the source terms. At the same time they are able to deal accurately with regions of zero porosity, where the water cannot flow. Several numerical test cases are used in order to verify the properties of the discretization schemes proposed. Copyright © 2009 John Wiley & Sons, Ltd.     

Solution of the shallow‐water equations using an adaptive moving mesh method

This paper extends an adaptive moving mesh method to multi‐dimensional shallow water equations (SWE) with source terms. The algorithm is composed of two independent parts: the SWEs evolution and the mesh redistribution. The first part is a high‐resolution kinetic flux‐vector splitting (KFVS) method combined with the surface gradient method for initial data reconstruction, and the second part is based on an iteration procedure. In each iteration, meshes are first redistributed by a variational principle and then the underlying numerical solutions are updated by a conservative‐interpolation formula on the resulting new mesh. Several test problems in one‐ and two‐dimensions with a general geometry are computed using the proposed moving mesh algorithm. The computations demonstrate that the algorithm is efficient for solving problems with bore waves and their interactions. The solutions with higher resolution can be obtained by using a KFVS scheme for the SWEs with a much smaller number of grid points than the uniform mesh approach, although we do not treat technically the bed slope source terms in order to balance the source terms and flux gradients. Copyright © 2004 John Wiley & Sons, Ltd.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.