Characterising strongly advected discharges in the initial dilution zone

Mean concentration fields of strongly advected non-buoyant discharges are characterised with a double-Gaussian assumption. Comparisons with experimental data show that the approximation provides a reasonable representation of the cross-sectional profiles. The self-similarity of these profiles enables their form to be represented by two additional parameters, one describing the relative separation of the peaks and the other the ratio of the cross-sectional spreads. Values for these additional parameters are determined from experimental data. This systematic approach to characterising the strongly advected flows provides a consistent framework for determining spreading rates and concentration ratios, such as the peak to centreline maximum and the peak to top hat. The double-Gaussian framework also provides a basis for comparisons with the CorJet and VisJet numerical models. In addition the double-Gaussian assumption is employed to interpret data obtained using the Light Attenuation technique. This is a relatively simple measuring system, which provides depth integrated concentration information. The data obtained using this technique is shown to be generally consistent with that from previous studies.     

Integral Model for Turbulent Buoyant Jets in Unbounded Stratified Flows. Part I: Single Round Jet

The mechanics of buoyant jet flows issuing with a general three-dimensional geometry into an unbounded ambient environment with uniform density or stable density stratification and under stagnant or steady sheared current conditions is investigated. An integral model is formulated for the conservation of mass, momentum, buoyancy and scalar quantities in the turbulent jet flow. The model employs an entrainment closure approach that distinguishes between the separate contributions of transverse shear (leading to jet, plume, or wake internal flow dynamics) and of azimuthal shear mechanisms (leading to advected momentum puff or thermal flow dynamics), respectively. Furthermore, it contains a quadratic law turbulent drag force mechanism as suggested by a number of recent detailed experimental investigations on the dynamics of transverse jets into crossflow. The model is validated in several stages: First, comparison with basic experimental data for the five asymptotic, self-similar stages of buoyant jet flows, i.e., the pure jet, the pure plume, the pure wake, the advected line puff, and the advected line thermal, support the choice and magnitude of the turbulent closure coefficients contained in the entrainment formulation. Second, comparison with many types of non-equilibrium flows support the proposed transition function within the entrainment relationship, and also the role of the drag force in the jet deflection dynamics. Third, a number of spatial limits of applicability have been proposed beyond which the integral model necessarily becomes invalid due to its parabolic formulation. These conditions, often related to the breakdown of the boundary layer nature of the flow, describe features such as terminal layer formation in stratification, upstream penetration in jets opposing a current, or transition to passive diffusion in a turbulent ambient shear flow. Based on all these comparisons, that include parameters such as trajectories, centerline velocities, concentrations and dilutions, the model appears to provide an accurate and reliable representation of buoyant jet physics under highly general flow conditions.     

Integral Model for Turbulent Buoyant Jets in Unbounded Stratified Flows Part 2: Plane Jet Dynamics Resulting from Multiport Diffuser Jets

An integral model for the plane buoyant jet dynamics resulting from the interaction of multiple buoyant jet effluxes spaced along a diffuser line is considered as an extension of the round jet formulation that was proposed in Part I. The receiving fluid is given by an unbounded ambient environment with uniform density or stable density stratification and under stagnant or steady sheared current conditions. Applications for this situation are primarily for submerged multiport diffusers for discharges of liquid effluents into ambient water bodies, but also for multiple cooling tower plumes and building air-conditioning. The CorJet model formulation describes the conservation of mass, momentum, buoyancy and scalar quantities in the turbulent jet flow in the plane jet geometry. It employs an entrainment closure approach that distinguishes between the separate contributions of transverse shear and of internal instability mechanisms, and contains a quadratic law turbulent pressure force mechanism. But the model formulation also includes several significant three-dimensional effects that distinguish actual diffuser installations in the water environment. These relate to local merging processes from the individual multiple jets, to overall finite length effects affecting the plume geometry, and to bottom proximity effects given by a “leakage factor” that measures the combined affect of port height and spacing in allowing the ambient flow to pass through the diffuser line in order to provide sufficient entrainment flow for the mixing downstream from the diffuser. The model is validated in several stages: First, comparison with experimental data for the asymptotic, self-similar stages of plane buoyant jet flows, i.e. the plane pure jet, the pure plume, the pure wake, the advected line puff, and the advected line thermal, support the choice of the turbulent closure coefficients contained in the entrainment formulation. Second, comparison with data for many types of non-equilibrium flows with a plane geometry support the proposed functional form of the entrainment relationship, and also the role of the pressure force in the jet deflection dynamics. Third, the observed behavior of the merging process from different types of multiport diffuser discharges in both stagnant and flowing ambient conditions and with stratification appears well predicted with the CorJet formulation. Fourth, a number of spatial limits of applicability, relating to terminal layer formation in stratification or transition to passive diffusion in a turbulent ambient shear flow, have been proposed. In sum, the CorJet integral model appears to provide a mechanistically sound, accurate and reliable representation of complex buoyant jet mixing processes, provided the condition of an unbounded receiving fluid is satisfied.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

