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.     

Laboratory measurements and multi-block numerical simulations of the mean flow and turbulence in the non-aerated skimming flow region of steep stepped spillways

We present and discuss the results of a comprehensive study addressing the non-aerated region of the skimming flow in steep stepped spillways. Although flows in stepped spillways are usually characterized by high air concentrations concomitant with high rates of energy dissipation, the non-aerated region becomes important in small dams and/or spillways with high specific discharges. A relatively large physical model of such spillway was used to acquire data on flow velocities and water levels and, then, well-resolved numerical simulations were performed with a commercial code to reproduce those experimental conditions. The numerical runs benefited from the ability of using multi-block grids in a Cartesian coordinate system, from capturing the free surface with the TruVOF method embedded in the code, and from the use of two turbulence models: the k-e{k{-}\varepsilon} and the RNGk-e{k{-}\varepsilon} models. Numerical results are in good agreement with the experimental data corresponding to three volumetric flow rates in terms of the time-averaged velocities measured at diverse steps in the spillway, and they are in very satisfactory agreement for water levels along the spillway. In addition, the numerical results provide information on the turbulence statistics of the flow. This work also discusses important aspects of the flow, such as the values of the exponents of the power-law velocity profiles, and the characteristics of the development of the boundary layer in the spillway.     

Numerical simulations of intrusive gravity currents interacting with a bottom-mounted obstacle in a continuously stratified ambient

In this study, the flow dynamics of intrusive gravity currents past a bottom-mounted obstacle were investigated using highly resolved numerical simulations. The propagation dynamics of a classic intrusive gravity current was first simulated in order to validate the numerical model with previous laboratory experiments. A bottom-mounted obstacle with a varying non-dimensional height of \(\tilde{D}=D/H\), where D is the obstacle height and H is the total flow depth, was then added to the problem in order to study the downstream flow pattern of the intrusive gravity current. For short obstacles, the intrusion re-established itself downstream without much distortion. However, for tall obstacles, the downstream flow was found to be a joint effect of horizontal advection, overshoot-springback phenomenon, and associated Kelvin-Helmholtz instabilities. Analysis of the numerical results show that the relationship between the downstream propagation speed and the obstacle height can be subdivided into three regimes: (1) a retarding regime (\(\tilde{D}\) \(\approx \) 0–0.3) where a 30 % increase in obstacle height leads to a 20 % reduction in propagation speed, simply due to the obstacle’s retarding effect; (2) an impounding regime (\(\tilde{D}\) \(\approx \) 0.3–0.6) where the additional 30 % increase in obstacle height only leads to a further (negligible) 5 % reduction in propagation speed, due to the accelerating effect of upstream impoundment and downstream enhanced mixing; and (3) a choking regime (\(\tilde{D}\) \(\approx \) 0.6–1.0) where the propagation speed is dramatically reduced due to the dominance of the obstacle’s blocking effect. The obstacle thickness was found to be irrelevant in determining the downstream propagation speed at least for the parameter range explored in this study. The present work highlights the significance of topographic effects in stratified flows with horizontal pressure forcing.     

Characterization of turbulence statistics on the non-aerated skimming flow over stepped spillways: a numerical study

Environmental Fluid Mechanics - In this work we address the mean flow and turbulence statistics in the non-aerated region of a stepped spillway by using two different numerical strategies in two...     

Flow aeration, cavity processes and energy dissipation on flat and pooled stepped spillways for embankments

The design floods of several reservoirs were recently re-evaluated and the revised spillway outflow could result in dam overtopping with catastrophic consequences for some embankment structures. Herein a physical study was performed on flat and pooled stepped spillways with a slope typical of embankments $(\uptheta = 26.6^{\circ })$ and four stepped configurations were tested: a stepped spillway with flat horizontal steps, a pooled stepped spillway, and two stepped spillways with in-line and staggered configurations of flat and pooled steps. The focus of the study was on the flow aeration, air–water flow properties, cavity flow processes, and energy dissipation performances. The results demonstrated the strong aeration of the flow for all configurations. On the in-line and staggered configurations of flat and pooled steps, the flow was highly three-dimensional. The residual head and energy dissipation rates at the stepped chute downstream end were calculated based upon the detailed air–water flow properties. The results showed that the residual energy was the lowest for the flat stepped weir. The data for the stepped spillway configuration with in-line and staggered configurations of flat and pooled steps showed large differences in terms of residual head in the transverse direction. Altogether the present results showed that, on a $26.6^{\circ }$ slope stepped chute, the designs with in-line and staggered configurations of flat and pooled steps did not provide any advantageous performances in terms of energy dissipation and flow aeration, but they were affected by three-dimensional patterns leading to some flow concentration.     

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.     

