Scour around various structures obstructing flow in an open channel is a common problem; therefore a better understanding of how turbulent flow affects sediment transport is needed. Additionally, is it the mean flow or the turbulence properties that are more important in contributing to bed shear stress? To this end, an experimental study was conducted in a fixed-bed flume containing a series of dikes. Turbulence intensities and Reynold's stresses were calculated from 3-D velocity measurements gathered with a microADV. Results showed that the maximum shear stress was nearly 12-20 times that of the approach flow, while maximum turbulence intensities were about 3-5 times those of the incoming flow. Highest magnitudes of both were seen at the tip of the second dike in the three-dike series. The mean velocity appeared to have no relation to the formation of scour near the tips of the dikes but the turbulence intensities did.

The unsteady nature of coastal hydrodynamics is associated with complex boundary layer dynamics and hence engineering predictions of shear stresses and sediment transport are difficult. This thesis explores some of the complex hydrodynamic problems and boundary layer behaviour in the coastal zone and seeks to provide new and improved modelling approaches. The latest experimental results are used to inform the model development process. New laboratory experiments carried out as part of this thesis illustrate the value of convolution integral calculations for both pressure and skin friction forces on particles on the bed. The experiments also highlight the importance of the phase differences between free stream velocity and boundary layer shear stresses. The use of a “bed” shear stress as a model input is found to be problematic whenever there is a large vertical gradient in the boundary layer shear stress. New experimental and modelling work has helped to improve our understanding of sheet flow boundary layer dynamics. This thesis builds on some of these new discoveries to propose a new simplified model framework for sheet flow sediment transport prediction using convolution integrals. This time domain technique has the advantage of simplicity while incorporating the most important physical processes from more detailed models. The new model framework could be incorporated into any depth averaged coastal hydrodynamic modelling software package. Boundary layer analysis techniques presented in the thesis provide an improved understanding of the effective roughness of mobile beds and can be used to calculate instantaneous shear stress profiles throughout the mobile bed boundary layer. New solutions for swash zone hydrodynamics are presented which illustrate the limitations of the previous benchmark analytical model for swash hydrodynamics. It is shown that real swash necessarily involves a much larger influx of mass and momentum than the analytical solution which was previously used by many in the swash sediment transport research community. Models for swash boundary layer development are also presented.

A 3-axis acoustic Doppler Velocimeter (ADV) with high sampling rate was used to measure the bed shear stress and turbulence under wave-current interaction conditions in this study. Experimental sites include laboratory tank, Love River in Kaohsiung and Howan coastal waters in Pingtung. Bed shear stress is estimated primarily by the inertial dissipation method, also by the turbulent kinetic energy method and eddy correlation method. Results of the laboratory experiments indicate that the bed shear stress increases as both the flow speed and wave height increase, and the flow speed is a more important factor. Field experiments can be divided into several types. The first type is under slow flows and calm waves. Love River is a typical condition of this type with turbid waters and a low flow speed. During the experimental period the ADV correlations reach 90% or more. Because the river flow is quite small, no significant bed shear stress is produced and u* is mostly less than 1 cm/s. As a result the deposition effect is much larger than erosion, thus a very thick layer of mud is formed on the riverbed. The observations in Howan in April 2010 also reveal the condition of slow flows and small waves, and the bed shear stress is also quite small. Due to the factors of clean coastal waters and weak turbulence in this season, the quality of ADV signals is poor. The second type is under large flows and small waves, as shown from the observations of Howan in April 2011, during which the maximum speed reached 25 cm/s and wave heights less than 20 cm. In this experiment the shear stress is large, the u* are mostly greater than 0.8 cm/s and the value of the drag coefficient is 0.0021; the ADV signals have good quality and the inertial sub-range is well defined. The third type is under weak flows and large waves. The observations of Howan in July 2011 show significant rainfall and maximum wave heights of 90 cm. In this case the u* are mostly centered around 1 cm/s. The acoustic backscatter intensity is positively correlated with the turbidity and wave height. Sizable bed shear stress induced by the orbital velocity of waves contributes a significant part to the total bed shear stress.

