Numerical simulation of flow data over two-dimensional hills

We have studied a neutrally-stratified flow over two-dimensional hills using a two-dimensional, non-hydrostatic version of the Regional Atmospheric Modeling System (RAMS). We have implemented three different turbulence closure models: the standardE- model, an Algebraic Reynolds Stress Model (ARSM) and a new model. Model predictions for the mean and turbulence flows using different closure schemes are compared with the data of a wind tunnel experiment containing isolated two-dimensional hills of varying slope. From the comparison, it is concluded that all three models predict the mean flow velocities equally well while only the new closure model accurately predicts the turbulence data statistics.The research reported in this paper was conducted while the first author held a National Research Council (NRC) Associateship.     

Numerical Simulation of Flow over Two-Dimensional Hills Using a Second-Order Turbulence Closure Model

We study turbulent flow over two-dimensional hills. The Reynolds stresses are represented by a second-order closure model, where advection, diffusion, production and dissipation processes are all accounted for. We solve a full set of primitive non-hydrostatic dynamic equations for mean flow quantities using a finite-difference numerical method. The model predictions for the mean velocity and Reynolds stresses are compared with the measured data from a wind-tunnel experiment that simulates the atmospheric boundary layer. The agreement is good. The performance of the second-order closure model is also compared withthat of lower level turbulence models, including the eddy-viscositymodel and algebraic Reynolds stress models. It is concluded that thepresent closure is a considerable improvement over the other modelsin representing various physical effects in flow over hills. Thefeasibility of running a finite-difference numerical simulationincorporating a full second-order closure model on an IBM workstationis also demonstrated.     

A wind tunnel study of turbulent flow over model hills

Detailed wind tunnel measurements have been made of mean flow and turbulence over a two-dimensional ridge and a circular hill, both having cosine-squared cross-section and maximum slope about 15 °. The measurements were made in an artificially thickened neutrally stratified boundary layer, and have been compared with results from linear models and rapid distortion theory as appropriate.Our study shows that linear theory gives generally good predictions of the mean flow on the upwind side of the hills, and especially of the flow speedup at the hill top, but that the turbulence is less well predicted. In particular, the measurements show a major increase in the vertical component of turbulence and in the shear stress on the upwind slope of both the two- and three-dimensional hills which is not predicted by either equilibrium or isotropic rapid-distortion theories, although this may be partly due to the effect of streamline curvature. Rapid-distortion theory is successful only in describing the streamwise component of turbulence in the outer region of the flow, while in the upper part of the inner region of the flow, the turbulence measurements show disagreement with both the equilibrium and the rapid-distortion theories. Our experiments also confirm that the equilibrium region is a very thin layer close to the surface, while above this region and below the outer region, there is a transitional region where all terms in the stress equation are important.The measurements over the three-dimensional hill suggest that the mean flow and turbulence are broadly similar to those over the two-dimensional ridge, but with reduced perturbation amplitudes. The major differences between the two cases are found on the upwind slope and in the wake where, respectively, horizontal divergence and convergence of the three-dimensional flow are most pronounced.     

Numerical Modelling of the Turbulent Flow Developing Within and Over a 3-D Building Array, Part I: A High-Resolution Reynolds-Averaged Navier—Stokes Approach

A study of the neutrally-stratified flow within and over an array of three-dimensional buildings (cubes) was undertaken using simple Reynolds-averaged Navier—Stokes (RANS) flow models. These models consist of a general solution of the ensemble-averaged, steady-state, three-dimensional Navier—Stokes equations, where the k-ε turbulence model (k is turbulence kinetic energy and ε is viscous dissipation rate) has been used to close the system of equations. Two turbulence closure models were tested, namely, the standard and Kato—Launder k-ε models. The latter model is a modified k-ε model designed specifically to overcome the stagnation point anomaly in flows past a bluff body where the standard k-ε model overpredicts the production of turbulence kinetic energy near the stagnation point. Results of a detailed comparison between a wind-tunnel experiment and the RANS flow model predictions are presented. More specifically, vertical profiles of the predicted mean streamwise velocity, mean vertical velocity, and turbulence kinetic energy at a number of streamwise locations that extend from the impingement zone upstream of the array, through the array interior, to the exit region downstream of the array are presented and compared to those measured in the wind-tunnel experiment. Generally, the numerical predictions show good agreement for the mean flow velocities. The turbulence kinetic energy was underestimated by the two different closure models. After validation, the results of the high-resolution RANS flow model predictions were used to diagnose the dispersive stress, within and above the building array. The importance of dispersive stresses, which arise from point-to-point variations in the mean flow field, relative to the spatially-averaged Reynolds stresses are assessed for the building array.     

