Water‐worked gravel beds in laboratory flumes – a natural analogue?

An investigation has been conducted to identify the key parameters that are likely to scale laboratory sediment deposits to the field scale. Two types of bed formation were examined: one where sediment is manually placed and screeded and the second where sediment is fed into a running flume. This later technique created deposits through sequential cycles of sediment transport and deposition. Detailed bed surface topography measurements have been made over a screeded bed and three fed beds. In addition, bulk subsurface porosity and hydraulic conductivity have been measured. By comparing the four beds, results revealed that certain physical properties of the screeded bed were clearly different from those of the fed beds. The screeded bed had a random organization of grains on both the surface and within the subsurface. The fed beds exhibited greater surface and subsurface organization and complexity, and had a number of properties that closely resembled those found for water‐worked gravel beds. The surfaces were water‐worked and armoured and there was preferential particle orientation and direction of imbrication in the subsurface. This suggested that fed beds are able to simulate, in a simplified manner, both the surface and subsurface properties of established gravel‐bed river deposits. The near‐bed flow properties were also compared. It revealed that the use of a screeded bed will typically cause an underestimation in the degree of temporal variability in the flow. Furthermore, time‐averaged streamwise velocities were found to be randomly organized over the screeded bed but were organized into long streamwise flow structures over the fed beds. It clearly showed that caution should be taken when comparing velocity measurements over screeded beds with water‐worked beds, and that the formation of fed beds offers an improved way of investigating intragravel flow and sediment–water interface exchange processes in gravel‐bed rivers at a laboratory scale. Copyright © 2008 John Wiley & Sons, Ltd.     

Examining the physical components of boundary shear stress for water‐worked gravel deposits

It is argued in this commentary that, in order to understand better the physical mechanisms that generate boundary shear stress over water‐worked gravel beds, flow velocity data should be re‐evaluated by spatial averaging the Reynolds equations to produce time‐ and space‐averaged (double‐averaged) momentum equations. A series of laboratory experiments were conducted in which the flow velocities were measured using a PIV system over two water‐worked gravel deposits. Combined with detailed data on the bed surface topography and vertical porosity, the physical components of shear stress were obtained. This enabled the various momentum transfer mechanisms present above, within and at the interface of a porous, fluvial deposit, to be quantified. This included the examination of the relevant contributions of temporal and spatial fluctuations in velocity and surface drag to the overall momentum transfer. It is demonstrated that double‐averaging represents a logical framework for assessing the fluid forces responsible for sediment entrainment and for investigating intragravel flow and sediment–water interface exchange mechanisms within the roughness layer in water‐worked gravel deposits. By considering the physical components of shear stress and their relative sizes it was possible to provide a physically based explanation for existing observations of enhanced mobility of gravel–sand mixtures and the transfer of solutes into porous, gravel deposits. This analysis reveals the importance of obtaining co‐located, high quality spatial data on the flow field and bed surface topography in order to gain a physical understanding of the mechanisms which generate boundary shear stress. Copyright © 2010 John Wiley & Sons, Ltd.     

Time-averaged flow structure in the central region of a stream confluence

Previous process-oriented field studies of stream confluences have focused mainly on fluvial dynamics at or immediately downstream of the location where the confluent flows enter the downstream channel. This study examines in detail the spatial evolution of the time-averaged downstream velocity, cross-stream velocity, and temperature fields between the junction apex, where the flows initially meet, and the entrance to the downstream channel. A well-defined, vertically oriented mixing interface exists within this portion of the confluence, suggesting that lateral mixing of the incoming flows is limited. The downstream velocity field near the junction apex is characterized by two high-velocity cores separated by an intervening region of low-velocity or recirculating fluid. In the downstream direction, the high-velocity cores move inwards towards the mixing interface and high-velocity fluid progressively extends downwards into a zone of scour, resulting in an increase in flow velocity in the centre of the confluence. The cross-stream velocity field is dominated by flow convergence, but also includes a component associated with a consistent pattern of secondary circulation. This pattern is characterized by two surface-convergent helical cells, one on each side of the mixing interface. The helical cells appear to be the mechanism by which high-momentum fluid near the surface is advected downwards into the zone of scour. For transport-ineffective flows, the dimensions and intensities of the cells are controlled by the momentum ratio of the confluent streams and by the extant bed morphology within the confluence. Although the flow structure of formative events was not measured directly in this study, documented patterns of erosion and deposition within the central region of the confluence suggest that these events are dynamically similar to the measured flows, except for the fact that formative flows are not constrained by, but can reshape, the bed morphology. The results of this investigation are consistent with and augment previous findings on time-averaged flow structure in the downstream portion of the confluence. © 1998 John Wiley & Sons, Ltd.     

