An evaluation of force terms in the lattice Boltzmann models in simulating shallow water flows over complex topography

Proper approximation of the force terms, especially the bed slope term, is of crucial importance to simulating shallow water flows in lattice Boltzmann (LB) models. However, there is little discussion on the schemes of adding force terms to LB models for shallow water equations (SWEs). In this study, we evaluate the performance of forcing schemes coupled with different LB models (LABSWE and MLBSWE) in simulating shallow water flows over complex topography and try to find out their intrinsic characteristics and applicability. Three cases are adopted for evaluation, including a stationary case, a one-dimensional tidal wave flow over an irregular bed, and a steady flow over a two-dimensional seamount. The simulating results are compared with analytical solutions or the results produced by the finite difference method. For LABSWE, all the forcing schemes, except for the weighting factor method, fail to produce accurate solutions for the test cases; this is probably due to the mismatch between the bed slope term in source terms and the quadratic depth term of the equilibrium distribution functions in these forcing schemes. For MLBSWE, all the forcing schemes are capable of simulating flows over the complex topography accurately; furthermore, those schemes taking into account the collision effect τ to eliminate the momentum induced by forces provide more accurate solutions with quicker convergence as the lattice size decreases. In this view, MLBSWE can bring more flexibility in treating the force terms and thus can be a better tool to simulate shallow water flows over complex topography in practical application.     

A coupled lattice Boltzmann model for the shallow water‐contamination system

A coupled numerical method for the direct simulation of shallow water dynamics and pollutant transport is formulated and implemented. The conservation equations of shallow water dynamics equations and the convection–diffusion equations are solved using the lattice Boltzmann (LB) method. The local equilibrium distribution of the pollutant has no terms of second order in flow velocity. And the relaxation time of the pollutant deviates from a constant for the flows with variable free surface water depth. The numerical tests show that this scheme strictly obeys the conservation law of mass and momentum. Excellent agreement is obtained between numerical predictions and analytical solutions in the pure diffusion problem and convection–diffusion problem. Furthermore, the influences on the accuracy of the lattice size and the diffusivity are also studied. The results indicate that the variation in the free surface water depth cannot affect the conservation of the model, and the model has the ability to simulate the complex topography problem. The comparison shows that the LB scheme has the capacity to solve the complex convection–diffusion problem in shallow water. Copyright © 2008 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of shallow water flows over discontinuous beds

Numerical modeling of shallow water flows over discontinuous beds is presented. The flows are described with the shallow water equations and the equations are solved using the lattice Boltzmann method (LBM) with single relaxation time (Bhatnagar–Gross–Krook‐LBM (BGK‐LBM)) and the multiple relaxation time (MRT‐LBM). The weighted centered scheme for force term together with the bed height for a bed slope is described to improve simulation of flows over discontinuous bed. Furthermore, the resistance stress is added to include the local head loss caused by flow over a step. Four test cases, one‐dimensional tidal over regular bed and steps, dam‐break flows, and two‐dimensional shallow water flow over a square block, are considered to verify the present method. Agreements between predictions and analytical solutions are satisfactory. Furthermore, the performance and CPU cost time of BGK‐LBM and MRT‐LBM are compared and studied. The results have shown that the lattice Boltzmann method is simple and accurate for simulating shallow water flows over discontinuous beds. This demonstrates the capability and applicability of the lattice Boltzmann method in modeling shallow water flows on bed topography with a discontinuity in practical hydraulic engineering. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Bathymetry reconstruction based on the zero-inertia shallow water approximation

