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.     

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.     

Direct numerical simulations of boundary condition effects on the propagation of density current in wall-bounded and open channels

The propagation of density current under different boundary conditions is investigated using high resolution direct numerical simulations (DNS). A revised Kleiser and Schumann influence-matrix method is used to treat the general Robin type velocity boundary conditions and the related “tau” error corrections in the numerical simulations. Comparison of the simulation results reveals that the boundary conditions change the turbulent flow field and therefore the propagation of the front. This paper mainly focuses on the effects of boundary conditions and initial depth of the dense fluid. The differences in energy dissipation and overall front development in wall-bounded and open channels are examined. Through DNS simulations, it is evident that with the decrease of initial release depth ratio ( $D/H$ ), the effect of the top boundary becomes less important. In wall-bounded channels, there are three distinctive layers in the vertical distribution of energy dissipation corresponding to the contributions from bottom wall, interface, and top wall, respectively. In open channels, there are only two layers with the top one missing due to the shear free nature of the boundary. It is found that the energy dissipation distribution in the bottom layer is similar for cases with the same $D/H$ ratio regardless the top boundary condition. The simulation results also reveal that for low Reynolds number cases, the energy change due to concentration diffusion cannot be neglected in the energy budget. To reflect the real dynamics of density current, the dimensionless Froude number and Reynolds number should be defined using the release depth $D$ as the length scale.     

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.     

A closure approach with multi-scale mixing length,and its application to wind and shear stress profiles within the atmospheric boundary layer

A method to determine flow specific first-order closure for the turbulent flux of momentum in the atmospheric boundary layer (ABL) is presented. This is based on the premise that eddy viscosity is a flow rather than a fluid property, and the physically more realistic assumption that the transfer of momentum and other scalar quantities in a turbulent flow takes place by a large, but finite number of length scales, than the often used single length scale, the ‘mixing length’. The resulting eddy viscosity is flow specific and when applied to the study of the ABL, yields the vertical profiles of shear stress and mean wind velocity in good agreement with observations. The method may be extended to other types of turbulent flows, however it should be recognized that each type of flow may yield a different eddy viscosity profile. Using the derived eddy viscosity the paper presents simple analytical solutions of the ABL equations to determine observationally consistent wind speed and shear stress profiles in the ABL for a variety of practical applications including air pollution modelling.     

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.     

The Structure of the Shear Layer in Flows over Rigid and Flexible Canopies

Flume experiments were conducted with rigid and flexible model vegetation to study the structure of coherent vortices (a manifestation of the Kelvin–Helmholtz instability) and vertical transport in shallow vegetated shear flows. The vortex street in a vegetated shear layer creates a pronounced oscillation in the velocity profile, with the velocity near the top of a model canopy varying by a factor of three during vortex passage. In turn, this velocity oscillation drives the coherent waving of flexible canopies. Relative to flows over rigid vegetation, the oscillation in canopy geometry has the effect of decreasing the amount of turbulent vertical momentum transport in the shear layer. Using a waving plant to determine phase in the vortex cycle, each vortex is shown to consist of a strong sweep at its front (during which the canopy is most deflected), followed by a weak ejection at its rear (when the canopy height is at a maximum). Whereas in unobstructed mixing layers the vortices span the entire layer, they encompass only 70% of the flexibly obstructed shear layer studied here.     

Analysis of the velocity field in a large rectangular channel with lateral shockwave

In this work the authors describe the main characteristics of the velocity field of hydraulic jumps in a very large channel where lateral shockwaves occur. Experiments were carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). Extensive flow velocity measurements were investigated in order to have a clearer understanding of both hydraulic jump development and lateral shockwave formation in a very large channel. Eight experiments were performed in a 4m wide rectangular channel; the experiments differed in the inlet Froude number F ₀ and the jump type. Seven tests were carried out with undular jumps and one with a roller jump. The flow velocity and the flow free surface measurements were taken using a two-dimensional Acoustic Doppler Velocimeter (ADV) and an ultrasonic profiler, respectively. The experimental results can be summarized as follow: (i) the formation of well developed lateral shockwaves similar to those of oblique jumps were observed; (ii) the comparison of the experimental and theoretical data shows that the classic shockwave theory is sufficiently confirmed in the analyzed range of Reynolds number, taking into account the experimental errors and the difference between the theoretical and experimental assumptions; (iii) the transversal flow velocity profiles in the recirculating zone show a good agreement with the numerical simulations presented in literature in the case of a separated turbulent boundary layer over a flat plate. This conclusion enables us to confirm the hypothesis that the lateral shockwaves in the channel are the result of a boundary layer which, as observed, forms on the channel sidewalls.     