Estimation of Roughness Parameters Within Sparse Urban-Like Obstacle Arrays

We conduct wind-tunnel experiments on three different uniform roughness arrays composed of sparsely distributed rectangular cylinders for the estimation of surface parameters. Roughness parameters such as the roughness length z ₀ and zero-plane displacement d are extracted using a best-fit approximation of the measured wind velocity. We also perform a large-eddy simulation (LES) to confirm that four sampling points are sufficient to surrogate a space average above the canopy layer of the sparse roughness arrays. We propose a new morphological model from a systematic analysis of experimental data on the arrays. The friction velocity predicted by the proposed model agrees well with the peak value of the measured Reynolds shear stress ${(-\left<\overline{u'w'}\right>)^{0.5}}${(-\left<\overline{u'w'}\right>)^{0.5}}. The proposed model is further validated in an additional wind-tunnel experiment conducted on a scaled configuration of a real urban area exposed to four wind directions. The proposed model is found to perform very well particularly in the estimation of the friction velocity, readily leading to a better estimation of turbulence, which is essential for an accurate prediction of pollutant dispersion.     

A two-dimensional trajectory-simulation model for non-Gaussian,inhomogeneous turbulence within plant canopies

We report a two-dimensional (alongwind u, vertical w) trajectory-simulation model, consistent with Thomson's (1987) well-mixed criteria, that allows for the non-Gaussian turbulence typical of flow within a plant canopy. The effect of non-Gaussian turbulence was examined by formulating a non-Gaussian u, w joint probability density function (PDF) as the sum of two Gaussian joint-PDFs. The resultant PDF reproduced the desired means, variances, skewnesses, and kurtoses, and the correct covariance. In prediction of the location of maximum concentration downwind of a line source in homogeneous, slightly non-Gaussian turbulence, it proved advantageous to incorporate skewness and kurtosis. However, in the case of inhomogeneous, highly non-Gaussian turbulence, the addition of skewness and kurtosis in the model resulted in substantially worse agreement with measurements than the results of the model using Gaussian PDFs. This may be due to inaccuracy in our PDF formulation. Dispersion predictions from the model with Gaussian PDFs were generally not statistically different from measurements. These results indicate that a two-dimensional Gaussian trajectory-simulation approach is adequate to predict mean concentrations and fluxes resulting from sources within plant canopies.     

Influence of non-flat terrain and wind direction shear on canopy turbulence

We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk _w and r _uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.     

Numerical studies of neutrally stratified planetary boundary-layer flow above gentle topography

A numerical model of planetary boundary-layer flow above two-dimensional gentle topography is developed as an extension of the surface-layer model described by Taylor (1977). Comparisons are made with surface-layer predictions for flow over Gaussian hills; and the flow at various angles above hills, valleys and escarpments is modelled. Some simple case studies of the influence of gentle two-dimensional topography on pollutant dispersion are made which indicate relatively minor effects on surface pollutant concentrations in comparison with results for dispersion above a plane surface.     

A Flume Experiment on the Adjustment of the Mean and Turbulent Statistics to a Transition from Short to Tall Sparse Canopies

Water-flume experiments are conducted to study the structure of turbulent flow within and above a sparse model canopy consisting of two rigid canopies of different heights. This difference in height specifies a two-dimensional step change from a rough to a rougher surface, as opposed to a smooth-to-rough transition. Despite the fact that the flow is in transition from a rough to a rougher surface, the thickness of the internal boundary layer scales as x ^4/5, consistent with smooth-to-rough boundary layer adjustment studies, where x is the downstream distance from the step change. However, the analogy with smooth-to-rough transitions no longer holds when the flow inside the canopy and near the canopy top is considered. Results show that the step change in surface roughness significantly increases turbulence intensities and shear stress. In particular, there is an adjustment of the mean horizontal velocity and shear stress as the flow passes over the rougher canopy, so that their vertical profiles adjust to give maximum values at the top of this canopy. We also observe that the magnitude and shape of the inflection in the mean horizontal velocity profile is significantly affected by the transition. The horizontal and vertical turbulence spectra compare well with Kolmogorov’s theory, although a small deviation at high frequencies is observed in the horizontal spectrum within the canopy. Here, for relatively low leaf area index, shear is found to be a more effective mechanism for momentum transfer through the canopy structure than vortex shedding.     