The knowledge of the channel bed topography is paramount in modeling the hydrodynamics of open channel flows. Indeed, flow models based on the Shallow Water Approximation require prior information on the channel bed topography to accurately capture the flow features in natural rivers, estuaries, and flood plains. We present here a numerical technique for reconstructing the channel bed topography from given free surface elevation data for steep open channel flows for which the zero-inertia shallow water approximation holds. In this context, the shallow water equations are modified by neglecting inertia terms while retaining the effects of the bed slope and friction terms. We show in this work that by algebraic manipulation, we can recast the governing equations into a single first-order partial differential equation which describes the inverse problem which consists in finding the bed topography from known free surface elevation data. Interestingly, the analysis shows that the inverse problem does not require the knowledge of the bed roughness. The forward problem is solved using MacCormack’s explicit numerical scheme by considering unsteady modified shallow water equations. However, the inverse problem is solved using the method of characteristics. The results of the inverse and the forward problem are successfully tested against each other on two different test cases.     

A correction for balancing discontinuous bed slopes in two‐dimensional smoothed particle hydrodynamics shallow water modeling

In this paper, a smoothed particle hydrodynamics (SPH) numerical model for the shallow water equations (SWEs) with bed slope source term balancing is presented. The solution of the SWEs using SPH is attractive being a conservative, mesh‐free, automatically adaptive method without special treatment for wet‐dry interfaces. Recently, the capability of the SPH–SWEs numerical scheme with shock capturing and general boundary conditions has been used for predicting practical flooding problems. The balance between the bed slope source term and fluxes in shallow water models is desirable for reliable simulations of flooding over bathymetries where discontinuities are present and has received some attention in the framework of Finite Volume Eulerian models. The imbalance because of the source term resulting from the calculation of the the water depth is eradicated by means of a corrected mass, which is able to remove the error introduced by a bottom discontinuity. Two different discretizations of the momentum equation are presented herein: the first one is based on the variational formulation of the SWEs in order to obtain a fully conservative formulation, whereas the second one is obtained using a non‐conservative form of the free‐surface elevation gradient. In both formulations, a variable smoothing length is considered. Results are presented demonstrating the corrections preserve still water in the vicinity of either 1D or 2D bed discontinuities and provide close agreement with 1D analytical solutions for rapidly varying flows over step changes in the bed. The method is finally applied to 2D dam break flow over a square obstacle where the balanced formulation improves the agreement with experimental measurements of the free surface. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical model for the flooding and drying of irregular domains

A numerical technique for the modelling of shallow water flow in one and two dimensions is presented in this work along with the results obtained in different applications involving unsteady flows in complex geometries. A cell‐centred finite volume method based on Roe's approximate Riemann solver across the edges of both structured and unstructured cells is presented. The discretization of the bed slope source terms is done following an upwind approach. In some applications a problem arises when the flow propagates over adverse dry bed slopes, so a special procedure has been introduced to model the advancing front. It is shown that this modification reproduces exactly steady state of still water in configurations with strong variations in bed slope and contour. The applications presented are mainly related with unsteady flow problems. The scheme is capable of handling complex flow domains as will be shown in the simulations corresponding to the test cases that are going to be presented. Comparisons of experimental and numerical results are shown for some of the tests. Copyright © 2002 John Wiley & Sons, Ltd.     

An unstructured finite‐volume method for coupled models of suspended sediment and bed load transport in shallow‐water flows

The aim of this work is to develop a well‐balanced finite‐volume method for the accurate numerical solution of the equations governing suspended sediment and bed load transport in two‐dimensional shallow‐water flows. The modelling system consists of three coupled model components: (i) the shallow‐water equations for the hydrodynamical model; (ii) a transport equation for the dispersion of suspended sediments; and (iii) an Exner equation for the morphodynamics. These coupled models form a hyperbolic system of conservation laws with source terms. The proposed finite‐volume method consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms. The gradient fluxes are discretized using a modified Roe's scheme using the sign of the Jacobian matrix in the coupled system. A well‐balanced discretization is used for the treatment of source terms. In this paper, we also employ an adaptive procedure in the finite‐volume method by monitoring the concentration of suspended sediments in the computational domain during its transport process. The method uses unstructured meshes and incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep sediment concentrations and bed load gradients that may form in the approximate solutions. Details are given on the implementation of the method, and numerical results are presented for two idealized test cases, which demonstrate the accuracy and robustness of the method and its applicability in predicting dam‐break flows over erodible sediment beds. The method is also applied to a sediment transport problem in the Nador lagoon.Copyright © 2013 John Wiley & Sons, Ltd.     

