Sensitivity of the 1D shallow water equations with source terms: Solution method for discontinuous flows

A finite volume‐based numerical technique is presented concerning the sensitivity of the solution of the one‐dimensional Shallow Water Equations with scalar transport. An approximate Riemann solver is proposed for direct sensitivity calculation even in the presence of discontinuous solutions. The Shallow Water Sensitivity Equations are first derived as well as the expressions of the sensitivity source terms, initial and boundary conditions. The numerical technique is then detailed and application examples are provided to assess the method's efficiency in estimating the sensitivity to different parameters (friction coefficient and initial and boundary conditions). The application of the dam‐break problem to a trapezoidal channel is also provided. The comparison with the analytical solution and the classical empirical approach illustrates the usefulness of the direct sensitivity calculation. Copyright © 2010 John Wiley & Sons, Ltd.     

Solution of shallow water equations using fully adaptive multiscale schemes

The concept of fully adaptive multiscale finite volume methods has been developed to increase spatial resolution and to reduce computational costs of numerical simulations. Here grid adaptation is performed by means of a multiscale analysis based on biorthogonal wavelets. In order to update the solution in time we use a local time stepping strategy that has been recently developed for hyperbolic conservation laws. The adaptive multiresolution scheme is now applied to two‐dimensional shallow water equations with source terms. The efficiency of the scheme is demonstrated on several problems with a general geometry, including circular damp breaks, oblique hydraulic jump, supercritical channel flows encountering sudden change in cross‐section, and, finally, the bore wave and its interactions. Copyright © 2005 John Wiley & Sons, Ltd.     

A numerical study of the motion of a wave running up a beach

The main difficulty for the numerical calculation of the wave running up a beach is the treatment of its moving water boundary. In this paper a scheme of turning the free boundary problem into a fixed boundary problem is designed. The calculated run-up height is consistent with the experiments. Some interesting wave phenomena are also found.     

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.     

Similarities between the quasi‐bubble and the generalized wave continuity equation solutions to the shallow water equations

Two common strategies for solving the shallow water equations in the finite element community are the generalized wave continuity equation (GWCE) reformulation and the quasi‐bubble velocity approximation. The GWCE approach has been widely analysed in the literature. In this work, the quasi‐bubble equations are analysed and comparisons are made between the quasi‐bubble approximation of the primitive form of the shallow water equations and a linear finite element approximation of the GWCE reformulation of the shallow water equations. The discrete condensed quasi‐bubble continuity equation is shown to be identical to a discrete wave equation for a specific GWCE weighting parameter value. The discrete momentum equations are slightly different due to the bubble function. In addition, the dispersion relationships are shown to be almost identical and numerical experiments confirm that the two schemes compute almost identical results. Analysis of the quasi‐bubble formulation suggests a relationship that may guide selection of the optimal GWCE weighting parameter. Copyright © 2004 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.     

High‐order compact finite difference schemes for the vorticity–divergence representation of the spherical shallow water equations

This work is devoted to the application of the super compact finite difference method (SCFDM) and the combined compact finite difference method (CCFDM) for spatial differencing of the spherical shallow water equations in terms of vorticity, divergence, and height. The fourth‐order compact, the sixth‐order and eighth‐order SCFDM, and the sixth‐order and eighth‐order CCFDM schemes are used for the spatial differencing. To advance the solution in time, a semi‐implicit Runge–Kutta method is used. In addition, to control the nonlinear instability, an eighth‐order compact spatial filter is employed. For the numerical solution of the elliptic equations in the problem, a direct hybrid method, which consists of a high‐order compact scheme for spatial differencing in the latitude coordinate and a fast Fourier transform in longitude coordinate, is utilized. The accuracy and convergence rate for all methods are verified against exact analytical solutions. Qualitative and quantitative assessments of the results for an unstable barotropic mid‐latitude zonal jet employed as an initial condition are addressed. It is revealed that the sixth‐order and eighth‐order CCFDMs and SCFDMs lead to a remarkable improvement of the solution over the fourth‐order compact method. Copyright © 2015 John Wiley & Sons, Ltd.     

Fourier analysis of a class of upwind schemes in shallow water systems for gravity and Rossby waves

A Fourier analysis has been performed for a class of upwind finite volume schemes, including the study of phase speed, group velocity, damping and dispersion. In the first part, pure gravity waves are investigated. As expected, most upwind schemes lead to a significant damping, but they exhibit a better phase behavior than most centered schemes. In the second part, the Coriolis parameter is considered and the Rossby modes are studied. In this case, all selected upwind schemes lead to a severe damping. The numerical results are also compared with those obtained by using a slope limiter approach. It is concluded that most upwind schemes with or without slope limiters present poor results for an accurate calculation of the Rossby modes. Copyright © 2008 John Wiley & Sons, Ltd.     

An implicit time‐marching algorithm for shallow water models based on the generalized wave continuity equation

