Markovian–Octant analysis based stable turbulent shear stresses in near-bed bursting phenomena of vortex settling chamber

The coherent structure in near-bed turbulent boundary layer of vortex chamber, particularly the bursting events and their associated shear stresses play the main role in sediment flushing process and consequently the trap efficiency of the vortex settling chamber. Hence, three-dimensional velocity measurements were made at 48 points near the bed of physical model of vortex chamber by using Micro-ADV. The pattern of sediment deposition at the bed of vortex settling chamber reveals three separate regions formed by three predominant currents of inlet flow, flushing flow and outlet over flow. Additionally, due to the instability and three-dimensional nature of the bursting events near the bed of chamber, the new method of Markovian–Octant analysis was applied to study the different classes of near-bed stable shear stresses of vortex chamber in three dimensions. Moreover, the role of each class of stable shear stresses on Sediment transport mechanism at the bed of vortex chamber is investigated.     

Evaluating the dimensionality and significance of “active periods” in turbulent environmental flows defined using Lipshitz/Hölder regularity

Active periods within perturbed boundary-layer flows are considered in terms of the local roughness of measured velocity time series and defined in terms of Hölder/Lipshitz exponents. Such events are associated with the passage of energetic, coherent flow structures and are responsible for exerting high turbulent stresses because of the rapid changes in velocity that occur at such times. A method is proposed for assessing the effective dimensionality of such active periods, as well as their significance to the flow field, for a particular choice of flow metric. The method is applied to the turbulent flow through a confluence flow geometry, with velocity samples acquired close to the bed of the channel in a zone of complex mixing. The dimensionality of the active periods is consistent with the observed patterns of sediment entrainment from the bed, with the significance of the active periods decaying away from the erosional zone.     

Pier scour within long contraction in cohesive sediment bed

In the present study, an attempt has been made to study the bridge pier scour embedded long contraction for the bed sediment prepared by a mixture of clay and fine sand with varying proportions having a specific range of antecedent moisture content. Results particularly focused on the clay–sand mixed cohesive bed at varying clay fractions of the sediment bed, approach flow velocity, contraction ratio and different pier shapes, on maximum equilibrium scour depth for pier scour within long contraction. Further, regression based equations for the estimation of non-dimensional maximum scour depth for piers within long contraction in clay–sand mixed cohesive bed embedded were proposed.     

Turbulent burst-sweep events around fully submerged vertical square cylinder over plane bed

Environmental Fluid Mechanics - The subtle variations of the turbulent shear stress affect the entrainment of bed sediments into the flow as well as sediment transport around submerged...     

Characterization of bedload intermittency near the threshold of motion using a Lagrangian sediment transport model

At the smallest scales of sediment transport in rivers, the coherent structures of the turbulent boundary layer constitute the fundamental mechanisms of bedload transport, locally increasing the instantaneous hydrodynamic forces acting on sediment particles, and mobilizing them downstream. Near the critical threshold for initiating sediment motion, the interactions of the particles with these unsteady coherent structures and with other sediment grains, produce localized transport events with brief episodes of collective motion occurring due to the near-bed velocity fluctuations. Simulations of these flows pose a significant challenge for numerical models aimed at capturing the physical processes and complex non-linear interactions that generate highly intermittent and self-similar bedload transport fluxes. In this investigation we carry out direct numerical simulations of the flow in a rectangular flat-bed channel, at a Reynolds number equal to Re = 3632, coupled with the discrete element method to simulate the dynamics of spherical particles near the bed. We perform two-way coupled Lagrangian simulations of 48,510 sediment particles, with 4851 fixed particles to account for bed roughness. Our simulations consider a total of eight different values of the non-dimensional Shields parameter to study the evolution of transport statistics. From the trajectory and velocity of each sediment particle, we compute the changes in the probability distribution functions of velocities, bed activity, and jump lengths as the Shields number increases. For the lower shear stresses, the intermittency of the global bedload transport flux is described by computing the singularity or multifr actal spectrum of transport, which also characterizes the widespread range of transport event magnitudes. These findings can help to identify the mechanisms of sediment transport at the particle scale. The statistical analysis can also be used as an ingredient to develop larger, upscaled models for predicting mean transport rates, considering the variability of entrainment and deposition that characterizes the transport near the threshold of motion.     