This research is part of the extension of the SRICOS-EFA method for predicting the maximum scour depth history around the bridge abutment. The basic objective is to establish the equation for predicting the maximum bed shear stress around the abutment at the initial condition of scouring. CHEN3D (Computerized Hydraulic ENgineering program for 3D flow) program is utilized to perform numerical simulations and predict bed shear stress before scouring. The Chimera technique incorporated in CHEN3D makes the program capable of simulating all kinds of complex geometry and moving boundary. CHEN3D program has been proven to be an accurate method to predict flow field and boundary shear stress in many fields and used in bridge scour study in cohesive soils for more than ten years. The maximum bed shear stress around abutment in open rectangular channel is studied numerically and the equation is proposed. Reynolds number is the dominant parameter, and the parametric studies have been performed based on the dimensional analysis. The influence of channel contraction ratio, abutment aspect ratio, water depth, abutment shape, and skew angle has been investigated, and the corresponding correction This research is part of the extension of the SRICOS-EFA method for predicting the maximum scour depth history around the bridge abutment. The basic objective is to establish the equation for predicting the maximum bed shear stress around the abutment at the initial condition of scouring. CHEN3D (Computerized Hydraulic ENgineering program for 3D flow) program is utilized to perform numerical simulations and predict bed shear stress before scouring. The Chimera technique incorporated in CHEN3D makes the program capable of simulating all kinds of complex geometry and moving boundary. CHEN3D program has been proven to be an accurate method to predict flow field and boundary shear stress in many fields and used in bridge scour study in cohesive soils for more than ten years. The maximum bed shear stress around abutment in open rectangular channel is studied numerically and the equation is proposed. Reynolds number is the dominant parameter, and the parametric studies have been performed based on the dimensional analysis. The influence of channel contraction ratio, abutment aspect ratio, water depth, abutment shape, and skew angle has been investigated, and the corresponding correction factors have been proposed. The study of the compound channel configuration is conducted further to extend the application of the proposed equation. Numerical simulations of overtopping flow in straight rectangular channel, straight compound channel and channel bend have been conducted. The bridge deck is found to be able to change the flow distribution and the bed shear stress will increase significantly once overtopping. The influence of the channel bend curvature, abutment location in the channel bend, and the abutment shape is also investigated. The corresponding variation of the bed shear stress has been concluded. The scour models, including the erosion rate function, roughness effect, and the turbulence kinetic energy, have been proposed and incorporated into the CHEN3D program. One flume test case in NCHRP 24-15(2) has been simulated to determine the parameters for the roughness and the turbulence kinetic energy. The prediction of the maximum scour depth history with the proposed model is in good agreement with the measurement for most cases. The influence of overtopping flow on the abutment scour development is also studied and the corresponding correction factor is proposed.

The development and testing of a shear cell for the purpose of measuring swash zone bed shear stress is presented. Direct measurements of bed shear stress were subsequently obtained using the shear plate in small, medium, and large-scale laboratory facilities. Measurements from both dam- break and bore-driven swash experiments are considered, covering a wide range of hydrodynamics and bed roughness. The dam-break problem is of interest here due to the theoretical analogy with the run-up of a solitary bore on a beach. Estimates of the flow velocities through the full swash cycle were obtained through numerical modelling and verified against measured velocity data. In conjunction, these data are used to calculate skin friction coefficients. The measurements indicate strong temporal and spatial variation in bed shear stress throughout the swash cycle, and a clear distinction between the uprush and backwash phase. For a single swash event, the maximum uprush bed shear stresses occur in the lower swash zone, within the range 0<x/Rx<0.3. The maximum backwash bed shear stresses also occur in the lower swash zone, and extend seaward of the initial bore collapse location. For a given cross-shore location the peak uprush bed shear stress is typically greater than the peak backwash bed shear stress by at least a factor two and up to a factor four. Local skin friction coefficients also indicate strong temporal and spatial variation. Furthermore, the behaviour of the local skin friction coefficient (back calculated from the measured bed shear stress using predicted, depth-averaged, flow velocities) over the swash cycle is inconsistent with the classical behaviour that is expected on the basis of the low Reynolds number flow. Smooth bed dam break and swash uprush friction coefficients appear to follow the general behaviour observed for smooth, turbulent open channel flow for an increasing Reynolds number. However, for a decreasing Reynolds number the behaviour of Cf differs from the steady flow relation. This is attributed the unsteady swash flow regime and flow history effects. It is expected that differences in flow history between the uprush and backwash have implications in terms of swash boundary layer growth and the resulting bed shear stress. A Lagrangian model for the swash boundary layer development is presented to consider these flow history effects. The model is based on the momentum integral approach for steady, turbulent, flat-plate boundary layers, with appropriate modifications to account for the unsteady flow regime. Fluid particle trajectories and velocity are computed and the boundary layer growth across the entire swash zone is estimated. Predictions of the bed shear stress agree well with the direct bed shear stress measurements and show a bias toward uprush sediment transport which has consistently been observed in measurements.