Water tank and numerical model studies of flow over steep smooth two-dimensional hills

The present work investigates the role of different treatments of the lower boundary condition on the numerical prediction of flows over two-dimensional, smooth, steep hills. Four different law of the wall formulations are tested when a large recirculating region is formed on the lee side of the hill. Numerical implementation of the near-wall functions was made through a finite elements code. The standard κ–ε model was used to close the averaged Navier–Stokes equations. Results are validated through original data obtained in a water tank. Measurements resorted to laser Doppler anemometry. The experiment provide detailed data for the characterization of the reverse flow in the region between the separation and the reattachment points, with emphasis on the near wall region. The experimental wall shear stress distribution is compared with the results provided by the different laws of the wall showing good agreement. The numerical predictions are shown to vary markedly between different formulations.     

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.     

The onset of separation in neutral,turbulent flow over hills

The onset of separation in turbulent, neutrally stratified, boundary-layer flow over hills is considered. Since the flows are fully turbulent, the occurrence of intermittent separation, in the sense of any reversal of near surface flow, will depend strongly on the detailed structure and behaviour of the turbulent eddies. Very little is known about such intermittent separation and the phenomenon cannot be studied with numerical models employing standard turbulence closures; eddy-resolving models are required. Therefore, here, as elsewhere in the literature, the arguably less physically significant process of mean flow separation is studied. Numerical simulations of flow over idealised two- and three-dimensional hills are examined in detail to determine the lowest slope, _crit, for which the mean flow separates.Previous work has identified this critical slope as that required to produce a zero surface stress somewhere over the hill. This criterion, when a mixing-length turbulence closure is applied, reduces to requiring the near-surface vertical velocity shear to vanish at some point on the hill's surface. By applying results from a recent linear analysis for the flow perturbations to this condition, a new expression for _crit is obtained. The expression is approximate but its relative simplicity makes it practically applicable without the need for use of a computer or for detailed mapping of the hill. The approach suggested differs from previous ones in that it applies linear results to a non-linear expression for the surface stress. In the past, a linear expression for the surface stress has been used. The proposed expression for _crit leads to critical angles that are about twice previous predictions. It is shown that the present expression gives good agreement with the numerical results presented here, as well as with other numerical and experimental results. It is also consistent with atmospheric observations.     

The Effects of Probe-Induced Flow Distortion on Velocity Covariances: Field Observations

We carried out measurements to test a simple theory of the effect of probe-induced flow distortion on turbulence measurements. We used two three-component sonic anemometers mounted 1.8m apart at a height of 6.7 m. Behind one was a horizontal circular cylinder of radius 0.15 m and length 1.2 m, chosen to model two-dimensional probe-induced flow distortion in the limit where the scale of the turbulence is very large compared to the scale of the probe. The second sonic anemometer measured the undistorted flow. The measured flow-distortion effects on the Reynolds shearing stress and the variances of streamwise and vertical velocity agree well with the theory.     

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.     

A group-kinetic theory with a closure by memory loss for modeling turbulence in the atmosphere and the oceans

Summary The principle of the group-kinetic method is elucidated. This method of renormalization serves as the basis for analyzing the spectral structure of turbulence. The spectral distributions include the Kolmogoroff lawk ^–5/3 for isotropic turbulence, the power lawk ^–1 for shear turbulence, the spectrum for stratified turbulence not in the power law form, the power lawk ^–3 for two-dimensional geostrophic turbulence, and the power lawsk ^–3,k ^–2 andk ^–5 for two-dimensional Rossby wave turbulence with uniform and differential rotations. We discuss a spectrum-dependent modeling in reference to the problems of the universal functions and parameters in the similarity theories for the atmospheric surface layer and the planetary boundary layer. A renormalization-based modeling of atmospheric turbulence is proposed.     

