Impact of the Diurnal Cycle of the Atmospheric Boundary Layer on Wind-Turbine Wakes: A Numerical Modelling Study

The wake characteristics of a wind turbine for different regimes occurring throughout the diurnal cycle are investigated systematically by means of large-eddy simulation. Idealized diurnal cycle simulations of the atmospheric boundary layer are performed with the geophysical flow solver EULAG over both homogeneous and heterogeneous terrain. Under homogeneous conditions, the diurnal cycle significantly affects the low-level wind shear and atmospheric turbulence. A strong vertical wind shear and veering with height occur in the nocturnal stable boundary layer and in the morning boundary layer, whereas atmospheric turbulence is much larger in the convective boundary layer and in the evening boundary layer. The increased shear under heterogeneous conditions changes these wind characteristics, counteracting the formation of the night-time Ekman spiral. The convective, stable, evening, and morning regimes of the atmospheric boundary layer over a homogeneous surface as well as the convective and stable regimes over a heterogeneous surface are used to study the flow in a wind-turbine wake. Synchronized turbulent inflow data from the idealized atmospheric boundary-layer simulations with periodic horizontal boundary conditions are applied to the wind-turbine simulations with open streamwise boundary conditions. The resulting wake is strongly influenced by the stability of the atmosphere. In both cases, the flow in the wake recovers more rapidly under convective conditions during the day than under stable conditions at night. The simulated wakes produced for the night-time situation completely differ between heterogeneous and homogeneous surface conditions. The wake characteristics of the transitional periods are influenced by the flow regime prior to the transition. Furthermore, there are different wake deflections over the height of the rotor, which reflect the incoming wind direction.     

The Effects of Building Representation and Clustering in Large-Eddy Simulations of Flows in Urban Canopies

We perform large-eddy simulations of neutral atmospheric boundary-layer flow over a cluster of buildings surrounded by relatively flat terrain. The first investigated question is the effect of the level of building detail that can be included in the numerical model, a topic not yet addressed by any previous study. The simplest representation is found to give similar results to more refined representations for the mean flow, but not for turbulence. The wind direction on the other hand is found to be important for both mean and turbulent parameters. As many suburban areas are characterised by the clustering of buildings and homes into small areas separated by surfaces of lower roughness, we look at the adjustment of the atmospheric surface layer as it flows from the smoother terrain to the built-up area. This transition has unexpected impacts on the flow; mainly, a zone of global backscatter (energy transfer from the turbulent eddies to the mean flow) is found at the upstream edge of the built-up area.     

Optimization of Upstream Profiles in Modelled Flow Over Complex Terrain

This paper describes an adjoint method for data assimilation intoupstream boundary conditions of numerical modelsusing optimal control theory. Mathematical formalisms are given along with the numerical implementation of the schemein a column model of the atmospheric boundary layer. The optimized mean and turbulence profiles are used as an upstream solutionin a model of turbulent flow in complex terrain. To contrast thiswith other methods, two solutions for flow over an isolatedhill are calculated, one with an optimized upstream solution andone with a simple surface-layer formulation for the upstream solution.These two solutions are compared to observations and analytical theory. The adjoint optimization method is shown to producesolutions of flow in complex terrain that are substantively differentat the two solutions, with the optimized solution giving more accurate results.     

Pollutant Plume Dispersion in the Atmospheric Boundary Layer over Idealized Urban Roughness

The Gaussian model of plume dispersion is commonly used for pollutant concentration estimates. However, its major parameters, dispersion coefficients, barely account for terrain configuration and surface roughness. Large-scale roughness elements (e.g. buildings in urban areas) can substantially modify the ground features together with the pollutant transport in the atmospheric boundary layer over urban roughness (also known as the urban boundary layer, UBL). This study is thus conceived to investigate how urban roughness affects the flow structure and vertical dispersion coefficient in the UBL. Large-eddy simulation (LES) is carried out to examine the plume dispersion from a ground-level pollutant (area) source over idealized street canyons for cross flows in neutral stratification. A range of building-height-to-street-width (aspect) ratios, covering the regimes of skimming flow, wake interference, and isolated roughness, is employed to control the surface roughness. Apart from the widely used aerodynamic resistance or roughness function, the friction factor is another suitable parameter that measures the drag imposed by urban roughness quantitatively. Previous results from laboratory experiments and mathematical modelling also support the aforementioned approach for both two- and three-dimensional roughness elements. Comparing the UBL plume behaviour, the LES results show that the pollutant dispersion strongly depends on the friction factor. Empirical studies reveal that the vertical dispersion coefficient increases with increasing friction factor in the skimming flow regime (lower resistance) but is more uniform in the regimes of wake interference and isolated roughness (higher resistance). Hence, it is proposed that the friction factor and flow regimes could be adopted concurrently for pollutant concentration estimate in the UBL over urban street canyons of different roughness.     

