Modelling dam-break flows over mobile beds using a 2D coupled approach

Dam-break flows usually propagate along rivers and floodplains, where the processes of fluid flow, sediment transport and bed evolution are closely linked. However, the majority of existing two-dimensional (2D) models used to simulate dam-break flows are only applicable to fixed beds. Details are given in this paper of the development of a 2D morphodynamic model for predicting dam-break flows over mobile beds. In this model, the common 2D shallow water equations are modified, so that the effects of sediment concentrations and bed evolution on the flood wave propagation can be considered. These equations are used together with the non-equilibrium transport equations for graded sediments and the equation of bed evolution. The governing equations are solved using a matrix method, thus the hydrodynamic, sediment transport and morphological processes can be jointly solved. The model employs an unstructured finite volume algorithm, with an approximate Riemann solver, based on the Roe-MUSCL scheme. A predictor–corrector scheme is used in time stepping, leading to a second-order accurate solution in both time and space. In addition, the model considers the adjustment process of bed material composition during the morphological evolution process. The model was first verified against results from existing numerical models and laboratory experiments. It was then used to simulate dam-break flows over a fixed bed and a mobile bed to examine the differences in the predicted flood wave speed and depth. The effects of bed material size distributions on the flood flow and bed evolution were also investigated. The results indicate that there is a great difference between the dam-break flow predictions made over a fixed bed and a mobile bed. At the initial stage of a dam-break flow, the rate of bed evolution could be comparable to that of water depth change. Therefore, it is often necessary to employ the turbid water governing equations using a coupled approach for simulating dam-break flows.     

A depth‐averaged two‐phase model for debris flows over erodible beds

Mass exchange between debris flow and the bed plays a vital role in debris flow dynamics. Here a depth‐averaged two‐phase model is proposed for debris flows over erodible beds. Compared to previous depth‐averaged two‐phase models, the present model features a physical step forward by explicitly incorporating the mass exchange between the flow and the bed. A widely used closure model in fluvial hydraulics is employed to estimate the mass exchange between the debris flow and the bed, and an existing relationship for bed entrainment rate is introduced for comparison. Also, two distinct closure models for the bed shear stresses are evaluated. One uses the Coulomb friction law and Manning's equation to determine the solid and fluid resistances respectively, while the other employs an analytically derived formula for the solid phase and the mixing length approach for the fluid phase. A well‐balanced numerical algorithm is applied to solve the governing equations of the model. The present model is first shown to reproduce average sediment concentrations in steady and uniform debris flows over saturated bed as compared to an existing formula underpinned by experimental datasets. Then, it is demonstrated to perform rather well as compared to the full set of USGS large‐scale experimental debris flows over erodible beds, in producing debris flow depth, front location and bed deformation. The effects of initial conditions on debris flow mass and momentum gain are resolved by the present model, which explicitly demonstrates the roles of the wetness, porosity and volume of bed sediments in affecting the flow. By virtue of extended modeling cases, the present model produces debris flow efficiency that, as revealed by existing observations and empirical relations, increases with initial volume, which is enhanced by mass gain from the bed. Copyright © 2017 John Wiley & Sons, Ltd.     

Cutting management of riparian vegetation by using hydrodynamic model simulations

This study numerically investigates effects of cutting riparian vegetation on flow characteristics by using a two-dimensional numerical model. The numerical model is based on depth-averaging the time- and volume-averaged Navier–Stokes equation with turbulent effects determined by the standard k–ε turbulence model. Drag forces exerted by the flow on vegetation are considered by adding source terms into momentum equations. In a rectangular channel and compound channel with vegetation along one side, numerical predictions show are in good agreement with those of previous studies. Five cutting scenarios, including the original, cutting along the main channel side, cutting along the bank side, alternative cutting, and reducing vegetative density, are analyzed in this study. The influences of the cutting scenarios on hydrodynamic behaviors are evaluated via numerical simulations. Simulation results suggest that cutting along the main channel side is the most effective scenario for reducing water depth and flow velocities.     

A robust well-balanced finite volume model for shallow water flows with wetting and drying over irregular terrain

