Large-Eddy Simulation of Turbulent Flow Over a Rough Surface

A family of wall models is proposed that exhibits moresatisfactory performance than previousmodels for the large-eddy simulation (LES) of the turbulentboundary layer over a rough surface.The time and horizontally averaged statistics such asmean vertical profiles of windvelocity, Reynolds stress, turbulent intensities, turbulentkinetic energy and alsospectra are compared with wind-tunnel experimental data.The purpose of the present study is to obtain simulatedturbulent flows that are comparable with wind-tunnelmeasurements for use as the wind environment for thenumerical prediction by LES of source dispersion in theneutral atmospheric boundary layer.     

Large-Eddy Simulation Of Turbulent Separated Flow Over Rough Hills

Large-eddy simulations (LES) have been performed ofneutral turbulent flow over two-dimensional ridges steepenough to cause separation. Both periodic and isolated ridges havebeen considered. The results are compared with wind-tunnel observations and with the predictions of various turbulence closure models.For the periodic case the LES results are qualitatively reasonable,although the depth of the separated region appears to besensitive to the use of a distributed drag near the lower boundary.The isolated ridge results compare very favourably with the experimentaldata, with the LES performance appearing to be at least as good as that ofthe closure models.     

A Large-Eddy Simulation Study of Turbulent Flow Over Multiscale Topography

Most natural landscapes are characterized by multiscale (often multifractal) topography with well-known scale-invariance properties. For example, the spectral density of landscape elevation fields is often found to have a power-law scaling behaviour (with a −2 slope on a log–log scale) over a wide span of spatial scales, typically ranging from tens of kilometres down to a few metres. Even though the effect of topography on the atmospheric boundary layer (ABL) has been the subject of numerous studies, few have focussed on multiscale topography. In this study, large-eddy simulation (LES) is used to investigate boundary-layer flow over multiscale topography, and guide the development of parametrizations needed to represent the effects of subgrid-scale (SGS) topography in numerical models of ABL flow. Particular emphasis is placed on the formulation of an effective roughness used to account for the increased aerodynamic roughness associated with SGS topography. The LES code uses the scale-dependent Lagrangian dynamic SGS model for the turbulent stresses and a terrain-following coordinate transformation to explicitly resolve the effects of the topography at scales larger than the LES resolution. The terrain used in the simulations is generated using a restricted solid-on-solid landscape evolution model, and it is characterized by a −2 slope of the elevation power spectrum. Results from simulations performed using elevation fields band-pass filtered at different spatial resolutions indicate a clear linear relation between the square of the effective roughness and the variance of elevation.     

Numerical Simulation of Strongly Stratified Flow Over a Three-Dimensional Hill

A numerical study of stably stratified flow over a three-dimensional hill is presented. Large-eddy simulation is used here to examine in detail the laboratory experimental flows described in the landmark work of Hunt and Snyder about stratified flow over a hill. The flow is linearly stratified and U/Nh is varied from 0.2 to 1.0. Here N and U are the buoyancy frequency and freestream velocity respectively, and h is the height of the hill. The Reynolds number based on the hill height is varied from 365 to 2968. The characteristic flow patterns at various values of U/Nh have been obtained and they are in good agreement with earlier theoretical and experimental results. It is shown that the flow field cannot be predicted by Drazin's theory when recirculation exists at the leeside of the hill even at UNh 1. The wake structure agrees well with a two-dimensional wake assumption when U/Nh 1 but lee waves start to influence the wake structure as U/Nh increases. The dividing-streamline heights obtained in the simulation are in accordance with experimental results and Sheppard's formula. The energy loss along the dividing streamline due to friction/turbulence approximately offsets the energy gained from pressure field. When lee waves are present, linear theory always underestimates the amplitude and overestimates the wavelength of three-dimensional lee waves. The simulated variations of drag coefficients with the parameterK (=ND/ U) are qualitatively consistent with experimental data and linear theory. Here D is the depth of the tank.     

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.     

Large-Eddy Simulation of Flow Over a Vegetation-Like Canopy Modelled as Arrays of Bluff-Body Elements