Yes / The bed slope and the tailwater depth are two important ones among the factors that affect the propagation of the dam-break flood and Favre waves. Most previous studies have only focused on the macroscopic characteristics of the dam-break flows or Favre waves under the condition of horizontal bed, rather than the internal movement characteristics in sloped channel. The present study applies two numerical models, namely, large eddy simulation (LES) and shallow water equations (SWEs) models embedded in the CFD software package FLOW-3D to analyze the internal movement characteristics of the dam-break flows and Favre waves, such as water level, the velocity distribution, the fluid particles acceleration and the bed shear stress, under the different bed slopes and water depth ratios. The results under the conditions considered in this study show that there is a flow state transition in the flow evolution for the steep bed slope even in water depth ratio α = 0.1 (α is the ratio of the tailwater depth to the reservoir water depth). The flow state transition shows that the wavefront changes from a breaking state to undular. Such flow transition is not observed for the horizontal slope and mild bed slope. The existence of the Favre waves leads to a significant increase of the vertical velocity and the vertical acceleration. In this situation, the SWEs model has poor prediction. Analysis reveals that the variation of the maximum bed shear stress is affected by both the bed slope and tailwater depth. Under the same bed slope (e.g., S0 = 0.02), the maximum bed shear stress position develops downstream of the dam when α = 0.1, while it develops towards the end of the reservoir when α = 0.7. For the same water depth ratio (e.g., α = 0.7), the maximum bed shear stress position always locates within the reservoir at S0 = 0.02, while it appears in the downstream of the dam for S0 = 0 and 0.003 after the flow evolves for a while. The comparison between the numerical simulation and experimental measurements shows that the LES model can predict the internal movement characteristics with satisfactory accuracy. This study improves the understanding of the effect of both the bed slope and the tailwater depth on the internal movement characteristics of the dam-break flows and Favre waves, which also provides a valuable reference for determining the flood embankment height and designing the channel bed anti-scouring facility. / National Natural Science Foundation of China (Grant No: 51879179, 52079081), the Open Fund from the State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University (SKHL1809) and the Sichuan Science and Technology Program (No. 2019JDTD0007)

京都大学 / 新制・課程博士 / 博士(農学) / 甲第23528号 / 農博第2475号 / 新制||農||1087(附属図書館) / 学位論文||R3||N5359(農学部図書室) / 京都大学大学院農学研究科地域環境科学専攻 / (主査)教授 藤原 正幸, 教授 中村 公人, 准教授 藤澤 和謙 / 学位規則第4条第1項該当 / Doctor of Agricultural Science / Kyoto University / DGAM

This thesis reports on investigations carried out to study of the effect of horizontal wall jets on rough, fixed and mobile beds in open channel flow. Experimental tests were carried out, using fixed and mobile sediment beds. Computer simulation models for the flow within the jet and resulting sediment transport were developed and their results analysed in this study. In the experimental phase, tests were carried out with both fixed and mobile sediment beds. The shape of the water surface, numerous point velocity measurements and measurements of the evolving scour hole shape were made. Detailed descriptions of the turbulent flow field over a fixed rough bed and for scour holes at equilibrium were obtained for a range of initial jet conditions. Fully turbulent, multiphase flow was modelled using the Fluent Computational Fluid Dynamics software. This was used to analyze the flow caused by a jet in a rectangle open-channel with a rough bed, and also the flow pattern in a channel with a local scour hole. The volume of fluid (VOF) multiphase method and K- model was used to model the fluid flow in both cases. The model predictions of velocity and shear stress were compared against experimental observations. The experimental data was used to develop new empirical relationships to describe the pattern of boundary shear stress caused by a wall jet over fixed beds and in equilibrium scour holes. These relationships were linked with existing bed-load transport rate models in order to predict the temporal evolution of scour holes. An analytical model describing the relationship between the wall jet flow and the development of a local scour hole shape was reported and its predictions compared with experimental data.

