Sediment transport is ubiquitous in aquatic environments, and the study of sediment transport is important for both engineering and environmental reasons. However, the understanding and prediction of sediment transport are hindered by its complex dynamics and regimes. In this dissertation, the open-source solver SediFoam is developed for high-fidelity particle-resolving simulations of various sediment transport problems based on open-source solvers OpenFOAM and LAMMPS. OpenFOAM is a CFD toolbox that can perform three-dimensional flow simulations on unstructured mesh; LAMMPS is a massively parallel DEM solver for molecular dynamics. To enable the particle-resolving simulation of sediment transport on an arbitrary mesh, a diffusion-based algorithm is used in SediFoam to obtain the averaged Eulerian fields from discrete particle data. The parallel interface is also implemented for the communication of the two open-source solvers. Extensive numerical simulations are performed to validate the capability of SediFoam in the modeling of sediment transport problems. The predictions of various sediment transport regimes, including `flat bed in motion', `small dune' and `vortex dune', are in good agreement of with the experimental results and those obtained by using interface resolved simulations. The capability of the solver in the simulation of sediment transport in the oscillatory boundary layer is also demonstrated. Moreover, this well-validated high-fidelity simulation tool has been used to probe the physics of particle dynamics in self-generated bedforms in various hydraulic conditions. The results obtained by using SediFoam not only bridge the gaps in the experimental results but also help improve the engineering practice in the understanding of sediment transport. By using the particle-resolving simulation results and the insights generated therein, the closure terms in the two-fluid models or hydro-morphodynamic models can be improved, which can contribute to the numerical modeling of sediment transport in engineering scales. / Ph. D. / The study and prediction of sediment transport are important for both engineering and environmental reasons. However, the understanding of sediment transport is hindered by the complex dynamics of sediment particles in turbulent flow. In this dissertation, the open-source solver SediFoam is developed for the simulations of various sediment transport problems. Both turbulent flow and particle motions can be resolved by using SediFoam, and thus high-fidelity predictions can be provided. The SediFoam is validated extensively with respect to various sediment transport applications, including “flat bed in motion”, “dune generation and migration”, and “sediment transport in oscillatory flow”. The results obtained by using SediFoam are in good agreement of with available data in the literature. By using this well-validated high-fidelity simulation tool, the physics of particle dynamics in sediment bed and self-generated dunes are investigated. Physical insights of sediment transport that have not been captured by experimental measurements are provided by the high-fidelity simulations. Although the domain length in high-fidelity simulations is only 0.1 m, the results can also be used to improve low-fidelity numerical modeling in macro-scale engineering problems.

Thesis advisor: Noah P. Snyder / Reservoir sedimentation is a significant issue not only because it limits a reservoir’s water storage capacity and threatens its ability to meet environmental and societal needs, but also because it reduces the amount of sediment reaching downstream coastal ecosystems where sediment loading sustains critical habitat for wildlife and fisheries. Reservoir deltas in particular can decrease channel capacity and lead to an increased flood risk for populations living in the alluvial plain upstream of dams; however, studies of these landforms are underrepresented in the literature. To address that knowledge gap, this project examines reservoir sedimentation in Lake Seminole, a 123 km² surface-water impoundment created in 1954 and located at the junction of the Chattahoochee and Flint Rivers in Florida and Georgia. Where it enters Lake Seminole, the Chattahoochee River has a large subaerial delta that is actively prograding. High-resolution topographic and bathymetric datasets and historical cross section data were analyzed to measure the evolution of this delta and characterize subaerial and subaqueous sedimentation in the reservoir more broadly. In addition, a comparative land-cover change analysis was conducted for the three watersheds that drain to Lake Seminole to explore potential links between the suspended sediment generated by these surface disturbances and sedimentation patterns in each arm of the reservoir. Across the entire reservoir body, subaqueous sediment accumulated at a rate of 0.81 cm/yr between 1957–1976 and 0.12 cm/yr between 1976–2009, with the highest sedimentation rates occurring shortly after dam construction. Within the Chattahoochee arm, subaqueous sediment accumulated at a rate of 2.79 cm/yr between 1957–1976 and 0.68 cm/yr between 1976–2009, whereas, in the Flint arm, subaqueous sedimented eroded at a rate of -0.21 cm/yr between 1957–1976 and accumulated at a rate of 0.08 cm/yr between 1976–2009. Qualitative observations indicate that the erosional signal in the Flint is focused in the upstream-most portion of the arm and a depositional signal emerges farther downstream. On the Chattahoochee River delta, subaerial sediment accumulated at a rate of 1.46 ± 0.48 cm/yr between 2007–2018, with the most rapid areas of aggradation (> 4.5 cm/yr) located around the exterior edges of the delta islands. No strong link was found between the percentage of each watershed that underwent land cover change and the sedimentation rates in each reservoir arm. This analysis will provide new insights into the physical processes of reservoir sedimentation that can be used to inform river management practices and decrease the negative impacts of sediment trapping not only within reservoirs, but also upstream and downstream of dams. / Thesis (MS) — Boston College, 2024. / Submitted to: Boston College. Graduate School of Arts and Sciences. / Discipline: Earth and Environmental Sciences.

