Physical modeling of complex airflows developing above rural terrains

The rural atmospheric boundary layer (ABL) flow was reproduced in a wind tunnel at three different simulation length scales to investigate possible effects of the simulation length scale on flow characteristics. Performance of truncated vortex generators developed for part-depth ABL wind-tunnel simulations was tested in rural terrain exposure against the full-size Counihan vortex generators. A procedure to design the ABL developing above rural type terrain has been described. The 1:395 and 1:236 simulations were created as full-depth simulations, i.e., wind characteristics throughout an entire ABL were reproduced in the wind tunnel. The 1:208 simulation was a part-depth simulation, i.e., only a lower 70% of the ABL was experimentally modelled. The projected scaled-up ABL thicknesses are 395, 354, and 416 m full-scale in the 1:395, 1:236, and 1:208 simulations, respectively. Experimental results show similar trends in all three configurations not depending on the simulation length scale factor. This clearly indicates a possibility to physically, in the wind tunnel, reproduce the same rural atmospheric airflows at different simulation length scales.     

Wind-induced pressure at a tunnel portal

In order to properly size the mechanical ventilation system of a tunnel, it is essential to estimate the wind-driven pressure difference that might rise between its two portals. In this respect, we explore here the pressure distribution over a tunnel portal under the influence of an incident atmospheric boundary layer and, in particular, its dependency on wind direction and on tunnel geometry. Reduced scale models of generic configurations of a tunnel portal are studied in an atmospheric wind tunnel. Pressure distributions over the front section of different open cavities are measured with surface taps, which allows us to infer the influence of the tunnel aspect ratio and wind direction on a pressure coefficient \(C_{P}\), defined as a spatially and time averaged non-dimensional pressure. Experiments reveal that the magnitude of the coefficient \(C_{P}\), as a function of the wind direction, is significantly influenced by the portal height-to-width ratio and almost insensitive to its length. The experimental data set is completed by hot-wire anemometry measurements providing vertical distribution of velocity statistics. The same configurations are simulated by numerically solving the Reynolds-averaged Navier–Stokes equations, adopting the standard \(k - \varepsilon\) turbulence model. Despite some discrepancies between numerical and experimental estimates of some flow parameters (namely the turbulent kinetic energy field), the numerical estimates of the pressure coefficients \(C_{P}\) show very good agreement with experimental data. The latter is also compared to the predictions of an analytical model, based on the estimate of a spatially averaged velocity within an infinitely long street canyon. The results of the model, which takes into account varying canyon aspect ratios, are in reasonable agreement with experimental data for all cases studied. Notably, its predictions are significantly better than those provided by the simple analytical relations usually adopted as a reference in tunnel ventilation studies.     

Evaluation of CFD Simulation using RANS Turbulence Models for Building Effects on Pollutant Dispersion

CFD evaluations were performed to examine the applicability of the RANS methods in simulating pollutant dispersion near, within and over three typical building configurations: (1) an isolated building, (2) a building array and (3) an urban intersection. The CFD results are compared with values obtained from wind tunnel tests. In some situations major differences between the wind tunnel tests and the CFD results were observed. The main source of difference between the CFD and wind tunnel results was inadequate modelling of local flow patterns using the RANS turbulence models. Also inappropriate evaluation of high intermittent turbulent mixing in the RANS approach may lead to either over-prediction or under-prediction of the concentration level, by up to a factor of 10, depending on the case investigated.     

Experimental surface flow visualization and numerical investigation of flow structure around an oblong stockpile

Emission factors are largely used to quantify particle emissions from industrial open storage piles. These factors are based on the knowledge of velocity distribution and flow patterns over the stockpile surface which still requires further research. The aim of the present work is to investigate the airflow characteristics over a single typical oblong pile and in its near-ground surroundings for various wind flow directions. Wind tunnel experiments using an oil-film surface coating technique were carried out for near-wall flow visualization. Numerical simulation results, favorably compared to PIV measurements, were used to allow comparison analysis of flow features. For the stockpile oriented 90° to the wind main direction, typical topology of flow around wall-mounted obstacles were observed, notably a wake zone downstream the pile including two main counter-rotating vortices. Further analysis of numerical wall shear stress distribution and streamlines indicates that two complex three-dimensional vortical flow structures develop downstream the pile. For other incoming wind flow directions (30 and 60°), the flow characteristics over the storage pile greatly differ as a single helical main vortex develops from the pile’s crest. Corresponding high values of wall shear stress are noticed downstream the storage pile. For each configuration studied, downwash and upwash zones are induced by the vortical structures developed. This near-wall flow topology combined with areas of high friction levels may be linked to potential dust emission from the ground surface surrounding industrial stockpiles.     