Estimation of maximum scour depths at upstream of front and rear piers for two in-line circular columns

Previous investigations indicate that scour around bridge piers is one of the most important factors for the failure of waterway bridges. Hence, it is essential to determine the accurate scour depth around the bridge piers. Most of the previous studies were based on scour around a single pier; however, in practice, new bridges are usually wide and then piers comprise two circular piers aligned in the flow direction that together support the loading of the structure. In this study, the effect on maximum scour depth of the spacing between two piers aligned in the flow direction was investigated experimentally under clear water scour conditions. The results show that the maximum scour depth at upstream of the front pier occurs when the spacing between the two piers is 2.5 times the diameter of the pier. Two semi empirical equations have been developed to predict the maximum scour depth at upstream of both front and rear piers as a function of the spacing between the piers, in terms of a pier-spacing factor. If the new equations for the pier-spacing factor are used with some of the existing equations for scour at a single pier, the predicted scouring depths are in good agreement with observed results. The S/M equation exhibited the best performance among the various equations tested and was recommended for use in prediction of the equilibrium scour depth. The findings of this study can be used to facilitate the positioning of piers when scouring is a design concern.     

Variations of bed elevations due to turbulence around submerged cylinder in sand beds

This paper presents the spatio-temporal variations in bed elevations and the near-bed turbulence statistics over the deformed bed generated around the submerged cylindrical piers embedded vertically on loose sediment bed at a constant flow discharge. Experiments were carried out in a laboratory flume for three blockage ratios in the range of 0.04–0.06 using three different sizes of submerged cylinders individually placed vertically at the centerline of the flume. Clear-water experimental conditions were maintained over the smooth sediment bed surface with a constant flow discharge (\(Q = 0.015\,{\rm m}^3/{\rm sec}\)), thereby giving three different cylinder Reynolds numbers \(Re_{D_c} = \frac{U_mD_c}{\nu }\) (=10200, 12750, 15300) away from the cylinder locations, where \(U_m\) is the maximum mean velocity, \(D_c\) is the cylinder diameter and \(\nu\) is the kinematic viscosity of fluid. Instantaneous sand bed elevations around the cylinders were recorded using a SeaTek 5MHz ultrasonic ranging system of net 24 transducers to estimate bed form migration, and the near-bed velocity data at transducer locations over the stable deformed bed around the pier-like structures were collected using down-looking three-dimensional (3D) Micro-acoustic Doppler velocimeter to estimate the bottom Reynolds shear stresses and the contributions of bursting events to the dominant shear stress component. The flow perturbation generated due to relatively lower flow blockage ratio favored to achieve the stable bed condition more rapidly than the others, and larger upstream scour-depth and deformed areas were noticed for greater flow blockage ratio due to larger cylinder diameter. For larger blockage ratio in the upstream of scour-hole near the bed, occurrences of probabilities of both boundary-ward interactions (Q1 and Q3) were the dominant; whereas in the downstream of the scoured region, occurrences of probabilities of second and third quadrant events (Q2 and Q4) were dominant. On the other hand, for the lower blockage ratio, quadrant (Q2) was dominant over Q4 in the downstream of scour-hole, and in the upstream of scour-hole, quadrant Q4 was the dominant.     

Hydrodynamic characteristics and related mass-transfer properties in open-channel flows with rectangular embayment zone