This thesis reports on investigations carried out to study of the effect of horizontal wall jets on rough, fixed and mobile beds in open channel flow. Experimental tests were carried out, using fixed and mobile sediment beds. Computer simulation models for the flow within the jet and resulting sediment transport were developed and their results analysed in this study. In the experimental phase, tests were carried out with both fixed and mobile sediment beds. The shape of the water surface, numerous point velocity measurements and measurements of the evolving scour hole shape were made. Detailed descriptions of the turbulent flow field over a fixed rough bed and for scour holes at equilibrium were obtained for a range of initial jet conditions. Fully turbulent, multiphase flow was modelled using the Fluent Computational Fluid Dynamics software. This was used to analyze the flow caused by a jet in a rectangle open-channel with a rough bed, and also the flow pattern in a channel with a local scour hole. The volume of fluid (VOF) multiphase method and K- model was used to model the fluid flow in both cases. The model predictions of velocity and shear stress were compared against experimental observations. The experimental data was used to develop new empirical relationships to describe the pattern of boundary shear stress caused by a wall jet over fixed beds and in equilibrium scour holes. These relationships were linked with existing bed-load transport rate models in order to predict the temporal evolution of scour holes. An analytical model describing the relationship between the wall jet flow and the development of a local scour hole shape was reported and its predictions compared with experimental data.