Evaluation of wind environment around a residential complex using a PIV velocity field measurement technique

Wind-tunnel simulations were employed to evaluate the wind environment around a tested residential area located near industrial complexes. The scaled-down geomorphological model of the test area was placed in the test section of a boundary layer wind tunnel. Particle image velocimetry (PIV) measurements were made in five vertical planes and one horizontal plane around the test area for two prevailing wind directions. The results showed that the wind speed decreased in the near surface layer and the velocity fluctuations increased in the upper region due to the presence of hills and high-rise buildings around the test area. Regions of flow separation and low-speed flow were found inside the test area for both the wind directions. The result suggests that the high-rise buildings should be well arranged with respect to the main wind directions to increase the natural ventilation inside the residential complex at the initial design stage.     

Evaluation of the QUIC-URB fast response urban wind model for a cubical building array and wide building street canyon

This paper describes the QUIC-URB fast response urban wind modeling tool and evaluates it against wind tunnel data for a 7 × 11 cubical building array and wide building street canyon. QUIC-URB is based on the Röckle diagnostic wind modeling strategy that rapidly produces spatially resolved wind fields in urban areas and can be used to drive urban dispersion models. Röckle-type models do not solve transport equations for momentum or energy; rather, they rely heavily on empirical parameterizations and mass conservation. In the model-experiment comparisons, we test two empirical building flow parameterizations within the QUIC-URB model: our implementation of the standard Röckle (SR) algorithms and a set of modified Röckle (MR) algorithms. The MR model attempts to build on the strengths of the SR model and introduces additional physically based, but simple parameterizations that significantly improve the results in most regions of the flow for both test cases. The MR model produces vortices in front of buildings, on rooftops and within street canyons that have velocities that compare much more favorably to the experimental results. We expect that these improvements in the wind field will result in improved dispersion calculations in built environments.     

Development and testing of a micro wind tunnel for on-site wind erosion simulations

Wind erosion processes affect soil surfaces across all land uses worldwide. Understanding the spatial and temporal scales of wind erosion is a challenging undertaking because these processes are diverse and highly variable. Wind tunnels provide a useful tool as they can be used to simulate erosion at small spatial scales. Portable wind tunnels are particularly valued because erosion can be simulated on undisturbed soil surfaces in the field. There has been a long history of use of large portable wind tunnels, with consensus that these wind erosion simulation tools can meet real world aerodynamic criteria. However, one consequence of striving to meet aerodynamic reality is that the size of the tunnels has increased, making them logistically difficult to work with in the field and resulting in a tendency to homogenise naturally complex soil surfaces. This homogenisation is at odds with an increasing awareness of the importance that small scale processes have in wind erosion. To address these logistical and surface homogenisation issues we present here the development and testing of a micro wind tunnel (MWT) designed to simulate wind erosion processes at high spatial resolution. The MWT is a duct-type design—0.05 m tall 0.1 m wide and with a 1.0 m working section. The tunnel uses a centrifugal motor to suck air through a flow‐conditioning section, over the working section and then through a sediment collection trap. Simulated wind velocities range from 5 to 18 m s^?1, with high reproducibility. Wind speeds are laterally uniform and values of u _* at the tunnel bed (calculated by measuring the pressure gradients within the MWT) are comparable with those of larger tunnels in which logarithmic profiles can be developed. Saltation sediment can be added. The tunnel can be deployed by a single person and operated on slopes ranging from 0 to 10°. Evidence is presented here that the MWT provides new and useful understanding of the erodibility of rangelands, claypans and ore stockpiles.     

Effect of the topography of an industrial site on dust emissions from open storage yards