Flow structure over square bars at intermediate submergence: Large Eddy Simulation study of bar spacing effect

Turbulent open-channel flow over 2D roughness elements is investigated numerically by Large Eddy Simulation (LES). The flow over square bars for two roughness regimes (k-type roughness and transitional roughness between d-type and k-type) at a relative submergence of H/k = 6.5 is considered, where H is the maximum water depth and k is the roughness height. The selected roughness configurations are based on laboratory experiments, which are used for validating numerical simulations. Results from the LES, in turn, complement the experiments in order to investigate the time-averaged flow properties at much higher spatial resolution. The concept of the double-averaging (DA) of the governing equations is utilized to quantify roughness effects at a range of flow properties. Double-averaged velocity profiles are analysed and the applicability of the logarithmic law for rough-wall flows of intermediate submergence is evaluated. Momentum flux components are quantified and roughness effect on their vertical distribution is assessed using an integral form of the DA-equations. The relative contributions of pressure drag and viscous friction to the overall bed shear stress are also reported.     

Bed morphology in vegetated estuarine river with mild-curved meander bend

ABSTRACT

In this study, the effect of single and double row piles for reducing scouring in a mild-curved river meander was studied experimentally. The experimental study focused on the effect of vegetation on bed topography in a mild-curved meander bend. The experimental tests were conducted in a laboratory flume under clear water flow conditions. A series of experimental tests were carried out with a fixed bed and non-vegetated and vegetated moveable beds with different vegetation patterns. Analysis of the flow characteristics indicated that when the bed was mobile with vegetation on the inner bank, the core of maximum streamwise velocity shifted towards the centreline of the bend. Additionally, the cross-sectional kinetic energy increased from 0.05% for the fixed-bed test to 4.30% for the test with a double row of vegetation. Furthermore, the presence of vegetation was found to increase the uniformity of the distribution of turbulence intensity and to reduce the Reynolds shear stress along the test section. Also, the mass fluxes increased from the outer bank to the inner bank and from the upstream towards the downstream of the bend. Finally, comparison of bed topography in vegetated and non-vegetated channels showed that the maximum scour depth at the bend apex was reduced by 77% and 62% for the cases with one row and two rows of vegetation, respectively. The results of this study were compared with previously proposed models for predicting the vertical distribution of the streamwise velocity at the bend apex. It was found that Johannesson and Parker’s model (JPM) gave the lowest value of standard error. The above findings are useful in river training works and, in particular, for restoration of meandering rivers.

Editor D. M.C. Acreman; Associate editor C. Cudennec     

Does the permeability of gravel river beds affect near‐bed hydrodynamics?