Application of a second‐order Runge–Kutta discontinuous Galerkin scheme for the shallow water equations with source terms

The present work addresses the numerical prediction of discontinuous shallow water flows by the application of a second‐order Runge–Kutta discontinuous Galerkin scheme (RKDG2). The unsteady flow of water in a one‐dimensional approach is described by the Saint Venant's model which incorporates source terms in practical applications. Therefore, the RKDG2 scheme is reformulated with a simple way to integrate source terms. Further, an adequate boundary conditions handling, by the theory of characteristics, was overviewed to be adapted to the external points of the mesh, as well as to some points of local invalidity of the Saint Venant's model. To validate the proposed technique, steady and transient test problems (all having a reference solution) were considered and computed by means of the overall method. The results were illustrated jointly with the reference solution and the results carried out by a traditional second‐order finite volume (FV2) scheme implemented with the same techniques as the RKDG2. The proposed method has proven its practical consideration when solving discontinuous shallow water flow involving: non‐prismatic channels, various cross‐sections, smoothly varying bed topography and internal boundary conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry

The lattice‐Boltzmann (LB) method, derived from lattice gas automata, is a relatively new technique for studying transport problems. The LB method is investigated for its accuracy to study fluid dynamics and dispersion problems. Two problems of relevance to flow and dispersion in porous media are addressed: (i) Poiseuille flow between parallel plates (which is analogous to flow in pore throats in two‐dimensional porous networks), and (ii) flow through an expansion–contraction geometry (which is analogous to flow in pore bodies in two‐dimensional porous networks). The results obtained from the LB simulations are compared with analytical solutions when available, and with solutions obtained from a finite element code (FIDAP) when analytical results are not available. Excellent agreement is found between the LB results and the analytical/FIDAP solutions in most cases, indicating the utility of the lattice‐Boltzmann method for solving fluid dynamics and dispersion problems. Copyright © 1999 John Wiley & Sons, Ltd.     

A modified lattice Boltzmann model for shallow water flows over complex topography

A modified lattice Boltzmann model is proposed to describe shallow water flows over complex topography. In the proposed model, the quadratic depth term is excluded from the equilibrium distribution functions (EDFs), and the hydrostatic pressure term is combined with the bed slope term to be treated as a part of the sourcing term in the lattice Boltzmann equation (LBE). Therefore, it is unnecessary to match the coefficients of the quadratic depth term in the EDFs with those of the bed slope term in the sourcing terms in the LBE. This would bring more flexibility to the treatment of the sourcing terms in the LBE. In order to recover the shallow water equations (SWEs), the basic constraints are redefined, and under these constraints, the coefficients of the EDFs are derived afterwards. Several benchmark problems are used to validate the proposed model, including stationary case, steady flows over a two‐dimensional bump and tidal wave flows over irregular bed elevation. The computed results are in excellent agreement with the results of the other numerical methods and the analytical solutions, indicating that the proposed model is capable of simulating shallow water flows over complex bathymetry. It also proves that the proposed model has potential to produce competitive solutions to shallow water flows over complex bed topography. Copyright © 2014 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 robust high‐resolution finite volume scheme for the simulation of long waves over complex domains

The propagation, runup and rundown of long surface waves are numerically investigated, initially in one dimension, using a well‐balanced high‐resolution finite volume scheme. A conservative form of the nonlinear shallow water equations with source terms is solved numerically using a high‐resolution Godunov‐type explicit scheme coupled with Roe's approximate Riemann solver. The scheme is also extended to handle two‐dimensional complex domains. The numerical difficulties related to the presence of the topography source terms in the model equations along with the appearance of the wet/dry fronts are properly treated and extended. The resulting numerical model accurately describes breaking waves as bores or hydraulic jumps and conserves volume across flow discontinuities. Numerical results show very good agreement with previously presented analytical or asymptotic solutions as well as with experimental benchmark data. Copyright © 2007 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.     

