The improved space–time conservation element and solution element scheme for two‐dimensional dam‐break flow simulation

A new numerical scheme, namely space–time conservation element and solution element (CE/SE) method, has been used for the solution of the two‐dimensional (2D) dam‐break problem. Distinguishing from the well‐established traditional numerical methods (such as characteristics, finite difference, finite element, and finite‐volume methods), the CE/SE scheme has many non‐traditional features in both concept and methodology: space and time are treated in a unified way, which is the most important characteristic for the CE/SE method; the CEs and SEs are introduced, both local and global flux conservations in space and time rather than space only are enforced; an explicit scheme with a stagger grid is adopted. Furthermore, this scheme is robust and easy to implement. In this paper, an improved CE/SE scheme is extended to solve the 2D shallow water equations with the source terms, which usually plays a critical role in dam‐break flows. To demonstrate the accuracy, robustness and efficiency of the improved CE/SE method, both 1D and 2D dam‐break problems are simulated numerically, and the results are consistent with either the analytical solutions or experimental results. Copyright © 2011 John Wiley & Sons, Ltd.     

HLLC scheme with novel wave‐speed estimators appropriate for two‐dimensional shallow‐water flow on erodible bed

This paper presents a first‐order HLLC (Harten‐Lax‐Van Leer with contact discontinuities) scheme to solve the Saint‐Venant shallow‐water equations, including morphological evolution of the bed by erosion and deposition of sediments. The Exner equation is used to model the morphological evolution of the bed, while a closure equation is needed to evaluate the rate of sediment transport. The system of Saint‐Venant–Exner equations is solved in a fully coupled way using a finite‐volume technique and a HLLC solver for the fluxes, with a novel wave‐speed estimator adapted to the Exner equation. Wave speeds are usually derived by computing the eigenvalues of the full system, which is highly time‐consuming when no analytical expression is available. In this paper, an eigenvalue analysis of the full system is conducted, leading to simple but still accurate wave‐speed estimators. The new numerical scheme is then tested in three different situations: (1) a circular dam‐break flow over movable bed, (2) an one‐dimensional bed aggradation problem simulated on a 2D unstructured mesh and (3) the case of a dam‐break flow in an erodible channel with a sudden enlargement, for which experimental measurements are available. Copyright © 2010 John Wiley & Sons, Ltd.     

An unstructured finite volume model for dam‐break floods with wet/dry fronts over complex topography

A robust, well‐balanced, unstructured, Godunov‐type finite volume model has been developed in order to simulate two‐dimensional dam‐break floods over complex topography with wetting and drying. The model is based on the nonlinear shallow water equations in hyperbolic conservation form. The inviscid fluxes are calculated using the HLLC approximate Riemann solver and a second‐order spatial accuracy is achieved by implementing the MUSCL reconstruction technique. To prevent numerical oscillations near shocks, slope‐limiting techniques are used for controlling the total variation of the reconstructed field. The model utilizes an explicit two‐stage Runge–Kutta method for time stepping, whereas implicit treatments for friction source terms. The novelties of the model include the flux correction terms and the water depth reconstruction method both for partially and fully submerged cells, and the wet/dry front treatments. The proposed flux correction terms combined with the water depth reconstruction method are necessary to balance the bed slope terms and flux gradient in the hydrostatical steady flow condition. Especially, this well‐balanced property is also preserved in partially submerged cells. It is found that the developed wet/dry front treatments and implicit scheme for friction source terms are stable. The model is tested against benchmark problems, laboratory experimental data, and realistic application related to dam‐break flood wave propagation over arbitrary topography. Numerical results show that the model performs satisfactorily with respect to its effectiveness and robustness and thus has bright application prospects. Copyright © 2010 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.     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of dam‐ and dyke‐break hydrodynamics on dynamically adaptive quadtree grids

Flooding due to the failure of a dam or dyke has potentially disastrous consequences. This paper presents a Godunov‐type finite volume solver of the shallow water equations based on dynamically adaptive quadtree grids. The Harten, Lax and van Leer approximate Riemann solver with the Contact wave restored (HLLC) scheme is used to evaluate interface fluxes in both wet‐ and dry‐bed applications. The numerical model is validated against results from alternative numerical models for idealized circular and rectangular dam breaks. Close agreement is achieved with experimental measurements from the CADAM dam break test and data from a laboratory dyke break undertaken at Delft University of Technology. Copyright © 2004 John Wiley Sons, Ltd.     

A Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Two‐dimensional,unstructured mesh generation for tidal models

The successful implementation of a finite element model for computing shallow‐water flow requires the identification and spatial discretization of a surface water region. Since no robust criterion or node spacing routine exists, which incorporates physical characteristics and subsequent responses into the mesh generation process, modelers are left to rely on crude gridding criteria as well as their knowledge of particular domains and their intuition. Two separate methods to generate a finite element mesh are compared for the Gulf of Mexico. A wavelength‐based criterion and an alternative approach, which employs a localized truncation error analysis (LTEA), are presented. Both meshes have roughly the same number of nodes, although the distribution of these nodes is very different. Two‐dimensional depth‐averaged simulations of flow using a linearized form of the generalized wave continuity equation and momentum equations are performed with the LTEA‐based mesh and the wavelength‐to‐gridsize ratio mesh. All simulations are forced with a single tidal constituent, M₂. Use of the LTEA‐based procedure is shown to produce a superior (i.e., less error) two‐dimensional grid because the physics of shallow‐water flow, as represented by discrete equations, are incorporated into the mesh generation process. Copyright © 2001 John Wiley & Sons, Ltd.     

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

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

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.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.     