Characteristics of clustered particles in skimming flows on a stepped spillway

Air–water flows at hydraulic structures are commonly observed and called white waters. The free-surface aeration is characterised by some intense exchanges of air and water leading to complex air–water structures including some clustering. The number and properties of clusters may provide some measure of the level of particle-turbulence and particle–particle interactions in the high-velocity air–water flows. Herein a re-analysis of air–water clusters was applied to a highly aerated free-surface flow data set (Chanson and Carosi, Exp Fluids 42:385–401, 2007). A two-dimensional cluster analysis was introduced combining a longitudinal clustering criterion based on near-wake effect and a side-by-side particle detection method. The results highlighted a significant number of clustered particles in the high-velocity free-surface flows. The number of bubble/droplet clusters per second and the percentage of clustered particles were significantly larger using the two-dimensional cluster analysis than those derived from earlier longitudinal detection techniques only. A number of large cluster structures were further detected. The results illustrated the complex interactions between entrained air and turbulent structures in skimming flow on a stepped spillway, and the cluster detection method may apply to other highly aerated free-surface flows.     

Improving accuracy of LiDAR-derived digital terrain models for saltmarsh management

Accurate digital elevation models of saltmarshes are crucial for both conservation and management goals. Light detection and ranging (LiDAR) is increasingly used for topographic surveys due to the ability to acquire high resolution data over spatially-extensive areas. This capability is ideally suited to saltmarsh environments, which are often vast, inaccessible systems where topographic variations can be very subtle. Derivation of surface (DSMs) (ground elevation plus vegetation) versus terrain (bare ground elevation) models (DTMs) relies on the ability of the LiDAR sensor to accurately record multiple returns. In saltmarshes however, the dense stands of low (< 1 m) vegetation commonly found precludes the acquisition of more than one return, and the resulting DTM is not different to the DSM. Establishing the offset between ground and vegetation surface in order to correct the LiDAR-derived DTM can be challenging due to the spatial variability in saltmarsh habitats. Here we show the development and application of a habitat-specific correction factor (HSCF) for the Odiel Saltmarshes using a combination of habitat object-based classification (82% overall accuracy) and ground control surveys that reduces the DTM error to within that associated with the LiDAR sensor (average error 0.1 m). We also show that the true accuracy of supplied (unmodified) DTMs can be >0.5 m in saltmarshes dominated by dense vegetation such as Spartina densiflora. In particular, global projections of sea-level rise across the next 80 years (0.18–0.59 m) significantly overlaps this accuracy margin, implying that assessments and modelling of sea-level impacts in saltmarsh systems will likely be erroneous if based on Lidar-derived DTMs. Erroneous assumptions and conclusions can result if the real accuracy of DTMs (bare ground) on vegetated saltmarshes is not considered, and the consequences of the propagation of this misinformation through to management decisions should not be over-looked.     

Material dispersion by oceanic internal waves

Internal gravity waves that are generated in the open ocean have a universal frequency spectrum, called Garrett–Munk spectrum. By initializing internal waves that satisfy the Garrett–Munk spectrum in a non-hydrostatic numerical model, we investigate the material dispersion produced by these internal waves. Three numerical experiments are designed: Exp.-1 uses a linearly stratified fluid, Exp.-2 has an upper mixed layer, and Exp.-3 incorporates a circular front into the upper mixed layer. Resorting to neutrally buoyant particles, we investigate the dispersion in terms of metrics of the relative dispersion and finite-scale Lyapunov exponent (FSLE). Exp.-1 shows that the dispersion regime produced by these internal waves is between ballistic and diffusive based on relative dispersion, and is however ballistic according to FSLE. The maximum FSLE at scales of 100 m is about 5 day\(^{-1}\), which is comparable to that calculated using ocean drifters. Exp.-2 demonstrates that internal waves can generate flows and material dispersion in an upper mixed layer. However, when mixed layer eddies are present, as in Exp.-3, the dispersion in the mixed layer is controlled by the eddies. In addition, we show that inertial oscillations do not affect the relative dispersion, but impact FSLE at scales of inertial oscillations.     

Characterization and speciation of mercury in mosses and lichens from the high-altitude Tibetan Plateau

