Combined effects of bed friction and emergent cylinder drag in open channel flow

Environmental Fluid Mechanics - Open channel flows subjected to a longitudinal transition in roughness, from bed friction to emergent cylinder drag and vice versa, are investigated experimentally...     

One-dimensional numerical modelling of dam-break waves over movable beds: application to experimental and field cases

This paper reports a numerical study on dam-break waves over movable beds. A one-dimensional (1-D) model is built upon the Saint-Venant equations for shallow water waves, the Exner equation of sediment mass conservation and a spatial lag equation for non-equilibrium sediment transport. The set of governing equations is solved using an explicit finite difference scheme. The model is tested in various idealized experimental cases, with fairly good agreement between the numerical predictions and measurements. Discrepancies are observed at the earlier stage of the dam-break wave and around the dam location due to no vertical velocity component being taken into account. Sensitivity tests confirm that the friction coefficient is an important parameter for the evaluation of sediment transport processes operating during a dam-break wave. The influence of the non-equilibrium adaptation length (or the lag distance) is negligible on the wavefront celerity and weak on the free surface and bed profiles, which indicates that one may ignore the spatial lag effect in dam-break wave studies. Finally, the simulation of the Lake Ha!Ha! dyke-break flood event shows that the model can provide relevant results if a convenient formula for computing the sediment transport capacity and an appropriate median grain diameter of riverbed material are selected.     

Automatic sensitivity analysis of a finite volume model for two-dimensional shallow water flows

Given a numerical model for solving two-dimensional shallow water equations, we are interested in the robustness of the simulation by identifying the rate of change of the water depths and discharges with respect to a change in the bottom friction coefficients. Such a sensitivity analysis can be carried out by computing the corresponding derivatives. Automatic differentiation (AD) is an efficient numerical method, free of approximation errors, to evaluate derivatives of the objective function specified by the computer program, Rubar20 for example. In this paper AD software tool Tapenade is used to compute forward derivatives. Numerical tests were done to show the robustness of the model and to demonstrate the efficiency of these AD-derivatives.     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.