Measurements and modeling of open-channel flows with finite semi-rigid vegetation patches

The hydrodynamics of flows through a finite length semi-rigid vegetation patch (VP) were investigated experimentally and numerically. Detailed measurements have been carried out to determine the spatial variation of velocity and turbulence profiles within the VP. The measurement results show that an intrusion region exists in which the peak Reynolds stress remains near the bed. The velocity profile is invariant within the downstream part of the VP while the Reynolds stress profile requires a longer distance to attain the spatially invariant state. Higher vegetation density leads to a shorter adjustment length of the transition region, and a higher turbulence level within the VP. The vegetation density used in the present study permits the passing through of water and causes the peak Reynolds stress and turbulence kinetic energy each the maximum at the downstream end of the patch. A 3D Reynolds-averaged Navier–Stokes model incorporating the Spalart–Allmaras turbulence closure was employed subsequently to replicate the flow development within the VP. The model reproduced transitional flow characteristics well and the results are in good agreement with the experimental data. Additional numerical experiments show that the adjustment length can be scaled by the water depth, mean velocity and maximum shear stress. Empirical equations of the adjustment lengths for mean velocity and Reynolds stress were derived with coefficients quantified from the numerical simulation results.     

Model of turbulence penetration into a suspension layer on a sediment bed and effect on vertical solute transfer

The vertical diffusional mass (solute) transfer through a suspended sediment layer, e.g. at the bottom of a lake, reservoir or estuary, by the propagation of velocity fluctuations from above was investigated. The attenuation of the velocity fluctuations in the suspension layer and the associated effect on solute transfer through the suspension layer was simulated. To represent large eddies traveling downstream in water over a high-concentration suspended sediment layer, a streamwise velocity fluctuation moving in downstream direction was imposed along the upper boundary of the suspension layer. Velocity fluctuations and downstream velocity were normalized by the shearvelocity (U_*) at the top of the suspension layer. Streamwise and vertical velocity components inside the suspension layer, were obtained from the 2-D continuity and the Navier–Stokes equations. The persistence of turbulence with depth—as it penetrates from the overlying water into the suspension layer—was found to depend on its amplitude, its period, and on the apparent viscosity of the suspension. The turbulence was found to propagate efficiently into the suspension layer when its frequency is low, and the apparent viscosity of the suspension is high. Effects on vertical mass transfer were parameterized by penetration depth and effective diffusion coefficient, and related to apparent viscosity of the suspension, Schmidt number and shear velocity on top of the suspension layer. The enhancement of turbulence penetration by viscosity is similar to the flow near an oscillating flat plate (Stokes’ second problem), but is opposite to turbulence penetration into a stationary porous and permeable sediment bed. The information is applicable to water quality modeling mear the sediment/water interface of lakes, river impoundments and estuaries.     

Spatial–temporal structure of mixing interface turbulence at two large river confluences

Converging flows at stream confluences often produce highly turbulent conditions. The shear layer/mixing interface that develops within the confluence hydrodynamic zone (CHZ) is characterized by complex patterns of three-dimensional flow that vary both spatially and temporally. Previous research has examined in detail characteristics of mean flow and turbulence along mixing interfaces at small stream confluences and laboratory junctions; however few, if any, studies have examined these characteristics within mixing interfaces at large river confluences. This study investigates the structure of mean velocity profiles as well as spatial and temporal variations in velocity, backscatter intensity, and temperature within the mixing interfaces of two large river confluences. Velocity, temperature, and backscatter intensity data were obtained at stationary locations within the mixing interfaces and at several cross sections within the CHZ using acoustic Doppler current profilers. Results show that mean flow within the mixing interfaces accelerates over distance from the junction apex. Turbulent kinetic energy initially increases rapidly over distance, but the rate of increase diminishes downstream. Hilbert–Huang transform analysis of time series data at the stationary locations shows that multiple dominant modes of fluctuations exist within the original signals of velocity, backscatter intensity, and temperature. Frequencies of the largest dominant modes correspond well with predicted frequencies for shallow wake flows, suggesting that mixing-interface dynamics include wake vortex shedding—a finding consistent with spatial patterns of depth-averaged velocities at measured cross sections. Spatial patterns of temperature and backscatter intensity show that the converging flows at both confluences do not mix substantially, indicating that turbulent structures within the mixing interfaces are relatively ineffective at producing mixing of the flows in the CHZ.     

Canopy-wake dynamics and wind sheltering effects on Earth surface fluxes