The accumulation and species of mercury (Hg) in mosses and lichens collected from high-altitude Tibetan Plateau were studied. The altitudes of the sampling sites spanned from 1983 to 5147 m, and a total of 130 mosses and 52 lichens were analyzed. The total mercury (THg) contents in mosses and lichens were in the ranges of 13.1–273.0 and 20.2–345.9 ng/g, respectively. The average ratios of methylmercury (MeHg) in THg in mosses and lichens were 2.4 % (0.3–11.1 %) and 2.7 % (0.4–9.6 %), respectively, which were higher than those values reported in other regions. The contents of THg in both mosses and lichens were not correlated with the THg in soils (p > 0.05). The lipid contents displayed a significantly positive correlation with concentrations of MeHg in mosses (r = 0.461, p < 0.01, n = 90), but not in lichens. The correlations between Hg contents in mosses and the altitudes, latitudes and longitudes of sampling sites indicated the mountain trapping and spatial deposition of Hg in the Tibetan Plateau.     

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.     

Quantifying the effect of wind on internal wave resonance in Lake Villarrica,Chile

Lake Villarrica, located in south central Chile, has a maximum depth of 167 m and a maximum fetch of about 20 km. The lake is monomictic, with a seasonal thermocline located at a depth of approximately 20 m. Field data show the presence of basin-scale internal waves that are forced by daily winds and affected by Coriolis acceleration. A modal linear and non-linear analysis of internal waves has been used, assuming a two-layer system. The numerical simulations show good agreement with the internal wave field observations. The obtained modes were used to study the energy dissipation within the system, which is necessary to control the amplitude growth. Field data and numerical simulations identify (1) the occurrence of a horizontal mode 1 Kelvin wave, with a period of about a day that coincides with the frequency of daily winds, suggesting that this mode of the Kelvin waves is in a resonant state (subject to damping and controlled by frictional effects in the field) and (2) the presence of higher-frequency internal waves, which are excited by non-linear interactions between basin-scale internal waves. The non-linear simulation indicates that only 10 % of the dissipation rate of the Kelvin wave is because of bottom friction, while the rest 90 % represents the energy that is radiated from the Kelvin wave to other modes. Also, this study shows that modes with periods between 5 and 8 h are excited by non-linear interactions between the fundamental Kelvin wave and horizontal Poincaré-type waves. A laboratory study of the resonant interaction between a periodic forcing and the internal wave field response has also been performed, confirming the resonance for the horizontal mode 1 Kelvin wave.     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.     

Evolution of a beach nourishment project using dredged sand from navigation channel,Dunkirk, northern France

The largest beach replenishment project ever in France was completed in February 2014 in Dunkirk on the coast of northern France. A volume of 1.5 × 10⁶ m³ of sand extracted from a navigation channel was placed on the beach to build up a 150 to 300 m wide supratidal platform in front of a dike, called « Digue des Alliés », which protects several residential districts of Dunkirk from marine flooding. High resolution topographic surveys were carried out during 2½ years to monitor beach morphological changes, completed by a hydrodynamic field experiment conducted in February 2016. Approximately ?138,200 m³ of sand, corresponding to 9.2% of the initial nourishment volume, were eroded over the nourishment area in about 2 years. An obvious decrease in erosion eastward with a shift from erosion to accumulation was observed, suggesting an eastward redistribution of sand. This longshore sand drift is beneficial for the eastward beach of Malo-les-Bains where most of the recreational activities are concentrated. Hydrodynamic measurements showed that waves and wave-induced currents play a major role on the longshore sand redistribution compared to tidal flows. Strong relationships were observed between cumulative offshore wave power and beach volume change during distinct beach survey periods (R² = 0.79 to 0.87), with more significant correlations for northerly waves. A slight decrease in erosion during the second year compared to the first year after nourishment suggests that the loss of sand should decrease after an initial phase of rapid readjustment of the beach shape towards equilibrium.     

Graphene quantum dots–eggshell nanocomposite to extract polycyclic aromatic hydrocarbons in water

In solid-phase extraction, the sorbent controls selectivity, affinity, and loading capacity. Many sorbents can be used, such as functionalized resin, naphthalene, activated carbon, molecular-imprinted polymers, and nanometer-sized materials. Here, we prepared a nanocomposite made of graphene quantum dots and eggshell to extract polycyclic aromatic hydrocarbons from aqueous solutions. We studied the effect of sample volume, solution pH, amount of nanosorbent, and sample loading flow rate on the extraction efficiency. We obtained a wide dynamic linear range from 0.1 to 200 ng m/L, good linearity with r ² higher than 0.99, and low-detection limits of 5–75 pg m/L. We measured PAH in spiked water samples.     

Eddy diffusivity: a single dispersion analysis of high resolution drifters in a tidal shallow estuary

