Large-scale turbulent structure of uniform shallow free-surface flows

The hydrodynamics of super- and sub-critical shallow uniform free-surface flows are assessed using laboratory experiments aimed at identifying and quantifying flow structure at scales larger than the flow depth. In particular, we provide information on probability distributions of horizontal velocity components, their correlation functions, velocity spectra, and structure functions for the near-water-surface flow region. The data suggest that for the high Froude number flows the structure of the near-surface layer resembles that of two-dimensional turbulence with an inverse energy cascade. In contrast, although large-scale velocity fluctuations were also present in low Froude number flow its behaviour was different, with a direct energy cascade. Based on our results and some published data we suggest a physical explanation for the observed behaviours. The experiments support Jirka’s [Jirka GH (2001) J Hydraul Res 39(6):567–573] hypothesis that secondary instabilities of the base flow may generate large-scale two-dimensional eddies, even in the absence of transverse gradients in the time-averaged flow properties.     

Fluid dynamics at the margin of rotational control

Geophysicists have developed a perspective on the dynamical influence of Earth’s rotation, while most other areas of fluid dynamics can safely disregard rotation. In this essay the dominant turbulence and wave behaviors in the rotating and non-rotating fluid-dynamical realms are described, and particular attention is given to their borderlands, where rotational influences are significant but not dominant. Contrary to the inverse energy cascade of geostrophic turbulence toward larger scales, a forward energy cascade toward microscale dissipation develops within the borderlands from the breakdown of diagnostic force balances, frontogenesis, and frontal instabilities, and then it continues further through the small-scale non-rotating realm until it dissipates at the microscale.     

Experimental investigation of meandering jets in shallow reservoirs

Meandering flows in rectangular shallow reservoirs were experimentally investigated. The characteristic frequency, the longitudinal wave length and the mean lateral extension of the meandering jet were extracted from the first paired modes, obtained by a proper orthogonal decomposition of the surface velocity field measured by large scale PIV. The depth-normalised characteristic lengths and the Strouhal number were then compared to the main dimensionless numbers characterizing the experiments: Froude number, friction number and reservoir shape factor. The normalised wave length and mean lateral extension of the meandering jet are neither correlated with the Froude number nor with the reservoir shape factor; but a clear relationship is found with the friction number. Similarly, the Strouhal number is found proportional to a negative power of the friction number. In contrast, the Froude number and the reservoir shape factor enable to predict the occurrence of a meandering flow pattern: meandering jets occur for Froude number greater than 0.21 and for a shape factor smaller than 6.2.     

Far-field vortex structures in mountain wakes

Laboratory experiments on scaled mountain wakes have provided the first demonstration that indicates the existence of two distinct mechanisms for the formation of the long-lived vortex structures that have been frequently observed in satellite imagery. In a series of laboratory experiments the well-known von Kármán vortex shedding mechanism has been observed at low values of Froude number (corresponding to low speed flow and strong stratification). At high Froude numbers the initial highly turbulent wake evolves into a vortex pattern similar to the far-wake dipole eddies previously observed in laboratory studies of submerged stratified wakes. In satellite imagery vortex structures have been observed to persist for hundreds of kilometers. In the present laboratory experiments far-wake dipole eddies were observed to persist for at least 200 diameters downstream, which corresponds to ~4,000?km when scaled to actual mountains. When generated at multiple sites along a continental mountain chain, these eddy structures could play a significant role in the development and variability of multi-scale weather and climate patterns.     

Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow

In this study, a three-dimensional model was used to numerically study the buoyant flow, along with its mixing characteristics, of heated water discharged from the surface and submerged side outfalls in shallow and deep water with a cross flow. Hydraulic experimental data were used to evaluate the applicability of the model. The simulation results agree well with the experimental results, particularly, the jet trajectories, the dimensions of the recirculating zone, and the distribution of the dimensionless excess temperature. The level of accuracy of the simulation results of the present study is nearly identical to that of the results conducted by McGuirk and Rodi (1978). If the heated water is discharged into shallow water where the momentum flux ratio and the discharge densimetric Froude number are high, the submerged discharge method is better than the surface discharge method in terms of the scale of the recirculating zone and the minimum dilution. In deep water, where the momentum flux ratio and discharge densimetric Froude number are low, however, the submerged discharge method had few advantages. In shallow water, the discharge jet is deflected by the ambient cross flow, while forcing the ambient flow to bend towards the far bank for the full depth. For a submerged discharge in shallow water, the recirculating zone is the largest in the lowest layer but becomes smaller in the upper layer. As the water depth increases, the ambient flow goes over the jet and diminishes the blocking effect, thereby decreasing the bending of the jet.     

