Flow Over Partially Forested Ridges

Numerical simulations of flow over hills that are partially covered with a forest canopy are performed. This represents a much more realistic situation than previous studies that have generally concentrated on hills that are fully-forested. The results show that the flow over the hill is sensitive to where on the hill the forest is positioned. In particular, for low slopes flow separation is predominantly located within the forest on the lee slope. This has implications for the transport of scalars in the forest canopy. For large hills the results show more variability in scalar concentrations within the canopy compared to either a fully-forested hill or a patch of forest over flat terrain. These results are likely to have implications for a range of applications including the siting and interpretation of flux measurements over forests in complex terrain, predicting wind damage to trees and wind-farm developments. Calculation of the hill-induced pressure drag and canopy-plus-surface stress shows a strong sensitivity to the position of the forest relative to the hill. Depending on the position of the forest the individual drag terms may be strongly enhanced or reduced and may even change sign. The net impact is generally to reduce the total drag compared to an equivalent fully-forested hill, but the amount of the reduction depends strongly on the position of the forest canopy on the hill. In many cases with large, wide hills there is a clear separation of scales between the adjustment of the canopy to a forest edge (of order 6 ? 8L _c, where L _c is the canopy adjustment length scale) and the width of the hill. This separation means that the hill-induced pressure and flow fields and the forest-edge induced pressure and flow fields can in some sense be considered as acting separately. This provides a means of explaining the combined effects of partial forestation and terrain. It also offers a simple method for modelling the changes in drag over a hill due to partial forest cover by considering the impact of the hill and the partial canopy separately. Scaling arguments based on this idea successfully collapse the modelled drag over a range of different hill widths and heights and for different canopy parameters. This offers scope for a relatively simple parametrization of the effects of partial forest cover on the drag over a hill.     

Large-eddy Simulations of Flow Over Forested Ridges

Large-eddy simulations (LES) of flow over a series of small forested ridges are performed, and compared with numerical simulations using a one-and-a-half order mixing length closure scheme. The qualitative and quantitative similarity between these results provides some confidence in the results of recent analytical and numerical studies of flow over forested hills using first-order mixing length schemes. Time series of model velocities at various locations within the canopy allow the application of various experimental techniques to study the turbulence in the LES. The application of conditional analysis shows that the structure of the turbulence over a forested hill is broadly similar to that over flat ground, with sweeps and ejections dominating. Differences are seen across the hill, particularly associated with regions of mean flow separation and recirculation near the summit and in the lee of the hill. Detailed comparison of derived mixing lengths from the LES with the assumed values used in mixing-length closure schemes show that the mixing length varies with location across the hill and with height in the canopy. This is consistent with previous wind-tunnel measurements, and demonstrates that a constant mixing-length assumption is not strictly valid within the canopy. Despite this, the first-order mixing-length schemes do give similar results both for the mean flow and the turbulence in such situations.     

Large-Eddy Simulation Of Neutral Turbulent Flow Over Rough Sinusoidal Ridges

The requirements for a credible large-eddy simulation of neutral, turbulent flow over hills with an aerodynamically rough surface are discussed, in order to select a suitable case for simulation. As well as providing adequate resolution within the dynamically important inner region, obtaining a realistic upstream or undisturbed mean velocity profile is also of critical importance. A distributed drag canopy formulation has been introduced to the model to allow it to obtain such a profile while resolving very close to the rough surface. Simulations have then been performed of flow over ridges of varying heights. The results from the steepest case, which is just on the verge of separation, are compared with wind-tunnel observations. It is shown that the large-eddy simulation results are in much better agreement with the experimental data than are the results from a simple first-order mixing-length closure model. An encouraging lack of sensitivity of the simulation results to numerical resolution is also demonstrated.     

Turbulent Pressure and Velocity Perturbations Induced by Gentle Hills Covered with Sparse and Dense Canopies