It is demonstrated that the method for locating a velocity profile origin, or plane of zero velocity, by fitting log profiles to streamwise velocity measurements is applicable to a larger range of roughness scales than previously expected. Five different sets of detailed, experimental velocity measurements were analyzed encompassing sediment-scale roughness elements, roughness caused by rigid vegetation, and large-scale roughness elements comprised of mobile bedforms. The method resulted in similar values of normalized zero-plane displacement for all roughness types considered. The ratios of zero-plane displacement, dh, to roughness height, ks, were 0.20 and 0.26 for the sediment- and vegetation-scale experiments, respectively. The results for the two experiments with bedform dominated roughness were 0.34 and 0.41. An estimate of dh/ks ranging from 0.2 to 0.4 is therefore recommended for a range of roughness types with the higher end of the range being more appropriate for the larger, bedform-scale roughness elements, and the lower end for the sediment-scale roughness elements. In addition, it is demonstrated that the location of the plane of zero velocity is temporally constant even when the bed height is not. The effects of roughness element packing density were also examined with the identification of a possible threshold at 4%, above which zero-plane displacement is independent of packing density. The findings can be applied to field velocity measurements under mobile bed conditions, facilitating the calculation of turbulence parameters such as shear velocity, by using point measurements and providing guidelines for the estimation of an appropriate value for zero-plane displacement. / Master of Science

The U.S. Army Corps of Engineers is investigating ways to improve lockage efficiency at the Marmet Locks on the Kanawha River, West Virginia. These improvements are needed to meet projected future river traffic demands. A physical effects prediction model has been developed to evaluate the changes on the river due to changes in future traffic patterns and tow configurations. The model has been verified based on the field measured data which cover a wide range of traffic conditions. The model consists of selected equations which simulate ambient velocities; wind waves; diverging waves; squat and backwater generated by the tows; and propeller jet velocities and its associated entrainment velocities. Based on the calculated velocities, concentrations of suspended sediment induced by the disturbance are estimated, at the time of disturbance and after suspended sediments have been allowed to settle. Field data on tow induced velocity changes in a back channel and a tributary are also analyzed. The effects of passing recreational boats has also been measured in the main channel. / Master of Science