Wave equation models currently discretize the generalized wave continuity equation with a three‐time‐level scheme centered at k and the momentum equation with a two‐time‐level scheme centered at k+1/2; non‐linear terms are evaluated explicitly. However in highly non‐linear applications, the algorithm becomes unstable at even moderate Courant numbers. This paper examines an implicit treatment of the non‐linear terms using an iterative time‐marching algorithm. Depending on the domain, results from one‐dimensional experiments show up to a tenfold increase in stability and temporal accuracy. The sensitivity of stability to variations in the G‐parameter (a numerical weighting parameter in the generalized wave continuity equation) was examined; results show that the greatest increase in stability occurs with G/τ=2–50. In the one‐dimensional experiments, three different types of node spacing techniques—constant, variable, and LTEA (Localized Truncation Error Analysis)—were examined; stability is positively correlated to the uniformity of the node spacing. Lastly, a scaling analysis demonstrates that the magnitudes of the non‐linear terms are positively correlated to those that most influence stability, particularly the term containing the G‐parameter. It is evident that the new algorithm improves stability and temporal accuracy in a cost‐effective manner. Copyright © 2001 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.     

A two‐dimensional numerical scheme of dry/wet fronts for the Saint‐Venant system of shallow water equations

We propose a new two‐dimensional numerical scheme to solve the Saint‐Venant system of shallow water equations in the presence of partially flooded cells. Our method is well balanced, positivity preserving, and handles dry states. The latter is ensured by using the draining time step technique in the time integration process, which guarantees non‐negative water depths. Unlike previous schemes, our technique does not generate high velocities at the dry/wet boundaries, which are responsible for small time step sizes and slow simulation runs. We prove that the new scheme preserves ‘lake at rest’ steady states and guarantees the positivity of the computed fluid depth in the partially flooded cells. We test the new scheme, along with another recent scheme from the literature, against the analytical solution for a parabolic basin and show the improved simulation performance of the new scheme for two real‐world scenarios. Copyright © 2014 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.     

Finite volume solution to integrated shallow surface–saturated groundwater flow

This paper describes development of an integrated shallow surface and saturated groundwater model (GSHAW5). The surface flow motion is described by the 2‐D shallow water equations and groundwater movement is described by the 2‐D groundwater equations. The numerical solution of these equations is based on the finite volume method where the surface water fluxes are estimated using the Roe shock‐capturing scheme, and the groundwater fluxes are computed by application of Darcy's law. Use of a shock‐capturing scheme ensures ability to simulate steady and unsteady, continuous and discontinuous, subcritical and supercritical surface water flow conditions. Ground and surface water interaction is achieved by the introduction of source‐sink terms into the continuity equations. Two solutions are tightly coupled in a single code. The numerical solutions and coupling algorithms are explained. The model has been applied to 1‐D and 2‐D test scenarios. The results have shown that the model can produce very accurate results and can be used for simulation of situations involving interaction between shallow surface and saturated groundwater flows. Copyright © 2005 John Wiley & Sons, Ltd.     

Third- and fourth-order well-balanced schemes for the shallow water equations based on the CWENO reconstruction

High-order finite volume schemes for conservation laws are very useful in applications, due to their ability to compute accurate solutions on quite coarse meshes and with very few restrictions on the kind of cells employed in the discretization. For balance laws, the ability to approximate up to machine precision relevant steady states allows the scheme to compute accurately, also on coarse meshes, small perturbations of such states, which are very relevant for many applications. In this paper, we propose third- and fourth-order accurate finite volume schemes for the shallow water equations. The schemes have the well-balanced property thanks to a path-conservative approach applied to an appropriate nonconservative reformulation of the equations. High-order accuracy is achieved by designing truly two-dimensional (2D) reconstruction procedures of the central WENO (CWENO ) type. The novel schemes are tested for accuracy and well-balancing and shown to maintain positivity of the water height on wet/dry transitions. Finally, they are applied to simulate the Tohoku 2011 tsunami event.     

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.     

The prediction of spatially periodic flows using a finite‐volume model

A periodic boundary condition has been developed that can be used in conjunction with a specified flow rate to produce accurate results in spatially periodic geometries. This condition is useful in situations where the flow rate is known, or more importantly, in cases where the pressure gradient is not known a priori, such as in countercurrent flows. Using the present condition, the flow rate is imposed at the inlet in terms of a bulk velocity, but the velocity field evolves as part of the solution. The condition is formulated to be suitable for both fixed and moving periodic domains. For the case of a moving domain, a correction is introduced to account for changes in the instantaneous velocity through the periodic edges. Under periodic conditions, these corrections integrate to zero over a complete (temporal) period. The new periodic condition is shown to produce accurate results for flat and wavy‐walled channels under both induced flow and countercurrent conditions. Copyright © 2003 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.     

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 second-order radiation boundary condition for the shallow water wave equations on two-dimensional unstructured finite element grids

A second-order radiation boundary condition (RBC) is derived for 2D shallow water problems posed in ‘wave equation’ form and is implemented within the Galerkin finite element framework. The RBC is derived by matching the dispersion relation for the interior wave equation with an approximate solution to the exterior problem for outgoing waves. The matching is correct to second order, accounting for curvature of the wave front and the geometry. Implementation is achieved by using the RBC as an evolution equation for the normal gradient on the boundary, coupled through the natural boundary integral of the Galerkin interior problem. The formulation is easily implemented on non-straight, unstructured meshes of simple elements. Test cases show fidelity to solutions obtained on extended meshes and improvement relative to simpler first-order RBCs.