Large-Eddy Simulations of Atmospheric Flows Over Complex Terrain Using the Immersed-Boundary Method in the Weather Research and Forecasting Model

Atmospheric flow over complex terrain, particularly recirculation flows, greatly influences wind-turbine siting, forest-fire behaviour, and trace-gas and pollutant dispersion. However, there is a large uncertainty in the simulation of flow over complex topography, which is attributable to the type of turbulence model, the subgrid-scale (SGS) turbulence parametrization, terrain-following coordinates, and numerical errors in finite-difference methods. Here, we upgrade the large-eddy simulation module within the Weather Research and Forecasting model by incorporating the immersed-boundary method into the module to improve simulations of the flow and recirculation over complex terrain. Simulations over the Bolund Hill indicate improved mean absolute speed-up errors with respect to previous studies, as well an improved simulation of the recirculation zone behind the escarpment of the hill. With regard to the SGS parametrization, the Lagrangian-averaged scale-dependent Smagorinsky model performs better than the classic Smagorinsky model in reproducing both velocity and turbulent kinetic energy. A finer grid resolution also improves the strength of the recirculation in flow simulations, with a higher horizontal grid resolution improving simulations just behind the escarpment, and a higher vertical grid resolution improving results on the lee side of the hill. Our modelling approach has broad applications for the simulation of atmospheric flows over complex topography.     

The dispersion of chemically reactive species in the atmospheric boundary layer

Summary The role of turbulence in the dispersion of atmospheric pollutants that react with linear (decay) and nonlinear (second-order) chemical reactions is examined. The most relevant processes that drive the reactivity of species emitted in a surface area or released by a point source are studied by deriving the dimensionless scaling numbers from equations for the atmospheric turbulent reacting flow. The first number is the ratio of the time scale of turbulence to the time scale of the chemical reaction, namely the Damköhler number. The second number is the ratio of the concentrations of the species present in the chemical transformations. In this paper, model results and experimental studies of turbulent reacting flows in the atmospheric boundary layer are presented to show the modifications and control exerted by turbulence on the atmospheric chemistry as a function of these numbers and processes. We also discuss how the chemical transformation is affected when species are in a state of chemical equilibrium.By studying the plume dispersion of a reactant, that decays with a simple chemical reaction, one can analyse the dependence of concentration fluctuations on the Damköhler number. The study is extended to plumes that react nonlinearly. In such reacting systems, the large gradients and segregation of species result in a significant reduction in the reaction rates. Because of this modification, the chemistry of species related to NO_x and HO_x can be very different from the chemistry in conditions where the species are uniformly mixed. The lack of complete observational evidence is hampering our understanding of these processes and our evaluation of numerical modelling results. Finally, we discuss briefly how to represent, in the form of a parameterization, the effect that turbulence can have on the reactivity of species emitted by a point source or an area source.     

A 3-dimensional nonhydrostatic dispersion modeling system for modeling of atmospheric transport and diffusion over coastal complex terrain in the Hongkong-Shenzhen area