A numerical model of wave- and current-driven nutrient uptake by coral reef communities

We developed a numerical model capable of simulating the spatial zonation of nutrient uptake in coral reef systems driven by hydrodynamic forcing (both from waves and currents). Relationships between nutrient uptake and bed stress derived from flume and field studies were added to a four-component biogeochemical model embedded within a three-dimensional (3-D) hydrodynamic ocean model coupled to a numerical wave model. The performance of the resulting coupled physical-biogeochemical model was first evaluated in an idealized one-dimensional (1-D) channel for both a pure current and a combined wave-current flow. Waves in the channel were represented by an oscillatory flow with constant amplitude and frequency. The simulated nutrient concentrations were in good agreement with the analytical solution for nutrient depletion along a uniform channel, as well as with existing observations of phosphate uptake across a real reef flat. We then applied this integrated model to investigate more complex two-dimensional (2-D) nutrient dynamics, firstly to an idealized coral reef-lagoon morphology, and secondly to a realistic section of Ningaloo Reef in Western Australia, where nutrients were advected into the domain via alongshore coastal currents. Both the idealized reef and Ningaloo Reef simulations showed similar patterns of maximum uptake rates on the shallow forereef and reef crest, and with nutrient concentration decreasing as water flowed over the reef flat. As a result of the cumulative outflow of nutrient-depleted water exiting the reef channels and then being advected down the coast by alongshore currents, both reef simulations exhibited substantial alongshore variation in nutrient concentrations. The coupled models successfully reproduced the observed spatial-variability in nitrate concentration across the Ningaloo Reef system.     

Removing the boundary influence on negatively buoyant jets

A comprehensive laboratory study of negatively buoyant discharges is presented. Unlike previous studies, here the focus is on generating data sets where influences of the bottom boundary have been eliminated. There are significant discrepancies in the published dilution data for these flows and a contributing factor is the large variation in the bottom boundary condition. A Laser-induced Fluorescence system is employed to gather flow spread, peak concentration (minimum dilution) and trajectory data for a wide range of densimetric Froude numbers and initial discharge angles. Data from these experiments are compared with previously published data, along with predictions from integral models and a revised form of the previously published semi-analytical solutions. The new data sets are not distorted by mixing processes associated with the bottom boundary and therefore provide the basis for more meaningful assessments of the predictive capabilities of existing models, given that the influences of the bottom boundary on contaminant mixing are not incorporated into these models. In general the models assessed are able to predict key geometric quantities with reasonable accuracy, but their minimum dilution predictions are conservative. Importantly dilution at the return point shows a strong dependence on the initial discharge angle and this could have important implications for the design of discharge systems.     

Formation of tidal starting-jet vortices through idealized barotropic inlets with finite length

This paper presents a surface particle image velocimetry study to investigate the dynamics of shallow starting-jet dipoles formed by tidal flow through inlets and their interaction with vorticity formed at the inlet channel lateral boundaries. Vortical structure in the flow field is identified using a local swirl strength criterion evaluated from the two-dimensional flow field. The starting jet dipole vortices and vortices formed as the lateral boundary layers are expelled during flow reversal are characterized by their trajectory, size, and circulation. Using these quantities, a model is developed to predict the size and strength of the expelled lateral boundary layer vortices based on the inlet velocity, channel length, and width of the lateral boundary layer. The expelled boundary layer vortices are found to disrupt the formation of the primary tidal jet dipole through two mechanisms. First, because the boundary layer vortices themselves form a dipole with each half of the starting-jet dipole, the starting-jet vortices are pulled apart and advected away from the inlet mouth early in the tidal cycle, resulting in a reduction in the spin-up time and the amount of vorticity input during starting-jet vortex formation. Second, the advection of each dipole away from the inlet disconnects each starting-jet vortex from the starting jet; hence, the vortices are not fed by fluid in the jet or energized by shear in the jet boundary layers. These influences of the lateral boundary layer on the starting-jet vortices’ formation and propagation are found to be a function of the channel length L, maximum velocity U, and tidal period T, resulting in a predictive value to characterize their trajectory, strength, and evolution.     