How the spatial perturbations of the first and second moments of the velocity and pressure fields differ for flow over a train of gentle hills covered by either sparse or dense vegetation is explored using large-eddy simulation (LES). Two simulations are investigated where the canopy is composed of uniformly arrayed rods each with a height that is comparable to the hill height. In the first simulation, the rod density is chosen so that much of the momentum is absorbed within the canopy volume yet the canopy is not dense enough to induce separation on the lee side of the hill. In the second simulation, the rod density is large enough to induce recirculation inside the canopy on the lee side of the hill. For this separating flow case, zones of intense shear stress originating near the canopy-atmosphere interface persist all the way up to the middle layer, ‘contaminating’ much of the middle and outer layers with shear stress gradients. The implications of these persistent shear-stress gradients on rapid distortion theory and phase relationships between higher order velocity statistics and hill-induced mean velocity perturbations (Δu) are discussed. Within the inner layer, these intense shear zones improve predictions of the spatial perturbation by K-theory, especially for the phase relationships between the shear stress (~ ?Δu/?z) and the velocity variances, where z is the height. For the upper canopy layers, wake production increases with increasing leaf area density resulting in a vertical velocity variance more in phase with Δu than with ?Δu/?z. However, background turbulence and inactive eddies may have dampened this effect for the longitudinal velocity variance. The increase in leaf area density does not significantly affect the phase relationship between mean surface pressure and topography for the two simulations, though the LES results here confirm earlier findings that the minimum mean pressure shifts downstream from the hill crest. The increase in leaf area density and associated flow separation simply stretches this difference further downstream. This shift increases the pressure drag, the dominant term in the overall drag on the hill surface, by some 15%. With regards to the normalized pressure variance, increasing leaf area density increases ${\sigma_p/u_{*}^{2}}$ near the canopy top, where u _* is the longitudinally averaged friction velocity at the canopy top and σ _p is the standard deviation of the pressure fluctuations. This increase is shown to be consistent with a primitive scaling argument on the leading term describing the mean-flow turbulent interaction. This scaling argument also predicts the spatial variations in σ _p above the canopy reasonably well for both simulations, but not inside the canopy.     

Comparison of Speed-Up Over Hills Derived from Wind-Tunnel Experiments,Wind-Loading Standards,and Numerical Modelling

We evaluate the accuracy of the speed-up provided in several wind-loading standards by comparison with wind-tunnel measurements and numerical predictions, which are carried out at a nominal scale of 1:500 and full-scale, respectively. Airflow over two- and three-dimensional bell-shaped hills is numerically modelled using the Reynolds-averaged Navier–Stokes method with a pressure-driven atmospheric boundary layer and three different turbulence models. Investigated in detail are the effects of grid size on the speed-up and flow separation, as well as the resulting uncertainties in the numerical simulations. Good agreement is obtained between the numerical prediction of speed-up, as well as the wake region size and location, with that according to large-eddy simulations and the wind-tunnel results. The numerical results demonstrate the ability to predict the airflow over a hill with good accuracy with considerably less computational time than for large-eddy simulation. Numerical simulations for a three-dimensional hill show that the speed-up and the wake region decrease significantly when compared with the flow over two-dimensional hills due to the secondary flow around three-dimensional hills. Different hill slopes and shapes are simulated numerically to investigate the effect of hill profile on the speed-up. In comparison with more peaked hill crests, flat-topped hills have a lower speed-up at the crest up to heights of about half the hill height, for which none of the standards gives entirely satisfactory values of speed-up. Overall, the latest versions of the National Building Code of Canada and the Australian and New Zealand Standard give the best predictions of wind speed over isolated hills.     

Dispersion of a Passive Scalar Within and Above an Urban Street Network

The transport of a passive scalar from a continuous point-source release in an urban street network is studied using direct numerical simulation (DNS). Dispersion through the network is characterized by evaluating horizontal fluxes of scalar within and above the urban canopy and vertical exchange fluxes through the canopy top. The relative magnitude and balance of these fluxes are used to distinguish three different regions relative to the source location: a near-field region, a transition region and a far-field region. The partitioning of each of these fluxes into mean and turbulent parts is computed. It is shown that within the canopy the horizontal turbulent flux in the street network is small, whereas above the canopy it comprises a significant fraction of the total flux. Vertical fluxes through the canopy top are predominantly turbulent. The mean and turbulent fluxes are respectively parametrized in terms of an advection velocity and a detrainment velocity and the parametrization incorporated into a simple box-network model. The model treats the coupled dispersion problem within and above the street network in a unified way and predictions of mean concentrations compare well with the DNS data. This demonstrates the usefulness of the box-network approach for process studies and interpretation of results from more detailed numerical simulations.     

The influence of geometry on recirculation and CO2 transport over forested hills

Flow distortion over a forested hill is asymmetric, forming a recirculation region on the lee slope that increases the complexity in understanding atmosphere–biosphere interaction. To understand the complexity, we examine the effect of the geometry of forested hills on recirculation formation, structure, and related CO₂ transport by performing numerical simulations over double-forested hills. The ratio (0.8) of hill height (H) to half length (L) is a threshold value of flow patterns in the recirculation region: below 0.8, sporadic reversed flow occurs; at 0.8, one vortex is formed; and above 0.8, a pair of counter-rotating vortices is formed. The depth of recirculation increases with increasing H/L. The contribution of advection to the CO₂ budget is non-negligible and topographic-dependent. Vertical advection is opposite in sign to horizontal advection but cannot exactly offset in magnitude. Height-integrated advection shows significant variation in fluxes across hills. Gentle slopes can cause larger advection error. However, the relative importance of advection to CO₂ budget is slope-independent.     