Estimation of far-field horizontal and vertical turbulent diffusion coefficients from the concentration field of a wastewater plume near the Akashi Strait

We carried out a field study of the plume discharged by a near-shore wastewater outfall near the Akashi Strait, Japan. Using an Acoustic Doppler Current Profiler and a tow-body CTD, we measured the near-surface salinity and temperature fields in the region throughout an M₂ tidal cycle. We filtered the data in T–S space to remove water masses other than the wastewater, and then used the adiabatic mixing assumption to calculate the concentration of wastewater in the far field of this plume. Averaging the T–S fields of repeated surveys over a time period during which the tidal regime did not change substantially, allowed comparison of the time-averaged plume with the analytical solution for a plume diffusing in both the horizontal and vertical dimensions. The resulting vertical turbulent diffusion coefficients agreed well with those resulting from Thorpe scales determined via a vertically-profiling CTD, as well as with the canonical value for open channel flow of D _z = 0.067hu _*. The corresponding horizontal turbulent diffusion coefficients, however, were two orders of magnitude larger than those typically observed in straight channels, and an order of magnitude larger than those observed in meandering rivers. This is likely a result of enhanced horizontal mixing due to barotropic eddies generated by the interaction of strong tidal flow with headlands and levees, as well as due to the time-varying nature of tidal flow, and baroclinic spreading of the buoyant wastewater plume.     

Numerical Simulation of Interaction of the Heavy Gas Cloud with the Atmospheric Surface Layer

The numerical time-dependent three-dimensional model [Kovalets, I.V. and Maderich, V.S.: 2001, Int. J. Fluid Mech. Res. 30, 410–429] of the heavy gas dispersion in the atmospheric boundary layer has been improved by parameterizing momentum and heat fluxes on the surface of Earth using Monin–Obukhov similarity theory. Three parameterizations of heat exchange with the surface of Earth were considered: (A) formula of Yaglom A.M. and Kader B.A. [1974, J. Fluid Mech. 62, 601–623] for forced convection, (B) interpolation formula for mixed convection and (C) similarity relationship for mixed convection [Kader, B.A. and Yaglom, A.M.: 1990, J. Fluid Mech. 212, 637–662]. Two case studies were considered. In the first study based on experiment of Zhu et al., J. Hazard Mater 62, 161–186], the interaction of an isothermal heavy gas plume with an atmospheric surface layer was simulated. It was found that stable stratification in the cloud essentially suppresses the turbulence in the plume, reducing the turbulent momentum flux by a factor of down to 1/5 in comparison with the undisturbed value. This reduction essentially influences velocities in the atmospheric boundary layer above the cloud, increasing the mean velocity by a factor of up to 1.3 in comparison with the undisturbed value. A simulation of cold heavy gas dispersion was carried out in the second case based on field experiment BURRO 8. It was shown that both forced and free convections under moderate wind speeds significantly influence the plume. The relative rms and bias errors of prediction the plume’s height were σ_H ≈ 30% and ɛ_H = − 10%, respectively, for parameterization B, while for A and C the errors were σ_H ≈ 80% and ɛ_H ≈ − 65%. It is therefore advised to use the simple parameterization B in dense gas dispersion models.     

Flow,flapping, and photosynthesis ofNereocystis leutkeana: a functional comparison of undulate and flat blade morphologies

