Estimation of the surface stress from the streamwise pressure gradient: The Kármán Integral Momentum Equation revisited

A method is developed to estimate the stress at the surface in a portable wind tunnel for wind erosion studies. The boundary layer height and the pressure gradient are used in a simple expression from the Kármán Integral Momentum Equation. Values of friction velocity u _* are within 10% of experimental values obtained through correlation techniques, including measurements of differential pressures with the Murdoch Turbulence Probe MTP and the X-wire, hot-wire anemometer XWA. Wind velocity and stress profiles reveal logarithmic trends and a constant stress layer near the surface in the DAWA portable wind tunnel. Realignment of the statistics with the mean wind is essential.     

Near-boundary velocity and turbulence in depth-varying stream flows

This research examined the temporal distribution of turbulent structure near a streambank toe through the progression of a flood wave in West Run (Morgantown, WV, USA). Three-dimensional velocities and water depths were measured through a 17-h flood event. Turbulence characteristics were examined: Reynolds stresses, turbulent kinetic energy, and turbulence intensities. On average, near-boundary velocity during the rising stage was less than the falling stage, likely due to the measurement location and local roughness. The velocity vectors shifted from towards bed before the flood wave to toward the streambank during progression of the flood wave. Turbulent kinetic energy increased with increasing water depth during the rising stage. Reynolds stress, τ_xz, increased with increasing water depth during the rising stage, but the majority of the stresses were negative through the storm event. Reynolds stress, τ_xy, was positive throughout the event and did not vary with depth. This work is among the first to evaluate turbulence during depth-varying flows in the field.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

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.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Turbulent velocity profile in fully-developed open channel flows

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

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.     

Advection and resulting CO2 exchange uncertainty in a tall forest in central Germany

Potential losses by advection were estimated at Hainich Forest, Thuringia, Germany, where the tower is located at a gentle slope. Three approaches were used: (1) comparing nighttime eddy covariance fluxes to an independent value of total ecosystem respiration by bottom-up modeling of the underlying processes, (2) direct measurements of a horizontal CO2 gradient and horizontal wind speed at 2 m height in order to calculate horizontal advection, and (3) direct measurements of a vertical CO2 gradient and a three-dimensional wind profile in order to calculate vertical advection. In the first approach, nighttime eddy covariance measurements were compared to independent values of total ecosystem respiration by means of bottom-up modeling of the underlying biological processes. Turbulent fluxes and storage term were normalized to the fluxes calculated by the bottom-up model. Below a u(*) threshold of 0.6 m/s the normalized turbulent fluxes decreased with decreasing u(*), but the flux to the storage increased only up to values less than 20% of the modeled flux at low turbulence. Horizontal advection was measured by a horizontal CO2 gradient over a distance of 130 m combined with horizontal wind speed measurements. Horizontal advection occurred at most of the evenings independently of friction velocity above the canopy. Nevertheless, horizontal advection was higher when u(*) was low. The peaks of horizontal advection correlated with changes in temperature. A full mass balance including turbulent fluxes, storage, and horizontal and vertical advection resulted in an increase of spikes and scatter but seemed to generally improve the results from the flux measurements. The comparison of flux data with independent bottom-up modeling results as well as the direct measurements resulted in strong indications that katabatic flows along the hill slope during evening and night reduces the measured apparent ecosystem respiration rate. In addition, anabatic flows may occur during the morning. We conclude that direct measurements of horizontal and vertical advection are highly necessary at sites located even on gentle hill slopes.     

The Modelling of Turbulence from Traffic in Urban Dispersion Models — Part II: Evaluation Against Laboratory and Full-Scale Concentration Measurements in Street Canyons

The paper addresses the problem of the parameterisation of traffic induced turbulent motion in urban dispersion models. Results from a variety of full-scale and wind-tunnel studies are analysed and interpreted within a modelling framework based on scaling considerations. The combined effects of traffic and wind induced dispersive motions are quantified for different traffic situations (variable traffic densities, vehicle velocities and vehicle types) and incorporated into the developed parameterisations. A new dispersive velocity scale is formulated and recommendations regarding its application in urban dispersion models are given. The necessity of accounting for traffic induced air motions in predictions of street-canyon pollution levels is demonstrated. Further research is needed to verify the empirical constants in the proposed parameterisations and to generalize the developed approach for a broader range of urban building configurations, meteorological conditions, and traffic situations.     

Turbulence Perturbations in the Neutrally Stratified Surface Layer due to the Interaction of a Katabatic Flow with a Steep Ridge