Turbulent flow over a vegetation canopy under neutral atmospheric conditions is investigated using large-eddy simulation. Each model tree, which consists of a sphere-shaped tree crown and a cylindrical trunk, is fully resolved. The resulting turbulence statistics and the drag force on the vegetation agree well with measurements from the corresponding wind-tunnel experiment described by Böhm et al. (Boundary-Layer Meteorol, 146:393–419, 2013). Statistically, this kind of model canopy exhibits both vegetation and bluff-body-flow characteristics. The time-averaged flow skims over the top of the underlying canopy, forming a low-momentum recirculation zone on the lee-side of the bluff elements, which causes significant dispersive stress within the canopy layer. Two other numerical representations of vegetation canopies, referred to as the drag-element and drag-crown approaches, have also been developed to assess the performance of simulations. Turbulence statistics suggest that the canopy shear layer interferes with wakes behind stems and crowns. The drag-crown approach yields better agreement between numerical results and experimental measurements than does the traditional drag-element approach, thus providing a promising numerical model for simulating canopy turbulence.     

Large-Eddy Simulation of Atmospheric Boundary-Layer Flow Over Fluvial-Like Landscapes Using a Dynamic Roughness Model

A recently developed dynamic surface roughness model (Anderson and Meneveau, J Fluid Mech 679:288–314, 2011) for large-eddy simulation (LES) of atmospheric boundary-layer flow over multi-scale topographies is applied to boundary-layer flow over several types of fluvial-like landscapes. The landscapes are generated numerically with simulation of a modified Kardar–Parisi–Zhang equation (Passalacqua et al., Water Resour Res 42:WOD611, 2006). These surfaces possess the fractal-like channel network and anisotropic features often found in real terrains. The dynamic model is shown to lead to accurate flow predictions when the surface-height distributions exhibit power-law scaling (scale invariance) in the prevalent mean flow direction. In those cases, the LES provide accurate predictions (invariant to resolution) of mean velocity profiles. Conversely, some resolution dependence is found for applications in which the landscape’s streamwise spectra do not exhibit pure power-law scaling near wavenumbers corresponding to the LES grid resolution.     

Large-Eddy Simulation of Windbreak Flow

A large-eddy simulation has been performed of turbulent flow around multiple windbreaks set within a wheat canopy under neutral stability conditions. The simulation is validated against a wind tunnel data set taken under similar conditions. Velocity profiles and second-order statistics are presented and compared to those found in the wind tunnel. From the numerical simulation, we discuss spatial distributions of instantaneous velocity fields and pressure statistics, which are important and telling features of the flow that are difficult to measure experimentally. We present a discussion of the momentum balance at various locations with respect to the windbreak, and similarly, we introduce the budget of a passive scalar. These discussions show the importance of the terms in each budget equation as they vary upstream and downstream of the windbreak.     

Large-Eddy Simulation of Coherent Flow Structures within a Cubical Canopy

Instantaneous flow structures “within” a cubical canopy are investigated via large-eddy simulation. The main topics of interest are, (1) large-scale coherent flow structures within a cubical canopy, (2) how the structures are coupled with the turbulent organized structures (TOS) above them, and (3) the classification and quantification of representative instantaneous flow patterns within a street canyon in relation to the coherent structures. We use a large numerical domain (2,560 m × 2,560 m × 1,710 m) with a fine spatial resolution (2.5 m), thereby simulating a complete daytime atmospheric boundary layer (ABL), as well as explicitly resolving a regular array of cubes (40 m in height) at the surface. A typical urban ABL is numerically modelled. In this situation, the constant heat supply from roof and floor surfaces sustains a convective mixed layer as a whole, but strong wind shear near the canopy top maintains the surface layer nearly neutral. The results reveal large coherent structures in both the velocity and temperature fields “within” the canopy layer. These structures are much larger than the cubes, and their shapes and locations are shown to be closely related to the TOS above them. We classify the instantaneous flow patterns in a cavity, specifically focusing on two characteristic flow patterns: flushing and cavity-eddy events. Flushing indicates a strong upward motion, while a cavity eddy is characterized by a dominant vortical motion within a single cavity. Flushing is clearly correlated with the TOS above, occurring frequently beneath low-momentum streaks. The instantaneous momentum and heat transport within and above a cavity due to flushing and cavity-eddy events are also quantified.     

Dispersive Fluxes Within and Over a Real Urban Canopy: A Large-Eddy Simulation Study