The atmospheric boundary layer adjustment at the abrupt transition from a canopy (forest) to a flat surface (land or water) is investigated in a wind tunnel experiment. Detailed measurements examining the effect of canopy turbulence on flow separation, reduced surface shear stress and wake recovery are compared to data for the classical case of a solid backward-facing step. Results provide new insights into the interpretation for flux estimation by eddy-covariance and flux gradient methods and for the assessment of surface boundary conditions in turbulence models of the atmospheric boundary layer in complex landscapes and over water bodies affected by canopy wakes. The wind tunnel results indicate that the wake of a forest canopy strongly affects surface momentum flux within a distance of 35–100 times the step or canopy height, and mean turbulence quantities require distances of at least 100 times the canopy height to adjust to the new surface. The near-surface mixing length in the wake exhibits characteristic length scales of canopy flows at the canopy edge, of the flow separation in the near wake and adjusts to surface layer scaling in the far wake. Components of the momentum budget are examined individually to determine the impact of the canopy wake. The results demonstrate why a constant flux layer does not form until far downwind in the wake. An empirical model for surface shear stress distribution from a forest canopy to a clearing or lake is proposed.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

Experimental study of recirculating flows generated by lateral shock waves in very large channels

Due to the lack of data on hydraulic-jump dynamics in very large channels, the present paper describes the main characteristics of the velocity field and turbulence in a large rectangular channel with a width of 4 m. Although a hydraulic jump is always treated as a wave that is transversal to the channel wall, in the case of this study it has a trapezoidal front shape, first starting from a point at the sidewalls and then developing downstream in an oblique manner, finally giving rise to a trapezoidal shape. The oblique wave front may be regarded as a lateral shockwave that arises from a perturbation at a certain point of the lateral wall and travels obliquely toward the centreline of the channel. The experimental work was carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). In addition to the hydraulic jump formation, a large recirculating flow zone starts to develop from the separating point of the lateral shock wave and a separate boundary layer occurs. Intensive measurements of the streamwise and spanwise flow velocity components along one-half width of the channel were taken using a bidimensional Acoustic Doppler Velocimeter (ADV). The water surface elevation was obtained by means of an ultrasonic profiler. Velocity vectors, transversal velocity profiles, turbulence intensities and Reynolds shear stresses were all investigated. The experimental results of the separated boundary layer were compared with numerical predictions and related work presented in literature and showed good agreement. The transversal velocity profiles indicated the presence of adverse pressure gradient zones and the law of the wall appears to govern the region around the separated boundary layer.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

Kelvin–Helmholtz instability in the presence of variable viscosity for mudflow resuspension in estuaries

The temporal stability of a parallel shear flow of miscible fluid layers of different density and viscosity is investigated through a linear stability analysis and direct numerical simulations. The geometry and rheology of this Newtonian fluid mixing can be viewed as a simplified model of the behavior of mudflow at the bottom of estuaries for suspension studies. In this study, focus is on the stability and transition to turbulence of an initially laminar configuration. A parametric analysis is performed by varying the values of three control parameters, namely the viscosity ratio, the Richardson and Reynolds numbers, in the case of initially identical thickness of the velocity, density and viscosity profiles. The range of parameters has been chosen so as to mimic a wide variety of real configurations. This study shows that the Kelvin–Helmholtz instability is controlled by the local Reynolds and Richardson numbers of the inflection point. In addition, at moderate Reynolds number, viscosity stratification has a strong influence on the onset of instability, the latter being enhanced at high viscosity ratio, while at high Reynolds number, the influence is less pronounced. In all cases, we show that the thickness of the mixing layer (and thus resuspension) is increased by high viscosity stratification, in particular during the non-linear development of the instability and especially pairing processes. This study suggests that mud viscosity has to be taken into account for resuspension parameterizations because of its impact on the inflection point Reynolds number and the viscosity ratio, which are key parameters for shear instabilities.     

Turbulent entrainment into fluid mud gravity currents

