Self-aeration in the rapidly- and gradually-varying flow regions of steep smooth and stepped spillways

In high-velocity chute flows, free-surface aeration is often observed. The phenomenon is called self-aeration or white waters. When the turbulent shear stresses next to the free-surface are large enough, air bubbles are entrained throughout the entire air–water column. A rapidly-varied flow region is observed immediately downstream of the inception point of free-surface aeration. An analytical solution of the air diffusion equation is proposed and the results compare well with new experimental data. Both experiments and theory indicate that the flow bulking spans over approximately 3–4 step cavities downstream of the inception point of free-surface aeration on a stepped chute. Further downstream the void fraction distributions follow closely earlier solutions of the air diffusion equation. The application of the diffusion equation solution to prototype and laboratory data shows air bubble diffusivities typically larger than the momentum transfer coefficient. The result highlights however a marked decrease in the ratio of air bubble diffusivity to eddy viscosity with increasing Reynolds number. The finding might indicate some limitation of laboratory investigations of air bubble diffusion process in self-aerated flows and of their extrapolation to full-scale prototype applications.     

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.     

Advanced post-processing and correlation analyses in high-velocity air–water flows

The interest in air–water flows has not diminished in recent years, but it is accompanied by frequent citations of early, sometimes outdated articles. A basic issue is the inadequate, incomplete interpretation of air–water flow instrumentation by hydraulic engineers and researchers. This article comments on high-velocity air–water flow measurements by means of intrusive phase detection probes. This article focus on the bubbly flow structure of high-velocity air–water flow based upon measurements by means of intrusive phase detection probes. It is shown that some advanced post-processing techniques may yield expanded information on the air–water turbulent flow properties and bubbly flow structures. The outcomes demonstrate simple techniques in high-velocity air–water flow analysis.     

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.     

Surface waves and roughness in self-aerated supercritical flow

In high-velocity open channel flows, free-surface aeration is commonly observed. The effects of surface waves on the air–water flow properties are tested herein. The study simulates the air–water flow past a fixed-location phase-detection probe by introducing random fluctuations of the flow depth. The present model yields results that are close to experimental observations in terms of void fraction, bubble count rate and bubble/droplet chord size distributions. The results show that the surface waves have relatively little impact on the void fraction profiles, but that the bubble count rate profiles and the distributions of bubble and chord sizes are affected by the presence of surface waves.     

Energy dissipation, flow resistance and gas-liquid interfacial area in skimming flows on moderate-slope stepped spillways

With the re-evaluation and revision of a number of design floods, several embankment overtopping protection systems have been developed and a common technique is the construction of a stepped spillway on the downstream slope. For such moderate slope stepped channels, detailed air–water flow measurements were performed in a large facility with a focus on the rate of energy dissipation, flow resistance, air–water interfacial areas and re-aeration rates. Past and present experimental results showed a significant aeration of the flow. The median dimensionless residual head was about 3 × d_c for the 21.8° sloping chute and smaller than that for flatter slopes (θ = 3.4° and 15.9°). The flow resistance results yielded an equivalent Darcy friction factor of about 0.25 implying a larger flow resistance for the 21.8° slope angle than for smaller slope angles. The re-aeration rate was deduced from the integration of the mass transfer equation using measured air–water interfacial areas and air–water flow velocities. The results suggested an increasing re-aeration rate with increasing rate of energy dissipation. The stepped invert contributed to intense turbulence production, free-surface aeration and large interfacial areas. The experimental data showed however some distinctive seesaw pattern in the longitudinal distribution of air–water flow properties with a wave length of about two step cavities. While these may be caused by the interactions between successive adjacent step cavities and their interference with the free-surface, the existence of such “instabilities” implies that the traditional concept of normal flow might not exist in skimming flows above moderate-slope stepped spillways.     

Bayesian model-based cluster analysis for predicting macrofaunal communities