An Evaluation Of Local Similarity At The Top Of The Mixed Layer Based On Large-Eddy Simulations

Local similarity, referred to as type II similarity,in the interfacial, stably-stratified layer at thetop of the atmospheric (or oceanic) mixed layer isdiscussed. Type II scales for scalars are based onthe local values of scalar gradients. Similaritypredictions are derived from the second-orderclosure model of Yamada and Mellor, and also fromsimilarity arguments. The obtainedformulation is verified for active and passive scalarsbased on the large-eddy simulation model.     

The Effect of Source Distribution on Bulk Scalar Transfer between a Rough Land Surface and the Atmosphere

We investigate the effect of source distribution on the bulk transfer of passive scalars between rough, vegetated land surfaces and the atmosphere, using data from a wind-tunnel experiment in which passive heat was emitted from both the underlying surface and canopy elements of a three-dimensional regular bluff-body array. The experimental results are compared with a simple one-dimensional, two-source model for scalar transfer. We find that: (1) the observed scalar transfer resistance across the boundary layer at the underlying surface is simply related to flat-plate theory by a constant of 0.62, despite the complexity of the turbulent flow within the wind-tunnel canopy; (2) one-dimensional gradient-transfer theory, even with extensions to account for the non-local nature of turbulent transfer within the canopy, does not describe the observed details of scalar concentration gradients in the highly three-dimensional canopy flow, but does provide a reasonable framework for bulk scalar transfer between the composite ground-canopy surface and the flow above the canopy; (3) the kB ⁻¹ parameter (which accounts for bulk excess resistance to scalar transfer over momentum transfer) is highly sensitive to scalar source partition between ground and canopy.     

The Effects of Thermal Stratification on Clustering Properties of Canopy Turbulence

The intermittent structure of turbulence within the canopy sublayer (CSL) is sensitive to the presence of foliage and to the atmospheric stability regime. How much of this intermittency originates from amplitude variability or clustering properties remains a vexing research problem for CSL flows. Using a five-level set of measurements collected within a dense hardwood canopy, the clustering properties of CSL turbulence and their dependence on atmospheric stability are explored using the telegraphic approximation (TA). The binary structure of the TA removes any amplitude variability from turbulent excursions but retains their zero-crossing behaviour, and thereby isolating the role of clustering in intermittency. A relationship between the spectral exponents of the actual and the TA series is derived across a wide range of atmospheric stability regimes and for several flow variables. This relationship is shown to be consistent with a relationship derived for long-memory and monofractal processes such as fractional Brownian motion (fBm). Moreover, it is demonstrated that for the longitudinal and vertical velocity components, the vegetation does not appreciably alter fine-scale clustering but atmospheric stability does. Stable atmospheric stability conditions is characterized by more fine scale clustering when compared to other atmospheric stability regimes. For scalars, fine-scale clustering above the canopy is similar to its velocity counterpart but is significantly increased inside the canopy, especially under stable stratification. Using simplified scaling analysis, it is demonstrated that clustering is much more connected to space than to time within the CSL. When comparing intermittency for flow variables and their TA series, it is shown that for velocity, amplitude variations modulate intermittency for all stability regimes. However, amplitude variations play only a minor role in scalar intermittency. Within the crown region of the canopy, a ‘double regime’ emerges in the inter-pulse duration probability distributions not observed in classical turbulence studies away from boundaries. The double regime is characterized by a power-law distribution for shorter inter-pulse periods and a log-normal distribution for large inter-pulse periods. The co-existence of these two regimes is shown to be consistent with near-field/far-field scaling arguments. In the near-field, short inter-pulse periods are controlled by the source strength, while in the far-field long inter-pulse periods are less affected by the precise source strength details and more affected by the transport properties of the background turbulence.     

Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

Modelling Street-Scale Flow and Dispersion in Realistic Winds—Towards Coupling with Mesoscale Meteorological Models