Velocity and Surface Shear Stress Distributions Behind a Rough-to-Smooth Surface Transition: A Simple New Model

A simple new model is proposed to predict the distribution of wind velocity and surface shear stress downwind of a rough-to-smooth surface transition. The wind velocity is estimated as a weighted average between two limiting logarithmic profiles: the first log law, which is recovered above the internal boundary-layer height, corresponds to the upwind velocity profile; the second log law is adjusted to the downwind aerodynamic roughness and local surface shear stress, and it is recovered near the surface, in the equilibrium sublayer. The proposed non-linear form of the weighting factor is equal to ln(z/z ₀₁)/ln(δ _i /z ₀₁), where z, δ _i and z ₀₁ are the elevation of the prediction location, the internal boundary-layer height at that downwind distance, and the upwind surface roughness, respectively. Unlike other simple analytical models, the new model does not rely on the assumption of a constant or linear distribution for the turbulent shear stress within the internal boundary layer. The performance of the new model is tested with wind-tunnel measurements and also with the field data of Bradley. Compared with other existing analytical models, the proposed model shows improved predictions of both surface shear stress and velocity distributions at different positions downwind of the transition.     

Wind flow over ridges in simulated atmospheric boundary layers

The flows over four two-dimensional triangular hills and three two-dimensional bell-shaped hills have been investigated in a simulated rural atmospheric boundary layer modelled to a scale of 1:300: Further measurements were made over two of the triangular hills in a simulated rural boundary layer of 1: 3000 scale and in a simulated urban boundary layer modelled to a scale of 1:400. The effect of the model hill surface roughness was also investigated. Flow measurements were restricted to the mean velocity U, RMS velocity fluctuations u and the energy spectra for the streamwise velocity component Measurements were made at a number of longitudinal positions in the approach flow, over the model hills and downstream of the model hills. For each model hill, the crest was the region of largest mean velocity and smallest velocity fluctuations. The largest mean velocities over the model hills occurred for hills of intermediate slope rather than for the steepest hills. A decrease in the scale of the simulated atmospheric boundary layer led to a reduction in the amplification factors at the hill crests, whereas an increase in the surface roughness of the approach flow resulted in increased amplification factors at the hill crests.     

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.     

The Use of Weather Forecasts to Characterise Near-Surface Optical Turbulence

The propagation of optical and electromagnetic waves is affected by small-scale atmospheric turbulence, quantified by the structure parameter of the refractive index. In the atmospheric surface layer, the mean structure parameter C_n²{C_{n}^{2}} , as averaged over the large-scale turbulence, relates to meteorological forcings through well-documented relationships. Present-day numerical weather forecast models routinely produce these forcings at the global scale. This study introduces a method where the products of such a model are used to calculate the mean optical turbulence near the surface. The method is evaluated against scintillometry measurements over climatologically distinct sites in Western Europe. The diurnal cycle modulation, and regional and seasonal contrasts, are all reproduced by our predictions. Hence, the present method explains and predicts some essential aspects of the meteorological variability of C_n²{C_{n}^{2}} near the surface. The noted discrepancies combine instrumental limitations, site peculiarities, differences related to distinct averaging procedures, and model errors, notably from weather forecasts. The minute-scale fluctuations of the measured scintillation rate are also analysed in the light of the forecast weather conditions. Fair-weather daytime periods consistently show a small short-term variability compared to the nighttime and perturbed weather periods. Thus, this short-term variability appears to have a predictable component.     

Analysis of wind tunnel dispersion data over two-dimensional obstacles

An analysis of wind tunnel data of dispersion from elevated point sources over a flat floor and gently sloping, 2-D hills is performed. (The data were obtained by Khurshudyan et al., 1981.) Formulas for the mixing lengths and the Lagrangian time scales are tested, suitable for use in various dispersion models. Some expressions for the vertical first moments of the plume concentration distribution suitable for shear flow (Hunt, 1985) are also tested.Then, a normalization is suggested, based on the source mean flow and turbulence parameters, for the ground-level concentration maximum value and position along the plume centerline. Using this normalization, the maximum position is almost constant, regardless of source height variation and of whether the hill is present or not, at least for source positions upstream or over the top.The maximum values allow the determination of normalized terrain amplification factors TAF_N, which are shown to be in most cases closer to one than the corresponding TAF obtained without normalization.