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.     

Large eddy simulation of turbulent shallow water flows using multi‐relaxation‐time lattice Boltzmann model

In this paper, the standard Smagorinsky's algorithm is embedded into the multiple relaxation time (MRT) lattice Boltzmann model (LBM) for large eddy simulation (LES) of turbulent shallow water flows (MRT‐LABSWE^TM). The model is based on the two‐dimensional nonlinear shallow water equations, giving the depth‐averaged features. It is verified by applying the model in three typical cases in engineering with turbulence: (i) the flow around a square cylinder, (ii) plane cavity flow, and (iii) flows in a junction of 90°. The results obtained by the MRT‐LABSWE^TM are compared with BGK‐LABSWE^TM results and experimental data. The objectives of this study are to validate the MRT‐LABSWE^TM in a turbulence simulation and perform a comparative analysis between the results of BGK‐LABSWE^TM and MRT‐LABSWE^TM. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of a generalized multi‐layer model for 3‐D simulation of free surface flows

A mathematical model was developed for three‐dimensional (3‐D) simulation of free surface flows. In this model, the flow depth is divided into a number of layers and shallow water equations are integrated in each layer to derive the hydrodynamic equations. To give a general form to this model, each layer is assumed to be non‐horizontal with varying thickness in the flow domain. A non‐orthogonal curvilinear coordinate system is employed in the model, to allow for flexibility in dealing with the irregular geometry of natural watercourses. Due to the similarity in governing equations, two‐dimensional (2‐D) depth averaged programs can be developed into a multi‐layer model. The development for a depth averaged program and its numerical scheme is described in this paper. Experimental data and semi‐analytical solutions are used to evaluate the performance of the model. Three different cases of open channel flow are tested: 1‐flow in a straight open channel, 2‐the flow development region in a channel, and 3‐flow in a meandering channel. It is shown that the model has the capability to predict velocity distribution and secondary flows in complex 3‐D flow conditions. Copyright © 2004 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.     

Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 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.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐dimensional dam break flow simulation

Numerical modelling of shallow water flow in two dimensions is presented in this work with the results obtained in dam break tests. Free surface flow in channels can be described mathematically by the shallow‐water system of equations. These equations have been discretized using an approach based on unstructured Delaunay triangles and applied to the simulation of two‐dimensional dam break flows. A cell centred finite volume method based on Roe's approximate Riemann solver across the edges of the cells is presented and the results are compared for first‐ and second‐order accuracy. Special treatment of the friction term has been adopted and will be described. The scheme is capable of handling complex flow domains as shown in the simulation corresponding to the test cases proposed, i.e. that of a dam break wave propagating into a 45° bend channel (UCL) and in a channel with a constriction (LNEC‐IST). Comparisons of experimental and numerical results are shown. Copyright © 2000 John Wiley & Sons, Ltd.     

A depth‐integrated non‐hydrostatic finite element model for wave propagation

This paper presents a two‐dimensional finite element model for simulating dynamic propagation of weakly dispersive waves. Shallow water equations including extra non‐hydrostatic pressure terms and a depth‐integrated vertical momentum equation are solved with linear distributions assumed in the vertical direction for the non‐hydrostatic pressure and the vertical velocity. The model is developed based on the platform of a finite element model, CCHE2D. A physically bounded upwind scheme for the advection term discretization is developed, and the quasi second‐order differential operators of this scheme result in no oscillation and little numerical diffusion. The depth‐integrated non‐hydrostatic wave model is solved semi‐implicitly: the provisional flow velocity is first implicitly solved using the shallow water equations; the non‐hydrostatic pressure, which is implicitly obtained by ensuring a divergence‐free velocity field, is used to correct the provisional velocity, and finally the depth‐integrated continuity equation is explicitly solved to satisfy global mass conservation. The developed wave model is verified by an analytical solution and validated by laboratory experiments, and the computed results show that the wave model can properly handle linear and nonlinear dispersive waves, wave shoaling, diffraction, refraction and focusing. Copyright © 2013 John Wiley & Sons, Ltd.     

A meshless method for numerical simulation of depth‐averaged turbulence flows using a k‐ϵ model

A three‐dimensional numerical model is developed to analyze free surface flows and water impact problems. The flow of an incompressible viscous fluid is solved using the unsteady Navier–Stokes equations. Pseudo‐time derivatives are introduced into the equations to improve computational efficiency. The interface between the two phases is tracked using a volume‐of‐fluid interface tracking algorithm developed in a generalized curvilinear coordinate system. The accuracy of the volume‐of‐fluid method is first evaluated by the multiple numerical benchmark tests, including two‐dimensional and three‐dimensional deformation cases on curvilinear grids. The performance and capability of the numerical model for water impact problems are demonstrated by simulations of water entries of the free‐falling hemisphere and cone, based on comparisons of water impact loadings, velocities, and penetrations of the body with experimental data. For further validation, computations of the dam‐break flows are presented, based on an analysis of the wave front propagation, water level, and the dynamic pressure impact of the waves on the downstream walls, on a specific container, and on a tall structure. Extensive comparisons between the obtained solutions, the experimental data, and the results of other numerical simulations in the literature are presented and show a good agreement. Copyright © 2015 John Wiley & Sons, Ltd.     