Boundary-Layer Meteorology - Large-eddy simulations (LES) are conducted to study the transport of momentum and passive scalar within and over a real urban canopy in the City of Boston, USA. This...     

Three-Dimensional Scalar Microfront Systems in a Large-Eddy Simulation of Vegetation Canopy Flow

Large-eddy simulation is used to reproduce neutrallystratified airflow inside and immediately above a vegetation canopy. A passive scalaris released from the canopy and the evolution of scalar concentration above the canopyis studied. The most significant characteristic of the scalar concentration is the repeatedformation and dissipation of scalar microfronts, a phenomenon that has been observedin nature. These scalar microfronts consist of downstream-tilted regions of highscalar concentration gradients. Computer visualization tools and a conditional samplingand compositing technique are utilized to analyze these microfronts. Peaks in positivepressure perturbation exceeding an experimental threshold are found to be effectiveindicators of scalar microfronts. Convergence of the streamwise velocity componentand divergence of the cross-stream velocity component are observed in the immediatevicinity of scalar microfronts, which helps explain their relatively longlifetimes. Many of these three-dimensional features have been observedin previous field studies of canopy flow.     

Large-Eddy Simulation Comparison of Neutral Flow Over a Canopy: Sensitivities to Physical and Numerical Conditions,and Similarity to Other Representations

The representation of a neutral atmospheric flow over roughness elements simulating a vegetation canopy is compared between two large-eddy simulation models, wind-tunnel data and recently updated empirical flux-gradient relationships. Special attention is devoted to the dynamics in the roughness sublayer above the canopy layer, where turbulence is most intense. By demonstrating that the flow properties are consistent across these different approaches, confidence in the individual independent representations is bolstered. Systematic sensitivity analyses with the Dutch Atmospheric Large-Eddy Simulation model show that the transition in the one-sided plant-area density from the canopy layer to unobstructed air potentially alters the flow in the canopy and roughness sublayer. Anomalously induced fluctuations can be fully suppressed by spreading the transition over four steps. Finer vertical resolutions only serve to reduce the magnitude of these fluctuations, but do not prevent them. To capture the general dynamics of the flow, a resolution of 10 % of the canopy height is found to suffice, while a finer resolution still improves the representation of the turbulent kinetic energy. Finally, quadrant analyses indicate that momentum transport is dominated by the mean velocity components within each quadrant. Consequently, a mass-flux approach can be applied to represent the momentum flux.     

Non-Linear, Microscale Modelling of the Flow Over Askervein Hill

The flow solver “3DWind” is used to explore new aspects of the Askervein hill flow case. Previous work has investigated sensitivities to the grid, the inflow boundary profile, the roughness and the turbulence model. Several different linear and non-linear numerical models have also been validated by means of the Askervein hill case. This analysis focuses on the flow sensitivity to the grid spacing, the incident wind direction and the vertical resolution of topographic data. The horizontal resolution is found to be fine enough to cause only minor differences compared to a grid where every second node is removed. The vertical resolution dependence is mainly attributed to the wall functions. Simulations are performed for wind directions 200°, 205°, 210° and 215° at the reference station. The smallest directional biases compared to experimental values along a line through the hilltop are found for the directions 200° and 205°. There are larger wind direction changes along this line through the hilltop in the 200° case than in the 215° case. Still the simulation results give less veering than found in the experimental results, and this is maybe caused by a slightly stable atmosphere. The sensitivity to the vertical resolution of the topographical data is found to be particularly high close to the ground at the top of the hill; this is where the speed-up is most important. Differences decrease with the height from the ground. At higher levels the speed-ups are smaller and caused by terrain formations with larger scales.     

Large-Eddy Simulation of Thermally Stratified Atmospheric Boundary-Layer Flow Using a Minimum Dissipation Model

A generalized form of a recently developed minimum dissipation model for subfilter turbulent fluxes is proposed and implemented in the simulation of thermally stratified atmospheric boundary-layer flows. Compared with the original model, the generalized model includes the contribution of buoyant forces, in addition to shear, to the production or suppression of turbulence, with a number of desirable practical and theoretical properties. Specifically, the model has a low computational complexity, appropriately switches off in laminar and transitional flows, does not require any ad hoc shear and stability corrections, and is consistent with theoretical subfilter turbulent fluxes. The simulation results show remarkable agreement with well-established empirical correlations, theoretical predictions, and field observations in the atmosphere. In addition, the results show very little sensitivity to the grid resolution, demonstrating the robustness of the model in the simulation of the atmospheric boundary layer, even with relatively coarse resolutions.     