A number of species of macroalagae possess a flat, strap-like blade morphology in habitats exposed to rapidly-moving water whereas those at protected sites have a wider, undulate blade shape. We have explored the functional consequences of flat, narrow vs. wide, undulate blade morphologies in the giant bull kelpNereocystis luetkeana. Our study focused on the behavior of blades in ambient water currents and the consequences of that behavior to breakage and to photosynthesis. In flowing water, the narrow, flat blades flap with lower amplitude and collapse together into a more streamlined bundle than do wide, undulate blades, and hence experience lower drag per blade area at a given flow velocity. If the algae at current-swept sites had ruffled blades, drag forces would sometimes be sufficient to break the stipes. However, flat blades in a streamlined bundle experience more self-shading than do undulate blades, which remain spread out in water currents. Thus, there is a morphological trade-off between reducing drag and reducing self-shading. Photosynthetic¹⁴C-HCO₃ uptake rates decrease in slow flow when the boundary layer along the blade surface across which diffusion takes place is relatively thick. However, blade flapping, which stirs water near the blade surface, enhances carbon uptake rates in slow water currents for both the undulate and the flat morphologies.     

A Numerical Study for Turbulent Flow and Thermal Influence over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier–Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052–2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349–377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369–398)’s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243–258)’s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.     

Large-eddy simulation of low-level jet-like flow in a canopy

The computational method of Large-Eddy Simulations has been used to study the weak, neutrally stable drainage flow within tree canopies. The computational results show that a secondary velocity maximum that resembles a jet is formed within the canopy under the nocturnal flow conditions. This jet-like flow is important in the analysis and measurements of the net ecosystem-atmosphere exchange (NEE) for carbon dioxide (CO₂). A uniformly distributed, plane source was placed within the canopy in order to simulate the nocturnal production of CO₂. The NEE is calculated as the sum of the integration of the rate of change of the concentration of CO₂ over the computational domain, the vertical turbulent flux measured directly by eddy-covariance (EC) method, and the advection terms, which are not taken into account in the EC method. Numerical results of the velocity and concentration fields, within and above the canopy, are presented and their impact on the CO₂ transport is investigated in detail. The computational results show that 15–20% of NEE is drained out by the advection process under the canopy. The results also show that the turbulent fluctuations in the lateral direction are also significant and may result in 2–5% CO₂ transport.     

Energy dissipation estimates in oscillating grid setup: LDV and PIV measurements

The dissipation of turbulent kinetic energy has been increasingly used as a scaling parameter to integrate microbiological accrual and metabolic rates with fluid-flow motion in natural and engineered aquatic ecosystems. The estimation of turbulent kinetic energy under field conditions and the generation of energy dissipation rates under controlled laboratory conditions with microbiological organisms are necessities required to integrate environmental/ecological laboratory protocols with a moving fluid in the environment. Turbulent fluid-flow conditions were generated in an oscillating grid setup, and turbulence variables were quantified using laser-Doppler velocimetry (LDV) and particle image velocimetry (PIV) measuring techniques. The rate of dissipation of the turbulent kinetic energy in the setup ranged from 10⁻⁹ to 10⁻⁴ m²/s³ and was similar to the levels of energy dissipation commonly reported in engineered and natural aquatic ecosystems. Time-averaged velocities were close to zero with the root-mean-square velocity ratios about 1, indicating nearly isotropic fluid-flow conditions in the setup. The velocity spectra, obtained by stationary LDV measurements for the vertical and horizontal velocity components across the setup revealed the existence of inertial subrange with the frequency power scaling law of “ω ^−5/3.” The estimated Eulerian frequency spectrum followed the theoretical functional relation and confirmed the applicability of inertial dissipation method for the estimation of turbulent kinetic energy dissipation rates. PIV was used for a direct estimation of dissipation by evaluating spatially distributed velocity gradients. The direct dissipation estimate in conjunction with the estimated Eulerian frequency spectrum provided evaluation of a “universal” constant, α, commonly used for the estimation of an energy dissipation rate over the inertial subrange of the Eulerian spectrum. The results demonstrated a range of values, rather than a universal constant, of α with a lognormal probability distribution for vertical and horizontal velocity components. In order to encompass a 0.955 probability range under the lognormal distribution the universal constant, α, should be in the range 2.91 ≥ α _u ≥ 0.43 and 4.44 ≥ α _w ≥ 0.42 for horizontal and vertical velocity components, respectively.     

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.     