Consecutive groynes and embayments form dead water zones, where sedimentation and high concentrations of pollutants are often observed. It is thus very important to understand the mass and momentum exchange between the main channel and side cavities in rivers and hydraulic engineering structures. The spanwise gradient of the streamwise velocity near the junction produces small-scale turbulent vortices because of shear instability. Furthermore, large-scale horizontal circulation is also generated in the cavity zone. These coherent turbulent structures play a significant role in mass and sediment transfer at the boundary between the mainstream and embayment. However, the relation between turbulence and mass transfer is poorly understood. In this study, we performed particle image velocity and laser-induced fluorescence experiments using a laboratory flume, laser light sheets and a high-speed CMOS camera. We examined the exchange properties of a dye as a function of bed configuration and sedimentation effect. Both primary and secondary gyres were observed in the flat bed and downward-sloping bed, whereas the primary gyre was prevalent in the upward-sloping bed. Moreover, the horizontal circulation strongly affected the mass-transfer properties between the mainstream and side cavity.     

Simulation-based optimization of in-stream structures design: rock vanes

We employ a three-dimensional coupled hydro-morphodynamic model, the Virtual Flow Simulator (VFS-Geophysics) in its Unsteady Reynolds Averaged Navier–Stokes mode closed with \(k-\omega\) model, to simulate the turbulent flow and sediment transport in large-scale sand and gravel bed waterways under prototype and live-bed conditions. The simulation results are used to carry out systematic numerical experiments to develop design guidelines for rock vane structures. The numerical model is based on the Curvilinear Immersed Boundary approach to simulate flow and sediment transport processes in arbitrarily complex rivers with embedded rock structures. Three validation test cases are conducted to examine the capability of the model in capturing turbulent flow and sediment transport in channels with mobile-bed. Transport of sediment materials is handled using the Exner equation coupled with a transport equation for suspended load. Two representative meandering rivers, with gravel and sand beds, respectively, are selected to serve as the virtual test-bed for developing design guidelines for rock vane structures. The characteristics of these rivers are selected based on available field data. Initially guided by existing design guidelines, we consider numerous arrangements of rock vane structures computationally to identify optimal structure design and placement characteristics for a given river system.     

Spatial evolution of coherent motions in finite-length vegetation patch flow

A number of experimental studies on submerged canopy flows have focused on fully-developed flow and turbulent characteristics. In many natural rivers, however, aquatic vegetation occurs in patches of finite length. In such vegetated flows, the shear layer is not formed at the upstream edge of the vegetation patch and coherent motions develop downstream. Therefore, more work is neededz to reveal the development process for large-scale coherent structures within vegetation patches. For this work, we considered the effect of a limited length vegetation patch. Turbulence measurements were intensively conducted in open-channel flows with submerged vegetation using Particle Image Velocimetry (PIV). To examine the transition from boundary-layer flow upstream of the vegetation patch to a mixing-layer-type flow within the patch, velocity profiles were measured at 33 positions in a longitudinal direction. A phenomenological model for the development process in the vegetation flow was developed. The model decomposed the entire flow region into four zones. The four zones are the following: (i) the smooth bed zone, (ii) the diverging flow zone, (iii) the developing zone and (iv) the fully-developed zone. The PIV data also confirmed the efficiency of the mixing-layer analogy and provided insight into the spatial evolution of coherent motions.     

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.     

Large eddy simulation of flow and scalar transport in a vegetated channel

Predicting flow and mass transport in vegetated regions has a broad range of applications in ecology and engineering practice. This paper presents large eddy simulation (LES) of turbulent flow and scalar transport within a fully developed open-channel with submerged vegetation. To properly represent the scalar transport, an additional diffusivity was introduced within the canopy to account for the contribution of stem wakes, which were not resolved by the LES, to turbulent diffusion. The LES produced good agreement with the velocity and concentration fields measured in a flume experiment. The simulation revealed a secondary flow distributed symmetrically about the channel centerline, which differed significantly from the circulation in a bare channel. The secondary circulation accelerated the vertical spread of the plume both within and above the canopy layer. Quadrant analysis was used to identify the form and shape of canopy-scale turbulent structures within and above the vegetation canopy. Within the canopy, sweep events contributed more to momentum transfer than ejection events, whereas the opposite occurred above the canopy. The coherent structures were similar to those observed in terrestrial canopies, but smaller in scale due to the constraint of the water surface.     