Summary A dispersion modeling system consisting of a three-dimensional PBL model NHECM (non-hydrostaticE- closure model) and SLPTDM (seven-level puff transport and diffusion model) is developed to simulate the transport and dispersion of pollutant over coastal complex terrain. As an application of the system, the transport and dispersion of SO₂ released from an elevated point source are simulated during typical sea-land breeze circulation in the Hongkong-Shenzhen area of China. The NHECM provides time-varying, three-dimensional distributions of wind and turbulence fields to the SLPTDM. The NHECM predictions compare well with observational data. Reflection of both the ground and the mixing layer top and penetration of the mixing layer top are improved in the SLPTDM. Results obtained from the system indicate that temporal variation and nonuniformity of airflow and turbulence obviously affect the concentration distributions, especially during the sea-land breeze transition period. A diurnal cycle of the GLC (ground-level concentration) is discussed. The results are compared with those estimated using a Gaussian model. The study's results illustrate the complexity of the dispersion patterns due to diurnal effects and mesoscale circulations, and demonstrate the potential of the mesoscale atmospheric dispersion modeling system for studies of air quality in complex terrain.With 8 Figures     

Test Of Turbulence Models For Wind Flow Over Terrain With Separation And Recirculation

Several two-equation turbulence models using isotropic eddy viscosity and wall functions are assessed by solution of the neutral atmospheric boundary layer over a flat surface and wind flow over two- and three-dimensional models and real terrain. Calculations are presented for wind flow over the Sirhowy Valley in Wales, an embankment along the Rhine in Germany and the Askervein Hill in Scotland. Comparisons of predictions with previous work, and laboratory and field data, show that the RNG-based k– model gives the best agreement with respect to the flow profiles and length of the separated flow region. The results of this model are analyzed with a non-linear stress-strain relation to gauge the potential effect of turbulence anisotropy.     

Structure of mean winds and turbulence in the planetary boundary layer over rural terrain

A detailed analysis has been carried out of the temporal and spatial structure of mean winds and turbulence in the neutrally-stable planetary boundary layer over typically rural terrain. The data were obtained from a horizontal array of tower-mounted propeller anemometers (z = 11 m) during a five-hour period for which the mean wind direction was virtually perpendicular to the main span of the array. Various turbulence characteristics have been obtained for all three components of velocity and have been compared with idealized models for such a flow and with some of the other available atmospheric results.Considerable tower-to-tower and block-to-block variability has been observed in many of the measured results, particularly in those for the horizontal-component integral scales. Surface shear stress, roughness length and turbulence intensities were in good agreement with expected values for such a site. Power spectra for all components displayed significantly more energy at middle and lower frequencies than that observed by Kaimal et al. (1972) over flat, relatively featureless terrain. This is felt to be a result of the generally rougher gross features of the terrain in the present case and has led to the development of a modified version of the Kaimal-spectral model which fits the observed data better than either the original Kaimal model or the von Kármán model. It is suggested that it may in future be possible to represent power spectra over a wide range of terrain types by using such a modified spectral model.Integral scales of turbulence were calculated by three different techniques and in most cases displayed a strong dependence on the technique used. Averaged values of scale showed reasonable agreement with most of the available atmospheric data and with the values suggested by ESDU (1975). The anticipated elongation of turbulent eddies in the longitudinal direction was confirmed for all three velocity components, although it was found to be not as large as some other observations.     

Operating ranges of mesoscale numerical models and meteorological wind tunnels for the simulation of sea and land breezes

The operating ranges of mesoscale numerical models and meteorological wind tunnels for sea- and land-breeze simulations are defined in this paper based on a review of the theoretical and practical limitations of these two approaches. Numerical-model operating ranges are limited by the choice of governing equations, the numerical methods used to solve the governing equations, the scales of the surface or atmospheric forcing and the atmospheric response, the specified grid resolution and domain size, and the available computer resources. Wind-tunnel operating ranges are limited by the dimensions of the simulated circulations and of the tunnel itself, the tunnel flow speed and turbulence characteristics, the temperature gradients within the tunnel, the lack of Coriolis force and moist processes, and the characteristics of the measurement instrumentation. The operating ranges of these two simulation methods are shown to overlap. In this common range, results of simulations from both approaches can be compared so as to strengthen the validity of the results and to help in the development and improvement of parameterizations of physical processes in numerical models. In addition, the coupling of meteorological wind tunnels and mesoscale numerical models offers a larger range of operating conditions than can be achieved by either approach alone. Together, they can be used in a hybrid form to predict atmospheric conditions at the scale of a few meters for complex terrain (e.g., buildings, hills, etc.) within larger mesoscale atmospheric flow regimes. In the case of sea and land breezes, the sea-land transition zone and coastal internal boundary layer can be studied using both approaches.     