To predict macrofaunal community composition from environmental data a two-step approach is often followed: (1) the water samples are clustered into groups on the basis of the macrofauna data and (2) the groups are related to the environmental data, e.g. by discriminant analysis. For the cluster analysis in step 1 many hard, seemingly arbitrary choices have to be made that nevertheless influence the solution (similarity measure, clustering strategy, number of clusters). The stability of the solution is often of concern, e.g. in clustering by the program. In the discriminant analysis of step 2 it can occur that a water sample is misclassified on the basis of the environmental data but on further inspection happens to be a borderline case in the cluster analysis. One would then rather reclassify such a sample and iterate the two steps. Bayesian latent class analysis is a flexible, extendable model-based cluster analysis approach that recently has gained popularity in the statistical literature and that has the potential to address these problems. It allows the macrofauna and environmental data to be modelled and analyzed in a single integrated analysis. An exciting extension is to incorporate in the analysis prior information on the habitat preferences of the macrofauna taxa such as is available in lists of indicator values. The output of the analysis is not a hard assignment of water samples to clusters but a probabilistic (fuzzy) assignment. The number of clusters is determined on the basis of the Bayes factor. A standard feature of the Bayesian method is to make predictions and to assess their uncertainty. We applied this approach to a data set consisting of 70 water samples, 484 macrofauna taxa and four environmental variables for which previously a five cluster solution had been proposed. The standard for Bayesian estimation, the Gibbs sampler, worked fine on a subset with only 12 selected taxa but did not converge on the full set with 484 taxa. This is due to many configurations in which the assignment probabilities are all very close to either 0 or 1. This convergence problem is comparable with the local optima problem in classical cluster optimization algorithms, including the EM algorithm used in Latent Gold, a Windows program for latent class analysis. The convergence problem needs to be solved before the benefits of Bayesian latent class analysis can come to fruition in this application. We discuss possible solutions.     

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.     

Factors influencing clustering in the intertidal snail Cerithium moniliferum

The snail Cerithium moniliferum Kiener at Heron Island (southern end of the Great Barrier Reef) clusters at outgoing tide and disperses at incoming tide. In containers with a constant water level, but exposed to the normal day and night rhythm, the snails cluster and disperse rhythmically for up to several days. Graphical simulation and spectral analysis, preceded by first difference filtering, showed two components, one corresponding to the tidal, and the other to the diurnal rhythm. Snails kept under continuous light also exhibit a two-component rhythm of clustering and dispersal. In containers with high (40 cm) water level and under continuous light or normal day and night rhythm, the two-component rhythm may—for some unknown reason—be replaced by long fluctuations of dispersal, with clustering occurring in the afternoon only in snails kept under normal alternation of day and night. Snails kept for 24 h in dryness prior to putting them into containers with water, remain clustered (sometimes with an initial dispersal “shock”) for a long time, with subsequent slow fluctuations of clustering and dispersal, probably induced by the light conditions. Factors contributing to clustering and dispersal are: (1) increased rate of dropping from container walls at outgoing tide; (2) lack of locomotory activity at low tide; (3) negative geotaxis at incoming tide; (4) positive geotaxis at outgoing tide; (5) attraction of snails by other snails (positive chemotaxis) throughout the rhythm. Sudden exposure to light increases the rate of dropping from the container walls and may be a contributing factor to stronger clustering during the day. Factors 1–4 show endogenous periodicities.     

Ecology of subtropical hermit crabs in S.W. Madagascar: cluster structure and function

The Madagascar coast (both) has a higher density and diversity of hermit crabs than is known from any other locality in the western Indian Ocean. Of the 20 species occurring at Anakao (S.W. Madagascar), 11 aggregated into clusters, including all but one of the species above the subtidal zone. The mean number of hermit crabs and species in clusters varied with several spatial parameters and time. Over 80% of the community clustered in certain habitats in particular tidal zones, whilst as low as 3% clustered in others. The highest intensity of clustering with shore zone was coincident with peak numbers of hermit crabs. The initiation and duration of clusters of hermit crabs above the eulittoral was driven by circadian rhythms, whilst those in the eulittoral were governed by tidal state. Clusters above and below the eulittoral were longer in duration, and those in the subtidal were more temporally variable than those above it. Certain pairs of species showed positive correlations of occurrence and (more rarely) of abundance, and all the correlations (of occurrence) of one, Clibanarius eurysternus, were negative. Positive correlations of occurrence were related to the degree of shell-use commonality between species pairs. Eulittoral species clustered with other individuals of approximately similar size and exchanged shells upon cluster disintegration. There was evidence of a dehydration-reduction function to shell clustering in addition to shell-exchange facilitation. The highly variable species-specific strategies of clustering may be important in alleviating both intra- and interspecific competition in assemblages of similar and highly abundant species.     

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.     

