A dual‐weighted trust‐region adaptive POD 4D‐VAR applied to a finite‐element shallow‐water equations model

We consider a limited‐area finite‐element discretization of the shallow‐water equations model. Our purpose in this paper is to solve an inverse problem for the above model controlling its initial conditions in presence of observations being assimilated in a time interval (window of assimilation). We then attempt to obtain a reduced‐order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4‐D VAR. Different approaches of POD implementation of the reduced inverse problem are compared, including a dual‐weighed method for snapshot selection coupled with a trust‐region POD approach. Numerical results obtained point to an improved accuracy in all metrics tested when dual‐weighing choice of snapshots is combined with POD adaptivity of the trust‐region type. Results of ad‐hoc adaptivity of the POD 4‐D VAR turn out to yield less accurate results than trust‐region POD when compared with high‐fidelity model. Directions of future research are finally outlined. Copyright © 2009 John Wiley & Sons, Ltd.     

An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations

An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd.     

An improved algorithm for the shallow water equations model reduction: Dynamic Mode Decomposition vs POD

We propose an improved framework for dynamic mode decomposition (DMD) of 2‐D flows for problems originating from meteorology when a large time step acts like a filter in obtaining the significant Koopman modes, therefore, the classic DMD method is not effective. This study is motivated by the need to further clarify the connection between Koopman modes and proper orthogonal decomposition (POD) dynamic modes. We apply DMD and POD to derive reduced order models (ROM) of the shallow water equations. Key innovations for the DMD‐based ROM introduced in this paper are the use of the Moore–Penrose pseudoinverse in the DMD computation that produced an accurate result and a novel selection method for the DMD modes and associated amplitudes and Ritz values. A quantitative comparison of the spatial modes computed from the two decompositions is performed, and a rigorous error analysis for the ROM models obtained is presented. Copyright © 2015 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.     

Comparison of POD reduced order strategies for the nonlinear 2D shallow water equations

This paper introduces tensorial calculus techniques in the framework of POD to reduce the computational complexity of the reduced nonlinear terms. The resulting method, named tensorial POD, can be applied to polynomial nonlinearities of any degree p. Such nonlinear terms have an online complexity of , where k is the dimension of POD basis and therefore is independent of full space dimension. However, it is efficient only for quadratic nonlinear terms because for higher nonlinearities, POD model proves to be less time consuming once the POD basis dimension k is increased. Numerical experiments are carried out with a two‐dimensional SWE test problem to compare the performance of tensorial POD, POD, and POD/discrete empirical interpolation method (DEIM). Numerical results show that tensorial POD decreases by 76× the computational cost of the online stage of POD model for configurations using more than 300,000 model variables. The tensorial POD SWE model was only 2 to 8× slower than the POD/DEIM SWE model but the implementation effort is considerably increased. Tensorial calculus was again employed to construct a new algorithm allowing POD/DEIM SWE model to compute its offline stage faster than POD and tensorial POD approaches. Copyright © 2014 John Wiley & Sons, Ltd.     

Development and testing of a simple 2D finite volume model of sub‐critical shallow water flow

Many environmental applications of shallow water flow modelling can be characterized as only slowly varying and everywhere sub‐critical. A simplified finite volume model is therefore developed that is capable of describing pertinent shallow water flow processes more efficiently than the usual Godunov/ Riemann characteristics approaches. The model is tested against a number of analytical and numerical solutions to the governing equations. The model reproduces accurately flow round a circular bend, flow over topography, flow up an initially dry beach and floodwave propagation down a meandering river reach, with mass conservative solutions. Copyright © 2004 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 study of reduced‐order 4DVAR with a finite element shallow water model