An unstructured Godunov-type finite volume model is developed for the numerical simulation of geometrically challenging two-dimensional shallow water flows with wetting and drying over convoluted topography. In the framework of sloping bottom model, a modified formulation of shallow water equations is used to preserve mass conservation during flooding and recession. The key ingredient of the model is the use of this combination of the sloping bottom model and the modified shallow water equations to provide a robust technique for wet/dry fronts tracking and, together with centered discretization of the bed slope source term, to exactly preserve the static flow on irregular topographies. The variable reconstruction technique ensures nonnegative reconstructed water depth and reasonable reconstructed velocity, and the friction terms are solved by semi-implicit scheme that does not invert the direction of velocity components. The robustness and accuracy of the proposed model are assessed by comparing numerical and reference results of extensive test cases. Moreover, the results of a dam-break flooding over real topography are presented to show the capability of the model on field-scale application.     

Simulation of dam-break waves on movable beds using a multi-stage centered scheme

This paper presents the application of the multi-stage first-order centered scheme GMUSTA to solve a two-phase flow model with four equations for simulating dam-break floods without and with sediment transport.Computation of generalized Riemann invariants can be particularly complex and costly in simulating dam-break floods with sediment transport.GMUSTA numerical scheme,which does not require complete knowledge of the eigenstructure of the hyperbolic mathematical model,offers a suitable and attractive option.The quality of the dam-break flood simulations with GMUSTA scheme is verified by comparing the results against laboratory tests and some experimental data available in the literature,on fixed and mobile bed conditions,with different bed materials and flush or stepped bottoms.The numerical results reproduce quite well the experimental evidence,proving that the model is capable of predicting the temporal evolution of the free-surface and the bed.The GMUSTA scheme,which is not only simple to implement but also both accurate and computationally efficient,is proposed as an appropriate tool for integrating non-equilibrium sediment-transport models.     

Modeling depth-averaged velocity and bed shear stress in compound channels with emergent and submerged vegetation

This paper presents an approach to modeling the depth-averaged velocity and bed shear stress in compound channels with emergent and submerged vegetation. The depth-averaged equation of vegetated compound channel flow is given by considering the drag force and the blockage effect of vegetation, based on the Shiono and Knight method (1991) [40]. The analytical solution to the transverse variation of depth-averaged velocity is presented, including the effects of bed friction, lateral momentum transfer, secondary flows and drag force due to vegetation. The model is then applied to compound channels with completely vegetated floodplains and with one-line vegetation along the floodplain edge. The modeled results agree well with the available experimental data, indicating that the proposed model is capable of accurately predicting the lateral distributions of depth-averaged velocity and bed shear stress in vegetated compound channels with secondary flows. The secondary flow parameter and dimensionless eddy viscosity are also discussed and analyzed. The study shows that the sign of the secondary flow parameter is determined by the rotational direction of secondary current cells and its value is dependent on the flow depth. In the application of the model, ignoring the secondary flow leads to a large computational error, especially in the non-vegetated main channel.     

Numerical calculation for bed variation in compound-meandering channel using depth integrated model without assumption of shallow water flow

In a compound meandering channel, patterns of flow structures and bed variations change with increasing water depth owing to complex momentum exchange between high-velocity flow in a main channel and low-velocity flows in flood plains. We have developed a new quasi-three-dimensional model without the shallow water assumption, i.e., hydrostatic pressure distribution; our method is known as the general bottom velocity computation (BVC) method. In this method, a set of depth-integrated equations, including depth-integrated momentum and vorticity equations, are prepared for evaluating bottom velocity and vertical velocity distributions. The objective of this study is to develop a bed variation calculation method for both single and compound meandering channels by using the BVC method coupled with a sediment transport model. This paper shows that the BVC method can reproduce the pattern change of bed variation in a compound meandering channel flow with increasing relative depth. The variation in sediment transport rate due to overbank flow is explained by experimental and computational results.     

Two-dimensional numerical modeling of dam-break flows over natural terrain using a central explicit scheme

A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases.     

A 2D weakly-coupled and efficient numerical model for transient shallow flow and movable bed

