No abstract

Rivers typically present heterogeneous bed material, but the effects of sediment nonuniformity on river bar characteristics are still unclear. This work investigates the impact of sediment size heterogeneity on alternate bars with a morphodynamic numerical model. The model is first used to reproduce a laboratory experiment showing alternate bar formation with nonuniform bed material. Subsequently, the influence of sediment size heterogeneity on alternate bars is investigated distinguishing hybrid from free bars, definition based on the presence/absence of morphodynamic forcing, considering the results of nine scenarios. In four of them, a transverse obstacle is used to generate forcing. The computations are carried out with the Telemac‐Mascaret system solving the two‐dimensional shallow‐water equations with a finite element approach, accounting for horizontal and vertical sediment sorting processes. The results show that sediment heterogeneity affects free migrating and hybrid bars in a different way. The difference lies in the presence/absence of a migration front, so that distinct relations between bed topography, bed shear stress, and sediment sorting are obtained. Sediment sorting and associated planform redistribution of bed roughness only slightly modify free migrating bar morphodynamics, whereas hybrid bars are greatly impacted, with decreased amplitude and increased wavelength. Increased sediment size heterogeneity increases the degree of sediment sorting, while the sorting pattern remains the same for both free and hybrid bars. Moreover, it produces averagely higher, longer, and faster free bars, while in the case of hybrid bars their wavelength is increased but no general trend can be determined for their amplitude.

The study of the relationship between flow structure and morphodynamic of bars in a channel expansion/contraction is essential to better understand the processes that control the evolution of rivers. Thus, multibeam echosoundings and Acoustic Doppler Profiler (ADP) measurements were performed with a high temporal resolution in an expansion/contraction zone of the Loire River (France) occupied by bars. During the monitoring period, the macroforms presented successively an alternate, a lateral and a transverse configuration. Field data were analyzed to study how the primary and secondary velocities, the flow directions, the bed shear stresses, and the bed roughnesses (associated to dunes) evolve as a function of the water discharge and bars configuration. The bars modify the flow structure imposed by the channel width variations. In fact, the bars induce a topographic forcing which enables the separation and reducing of the mixing of two currents formed in the upstream channel expansion. This forcing is enhanced by the turbulence formed by the large dunes superimposed on the bars. Therefore, the bars promote a nonuniform flow in the channel. In turn, in the channel expansion/contraction, the migrating bars' morphodynamics are affected by the downstream channel narrowing which stops their downstream migration and forced the bars in the system. Then the nonuniformity of the flow encourages the lateral migration of the macroforms until they reach a bank and become nonmigrating. Finally, the nonmigrating bars are eroded by the flow deflected during the migration of a new bar in the channel expansion/contraction.

The flow and sediment transport processes near steep streambanks, which are commonly found in meandering, braided, and anastomosing stream systems, exhibit complex patterns that produce intricate interactions between bed and bank morphologic adjustment. Increasingly, multi-dimensional computer models of riverine morphodynamics are used to aid in the study of these processes. A number of depth-averaged two-dimensional models are available to simulate morphologic adjustment of both bed and banks. Unfortunately, these models use overly simplified conceptual models of riverbank erosion, are limited by inflexible structured mesh systems, or are unable to accurately account for the flow and sediment transport adjacent to streambanks of arbitrary geometry. A new, nonlinear model is introduced that resolves these limitations.

The Bulle effect is a phenomenon in which a disproportionately higher amount of near‐bed sediment load at a fluvial diversion moves into the diverted channel, even for cases in which the proportion of water (with respect to the main flow) entering the diversion channel is relatively small. This phenomenon has wide‐ranging implications for both engineered and natural systems: from efficient design of channels to redirect water and sediment for reclaiming sinking deltas, designing navigational channels that do not need frequent dredging, to morphological evolution of river bifurcations. The first ever, and one of the most extensive set of experiments conducted to explore this phenomenon, were conducted by Bulle in . In the current study the experiments conducted by Bulle have been simulated using an open‐source, free‐surface finite‐element‐based hydrodynamic solver. The main objectives were to explore to what extent the complex phenomenon of the Bulle effect at the scale of a laboratory experiment can be simulated accurately using Reynolds‐averaged Navier–Stokes (RANS)‐based hydrodynamic solver, and to understand the details of the hydrodynamics that Bulle could not analyze through his experiments. The hydrodynamics captured by the simulations were found to match the observations made by Bulle through his experiments, and the distributions of sediment at the diversion predicted by the numerical simulations were found to match the general trend observed in the laboratory experiments. The results from the numerical simulations were also compared with existing one‐dimensional models for sediment distribution at bifurcations, and the three‐dimensional numerical model was found to perform appreciably better. This is expected due to the complex flow features at the diversion, which can only be captured satisfactorily using a three‐dimensional hydrodynamic model. Copyright © 2017 John Wiley & Sons, Ltd.