Streamwise velocity profiles in coastal currents

The ability to model marine currents can be a powerful device for many planning activities, for which the knowledge of the velocity field is of pivotal importance, such as the evaluation of current-induced loading on maritime structures or the diffusion and dispersion of polluted flow discharges. Observations of time-averaged velocity profiles, taken with a vessel mounted acoustic Doppler current profiler during a monitoring survey program in the seas of Southern Italy, are analysed in this paper. The measurements were taken under non-breaking conditions, offshore the surf zone, with the aim of reproducing the vertical trends of the streamwise velocity by means of standard theoretical laws. To evaluate also the possible influence of stratification on the current velocity profile shape, together with velocity measurements water temperature and salinity were also measured at the same time and locations, by means of a CTD recorder. The examined surveys referred to different time periods and sites, to guarantee a general validity of deductions. On the basis of the experiments, we verified the actual existence of a log-layer and concluded that the upper limit of the region in which the log law is applicable extends well beyond the inner region. Moreover, the deviations of the measured velocity from the logarithmic profiles above the height of the log layer is consistent with the effects of stratification. The parameters of the log law were estimated, depending on both flow dynamics and stratification in the target area. As a second step, in the most superficial and stratified layer, the velocity profiles were modelled by means of a power law, which fitted the measured data well. According to previous studies, the power law parameters result Reynolds number dependent by means of a new proposed formulation. Finally, the bottom stress and the bottom drag coefficient were investigated.     

Comments on “turbulent velocity profile in fully-developed open channel flows”

In a recent paper published in this journal, Bonakdari et al. (Environ Fluid Mech 8:1–17, 2008) presented a new formulation of the vertical velocity profile in the central portion of steady fully developed turbulent open-channel flows which is based on an analysis of the Navier–Stokes equations. The predicted mean streamwise velocity profiles presented by the authors represent dip phenomenon (maximum velocity below the free surface). The discusser would like to point out some contradictions and an error in the main demonstration. This error appears in the right-hand side term of the main Eq. 22 (Environ Fluid Mech 8:1–17, 2008) in the parameter α. Our corrected demonstration shows that α is equal to and not to as obtained by the authors (Environ Fluid Mech 8:1–17, 2008). The values of α used by the authors are under-estimated by 1, this difference has a significant effect on velocities and therefore will involve a difference between measured and predicted velocity profiles. Finally, the assumption of at the free surface is in contradiction with predicted velocity profiles.     

Different proxies for the reactivity of aquatic sediments towards oxygen: A model assessment

Information on benthic carbon mineralization rates is often derived from the analysis of oxygen microprofiles in sediments. To enable a direct comparison of different sediment environments, it is often desirable to characterize sediments by a single proxy that expresses their “reactivity” towards oxygen. For this, there are three commonly used proxies: the oxygen penetration depth (OPD), the oxygen flux at the sediment-water interface (DOU), and the maximum volumetric oxygen consumption rate (R_max). The OPD can be directly determined from the oxygen depth profile, while the DOU is usually obtained by a linear fit to the oxygen gradient either in diffusive boundary layer. The oxygen consumption rate R_max requires the fitting of a reactive-transport model to the data profile. This article shows that the OPD alone is a suboptimal proxy, because it shows a strong dependence on the half-saturation constant K_s, and secondly, because it is sensitive to the particular re-oxidation conditions right above the oxic-anoxic interface. Similarly, the volumetric oxygen consumption rate R_max is rather strongly dependent on the kinetic model formulation employed. To show this we fitted three different (Bouldin, Blackman and Monod) kinetics to the same oxygen data profiles. When fitting these models, the R_max values obtained differed by 20% for exactly the same oxygen profile. Accordingly, if one reports R_max values, it is crucial to specify the kinetic model alongside. Overall, DOU emerges as sediment reactivity proxy which is the least model dependent.