Deep-marine channels are conduits for turbidity currents that transport sediment from the continental shelf to the deep ocean floor. Along the way, some of the sediment suspended in the current is deposited within the channel whereas the rest of the sediment makes its way to deeper water to construct basin-floor lobes. Slope channel deposits provide information about past climates and other perturbations to Earth systems; the transfer of nutrients, pollutants, and organic carbon; and can be important subsurface reservoirs (hydrocarbon, groundwater, etc.). However, studying the processes associated with deep-marine channels is challenging and difficult to monitor directly. The Cretaceous Tres Pasos Formation in southern Chile is exposed in an extensive and well-documented outcrop belt where preserved slope channels and their sedimentary fill can be studied. Previous work has focused on these channels at a macroscopic scale (meters to kilometers) but very few studies have been conducted on the microscopic scale. For example, little is known about how the sedimentary texture, specifically grain size and sorting, of channel-fill sandstones vary across the width and along the length of a channel. Here, we quantify grain size and sorting by hand using petrographic analysis of 57 sandstone samples from eight channel-fill transects. Our results show that there are no systematic relationships of deep-water sandstone size and sorting within channel architectural transects (e.g., from axis to margin). We also show that there are no textural trends throughout the entire study area of the Figueroa unit (~25 km length of the slope system). We also find that these sandstones are poorly to very poorly sorted, which has implications to understanding the depositional processes. Finally, we also identify that these rocks should be classified as silty sandstones rather than sandstones according to Folk (1954). Based on our results we can predict that these flows were probably well fit to the channels as well as having similar velocities, concentrations, carrying capacity for grain size in suspension, etc. across multiple flows. We can also make predictions about these samples having a detrital vs. diagenetic matrix based on their volcanic lithic rich composition and looking at hybrid event beds (HEBs) to see if similar processes could be going on at this location. We hope that this human measured dataset can be used as a baseline and for the improvement of automative measuring processes that have started to develop. / Master of Science / Sediments that come from land get transported to the ocean by way of turbidity currents which travel through submarine channels bringing those sediments with them. Along the way, sediments will get deposited in the channel before the current reaches the ocean floor giving us channel-fill sandstones. The sediments that get deposited in these channels can provide us with information on things such as past climates; transfer of nutrients, pollutants, and organic carbon; and can be important subsurface reservoirs. Studying these channels directly in the ocean can be hard to do so we study outcrops which allow us to study sediments that have been buried, lithified, and uplifted onto land. This study looks at the Cretaceous Tres Pasos Formation in southern Chile where we have lots of exposure to these antient submarine channels. Lots of previous work in this area has focused on large scale dynamics (macroscale) but there has been very little study looking at these small-scale dynamics (microscale). For example, little is known about how sedimentary texture, grain size and sorting, of channel-fill sandstones might change across the width and along the length of a channel. In this study we measure the grain size and sorting by hand using petrographic analysis of 57 sandstone samples from eight channel-fill transects. Our results show no trends in deep-water sandstone size and sorting within different transects (e.g., from the middle to outside) or across the whole study area (~25 km). We also find that these sandstones are poorly to very poorly sorted and that the classification of these rocks according to Folk (1954) is silty sandstones rather than sandstones. Based on these results we can make predictions about the flow of turbidity currents in relation to submarine channels such as flows having similar velocities, concentrations, carrying capacity for grain size in suspension, etc. across multiple flows. We can also make predictions about where the matrix of these rocks might have originated based on the composition of the rocks and by looking at another type of flow, hybrid event beds (HEBs), to see if similar processes could be going on at this location. We hope that this human measured dataset can be used as a baseline and for the improvement of automative measuring processes that have started to develop.

