Idealized study on cohesive sediment flux by tidal asymmetry

The flux of cohesive sediment in an estuary is determined by many factors, including tidal asymmetry, wave effect, fluvial influence, phase difference between tidal velocity and tidal level fluctuations, sediment properties, flocculation, bed erodibility, bathymetry effect and other nonlocal effects. Our capability in predicting sediment fluxes in tide-dominant environments is critical to the morphodynamics and water quality of estuaries. Due to the difficulties in carrying out detailed measurement of sediment flux with high spatial and temporal resolutions, an one-dimensional-vertical (1DV) numerical model for cohesive sediment transport, previously verified and calibrated with field measured cohesive sediment concentration data, is utilized here to study some of the aforementioned factors in affecting tidal-driven sediment fluxes in idealized condition. Tidal-averaged sediment flux is shown to be correlated with tidal velocity skewness with a linear relationship. This linear relationship is different from that of non-cohesive sediment and it is demonstrated here to be mainly due to variable critical shear stress implemented for the mud bed in order to parameterize consolidation. The reason that tidal velocity skewness causes tidal-averaged residual sediment transport is shown to be due to nonlinear intra-tidal interactions between flow velocity and sediment concentration. Moreover, the effects of nonlinear intra-tidal interaction between tidal velocity and tidal level fluctuations is shown to mainly cause seaward transport, which is the most significant under progressive wave system (phase difference 90°) and almost negligible for standing wave system (phase difference 0°).     

Streamwise velocity profiles in coastal currents

The ability to model marine currents can be a powerful device for many planning activities, for which the knowledge of the velocity field is of pivotal importance, such as the evaluation of current-induced loading on maritime structures or the diffusion and dispersion of polluted flow discharges. Observations of time-averaged velocity profiles, taken with a vessel mounted acoustic Doppler current profiler during a monitoring survey program in the seas of Southern Italy, are analysed in this paper. The measurements were taken under non-breaking conditions, offshore the surf zone, with the aim of reproducing the vertical trends of the streamwise velocity by means of standard theoretical laws. To evaluate also the possible influence of stratification on the current velocity profile shape, together with velocity measurements water temperature and salinity were also measured at the same time and locations, by means of a CTD recorder. The examined surveys referred to different time periods and sites, to guarantee a general validity of deductions. On the basis of the experiments, we verified the actual existence of a log-layer and concluded that the upper limit of the region in which the log law is applicable extends well beyond the inner region. Moreover, the deviations of the measured velocity from the logarithmic profiles above the height of the log layer is consistent with the effects of stratification. The parameters of the log law were estimated, depending on both flow dynamics and stratification in the target area. As a second step, in the most superficial and stratified layer, the velocity profiles were modelled by means of a power law, which fitted the measured data well. According to previous studies, the power law parameters result Reynolds number dependent by means of a new proposed formulation. Finally, the bottom stress and the bottom drag coefficient were investigated.     

Mixing in the three-dimensional wake of an experimental modelled vehicle

Experimental results on tracer gas diffusion within the near wake of a simplified model car (Ahmed model with a rear slant angle of 25°) are presented. Pollutant emission is simulated using heated air injected through a small pipe at one side of the model base. Fine cold wire thermometry is used to measure instantaneous temperature excess and variance of temperature gradient in the near wake. Measurements of the three mean velocity components were made using a laser Doppler anemometers system. Characteristics of the mean and fluctuating temperature fields, time-averaged flow streamlines and scalar dissipation measurements are presented and discussed. The local mixing time scale is determined from the measured mean dissipation rate of temperature variance. Its value shows that micromixing is not a limiting phenomenon for chemical reactions in the near wake.     

Numerical simulation of water mixing and renewals in the Barcelona harbour area: the winter season

In the present paper, we use numerical simulation to investigate currents, mixing and water renewal in Barcelona harbour under typical conditions of wind forcing for the winter season. This site is of particular importance due to the interplay between touristic and commercial activities, requiring detailed and high-definition studies of water quality within the harbour. We use Large Eddy Simulation (LES) which directly resolves the anisotropic and energetic large scales of motion and parametrizes the small, dissipative, ones. Small-scale turbulence is modelled by the anisotropic Smagorinsky model (ASM) to be employed in presence of large cell anisotropy. The complexity of the harbour is modelled using a combination of curvilinear, structured, non-staggered grid and the immersed boundary method. Boundary conditions for wind and currents at the inlets of the port are obtained from in-situ measurements. Analysis of the numerical results is carried out based on both instantaneous and time-averaged velocity fields. First- and second-order statistics, such as turbulent kinetic energy and horizontal and vertical eddy viscosities, are calculated and their spatial distribution is discussed. The study shows the presence of intense current in the narrow and elongated part of the harbour together with sub-surface along-shore elongated rolling structures (with a time scale of a few hours), and they contribute to the vertical water mixing. Time-averaged velocity field reveals intense upwelling and downwelling zones along the walls of the harbour. The analysis of second-order statistics shows strong inhomogeneity of turbulent kinetic energy and horizontal and vertical eddy viscosities in the horizontal plane, with larger values in the regions characterized by stronger currents. The water renewal within the port is quantified for particular sub-domain regions, showing that the complexity of the harbour is such that certain in-harbour basins have a water renewal of over five days, including the yacht marina area. The LES solution compares favourably with available current-meter data. The LES solution is also compared with a RANS solution obtained in literature for the same site under the same forcing conditions, the comparison demonstrating a large sensitivity of properties to model resolution and frictional parametrization.     