The erosion, transport and subsequent deposition of particles by wind has been identified as an important environmental problem worldwide. Investigations assessing fugitive dust emissions represent a relevant subject for various applications such as prevention of desert expansion, farmland erosion, quantification of dust emissions due to wind erosion from mining zones or open storage yards of granular materials. The present study, within this latter framework, aims to improve the methods of quantification of fugitive emissions from granular materials stockpiles on industrial sites. The methodology mostly used by the industrial operators to quantify dust emissions from their plants is based on the wind flow analysis and application of emission factors into such pattern. In this paper, three-dimensional numerical simulations were done to simulate wind flow over a real configuration of a power plant to highlight the influence of the surroundings buildings on the real exposure of the granular materials stockpiles. This study emphasizes the necessity to take into consideration the presence of surrounding buildings for the calculation of the particle emissions rate on stockpiles in order to obtain a more accurate and relevant evaluation of fugitive dust emissions from open storage systems on industrial sites.     

Concentration fluctuations in a downtown urban area. Part II: analysis of Joint Urban 2003 wind-tunnel measurements

An analysis of concentration time series measured in a boundary-layer wind tunnel at the University of Hamburg is presented. The measurements were conducted with a detailed aerodynamic model of the Oklahoma City (OKC) central business district (CBD) at the scale of 1:300 and were part of the Joint Urban 2003 (JU2003) project. Concentration statistics, as well as concentration probability density (PDF) and exceedance probability (EDF) functions were computed for street- and roof-level sites for three different wind directions. Taking into account the different length scales and wind speeds in the wind-tunnel (WT) and full-scale experiments, dimensionless concentrations and a dimensionless time scale are computed for the comparison with data from the JU2003 full-scale tracer experiments, conducted in OKC in 2003. Using such dimensionless time, the WT time series cover a ~20 times longer time span than the JU2003 full-scale time series, which are analysed in detail in an accompanying, first part of this paper. The WT time series are thus divided into 20 consecutive blocks of equal length and the statistical significance of parameters based on relatively short records is assessed by studying the variability of the concentration statistics and probability functions for the different blocks. In particular at sites closer to the plume edge, the results for the individual blocks vary significantly and at such sites statistics from short records are not very representative. While the location of three sampling sites in the WT closely matched the sites during the full-scale experiments, the prevailing wind directions during the JU2003 releases were not exactly matched. The comparison between full-scale and WT concentration parameters should thus primarily be interpreted in a qualitative rather than direct quantitative sense. Given the differences in mean wind directions and concerns about the representativeness of full-scale concentration statistics, the WT and full-scale results compared well. The 98 percentile concentrations for almost all full-scale releases analyzed are within the scatter of the percentiles observed in the block analysis of the WT time series. Furthermore, the concentration percentiles appear linearly correlated with the fluctuation intensities and the linear relationships determined in the wind tunnel agree well with full-scale results.     

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 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.     

Canopy-wake dynamics and wind sheltering effects on Earth surface fluxes

The atmospheric boundary layer adjustment at the abrupt transition from a canopy (forest) to a flat surface (land or water) is investigated in a wind tunnel experiment. Detailed measurements examining the effect of canopy turbulence on flow separation, reduced surface shear stress and wake recovery are compared to data for the classical case of a solid backward-facing step. Results provide new insights into the interpretation for flux estimation by eddy-covariance and flux gradient methods and for the assessment of surface boundary conditions in turbulence models of the atmospheric boundary layer in complex landscapes and over water bodies affected by canopy wakes. The wind tunnel results indicate that the wake of a forest canopy strongly affects surface momentum flux within a distance of 35–100 times the step or canopy height, and mean turbulence quantities require distances of at least 100 times the canopy height to adjust to the new surface. The near-surface mixing length in the wake exhibits characteristic length scales of canopy flows at the canopy edge, of the flow separation in the near wake and adjusts to surface layer scaling in the far wake. Components of the momentum budget are examined individually to determine the impact of the canopy wake. The results demonstrate why a constant flux layer does not form until far downwind in the wake. An empirical model for surface shear stress distribution from a forest canopy to a clearing or lake is proposed.     

冬季降雪过程对城市大气气态汞污染的影响