Environmental Impact of Undular Tidal Bores in Tropical Rivers

A tidal bore impacts significantly on the estuarine ecosystem, although little is known on the flow field, mixing and sediment motion beneath tidal bores. In the absence of detailed systematic field measurements, a quasi-steady flow analogy was applied to investigate undular tidal bores with inflow Froude numbers between 1.25 and 1.6. Experimental results indicated that rapid flow redistributions occur beneath the free-surface undulations, with significant variations in bed shear stress between wave crests and troughs. Dynamic similarity was used to predict detailed flow characteristics of undular tidal bores. The effects of periodic loading on river sediments, scour of river bed and flow mixing behind the bore are discussed. A better understanding of these processes will contribute to better management practices in tidal bore affected rivers, including the Styx and Daly rivers in tropical Australia.     

Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.     

The generation of internal waves by tidal flow over continental shelf/slope topography

Experiments were performed to investigate the generation of internal waves induced by a barotropic tidal flow over continental shelf/slope topography in a two-layer stratified fluid. The interaction of the barotropic tide with the continental slope resulted in the formation of both linear and nonlinear waveforms, ranging from a linear wave of depression to a highly nonlinear internal bolus. The type of wave response was strongly dependent on the tidal forcing intensity and the position of the density interface relative to the shelf depth. Based on the values of the internal Froude number and the layer depth ratio, we delineate four distinct generation regimes, each with a distinct wave response.     

Calculation of Entrainment in Density Jumps

The density jump in a two-layer channel flow of miscible fluids, in which one of the layers is infinitely deep and at rest, is analyzed using the momentum flux and mass flux conservation equations. The analysis yields simple equations relating the heights upstream and downstream of the jump with the upstream Froude number and the rate of entrainment into the moving layer, as well as a relation between the maximum possible entrainment and the upstream Froude number. The analysis also shows that when the flow down-stream of the jump is controlled by an obstruction or a contraction, the entrainment rate into the jump depends solely on the dimensionless obstruction height contraction ratio and the upstream Froude number.     

Gravity currents in rotating,wedge-shaped,adverse channels

Results are presented from a series of parametric experimental and analytical studies of the behaviour of dense gravity currents along rotating, up-sloping, wedge-shaped channels. High resolution density profile measurements at fixed cross- and along-channel locations reveal the outflowing bottom gravity currents to adjust to quasi-steady, geostrophically-balanced conditions along the channels, with the outflow layer thickness and cross-channel interface slope shown to scale with the inlet Burger number for all experimental conditions tested. A general analytical solution to the classic rotating hydraulics problem has been developed under the assumption of inviscid, zero-potential-vorticity conditions to model dense water flow through a triangular constriction and thus simulate the vee-channel configurations under consideration. Predictions from this zero-PV model are shown to provide good overall quantitative agreement with experimental measurements obtained both under hydraulically-controlled conditions at the channel exit and for subcritical conditions generated along the channel length. Quantitative discrepancies between measurements and analytical predictions are attributed primarily to assumptions and limitations associated with the zero-PV modelling approach adopted, as well as the to the rapid adjustment in outflow characteristics as the channel exit is approached, as characterised by the along-channel variation in densimetric Froude number for the outflows.     

Froude scaling limitations in modeling of turbidity currents

The scaling problem associated with the modeling of turbidity currents has been recognized but is yet to be explored systematically. This paper presents an analysis of the dimensionless governing equations of turbidity currents to investigate the scale effect. Three types of flow conditions are considered: (i) conservative density current; (ii) purely depositional turbidity current; and (iii) mixed erosional/depositional turbidity current. Two controlling dimensionless numbers, the Froude number and the Reynolds number, appear in the non-dimensional governing equations. When densimetric Froude similarity is satisfied, the analysis shows that the results would be scale-invariant for conservative density current under the rough turbulent condition. In the case of purely depositional flows, truly scale-invariant results cannot be obtained, as the Reynolds-mediated scale effects appear in the bottom boundary conditions of the flow velocity and sediment fall velocity. However, the scale effect would be relatively modest. The Reynolds effect becomes more significant for erosional or mixed erosional/depositional turbidity currents as Reynolds-mediated scale effects also appear in the sediment entrainment relation. Numerical simulations have been conducted at three different scales by considering densimetric Froude scaling alone as well as combined densimetric Froude and Reynolds similarity. Simulation results confirm that although the scaling of densimetric Froude number alone can produce scale-invariable results for conservative density currents, variations occur in the case of turbidity currents. The results become scale invariant when densimetric Froude and Reynolds similarities are satisfied simultaneously.     