Turbulence measurements were performed in Antarctica, on the Nansen Ice Sheet, dominated by westerly katabatic winds. The measurements were taken at two sites aligned with the katabatic wind fall-line. The measuring stations were located in the middle of a wide, flat iced area at a distance of 14 km from the base of a sloping surface and at the top of a steep ridge (Inexpressible Island). The aim was to investigate the perturbation of turbulence close to the ground generated by the interaction of the flow with the ridge. We present an analysis comparing the data measured at the upstream unperturbed station with those at the top of the obstacle. Moments and spectra of velocity components have been calculated for almost steady periods. The topography and roughness change produce a combined effect on the flow acceleration (of the opposite sign) and on the turbulent stresses (of the same sign). Spectra of velocity components measured at the top of the ridge and scaled by unperturbed quantities evidence an increment of energy in the high frequency subrange with respect to the up-stream flow. Moreover, the horizontal velocity components display a shift in turbulence maximum towards higher frequencies. The vertical velocity spectrum exhibits an energy increment at low frequencies with respect to the upstream spectrum.     

The Modelling of Turbulence from Traffic in Urban Dispersion Models — Part I: Theoretical Considerations

The modelling of pollutant dispersion at the street scale in an urban environment requires the knowledge of turbulence generated by the traffic motion in streets. In this paper, a theoretical framework to estimate mechanical turbulence induced by traffic in street canyons at low wind speed conditions is established. The standard deviation of the velocity fluctuations is adopted as a measure of traffic-produced turbulence (TPT). Based on the balance between turbulent kinetic energy production and dissipation, three different parameterisations for TPT suitable for different traffic flow conditions are derived and discussed. These formulae rely on the calculations of constants that need to be estimated on the basis of experimental data. One such estimate has been made with the help of a wind tunnel data set corresponding to intermediate traffic densities, which is the most common regime, with interacting vehicle wakes.     

Turbulence closure models and their application in RAMS

Turbulence closures are fundamental for modelling the atmospheric diffusion in numerical codes and the resulting eddy diffusivities are key parameters in describing the transport and dispersion in the boundary layer. In this work, four turbulence closure schemes have been applied for reproducing a neutral flow over schematic complex terrain using the meteorological model RAMS. Two of the closures, a one-equation (E-l) and a two-equations (E-) model, have been implemented in RAMS in alternative to the ones originally available. In these cases, an analytical method based on the similarity theory for the atmospheric surface layer and boundary layer is adopted to calculate the empirical constants of the turbulence closures. Some examples of numerical studies performed to simulate the flow and turbulence over a 3-D hill in wind-tunnel experiment in neutral stratification are presented and discussed. An intercomparison of simulations related to different closures is considered by analysing the main features of the flow over the hill and by comparing calculated vertical profiles of turbulent kinetic energy with measured ones.     

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.     

An experimental analysis of lunule function in the sand dollar Mellita quinquiesperforata

Coefficients of lift and drag for fifty lunulate sand dollars [Mellita quinquiesperforata (Leske 1778)], poised over a flat, sandy surface, were determined at seven different air stream velocities. Coefficients of static friction were obtained experimentally on level sand in standing water. Critical velocity, at which a sand dollar would be dislodged, is defined as that current velocity in which the force of drag is equal to the product of residual weight (i.e. weight-lift) and the coefficient of static friction. Small individuals have lower critical velocities than large ones: 26 cm s^?1 and 43 cm s^?1 for individuals 48 and 94 mm in diameter, respectively. Measurements were repeated for 20 of the controls, their lunules smoothly filled with modelling clay, and 20 from which lunule spines were removed. In every case, blocking the lunules decreased critical velocity, by an average of 15%; this was shown to be statistically significant (P>0.01). Removal of spines, which partially occlude the lunules, significantly (P>0.01) increased critical velocity by an average of 13%. Direct measurements in a water tunnel showed that calculated values of critical velocity were slightly too high, but otherwise confirmed the role of the lunules. Allometric analysis showed that as sand dollars increase in size, the distances of the lunules from the center of the disk increase less rapidly than the radii. A principal components' analysis of profile measurements of nine species indicated that lunulate and non-lunulate sand dollars are quite different in form. Lunulate species tend to be thinedged, flat-domed and larger in diameter. Some alternative ideas of lunule function are discussed and it is shown that there are no allometric or experimental data to support any of the suggested feeding hypotheses.     

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.     

A Computational Fluid Dynamics Approach for Urban Area Transport and Dispersion Modeling