2009年降雪和非降雪期间对北京西北城区的气态总汞浓度进行了连续采样,比较了降雪期间、非降雪期间的气态总汞浓度日变化过程;降雪期间气态总汞浓度的降低和恢复过程。结果表明,降雪和非降雪期间大气气态总汞浓度的日均值有显著差异,降雪期间气态总汞的平均浓度为5．64ng·m^-3,非降雪期间的平均浓度为7．43ng·m^-3,前者约为后者的70％。降雪后约7h气态总汞浓度恢复到降雪前水平。研究中分析了气象因素（气压、风速、阵风速度、气温和相对湿度）对于气态总汞浓度的影响,结果表明：降雪期间主要受到风速（r=-0．527）和阵风速度（r=-0．574）的影响;非降雪期间主要受到风速（r=-0．691）,阵风速度（r=-0．726）和相对湿度（r=0．692）的影响,并且相对湿度的影响与风速的影响相近。降雪和非降雪期间气态总汞的日变化有所差异：非降雪期间气态总汞浓度在午夜和清晨较高,日变化趋势与相对湿度一致;降雪期间气态总汞的日变化没有明显规律。     

Effects of sand movement by wind on nematodes and soil-borne fungi in coastal foredunes

In stabilized dunesAmmophila arenaria (marram grass) degenerates due to a process involving soil-borne pathogens and parasites. This leads to exposure of the sand surface so that wind erosion may create blowouts.Ammophila rejuvenates on the edges of the blowouts, where the sand has accumulated. We tested the hypothesis that such rejuvenation of plants may be related to a reduction of the plant-parasitic nematodes and fungal propagules during the process of wind-driven transport. Field measurements in blowouts during storm events indicated that the drifted sand contained relatively low numbers of plant pathogenic fungi and plant-parasitic nematodes. A wind tunnel experiment showed that drifting sand may indeed reduce the numbers of fungi and nematodes. Although most fungi were attached to sand particles, they were also affected by the wind-borne sand movement. Sand that had been deposited by wind was made up of a larger proportion of large-sized particles. In our experiment the relatively small particles were lost during transport. Stirring the soil (part of the forces of natural winds) by mixing for 15 min. with a propeller mixer at 1500 rpm significantly reduced the number of nematodes and fungi. Both sand movement in the wind tunnel and intensive stirring of the sand enhanced the growth ofAmmophila test plants in a bioassay. It was concluded that in wind-blown sand the pathogen inoculum is reduced. Therefore, serious consideration should be given to allowing controlled reactivation of blowouts to rejuvenate decliningAmmophila in stabilized foredunes. Implications for dune management are briefly discussed. Nomenclature: van der Meijden (1990) for vascular plants. Nematodes were identified to the genus level according to Bongers (1988). The allocation of nematodes to feeding groups was according to Yeates et al. (1993). Nomenclature of fungi according to Domsch et al. (1980) was used throughout this study except forFusarium. Species of the latter genus were identified according to Nelson et al. (1983).     

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.     

RANS simulation of neutral atmospheric boundary layer flows over complex terrain by proper imposition of boundary conditions and modification on the k-ε model

In this study, a modelling methodology is proposed for RANS simulations of neutral Atmospheric Boundary Layer (ABL) flows on the basis of the standard k-ε model, which allows the adoption of an arbitrary shear stress model. This modelling methodology is first examined in the context of an open flat terrain in an empty domain to ascertain there are no substantial changes in the prescribed profiles. The results show that relatively good homogeneity can be achieved with this modelling methodology for various sets of inflow boundary profiles. In addition, to extend the solutions derived from the standard k-ε model to RNG k-ε model, the RNG k-ε model is in detail assembly and tuned. Finally, the topographic effects on surface wind speeds over a complex terrain are assessed with the combined use of the proposed methodology and the modified RNG model. The numerical results are in good agreement with wind tunnel testing results and long-term field observations. A discussion of the effects of horizontal homogeneity and turbulence models on the simulated wind flows over a complex terrain is also given.     

Nest selection by snow petrels Pagodroma nivea in East Antarctica: Validating predictive habitat selection models at the continental scale

