Three-dimensional numerical modeling of cohesive sediment transport and wind wave impact in a shallow oxbow lake

It was observed that in some closed inland lakes sediment transport was dominated by wind-induced currents, and the sediment resuspension was primarily driven by wind-induced waves. This paper presents the development and application of a three-dimensional numerical model for simulating cohesive sediment transport in water bodies where wind-induced currents and waves are important. In the model, the bottom shear stresses induced by currents and waves were calculated, and the processes of resuspension (erosion), deposition, settling, etc. were considered. This model was first verified by a simple test case consisting of the movement of a non-conservative tracer in a prismatic channel with uniform flow, and the model output agreed well with the analytical solution. Then it was applied to Deep Hollow Lake, a small oxbow lake in Mississippi. Simulated sediment concentrations were compared with available field observations, with generally good agreement. The transport and resuspension processes of cohesive sediment due to wind-induced current and wave in Deep Hollow Lake were also discussed.     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

LES-DEM simulations of sediment transport

In this work, a fully-coupled Computational Fluid Dynamics (CFD) model and Discrete Element Method (DEM) are used to simulate a unidirectional turbulent open-channel flow over the full range of sediment transport regimes. The fluid and particles are computed on separate grids using a dual-grid formulation to maintain consistency and avoid instability issues. The results of coupling the dispersed phase to a multiphase flow solver that uses volume-averaged Navier-Stokes equations are compared to those obtained from coupling through drag to a single flow solver. The current work also examines the applicability and limitations of lumping particles as a representative particle to reduce the cost of simulations. Insight to the impact of different turbulent events to the entrainment of particles is also given. The simulation results of sediment transport from both coupling techniques show good agreement with empirical formulas in the bedload regime, but under-predict sediment transport in the suspended load regime. In the suspended load regime, using partial coupling, the rate of sediment transport was found to be under-predicted as compared to full-coupling. The deviation in results in the suspended load regime was found to increase with increases in the applied shear stress. Both coupling methods revealed the same effect on the friction factor where friction increases in the bedload regime and decreases in the suspended load regime reaching a maximum at the transition between regimes. This result is contrary to past studies which have shown a discrete jump in the friction factor at the transition. Lumping particles as representative particles is shown to reduce the simulation cost by more than a factor of 5 when using a scaling factor of 2. By doing a quadrant analysis on information obtained from particle and flow field results, it was found that most of the particles are entrained by more frequent sweep events.     

Soil erosion due to rainfall impact with no inflow: a numerical solution with spatial and temporal effects of sediment settling velocity characteristics

Dynamic changes take place in the nature of sediment eroded from bare soil at low slopes by rainfall impact when there is no inflow of water at the top of the eroding slope. This relates initially to fine soil sediment not settling back onto the soil after the rainfall impact. Coupled partial differential equations describing such dynamic changes have been solved numerically for a bed of soil, bounded at its upper end, and subject to a constant rainfall rate. This solution allows prediction of the change with time and downslope distance in the concentration and settling velocity (or size) characteristics of eroding sediment, allowing critical evaluation of the assumption of space-independent sediment characteristics made in prior approximate analytical solutions of the equations involved. Following the determination of as yet unpredictable soil-related parameters in the equations, the solution was tested by comparison with experimented data on two soils of contrasting structural stability, namely a vertosol [The Australian Soil Classification (1996)] and a aridisol. Investigations included the determination of a minimum number of sediment size classes required to adequately describe the settling velocity characteristics, based on the shape of the underlying basic settling velocity characteristic, which is used to predict the dynamics of sediment deposition. The effect on the solution of observed structural breakdown in soil aggregation due to rainfall impact was investigated, leading to more accurate predictions of the settling velocity characteristics of eroded sediment. Other sources of discrepancy between theory and observation remain to be determined.     

Suspended sediment transport in the Gulf of Lions (NW Mediterranean): Impact of extreme storms and floods

In situ observations were combined with 3D modeling to gain understanding of and to quantify the suspended sediment transport in the Gulf of Lions (NW Mediterranean Sea). The outputs of a hydrodynamic–sediment transport coupled model were compared to near-bottom current and suspended sediment concentration measurements collected at the head of seven submarine canyons and at a shallow shelf site, over a 6-month period (November 2003–May 2004). The comparisons provide a reasonable validation of the model that reproduces the observed spatial and time variations. The study period was marked by an unusual occurrence of marine storms and high river inputs. The major water and sediment discharges were supplied by the Rhone, the largest Mediterranean river, during an exceptional flood accompanying a severe marine storm in early December 2003. A second major storm, with moderate flooding, occurred in February 2004. The estimate of river input during the studied period was 5.9 Mt. Our study reveals (i) that most of the particulate matter delivered by the Rhone was entrapped on the prodelta, and (ii) that marine storms played a crucial role on the sediment dispersal on the shelf and the off-shelf export. The marine storms occurring in early December 2003 and late February 2004 resuspended a very large amount of shelf sediment (>8 Mt). Erosion was controlled by waves on the inner shelf and by energetic currents on the outer shelf. Sediment deposition took place in the middle part of the shelf, between 50 and 100 m depth. Resuspended sediments and river-borne particles were transported to the southwestern end of the shelf by a cyclonic circulation induced by these onshore winds and exported towards the Catalan shelf and into the Cap de Creus Canyon which incises the slope close to the shore. Export taking place mostly during marine storms was estimated to reach 9.1 Mt during the study period.     