Testing spatial cluster occurrence in maps equipped with environmentally defined structures

We propose a novel tool for testing hypotheses concerning the adequacy of environmentally defined factors for local clustering of diseases, through the comparative evaluation of the significance of the most likely clusters detected under maps whose neighborhood structures were modified according to those factors. A multi-objective genetic algorithm scan statistic is employed for finding spatial clusters in a map divided in a finite number of regions, whose adjacency is defined by a graph structure. This cluster finder maximizes two objectives, the spatial scan statistic and the regularity of cluster shape. Instead of specifying locations for the possible clusters a priori, as is currently done for cluster finders based on focused algorithms, we alter the usual adjacency induced by the common geographical boundary between regions. In our approach, the connectivity between regions is reinforced or weakened, according to certain environmental features of interest associated with the map. We build various plausible scenarios, each time modifying the adjacency structure on specific geographic areas in the map, and run the multi-objective genetic algorithm for selecting the best cluster solutions for each one of the selected scenarios. The statistical significances of the most likely clusters are estimated through Monte Carlo simulations. The clusters with the lowest estimated p-values, along with their corresponding maps of enhanced environmental features, are displayed for comparative analysis. Therefore the probability of cluster detection is increased or decreased, according to changes made in the adjacency graph structure, related to the selection of environmental features. The eventual identification of the specific environmental conditions which induce the most significant clusters enables the practitioner to accept or reject different hypotheses concerning the relevance of geographical factors. Numerical simulation studies and an application for malaria clusters in Brazil are presented.     

Turbulent air–water flows in hydraulic structures: dynamic similarity and scale effects

In hydraulic structures, free-surface aeration is commonly observed: i.e., the white waters. The air bubble entrainment may be localised (hydraulic jumps, plunging jets) or continuous along an interface (water jets, chutes). Despite recent advances, there are some basic concerns about the extrapolation of laboratory results to large size prototype structures. Herein the basic air bubble entrainment processes are reviewed and the relevant dynamic similarities are discussed. Traditionally, physical studies are conducted using a Froude similitude which implies drastically smaller laboratory Reynolds numbers than in the corresponding prototype flows. Basic dimensional analyses are developed for both singular and interfacial aeration processes. The results are discussed in the light of systematic investigations and they show that the notion of scale effects is closely linked with the selection of relevant characteristic air–water flow properties. Recent studies of local air–water flow properties highlight that turbulence levels, entrained bubble sizes and interfacial areas are improperly scaled based upon a Froude similitude even in large-size models operating with the so defined Reynolds numbers ρ _w × q _w/μ _w up to 5 E+5. In laboratory models, the dimensionless turbulence levels, air–water interfacial areas and mass transfer rates are drastically underestimated.     

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 abundance from presence–absence maps using cluster models

A presence–absence map consists of indicators of the occurrence or nonoccurrence of a given species in each cell over a grid, without counting the number of individuals in a cell once it is known it is occupied. They are commonly used to estimate the distribution of a species, but our interest is in using these data to estimate the abundance of the species. In practice, certain types of species (in particular flora types) may be spatially clustered. For example, some plant communities will naturally group together according to similar environmental characteristics within a given area. To estimate abundance, we develop an approach based on clustered negative binomial models with unknown cluster sizes. Our approach uses working clusters of cells to construct an estimator which we show is consistent. We also introduce a new concept called super-clustering used to estimate components of the standard errors and interval estimators. A simulation study is conducted to examine the performance of the estimators and they are applied to real data.     

Two-phase flow insights into open-channel flows with suspended particles of different densities

The effect of particle density on the turbulent open-channel flow carrying dilute particle suspensions is investigated using two specific gravities and three concentrations of solid particles. The particles, identical in size and similar in shape, were natural sand and a neutrally buoyant plastic. The particles were fully suspended, and formed no particle streaks on the channel’s bed. Accordingly, the changes in the flow are attributed to the interactions between suspended particles and flow turbulence structures. Measurements were obtained by means of image velocimetry enabling simultaneous, but distinct, measurement of liquid and particle velocities. The experimental results show that, irrespective of particle specific gravity, particle suspension influences bulk velocity of flow and the Kármán coefficient, while friction velocity essentially remains constant. The results also show that particles in suspension modify local water turbulence over the flow depth, but in ways not accurately predicted using the customary parameters for characterizing turbulence modification.     

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.     