The entrainment of ambient water into non-Newtonian fluid mud gravity currents was investigated in this study. Constant volume release gravity currents were generated in a lock-exchange tank for a wide range of experimental conditions. A technique similar to the so-called light attenuation technique was used to find the boundary of the current, allowing for the calculation of both temporal and bulk entrainment parameters (in terms of the temporal and bulk entrainment velocities, respectively). It was found that the temporal entrainment velocity is dependent on different parameters in the different propagation phases. The slumping phase begins with an adjustment zone (henceforth, non-established zone) in which the temporal entrainment velocity is not a function of the current front velocity, followed by the established zone in which the temporal entrainment velocity is a function of the current front velocity. This dependence of the temporal entrainment velocity on the current front velocity carries through to the inertia-buoyancy phase. As expected, temporal entrainment velocity in the viscous-buoyancy phase was negligible in comparison to average entrainment velocity in the other phases. It is observed that the temporal entrainment characteristics in the non-established zone is governed by the competition between the entrainment-inhibiting density stratification effects and the entrainment-favouring effects of the Kelvin–Helmholtz billows that are quantified by the Richardson number and the Reynolds number of the gravity current, respectively. In the established zone, Reynolds number effects were observed to dominate over Richardson number effects in dictating temporal entrainment characteristics. A parameterization for the temporal entrainment velocity for non-Newtonian fluid mud gravity currents is developed based upon the experimental observations. This study also found that the bulk entrainment characteristics for the non-Newtonian fluid mud gravity currents can be parameterized by the Newtonian bulk entrainment parameterizations that rely solely on a bulk Richardson number. Interestingly, it was found that the non-Newtonian characteristics of the gravity current have little to no effect on the entrainment of the Newtonian ambient fluid.     

Stability and Mixing of a Vertical Axisymmetric Buoyant Jet in Shallow Water

The stability, mixing and effect of downstream control on axisymmetric turbulent buoyant jets discharging vertically into shallow stagnant water is studied using 3D Reynolds-averaged Navier–Stokes equations (RANS) combined with a buoyancy-extended k –ε model. The steady axisymmetric turbulent flow, temperature (or tracer concentration) and turbulence fields are computed using the finite volume method on a high resolution grid. The numerical predictions demonstrate two generic flow patterns for different turbulent heated jet discharges and environmental parameters (i) a stable buoyant discharge with the mixed fluid leaving the vertical jet region in a surface warm water layer; and (ii) an unstable buoyant discharge with flow recirculation and re-entrainment of heated water. A stratified counterflow region always appears in the far-field for both stable and unstable buoyant discharges. Provided that the domain radius L exceeds about 6H, the near field interaction and hence discharge stability is governed chiefly by the jet momentum length scale to depth ratio l_M/H, regardless of downstream control. The near field jet stability criterion is determined to be l_M/H = 3.5. A radial internal hydraulic jump always exists beyond the surface impingement region, with a 3- to 6-fold increase in dilution across the jump compared with vertical buoyant jet mixing. The predicted stability category, velocity and temperature/concentration fields are well-supported by experiments of all previous investigators.     

Embedded large eddy simulation approach for pollutant dispersion around a model building in atmospheric boundary layer

In the present article, the potential of embedded large eddy simulation (ELES) approach to reliably predict pollutant dispersion around a model building in atmospheric boundary layer is assessed. The performance of ELES in comparison with large eddy simulation (LES) is evaluated in several ways. These include a number of qualitative and quantitative comparisons of time-averaged and instantaneous results with wind tunnel measurements supplemented by statistical data analyses using scatter plots and standard evaluation metrics. Results obtained by both LES and ELES approaches show very good agreement with the experiment. However, addition of turbulence to mean flow at Reynolds averaged Navier–Stokes (RANS)–LES interface in ELES approach not only increases the turbulence intensity, it also results in larger values of turbulent kinetic energy (TKE) as well as a shorter reattachment length in the wake region. Accordingly, higher levels of TKE predicted by ELES increase the local intensity of concentration leading to shorter plume shapes as compared with LES. In general, ELES shows better agreement with experiment on the surfaces of model building and also in the downstream wake region. In terms of computational costs, the CPU time required to obtain statistical values in ELES is about 49 % lower than that of LES and the number of iterations per time step is also reduced by 55 % as compared with LES.     

Simulation of Ekman Boundary Layers by Large Eddy Model with Dynamic Mixed Subfilter Closure