Multigene genetic programming for sediment transport modeling in sewers for conditions of non-deposition with a bed deposit

It is known that construction of large sewers based on consideration of flow with non-deposition without a bed deposit is not economical. Sewer design based on consideration of flow with non-deposition with a bed deposit reduces channel bed slope and construction cost in which the presence of a small depth of sediment deposition on the bed increases the sediment transport capacity of the flow. This paper suggests a new Pareto-optimal model developed by the multigene genetic programming (MGGP) technique to estimate particle Froude number (Fr_p) in large sewers with conditions of sediment deposition on the bed. To this end, four data sets including wide ranges of sediment size and concentration, deposit thickness, and pipe size are used. On the basis of different statistical performance indices, the efficiency of the proposed Pareto-optimal MGGP model is compared to those of the best MGGP model developed in the current study as well as the conventional regression models available in the literature. The results indicate the higher efficiency of the MGGP-based models for Fr_p estimation in the case of no additional deposition onto a bed with a sediment deposit. Inasmuch as the Pareto-optimal MGGP model utilizes a lower number of input parameters to yield comparatively higher performance than the conventional regression models, it can be used as a parsimonious model for self-cleansing design of large sewers in practice.     

Physical and coupled fully three-dimensional numerical modeling of pressurized bottom outlet flushing processes in reservoirs

Sediment deposition in reservoirs is an important research topic in engineering practice. Reservoir sedimentation has the potential to affect ood levels, drainage for agricultural land, pump station and hydropower operation as well as navigation. This paper describes the development of a coupled fully three-dimensional (3D) numerical model for the prediction of the local sediment ushing scour upstream of the bottom outlet. The presented numerical model solves the Navier-Stokes equations in conjunction with the k- turbulence model which includes both sediment transport and hydrodynamic parameters. The proposed coupled fully 3D numerical model is used to simulate experimental tests based on non-cohesive sediment. The geometric features of the scour hole (temporal and spatial hole devel- opment) upstream of the bottom outlet were reasonably well predicted compared to the experimental data. Furthermore, the velocity eld upstream of the bottom outlet was in good agreement with mea- surements. The proposed numerical model for bottom outlet ushing was, therefore, validated because of its ability to accurately predict the scour hole development during the ushing process. The proposed numerical model can be considered reliable provided that the model is correctly calibrated and set up to re ect the conditions of a particular case study.     

Hydrologic-hydraulic modeling of sediment transport along the main stem of a watershed: role of tributaries and channel geometry

ABSTRACT

Current estimations of sediment transport at the watershed scale are limited by the difficulty of accurately simulating the sediment transfer along the main stem. The typical approach to simulating watershed sediment transport involves the adoption of hydrologic sediment routing schemes that do not fully capture the contribution and timing of side tributaries, and the inclusion of a simplified channel geometry that does not include its hydraulic feedback. In this paper, we present the results of a coupled hydrologic-hydraulic model of sediment transport applied to a small watershed of Iowa. The model was developed to simulate both the hydrologic network and non-equilibrium sediment transport that occur during a flood. The model results highlight the importance of including side tributaries in order to capture a realistic duration of shear stress that ultimately affects sediment transport. Comparisons with bank erosion measurements indicate that the presented approach is also promising to estimate sediment sources along the main stem.     

Laboratory experiments of sediment transport from bare soil with a rill

Abstract

Mathematical models developed for quantification of sediment transport in hydrological watersheds require data collected through field or laboratory experiments, but these are still very rare in the literature. This study aims to collect such data at the laboratory scale. To this end, a rainfall simulator equipped with nozzles to spray rainfall was constructed, together with an erosion flume that can be given longitudinal and lateral slopes. Eighty experiments were performed, considering microtopographical features by pre-forming a rill on the soil surface before the start of each experiment. Medium and fine sands were used as soil, and four rainfall intensities (45, 65, 85 and 105 mm h^-1) were applied in the experiments. Rainfall characteristics such as uniformity, granulometry, drop velocity and kinetic energy were evaluated; flow and sediment discharge data were collected and analysed. The analysis shows that the sediment transport rate is directly proportional to rainfall intensity and slope. In contrast, the volumetric sediment concentration stays constant and does not change with rainfall intensity unless the slope changes. These conclusions are restricted to the conditions of experiments performed under rainfall intensities between and 105 mm h^-1 for medium and fine sands in a 136-cm-wide, 650-cm-long and 17-cm-deep erosion flume with longitudinal and lateral slopes varying between 5 and 20%.

Editor Z.W. Kundzewicz; Associate editor G. Mahé

Citation Aksoy, H., Unal, N.E., Cokgor, S., Gedikli, A., Yoon, J., Koca, K., Inci, S.B., Eris, E., and Pak, G., 2013. Laboratory experiments of sediment transport from bare soil with a rill. Hydrological Sciences Journal, 58 (7), 1505–1518.     

Three‐dimensional numerical modeling of the Bulle effect: the nonlinear distribution of near‐bed sediment at fluvial diversions

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 1926 . 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.