The permeability of river beds is an important control on hyporheic flow and the movement of fine sediment and solutes into and out of the bed. However, relatively little is known about the effect of bed permeability on overlying near‐bed flow dynamics, and thus on fluid advection at the sediment–water interface. This study provides the first quantification of this effect for water‐worked gravel beds. Laboratory experiments in a recirculating flume revealed that flows over permeable beds exhibit fundamental differences compared with flows over impermeable beds of the same topography. The turbulence over permeable beds is less intense, more organised and more efficient at momentum transfer because eddies are more coherent. Furthermore, turbulent kinetic energy is lower, meaning that less energy is extracted from the mean flow by this turbulence. Consequently, the double‐averaged velocity is higher and the bulk flow resistance is lower over permeable beds, and there is a difference in how momentum is conveyed from the overlying flow to the bed surface. The main implications of these results are three‐fold. First, local pressure gradients, and therefore rates of material transport, across the sediment–water interface are likely to differ between impermeable and permeable beds. Second, near‐bed and hyporheic flows are unlikely to be adequately predicted by numerical models that represent the bed as an impermeable boundary. Third, more sophisticated flow resistance models are required for coarse‐grained rivers that consider not only the bed surface but also the underlying permeable structure. Overall, our results suggest that the effects of bed permeability have critical implications for hyporheic exchange, fluvial sediment dynamics and benthic habitat availability. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Bedload transport resistance in rough open‐channel flows

During bedload movement by saltation, streamwise momentum is transferred from the ?ow to the saltating grains. When the grains collide with other grains on the bed or in the ?ow, streamwise momentum is reduced, and there is a decrease in streamwise ?ow velocity and an increase in ?ow resistance, herein termed bedload transport resistance f_bt. Based on experiments in two ?umes with ?xed and mobile plane beds and previously published data, an equation is developed that may be used to predict f_bt for both capacity and non‐capacity ?ows. The variables in this equation are identi?ed by dimensional analysis and the coef?cients are determined by non‐linear regression. This equation applies to rough turbulent open‐channel ?ows, where the relative submergence is between 1 and 20 and the entire sediment load moves by saltation. An investigation of the relative magnitudes of f_bt and grain resistance f_c suggests that where dimensionless shear stress θ is less than 1 and saltation is the dominant mode of bedload transport, f_bt/f_c increases with θ but never exceeds 1. Copyright © 2004 John Wiley & Sons, Ltd.     

Detection and reconstruction of large‐scale coherent flow structures in gravel‐bed rivers

Large‐scale flow structures (LSFS) in the streamwise direction are important features of gravel‐bed river flows, because they may contribute to sediment transport and gas exchange. In the present study, these structures are detected using Huang's empirical mode decomposition and reconstructed with phase‐averaging techniques based on a Hilbert transform of the velocity signal. The analysis is based on the fluctuating component of 15 quasi‐instantaneous velocity profiles measured with a three‐dimensional (3D) acoustic Doppler velocity profiler (ADVP) in an armoured gravel‐bed river with a low relative submergence of 2.9 (ratio between flow depth and bed grain diameter). LSFS were identified in most of the measured profiles and consistently showed similar features. We were able to characterize the geometry of these large‐scale coherent structures: the front has a vertical linear shift in the time domain and a vertical profile corresponding to a first quarter moon with the apex situated at z/h ≈ 0.4. In the vertical, the front scales with flow depth h, and in the streamwise direction, LSFS scale with three to seven times the mean flow depth. On the bed, the effect of LSFS is a periodic non‐linear variation of the friction velocity on average between 0.90 and 1.10 times the mean value. A model for the friction velocity cycle resulting from LSFS oscillation is presented. Copyright © 2014 John Wiley & Sons, Ltd.     

Sensitivity of interfacial hydraulics to the microtopographic roughness of water‐lain gravels