Further to our previous work—simulations of flow and dispersion in an oblique wind over the DAPPLE site (Xie and Castro, Atmos Environ 43:2174–2185, 2009)—large-eddy simulations of flows and dispersion over the same site in a wind perpendicular to Marylebone Road and the windward surfaces of most of the buildings were performed. The DAPPLE site is located at the intersection of Marylebone Road and Gloucester Place in central London. In order to investigate the effects of wind direction on flows and dispersion, the velocity and scalar fields in the perpendicular wind were compared with those in the oblique wind. Furthermore, realistic wind conditions measured on the BT Tower at 190 m above street level were processed and used to drive the numerical simulations of flows and dispersion at the DAPPLE site. This leads to significant predictive improvements of the dispersion compared with field measurements, which provides validation and confidence for coupling mesoscale meteorological models, e.g. the UK Met Office’s Unified Model and the NCAR’s Weather Research & Forecasting Model, with the street-scale large-eddy simulation of urban environments.     

Subgrid-Scale Dynamics of Water Vapour, Heat, and Momentum over a Lake

We examine the dynamics of turbulence subgrid (or sub-filter) scales over a lake surface and the implications for large-eddy simulations (LES) of the atmospheric boundary layer. The analysis is based on measurements obtained during the Lake-Atmosphere Turbulent EXchange (LATEX) field campaign (August–October, 2006) over Lake Geneva, Switzerland. Wind velocity, temperature and humidity profiles were measured at 20 Hz using a vertical array of four sonic anemometers and open-path gas analyzers. The results indicate that the observed subgrid-scale statistics are very similar to those observed over land surfaces, suggesting that the effect of the lake waves on surface-layer turbulence during LATEX is small. The measurements allowed, for the first time, the study of subgrid-scale turbulent transport of water vapour, which is found to be well correlated with the transport of heat, suggesting that the subgrid-scale modelling of the two scalars may be coupled to save computational resources during LES.     

A wind tunnel study of turbulent dispersion over two- and three-dimensional gentle hills from upwind point sources in neutral flow

A study of turbulent dispersion over hills for upstream, elevated sources was conducted, based on wind tunnel tracer gas (CO₂) experiments over a gentle 2-D ridge and a 3-D circular hill, both having a cosine-square cross-section. The concentration measurements were made with four different source locations for each hill case (2-D or 3-D), and the study focused on dispersion parameters under the influence of the presence of the hills in order to provide a better understanding of the mechanisms involved.The wind tunnel measurements show that, in the case of gentle hills, the topographic impact on turbulent dispersion from upstream sources is only moderate and is more pronounced for the 3-D than for the 2-D hill. The perturbation in mean flow introduced by the hills, including streamline divergence/convergence, is shown to dominate the changes in the dispersion due to the hills in this case. The plume spread, both in the lateral and the vertical, is enhanced over the upwind hill foot and reduced over the hill top in response to the mean flow slow-down and speed-up at these places, and is further enhanced or reduced due to streamline divergence/convergence in the vertical over the hills as well as in the horizontal over the 3-D hill. These results are also compared with cases of turbulent dispersion over more steep hills (Snyder and Britter, 1987).     

Dissimilarity of Scalar Transport in the Convective Boundary Layer in Inhomogeneous Landscapes

A land-surface model (LSM) is coupled with a large-eddy simulation (LES) model to investigate the vegetation-atmosphere exchange of heat, water vapour, and carbon dioxide (CO₂) in heterogeneous landscapes. The dissimilarity of scalar transport in the lower convective boundary layer is quantified in several ways: eddy diffusivity, spatial structure of the scalar fields, and spatial and temporal variations in the surface fluxes of these scalars. The results show that eddy diffusivities differ among the three scalars, by up to 10–12%, in the surface layer; the difference is partly attributed to the influence of top-down diffusion. The turbulence-organized structures of CO₂ bear more resemblance to those of water vapour than those of the potential temperature. The surface fluxes when coupled with the flow aloft show large spatial variations even with perfectly homogeneous surface conditions and constant solar radiation forcing across the horizontal simulation domain. In general, the surface sensible heat flux shows the greatest spatial and temporal variations, and the CO₂ flux the least. Furthermore, our results show that the one-dimensional land-surface model scheme underestimates the surface heat flux by 3–8% and overestimates the water vapour and CO₂ fluxes by 2–8% and 1–9%, respectively, as compared to the flux simulated with the coupled LES-LSM.     

Turbulence in Sparse, Organized Vegetative Canopies: A Large-Eddy Simulation Study