This dissertation aims to better understand how floc aggregate characteristics and behaviors are modified under different local conditions and how such alterations impact the floc settling velocity, which is one of the most crucial parameters influencing sediment transport modeling. A series of laboratory experiments were conducted to examine the impact of suspended sediment concentration, mixes of clay and silt, and resuspension process to equilibrium floc size and floc settling velocity. In order to observe floc size evolution, a new floc imaging acquisition was first developed. This new method allows flocs in suspended sediment concentration up to C = 400 mg/L can be imaged non intrusively. This new method was applied in all three individual studies, which are composed of this dissertation. The first chapter investigates the behaviors of flocs under constant and decay suspended sediment concentrations within a steady turbulent suspension. In the constant-concentration set of experiments, floc size time series were measured for 12 h for each of the concentration C = 15, 25, 50, 100, 200, 300, and 400 mg/L. In the decay-concentration experiments, clear water was introduced to the mixing tank, simultaneously the suspension was drained out of the mixing tank at the same rate to make the suspended sediment concentration reduce while the turbulent shear was remained unchanged. The data shows that the equilibrium floc size is a weak, positive function of concentration. For example, in order to increase 20% of floc size (approximate 22 um) the concentration needs to be increased by 700% (going from 50 to 400 mg/L). The data also illustrates that during the decrease of concentration from C = 400 to 50 mg/L, the floc size responses to the changes of concentration in the order of 10 min or less. The second chapter examines how silt particles and clay aggregates interact in a turbulent suspension. Floc sizes and settling velocity of three different suspensions, i.e., pure clay, pure silt, and a mixture of clay and silt, were monitored. The floc size data show that the presence of silt particles does not have significant impacts on clay aggregate sizes. Silt particles, however, get bound up within floc aggregates, which in turn increase the settling velocity of the floc by at least 50%. The third chapter examines whether any changes in floc properties during the deposition and resuspension processes. The floc sizes and shapes in a set of experiments with different consolidation times, concentrations, and shear patterns were measured. The conditions at which the flocs deposited or resuspended were maintained the same. The data reveal that floc size and shape of freshly deposited and after resuspended are unchanged. The erosion rate and concentration is a function of consolidation time and the applied shear stress during the deposition phase. Hence, there is a small reduction in resuspended concentration resulting in a slight decrease in resuspension floc size since floc size is also a function of concentration. / Ph. D. / Sediment transport is a narrative poem from mother nature telling us about the evolution of ancient and modern rivers, deltas, and estuaries. For thousands of years, mankind has been examining the coarser part of the poem, the gravel and sand. The finer part, the mud, has not been systematically investigated until the last 60 years. The key difference between sand and mud is the capability of mud to aggregate and form flocs which have sizes, densities, and shapes that are vastly different from the original constitutive particles. This flocculation process adds a layer of dynamics to the erosion, deposition, and transport of mud that is not present in the transport of sand. Therefore, the primary motivations for this dissertation are 1) to better understand the behavior of floc size under different conditions, e.g., in the estuaries, and 2) to provide high-quality data of floc characteristics and size evolution for model development, testing, and calibration purposes. Laboratory studies are conducted to measure the floc size and in some cases settling velocity, as a function of time under different turbulent, concentration, and sediment mixture. The findings in this dissertation help to fill the gaps of knowledge in cohesive sediment transport processes. This dissertation also suggests how floc behaviors should be accounted for under different conditions. Such information is valuable for projects such as management of sediment supplies, mitigation of land loss, restoration, and land-building diversions, e.g., on the Mississippi and Atchafalaya Rivers. Data associated with this dissertation are also available on GitHub under https://github.com/FluidSedDynamics.

Cross-shore bathymetric evolution in the nearshore zone often leads to threatening consequences such as beach erosion and shoreline retreat that concern the coastal community. A new, comprehensive cross-shore morphodynamic model was developed that can be used to describe and predict these phenomena. The study included both physical and numerical models that were designed to focus on the influence of sediment size characteristics on the cross-shore sediment transport process. For a profile equilibrium timescale, three types of beach profiles with different sediment mixtures were simulated in a small-scale, random-wave flume laboratory using erosive, storm, and accretive wave conditions. Dynamic relationships between the sediment grain sorting and beach profile changes were found to be evident as size-graded sediment fractions tended to relocate to different energetic zones along the cross-shore profiles. Existing phase-averaged wave and circulation models were utilized together with several new intra-wave modules for predicting important hydrodynamic parameters that were validated using the experimental data. A novel, multi-size sediment transport model was formulated to compute individual transport rates of size-graded sediment fractions while accounting for their interaction and non-linear size dependencies. The model was coupled with a new grain sorting model that resolves cross-shore grain sorting and vertical grain lamination. Compared to a traditional modeling approach, the new comprehensive model proved to offer superior modeling accuracy for both profile evolution and sediment grain size change. The use of the model is most advantageous for a condition with intensive grain sorting, a common scenario on a natural beach profile. Equilibrium beach profile is also better simulated by the model as size-graded fractions are predicted to relocate to different zones where they could withstand local hydrodynamics. Other new components that also help improve the modeling capability include the terms for wave-breaking and bed-slope effects, wave-crest sediment flux, and acceleration-induced bottom-shear stress. Besides superior profile modeling accuracy, sediment size characteristics and their spatial and temporal variations are also a useful set of information provided by the new model.