Little is known on the factors controlling distribution and abundance of snow petrels in Antarctica. Studying habitat selection through modeling may provide useful information on the relationships between this species and its environment, especially relevant in a climate change context, where habitat availability may change. Validating the predictive capability of habitat selection models with independent data is a vital step in assessing the performance of such models and their potential for predicting species’ distribution in poorly documented areas.From the results of ground surveys conducted in the Casey region (2002–2003, Wilkes Land, East Antarctica), habitat selection models based on a dataset of 4000 nests were created to predict the nesting distribution of snow petrels as a function of topography and substrate. In this study, the Casey models were tested at Mawson, 3800 km away from Casey. The location and characteristics of approximately 7700 snow petrel nests were collected during ground surveys (Summer 2004–2005). Using GIS, predictive maps of nest distribution were produced for the Mawson region with the models derived from the Casey datasets and predictions were compared to the observed data. Models performance was assessed using classification matrixes and Receiver operating characteristic (ROC) curves. Overall correct classification rates for the Casey models varied from 57% to 90%. However, two geomorphologically different sub-regions (coastal islands and inland mountains) were clearly distinguished in terms of habitat selection by Casey model predictions but also by the specific variations in coefficients of terms in new models, derived from the Mawson data sets. Observed variations in the snow petrel aggregations were found to be related to local habitat availability.We discuss the applicability of various types of models (GLM, CT) and investigate the effect of scale on the prediction of snow petrel habitats. While the Casey models created with data collected at the nest scale did not perform well at Mawson due to regional variations in nest micro-characteristics, the predictive performance of models created with data compiled at a coarser scale (habitat units) was satisfactory. Substrate type was the most robust predictor of nest presence between Casey and Mawson. This study demonstrate that it is possible to predict at the large scale the presence of snow petrel nests based on simple predictors such as topography and substrate, which can be obtained from aerial photography. Such methodologies have valuable applications in the management and conservation of this top predator and associated resources and may be applied to other Antarctic, Sub-Antarctic and lower latitudes species and in a variety of habitats.     

Effect of pollinator abundance on self-fertilization and gene flow: application to GM Canola.

Cross-pollination from fields of transgenic crops is of great public concern. Although cross-pollination in commercial canola (Brassica napus) fields has been empirically measured, field trials are expensive and do not identify the causes of cross-pollination. Therefore, theoretical models can be valuable because they can provide estimates of cross-pollination at any given site and time. We present a general analytical model of field-to-field gene flow due to the following competing mechanisms: the wind, bees, and autonomous pollination. We parameterize the model for the particular case of field-to-field cross-pollination of genetically modified (GM) canola via the wind and via bumble bees (Bombus spp.) and honey bees (Apis mellifera). We make extensive use of the large data set of bee densities collected during the recent U.K. Farm Scale Evaluations. We predict that canola approaches almost full seed set without pollinators and that autonomous pollination is responsible for > or = 25% of seed set, irrespective of pollinator abundance. We do not predict the relative contribution of bees vs. the wind in landscape-scale gene flow in canola. However, under model assumptions, we predict that the maximum field-to-field gene flow due to bumble bees is 0.04% and 0.13% below the current EU limit for adventitious GM presence for winter- and spring-sown canola, respectively. We predict that gene flow due to bees is approximately 3.1 times higher at 20% compared to 100% male-fertility, and due to the wind, 1.3 times higher at 20% compared to 100% male-fertility, for both winter- and spring-sown canola. Bumble bee-mediated gene flow is approximately 2.7 times higher and wind-mediated gene flow approximately 1.7 times lower in spring-sown than in winter-sown canola, regardless of the degree of male-sterility. The model of cross-pollination due to the wind most closely predicted three previously published observations: field-to-field gene flow is low; gene flow increases with the proportion of plants that are male-sterile; and gene flow is higher in winter- than in spring-sown canola. Our results therefore suggest that the wind, not bees, is the main vector of long-distance gene flow in canola.     

Numerical simulation of sand dune erosion

Erosion of sand or other granular material is a subject of utmost importance in several fields of practical interest, including industrial processes or environmental issues. Resulting from intricate interaction between the incident flow field and localized body forces responsible for the granular material cohesion, erosion is a particularly complex phenomenon. The present work addresses this problem, proposing a numerical method to compute the time evolution of a sand dune subjected to aeolian erosion, along with the associated entrainment and deposition fluxes. Turbulent fluid flow is computed through a three-dimensional Navier-Stokes solver based on a generalized coordinate system. A Lagrangian approach is adopted for tracking the trajectories of particles entrained in the saltation regime, thus allowing prediction of the corresponding deposition locations. Different models for saltation fluxes are tested, along with several formulations for the creeping-to-saltation flux ratio, creeping threshold and creeping distance. Comparison with results from wind tunnel experiments is very encouraging, stressing the relative importance of creeping in the erosion process for the presently studied conditions.