A high‐order PIC method for advection‐dominated flow with application to shallow water waves

A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the ‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.     

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.     

Calculation of curved open channel flow using physical curvilinear non‐orthogonal co‐ordinates

There are two main difficulties in numerical simulation calculations using FD/FV method for the flows in real rivers. Firstly, the boundaries are very complex and secondly, the generated grid is usually very non‐uniform locally. Some numerical models in this field solve the first difficulty by the use of physical curvilinear orthogonal co‐ordinates. However, it is very difficult to generate an orthogonal grid for real rivers and the orthogonal restriction often forces the grid to be over concentrated where high resolution is not required. Recently, more and more models solve the first difficulty by the use of generalized curvilinear co‐ordinates (ξ,η). The governing equations are expressed in a covariant or contra‐variant form in terms of generalized curvilinearco‐ordinates (ξ,η). However, some studies in real rivers indicate that this kind of method has some undesirable mesh sensitivities. Sharp differences in adjacent mesh size may easily lead to a calculation stability problem oreven a false simulation result. Both approaches used presently have their own disadvantages in solving the two difficulties that exist in real rivers. In this paper, the authors present a method for two‐dimensional shallow water flow calculations to solve both of the main difficulties, by formulating the governing equations in a physical form in terms of physical curvilinear non‐orthogonal co‐ordinates (s,n). Derivation of the governing equations is explained, and two numerical examples are employed to demonstrate that the presented method is applicable to non‐orthogonal and significantly non‐uniform grids. Copyright © 2004 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.     

Comparing analytical and numerical solution of nonlinear two and three‐dimensional hydrostatic flows

New test cases for frictionless, three‐dimensional hydrostatic flows have been derived from some known analytical solutions of the two‐dimensional shallow water equations. The flow domain is a paraboloid of revolution and the flow is determined by the initial conditions, the nonlinear advective terms, the Coriolis acceleration and by the hydrostatic pressure. Wetting and drying is also included. Some specific properties of the exact solutions are discussed under different hypothesis and relative importance of the forcing terms. These solutions are proposed for testing the stability, the accuracy and the efficiency of numerical models to be used for simulating environmental hydrostatic flows. The computed solutions obtained with a semi‐implicit finite difference—finite volume algorithm on unstructured grid are compared with the corresponding analytical solutions in both two and three space dimension. Excellent agreement are obtained for the velocity and for the resulting water surface elevation. Comparison of the computed inundation area also shows a good agreement with the analytical solution with degrading accuracy observed when the inundation area becomes relatively large and for long simulation time. Copyright © 2006 John Wiley & Sons, Ltd.     

A σ‐coordinate non‐hydrostatic model with embedded Boussinesq‐type‐like equations for modeling deep‐water waves

A σ‐coordinate non‐hydrostatic model, combined with the embedded Boussinesq‐type‐like equations, a reference velocity, and an adapted top‐layer control, is developed to study the evolution of deep‐water waves. The advantage of using the Boussinesq‐type‐like equations with the reference velocity is to provide an analytical‐based non‐hydrostatic pressure distribution at the top‐layer and to optimize wave dispersion property. The σ‐based non‐hydrostatic model naturally tackles the so‐called overshooting issue in the case of non‐linear steep waves. Efficiency and accuracy of this non‐hydrostatic model in terms of wave dispersion and nonlinearity are critically examined. Overall results show that the newly developed model using a few layers is capable of resolving the evolution of non‐linear deep‐water wave groups. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

2D numerical flow modeling in a macro‐rough channel

A 2D numerical flow model, developed at the Research unit of Hydrology, Applied Hydrodynamics and Hydraulic Constructions at ULg, has been applied to flows in a macro‐rough channel. The model solves the shallow water equations (SWE) with a two length scale, depth‐integrated k‐type approach for turbulence modeling. Data for the comparison have been provided by experiments conducted at the Laboratory of Hydraulic Constructions at EPFL. In the experiments with different non‐prismatic channel configurations, namely large‐scale cavities at the side walls, three different 2D flow characteristics could be observed in cavities. With the used numerical model features, especially regarding turbulence and friction modeling, a single set of bottom and side wall roughness could be found for a large range of discharges investigated in a prismatic channel. For the macro rough configurations, the numerical model gives an excellent agreement between experimental and numerical results regarding backwater curves and flow patterns if the side wall cavities have low aspect ratios. For configurations with high aspect ratios, the head loss generated by the preservation of important recirculation gyres in the cavities is slightly underestimated. The results of the computations reveal clearly that the separation of turbulence sources in the mathematical model is of great importance. Indeed, the turbulence related to 2D transverse shear effects and the 3D turbulence, generated by bed friction, can have very different amplitude. When separating these two effects in the numerical models, most of the flow features observed experimentally can be reproduced accurately. Copyright © 2009 John Wiley & Sons, Ltd.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.