Turbulence Intensity Parameters over a Very Complex Terrain

Detailed knowledge of turbulence structure is important for the understanding of atmospheric phenomena in the boundary layer, especially over complex terrain. In the present study, turbulence intensity parameters are analyzed for different conditions regarding stability, wind speed and wind direction over a mountainous region. The purpose of the analysis is to verify whether the observed parameters follow Monin–Obukhov similarity theory (MOST), despite the terrain heterogeneity. The dataset was collected during an experimental campaign at the Nova Roma do Sul site, in southern Brazil, with a micrometeorological tower located near a sharp slope, approximately 400 m high. The results show that the normalized standard deviations of the vertical velocity component as well as the normalized standard deviation of temperature follow Monin–Obukhov similarity for all stability regimes, regardless of the wind direction. However the normalized standard deviation of the horizontal components of the turbulent velocity obeys the similarity relationship only for a limited range of the stability parameters.     

Resolved-scale turbulence in the atmospheric surface layer from a large eddy simulation

Large eddy simulation has encountered difficulties in handling turbulence in the atmospheric surface layer due to deficiencies in sub-grid scale models. This paper addresses the possibility of resolving the turbulence in the upper part of the surface layer by a low-aspect ratio of grid spacing. Results show that resolved-scale shear stresses dominate over the sub-grid scale components so that effects due to the sub-grid scale model can be ignored in this region. The effects of the lower boundary condition on the resolved-scale turbulence in the upper part of the surface layer are discussed. It is concluded that the normalized mean velocity shear and resolved turbulence in the upper part of the surface layer are not affected by the specification of the lower boundary condition. In addition, the present work proposes a new independent model parameter, the Smagorinsky Reynolds Number (Re_SM), and demonstrates that this number determines the resolved turbulence in the upper part of the surface layer.     

Mean horizontal wind profiles measured in the atmospheric boundary layer about a simulated block building

Instrumented wind towers are used to measure the three components of wind about a simulated block building. The mean horizontal wind profiles over the building are compared with wind profiles measured in the absence of the building and the wind speed deficit in the wake of the building is correlated.Horizontal mean wind speeds measured in the natural atmospheric boundary layer with and without the presence of a simulated building show excellent reproducibility and agreement with fundamental concepts of fluid mechanics. The data possess all the characteristic features reported from wind-tunnel studies of building flows. In the present study the turbulence intensity is of the order of 20% in the undisturbed flow whereas the free stream turbulence intensity of wind-tunnel studies is generally not more than 5%. The effect of smaller averaging periods and the structure of the turbulence will be reported at a later time.The velocity profiles measured in the undisturbed flow zones support the representation of a neutrally stable atmospheric boundary layer with a logarithmic wind profile.     

Characteristics of canopy turbulence over a deciduous forest on complex terrain

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy.     

对流边界层中过山气流的数值模拟

采用ARPS4.0非静力中尺度气象模式模拟了对流边界层中气流过山引起的地形波,讨论了地形及大气条件改变对其的影响.模拟表明,当大气边界层是对流边界层时,气流过山引起的地形强迫,仍能在上部稳定层结中造成足够的垂直扰动,产生向上传播的重力内波,重力内波引起的波动阻力仍不可忽略.     

广州市近地层大气的湍流微结构和谱特征

本文利用UVW脉动风速仪资料分析了广州市区近地层大气的湍流强度、相关系数、尺度和速度谱,并获得了不同稳定度条件下的速度谱模式.结果表明,城市近地层大气湍流在惯性副区接近局地各向同性、速度谱符合Kolmogorov相似理论;气流方向上下垫面粗糙度的增加,使沿海城市近地层大气湍流能量(特别是铅直方向)比平坦、均匀下垫面上的增加.     

Turbulence intensity and turbulent kinetic energy parameters over a heterogeneous terrain of Loess Plateau