Flow within the interfacial layer of gravel‐bed rivers is poorly understood, but this zone is important because the hydraulics here transport sediment, generate flow structures and interact with benthic organisms. We hypothesized that different gravel‐bed microtopographies generate measurable differences in hydraulic characteristics within the interfacial layer. This was tested using a high density of spatially and vertically distributed, velocity time series measured in the interfacial layers above three surfaces of contrasting microtopography. These surfaces had natural water‐worked textures, captured in the field using a casting procedure. Analysis was repeated for three discharges, with Reynolds numbers between 165000 and 287000, to evaluate whether discharge affected the impact of microtopography on interfacial flows. Relative submergence varied over a small range (3.5 to 8.1) characteristic of upland gravel‐bed rivers. Between‐surface differences in the median and variance of several time‐averaged and turbulent flow parameters were tested using non‐parametric statistics. Across all discharges, microtopographic differences did not affect spatially averaged (median) values of streamwise velocity, but were associated with significant differences in its spatial variance, and did affect spatially averaged (median) turbulent kinetic energy. Sweep and ejection events dominated the interfacial region above all surfaces at all flows, but there was a microtopographic effect, with Q2 and Q4 events less dominant and structures less persistent above the surface with the widest relief distribution, especially at the highest Reynolds number flow. Results are broadly consistent with earlier work, although this analysis is unique because of the focus on interfacial hydraulics, spatially averaged ‘patch scale’ metrics and a statistical approach to data analysis. An important implication is that observable differences in microtopography do not necessarily produce differences in interfacial hydraulics. An important observation is that appropriate roughness parameterizations for gravel‐bed rivers remain elusive, partly because the relative contributions to flow resistance of different aspects of bed microtopography are poorly constrained. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

An experimental study on the ripple–dune transition

Flume experiments were conducted on different bed stages across the ripple–dune transition. As flow velocity increases, an initially flat bed surface (made of fairly uniform sandy material) is gradually transformed into a two‐dimensional rippled bed. With further increase in velocity, two‐dimensional ripples are replaced by irregular, linguoid ripples. As the average velocity necessary for the ripple–dune transition to occur is imposed on the bed surface, these non‐equilibrium linguoid ripples are further transformed into larger, two‐dimensional dunes. For each of these stages across the transition, a concrete mould of the bed was created and the flow structure above each fixed bed surface investigated. An acoustic Doppler velocimeter was used to study the flow characteristics above each bed surface. Detailed profiles were used along a transect located in the middle of the channel. Results are presented in the form of spatially averaged profiles of various flow characteristics and of contour maps of flow fields (section view). They clearly illustrate some important distinctions in the flow structure above the different bedform types associated with different stages during the transition. Turbulence intensity and Reynolds stresses gradually increase throughout the transition. Two‐dimensional ripples present a fairly uniform spatial distribution of turbulent flow characteristics above the bed. Linguoid ripples induce three‐dimensional turbulence structure at greater heights above the bed surface and turbulence intensity tends to increase steadily with height above bed surface in the wake region. A very significant increase in turbulence intensity and momentum exchange occurs during the transition from linguoid ripples to dunes. The turbulent flow field properties above dunes are highly dependent on the position along and above the bed surface and these fields present a very high degree of spatial variability (when compared with the rippled beds). Further investigations under natural conditions emphasizing sediment transport mechanisms and rates during the transition should represent the next step of analysis, together with an emphasis on quadrant analysis. Copyright © 2001 John Wiley & Sons, Ltd.     

A three dimensional model for water-sediment movement

A three-dimensional k-ε-Ap two-fluid turbulence model is proposed to study liquid-particle two-phase flow and bed deformation.By solving coupled liquid-phase and solid-phase governing equations in a finite-volume method,the model can calculate the movement of both water and sediment.The model was validated by water-sediment transport in a 180° channel bend with a movable bed.The validation concerns two-phase time-averaged velocities,bed deformation,water depth,depth-averaged streamwise velocity,cross-stream bed profiles,and two-phase secondary flow velocity vectors.The agreement between numerical results and experimental results was generally good.The comparisons of the numerical results of different models show that the three-dimensional k-ε-Ap two-fluid turbulence model has a relatively higher accuracy than one-fluid model.     

NONUNIFORMOPEN CHANNEL FLOW WITH UPWARD SEEPAGE THROUGH LOOSE BEDS

The Reynolds stress and bed shear stress are important parameters in fluvial hydraulics. Steady-nonuniform flow in open channels with streamwise sloping beds having upward seepage through loose beds is theoretically analyzed to estimate the Reynolds stress and bed shear stress. Equations of the Reynolds stress and bed shear stress are developed, assuming a modified logarithmic velocity distribution law due to upward seepage, and using the Reynolds and continuity equations of two-dimensional flow in open channels.     