A two‐dimensional finite volume morphodynamic model on unstructured triangular grids

We discuss the application of a finite volume method to morphodynamic models on unstructured triangular meshes. The model is based on coupling the shallow water equations for the hydrodynamics with a sediment transport equation for the morphodynamics. The finite volume method is formulated for the quasi‐steady approach and the coupled approach. In the first approach, the steady hydrodynamic state is calculated first and the corresponding water velocity is used in the sediment transport equation to be solved subsequently. The second approach solves the coupled hydrodynamics and sediment transport system within the same time step. The gradient fluxes are discretized using a modified Roe's scheme incorporating the sign of the Jacobian matrix in the morphodynamic system. A well‐balanced discretization is used for the treatment of source terms. We also describe an adaptive procedure in the finite volume method by monitoring the bed–load in the computational domain during its transport process. The method uses unstructured meshes, incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep bed gradients that may form in the approximate solution. Numerical results are shown for a test problem in the evolution of an initially hump‐shaped bed in a squared channel. For the considered morphodynamical regimes, the obtained results point out that the coupled approach performs better than the quasi‐steady approach only when the bed–load rapidly interacts with the hydrodynamics. Copyright © 2009 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.     

A robust and well‐balanced scheme for the 2D Saint‐Venant system on unstructured meshes with friction source term

In the following lines, we propose a numerical scheme for the shallow‐water system supplemented by topography and friction source terms, in a 2D unstructured context. This work proposes an improved version of the well‐balanced and robust numerical model recently introduced by Duran et al. (J. Comp. Phys., 235 , 565–586, 2013) for the pre‐balanced shallow‐water equations, accounting for varying topography. The present work aims at relaxing the robustness condition and includes a friction term. To this purpose, the scheme is modified using a recent method, entirely based on a modified Riemann solver. This approach preserves the robustness and well‐balanced properties of the original scheme and prevents unstable computations in the presence of low water depths. A series of numerical experiments are devoted to highlighting the performances of the resulting scheme. Simulations involving dry areas, complex geometry and topography are proposed to validate the stability of the numerical model in the neighbourhood of wet/dry transitions. Copyright © 2015 John Wiley & Sons, Ltd.     

A Godunov method for the computation of erosional shallow water transients

A Godunov method is proposed for the computation of open‐channel flows in conditions of rapid bed erosion and intense sediment transport. Generalized shallow water equations govern the evolution of three distinct interfaces: the water free‐surface, the boundary between pure water and a sediment transport layer, and the morphodynamic bottom profile. Based on the HLL scheme of Harten, Lax and Van Leer (1983), a finite volume numerical solver is constructed, then extended to second‐order accuracy using Strang splitting and MUSCL extrapolation. Lateralisation of the momentum flux is adopted to handle the non‐conservative product associated with bottom slope. Computational results for erosional dam‐break waves are compared with experimental measurements and semi‐analytical Riemann solutions. Copyright © 2003 John Wiley & Sons, Ltd.     

Two‐dimensional modeling for stability analysis of two‐phase stratified flow

The effect of wavelength and relative velocity on the disturbed interface of two‐phase stratified regime is modeled and discussed. To analyze the stability, a small perturbation is imposed on the interface. Growth or decline of the disturbed wave, relative velocity, and surface tension with respect to time will be discussed numerically. Newly developed scheme applied to a two‐dimensional flow field and the governing Navier–Stokes equations in laminar regime are solved. Finite volume method together with non‐staggered curvilinear grid is a very effective approach to capture interface shape with time. Because of the interface shape, for any time advancement, a new grid is performed separately on each stratified field, liquid, and gas regime. The results are compared with the analytical characteristics method and one‐dimensional modeling. This comparison shows that solving the momentum equation including viscosity term leads to physically more realistic results. In addition, the newly developed method is capable of predicting two‐phase stratified flow behavior more precisely than one‐dimensional modeling. It was perceived that the surface tension has an inevitable role in dissipation of interface instability and convergence of the two‐phase flow model. Copyright © 2009 John Wiley & Sons, Ltd.     

An improved technique for solving two‐dimensional shallow water problems

A numerical model for solving the 2D shallow water equations is proposed herewith. This model is based on a finite volume technique in a generalized co‐ordinate system, coupled with a semi‐implicit splitting algorithm in which a Helmholtz equation is used for the surface elevation. Several benchmark problems have proven the good accuracy of this method in complex geometries. Nevertheless, several numerical perturbations were noted in the surface elevation. After finding the origin, a new numerical technique is suggested, to avoid these perturbations. Several severe tests are proposed to validate this technique. Copyright © 1999 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.