Despite a breadth of work investigating the spatiotemporal details of fluid dynamics over bedforms and bedload transport dynamics over flat beds, there remains a relative dearth of investigations into the spatiotemporal details of bedload transport over bedforms and on a sub-bedform scale. To address this, we conducted two sets of flume experiments focused on the two fundamental regions of flow associated with bedforms: flow separation/reattachment on the lee side of the bedform (Chapter 1; backward facing-step) and flow reacceleration up the stoss side of the next bedform (Chapter 2; two-dimensional bedform). Using Laser and Acoustic Doppler Velocimetry to record fluid turbulent events and manual particle tracking of high-speed imagery to record bedload transport dynamics, we identified the existence and importance of “permeable splat events” in the region proximal to flow reattachment.
These coupled turbulent and sediment transport events are integral to the spatiotemporal pattern of bedload transport over bedforms. Splat events are localized, high magnitude, intermittent flow features in which fluid impinges on the bed, infiltrates the top portion of bed, and then exfiltrates in all directions surrounding the point of impingement. This initiates bedload transport in a radial pattern. These turbulent structures are primarily associated with quadrant 1 and 4 turbulent structures (i.e. instantaneous fluid fluctuations in the streamwise direction that bring fluid down into the bed in the case of quadrant 1 events, or up away from the bed in the case of quadrant 4 events) and generate a distinct pattern of bedload transport compared to transport dynamics distal to flow reattachment. Distal to flow reattachment, bedload transport is characterized by relatively unidirectional transport. The dynamics of splat events, specifically their potential for inducing significant magnitudes of cross-stream transport, has important implications for the evolution of bedforms from simple, two dimensional features to complex, three-dimensional features.
New advancements in sonar technology have enabled more detailed quantification of bedload transport on the reach scale, a process paramount to the effective management of rivers with sand or gravel-dominated bed material. However, a practical and scalable field methodology for reliably estimating bedload remains elusive. A popular approach involves calculating transport from the geometry and celerity of migrating bedforms, extracted from time-series of bed elevation profiles (BEPs) acquired using echosounders. Using two sets of repeat multibeam sonar surveys from the Diamond Creek USGS gage station in Grand Canyon National Park with large spatio-temporal resolution and coverage, we compute bedload using three field techniques for acquiring BEPs: repeat multi-, single-, and multiple single-beam sonar. Significant differences in flux arise between repeat multibeam and single beam sonar. Mulitbeam and multiple single beam sonar systems can potentially yield comparable results, but the latter relies on knowledge of bedform geometries and flow that collectively inform optimal beam spacing and sampling rate. These results serve to guide design of optimal sampling, and for comparing transport estimates from different sonar configurations.
The numerical techniques and methodologies are first demonstrated on a homogeneous fluid contained within a three-dimensional lid-driven cavity. The edge state technique and linear stability analysis through Arnoldi iteration are used to resolve the complex dynamics of the canonical shear-driven benchmark problem. The techniques here lead to a dynamical description of an instability mechanism, and the work serves as a basis for the remainder of the dissertation.
A novel CFD algorithm called LEAP is currently being developed by the Kasbaoui Research Group (KRG) using the Immersed Boundary Method (IBM) to describe complex geometries. To validate the algorithm, this research project focused on testing the algorithm in three dimensions by simulating a sphere placed in a moving fluid. The simulation results were compared against the experimentally derived Schiller-Naumann Correlation. Over the course of 36 trials, various spatial and temporal resolutions were tested at specific Reynolds numbers between 10 and 300. It was observed that numerical errors decreased with increasing spatial and temporal resolution. This result was expected as increased resolution should give results closer to experimental values. Having shown the accuracy and robustness of this method, KRG will continue to develop this algorithm to explore more complex geometries such as aircraft engines or human lungs.