L’État haïtien, pour faire face à l’insécurité alimentaire que connaît sa population (à croissance rapide et à faible revenu), s’appuie particulièrement sur la Vallée du département de l’Artibonite qui constitue depuis toujours le véritable grenier agricole du pays, en fournissant, à elle seule, plus de 80% de la production rizicole nationale. Cette production agricole assure non seulement les besoins alimentaires de la population locale, mais également ceux particuliers des départements voisins (Ouest, Nord et Centre).L’irrigation de cette vaste plaine agricole (32400 ha de terres agricoles irrigables) est garantie, à plus de 75%, par le réseau d’irrigation établi en rive Gauche de l'ouvrage-partiteur de Canneau, alimenté par un Canal principal (CMRG) ayant un débit nominal de 40 m3/s.Néanmoins, ce Partiteur, source d’approvisionnement exclusive du CMRG, est lui-même régulé par le Barrage-réservoir de Péligre qui se trouve à 70 km plus en amont sur le fleuve de l’Artibonite. En termes hydrauliques, nos travaux ont d'abord examiné la situation du réservoir de Péligre. Alors qu'au moment de la construction il était initialement prévu d’y stocker 607 Mm3, il ne reste plus aujourd'hui qu'à peine 40% de cette capacité utile, en raison des dépôts sédimentaires qui se sont constitués année après année derrière le barrage.Cette sédimentation spectaculaire (due à l’érosion des bassins versants amont fort dégradés), combinée aux déficits pluviométriques, provoque en période d’étiage une rareté d’eau, des lâchers insuffisants et donc des déficits en eau utile envoyée en amont de ce Partiteur de Canneau. Cela impacte significativement le réseau d’irrigation aval, dont le CMRG. Nos travaux ont montré que la situation est d'autant plus critique que les débits eux-mêmes, en amont comme en aval du Partiteur, sont en fait très mal connus et devraient faire l'objet d'approches méthodologiques plus rigoureuses que celles déployées sur site actuellement.Aussi, face à ces constats, de nouvelles règles de distribution de l’eau s’imposent, pour continuer à irriguer et espérer obtenir un rendement agricole acceptable (souhaitable).Notre travail s'est ainsi donné pour objectif de contribuer à la mise en place des nouvelles règles de gestion de l’eau (répartition) au sein du réseau d’irrigation alimenté par le CMRG, pour continuer à fournir l’eau à l’irrigant en quantité acceptable (et connue). Cette amélioration de gestion a été envisagée à la fois en amont et en aval de ce réseau d’irrigation, particulièrement en ses différents nœuds-clés (Canal principal et Canal secondaire).La démarche méthodologique adoptée pour relever ce défi majeur s’appuie notamment sur un système d’information hydro-morpho-sédimentaire actualisé et de qualité. Il s’agit d’une base de données, riche en observations de hauteurs d’eau (lues aux stations limnimétriques), vitesses de surface au flotteur, champs de vitesse explorés au moyen d’un courantomètre et en données bathymétriques et granulométriques des tronçons des canaux étudiés, appréciées respectivement au moyen d’un GPS différentiel et du tamisage à sec.Les résultats fort encourageants obtenus permettent d'acquérir une meilleure compréhension du système et une amélioration particulière du réseau d’irrigation en rive gauche du Partiteur de Canneau. En s’appuyant sur les historiques de sédimentation du Lac-réservoir de Péligre (de 1960 à 2016), nous présentons un document de synthèse sur la sédimentation du Lac-réservoir de Péligre. Ce document met notamment en exergue le taux de sédimentation sévère de ce dernier (5.47 Mm3/an), qui continue d’augmenter encore aujourd’hui, ainsi que les conséquences de celui-ci sur les débits turbinés et la disponibilité de l’eau en amont du Partiteur de Canneau.Nous mettons également en évidence les formes irrégulières (Lit-non prismatique) des tronçons des canaux étudiés, via une vue axonométrique des profils en travers (issus de l’étude bathymétrique) des canaux d’irrigation en terre battue étudiés. Puis, nous présentons de manière détaillée le caractère très hétérogène des dépôts sédimentaires de ces derniers, à partir d’une analyse des représentations en Log-Probit des résultats du tamisage, construites au moyen du logiciel GrandPlots.En nous appuyant sur les mesures expérimentales des contraintes de Reynolds et des profils instantanés de vitesse (pris à intervalle de 64 ms), tirés de la base de données EPFL, nous avons montré qu’il faut absolument travailler dans les 18% inférieurs de la colonne d’eau (z/h<0.18) et en mode déficitaire, dans un écoulement turbulent comme celui-là, pour extraire de façon représentative et pertinente une pente expérimentale u*/, comme indicateur de u*.À l’issue d’un examen détaillé de la distribution verticale de la vitesse au canal secondaire FNE, nous validons un DMLWL (Dip-Modified-log-wake-law) à la fois en amont et en aval du réseau. Nous montrons que ceci permet de modéliser le Dip-phenomenon observé systématiquement au sein des profils explorés in situ. Nous proposons une relation entre le coefficient d’inégale répartition de la vitesse à la verticale αv (de Prony) et l’aspect ratio (W/h) pour tout le réseau d’irrigation étudié ;ceci afin d’obtenir une vitesse débitante (Ū), simplement à partir d’une prise de vitesse au flotteur, dans l’axe central d’écoulement.À partir des débits quantifiés à la section de référence du CMRG, via l’équation de continuité (Q=AŪ), nous fournissons un Abaque, diagramme à 3 entrées (débit (Qp), charge amont (H0) et ouverture de vanne (hv)), permettant aux vanniers de connaitre les débits au pont de fer correspondant aux différentes ouvertures de vanne et celui pour lequel le trop-plein (retour des eaux excédentaires vers le fleuve de l’Artibonite) commence à fonctionner.À l’égard des opérateurs locaux et gestionnaires du système, nous mettons enfin à disposition, des méthodes/outils simples et efficaces leur permettant de quantifier finement le débit au Canal principal en amont ainsi qu’au canal secondaire en aval, simplement à partir d’une mesure de hauteur d’eau (h) et de vitesse de surface au flotteur (Us). / Doctorat en Sciences agronomiques et ingénierie biologique / info:eu-repo/semantics/nonPublished