The Clean Water Act requires the establishment of Total Maximum Daily Loads (TMDLs) for quantifying allowable pollutant loads for stream reaches in which the biological integrity of the stream is threatened. Sediment TMDLs in urban streams are particularly difficult to establish because they require (1) reliable measurement of sediment loads and (2) the ability to locate sediment sources. This research has attempted to address these challenges through a field study of North Peachtree Creek located in DeKalb County, Georgia, which has been sampled at the Century Boulevard crossing through automatic point sampling and depth-integrated sampling. Storm events from October 2003 through October 2004 provided a field record of sediment concentration and turbidity data over a wide range of storm events. Bed and bank sediment samples were collected for comparison with the point samples and depth-integrated samples. A methodology is presented whereby point sampling is used to calculate suspended sediment discharge and turbidity analysis is used to locate and characterize sediment sources. Point samples provide the boundary condition in the Rouse solution for the vertical distribution of suspended sediment to obtain suspended sediment discharge, which is then calibrated through comparison with depth-integrated sampling. The computer model HEC-RAS (U.S. Army Corps of Engineers, 1998) was applied to the stream reach to calculate the energy grade line slope throughout each storm event for input into the sediment discharge calculations. A favorable relationship between turbidity and suspended fine sediment was found at the sampling cross-section and, through comparison with bed and bank sediment samples, was used to identify the contribution of eroded bank sediment to the total sediment discharge.