Simulation of shallow flows over variable topographies using unstructured grids

Simulation of shallow flows over variable topographies is a challenging case for most available shock‐capturing schemes. This problem arises because the source terms and flux gradients are not balanced in the numerical computations. Treatments for this problem generally work well on structured grids, but they are usually too expensive, and most of them are not directly applicable to unstructured grids. In this paper we propose two efficient methods to treat the source terms without upwinding and to satisfy the compatibility condition on unstructured grids. In the first method, the calculation of the bed slope source term is performed by employing a compatible approximation of water depth at the cell interfaces. In the second one, different components of the bed slope term are considered separately and a compatible discretization of the components is proposed. The present treatments are applicable for most schemes including the Roe's method without changing the performance of the original scheme for smooth topographies. Copyright © 2006 John Wiley & Sons, Ltd.     

Flux and source term discretization in two‐dimensional shallow water models with porosity on unstructured grids

Two‐dimensional shallow water models with porosity appear as an interesting path for the large‐scale modelling of floodplains with urbanized areas. The porosity accounts for the reduction in storage and in the exchange sections due to the presence of buildings and other structures in the floodplain. The introduction of a porosity into the two‐dimensional shallow water equations leads to modified expressions for the fluxes and source terms. An extra source term appears in the momentum equation. This paper presents a discretization of the modified fluxes using a modified HLL Riemann solver on unstructured grids. The source term arising from the gradients in the topography and in the porosity is treated in an upwind fashion so as to enhance the stability of the solution. The Riemann solver is tested against new analytical solutions with variable porosity. A new formulation is proposed for the macroscopic head loss in urban areas. An application example is presented, where the large scale model with porosity is compared to a refined flow model containing obstacles that represent a schematic urban area. The quality of the results illustrates the potential usefulness of porosity‐based shallow water models for large scale floodplain simulations. Copyright © 2005 John Wiley & Sons, Ltd.     

An efficient and flexible solver for the simulation of the morphodynamics of fast evolving flows on coarse sediment beaches

A well‐balanced total variation diminishing–McCormack scheme is used to simulate the fast evolving flow on a mobile coarse sediments bed. The scheme is chosen because of its shock capturing capabilities and its relative simplicity, which allow different sediment transport formulae to be slotted in easily. A typical example of the kind of flows treated here is bore‐driven wave run‐up. The analogy with a dam‐break on a mobile bed is used here to analyze the performance of the model. The model solves the nonlinear shallow water equations coupled with the Exner sediment balance equation for the mobile bed. Quasi‐analytical solutions to this problem for different expressions for instantaneous sediment discharge formulae are used to test the performance of the scheme. Together with the existing solution for the Grass formula, a further solution is obtained for a different formula. Numerical tests were also carried out for a further formula that is an industry standard. The agreement of the results with the solutions is very good and consistent results were obtained for all the formulae tested. Copyright © 2011 John Wiley & Sons, Ltd.     

用离散速度方法计算浅水长波方程

用离散速度法计算浅水波方程，将空气动力学方程和浅水波方程作了比较，用Nadiga提出的近平衡流动方法模拟浅水波方程的连续和间断解。计算了一维的溃坝波问题和Thacker提出的连续解问题，结果与精确解作了比较，并且计算了水流跃过障碍物的问题。     

Propagation of discontinuouswaves along a dry bed

The possibility of simulating the process of propagation of discontinuous waves along a dry bed on the basis of the equations of the first approximation of shallow water theory is studied. It is shown that the consistent losses of the free-stream total momentum and energy can be found from the mass conservation law within the framework of the shallow water equations. As an example, solutions of the problem of dam break with a dry bed in the lower pool and the problem of impingement of a discontinuous wave on a coastal shelf are constructed. These exact solutions are compared with the results of laboratory experiments.