A simulation tool has been developed to model the wind fields, turbulence fields, and the dispersion of Chemical, Biological, Radiological and Nuclear (CBRN) substances in urban areas on the building to city blocks scale. A Computational Fluid Dynamics (CFD) approach has been taken that naturally accounts for critical flow and dispersion processes in urban areas, such as channeling, lofting, vertical mixing and turbulence, by solving the steady-state, Reynolds-Averaged Navier–Stokes (RANS) equations. Rapid generation of high quality cityscape volume meshes is attained by a unique voxel-based model generator that directly interfaces with common Geographic Information Systems (GIS) file formats. The flow and turbulence fields are obtained by solving the steady-state RANS equations using a collocated, pressure-based approach formulated for unstructured and polyhedral mesh elements. Turbulence modeling is based upon the Renormalization Group variant of the k–ε model (k–ε RNG). Neutrally buoyant simulations are made by prescribing velocity boundary condition profiles found by a power–law relationship, while turbulence quantities boundary conditions are defined by a prescribed mixing length in conjunction with the assumption of turbulence equilibrium. Dispersion fields are computed by solving an unsteady transport equation of a dilute gas, formulated in a Eulerian framework, using the velocity and turbulence fields found from the steady-state RANS solution. In this paper the model is explained and detailed comparisons of predicted to experimentally obtained velocity, turbulence and dispersion fields are made to neutrally stable wind tunnel and hydraulic flume experiments.     

OHP–APT 2002 gravity wave campaign: waves,turbulence and forecasts

In July 2002, a multi-sensor campaign was conducted in southern France to investigate the hypothesized connection between gravity waves and optical turbulence. A generalized scidar (GS) was mounted on the 1.93 m diameter telescope at l’Observatoire de Haute-Provence (OHP). The GS provides continuous profiles of optical turbulence with 300 m vertical resolution from telescope altitude up to 25 km. Thermosondes, which provide in situ measurement of optical turbulence by measuring temperature variance, were launched at the OHP site and at a site approximately 20 km west-northwest of OHP. Gravity wave activity was deduced from temperature and wind velocity measured by radiosondes, which are part of the thermosonde system. In this paper, gravity waves were analyzed using techniques for simple two-dimensional mountain waves, with only fair results. Mesoscale models were run at moderately high resolution for the period. The forecasts were analyzed for wave activity, and a post processor model was used to diagnose the optical turbulence. Mountain waves were evident in the forecast, but quantitative comparison showed the forecast to be inadequate in predicting wave strength. The forecast optical turbulence was in fair agreement with measurements with notable exceptions. This paper is declared a work of the U.S. Government and is not subject to copyright protection in the United States. The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged. Presented at the Turbulence and Waves in Stably Stratified Atmospheric Shear Flows: Measurements, High Resolution Simulations and Numerical Prediction Challenges Conference, 13–15 September 2004, Lighthill Institute of Mathematical Sciences, University College London, London, UK.     

Dynamics of Turbulence in Low-Speed Oscillating Bottom-Boundary Layers of Stratified Basins

This paper focuses on the impact of an oscillating low-speed current on the structure and dynamics of the bottom-boundary layer (BBL) in a small stratified basin. A set of high-resolution current profile measurements in combination with temperature-microstructure measurements were collected during a complete cycle of the internal oscillation (`seiching') in the BBL of Lake Alpnach, Switzerland. It was found that even a relatively long seiching period of 24 hours significantly changed the form of the near-bottom current profiles as well as the dynamics of the turbulent dissipation rate compared to the steady-state law-of-the-wall. A logarithmic fit to the measured current profiles starting at a distance of 0.5 m above the sediment led to inconsistent estimates of both friction velocity and roughness length. Moreover, a phase lag between the current and the turbulent dissipation of 1.5 hours and a persistent maximum in the current profile at a height of 2.5 to 3 m above the sediment were observed. The experimental findings were compared to the results of a k- turbulence model and showed good agreement in general. Specifically, the inconsistent logarithmic fitting results and the observed phase lag were reproduced well by the model.     

Numerical and wind tunnel modeling on the windbreak effectiveness to control the aeolian erosion of conical stockpiles

There are many industrial sites where open aggregate conical piles exist to store granular materials, like coal, industrial residuals, or other minerals. Usually these storage piles are placed in open areas, making them susceptible to wind erosion, which can create health, environmental, and/or economical concerns. It is common to minimize the dust emission through the placement of windbreaks in the vicinity of the storage piles, which reduces the wind speed in the vicinity of the pile’s surface. In this work, some experimental results from a wind tunnel study on the erosion of a conical sand pile, exposed or protected by a fence with porosities of 0, 70, and 83%, are shown. For the sheltered cases, the windbreak was placed at several distances from the pile’s leading edge, ranging from H to 4H, where H is the initial height of the non-eroded pile. The evolution of the shape of the sand pile is shown, at different instances in time, and the pile deformation quantified, using a novel experimental setup developed for wind erosion studies. This information might be regarded as a useful dataset for the benchmark of computational models aiming to produce the transient simulation of the aeolian erosion of stockpiles. The CFD results are comprised of the modeling of several experimental scenarios. The computational results for the surface wind velocity show a good correlation with the initial deformation of the pile. Based on the results, the isocontours of (u _s /u _r ) presented might be regarded as a good basis for the estimation of the pile shear velocity.