A large-eddy simulation study was performed to characterize turbulence in sparse, row-oriented canopies. This was accomplished by simulating a set of heterogeneous row-oriented canopies with varying row vegetation density and spacing. To determine the effects of heterogeneity, results were compared to horizontally homogeneous canopies with an equivalent ‘effective’ leaf area index. By using a proper effective leaf area index, plane-averaged mean velocities and bulk scaling parameters contained only small errors when heterogeneity was ignored. However, many cases had significantly larger second- and third-order velocity moments in the presence of heterogeneity. Some heterogeneous canopies also contained dispersive fluxes in the lower canopy that were over 20 % as large as the turbulent flux. Impacts of heterogeneity were most pronounced in the cases of large row leaf area density and widely spaced rows. Despite the substantial amount of open space in the sparse canopies, vertical velocity skewness and quadrant-hole analysis indicated that the flow behaved predominantly as a canopy layer even though integral length scales at the canopy top no longer followed mixing-layer scaling. This was supported by the fact that similar composite-averaged coherent structures could be readily identified in both the heterogeneous and homogeneous canopies. Heterogeneity had an effect on coherent structures, in that structure detection events were most likely to occur just upwind of the vegetation rows. In simulations with large row spacing, these structures also penetrated deeper into the canopy when compared to the equivalent homogeneous canopy.     

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.     

Coherent Structures and the Dissimilarity of Turbulent Transport of Momentum and Scalars in the Unstable Atmospheric Surface Layer

Atmospheric stability effects on the dissimilarity between the turbulent transport of momentum and scalars (water vapour and temperature) are investigated in the neutral and unstable atmospheric surface layers over a lake and a vineyard. A decorrelation of the momentum and scalar fluxes is observed with increasing instability. Moreover, different measures of transport efficiency (correlation coefficients, efficiencies based on quadrant analysis and bulk transfer coefficients) indicate that, under close to neutral conditions, momentum and scalars are transported similarly whereas, as the instability of the atmosphere increases, scalars are transported increasingly more efficiently than momentum. This dissimilarity between the turbulent transport of momentum and scalars under unstable conditions concurs with, and is likely caused by, a change in the topology of turbulent coherent structures. Previous laboratory and field studies report that under neutral conditions hairpin vortices and hairpin packets are present and dominate the vertical fluxes, while under free-convection conditions thermal plumes are expected. Our results (cross-stream vorticity variation, quadrant analysis and time series analysis) are in very good agreement with this picture and confirm a change in the structure of the coherent turbulent motions under increasing instability, although the exact structure of these motions and how they are modified by stability requires further investigation based on three-dimensional flow data.     

Large-eddy Simulation of Turbulent Flow across a Forest Edge. Part II: Momentum and Turbulent Kinetic Energy Budgets

Momentum and turbulent kinetic energy (TKE) budgets across a forest edge have been investigated using large-eddy simulation (LES). Edge effects are observed in the rapid variation of a number of budget terms across this vegetation transition. The enhanced drag force at the forest edge is largely balanced by the pressure gradient force and by streamwise advection of upstream momentum, while vertical turbulent diffusion is relatively insignificant. For variance and TKE budgets, the most important processes at the forest edge are production due to the convergence (or divergence) of the mean flow, streamwise advection, pressure diffusion and enhanced dissipation by canopy drag. Turbulent diffusion, pressure redistribution and vertical shear production, which are characteristic processes in homogeneous canopy flow, are less important at the forest transition. We demonstrate that, in the equilibrated canopy flow, a substantial amount of TKE produced in the streamwise direction by the vertical shear of the mean flow is redistributed in the vertical direction by pressure fluctuations. This redistribution process occurs in the upper canopy layers. Part of the TKE in the vertical velocity component is transferred by turbulent and pressure diffusion to the lower canopy levels, where pressure redistribution takes place again and feeds TKE back to the streamwise direction. In this TKE cycle, the primary source terms are vertical shear production for streamwise velocity variance and pressure redistribution for vertical velocity variance. The evolution of these primary source terms downwind of the forest edge largely controls the adjustment rates of velocity variances.     

A comparison of wind-tunnel experiments and numerical simulations of neutral and stratified flow over a hill

A comparison is made of numerical and experimental results for flow over two-dimensional hills in both neutral and stably stratified flow. The numerical simulations are carried out using a range of one-and-a-half order and second-order closure schemes. The performance of the various turbulence schemes in predicting both the mean and turbulent quantities over the hill is assessed by comparing the results with new wind-tunnel measurements. The wind-tunnel experiments include both neutral and stably stratified flow over two different hills with different slopes, one of which is steep enough to induce flow separation. The dataset includes measurements of the mean and turbulent parts of the flow using laser Doppler anemometry. Pressure measurements are also made across the surface of the hill. These features make the dataset an excellent test of the model performance. In general second-order turbulence schemes provide the best agreement with the experimental data, however, they can be numerically unstable for steep hills. Some modifications can be made to the standard one-and-a-half order closure scheme, which results in improved performance at a fraction of the computation cost of the second-order schemes.