Theoretical analysis of boundary layer turbulence has suggested a feasibility of sufficiently accurate turbulence resolving simulations at relatively coarse meshes. However, large eddy simulation (LES) codes, which employ traditional eddy-viscosity turbulence closures, fail to provide adequate turbulence statistics at coarse meshes especially within a surface layer. Manual tuning of parameters in these turbulence closures may correct low order turbulence statistics but severely harms spectra of turbulence kinetic energy (TKE). For more than decade, engineering LES codes successfully employ dynamic turbulence closures. A dynamic Smagorinsky turbulence closure (DSM) has been already tried in environmental LES. The DSM is able to provide adequate turbulence statistics at coarse meshes but it is not completely consistent with the LES equations. This paper investigates applicability of an advanced dynamic mixed turbulence closure (DMM) to simulations of Ekman boundary layers of high Reynolds number flows. The DMM differs from the DSM by explicit calculation of the Leonard term in the turbulence stress tensor. The Horizontal Array Turbulence Study (HATS) field program has revealed that the Leonard term is indeed an important component of the real turbulence stress tensor. This paper presents validation of a new LES code LESNIC. The study shows that the LES code with the DMM provides rather accurate low order turbulence statistics and the TKE spectra at very coarse meshes. These coarse LES maintain more energetic small scale fluctuations of velocity especially within the surface layer. This is critically important for success of simulations. Accurate representation of higher order turbulence statistics, however, requires essentially better LES resolution. The study also shows that LES of the Ekman boundary layer cannot be directly compared with conventionally neutral atmospheric boundary layers. The depth of the boundary layer is an important scaling parameter for turbulence statistics.     

Turbulence characteristics within sparse and dense canopies

Boundary layer interactions with canopies control various environmental processes. In the case of dense and homogeneous canopies, the so-called mixing layer analogy is most generally used. When the canopy becomes sparser, a transition occurs between the mixing layer and the boundary layer perturbed by interactions between element wakes. This transition has still to be fully understood and characterized. The experimental work presented here deals with the effect of the canopy density on the flow turbulence and involves an artificial canopy placed in a fully developed turbulent boundary layer. One and two-component velocity measurements are performed, both within and above the canopy. The influence of the spacing between canopy elements is studied. Longitudinal velocity statistical moments and Reynolds stresses are calculated and compared to literature data. For spacings greater than the canopy height, evidences of this transition are found in the evolution of the skewness factor, shear length scale and mixing length.     

Turbulence investigation in the roughness sub-layer of a near wall flow

The turbulence behaviour along a wall roughened by pyramidal elements was analysed in the region extending from the apex of the roughness elements up to the external limit of the roughness sub-layer. The data used for the analysis were obtained by particle image velocimetry technique. The rough wall turbulent boundary layer flow is characterized by a relatively low Reynolds number. All the results on the rough wall were compared with data referring to the canonical flow on a smooth wall turbulent boundary layer. Mean values and turbulence quantities for the two flows collapse when approaching the external limit of the roughness sublayer. The quadrant analysis of the Reynolds shear stress, in the region near the surface, shows that the contribution of the sweep motions is about equivalent for the two flows (except for wall distances lower than 40 viscous units). The contribution of the ejection motions appears to be more important over the smooth wall than over the rough wall with increasing differences approaching the wall. The probability density functions of the streamwise fluctuating velocity field for the rough wall case appear to be positively skewed in the zone very close to the pyramid apex, in contrast with the behavior observed for the smooth wall case at corresponding distances from the wall. The integral and Taylor scales for the rough wall case appear to be strongly reduced by the presence of the roughness, while the Kolmogorov microscale shows higher values.     

Near-boundary velocity and turbulence in depth-varying stream flows

This research examined the temporal distribution of turbulent structure near a streambank toe through the progression of a flood wave in West Run (Morgantown, WV, USA). Three-dimensional velocities and water depths were measured through a 17-h flood event. Turbulence characteristics were examined: Reynolds stresses, turbulent kinetic energy, and turbulence intensities. On average, near-boundary velocity during the rising stage was less than the falling stage, likely due to the measurement location and local roughness. The velocity vectors shifted from towards bed before the flood wave to toward the streambank during progression of the flood wave. Turbulent kinetic energy increased with increasing water depth during the rising stage. Reynolds stress, τ_xz, increased with increasing water depth during the rising stage, but the majority of the stresses were negative through the storm event. Reynolds stress, τ_xy, was positive throughout the event and did not vary with depth. This work is among the first to evaluate turbulence during depth-varying flows in the field.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Turbulence and aeration in hydraulic jumps: free-surface fluctuation and integral turbulent scale measurements