Four‐dimensional variational data assimilation (4DVAR) is frequently used to improve model forecasting skills. This method improves a model consistency with available data by minimizing a cost function measuring the model–data misfit with respect to some model inputs and parameters. Associated with this type of method, however, are difficulties related to the coding of the adjoint model, which is needed to compute the gradient of the 4DVAR cost function. Proper orthogonal decomposition (POD) is a model reduction method that can be used to approximate the gradient calculation in 4DVAR. In this work, two ways of using POD in 4DVAR are presented, namely model‐reduced 4DVAR and reduced adjoint 4DVAR (RA‐4DVAR). Both techniques employ POD to obtain a reduced‐order approximation of the forward linear tangent operator. The difference between the two methods lies in the treatment of the forward model. Model‐reduced 4DVAR performs minimization entirely in the POD‐reduced space, thereby achieving very low computational costs, but sacrificing accuracy of the end result. On the other hand, the RA‐4DVAR uses POD to approximate only the adjoint model. The main contribution of this study is a comparative performance analysis of these 4DVAR methodologies on a nonlinear finite element shallow water model. The sensitivity of the methods to perturbations in observations and the number of observation points is examined. The results from twin experiments suggest that the RA‐4DVAR method is easy to implement and computationally efficient and provides a robust approach for achieving reasonable results in the context of variational data assimilation. Copyright © 2015 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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Fourth order finite volume solution to shallow water equations and applications

This study presents the fourth order accurate finite volume solution to shallow water equations. Fourth order accuracy in space was provided by using the Monotone Upstream‐centered Schemes for Conservation Laws–Total Variation Diminishing scheme, whereas fourth order accurate solution in time was achieved by using the third order predictor scheme of Adams–Basforth followed by the fourth order corrector scheme of Adams–Moulton. The applicability and accuracy of the solution algorithm were explored on complex flow conditions. These flow conditions cover a theoretical well‐known partial two‐dimensional dam break problems and an experimental flow in a compound channel with or without a bridge. The applicability limits of the solution algorithm were discussed. The overall performance of the solution was found to be reasonably good. Copyright © 2013 John Wiley & Sons, Ltd.     

An unstructured grid,three‐dimensional model based on the shallow water equations

A semi‐implicit finite difference model based on the three‐dimensional shallow water equations is modified to use unstructured grids. There are obvious advantages in using unstructured grids in problems with a complicated geometry. In this development, the concept of unstructured orthogonal grids is introduced and applied to this model. The governing differential equations are discretized by means of a semi‐implicit algorithm that is robust, stable and very efficient. The resulting model is relatively simple, conserves mass, can fit complicated boundaries and yet is sufficiently flexible to permit local mesh refinements in areas of interest. Moreover, the simulation of the flooding and drying is included in a natural and straightforward manner. These features are illustrated by a test case for studies of convergence rates and by examples of flooding on a river plain and flow in a shallow estuary. Copyright © 2000 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 finite element/volume method model of the depth‐averaged horizontally 2D shallow water equations

Analysis of surface water flows is of central importance in understanding and predicting a wide range of water engineering issues. Dynamics of surface water is reasonably well described using the shallow water equations (SWEs) with the hydrostatic pressure assumption. The SWEs are nonlinear hyperbolic partial differential equations that are in general required to be solved numerically. Application of a simple and efficient numerical model is desirable for solving the SWEs in practical problems. This study develops a new numerical model of the depth‐averaged horizontally 2D SWEs referred to as 2D finite element/volume method (2D FEVM) model. The continuity equation is solved with the conforming, standard Galerkin FEM scheme and momentum equations with an upwind, cell‐centered finite volume method scheme, utilizing the water surface elevation and the line discharges as unknowns aligned in a staggered manner. The 2D FEVM model relies on neither Riemann solvers nor high‐resolution algorithms in order to serve as a simple numerical model. Water at a rest state is exactly preserved in the model. A fully explicit temporal integration is achieved in the model using an efficient approximate matrix inversion method. A series of test problems, containing three benchmark problems and three experiments of transcritical flows, are carried out to assess accuracy and versatility of the model. Copyright © 2014 John Wiley & Sons, Ltd.