Bed roughness and grain sorting-an experimental study over fine to medium sand beds

The relation between grain-size distribution of the bed and in suspension was critically examined under a uniform flow velocity of 50 cm/s over two beds: one of mainly fine sands and the other of medium sands. Two sections – one 2.85 m downstream and the other 6.35 m downstream in the experimental channel-were selected for sampling to study the grain-sorting pattern in the vertical direction along the direction of transport. The shape and type of the grain-size distribution pattern were critically studied with height above the bed. The change in the distribution pattern has been attributed to the change of local bed roughness causing scouring against the protruded relatively coarse grains on the bed. Such trends are important to predict the nature of river bed topography. The sand of Bed-1 initially exhibits a log-skew-Laplace distribution at different heights of suspension. The distribution pattern, however, changes but this changing pattern is not consistent along the upstream side. For Bed-2, which initially exhibits a log-normal distribution, the same pattern persists from the height of suspension at 5 cm up to 20 cm. Such consistency in log-normality is also observed at the downstream points of measurement. It is generally expected that the mean grain-size would reduce with increases of suspension height but the results of the experiments, in some occasions, differ significantly from the gradual fining upward trend. This result has been attributed to local changes of bed roughness arising from the protruded relatively coarse grains causing eddies, scouring, and turbulent phenomena which moves coarse particles higher in suspension adding a coarse tail to the distribution increasing the mean grain-size.     

Flow and sedimentary processes in the meandering river South Esk,Glen Clova,Scotland

Channel geometry, flow and sedimentation in a meander bend of the River South Esk were studied from bankfull stages (January–February) to low water stages (May) in 1974. Bed topography varied little over the study period, showing a typical pool and ripple geometry. Variation of mean depth and velocity with discharge differed from section to section around the bend, due primarily to locally varying flow resistance with stage. The flow pattern for all stages was dominated by a single spiral over the point bar, with a development zone at the bend entrance. Deviation of bed shear stress from the mean flow direction was in general accord with theory, especially for high stages. The use of a uniform longitudinal water surface slope in the calculation of bed shear stress is not justified because of a complicated water surface topography, also such calculated shear may not represent effective bed shear on grains, as it accounts also for energy losses associated with secondary flows. Dunes covered much of the bar at high stages, with increasing proportions of ripples, sand ribbons and lower phase plane beds at low stages. Local flow resistance generally decreases from dunes, diminished and ripple-backed dunes, ripples, sand ribbons to plane beds, and bed forms are predicted quite well by the stream power-grain size scheme. Mean size, sorting and skewness of sediment over the bed changes little with stage. In general, size decreases, sorting improves and skewness changes from positive to negative from the talweg to the inner bank, and in the downstream direction. Allen's (1970a, b) force balance equation for moving bed load particles is supported for bankfull stage, with some reservations, and textural characteristics are explained by progressive sorting in the direction of sediment transport. Large-scale trough cross stratification (with some flat bedding) formed at high stage by dunes (and lower phase plane beds) dominates the point bar sediments. Alternations of fine-medium sand (often cross-laminated) and vegetation-rich layers result from periodic deposition on the grassed upper bar surface. Fining upwards sequences produced by lateral channel migration are modified by a coarsening upward subsequence in the upstream bar region where spiral flow is developing from the bend upstream.     

Reinvestigation on mixing length in an open channel turbulent flow

The present study proposes a model on vertical distribution of streamwise velocity in an open channel turbulent flow through a newly proposed mixing length, which is derived for both clear water and sediment-laden turbulent flows. The analysis is based on a theoretical consideration which explores the effect of density stratification on the streamwise velocity profile. The derivation of mixing length makes use of the diffusion equation where both the sediment diffusivity and momentum diffusivity are taken as a function of height from the channel bed. The damping factor present in the mixing length of sediment-fluid mixture contains velocity and concentration gradients. This factor is capable of describing the dip-phenomenon of velocity distribution. From the existing experimental data of velocity, the mixing length data are calculated. The pattern shows that mixing length increases from bed to the dip-position, having a larger value at dip-position and then decreases up to the water surface with a zero value thereat. The present model agrees well with these data sets and this behavior cannot be described by any other existing model. Finally, the proposed mixing length model is applied to find the velocity distribution in wide and narrow open channels. The derived velocity distribution is compared with laboratory channel data of velocity, and the comparison shows good agreement.     