Structural design of a natural windbreak using computational and experimental modeling

Computational and experimental approaches are employed for the structural design of a natural windbreak. It is intended to find the optimum tree shelterbelt to obviate the uneven wind speed distribution across the width dimension of a high-level competition rowing channel. The experimental results, obtained in a wind tunnel, and consisting of erosion-technique images and local wind-speed measurements, are used to benchmark the computational model. A good agreement between the two sets of results is obtained. Several windbreak configurations, considering one or two rows of different cross-sectional shape and porosity, are computationally modeled. For the shortest row a rectangular shape, with porosity of 35%, is considered; for the tallest row, which aims the modeling of a poplar tree, a porosity of 87% is assumed at the trunk level, and 60% at the crown. The optimum shelterbelt consists of two rows, composed by bamboo and poplar trees, which allows the attainment of a low and nearly uniform wind flow across the width of the channel.     

k-ε Predictions of the initial mixing of desalination discharges

Predictions from a k-ε model are compared with recently acquired experimental data from inclined negatively buoyant discharges. The k-ε model is part of a standard computational fluid dynamics package (CFX). Two approaches are taken when implementing the model. One involves using an essentially standard form of the model to predict flow behaviour. The other approach involves calibrating the model, through adjustment of the turbulent Schmidt number in the tracer transport equation, to achieve reasonable predictions for positively buoyant vertical discharges and then applying it to inclined negatively buoyant discharges. While the calibrated approach improves the predictions of some bulk parameters (notably the tracer spread and dilution) when compared to predictions from the standard model, the overall effect on the quality of the predictions is small. Comparisons with experimental data indicate that predictions from both the standard and calibrated simulations compare favourably with trajectory data, but integrated dilution predictions at the centreline maximum height are conservative (mean-integrated concentrations are over-predicted). The standard and calibrated k-ε predictions confirm the importance of buoyant instabilities on the lower (inner) side of the flow, the effects of which are clearly evident in the mean concentration profiles. However, these simulations have a tendency to overestimate the influence of stabilizing density gradients on the upper (outer) side of the flow and are unable to effectively predict the cross-sectional distribution of a tracer. In contrast to a previous study, the above comparisons indicate that predictions of bulk parameters from such models can be poor and indeed are no better than those obtained from relatively simple analytical solutions.     

Animal abundance estimation in independent observer line transect surveys

The theory of conventional line transect surveys is based on an essential assumption that 100% detection of animals right on the transect lines can be achieved. When this assumption fails, independent observer line transect surveys are used. This paper proposes a general approach, based on a conditional likelihood, which can be carried out either parametrically or nonparametrically, to estimate the abundance of non-clustered biological populations using data collected from independent observer line transect surveys. A nonparametric estimator is specifically proposed which combines the conditional likelihood and the kernel smoothing method. It has the advantage that it allows the data themselves to dictate the form of the detection function, free of any subjective choice. The bias and the variance of the nonparametric estimator are given. Its asymptotic normality is established which enables construction of confidence intervals. A simulation study shows that the proposed estimator has good empirical performance, and the confidence intervals have good coverage accuracy.     

Numerical modelling of horizontal sediment-laden jets

Sediment-laden turbulent flows are commonly encountered in natural and engineered environments. It is well known that turbulence generates fluctuations to the particle motion, resulting in modulation of the particle settling velocity. A novel stochastic particle tracking model is developed to predict the particle settling out and deposition from a sediment-laden jet. Particle velocity fluctuations in the jet flow are modelled from a Lagrangian velocity autocorrelation function that incorporates the physical mechanism leading to a reduction of settling velocity. The model is first applied to study the settling velocity modulation in a homogeneous turbulence field. Consistent with basic experiments using grid-generated turbulence and computational fluid dynamics (CFD) calculations, the model predicts that the apparent settling velocity can be reduced by as much as 30 % of the stillwater settling velocity. Using analytical solution for the jet mean flow and semi-empirical RMS turbulent velocity fluctuation and dissipation rate profiles derived from CFD predictions, model predictions of the sediment deposition and cross-sectional concentration profiles of horizontal sediment-laden jets are in excellent agreement with data. Unlike CFD calculations of sediment fall out and deposition from a jet flow, the present method does not require any a priori adjustment of particle settling velocity.     

Turbulence suppression by suspended sediment within a geophysical flow

