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.     

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.     

Monitoring and estimating the flow conditions and fish presence probability under various flow conditions at reach scale using genetic algorithms and kriging methods

The combination of current velocity and water depth influences stream flow conditions, and fish activities prefer particular flow conditions. This study develops a novel optimal flow classification method for identifying types of stream flow based on the current velocity and the water depth using a genetic algorithm. It is applied to the Datuan stream in northern Taiwan. Fish were sampled and their habitat investigated at the study site during the spring, summer, fall and winter of 2008-2009. The current velocity, water depth and maps of the presence probability of fish were estimated by ordinary and indicator kriging. The optimal classification results were compared with the classification results obtained using the Froude number and empirical methods. The flow classification results demonstrate that the proposed optimal flow classification method that considers depth-velocity and optimally identified criteria for classifying flow types, yields a current velocity and water depth of 0.32 (m/s) and 0.29 (m), respectively, and classifies the flow conditions in the study area as pool, run, riffle and slack. The variography results of the current velocity and the water depth data reveal that seasonal flows are not spatially stationary among seasons in the study area. Kriging methods and a two-dimensional hydrodynamic model (River 2D) with empirical and optimal flow classification methods are more effective than the Froude number method in classifying flow conditions in the study area. The flow condition classifications and probability maps were generated by River 2D, ordinary kriging and indicator kriging, to quantify the flow conditions preferred by Sicyopterus japonicus in the study area. However, the proposed optimal classification method with kriging and River 2D is an effective alternative method for mapping flow conditions and determining the relationship between flow and the presence probability of target fish in support of stream restoration.     

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.     

Mass-based depth and velocity scales for gravity currents and related flows

Gravity driven flows on inclines can be caused by cold, saline or turbid inflows into water bodies. Another example are cold downslope winds, which are caused by cooling of the atmosphere at the lower boundary. In a well-known contribution, Ellison and Turner (ET) investigated such flows by making use of earlier work on free shear flows by Morton, Taylor and Turner (MTT). Their entrainment relation is compared here with a spread relation based on a diffusion model for jets by Prandtl. This diffusion approach is suitable for forced plumes on an incline, but only when the channel topography is uniform, and the flow remains supercritical. A second aspect considered here is that the structure of ET’s entrainment relation, and their shallow water equations, agrees with the one for open channel flows, but their depth and velocity scales are those for free shear flows, and derived from the velocity field. Conversely, the depth of an open channel flow is the vertical extent of the excess mass of the liquid phase, and the average velocity is the (known) discharge divided by the depth. As an alternative to ET’s parameterization, two sets of flow scales similar to those of open channel flows are outlined for gravity currents in unstratified environments. The common feature of the two sets is that the velocity scale is derived by dividing the buoyancy flux by the excess pressure at the bottom. The difference between them is the way the volume flux is accounted for, which—unlike in open channel flows—generally increases in the streamwise direction. The relations between the three sets of scales are established here for gravity currents by allowing for a constant co-flow in the upper layer. The actual ratios of the three width, velocity, and buoyancy scales are evaluated from available experimental data on gravity currents, and from field data on katabatic winds. A corresponding study for free shear flows is referred to. Finally, a comparison of mass-based scales with a number of other flow scales is carried out for available data on a two-layer flow over an obstacle. Mass-based flow scales can also be used for other types of flows, such as self-aerated flows on spillways, water jets in air, or bubble plumes.     

Turbulent velocity profile in fully-developed open channel flows

The determination of velocity profile in turbulent narrow open channels is a difficult task due to the significant effects of the anisotropic turbulence that involve the Prandtl’s second type of secondary flow occurring in the cross section. With these currents the maximum velocity appears below the free surface that is called dip phenomenon. The well-known logarithmic law describes the velocity distribution in the inner region of the turbulent boundary layer but it is not adapted to define the velocity profile in the outer region of narrow channels. This paper relies on an analysis of the Navier–Stokes equations and yields a new formulation of the vertical velocity profile in the center region of steady, fully developed turbulent flows in open channels. This formulation is able to predict time averaged primary velocity in the outer region of the turbulent boundary layer for both narrow and wide open channels. The proposed law is based on the knowledge of the aspect ratio and involves a parameter C_Ar depending on the position of the maximum velocity (ξ_dip). ξ_dip may be derived, either from measurements or from an empirical equation given in this paper. A wide range of longitudinal velocity profile data for narrow open channels has been used for validating the model. The agreement between the measured and the computed velocities is rather good, despite the simplification used.     