Efficient numerical computation and experimental study of temporally long equilibrium scour development around abutment

For the abutment bed scour to reach its equilibrium state, a long flow time is needed. Hence, the employment of usual strategy of simulating such scouring event using the 3D numerical model is very time consuming and less practical. In order to develop an applicable model to consider temporally long abutment scouring process, this study modifies the common approach of 2D shallow water equations (SWEs) model to account for the sediment transport and turbulence, and provides a realistic approach to simulate the long scouring process to reach the full scour equilibrium. Due to the high demand of the 2D SWEs numerical scheme performance to simulate the abutment bed scouring, a recently proposed surface gradient upwind method (SGUM) was also used to improve the simulation of the numerical source terms. The abutment scour experiments of this study were conducted using the facility of Hydraulics Laboratory at Nanyang Technological University, Singapore to compare with the presented 2D SGUM–SWEs model. Fifteen experiments were conducted with their scouring flow durations vary from 46 to 546 h. The comparison shows that the 2D SGUM–SWEs model gives good representation to the experimental results with the practical advantage.     

Experimental assessment of characteristic turbulent scales in two-phase flow of hydraulic jump: from bottom to free surface

A hydraulic jump is a turbulent shear flow with a free-surface roller. The turbulent flow pattern is characterised by the development of instantaneous three-dimensional turbulent structures throughout the air–water column up to the free surface. The length and time scales of the turbulent structures are key information to describe the turbulent processes, which is of significant importance for the improvement of numerical models and physical measurement techniques. However, few physical data are available so far due to the complexity of the measurement. This paper presents an investigation of a series of characteristic turbulent scales for hydraulic jumps, covering the length and time scales of turbulent flow structures in bubbly flow, on free surface and at the impingement point. The bubbly-flow turbulent scales are obtained for Fr = 7.5 with 3.4 × 10⁴ < Re < 1.4 × 10⁵ in both longitudinal and transverse directions, and are compared with the free-surface scales. The results highlight three-dimensional flow patterns with anisotropic turbulence field. The turbulent structures are observed with different length and time scales respectively in the shear flow region and free-surface recirculation region. The bubbly structures next to the roller surface and the free-surface fluctuation structures show comparable length and time scales, both larger than the scales of vortical structures in the shear flow and smaller than the scales of impingement perimeter at the jump toe. A decomposition of physical signals indicates that the large turbulent scales are related to the unsteady motion of the flow in the upper part of the roller, while the high-frequency velocity turbulence dominates in the lower part of the roller. Scale effects cannot be ignored for Reynolds number smaller than 4 × 10⁴, mainly linked to the formation of large eddies in the shear layer. The present study provides a comprehensive assessment of turbulent scales in hydraulic jump, including the analyses of first data set of longitudinal bubbly-flow integral scales and transverse jump toe perimeter integral scales.     

Small-scale entrainment in inclined gravity currents

We investigate the effect of buoyancy on the small-scale aspects of turbulent entrainment by performing direct numerical simulation of a gravity current and a wall jet. In both flows, we detect the turbulent/nonturbulent interface separating turbulent from irrotational ambient flow regions using a range of enstrophy iso-levels spanning many orders of magnitude. Conform to expectation, the relative enstrophy isosurface velocity \(v_n\) in the viscous superlayer scales with the Kolmogorov velocity for both flow cases. We connect the integral entrainment coefficient E to the small-scale entrainment and observe excellent agreement between the two estimates throughout the viscous superlayer. The contribution of baroclinic torque to \(v_n\) is negligible, and we show that the primary reason for reduced entrainment in the gravity current as compared to the wall-jet are 1) the reduction of \(v_n\) relative to the integral velocity scale \(u_T\); and 2) the reduction in the surface area of the isosurfaces.     

Large Eddy Simulation of turbulence generated by a weak breaking tidal bore