The overall objective is to improve the one-dimensional numerical prediction of the fine and non-cohesive bed material load in alluvial rivers, especially during high intensity episodes during which sediment beds are strongly remobilized. For this reason, we attempt to reduce the major inaccuracy sources coming from the alluvial resistance and bed material load relations needed to close the mathematical system. Through a shared parameter called the control factor m, the interactions occurring in alluvial rivers are incorporated more deeply into the mathematical model and more particularly into the closure laws: bed material load (SVRD, Suction-Vortex Resuspension Dynamics) and the energy slope (Verbanck et al. 2007). The control factor m is assumedly related to the Rossiter resonance modes of the separated flow downstream the bed form crest. <p><p>To further improve the representation of the flow-sediment-morphology interactions, a fully coupled model approach has been naturally chosen. In this work the terminology fully coupled means that the three equations forming the system are solved synchronously and that the terms often neglected by more traditional decoupled models are kept. <p><p>The feasibility of the new closure methodology has been drawn up by reproducing numerically the silt-flushing experiment conducted by the Yellow River Conservancy Commission (Y.R.C.C.) in the Lower Yellow River (LYR) in Northern China. The objective of the silt flushing experiment is to reverse the aggradation trend of the Lower Yellow River which, in the last decades, has become a perched river. The numerical simulation specifically reproduces the silt-flush effects in a reach of LYR located in the meandering part of the river. This reach (around 100 km) is delimited by Aishan and Luokou hydrometric stations.<p><p>Since the SVRD formulation has been developed from flume observations, the law has first been confronted to river datasets. The confrontation has revealed that the SVRD law becomes less suitable for fine sediment fluxes (ratio of water depth over median particle size > 5000). Therefore, a modified equation SVRD-2 has been built to enlarge the validity range.<p><p>The suitability of the SVRD-2 equation to predict fine sediment fluxes has been tested on data available from several hydrometric stations located in the meandering reach of the LYR: historical observations and measures collected during the flushes. The SVRD-2 has also been compared with relations specifically calibrated for this configuration. The comparison has pointed out that the performance of the two formulas is similar, which is encouraging for the SVRD-2 approach as it has not been calibrated on those data. <p><p>The closed equation system has been written on its quasi-linear form and is solved by a Finite Volume Method combined with a linearized Riemann algorithm. The numerical model has been checked up on two test cases: deposition upstream of a dam and the aggradation experiment conducted by Soni 1975. <p><p>As it is not yet possible to predict dynamically the value of the control factor m, a possible solution would be to extract its value from the measured data at the inlet cross section. Unfortunately, the necessary data are not measured locally. Moreover, a uniform value of the control factor m may not suffice to reproduce the flow along the whole reach. Therefore, it has been proposed to work temporarily in the reverse way. <p><p>From the comparison between the numerical results and the experimental data, a time evolution of the control factor m has effectively been extracted and it has been shown that it varies along the reach. At Aishan, the evolution of the control factor m corresponds to the evolution expected from the data analysis previously conducted on other data sets: the value of the control factor m decreases during the flush as it tries to reach the optimal value m=1. The time evolution at Luokou behaves differently to the one at Aishan, but remains in agreement with m evolution patterns observed historically for the river section flowing round Jinan City walls. For Luokou, the highlighted differences may come from three dimensional effects coming from the meander bend upstream the station.<p><p>Generally, the results obtained for the hydraulics, the sediment transport and bed adaptation are encouraging but still need improvements and additional feeding from the experimental data. The results for the concentration and therefore the bed elevation are very sensitive to the value of the control factor m as it influences most of the terms of the bed material load equation (SVRD-2). <p><p>The major remaining difficulties are, firstly, to deal with the rapid transients for which the model is less suitable and, secondly, to improve the prediction of the value of control factor m. Before paying more attention into the transients, enhancements concerning the flow along the reach (initial condition and discharge rates during the first days of the flush) must be conducted in priority. Indeed as the prediction of the bed or the cross section evolutions depend directly on the quality of the prediction of the sediment concentration and the hydraulics, one should first improve these aspects. To perform this study, more information about the water levels or sediment concentrations is necessary at some intermediate stations. One solution is to lengthen the studied reach, upstream to Sunkou and downstream to Lijin, totaling a river length of 456 Km.<p><p>A more entire signal of the energy slopes and the associated bed configurations at different stations would enlighten how the control factor m evolves along the reach during the silt-flush events. <p><p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Soils overlying residuum on upland divides and interfluves that formed from transported material are common in the Virginia Piedmont. They are thought to occur on the oldest landscapes in the region. A study was initiated in Appomattox County and a small portion of Buckingham County encompassing an area of 238 square miles. The origin, age and characterization of these soils is studied. Mapping units comprised of red subsoil components and mapping units with non-red subsoil components are compared. Trend-surface analysis of the elevations at which they occur and chemical and physical data from twenty-four pedons in Appomattox County are used. The mapping units contain a complex mixture of taxonomic classifications that encompass pedons with and without palic clay distributions. Wetness due to perched water tables at variable depths also affects classifications. The red subsoil mapping units tend to occupy the older landscapes. Age estimates are derived from a comparison of trend-surface elevations between the transported soils and the present-day surface. These comparisons result in age estimates of 0.8 million years to 6.25 million years BP. Therefore, the oldest geomorphic surfaces in the south central Piedmont of Virginia may be estimated as late Pliocene to Miocene age landscapes. These soil materials were deposited through a process of landscape inversion dominated by subsidence and colluviation. / Master of Science