How do larvae attach to a solid in a laminar flow?

A hydrodynamic model explaining the mechanism of contact of marine larvae in vertical flows is presented. Two hydrodynamic factors—flow vorticity and larval self-propulsion—are the key components in the mathematical model. It is shown that flow vorticity causes a larva to rotate and change the direction of self-thrust, thus leading to its migration across the mean flow. The latter motion is of an oscillatory nature. Contact will be enabled only for sufficiently large amplitudes of oscillations. Simple expressions for the probability of initial contact are obtained for two-dimensional Couette and Poiseuille flows. The three-dimensional motion of a larva in a tube is studied using the Monte-Carlo simulations. It is shown that contact probability depends mainly on the ratio of the characteristic flow velocity and the larva’s swimming speed. The theoretical results compare favorably with available experimental data. Possible applications of the method and results presented here to the classical problem of larval attachment to bodies of general geometry are briefly discussed in the concluding section.     

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.     

Mechanism of the intermittent motion of two-phase debris flows

A typical two-phase debris flow exhibits a high and steep flow head consisting of rolling boulders and cobbles with intermittent or fluctuating velocity. The relative motion between the solid phase and the liquid phase is obvious. The motion of a two-phase debris flow depends not only on the rheological properties of the flow, but also on the energy transmission between the solid and liquid phases. Several models have been developed to study two-phase debris flows. An essential shortcoming of most of these models is the omission of the interaction between the two phases and identification of the different roles of the different materials in two-phase debris flows. The tracer particles were used for the velocity of solid phase and the velocity of liquid phase was calculated by the water velocity on the surface of the debris flow in the experiments. This paper analyzed the intermittent feature of two-phase debris flows based on videos of debris flows in the field and flume experiments. The experiments showed that the height of the head of the two-phase debris flow increased gradually in the initiation stage and reached equilibrium at a certain distance from the start of the debris flow. The height growth and the velocity of the flow head showed fluctuating characteristics. Model equations were established and the analyses proved that the average velocity of the two-phase debris flow head was proportional to the flood discharge and inversely proportional to the volume of the debris flow head.     

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.     

Flow and drag force around a free surface piercing cylinder for environmental applications

This paper investigates flows around a free surface piercing cylinder with Froude number F > 0.5 and Reynolds number around Re = 50,000. The aim of this work is to gain a better understanding of the flow behaviour in environmental systems such as fishways. The advances are based upon experimental and numerical results. Several flow discharges and slopes are tested to obtain both subcritical and supercritical flows. The drag force exerted on the cylinder is measured with the help of a torque gauge while the velocity field is obtained using particle velocimetry. For the numerical part, two URANS turbulence models are tested, the k-\(\omega\) SST and the RNG k-\(\varepsilon\) models using the OpenFOAM software suite for subcritical cases, and then compared with the corresponding experimental results. With fishways applications in mind, the changes in drag coefficient \(C_d\) versus Froude number and water depth are studied and experimental correlations proposed. We conclude that the most suitable URANS turbulence model for reproducing this kind of flow is the k-\(\omega\) SST model.     

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.     

Steady Shallow-Water Current and Solute Transport around a Semi-Conical Headland

Laboratory experiments have been carried out to investigate the effects of a sloping wall headland on the flow characteristics and the associated concentration distributions from a point source around the headland. A semi-conical headland with a slope of 1:2 was set up in a flow basin, 4.8 m long and 3.8 m wide. In this paper, the experimental results of a steady shallow-water current are reported. Three dimensional flow velocities in the basin were measured using Sontek-ADV instrument. The dye concentration levels in the basin were measured by two fluorometers. The experimental results showed a large-scale re-circulation region behind the semi-conical headland. The peak turbulence energy, at about 53% of the local kinetic flow energy, coincides with the region of high velocity gradient. Significant vertical flows were observed around the area near the downhill slope of the headland, with a maximum ratio of vertical to horizontal velocities being about 22%. Such relatively significant vertical scouring velocities, coupled with strong turbulence energy and high horizontal velocity gradients in the same region, could cause severe bed erosion. The experimental results have also been compared with the predicted results of a depth-averaged numerical model. The predicted eddy structure and the concentration distribution in the re-circulation area were found to compare favourably with the experimental results. However, the discrepancies in the flow velocities and the concentration levels near the headland were apparent. It was observed that the dye concentration continued to spread in the cross-stream direction after passing the headland, whereas only a limited extent of the lateral spreading was predicted by the numerical model further downstream of the headland.     