Recent advances in the simulation of free surface flows over mobile bed have shown that accurate and stable results in realistic problems can be provided if an appropriate coupling between the shallow water equations (SWE) and the Exner equation is performed. This coupling can be done if using a suitable Jacobian matrix. As a result, faithful numerical predictions are available for a wide range of flow conditions and empirical bed load discharge formulations, allowing to investigate the best option in each case study, which is mandatory in these type of environmental problems. When coupling the equations, the SWE are considered but including an extra conservation law for the sediment dynamics. In this way the computational cost may become unrealistic in situations where the application of the SWE over rigid bed can be used involving large time and space scales without giving up to the adequate level of mesh refinement. Therefore, for restoring the numerical efficiency, the coupling technique is simplified, not decreasing the number of waves involved in the Riemann problem but simplifying their definitions. The effects of the approximations made are tested against experimental data which include transient problems over erodible bed. The simplified model is formulated under a general framework able to insert any desirable discharge solid load formula.     

Transverse structure of tidal and residual flow and sediment concentration in estuaries

An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s^− 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.     

Spatially-averaged oscillatory flow over a rough bed

A rigorous framework involving flow decomposition and averaging is presented, within which the mechanics of rough-(e.g., rippled-) bed oscillatory flows can be better interpreted and understood. Spatiallyaveraged equations for conservation of fluid mass and momentum are developed for analyses of rapidly-changing bed conditions, e.g., for growing ripples. Where repeated observations of the changing bed conditions are available, the ensemble and spatially-averaged versions of these equations can be used for more detailed analyses of the flow dynamics. The double-averaged (in space and phase or time) equations of mass and momentum conservation are shown to be appropriate for analyses of flows over fixed rough beds and equilibrium ripples. The value of the present framework is highlighted herein by its application to PIV-measured oscillatory-flow velocities, stresses and vorticities over growing and equilibrium wave-induced intermediate-depth orbital-vortex ripples. In particular, discussions are provided regarding the mechanisms by which gravity-induced and pressure-gradient-induced momentum is transferred to the bed, with the analysis framework naturally and explicitly including the combination of the full range of fluid stresses and boundary form and skin friction drag that is important in defining the flow mechanics.     

2D stream hydrodynamic,sediment transport and bed morphology model for engineering applications

A 2D depth‐averaged hydrodynamic, sediment transport and bed morphology model named STREMR HySeD is presented. The depth‐averaged sediment transport equations are derived from the 3D dilute, multiphase, flow equations and are incorporated into the hydrodynamic model STREMR. The hydrodynamic model includes a two‐equation turbulence model and a correction for the mean flow due to secondary flows. The suspended sediment load can be subdivided into different size classes using the continuum (two‐fluid) approach; however, only one bed sediment size is used herein. The validation of the model is presented by comparing the suspended sediment transport module against experimental measurements and analytical solutions for the case of equilibrium sediment‐laden in a transition from a rigid bed to a porous bed where re‐suspension of sediment is prevented. On the other hand, the bed‐load sediment transport and bed evolution numerical results are compared against bed equilibrium experimental results for the case of a meander bend. A sensitivity analysis based on the correction for secondary flow on the mean flow including the effect of secondary flow on bed shear stresses direction as well as the downward acceleration effect due to gravity on transverse bed slopes is performed and discussed. In general, acceptable agreement is found when comparing the numerical results obtained with STREMR HySeD against experimental measurements and analytical solutions. Copyright © 2007 John Wiley & Sons, Ltd.     

Does the permeability of gravel river beds affect near‐bed hydrodynamics?

The permeability of river beds is an important control on hyporheic flow and the movement of fine sediment and solutes into and out of the bed. However, relatively little is known about the effect of bed permeability on overlying near‐bed flow dynamics, and thus on fluid advection at the sediment–water interface. This study provides the first quantification of this effect for water‐worked gravel beds. Laboratory experiments in a recirculating flume revealed that flows over permeable beds exhibit fundamental differences compared with flows over impermeable beds of the same topography. The turbulence over permeable beds is less intense, more organised and more efficient at momentum transfer because eddies are more coherent. Furthermore, turbulent kinetic energy is lower, meaning that less energy is extracted from the mean flow by this turbulence. Consequently, the double‐averaged velocity is higher and the bulk flow resistance is lower over permeable beds, and there is a difference in how momentum is conveyed from the overlying flow to the bed surface. The main implications of these results are three‐fold. First, local pressure gradients, and therefore rates of material transport, across the sediment–water interface are likely to differ between impermeable and permeable beds. Second, near‐bed and hyporheic flows are unlikely to be adequately predicted by numerical models that represent the bed as an impermeable boundary. Third, more sophisticated flow resistance models are required for coarse‐grained rivers that consider not only the bed surface but also the underlying permeable structure. Overall, our results suggest that the effects of bed permeability have critical implications for hyporheic exchange, fluvial sediment dynamics and benthic habitat availability. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Turbulence in mobile-bed streams