Stratified flow interactions with a suspended canopy

Field observations of the interactions between a stratified flow and a canopy suspended from the free surface above a solid boundary are described and analysed. Data were recorded in and around the canopy formed by a large long-line mussel farm. The canopy causes a partial blockage of the water flow, reducing velocities in the upper water column. Deceleration of the approaching flow results in a deepening of isopycnals upstream of the canopy. Energy considerations show that the potential for an approaching stratified flow to be diverted beneath a porous canopy is indicated by a densimetric Froude number. Strong stratification or low-velocities inhibit vertical diversion beneath the canopy, instead favouring a horizontal diversion around the sides. The effect on vertical mixing is also considered with a shear layer generated beneath the canopy and turbulence generated from drag within the canopy. In the observations, stratification is shown to be of sufficient strength to limit the effectiveness of the first mixing process, while the turbulence within the canopy is likely to enhance vertical exchange. Velocity and temperature microstructure measurements are used to investigate the effect of the canopy on turbulent dissipation and show that dissipation is enhanced within the canopy.     

Flow regimes in gaps within stands of flexible vegetation: laboratory flume simulations

Analyses of results from laboratory flume experiments are presented in which flow within gaps in canopies of flexible, submerged aquatic vegetation simulations is investigated. The aims of the work are (a) to identify the different flow regimes that may be found within such gaps, using Morris’ classical definitions of skimming flow, wake interference flow and isolated roughness flow as a template, (b) to determine the parameter space in which those flow regimes are most consistently delineated, and (c) to provide quantitative measurements of the loci of each flow regime within that parameter space for these experiments. The sedimentary and biological implications of each flow regime are also discussed. The results show that five flow regimes may be identified, expanding on Morris’ original set of three. The five are: (i) skimming flow; (ii) recirculation flow; (iii) boundary layer recovery; (iv) canopy through-flow; and (v) isolated roughness flow, the last being assumed to occur in some cases though it is not directly observed in these experiments. A Reynolds number based on the canopy overflow speed and the gap depth, and the gap aspect ratio are found to be the key parameters that determine these flow regimes, though a Froude number is found to be important for determining bed shear stress, and the length of leaves overhanging the gap from the upstream canopy is found to be important in determining the location of flow recirculation cells within the gap.     

Hydraulic modeling of mussel habitat at a bridge-replacement site, Allegheny River, Pennsylvania, USA

The Allegheny River in Pennsylvania supports a large and diverse freshwater-mussel community, including two federally listed endangered species, Pleurobema clava (Clubshell) and Epioblasma torulosa rangiana (Northern Riffleshell). It is recognized that river hydraulics and morphology play important roles in mussel distribution. To assess the hydraulic influences of bridge replacement on mussel habitat, metrics such as depth, velocity, and their derivatives (shear stress, Froude number) were collected or computed.The objectives of the project were to evaluate mussel and hydraulic data at a reference site and to compare those findings to a bridge-replacement site. The findings were used to support a statistical analysis, which establishes correlations between mussel count and hydraulics, and a numerical model to forecast habitat based on the statistics.ArcGIS was selected to manage the data and generate a grid to compute area statistics for 3319, 4.9-m × 4.9-m cells (cell) for total mussel count, depth, velocity, shear stress, and Froude number. The Wilcoxon Rank Sum test indicated no statistical significance between the total mussel count and the hydraulic variables; however, trellis graphs were used to account for the spatial variability in the data set. For the flow conditions measured, the total mussel count per cell is greatest at sections where (1) velocities range from 0.061 to 0.21 m/s, (2) shear stresses range from 0.48 to 3.8 dyne/cm², and (3) Froude numbers range from 0.006 to 0.04.Based on the statistical targets established, the hydraulic model results suggest that an additional 2428 m² or a 30-percent increase in suitable mussel habitat could be generated at the replacement-bridge site when compared to the baseline condition associated with the existing bridge at that same location. The study did not address the influences of substrate, acid mine drainage, sediment loads from tributaries, and surface-water/ground-water exchange on mussel habitat. Future studies could include methods for quantifying (1) channel-substrate composition and distribution using tools such as hydroacoustic echosounders specifically designed and calibrated to identify bed composition and mussel populations, (2) surface-water and ground-water interactions, and (3) a high-streamflow event.     