Similarity Scaling Over a Steep Alpine Slope

In this study, we investigate the validity of similarity scaling over a steep mountain slope (30–41 $^\circ $ ). The results are based on eddy-covariance data collected during the Slope Experiment near La Fouly (SELF-2010); a field campaign conducted in a narrow valley of the Swiss Alps during summer 2010. The turbulent fluxes of heat and momentum are found to vary significantly with height in the first few metres above the inclined surface. These variations exceed by an order of magnitude the well-accepted maximum 10 % required for the applicability of Monin–Obukhov similarity theory in the surface layer. This could be due to a surface layer that is too thin to be detected or to the presence of advective fluxes. It is shown that local scaling can be a useful tool in these cases when surface-layer theory breaks down. Under convective conditions and after removing the effects of self-correlation, the normalized standard deviations of slope-normal wind velocity, temperature and humidity scale relatively well with $z/\varLambda $ , where $z$ is the measurement height and $\varLambda (z)$ the local Obukhov length. However, the horizontal velocity fluctuations are not correlated with $z/\varLambda $ under all stability regimes. The non-dimensional gradients of wind velocity and temperature are also investigated. For those, the local scaling appears inappropriate, particularly at night when shallow drainage flows prevail and lead to negative wind-speed gradients close to the surface.     

Flow over Hills: A Large-Eddy Simulation of the Bolund Case

Simulation of local atmospheric flows around complex topography is important for several applications in wind energy (short-term wind forecasting and turbine siting and control), local weather prediction in mountainous regions and avalanche risk assessment. However, atmospheric simulation around steep mountain topography remains challenging, and a number of different approaches are used to represent such topography in numerical models. The immersed boundary method (IBM) is particularly well-suited for efficient and numerically stable simulation of flow around steep terrain. It uses a homogenous grid and permits a fast meshing of the topography. Here, we use the IBM in conjunction with a large-eddy simulation (LES) and test it against two unique datasets. In the first comparison, the LES is used to reproduce experimental results from a wind-tunnel study of a smooth three-dimensional hill. In the second comparison, we simulate the wind field around the Bolund Hill, Denmark, and make direct comparisons with field measurements. Both cases show good agreement between the simulation results and the experimental data, with the largest disagreement observed near the surface. The source of error is investigated by performing additional simulations with a variety of spatial resolutions and surface roughness properties.     

Edge Flow and Canopy Structure: A Large-Eddy Simulation Study

Sharp heterogeneities in forest structure, such as edges, are often responsible for wind damage. In order to better understand the behaviour of turbulent flow through canopy edges, large-eddy simulations (LES) have been performed at very fine scale (2 m) within and above heterogeneous vegetation canopies. A modified version of the Advanced Regional Prediction System (ARPS), previously validated in homogeneous conditions against field and wind-tunnel measurements, has been used for this purpose. Here it is validated in a simple forest-clearing-forest configuration. The model is shown to be able to reproduce accurately the main features observed in turbulent edge flow, especially the “enhanced gust zone” (EGZ) present around the canopy top at a few canopy heights downwind from the edge, and the turbulent region that develops further downstream. The EGZ is characterized by a peak in streamwise velocity skewness, which reflects the presence of intense intermittent wind gusts. A sensitivity study of the edge flow to the forest morphology shows that with increasing canopy density the flow adjusts faster and turbulent features such as the EGZ become more marked. When the canopy is characterized by a sparse trunk space the length of the adjustment region increases significantly due to the formation of a sub-canopy wind jet from the leading edge. It is shown that the position and magnitude of the EGZ are related to the mean upward motion formed around canopy top behind the leading edge, caused by the deceleration in the sub-canopy. Indeed, this mean upward motion advects low turbulence levels from the bottom of the canopy; this emphasises the passage of sudden strong wind gusts from the clearing, thereby increasing the skewness in streamwise velocity as compared with locations further downstream where ambient turbulence is stronger.