Generalized extreme value regression for ordinal response data

This paper introduces a flexible skewed link function for modeling ordinal response data with covariates based on the generalized extreme value (GEV) distribution. Commonly used probit, logit and complementary log-log links are prone to link misspecification because of their fixed skewness. The GEV link is flexible in fitting the skewness in the response curve with a free shape parameter. Using Bayesian methodology, it automatically detects the skewness in the response curve along with the model fitting. The flexibility of the proposed model is illustrated by its application to an ecological survey data about the coverage of Berberis thunbergii in New England. We employ the latent variable approach by Albert and Chib (J Am Stat Assoc 88:669–679, (1993) to develop computational schemes. For model selection, we employ the Deviance Information Criterion (DIC).     

On interacting bee/mite populations: a stochastic model with analysis using cumulant truncation

Parasitism by the Varroa mite has had recent drastic impact on both managed and feral bee colonies. This paper proposes a stochastic population dynamics model for interacting African bee colony and Varroa mite populations. Cumulant truncation procedures are used to obtain approximate transient cumulant functions, unconstrained by the usual assumption of bivariate Normality, for an assumed large-scale model. The apparent size of the variance and skewness functions suggest the importance of the proposed truncation procedure which retains some higher-order cumulants, but determining the accuracy of the approximations is problematical. A smaller-scale bee/Varroa mite model is hence proposed and investigated. The accuracy for the means is exceptional, for the second-order cumulants is moderate, and for some third-order cumulants is poor. Notwithstanding the poor accuracy of a skewness approximation, the saddlepoint approximations for the marginal transient population size distributions are excellent. The cumulant truncation methodology is very general, and research is continuing in its application to this new class of host-parasite models.     

A simple model for simulation of insect pheromone dispersion within forest canopies

Synthetic pheromones and other behavioral chemicals are used by land managers to prevent insect-caused tree mortality or crop failure in forest and agricultural systems. Currently, no method exists to continuously measure pheromone concentration or movement in real-time. To improve our understanding of pheromone fate and transport under different forest canopies, results from a set of surrogate pheromone (sulfur hexafluoride tracer) experimental trials were used to evaluate a simple, instantaneous, three-dimensional Lagrangian dispersion model. The model was designed to predict both instantaneous and time-averaged pheromone concentrations. Overall, the results from the model show simulated time-averaged arc maximum concentrations within a factor of two of the observed data. The model correctly matched the sharp peaks and narrow widths of the meandering plumes observed in the instantaneous data, however the magnitude of the instantaneous peaks was often under-estimated. This model and evaluation provide the basis for a tool that can be used to guide deployment of synthetic pheromones or other semiochemicals for monitoring, mass trapping, or disruption of mating or aggregation.     

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.     

Dressler’s theory for curved topography flows: iterative derivation,transcritical flow solutions and higher-order wave-type equations

The Dressler equations are a system of two non-linear partial differential equations for shallow fluid flows over curved topography. The theory originated from an asymptotic stretching method formulating the equations of motion in terrain-fitted curvilinear coordinates. Apparently, these equations failed to produce a transcritical flow profile changing from sub- to supercritical flow conditions. Further, wave-like motions over a flat bottom are excluded because the bed-normal velocity component is not accounted for. However, the theory was found relevant for several environmental flow problems including density currents over mountains and valleys, seepage flow in hillslope hydrology, the development of antidunes, the formation of geological deposits from hyper-concentrated flows, and shallow-water flow modeling in hydraulics. In this work, Dressler’s theory is developed in an alternative way by a systematic iteration of the stream and potential functions in terrain-fitted coordinates. The first iteration was found to be the former Dressler’s theory, whereas a second iteration of the governing equations results in velocity components generalizing Dressler’s theory to wave-like motion. Dressler’s first-order theory produces a transcritical flow solution over topography only if the total head is fixed by a minimum value of the specific energy at the transition point. However, the theory deviates from measurements under subcritical flow conditions, given that the bed-normal velocity component is significant. A second iteration to the velocity field was used to produce a second-order differential equation that resembles the cnoidal-wave theory. It accurately produces flow over an obstacle including the critical point and the minimum specific energy as part of the numerical solution. The new cnoidal-wave model compares well with the theory of a Cosserat surface for directed fluid sheets, whereas the Saint-Venant theory appears to be poor.     

Turbulence characteristics within sparse and dense canopies

Boundary layer interactions with canopies control various environmental processes. In the case of dense and homogeneous canopies, the so-called mixing layer analogy is most generally used. When the canopy becomes sparser, a transition occurs between the mixing layer and the boundary layer perturbed by interactions between element wakes. This transition has still to be fully understood and characterized. The experimental work presented here deals with the effect of the canopy density on the flow turbulence and involves an artificial canopy placed in a fully developed turbulent boundary layer. One and two-component velocity measurements are performed, both within and above the canopy. The influence of the spacing between canopy elements is studied. Longitudinal velocity statistical moments and Reynolds stresses are calculated and compared to literature data. For spacings greater than the canopy height, evidences of this transition are found in the evolution of the skewness factor, shear length scale and mixing length.     

Velocity measurements in inclined negatively buoyant jets

Experiments were performed with a particle tracking velocimetry system to investigate the behaviour of inclined negatively buoyant jets with source angles of 15°, 30°, 45°, 60°, 65°, 70°, and 75° in stationary ambient conditions. Velocities were measured in a plane aligned with the central axis of the flow and the experiments were designed such that the flow did not interact with boundaries in the region were the flow behaviour was measured. The results of this study complement previous research, which has largely focused on the mean geometric characteristics and the mean dilution of the discharged fluid. Geometric characteristics, spreading rates, and time-averaged (mean) centreline velocity results are compared with relevant experimental results from previous studies and integral model predictions. Axial and transverse mean velocity profiles at maximum height and the return point provide additional insights into the detrainment of discharged fluid due to the unstable density gradient on the inner side of the flow.     

A modification of the classical Ekman model accounting for the Stokes drift and stratification effects

A modification of the classical Ekman model of oceanic wind-driven currents including the Stokes drift and stratification effects is discussed. The modification is formulated as an application of turbulence mechanics accounting for the curvature effect of velocity fluctuation streamlines. It is shown that similar to the Stokes drift effect, the presence of a density jump layer (pycnocline) decreases the veering of the flow velocity vector at the surface from the direction of the wind stress. It is shown also that in the pycnocline the decrease of the norm of the velocity vector as well as its rotation with depth is smaller than in the regions adjacent to the pycnocline. If the Stokes drift and stratification effects are neglected, the model reduces to the classical Ekman solution with the coefficient of the turbulent shear viscosity replaced by an effective viscosity coefficient. The vertical distributions of velocity predicted by the modified model are compared with the velocity data measured in the Drake Passage and within the Long-Term Upper Ocean Study (LOTUS) in the North Atlantic.     

A Numerical and Field Investigation of Surface Heat Fluxes from Small Wind-Sheltered Waterbodies in Semi-Arid Western Australia

Shelterbelts are used for a variety of purposes in agricultural environments, primarily because of their ability to improve the downwind microclimate. Excessive evaporative losses from small, agricultural water supply reservoirs in semi-arid Western Australia motivated a combined numerical modelling and field investigation into the potential for using shelterbelts to reduce evaporation from these open waterbodies. A numerical model of the disturbed momentum and turbulence fields in the region modified by the wind-shelter was employed and accounted for the presence of a waterbody downwind. The model was coupled with conservation equations for heat and moisture and sensible and latent heat fluxes were estimated from the simulated momentum, temperature and humidity fields. The numerical simulations were tested against four days of field data from two experiments conducted in the agricultural districts of southwest Western Australia that measured boundary-layer evolution over a variety of small waterbodies protected by artifical and natural wind-shelters. The model provided good predictions of windspeed during neutral conditions, but inadequate specification of the upwind boundary during non-neutral stabilities resulted in the model failing to capture any sensitivity to atmospheric stability as seen in the field data. Despite this limitation, the temperature and humidity fields were adequately captured by the model, and evaporative mass flux predictions also agreed well with estimates taken from water-balance measurements. It is concluded that well-designed wind-shelters can reduce evaporation from open waterbodies by 20–30% as a result of reductions in the velocity scales responsible for removing moisture from the water surface. The model can be used to estimate the values of various shelterbelt design parameters (e.g., porosity, height) that could be applied in the field to provide optimum evaporation reductions.     

Analytical expression for predicting the reduced settling velocity of small particles in turbulence

The influence of turbulence on the settling velocity of small particles remains an inconclusive research subject. Both enhanced and retarded particle settling compared to quiescent settling have been reported in previous literature, and several theories have been proposed. Among the mechanisms that account for reduced settling velocities in turbulence, the loitering effect is an important one that generally exists in various conditions. This study focuses on an analytical prediction of reduced particle settling velocities due to the loitering effect. By considering the velocity autocorrelation function as a step function analogous to the free path theory, and using a modified integral time scale to capture the essence of the loitering effect, an analytical expression is derived for predicting the reduction in particle settling velocities in turbulent flows. Calculation results of the expression are then examined by comparing with the results of a random walk model, direct numerical simulations that conditionally captured the reduction of particle settling velocities due to the loitering effect, and several representative experiments. Major possible influencing factors on the prediction of the analytical expression and applicable conditions of the expression are then further discussed. The proposed analytical expression is shown to be suitable for predicting the reduced settling velocities of small particles with relatively weak inertia in turbulent flows and could provide a reasonable explanation for reported cases in which small to moderate reductions in particle settling velocities were observed.

     

基于遥感技术的太湖放养凤眼莲的生长模型

根据从卫星遥感影像或全球定位系统获得的凤眼莲面积变化数据,对在开放水域中大面积放养凤眼莲面积变化的研究,建立了放养凤眼莲生长的Logistic动力学模型来估算凤眼莲面积。在太湖竺山湾凤眼莲放养试验表明：从2009年8月4日到10月15日,凤眼莲面积从54 250.34 m2增加到135 525.66 m2,建立的估计凤眼莲面积的Logistic模型的环境容量k s为2 697 210.03,面积增长率rs为0.012 4,对试验数据的拟合优度检验R 2为0.975 3,F检验为78.93,均方根误差RMSE为2 705.55,表明该模型准确模拟了凤眼莲生长的动力学特征。     

Analytical solution of two-dimensional advection–dispersion equation with spatio-temporal coefficients for point sources in an infinite medium using Green’s function method

In the present study analytical solutions of a two-dimensional advection–dispersion equation (ADE) with spatially and temporally dependent longitudinal and lateral components of the dispersion coefficient and velocity are obtained using Green’s Function Method (GFM). These solutions describe solute transport in infinite horizontal groundwater flow, assimilating the spatio-temporal dependence of transport properties, dependence of dispersion coefficient on velocity, and the particulate heterogeneity of the aquifer. The solution is obtained in the general form of temporal dependence and the source term, from which solutions for instantaneous and continuous point sources are derived. The spatial dependence of groundwater velocity is considered non-homogeneous linear, whereas the dispersion coefficient is considered proportional to the square of spatial dependence of velocity. An asymptotically increasing temporal function is considered to illustrate the proposed solutions. The solutions are validated with the existing solutions derived from the proposed solutions in three special cases. The effect of spatially/temporally dependent heterogeneity on the solute transport is also demonstrated. To use the GFM, the ADE with spatio-temporally dependent coefficients is reduced to a dispersion equation with constant coefficients in terms of new position variables introduced through properly developed coordinate transformation equations. Also, a new time variable is introduced through a known transformation.     

Applying the discrete vortex method in environmental fluid mechanics: A study of the time-averaged near wake behind a circular cylinder

This work illustrates the discrete vortex method (DVM) as a tool for simulating environmental fluid mechanics problems involving transport in the wake of a bluff body. The DVM was used to model both the long-time-averaged and instantaneous features of flow past a circular cylinder. Simulations were performed for Re = 140, 000. Verification testing was accomplished by refining time-step and vortex element circulation. The DVM was validated through comparison with experimental data from Cantwell and Coles. Verification testing demonstrated that, while global convergence is not possible for an unsteady flow simulation, it is possible to have convergence to physical results. This verification entails identifying a range of parameters in which a flow containing coherent structures and matching physical conditions is attainable. Validation tests demonstrated excellent agreement between experimental and simulated results for time-averaged velocity and shear stress profiles, as well as Strouhal number.     

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.     

A proposed methodology for assessing the quality of wildlife habitat

Measurement of the quality and quantity of wildlife habitat is a necessity in planning and analyzing water resource projects. A methodology that provides a means of examining the following major components of habitat quality is proposed: (1) the quantity of land uses, (2) the degree of interspersion of land use, aand (3) the stage of land management and vegetative types. The model can be calibrated to a specific region by experienced wildlife biologists and applied by field technicians using aerial photographs and field evaluation of randomly located points. Biologists familiar with a region develop transformation curves relating variables that represent the above three major components to standardizing factors that range from 0 to 1; these transformation curves reflect the sensitivity of wildlife populations to determinants of habitat quality. The methodology also provides the means for biologists to weigh each components as to its relative importance to the wildlife group or species. A weighted geometric mean of the components provides an indication of the overall quality of the habitat. Changes in the habitat components that may occur because of alternative project actions can be used to evaluate the effectiveness of the alternatives. The change in the index of habitat quality will indicate the need for mitigation in water resource projects. The system then can be used to identify areas for mitigation and the effects of mitigating measures.