In an estuary, mixing and dispersion resulting from turbulence and small scale fluctuation has strong spatio-temporal variability which cannot be resolved in conventional hydrodynamic models while some models employs parameterizations large water bodies. This paper presents small scale diffusivity estimates from high resolution drifters sampled at 10 Hz for periods of about 4 h to resolve turbulence and shear diffusivity within a tidal shallow estuary (depth <3 m). Taylor’s diffusion theorem forms the basis of a first order estimate for the diffusivity scale. Diffusivity varied between 0.001 and 0.02 m²/s during the flood tide experiment. The diffusivity showed strong dependence (R² > 0.9) on the horizontal mean velocity within the channel. Enhanced diffusivity caused by shear dispersion resulting from the interaction of large scale flow with the boundary geometries was observed. Turbulence within the shallow channel showed some similarities with the boundary layer flow which include consistency with slope of 5/3 predicted by Kolmogorov’s similarity hypothesis within the inertial subrange. The diffusivities scale locally by 4/3 power law following Okubo’s scaling and the length scale scales as 3/2 power law of the time scale. The diffusivity scaling herein suggests that the modelling of small scale mixing within tidal shallow estuaries can be approached from classical turbulence scaling upon identifying pertinent parameters.     

Migration of aluminum from food contact materials to food—a health risk for consumers? Part III of III: migration of aluminum to food from camping dishes and utensils made of aluminum

Background

When cooking on a barbecue grill, consumers often use aluminum grill pans. For one, the pan catches the fats and oils that would drip into the embers causing the formation of potentially noxious smoke, and the pan also protects the food from being burned by direct heat from the coals. In addition, new aluminum products for use in ovens and grills are becoming increasingly popular. Due to their light weight and excellent heat transfer camping, utensils made of aluminum are, for example, often used by fishermen and mountain climbers. Preparing food in aluminum utensils can, however, result in migration of the aluminum to the foodstuffs.

Results/Conclusions

In this study presented here, it was found that the transfer limit of 5.00 mg/L for aluminum is not exceeded using simulants for oil or for tap water; however, with an aqueous solution of 0.5% citric acid, the limit is clearly exceeded at 638 mg/L. This means that the Tolerable Weekly Intake (TWI) is exceeded by 298% for a child weighing 15 kg and for an adult weighing 70 kg it is equivalent to 63.8% of the TWI, assuming a daily uptake of 10 mL marinade containing lemon juice over a period of 1 week. Preparation of a fish dish with a marinade containing lemon juice in camping dishes would result in the TWI being exceeded by 871% for a child weighing 15 kg and by 187% for an adult weighing 70 kg assuming a daily uptake of 250 g over a period of 1 week.

     

Modeling the annual thermal regime of lake Ohrid,Yugoslavia, using daily weather data

The time-dependent characteristics of the five weather variables which control the annual thermal response of Lake Ohrid are analyzed in detail. These are daily values for solar radiation, air temperature, humidity, wind speed and cloud cover. A simple numerical model of the lake's thermal response, forced by thermally driven density mixing, is developed and tested using observed thermal profiles for verification. The numerical model successfully reproduces the major features of the lake's thermal regime over a 6 y period from 1972 to 1977, the average root mean square value for the simulated profiles being 1.2°C with extremes of 2.2 and 0.3°C and a standard deviation of 0.4°C.     

Salt-marsh erosion and restoration in relation to flood protection on the Wadden Sea barrier island Terschelling

This paper explores the impact of erosion and restoration measures on habitat development and on wave damping by a small salt marsh nestled alongside a dike on the Wadden island of Terschelling. The aim is to advance knowledge about the benefits and possible side-effects of salt-marsh restoration. Analysis of a time series of aerial photographs from 1944 to 2010 indicates that the salt marsh decreased steadily in size after maintenance of accretion works was terminated. In the western part of the marsh, which is accessible to sheep, vegetation is low (5–15 cm) and dominated by Salicornia europaea and by Spartina anglica. In the most intensively grazed parts, vegetation is very scarce. The eastern, inaccessible part of the salt marsh is covered by dense patches of the shrubby perennial Atriplex portulacoides and Spartina anglica (15–25 cm in height). SWAN wave models show that wave height at this location is significantly affected by the areal extent of the salt marsh as well as by the vegetation. High or dense vegetation are in the models nearly as effective in damping waves (with an initial height of 0.15 and 0.5 m) as widening the salt-marsh area by 350 m. A low density of low plants, as observed in the grazed part of the marsh, has almost no wave-damping effect. Even under conditions of sea level rise, a broader salt marsh vegetated with high plants significantly affects modelled wave height. Therefore, salt-marsh restoration is an adaptation measure worth exploring, though an array of effect types must be considered.