Instabilities in non-rotating and rotating shallow shear flows

Numerical simulations for the wave radiation effect on the linear and nonlinear instabilities of rotating and non-rotating shallow flows are conducted using shallow-water equations. At a low convective Froude number, the results of the instabilities is a string of eddies. The coalescence between the neighbouring eddies decides the transverse mixing of the shallow shear flow. At a higher convective Froude number, the development of the shear flow is characterized by wave radiation and the production of shocklets. The radiation of waves in the non-rotating shallow flow is a phenomenon analogous to the radiation of sound in gas dynamics. In the rotating flow on the other hand, the shallow-flow instabilities are intensified due to rotational interference within a window of instability over a narrow range of Rossby numbers.     

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.     

Energy-economic theory and mathematical models for combining the systems of man and nature, case study: The urban region of Miami, Florida

This paper presents a study of the urban region of Miami, Fla. with consideration of energy flow and the relationship between energy theory and economics. Much of this work is based on the theories and work of Dr Howard T. Odum at the University of Florida. A theory of energy quality is presented which is an attempt to relate energies of different concentrations in their ability to do work. This theory allows comparison of the systems of man and nature. Another theory is proposed which seeks to describe the ability of a region to compete based on its flows of natural and fossil fuel energies.Economic, natural system, and energy data were compiled for the Miami urban region from 1950–1972. Cross-correlation of this data showed significant levels of correlation between the rate of change of fossil fuel use and the rates of change of population, budget, sales tax, income, building structure, and number of telephones. Calculation of several urban indicators for 1972 showed a fossil fuel energy density of 300 kcal/m²/day in the urbanized area, a per capita energy consumption of 53.8 × 10⁶ kcal/capita/year, a ratio of natural to fossil fuel energies of 0.25, a developed area of 260 miles² (673.4 km²), and a rate of development of 6.5 miles² (16.8 km²) per year.An overall model of Miami is presented with flows and storages quantified for 1972. Based on this model a simpler model was simulated on an analog computer. This model consisted of a system of first-order in time, non-linear differential equations which included fossil fuel energy flows, main economic flows, external price functions, building structure, natural energies, and population. This model was simulated for several linearly increasing future price functions and several sets of future energy functions.Natural energies within the region were calculated by determining the land areas associated with various ecosystem types. Estimating the productivities of these systems on a per area basis allowed calculation of total energy flows. The energies associated with winds, tides, waves, and fresh/salt water concentration gradients were also determined. It was found that the ratio of natural to fossil fuel energy changed from 1.77 in 1950 to 0.25 in 1972.     

Dominant features in three-dimensional turbulence structure: comparison of non-uniform accelerating and decelerating flows

The results are presented from an experimental study to investigate three-dimensional turbulence structure profiles, including turbulence intensity and Reynolds stress, of different non-uniform open channel flows over smooth bed in subcritical flow regime. In the analysis, the uniform flow profiles have been used to compare with those of the non-uniform flows to investigate their time-averaged spatial flow turbulence structure characteristics. The measured non-uniform velocity profiles are used to verify the von Karman constant κ and to estimate sets of log-law integration constant B_r and wake parameter П, where their findings are also compared with values from previous studies. From κ, B_r and П findings, it has been found that the log-wake law can sufficiently represent the non-uniform flow in its non-modified form, and all κ, B_r and П follow universal rules for different bed roughness conditions. The non-uniform flow experiments also show that both the turbulence intensity and Reynolds stress are governed well by exponential pressure gradient parameter β equations. Their exponential constants are described by quadratic functions in the investigated β range. Through this experimental study, it has been observed that the decelerating flow shows higher empirical constants, in both the turbulence intensity and Reynolds stress compared to the accelerating flow. The decelerating flow also has stronger dominance to determine the flow non-uniformity, because it presents higher Reynolds stress profile than uniform flow, whereas the accelerating flow does not.     

Experimental analysis of Froude number effect on air entrainment in the hydraulic jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves, and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in a rectangular horizontal flume with partially developed inflow conditions. The vertical distributions of the void fraction and the air bubbles count rate were recorded for inflow Froude number Fr ₁ in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction C _max increases with the increasing Fr ₁. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was derived.