Experiments are performed in a mixing box to evaluate the effect of suspended sediment on turbulence generated by an oscillating grid. Quartz-density sand of varying sizes and concentrations is used, and particle image velocimetry is employed to quantify only the fluid phase. Results show that (1) while a relatively large secondary flow field is present in the box, turbulence is a maximum near the grid and it decreases systematically toward the water surface; (2) relatively high concentrations of fine sediment can markedly alter this secondary flow field and significantly decrease both the time-mean and turbulent kinetic energy within the flow, yet these same sediment concentrations have little effect on the integral time and length scales derived for each velocity component; and (3) the overall turbulence suppression observed can be related to the transfer of energy from the fluid to the sediment and the maintenance of a suspended sediment load rather than commonly employed turbulence modulation criteria. These experimental data demonstrate unequivocally that the presence of a suspended sediment load can significantly reduce overall turbulent kinetic energy, and these results should be applicable to a range of sediment-laden geophysical flows.     

Enhancing transverse mixing in streams with submerged vanes

In order to maintain the water quality of moving streams, it is essential to know the process of pollutant mixing. The transverse mixing is very important which is needed to be modeled to understand mixing phenomenon. It was observed that transverse mixing is a strong function of secondary currents, thus, submerged vanes, which are aerofoil skewed at angle of 10°–40° with respect to flow, generate transverse circulations that can be utilized to induce secondary currents in the flow to enhance transverse mixing. Present study is an attempt to utilize submerged vanes as an instrument to enhance the transverse mixing by incorporating various vane configurations. In order to study the effect of vane generated circulations on transverse mixing, experimentations were conducted on three vane sizes and for various row arrangements of vanes attached to bed. An attempt is made to investigate the effect of submerged vane size and rows on transverse velocity, concentration profile and transverse mixing coefficient. It was observed by measurement of concentration profile that transverse mixing was more enhanced for submerged vanes of higher height. It was also observed that as the number of rows is proportional to the transverse mixing. By measuring the transverse velocity profile, it was observed that more and more fluid was advected in transverse direction for higher rows of vanes. By utilizing the observed transverse mixing coefficients, number of vane rows and relative height of vane, a predictor was derived to predict transverse mixing coefficient in the presence of submerged vane rows. It was observed that the derived predictor shows a fair amount of agreement in the result predicted.     

Double-Arranged Mound-Mounted Shelterbelts: Influence of Porosity on Wind Reduction between the Shelters

Using a two-dimensional LDA-System, measurements in an atmospheric boundary layer wind tunnel have been made in the intermediate field of double-arranged, mound-mounted shelterbelts. The porosity of the shelterbelts ranged from 0% to 52%. For each porosity, seven different distances of the two windbreaks between 5 and 25 times the windbreak height were investigated. The efficiency of the windbreak systems is assessed by means of protection parameters considering the change of horizontal momentum on bodies immersed in the flow, of pedestrian comfort conditions and of momentum exchange in vertical direction. Furthermore recirculation zones and the change of kinetic energy are shown. The measurement results are given in the form of contour plots indicating flow quantities and protection volumes between the two windbreaks. The shelter efficiency of the investigated windbreak systems is compared to that of adequate single line windbreaks. It is shown that for certain windbreak configurations the wind-protected area in the intermediate field of consecutively arranged windbreaks becomes larger in comparison to a single windbreak.     

Genetic structure in lagoons: the effects of habitat discontinuity and low dispersal ability on populations of <Emphasis Type="Italic">Atherina boyeri</Emphasis>

Atherina boyeri is a common inhabitant of the Mediterranean brackish and coastal waters. Several recent studies proposed the recognition of three distinct forms/species (one lagoon/freshwater and two marine) within an A. boyeri complex. In this study, we used 10 newly available microsatellite loci to better define the genetic structure of the lagoon form. A total of 428 fishes from 11 lagoons along the Italian and Albanian coasts and three Italian lakes were analyzed. A pronounced genetic spatial structure was confirmed for all considered geographic scales (near, medium and broad). A restricted gene flow was detected, suggesting reduced population connectivity. The comparison of present data and previous studies, based on mtDNA and RAPD markers, allowed a preliminary identification of the roles of recent factors (gene flow, habitat discontinuity and water currents) and past events (e.g., Pleistocene glaciations) in shaping the pattern of genetic differentiation.     

A little movement orientated to the geomagnetic field makes a big difference in strong flows