In this paper, a second order space discontinuous Galerkin (DG) method is presented for the numerical solution of inviscid shallow water flows over varying bottom topography. Novel in the implementation is the use of HLLC and kinetic numerical fluxes 1 in combination with a dissipation operator, applied only locally around discontinuities to limit spurious numerical oscillations. Numerical solutions over (non-)uniform meshes are verified against exact solutions; the numerical error in the L 2 -norm and the convergence of the solution are computed. Bore-vortex interactions are studied analytically and numerically to validate the model; these include bores as "breaking waves" in a channel and a bore traveling over a conical and Gaussian hump. In these complex numerical test cases, we correctly predict the generation of potential vorticity by non-uniform bores. Finally, we successfully validate the numerical model against measurements of steady oblique hydraulic jumps in a channel with a contraction. In the latter case, the kinetic flux is shown to be more robust.

Abstract:The problem of quantifying the effects of flexible plants on flow resistance and eddy viscosity by vegetated floodplains is first addressed with a one-dimensional (1D) approximation based upon the so-called lateral distribution method. The estimates so obtained are then tested with two-dimensional (2D) numerical simulations based on the full shallow water equations through the use of the computational code Telemac-2D. Data obtained on a physical model of the Besòs River (Spain), whose floodplains were covered with plastic ornamental plants to mimic the effect of flexible vegetation, is used for the validation of the numerical results. Additionally, the values of flow resistance estimated numerically with the 1D and 2D simulations are compared with values obtained in a rectangular flume under flow conditions (slope, water depth and artificial lining) similar to those used on the reduced model. It is then established that as more physical mechanisms are included in the mathematical model used to study the problem, the ratio between the floodplain and the main channel flow resistance coefficient increases. The approach demonstrates that whenever enough flow data is available, the lateral distribution method delivers values of flow resistance and eddy viscosity which are highly consistent with 2D numerical modelling. This finding could mean considerable savings in the burdensome task of specifying flow resistance and turbulence dissipation values for 2D modelling of large compound channel systems.

Alpine gravel‐bed rivers are dynamic systems that have been subjected to many anthropic alterations in the past centuries. Riparian vegetation development on previously bare sediment bedforms has been a common adjustment, raising important management issues in terms of flood risks and biodiversity. Many of these rivers are also channelized, and as a result present a pattern of alternate bars. Considering recent advances in numerical biomorphodynamic modeling, this study aims at exploring numerically the morphodynamics of alternate bars in the presence of riparian vegetation. To this end, a dynamic vegetation module has been implemented on top of an existing morphodynamic model, accounting for ecological processes of seed dispersal, seedling recruitment, growth, and mortality. Numerical simulations have been performed on a simplified reach of a gravel‐bed river with free migrating alternate bars at initial state. In this work 96 scenarios have been simulated, each representing 50 years of channel evolution, with different flood regimes characterized by various peak discharges and flood durations. Yearly peak discharge variability is explicitly modeled in 48 scenarios. Model outcomes present two possible equilibrium biomorphodynamic behaviors: stationary vegetated bars, or free migrating bars in the case of frequent vegetation removal during floods. This binary behavior holds true when the stochasticity of annual peak discharges is represented, and for a wide range of parameter values included in vegetation dynamic modeling. Transient mobility of vegetated bars is observed in specific configurations where large sediment deposits deflect the flow field, eroding bar heads. Modeled bar wavelengths are in the range of values predicted for free bars by linear bar theory, and remain far from the theoretical values of hybrid, steady bars. The shift from unvegetated migrating bars to steady vegetated bars seems to show that in these simulations vegetation constitutes a hydraulic forcing, leading to a shift from free bars to forced bars, with a final configuration largely inherited from the initial state. © 2019 John Wiley & Sons, Ltd.