Large-Eddy Simulation of Coastal Upwelling Flow

Large-eddy simulations were carried out to study laboratory-scale realizations of coastal upwelling in an annular rotating tank with a sloping bottom. A two-layer stratified fluid was set into rigid body motion with the tank and then driven by the relative rotation of a solid top lid. The simulation code developed in this work was a three-dimensional incompressible Navier-Stokes solver using the message passing interface. The simulation runs were performed on a distributed memory massively parallel computer, namely, the IBM SP2. The simulation results were able to reveal the evolution of the complex upwelling structures in detail. The results were used to compare with and to complement two relevant series of coastal upwelling experiments. A Rayleigh-Taylor type of instability took place in the top inversion layer due to the unstable stratification after establishment of the upwelling front. The primary upwelling front was unstable to azimuthal perturbations and developed large amplitude baroclinic waves. The frontal wave structure consists of cyclone/anticyclone pairs. Whether cyclonic eddies containing the lower-layer fluid pinch off from the front depends on the _* value. The non-dimensional parameter _*=gh ₀/u _* f_s, which was first introduced by Narimousa and Maxworthy, combines the effects of stratification, rotation and surface stress and can be used to characterize the upwelling flow field. Our studies show that the frontal instabilities are much more intense and the upwelling front itself displays strong unsteadiness and cyclonic eddies containing the lower-layer fluid pinch off from the front when _* is significantly less than 5.8. For _*=5.8, the frontal instabilities are less intense and no pinched-off process is observed. To separate these regimes, a critical value of _* of about 5.4 is consistent with Narimousa and Maxworthy's results.     

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.     

Impulsive waves caused by subaerial landslides

This paper presents the experimental results of impulsive waves caused by subaerial landslides. A wide range of effective parameters are considered and studied by performing 120 laboratory tests. Considered slide masses are both rigid and deformable. The effects of bed slope angle, water depth, slide impact velocity, geometry, shape and deformation on impulse wave characteristics have been inspected. The impulse wave features such as amplitude, period and also energy conversation are studied. The effects of slide Froude number and deformation on energy conversation from slide into wave are also investigated. Based on laboratory measured data an empirical equation for impulse wave amplitude and period have been presented and successfully verified using available data of previous laboratory works.     

Fourier analysis of the roll-up and merging of coherent structures in shallow mixing layers

The current study investigates the role of nonlinearity in the development of two-dimensional coherent structures (2DCS) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate temporal shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layer. The flow is periodic in the stream-wise direction and the transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. The numerical results are examined with the aid of Fourier decomposition. Results show that the previous success in applying local linear theory to shallow mixing layers does not imply that the flow is truly linear. Linear stability theory is confirmed to be only valid within a short distance from the inflow boundary. Downstream of this linear region, nonlinearity becomes important for the roll-up and merging of 2DCS. While the energy required for the merging of 2DCS is still largely provided by the velocity shear, the merging mechanism is one where nonlinear mode interaction changes the velocity field of the subharmonic mode and the gradient of the along-stream velocity profile which, in turn, changes the magnitude of the energy production of the subharmonic mode by the velocity shear implicitly. The nonlinear mode interaction is associated with energy up-scaling and is consistent with the inverse energy cascade which is expected to occur in shallow shear flows. Current results also show that such implicit nonlinear interaction is sensitive to the phase angle difference between the most unstable mode and its subharmonic. The bed friction effect on the 2DCS is relatively small initially and grows in tandem with the size of the 2DCS. The bed friction also causes a decrease in the velocity gradient as the flow develops downstream. The transition from unstable to stable flow occurs when the bed friction balances the energy production. Beyond this point, the bed friction is more dominant and the 2DCS are progressively damped and eventually get annihilated. The energy production by the velocity shear plays an important role from the upstream end all the way to the point of transition to stable flow. The fact that linear stability theory is valid only for a short distance from the inflow boundary suggests that some elements of nonlinearity is incorporated in the mean velocity profile in experiments by the averaging process. The implicit nature of nonlinear interaction in shallow mixing layers and the sensitivity of the nonlinear interaction to phase angle difference between the most unstable mode and its subharmonic allows local linear theory to be successful in reproducing features of the instability such as the dominant mode of the 2DCS and its amplitude.     

Free surface profiles of near-critical instabilities in open channel flows: undular hydraulic jumps

Environmental Fluid Mechanics - When the Froude number F of a free-surface flow ranges between 0.3 and 3, the flow is unstable and frequently characterised by free surface undulations, with the...