Experimental study on flow pattern and sediment transportation at a 90°open-channel confluence

In this paper,the evolutions of flow pattern and sediment transportation at a 90° open-channel confluence with different discharge ratios (q*) of the tributary flow to the total flow were studied.The e...     

Effect of different meander curvatures on spatial variation of coherent turbulent flow structure inside ingoing multi-bend river meanders

In this paper, the effect of different curvatures on the spatial variation of coherent flow structure inside two physical models with both strongly curved and mild multi-bend meanders is investigated. Three dimensional flow velocities at three sequential meanders were measured using an Acoustic Doppler Velocity meter (Micro-ADV). Three dimensions of flow velocity are classified into two major classes and eight different bursting events. The contribution probability and transition probability of each zone is calculated from experimental data. The results indicated that the effect of curvature in sequential bends was important particularly for strongly curved bends. The contribution probability of the events for strongly curved meanders with relative curvature (R_c/B) of 2.6 were found to be higher than for mild curved meanders with relative curvature (R_c/B) of 4.43. The minimum contribution probability was found in external inward interaction event. In addition, analysis of bursting events showed that the highest values of transition probabilities occurred in the stable organizations for both models. The influences of different curvatures on distributions of the Reynolds shear stress, the turbulent kinetic energy, the streamwise velocity and the vertical velocity were also shown to be in good agreement with eroded bed. The above results can be useful for finding meandering patterns inside rivers and also in river training works.     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

A three dimensional k-ε-A_p model for water-sediment movement

A three-dimensional k-ε-Ap two-fluid turbulence model is proposed to study liquid-particle two-phase flow and bed deformation. By solving coupled liquid-phase and solid-phase governing equations in a finite-volume method, the model can calculate the movement of both water and sediment. The model was validated by water-sediment transport in a 180° channel bend with a movable bed. The validation concerns two-phase time-averaged velocities, bed deformation, water depth, depth-averaged streamwise velocity, cross-stream bed profiles, and two-phase secondary flow velocity vectors. The agreement between numerical results and experimental results was generally good. The comparisons of the numerical results of different models show that the three-dimensional k-ε-Ap two-fluid turbulence model has a relatively higher accuracy than one-fluid model.     

Spatial patterns of water surface topography at a river confluence

Understanding flow structures in river confluences has largely been the product of interpretations made from measured flow velocity data. Here, we turn the attention to the investigation of the patterns of both the average and standard deviations of the micro‐topography of the water surface at an asymmetrical natural discordant confluence for different flow conditions. Water surface topography is measured using a total station to survey the position of a reflector mounted on a custom‐built raft. To limit error problems related to changes in the water level, measurements are taken and analysed by cross‐stream transects where five water surface profiles are taken before moving to the next transect. Three‐dimensional numerical simulations of the flow dynamics at the field site are used to examine predicted water surface topography for a steady‐state situation. The patterns are interpreted with respect to flow structure dynamics, visual observations of boils, and bed topography. Results indicate that coherent patterns emerge at the water surface of a discordant bed confluence for different flow conditions. The zone of stagnation and the mixing layer are characterized by super‐elevation, a lateral tilt is present at the edge of the mixing layer, and a zone of super‐elevation is present on the tributary side at the downstream junction corner. The latter seems associated with periodical upwelling and is not present in the numerical simulations that do not take into account instantaneous velocity fluctuations. Planform curvature, topographic steering related to the tributary mouth bar, and turbulent structures associated with the mixing layer all play a key role in the pattern of both the average and standard deviation of the water surface topography at confluences. Copyright © 2002 John Wiley & Sons, Ltd.