A tidal bore is a natural and fragile phenomenon, which is of great importance for the ecology of an estuary. The bore development is closely linked with the tidal range and the river mouth shape, and its existence is sensitive to any small change in boundary conditions. Despite their ecological and cultural value, little is known on the flow field, turbulent mixing and sediment motion beneath tidal bores. Indeed, some striking features can be highlighted in two-dimensional simulations, such as large velocity fluctuations and flow recirculation structures. Using Large Eddy Simulation method, the numerical results emphasised the complicated turbulent structures and their unsteadiness under a tidal bore.     

Turbulence and turbulent flux events in a small estuary

Relatively little systematic research has been conducted on the turbulence characteristics of small estuaries. In the present study, detailed measurements were conducted in a small subtropical estuary with a focus on turbulent flux events. Acoustic Doppler velocimeters were installed in the mid-estuary at fixed locations and sampled simultaneously and continuously for 50 h. A turbulent flux event analysis was performed for the entire data sets extending the technique of Narasimha et al. (Phil Trans R Soc Ser A 365:841–858, 2007) to the unsteady open channel flow motion and to turbulent sub-events. Turbulent bursting events were defined in terms of the instantaneous turbulent flux. The data showed close results for all ADV units. The very-large majority of turbulent events lasted between 0.04 and 0.3 s with an average of 1 to 4 turbulent events observed per second. A number of turbulent bursting events consisted of consecutive turbulent sub-events, with between 1 and 3 sub-events per main event on average. For all ADV systems, the number of events, event duration and event amplitude showed some tidal trends, with basic differences between high- and low-water periods. A comparison between the present estuary data and the atmospheric boundary layer results of Narasimha et al. (Phil Trans R Soc Ser A 365:841–858, 2007) showed a number of similarities and demonstrated the significance of turbulent events in environmental flows. A burstiness index of 0.85 was found for the present data.     

Simulation-based optimization of in–stream structures design: bendway weirs

Bendway weirs are one of the most practical in–stream rock structures utilized to protect the outer bend of meandering streams and rivers from erosion. We present development of a simulation-based paradigm for effective design of bendway weir structures to enhance meandering stream bank stability and control lateral migration. To do so, we employ the St. Anthony Falls Laboratory Virtual StreamLab (VSL3D) code to elucidate the flow and sediment transport phenomena induced by interaction of flow, mobile bed, and in–stream structures in large rivers under prototype conditions. We carried out numerous numerical experiments to systematically simulate various arrangements of bendway weir in two river test-beds and gaining insights into the physical mechanisms via which such bendway weirs modify turbulent flow, sediment transport and scour processes. The so-gained physical insights are then taken into account to develop a set of practical physics-based design criteria for optimal placement of bendway weirs in large rivers.     

Turbulence measurements in a small subtropical estuary under king tide conditions

In natural waterways and estuaries, the understanding of turbulent mixing is critical to the knowledge of sediment transport, stormwater runoff during flood events, and release of nutrient-rich wastewater into ecosystems. In the present study, some field measurements were conducted in a small subtropical estuary with micro-tidal range and semi-diurnal tides during king tide conditions: i.e., the tidal range was the largest for both 2009 and 2010. The turbulent velocity measurements were performed continuously at high-frequency (50Hz) for 60?h. Two acoustic Doppler velocimeters (ADVs) were sampled simultaneously in the middle estuarine zone, and a third ADV was deployed in the upper estuary for 12?h only. The results provided an unique characterisation of the turbulence in both middle and upper estuarine zones under the king tide conditions. The present observations showed some marked differences between king tide and neap tide conditions. During the king tide conditions, the tidal forcing was the dominant water exchange and circulation mechanism in the estuary. In contrast, the long-term oscillations linked with internal and external resonance played a major role in the turbulent mixing during neap tides. The data set showed further that the upper estuarine zone was drastically less affected by the spring tide range: the flow motion remained slow, but the turbulent velocity data were affected by the propagation of a transient front during the very early flood tide motion at the sampling site.