This study is devoted to quantify the near-bed turbulence parameters in mobile-bed flows with bed-load transport. A reduction in near-bed velocity fluctuations due to the decrease of flow velocity relative to particle velocity of the transporting particles results in an excessive near-bed damping in Reynolds shear stress (RSS) distributions. The bed particles are associated with the momentum provided from the flow to maintain their motion overcoming the bed resistance. It leads to a reduction in RSS magnitude over the entire flow depth. In the logarithmic law, the von Kármán coefficient decreases in presence of bed-load transport. The turbulent kinetic energy budget reveals that for the bed-load transport, the pressure energy diffusion rate near the bed changes sharply to a negative magnitude, implying a gain in turbulence production. According to the quadrant analysis, sweep events in mobile-bed flows are the principal mechanism of bed-load transport. The universal probability density functions for turbulence parameters given by Bose and Dey have been successfully applied in mobile-bed flows.     

Energy losses in compound open channels

This paper investigates energy losses in compound channel under non-uniform flow conditions. Using the first law of thermodynamics, the concepts of energy loss and head loss are first distinguished. They are found to be different within one sub-section (main channel or floodplain). Experimental measurements of the head within the main channel and the floodplain are then analyzed for geometries with constant or variable channel width. Results show that head loss differs from one sub-section to another: the classical 1D hypothesis of unique head loss gradient appears to be erroneous. Using a model that couple 1D momentum equations, called “Independent Sub-sections Method (ISM)”, head losses are resolved. The relative weights of head losses related to bed friction, turbulent exchanges and mass transfers between sub-sections are estimated. It is shown that water level and the discharge distribution across the channel are influenced by turbulent exchanges for (a) developing flows in straight channels, but only when the flow tends to uniformity; (b) flows in skewed floodplains and symmetrical converging floodplains for small relative flow depth; (c) flows in symmetrical diverging floodplains for small and medium relative depth. Flow parameters are influenced by the momentum flux due to mass exchanges in all non-prismatic geometries for small and medium relative depth, while this flux is negligible for developing flows in straight geometry. The role of an explicit modeling of mass conservation between sub-sections is eventually investigated.     

Flow analysis in a channel with flexible vegetation using double-averaging method

The paper addresses the problem of the resistance due to vegetation in an open channel flow, characterized by partially and fully submerged vegetation formed by colonies of bushes. The flow is characterized by significant spatial variations of velocity between vertical profiles that make the traditional approach based on time averaging of turbulent fluctuations inconvenient. A more useful procedure, based on time and spatial averaging (Double-Averaging Method) is applied for the flow field analysis and characterization. The vertical distribution of mean velocity and turbulent stresses at different spatial locations has been measured with a 3D Acoustic Doppler Velocimeter (ADV) for two different vegetation densities where fully submerged real bushes (salix pentandra) have been used. Velocity measurements were completed together with the measurements of drag exerted on the flow by bushes at different flow depths. The analysis of velocity measurements allows depicting the fundamental characteristics of both the mean flow field and turbulence. The experimental data show that the contribution of form-induced stresses to the momentum balance cannot be neglected. The mean velocity profiles and the spatially averaged turbulent intensity profiles allow inferring that the vegetation density is a driving parameter for the development of a mixing layer at the canopy top in the case of submerged vegetation. Moreover, the net upward turbulent momentum flux, evaluated with the methodology proposed by Lu and Willmarth (1973), appears to be damped for increased vegetation density; this finding can rationally explain the reduction of the suspended sediment transport capacity typically observed in free surface flows over a vegetated bed.     

Two-dimensional numerical model of two-layer shallow water equations for confluence simulation