Whilst a range of animals have been shown to respond behaviourally to components of the Earth’s magnetic field, evidence of the value of this sensory perception for small animals advected by strong flows (wind/ocean currents) is equivocal. We added geomagnetic directional swimming behaviour for North Atlantic loggerhead turtle hatchlings (Caretta caretta) into a high-resolution (1/4°) global general circulation ocean model to simulate 2,925-year-long hatchling trajectories comprising 355,875 locations. A little directional swimming (1–3 h per day) had a major impact on trajectories; simulated hatchlings travelled further south into warmer water. As a result, thermal elevation of hatchling metabolic rates was estimated to be between 63.3 and 114.5% after 220 days. We show that even small animals in strong flows can benefit from geomagnetic orientation and thus the potential implications of directional swimming for other taxa may be broad.     

Optimal investment in pollution control capital in a neoclassical growth context

A neoclassical growth model is employed to study the optimal division of investment between pollution control capital and directly productive capital while simultaneously treating the saving rate as a policy variable to be varied to the optimal. The model treats pollution as a flow positively related to aggregate output, negatively related to the stock of pollution control capital, and having a negative effect on utility. Under the assumptions chosen, the optimal policy is an unbalanced one of specializing investment first in directly productive capital and subsequently in pollution control capital.     

Suspended solids within sewers: an experimental study

A precise determination of the pollutant loads processed through the sewer networks is required for legal, environmental and management reasons. Such a determination requires a good evaluation of concentrations, which is not an easy task because samplers or sensors can perform measurements in a cross-sectional area of only a few locations, and this local information must then be extrapolated across the whole cross-sectional area of the sewer network. This latter step is problematical. A great deal of attention must be given to suspended solids because they are prone to settle and display heterogeneous concentration patterns. A two-dimensional sampler has been designed to measure the variation of suspended solids concentration throughout a cross-sectional area. This paper presents the principles upon which the conception of this experimental device was based and the building process. Experimental results are presented and discussed.     

Measurements and modeling of open-channel flows with finite semi-rigid vegetation patches

The hydrodynamics of flows through a finite length semi-rigid vegetation patch (VP) were investigated experimentally and numerically. Detailed measurements have been carried out to determine the spatial variation of velocity and turbulence profiles within the VP. The measurement results show that an intrusion region exists in which the peak Reynolds stress remains near the bed. The velocity profile is invariant within the downstream part of the VP while the Reynolds stress profile requires a longer distance to attain the spatially invariant state. Higher vegetation density leads to a shorter adjustment length of the transition region, and a higher turbulence level within the VP. The vegetation density used in the present study permits the passing through of water and causes the peak Reynolds stress and turbulence kinetic energy each the maximum at the downstream end of the patch. A 3D Reynolds-averaged Navier–Stokes model incorporating the Spalart–Allmaras turbulence closure was employed subsequently to replicate the flow development within the VP. The model reproduced transitional flow characteristics well and the results are in good agreement with the experimental data. Additional numerical experiments show that the adjustment length can be scaled by the water depth, mean velocity and maximum shear stress. Empirical equations of the adjustment lengths for mean velocity and Reynolds stress were derived with coefficients quantified from the numerical simulation results.     

Trophic structure of coastal Antarctic food webs associated with changes in sea ice and food supply

Predicting the dynamics of ecosystems requires an understanding of how trophic interactions respond to environmental change. In Antarctic marine ecosystems, food web dynamics are inextricably linked to sea ice conditions that affect the nature and magnitude of primary food sources available to higher trophic levels. Recent attention on the changing sea ice conditions in polar seas highlights the need to better understand how marine food webs respond to changes in such broad-scale environmental drivers. This study investigated the importance of sea ice and advected primary food sources to the structure of benthic food webs in coastal Antarctica. We compared the isotopic composition of several seafloor taxa (including primary producers and invertebrates with a variety of feeding modes) that are widely distributed in the Antarctic. We assessed shifts in the trophic role of numerically dominant benthic omnivores at five coastal Ross Sea locations. These locations vary in primary productivity and food availability, due to their different levels of sea ice cover, and proximity to polynyas and advected primary production. The delta15N signatures and isotope mixing model results for the bivalves Laternula elliptica and Adamussium colbecki and the urchin Sterechinus neumeyeri indicate a shift from consumption of a higher proportion of detritus at locations with more permanent sea ice in the south to more freshly produced algal material associated with proximity to ice-free water in the north and east. The detrital pathways utilized by many benthic species may act to dampen the impacts of large seasonal fluctuations in the availability of primary production. The limiting relationship between sea ice distribution and in situ primary productivity emphasizes the role of connectivity and spatial subsidies of organic matter in fueling the food web. Our results begin to provide a basis for predicting how benthic ecosystems will respond to changes in sea ice persistence and extent along environmental gradients in the high Antarctic.