Drivers of carbon dioxide emissions: an empirical investigation using hierarchical and non-hierarchical clustering methods

The mitigation of CO2 emissions requires a global effort with common but differentiated responsibilities. In this paper, we identify clusters of CO2 emissions across 72 countries. First, using the stochastic version of the IPAT and employing the dynamic common correlated effects technique, we identify three key determinants affecting CO2 emissions (non-renewables, population, and real GDP). In the second step, both hierarchical and non-hierarchical clustering methods are considered to identify the optimal number of clusters. We identify two to four clusters with different member countries, and in particular establish that in most cases, a 2-cluster solution appears to be optimal. The contents of clusters vary slightly according to the clustering methods for each period. The clustering results from using only the overall CO2 emissions indicate that the countries we consider form three clusters, with China and the USA each within a single member cluster. The remaining 70 countries form the third cluster. Our findings reflect the prominent roles of China and the USA in overall CO2 emissions. Analyses with sub-period and largest emitters reflect a different clustering structure. Some policy recommendations in setting emission reductions are made, considering different clusters across countries.     

Adaptive cluster sampling with clusters selected without replacement and stopping rule

Adaptive cluster sampling (ACS) has received much attention in recent years since it yields more precise estimates than conventional sampling designs when applied to rare and clustered populations. These results, however, are impacted by the availability of some prior knowledge about the spatial distribution and the absolute abundance of the population under study. This prior information helps the researcher to select a suitable critical value that triggers the adaptive search, the neighborhood definition and the initial sample size. A bad setting of the ACS design would worsen the performance of the adaptive estimators. In particular, one of the greatest weaknesses in ACS is the inability to control the final sampling effort if, for example, the critical value is set too low. To overcome this drawback one can introduce ACS with clusters selected without replacement where one can fix in advance the number of distinct clusters to be selected or ACS with a stopping rule which stops the adaptive sampling when a predetermined sample size limit is reached or when a given stopping rule is verified. However, the stopping rule breaks down the theoretical basis for the unbiasedness of the ACS estimators introducing an unknown amount of bias in the estimates. The current study improves the performance of ACS when applied to patchy and clustered but not rare populations and/or less clustered populations. This is done by combining the stopping rule with ACS without replacement of clusters so as to further limit the sampling effort in form of traveling expenses by avoiding repeat observations and by reducing the final sample size. The performance of the proposed design is investigated using simulated and real data.     

Drinking water biotic safety of particles and bacteria attached to fines in activated carbon process

In this paper, the drinking water biotic safety of particles and bacteria attached to fines in activated carbon process was investigated by actual treatment process and advanced treatment pilot trial with granular activated carbon. In the experiment, the particles were detected by IBR particle calculating instrument, the activated carbon fines were counted on the basis of the most probable number (MPN) with a microscope, the total number of bacteria was analyzed between the conventional agar culture medium and the one with R2A, and the bacteria attached to activated carbon fines was resolved by the homogenization technique. The experimental results showed that the average total number of particles was 205 CNT/mL in the activated carbon effluent during a filter cycle, of which the number of particles with sizes > 2 μm was 77 CNT/mL more than the present particle control criterion of the American drinking water product standard (50 CNT/mL). The backwash of low density and long duration lowered particle number in the effluent. The MPN of activated carbon fines in the effluent was between 400 and 600 CNT/L, which accounted for less than 5‰ of the total particles from activated carbon filtration for a poor relative level (R ² = 0.34). The microorganisms in activated carbon effluent consisted mostly of heterotrophic bacillus and the total bacteria number was five times as high as that of the inflow, i.e. the effluent from sand filter. The actual bacteria number may be truly indicated by the detection technique with R2A culture medium compared with the traditional agar cultivation. The inactivation efficiency of bacteria attached to activated carbon fines was less than 40% under 1.1 mg/L of chlorine contacting for 40 min. Results showed that the particles and bacteria attached to activated carbon fines may influence drinking water biotic safety, and that the effective control measures need to be further investigated.