This study presents a finite-volume explicit method to solve 2D two-layer shallow water equations. This numerical model is intended to describe two-layer shallow flows in which the superposed layers differ in velocity, density and rheology in a two-dimensional domain. The rheological behavior of mudflow or debris flow is called the Bingham fluid. Thus, the shear stress on rigid bed can be derived from the constitutive equation. The computational approach adopts the HLL scheme, a novel approach for the purpose of computing a Godunov flux and solving the Riemann problem approximately proposed by Harten, Lax and van Leer, as a basic building block, treats the bottom slope by lateralizing the momentum flux, and refines the scheme using the Strang splitting to manage the frictional source term. This study successfully performed 2D two-layer shallow water computations on a rigid bed. The proposed numerical model can describe the variety of depths and velocities of substances including water and mud, when the hyperconcentrated tributary flows into the main river. The analytical results in this study will be valuable for further advanced research and for designing or planning hydraulic engineering structures.     

Flow structure over square bars at intermediate submergence: Large Eddy Simulation study of bar spacing effect

Turbulent open-channel flow over 2D roughness elements is investigated numerically by Large Eddy Simulation (LES). The flow over square bars for two roughness regimes (k-type roughness and transitional roughness between d-type and k-type) at a relative submergence of H/k = 6.5 is considered, where H is the maximum water depth and k is the roughness height. The selected roughness configurations are based on laboratory experiments, which are used for validating numerical simulations. Results from the LES, in turn, complement the experiments in order to investigate the time-averaged flow properties at much higher spatial resolution. The concept of the double-averaging (DA) of the governing equations is utilized to quantify roughness effects at a range of flow properties. Double-averaged velocity profiles are analysed and the applicability of the logarithmic law for rough-wall flows of intermediate submergence is evaluated. Momentum flux components are quantified and roughness effect on their vertical distribution is assessed using an integral form of the DA-equations. The relative contributions of pressure drag and viscous friction to the overall bed shear stress are also reported.     

A 3-D phase-averaged model for shallow-water flow with waves in vegetated water

A 3-D shallow-water flow model has been developed to simulate the flow in coastal vegetated waters with short waves. The model adopts the 3-D phase-averaged shallow-water flow equations with radiation stresses induced by short waves. It solves the governing equations using an implicit finite volume method based on quadtree rectangular mesh in the horizontal plane and stretching mesh in the vertical direction. The flow model is coupled with a spectral wave deformation model called CMS-Wave. The wave model solves the spectral wave-action balance equation and provides wave characteristics to the flow model. The model considers the effects of vegetation on currents and waves by including the drag and inertia forces of vegetation in the momentum equations and the wave energy loss due to vegetation resistance in the wave-action balance equation. The model has been tested using several sets of laboratory experiments, including steady flows in a straight channel with submerged vegetation and in a compound channel with vegetated floodplain and random waves through a vegetated channel and on a vegetated beach slope. The calculated water levels, current velocities, and wave heights are in general good agreement with the measured data.     

Three-layer model for vertical velocity distribution in open channel flow with submerged rigid vegetation

Based on the detailed laboratory experiments and theoretical analysis, a new three-layer model is proposed to predict the vertical velocity distribution in an open channel flow with submerged vegetation. The time averaged velocity and turbulence behaviour of a steady uniform flow with fully submerged artificial rigid vegetation was measured using a 3D Micro ADV, and the vertical distribution of velocity and Reynolds shear stress at different vegetation height, vegetation density and measuring positions were obtained. The results show that the velocity profile consists of three hydrodynamic regimes (i.e. the upper non-vegetated layer, the outer and bottom layer within vegetation); accordingly different methods had been adopted to describe the vertical velocity distribution. For the upper non-vegetated layer, a modified mixing length theory combined with the concept of ‘the new vegetation boundary layer’ was adopted, and an analytical model was presented to predict the vertical velocity distribution in this region. For the bottom layer within vegetation, the depth average velocity was obtained by numerically solving the momentum equations. For the upper layer within vegetation, the analytical solution was presented by expressing the shear stress as a formula fitted to the experimental data. Finally, the analytical predictions of the vertical velocity over the whole flow depth were compared with the results obtained by other researchers, and the good agreement proved that the three-layer model can be used to predict the velocity distribution of the open channel flow with submerged rigid vegetation.