In an open channel, a change from a supercritical to subcritical flow is a strong dissipative process called a hydraulic jump. Herein some new measurements of free-surface fluctuations of the impingement perimeter and integral turbulent time and length scales in the roller are presented with a focus on turbulence in hydraulic jumps with a marked roller. The observations highlighted the fluctuating nature of the impingement perimeter in terms of both longitudinal and transverse locations. The results showed further the close link between the production and detachment of large eddies in jump shear layer, and the longitudinal fluctuations of the jump toe. They highlighted the importance of the impingement perimeter as the origin of the developing shear layer and a source of vorticity. The air–water flow measurements emphasised the intense flow aeration. The turbulent velocity distributions presented a shape similar to a wall jet solution with a marked shear layer downstream of the impingement point. The integral turbulent length scale distributions exhibited a monotonic increase with increasing vertical elevation within 0.2 < L_z/d₁ < 0.8 in the shear layer, where L_z is the integral turbulent length scale and d₁ the inflow depth, while the integral turbulent time scales were about two orders of magnitude smaller than the period of impingement position longitudinal oscillations.     

Turbulence structure and coherent vortices in open-channel flows with wind-induced water waves

When wind-induced water waves appear over the free-surface flows such as natural rivers and artificial channels, large amounts of oxygen gas and heat are transported toward the river bed through the interface between water and wind layers. In contrast, a bed region is a kind of turbulent boundary layer, in which turbulence generation and its transport is promoted by the production of bed shear stress. In particular, coherent hairpin vortices, together with strong ejection events toward the outer part of the layer, promote mass and momentum exchanges between the inner and outer layers. It is inferred that such a near-bed turbulence may be influenced significantly by these air–water interfacial fluctuations accompanied with free-surface velocity shear and wind-induced water waves. However, these wind effects on the wall-turbulence structure are less understood. To address these exciting and challenging topics, we conducted particle imagery velocimetry (PIV) measurements in open-channel flows combined with air flows, and furthermore the present measured data allows us to investigate the effects of air–water interactions on turbulence structure through the whole depth region.     

Mixing at a sediment concentration interface in turbulent open channel flow

In this work we address the role of turbulence on mixing of clear layer of fluid with sediment-laden layer of fluid at a sediment concentration interface. This process can be conceived as the entrainment of sediment-free fluid into the sediment-laden layer, or alternatively, as the transport of sediment into the top sediment-free flow. This process is governed by four parameters—Reynolds number of the flow \(Re_\tau\), non-dimensional settling velocity of the sediment (proxy for sediment size) \(\tilde{V}\), Richardson number \(Ri_\tau\) and Schmidt number Sc. For this work we have performed direct numerical simulations for fixed Reynolds and Schmidt numbers while varying the values of Richardson number and particle settling velocity. In the simple model considered here, the flow’s momentum and turbulence pre-exists over the entire layer of fluid, while the sediment is initially confined to a layer close to the bed. Mixing of sediment-free fluid with the sediment-laden layer is associated primarily with upward transport of sediment and buoyancy. There is no simultaneous upward transport of fluid momentum and turbulence into the sediment-free fluid layer, which is already in motion and turbulent. The analysis performed shows that the ability of the flow to transport a given sediment size decreases with the distance from the bottom, and thus only fine enough sediment particles are transported across the sediment concentration interface. For these cases, the concentration profiles evolve to a final steady state in good agreement with the well-known Rouse profile. The approach towards the Rouse profile happens through a transient self-similar state. This behavior of the flow is not seen for larger particles. Detailed analysis of the three dimensional structure of the sediment concentration interface shows the mechanisms by which sediment particles are lifted up by tongues of sediment-laden fluid with positive correlation between vertical velocity and sediment concentration. Finally, the mixing ability of the flow is addressed by monitoring the time evolution of the center of mass of the sediment-laden layer and the vertical location of the sediment-free/sediment-laden interface.     

Hydrodynamics and contaminant transport on a degraded bed at a 90-degree channel confluence

Channel confluences at which two channels merge have an important effect on momentum exchange and contaminant diffusion in both natural rivers and artificial canals. In this study, a three-dimensional numerical model, which is based on the Reynolds Averaged Navier–Stokes equations and Reynolds Stress Turbulence model, is applied to simulate and compare flow patterns and contaminant transport processes for different bed morphologies. The results clearly show that the distribution of contaminant concentrations is mainly controlled by the shear layer and two counter-rotating helical cells, which in turn are affected by the discharge ratio and the bed morphology. As the discharge ratio increases, the shear flow moves to the outer bank and the counter-clockwise tributary helical cell caused by flow deflection is enlarged, leading the mixing happens near the outer bank and the mixing layer distorted. The bed morphology can induce shrinkage of the separation zone and increase of the clockwise main channel helical cell, which is initiated by the interaction between the tributary helical cell and the main channel flow and strengthened by the deep scour hole. The bed morphology can also affect the distortion direction of the mixing layer. Both a large discharge ratio and the bed morphology could lead to an increase in mixing intensity.