A deep understanding of turbulence structure is important for investigating the characteristics of the atmospheric boundary layer, especially over heterogeneous terrain. In the present study, turbulence intensity and turbulent kinetic energy (TKE) parameters are analyzed for different conditions with respect to stability, wind direction and wind speed over a valley region of the Loess Plateau of China during December 2003 and January 2004. The purpose of the study is to examine whether the observed turbulence intensity and TKE parameters satisfy Monin-Obukhov similarity theory (MOST), and analyze the wind shear effect on, and thermal buoyancy function of, the TKE, despite the terrain heterogeneity. The results demonstrate that the normalized intensity of turbulence follows MOST for all stability in the horizontal and vertical directions, as well as the normalized TKE in the horizontal direction. The shear effect of the wind speed in the Loess Plateau region is strong in winter and could enhance turbulence for all stability conditions. During daytime, the buoyancy and shear effect together constitute the generation of TKE under unstable conditions. At night, the contribution of buoyancy to TKE is relatively small, and mechanical shearing is the main production form of turbulence.     

Effects of Topographical Slope Angle and Atmospheric Stratification on Surface-Layer Turbulence

The effect of topographical slope angle and atmospheric stratification on turbulence intensities in the unstably stratified surface layer have been parameterized using observations obtained from a three-dimensional sonic anemometer installed at 8 m height above the ground at the Seoul National University (SNU) campus site in Korea for the years 1999–2001. Winds obtained from the sonic anemometer are analyzed according to the mean wind direction, since the topographical slope angle changes significantly along the azimuthal direction. The effects of the topographical slope angle and atmospheric stratification on surface-layer turbulence intensity are examined with these data. It is found that both the friction velocity and the variance for each component of wind normalized by the mean wind speed decrease with increase of the topographical slope angle, having a maximum decreasing rate at very unstable stratification. The decreasing rate of the normalized friction velocity (u _* /U) is found to be much larger than that of the turbulence intensity of each wind component due to the reduction of wind shear with increase in slope angle under unstable stratification. The decreasing rate of the w component of turbulence intensity (σ _w /U) is the smallest over the downslope surface whereas that of the u component (σ _u /U) has a minimum over the upslope surface. Consequently, σ _w /u _* has a maximum increasing rate with increase in slope angle for the downslope wind, whereas σ _u /u _* has its maximum for the upslope wind. The sloping terrain is found to reduce both the friction velocity and turbulence intensity compared with those on a flat surface. However, the reduction of the friction velocity over the sloping terrain is larger than that of the turbulence intensity, thereby enhancing the turbulence intensity normalized by the friction velocity over sloping terrain compared with that over a flat surface.     

Etude de la couche limite convective sahelienne en presence de brumes seches (Expérience ECLATS)

The ECLATS experiment was conducted in order to investigate the influence of radiative processes on the dynamics of the atmospheric boundary layer during its diurnal evolution. This experiment was carried out over Niger, near Niamey, by measuring continuously the energy balance at ground level and by using an instrumented aircraft for turbulence, radiative fluxes and aerosol measurements in the boundary layer during dusty conditions (brumes sèches). This paper is restricted to an analysis of the turbulent structure in the homogeneous and stationary convective boundary layer. The turbulence moments for kinetic energy and the spectral characteristics of the vertical velocity are discussed. These results are compared with a set of data obtained for clear convective boundary layers. The differences observed are quite important and seem, at least in part, due to radiative processes (infrared radiative divergence in the surface layer and absorption of solar radiation in the boundary layer).     

Comparing mixing-length models of the diabatic wind profile over homogeneous terrain

Models of the diabatic wind profile over homogeneous terrain for the entire atmospheric boundary layer are developed using mixing-length theory and are compared to wind speed observations up to 300 m at the National Test Station for Wind Turbines at Høvsøre, Denmark. The measurements are performed within a wide range of atmospheric stability conditions, which allows a comparison of the models with the average wind profile computed in seven stability classes, showing a better agreement than compared to the traditional surface-layer wind profile. The wind profile is measured by combining cup anemometer and lidar observations, showing good agreement at the overlapping heights. The height of the boundary layer, a parameter required for the wind profile models, is estimated under neutral and stable conditions using surface-layer turbulence measurements, and under unstable conditions based on the aerosol backscatter profile from ceilometer observations.