Research Article  Open Access
FullScale Experimental Validation of LargeEddy Simulation of Wind Flows over Complex Terrain: The Bolund Hill
Abstract
Numerical simulation of local atmospheric flows around a complex topography is important for many wind energy applications, as it can help in locating and optimizing wind farms. The aim of this paper is to simulate the atmospheric flows over the challenging and complex site of Bolund Hill (in Denmark) using LargeEddy Simulations (LES) and to validate the simulation methodology against fullscale measurements. Wind prediction over a potential inland wind park is the main application of the validated LES methodology. For the first time, LESbased results from more than one wind direction of the Bolund case are reported and analyzed in detail herein. The LES results from each direction are compared with the Bolund field data in a quantitative way, for example, by using simulation error. According to the results comparison, it is found that the LES model provides by far the most accurate prediction for turbulent kinetic energy. The simulation error of this LES model for predicting turbulent kinetic energy is 19% smaller than that of all other models, whether experimental or numerical, employed previously over Bolund Hill.
1. Introduction
The production of wind energy is of everincreasing interest due to risks of global warming and nuclear energy production plant accidents. Nowadays, many wind farms are being constructed in areas with complex topographies containing hills, ridges, forests, and mountains, with the intention of reaching higher wind speed. Despite the advantages, establishing a turbine over a hill presents specific difficulties, such as a considerable surge in structural loads on wind turbine blades due to inherent flow variability, largescale unsteadiness, and flow separation, especially in complex terrains. In such terrains, the evaluation of wind resources becomes more and more challenging, and the prediction of windfarm profitability is directly affected. Therefore, an accurate modelling of atmospheric flows over a complicated terrain has been the subject of considerable research over the past few decades. A field measurement campaign is the most accurate way to carry out wind measurements since it accounts for all environmental parameters and weather fluctuations. However, field measurements are limited to only point locations that may not provide information over the entire site. This is particularly true in the case of complex terrains. Therefore, other alternatives should be taken into account for the scale of a large wind farm. To fulfill this requirement, several numerical and experimental models are being developed to predict the wind on a larger scale. In this way, models can provide valuable information for locating and optimizing wind farms in the absence of field data. Also, models can be employed to simulate the presence of wind turbines and their interaction with the flow. However, all these modelling approaches need to be validated against fullscale measurements.
Among the numerical models, the Reynoldsaveraged NavierStokes (RANS) approach with twoequation turbulence models has been widely used for wind prediction over complex terrains [1–9]. However, it must be mentioned that the RANS models provide solutions on the mean wind speed and the mean level of turbulent kinetic energy due to the Reynoldsaveraged approach [10]. Additionally, the models are dependent on a number of constants which are typically calibrated for simple, canonical flows, such as turbulent channel flow [11]. Therefore, the RANS approach still has some fundamental limitations when used to predict flows with complex phenomena. For example, as pointed out by Sumner et al. [12], the RANS formulation is inherently incapable of capturing unsteady effects, such as intermittent flow separation and eddy generation and transport, and generally gives difficulties while modelling turbulence in areas of strong separation [12]. On the other hand, the LargeEddy Simulation (LES) approach resolves the large (and thus the most important) turbulent structures, and it models only the smallscale motions which are more universal than the large eddies. As a result, LES is able to capture the important unsteadiness in flows. In contrast to RANS, LES typically requires a much higher spatial resolution, resulting in substantially longer computational time [11]. However, due to rapidly increasing computer technology, the use of LES towards atmospheric flow simulation is becoming increasingly common nowadays. For wind energy applications, LES is a promising tool, as it can adequately resolve the most important eddies of the flow which are relevant to wake evolution and power production [13].
The aim of this work is to simulate atmospheric flows over a natural complex terrain and to validate the LES methodology against real measurements. Wind prediction over a potential inland wind park is the main application of the validated LES methodology. However, unlike idealized terrains [14–17], the application of LES over realistic complex topography is not so common due to several numerical challenges, such as the realistic upstream boundary condition, heterogeneous nature of surface roughness, resolving nearwall turbulence structures, and highly turbulent wind flow. Particularly, reproducing the inflow boundary condition is one of the critical challenges, as it is generally not known before a numerical simulation. Further, the inflow condition in LES must contain turbulent fluctuations as a function of space and time, with a realistic energy distribution over the spatial directions and the simulated wavenumber range. In this study, the socalled recycling technique is employed for generating the instantaneous inflow boundary condition. Recently, the recycling technique was systematically studied and validated against the field measurements [18].
Several field experiments over natural terrains have been performed to provide fullscale measurements for validating flow modelling approaches. The field experiments have been carried out over Black Mountain [19], Ailsa Craig Island [20], Blashaval Hill [21], Askervein Hill [22], Kettles Hill [23], Hjardemal Escarpment [24], and so forth. Up to now, the Askervein Hill project [22] has been the most wellknown field campaign for validating numerical codes and evaluating turbulence models. In the past, several studies [3, 10, 25–28] reported successful validations of their numerical models chiefly over Askervein Hill. However, it must be mentioned that the steepness of Askervein Hill is below , and thus it may be characterized as a smooth terrain which presents almost twodimensional flow features [6]. Also, due to the gentle slope (), the Hjardemal site cannot truly be considered challenging, as the complex flow phenomena, such as flow separation, do not occur.
Recently, the Bolund field campaign [29–31] over Bolund Hill in Denmark has provided a new experimental dataset of the mean flow and the turbulence properties. In this work, to validate the LES methodology over a realistic complex terrain, we chose the challenging site of the Bolund. The hill is relatively smaller (12 m high) in scale than Askervein Hill, but it offers a challenging topography with steep slopes and a cliff (). The selection of Bolund Hill is justified by the threedimensional nature of the hill and by the flow separations induced by the cliff that make it more challenging to model than gentler terrains, such as Askervein Hill or Hjardemal Hill.
The first ever numerical and/or experimental (at laboratory scale) validation of the Bolund flows was initiated by Bechmann et al. [30] as a blind comparison of different microscale atmospheric flow models. They employed 57 models ranging from numerical to physical, including LES, RANS, and linearized models, as well as experimental models such as wind tunnel and water channel. From the blind comparison [30], they concluded that a RANS model with twoequation closure appeared to be the best model, due to the smallest simulation error for predicting the velocity speedup and turbulent kinetic energy (TKE), even compared to LES and experimental models. Such interesting feedback from the blind test is obviously encouraging for further investigation of the flow structures over the complicated Bolund site, especially by using LES model. After the blind test, there have been a number of attempts to simulate the Bolund flows by using numerical [6, 32–35] as well as experimental (windtunnel) [34, 36–38] models. However, it should be emphasized that all of the previous attempts, except the RANS simulations conducted by Prospathopoulos et al. [6], are limited to only onewinddirection results. To our knowledge, the present work is the first LESbased study presenting and analyzing the results from more than one inflow wind direction of the Bolund case. In addition, the previous studies have focused mostly on the mean flow analysis but have not focused on the turbulence statistics in enough detail. In fact, the information on turbulence prediction is one of the important parts for a complete model validation.
In the present work, LES are carried out to investigate the atmospheric flow over the Bolund Hill, and the LES model is validated against the full scale measurement [29, 30]. The LESbased results from two different wind directions, and , are reported and analyzed in detail for both the mean and the turbulence quantities. Following Bechmann et al. [30], the results are evaluated in quantitative manner, for example, by using simulation error. The results from the two different wind directions are presented together and discussed in order to learn the differences and/or similarities that emerge owing to a change in wind direction. In the final section, using the simulation error and the modelvalidation criteria, this paper evaluates and compares the performance of the LES model with other models previously employed over the same hill. Through the quantitative comparison, it is found that the present LES yields the most accurate prediction for the TKE, due to the smallest simulation error. The simulation error for the TKE increase is 19% smaller than in any other (experimental or numerical) models employed so far over the Bolund Hill.
2. The Bolund Experiment
Bolund Hill is a 12 m high, roughly 130 m long, and 75 m wide, and it is located on the north of DTU (Technical University of Denmark) in the city of Roskilde (Denmark). It is surrounded by sea, with the exception of the east side, which is connected to the mainland by a narrow isthmus. The field experiment was performed over a period of three months in 20072008, with the intention of providing a new fullscale dataset for validating atmospheric flow models in wind energy. Figure 1 shows the Bolund orography data with the locations of measuring masts (M0 to M9) installed on the hill during the field campaign. In this field campaign, 23 sonic and 12 cup anemometers were simultaneously employed to measure the wind velocity and variances from the 10 different masts. As it can be seen in Figure 1, the measuring masts (referred to hereafter as M1 to M8) were placed on two lines, LineA and LineB. In addition, two more masts, M0 and M9, were utilized to provide the inflow wind conditions for westerly and easterly winds, respectively. At each mast, the data were recorded for four different wind directions: , , , and . According to Berg et al. [29], the mean wind speed during the measurements was around 10 m/s, which leads to a Reynolds number of , and, thus, the flow can be considered to be independent of the . Also, the field data were collected only for the neutral stability condition. In addition to offering the latest fullscale database, the Bolund field experiment also provides an interesting analysis from the blind comparison performed by Bechmann et al. [30], which allows one to analyze the validation in great detail. More detailed information about the field experiment can be found in [29–31].
3. LES Methodology
3.1. LES Model and the Code
In LES, the turbulent eddies larger than the grid size, which are the most energetic ones and responsible for the majority of turbulent transport, are resolved directly, whereas the eddies smaller than the grid size are more isotropic and are modelled using a subgrid scale (SGS) model. The governing equations for LES can be obtained by filtering the original continuity and NavierStokes equations. The filtered governing equations for incompressible flow can be written aswhere and are the gridfiltered values of the velocity and the pressure, respectively, is the kinematic viscosity, is the fluid density, and is the subgrid scale stress defined as follows:Since the smaller eddies are not resolved in the simulations, their effects are modelled using the oneequation eddy viscosity SGS model proposed by Yoshizawa [39]:where , is the strain rate tensor, and , , where is the volume of the computational cell. and .
In this work, the unstructured finitevolume method based code, OpenFOAM [40], is employed to solve the governing equations. The LES calculations are carried out by using our own inhouse OpenFOAMbased LES solver [18, 33], recently implemented into OpenFOAM. The solver uses the classical fourthorder Runge–Kutta (RK4) scheme for time integration and a projection method where the velocity field is corrected using the pressure gradient inbetween the RK4 substeps. The pressure is then obtained by solving the Poisson equation. More information on the LES solver, including its validation, can be found in [18]. The oneequation model (3) is time integrated numerically using the secondorder backward implicit method, and it is solved after the fourth (last) RK4 substep. All the governing equations are discretized in space using the central differencing scheme. Following the recommendations of Berg et al. [29] and Bechmann et al. [30], the present simulations do not include the effects of the Coriolis force, and the atmosphere was considered to be neutrally stratified.
3.2. Numerical Domain and Grid Resolution
Figure 2 shows the computational setup of the simulations, including an explanation of the recycling approach, and Figure 3 shows the numerical model of the hill with a surface grid. As discussed previously, the simulations use the recycling technique to generate the timedependent inlet boundary conditions at the inflow plane (see Figure 2). For this reason, the length of the computational domain must be extended on the upstream side in order to ensure the fully developed flow before the reference mast M0. In the streamwise direction, the domain ranges from m to 400 m. The location of the outlet boundary was followed from the previous RANS study [6]. The vertical height of the domain is set to 120 m, and the boundarylayer depth is also 120 m, which is approximately 10 times higher than the height of Bolund Hill , that is, . The computational domain is also extended from −205 m to 185 m in the spanwise direction to avoid the periodic effects on the hill part, as the periodic boundary conditions are applied in that direction. The center of the hill is located at . In this way, the LES computational domain is of the size in , , and directions, respectively.
The computational domains are then discretized using a blockstructured hexahedron mesh with finer (order of cm) grid resolution on the hill surface. To select the grid resolutions in the present case, some guidelines were followed from the previous work [6, 35] dealing with the same hill. The computational domains were meshed using the commercial grid generator called Pointwise [41]. In both LES cases, the grid created was nonisotropic with the finest grid spacing of m and m being kept on the vertical escarpment (see Figure 3). On the rest of the hill part, the streamwise resolution varies from 0.31 m to 0.85 m and the spanwise from 0.40 m to 0.70 m depending on the local steepness of the hill. Outside the hill part, grid sizes in both horizontal directions are relaxed up to 4 m near the domain boundaries. The vertical grid spacing also varies gradually from 0.15 m to 4.97 m at the top boundary. In this way, the entire computational grid consists of hexahedron cells in each LES calculation.
3.3. Boundary Conditions
Reproducing the correct inflow velocity conditions is a challenge in LES. The inflow conditions must contain a realistic distribution of turbulence level throughout the boundary layer. Previously in several LES calculations [10, 25–27, 35], a separate precursor simulation over a flat and homogeneous surface has been done to generate inlet conditions for complex terrain simulations. This means that the simulation must be done in two parts (for two domains or grids), and then the data from the first (precursor) simulation is used in the second (main) simulation, that is, over complex terrain. As a consequence, the entire process (precursor and main simulations) takes a longer clock time for finalizing the simulation. Further, it requires large memory storage to store the data at each time step. Even more importantly, if one wants to make any changes in the main computational domain or grid, then the precursor simulation needs to be run again after introducing the new parameters. However, on the positive side, the precursor method may be the best option when several simulations are planned to be run using the same inflow conditions, as one can utilize the stored data in all those simulations. For the present case, the fully developed inflow boundary condition is obtained using the recycling (or mapping) method; see, for example, Chaudhari [18]. With the recycling technique, the precursor simulation is not required, and thus the entire simulation is performed in a single step; however, the user have to bear the extra CPU cost of the recycling section every time.
Using the recycling technique, the flow variables are sampled on a spanwise plane (recycling plane in Figure 2), which is sufficiently downstream of the inflow plane. Then, the sampled data are recycled back to the inflow plane at each time step, as explained in Figure 2. In this way, a recycling section is created between the inflow and the recycling planes, in which the flow becomes fully developed. Recently, Chaudhari [18] tested the suitability of the method under different recycling lengths , that is, the distance between the inflow and the recycling planes. It was concluded that should be at least three times the boundarylayer depth, that is, , in order to avoid artificial turbulence structures within the recycling section [18]. Hence, in this case, was used. Furthermore, the average velocity flux is fixed at the inflow boundary in order to maintain the same amount of volumetric flow throughout the simulation. The static pressure is fixed to a constant value on the outlet plane and the zerogradient (Neumann) [40] boundary condition is used for the rest of the flow variables. The slip boundary condition is used for all the flow variables at the top boundary. The periodic boundary conditions are set in the crosswind direction. A roughnesslength based wall function is utilized to determine the ground surface fluxes on the lower boundary. The implementation and the validation of the wall function can be found from [18, 33]. Following Berg et al. [29], two types of the roughnesslength values, m and m, are used in the simulations for the water and the hill surfaces, respectively.
In this work, two LES calculations for the two inflow wind directions ( and ) are carried out. The simulations are run until the physical time of s, in which all the flow quantities are timeaveraged over the last 600 s for calculating the flow statistics. The wallclock time for each of the calculations was about 40 days using 512 processors for simulating the physical time of 2700 s. The calculations were carried out using Finland’s most powerful supercomputer called Sisu (Cray XC30). Due to the high CPU requirement, a grid independent study of the present case is subject to large computational resources. We therefore rely on the agreement of the results (see Table 1) between the present LES and the previous numerical and experimental cases in order to evaluate the quality and the reliability of the present simulations. For the same reason (high CPU resources), the current study is based on the 10 min timeaveraged results. However, the averaging time taken here is more than 6.5 advection times along the entire computational domain, or 32 advection times along the hill length, and importantly it is 3 times as long as that used in the previous LES by Diebold et al. [35].

4. Results and Discussions
We begin this section by discussing the results at the reference mast M0. After this, the LES results at the rest of the masts (M1 to M8) are presented and discussed. The results from both wind directions are discussed simultaneously in order to illuminate the differences between the results from the two directions. Finally, simulation error is calculated for both wind directions, and the errors are compared against each other as well as against the literature values.
4.1. Results from the Reference Location
Figure 4 compares the LES profiles of the mean velocity magnitude and the turbulent kinetic energy (TKE) with the field measurements at the reference location M0. The results are normalized by the frictional velocity at M0. Here, the mean flow profile is very well reproduced by the recycling approach, as the simulated profile agrees well with the measurements at all heights. The TKE is also well predicted by the simulation at 5 m, but it is underestimated by about 25% at the height of 12 m. The total (modelled and resolved) TKE also underestimates the measured TKE at 12 m. The measurements show almost constant level of the TKE at M0, but in the simulation, the TKE decreases with height. Earlier, the LES runs submitted in the blind comparison also did not predict the constant TKE profile. This indicates that the deficiency could be very typical for LES modelling in general and thus occurs in the preset model, too. Since the modelled part of the TKE is significant in the area very close to the surface, we present the total TKE in all the results (Figures 6–13) hereafter. The error bars in the measurements show the standard deviations of the respective mean values. The standard deviation values were provided together with other statistics in the field data [29, 30].
(a)
(b)
4.2. Results on the Hill
Here, we first visualize the LES results on the instantaneous and mean flow structures around the hill. Figure 5 depicts the instantaneous velocity fields as well as the 10min timeaveraged ones. It can clearly be seen that the local flow field is strongly affected by the hill. In particular, the cliff is noted to enhance the turbulent fluctuations. The complexity occurring in the flow due to the cliff is also illustrated by the streamlines which highlight the mean (10min) flow patterns around the hill (Figure 5(b)). As expected, the flow is found to be separated at several places just after the cliff as well as on the lee sides. Particularly, the cliff is forming the small zone of the reversed flow, near the ground and along the entire cliff. Masts M2 and M6 are highly influenced with this reversed flow. However, the previous studies [6, 35] including those wind tunnels [34, 36–38] did not predict this flow separation nearby the cliff.
(a)
(b)
(a) Mean velocity
(b) TKE
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
Next, the results comparison between the simulation and the field data at different heights of all 10 masts (total 3738 locations) is discussed. In order to see the correlation of the results, the scatter plots of all the available field data against the LES predictions for the mean velocity magnitude and the TKE are shown in Figures 6(a) and 6(b), respectively. In the plots, an error interval of is also drawn for reference purpose. In Figure 6(a), it can be observed that 82% of the LESpredicted mean velocity in the case of is within the interval of . The velocity prediction is almost the same, 81%, in the case of . However, a very small deviation (about 9%) is observed between the two cases when examining the turbulence part. For the TKE, 92% and 83% of the predictions appear within the interval () in the cases of and , respectively. When one considers the interval of , the TKE prediction is still good, 80%, in both wind directions.
To describe the discrepancies between the simulations and the measurements in detail, more scatter plots on the velocity and TKE predictions are presented in Figures 7 and 8, respectively. In these figures, the results are sorted with respect to the measurement height above the ground level (Figures 7(a) and 8(a)) and with respect to the masts (Figures 7(b) and 8(b)). The locations of the masts can be found from Figure 1. According to Figure 7(a), the predicted mean velocity is scattered for the height below 2 m, reaching a maximum error of almost 128% at M2 in the case of . On the other hand, the largest error (at 2 m) is smaller by about 35% in the case of . Generally, in both cases at , masts M2, M6, and M8 have larger errors than the rest of the masts, as can also be seen in Figure 7(b). The large errors are due to the mast locations in which the flows are highly turbulent with flow separations (see Figure 5(b)). In both cases, the maximum error is reduced by about 100% within a distance of 3 m, from m to 5 m. Above the height of 5 m, the prediction is greatly improved, as the maximum error is only 11%, while the majority of the error is below 5%.
The overall prediction of the TKE is also scattered in both cases ( and ), although the level of maximum error is reduced by about 90% (Figure 8(a)) as compared to that of the velocity. The largest discrepancy is again found at M2, which is 70% (at 3.6 m) in the case of and 65% (at 2 m) in the case of . However, unlike the velocity, the scatter is relatively independent of the mast locations and measurement heights (Figure 8(b)). Despite having less error (as compared to velocity), the TKE prediction does not seem to be improved with the height from 2 m to 5 m (Figure 8(a)). Overall, the velocity prediction is better in the case of , whereas the TKE prediction is better in the case of . More quantitative comparison between the twowinddirection results is presented in Section 4.3.
Following Bechmann et al. [30], the numerical results are further analyzed in terms of velocity speedup and TKE increase . The speedup and the TKE increase are defined as follows [30]: where the subscript 0 indicates the values ( or ) from the reference mast M0. Vertical profiles of the speedup and the TKE increase at all the masts are shown in Figures 9 and 10, respectively.
Here, it can be observed that the LES prediction of the speedup is generally acceptable (error less than 10%; Figure 11(a)) at most (70%) of the places for both wind directions. Nevertheless, for the wind, the simulation overestimates the measurements at the lowest heights of M2, M5, and M6, while it underestimates them at of M8. The maximum speedup error at the lowest height (2 m) is 53% and it decreases to 25% at 5 m in the case of wind. Concerning the lowest heights of M2, M5, and M6, the speedup is predicted better in the case of , as the maximum error is only about 34%, which is at 2 m height of M2 (see Figure 11(b)). We note that the vertical profiles of the predicted speedup at M5 and M6 (below 5 m) differ significantly between the two directions. At the rest of the masts, there is no such a significant difference with respect to the wind direction. This means the flows in those regions (M5 and M6) are highly sensitive to the wind direction. According to the measurements, the speedup prediction in those areas is better in the case of .
The large discrepancies between the simulations and the measurements at M2, M5, M6, and M8 are justified by the complex nature of the flows locally, that is, highly turbulent and reversing with several flow separations generated by the steepness of the hill. Particularly, masts M2 and M6 were placed about 10 m away from the cliff and were strongly influenced by the very small flow separation (due to the cliff) along the edge of the cliff, as it can be seen in Figure 5(b). This particular separation is indeed much challenging to be even captured by numerical simulations because of its small size and location (near the cliff). Consequently, none of the previous attempts, including the windtunnel tests, have predicted this flow separation. Masts M8 and M4 were located on the lee side of the hill for and winds, respectively, and were also continuously influenced by flow separations. The simulation does predict most of the flow separations there, but the depth of the flow separation is somewhat compromised, although more at M8 than at M4. Previous studies [6, 30, 35–38] also reported large discrepancies at several masts, except for M1 and M7.
Due to the flow separations, the turbulence level is significantly increased in a few locations. As a result, a clear peak in the TKE increase profile is visible at M2 and M6 for both winds and at M5 for the wind (see Figures 10(a) and 10(b)). Importantly, the simulation appears to accurately predict the extreme TKE vertical gradient, especially in the case of , including the lower heights (Figure 10(a)), even though the velocity speedup is highly overestimated at the same locations (i.e., M2 and M6). With the exception of M2 and M6, the TKE increase is also reproduced satisfactorily in the case of direction. At 2 m height of M2 (for ), the simulation shows the largest error of 140% for predicting the TKE increase (Figure 11(b)). Most of the discrepancies between the simulation and the measurements, visible in Figures 9 and 10, and especially at M2, M6, and M8, may also be justified by the uncertainties in the measured statistics. The high standard deviations at the lowest heights of M2 and M6 indicate high variabilities of the measured data at those places.
4.3. Simulation Error and Validation Metric
To evaluate the performance of the LES model in a quantitative way, the simulation (modelling or prediction) error from the LES results is calculated by considering the field data as a reference. Originally, Bechmann et al. [30] used the following equations to calculate the simulation errors for speedup and TKE increase. The error of the simulated speedup, , is defined as follows:and the error of the simulated TKE, , iswhere is the speedup, , subscript zero denotes that the value is taken at the reference locations, and subscripts and denote the measured and simulated values, respectively. Figure 11 shows the simulation errors for the speedup and the TKE increase at each individual mast as well as at each possible height for the both wind directions. In the same way, Figure 12 displays the simulation errors, but they are averaged over the heights of each mast. The individual error from each height and/or location in Figure 11 may also be used for comparison with the results shown in the previous figures.
In the following, we turn our focus on the averaged error to discuss the results. According to Figure 12, the speedup error is less than 10% for half (50%) of the masts in both cases. The maximum error in the case of is 27%, which is found at M6. On the other hand, the maximum error for the wind is found at M5 and it is 18%. For the TKE increase, the error () is much higher than the speedup error (), reaching 59% (at M2) and 80% (at M6) in the cases of and winds, respectively. The high error for the TKE increase is obviously because the secondorder statistic, for example, TKE or the Reynolds stress, is challenging to predict with enough accuracy as compared to mean flow (firstorder statistic) in general. An interesting fact noted here is that the TKE increase error is almost the same at M2 for both winds, but a significant difference (75%) is observed at M6 due to the change in wind direction. Despite the difficulty of predicting the TKE compared to the mean wind, the error is found to be smaller than at M6 and M8 in the case of and at M5 in the case of . The high error at certain masts, such as M2 and M6, is due to their locations (near the cliff), where the aforementioned extremely complex flow behaviors occur. Indeed, virtually all previously reported models, including those wind tunnels, have failed to predict the flow accurately enough at masts M2 and M6. When excluding both M2 and M6, the averaged speedup error changes from 10.3% to 10.2%, which is insignificant. However, a high reduction from 24.1% to 13.5% is observed for the turbulence quantity. In addition, among all the masts, the standard deviation of the measured TKE is also at its highest at M2 and M6, reaching up to 24% in the case of , as could be expected from the error bars in Figures 9 and 10. The standard deviation or any other uncertainties of the measurements are, of course, not taken into account in the calculations of the simulation error.
Next, Figure 13 shows the (averaged) simulation error at three different heights. Here, the error is averaged over the measurement locations having the same height. For the wind direction, it can be seen that the errors for the speedup and for the TKE increase at the height of 2 m are nearly the same (21%). This indicates that, at the lowest height, the model performs equally for both the speedup and the TKE increase. Above 2 m, the speedup error decreases rapidly down to 2.5% at 9 m. However, this does not happen with the TKE increase. The TKE increase error remains in the same order from 2 m to 5 m, and it decreases nominally, only by 4%, at 9 m.
For the wind direction, the speedup error also decreases with increasing height, but not in the same trend as in the case. It reduces from 13.5% at 2 m to 5.5% at 9 m. The TKE error is much higher, about 41% at 2 m (Figure 13). At this height, the TKE error is almost double compared to that of the case. The high error is mostly due to the poor predictions at M2 and M6. At 5 m, the TKE error reduces greatly to 21%, about the same level as in the case; however, it does not reduce further at 9 m.
It should be emphasized that the measurements of the TKE at 2 m and 9 m heights of the reference location (M0) were missing in the field data (see Figure 4(b)). Therefore, by following [30], a constant value of the TKE is considered at all three heights (2 m, 5 m, and 9 m) to calculate the measured TKE increase. This assumption (constant TKE with height) can highly bias the errors especially at 2 m, as the LES results follow the heightdependent TKE values from the simulation. Moreover, it should be remembered that the flow within the surface layer, for example, m in the Bolund case, may be of little interest for a windfarm developer, as one is considering the wind condition within the rotor diameter. With the average over both the wind directions, the speedup error at 5 m is 9.3% and is reduced to 3.9% at 9 m. The error is not completely negligible, but the prediction could be considered accurate (less than 10% error) [30, 35].
Table 1 shows the mean errors from the two wind directions, and it compares the errors from the present case with those reported previously in the literature. The mean simulation error of the present LES is 10.3% and 24.1% for speedup and TKE increase, respectively. The LES prediction of the speedup is among the best results, as the mean speedup error from our case is only 0.7% higher than that of the best results [35]. Concerning the turbulence aspect, the present results show far smaller error in comparison with all previous results. Indeed, all prior models, except a RANS twoequation model in [30], are suffering a high TKE increase error. A couple of studies [6, 35] did not report the TKE increase error at all. The TKE increase error of 24.1% from this case is about 19% smaller than that of any other model. To our knowledge, this is thus the best prediction of the TKE over the Bolund Hill among all models employed over the hill so far.
Finally, the present LES model is evaluated here using three different validation metrics. First, the validation metric called the fraction of two (FAC2) shows the fraction of data that is within the interval:where is the predicted and is the observed value. For an ideal model prediction, one would have FAC2 = 1 [43]. Second, the Fractional Bias (FB) can be defined aswhere the angular bracket represents the average value over the data samples. FB indicates the relative amount of bias in the modelled results such that FB means no bias at all and FB means systematic underestimation [44]. Finally, a third metric, known as the Normalized Mean Square Error (NMSE), is defined as follows:NMSE is sensitive to large disagreements, in particular those occurring on the highest values [44]. The validationmetric values from the present results for both wind directions are listed in Table 2. According to the table, the present LES model passes all three validation tests, using the acceptance criteria for the model evaluation (in a rural area) proposed by Hanna and Chang [42]. In particular, the FAC2 value obtained by the present results is in general a notably high value in the evaluation of atmospheric flow models.

5. Conclusions
This paper has focused on simulating the atmospheric flows over natural complex terrain and on validating the LES model against fullscale measurements. Previously, several researchers utilized LESbased numerical models to validate their simulations, especially over Askervein Hill, which is characterized by a smooth terrain. However, in this paper, we have selected much challenging complex terrain, the Bolund Hill, for validating the LES methodology. Unlike most of the previous studies on the Bolund case, we have carried out LES calculations for two different wind directions. The results on the mean and turbulence (e.g., TKE) quantities from the two wind directions have been analyzed equally and reported in detail herein. The results from both wind directions have been presented simultaneously and discussed quantitatively in order to learn the differences and the similarities in the results due to the change in wind direction. Using the simulation error and the modelvalidation criteria, the paper has evaluated and compared the performance of the LES model with other models previously employed over the hill.
The recycling method reproduced the correct mean velocity profile at the reference location (M0) but showed some difficulties in the prediction of especially constant TKE profile. Unlike the field data, the LESpredicted TKE profile decreases gradually above m; however, it reproduced the correct (measured) value at 5 m. Bechmann et al. [30] mentioned that all the LES models that participated in the blind test did not reproduce the correct inflow (velocity and TKE) profiles. This indicates that the reproduction of the correct inflow profiles, and especially the constant TKE, in the Bolund case is a global challenge to LES modelling in general. However, in comparison with other prior LES models [30, 35], the prediction of the inflow conditions is better in this study thanks to the recycling technique.
Generally speaking, the overall predictions of the speedup and the TKE increase over the hill are in close agreement with the measurements at all the positions, with the averaged difference of 10.2% for the speedup and 13.5% for the TKE increase, excluding M2 and M6. Due to the complex locations of M2 and M6, the simulation suffers from a maximum error of 140% to predict the measured values at 2 m height. Indeed, the measurements also show high standard deviations (up to 24%) at the same places, and this indicates that the measured data have high variabilities at those places (M2 and M6). Above the height of 5 m, the LES prediction of the mean flow is convincingly accurate, as the maximum difference between the simulations and the measurements is only 11%. In particular, the successful prediction of the reversed flow at the escarpment edge is believed to be an excellent achievement of this particular LES model, as no other previous numerical simulations had predicted the separation. For the turbulence part, the simulation reports the best prediction of the TKE, with the averaged error of 24.1%, which is 19% smaller than any other models, including the windtunnel experiments, reported in the past. Hence, it can be concluded that the present LES model, using the recycling technique on the inlet side, is able to reproduce the complex turbulent structures of the local wind flow over a complicated terrain such as Bolund Hill. Thus, it is possible to employ the same LES methodology to analyze the wind structures over other real terrains.
Currently, this work is being extended with the inclusion of the full depth of the atmospheric boundary layer with capping inversion and the Earth’s rotation into the simulations. In the near future, the aim is to study the influence of capping inversion and the Coriolis force on the local atmospheric flows over this challenging Bolund site. The new results will then be compared with the existing results presented in this paper.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The computational resources for this work were provided by CSCFinnish IT Center for Science Ltd. in Espoo, Finland. This research has been a part of the RENEWTECH Project (2011–2014), which aimed at the development of wind power technology and business in SouthEast Finland. The work was funded by the European Regional Development Fund (ERDF). The authors acknowledge with thanks the support of the TU1304 action “WINERCOST.” The first author, A. Chaudhari, would like also to thank Lappeenranta University of Technology (LUT) for funding this research as a part of his postdoctoral research at LUT.
References
 A. Maurizi, J. M. L. M. Palma, and F. A. Castro, “Numerical simulation of the atmospheric flow in a mountainous region of the North of Portugal,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 74–76, pp. 219–228, 1998. View at: Publisher Site  Google Scholar
 H. G. Kim, V. C. Patel, and C. M. Lee, “Numerical simulation of wind flow over hilly terrain,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 87, no. 1, pp. 45–60, 2000. View at: Publisher Site  Google Scholar
 F. A. Castro, J. M. L. M. Palma, and A. S. Lopes, “Simulation of the Askervein flow. Part 1: Reynolds averaged NavierStokes equations ($k\text{}\in $ turbulence model),” BoundaryLayer Meteorol, vol. 107, no. 3, pp. 501–530, 2003. View at: Publisher Site  Google Scholar
 L. M. S. Paiva, G. C. R. Bodstein, and W. F. Menezes, “Numerical simulation of atmospheric boundary layer flow over isolated and vegetated hills using RAMS,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 97, no. 910, pp. 439–454, 2009. View at: Publisher Site  Google Scholar
 A. El Kasmi and C. Masson, “Turbulence modeling of atmospheric boundary layer flow over complex terrain: a comparison of models at wind tunnel and full scale,” Wind Energy, vol. 13, no. 8, pp. 689–704, 2010. View at: Publisher Site  Google Scholar
 J. M. Prospathopoulos, E. S. Politis, and P. K. Chaviaropoulos, “Application of a 3D RANS solver on the complex hill of Bolund and assessment of the wind flow predictions,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 107108, pp. 149–159, 2012. View at: Publisher Site  Google Scholar
 B. Blocken, A. van der Hout, J. Dekker, and O. Weiler, “CFD simulation of wind flow over natural complex terrain: case study with validation by field measurements for Ria de Ferrol, Galicia, Spain,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 147, pp. 43–57, 2015. View at: Publisher Site  Google Scholar
 L. Li, P. W. Chan, L. Zhang, and F. Hu, “Numerical simulation of a lee wave case over threedimensional mountainous terrain under strong wind condition,” Advances in Meteorology, vol. 2013, Article ID 304321, 13 pages, 2013. View at: Publisher Site  Google Scholar
 Y. Zheng, Y. Miao, S. Liu, B. Chen, H. Zheng, and S. Wang, “Simulating flow and dispersion by using WRFCFD coupled model in a builtup area of Shenyang, China,” Advances in Meteorology, vol. 2015, Article ID 528618, 15 pages, 2015. View at: Publisher Site  Google Scholar
 A. Bechmann and N. N. Sørensen, “Hybrid RANS/LES method for wind flow over complex terrain,” Wind Energy, vol. 13, no. 1, pp. 36–50, 2010. View at: Publisher Site  Google Scholar
 S. B. Pope, Turbulent Flows, Cambridge University Press, Cambridge, UK, 2nd edition, 2000. View at: Publisher Site  MathSciNet
 J. Sumner, C. S. Watters, and C. Masson, “CFD in wind energy: the virtual, multiscale wind tunnel,” Energies, vol. 3, no. 5, pp. 989–1013, 2010. View at: Publisher Site  Google Scholar
 B. Zhou and F. K. Chow, “Turbulence modeling for the stable atmospheric boundary layer and implications for wind energy,” Flow, Turbulence and Combustion, vol. 88, no. 12, pp. 255–277, 2012. View at: Publisher Site  Google Scholar
 T. Allen and A. R. Brown, “Largeeddy simulation of turbulent separated flow over rough hills,” BoundaryLayer Meteorology, vol. 102, no. 2, pp. 177–198, 2002. View at: Publisher Site  Google Scholar
 T. Tamura, S. Cao, and A. Okuno, “LES study of turbulent boundary layer over a smooth and a rough 2D hill model,” Flow, Turbulence and Combustion, vol. 79, no. 4, pp. 405–432, 2007. View at: Publisher Site  Google Scholar
 A. Chaudhari, A. Hellsten, O. Agafonova, and J. Hämäläinen, “Large eddy simulation of boundarylayer flows over twodimensional hills,” in Progress in Industrial Mathematics at ECMI 2012, M. Fontes, M. Günther, and N. Marheineke, Eds., vol. 19 of Mathematics in Industry, pp. 211–218, Springer, 2014. View at: Google Scholar
 O. Agafonova, A. Koivuniemi, A. Chaudhari, and J. Hämäläinen, “Limits of WAsP modeling in comparison with CFD for wind flow over twodimensional hills,” in Proceedings of the European Wind Energy Association Conference (EWEA '14), Barcelona, Spain, 2014. View at: Google Scholar
 A. Chaudhari, Largeeddy simulation of wind flows over complex terrains for wind energy applications [Ph.D. thesis], Lappeenranta University of Technology, 2014.
 E. F. Bradley, “An experimental study of the profiles of wind speed, shearing stress and turbulence at the crest of a large hill,” Quarterly Journal of the Royal Meteorological Society, vol. 106, no. 447, pp. 101–123, 1980. View at: Publisher Site  Google Scholar
 G. J. Jenkins, P. J. Mason, W. H. Moores, and R. I. Sykes, “Measurements of the flow structure around Ailsa Craig, a steep threedimensional, isolated hill,” Quarterly Journalof the Royal Meteorological Society, vol. 107, no. 454, pp. 833–851, 1981. View at: Publisher Site  Google Scholar
 P. J. Mason and J. C. King, “Measurements and predictions of flow and turbulence over an isolated hill of moderate slope,” Quarterly Journal of the Royal Meteorological Society, vol. 111, no. 468, pp. 617–640, 1985. View at: Publisher Site  Google Scholar
 P. A. Taylor and H. W. Teunissen, “The Askervein hill project: overview and background data,” BoundaryLayer Meteorology, vol. 39, no. 12, pp. 15–39, 1987. View at: Publisher Site  Google Scholar
 J. R. Salmon, H. W. Teunissen, R. E. Mickle, and P. A. Taylor, “The Kettles hill project: field observations, windtunnel simulations and numerical model predictions for flow over a low hill,” BoundaryLayer Meteorology, vol. 43, no. 4, pp. 309–343, 1988. View at: Publisher Site  Google Scholar
 S. Emeis, M. S. Courtney, J. Højstrup, and N. O. Jensen, “Hjardemaal experiment data report,” Tech. Rep., Risø National Laboratory, Roskilde, Denmark, 1993. View at: Google Scholar
 A. Bechmann, Largeeddy simulation of atmospheric flow over complex terrain [RisøPhD thesis], RisøPhD28(EN), 2006.
 A. Silva Lopes, J. M. L. M. Palma, and F. A. Castro, “Simulation of the Askervein flow. Part 2: largeeddy simulations,” BoundaryLayer Meteorology, vol. 125, no. 1, pp. 85–108, 2007. View at: Publisher Site  Google Scholar
 F. K. Chow and R. L. Street, “Evaluation of turbulence closure models for LargeEddy simulation over complex terrain: flow over Askervein Hill,” Journal of Applied Meteorology and Climatology, vol. 48, no. 5, pp. 1050–1065, 2009. View at: Publisher Site  Google Scholar
 J.C. Golaz, J. D. Doyle, and S. Wang, “Oneway nested largeeddy simulation over the Askervein hill,” Journal of Advances in Modeling Earth Systems, vol. 1, no. 3, p. 6, 2009. View at: Publisher Site  Google Scholar
 J. Berg, J. Mann, A. Bechmann, M. S. Courtney, and H. E. Jørgensen, “The Bolund experiment, part I: flow over a steep, threedimensional hill,” BoundaryLayer Meteorology, vol. 141, no. 2, pp. 219–243, 2011. View at: Publisher Site  Google Scholar
 A. Bechmann, N. N. Sørensen, J. Berg, J. Mann, and P.E. Réthoré, “The Bolund experiment, part II: blind comparison of microscale flow models,” BoundaryLayer Meteorology, vol. 141, no. 2, pp. 245–271, 2011. View at: Publisher Site  Google Scholar
 A. Bechmann, J. Berg, M. Courtney, H. Ejsing Jørgensen, J. Mann, and N. N. Sørensen, The Bolund Experiment: Overview and Background, Danmarks Tekniske Universitet, Risø National Laboratory for Sustainable Energy, 2009.
 A. Chaudhari, V. Vuorinen, O. Agafonova, A. Hellsten, and J. Hämäläinen, “Largeeddy simulation for atmospheric boundary layer flows over complex terrains with applications in wind energy,” in Proceedings of the 11th World Congress on Computational Mechanics (WCCM '14), 5th European Conference on Computational Mechanics (ECCM '14), and 6th European Conference on Computational Fluid Dynamics, (ECFD '14), pp. 5205–5216, Barcelona, Spain, 2014. View at: Google Scholar
 V. Vuorinen, A. Chaudhari, and J.P. Keskinen, “Largeeddy simulation in a complex hill terrain enabled by a compact fractional step OpenFOAM® solver,” Advances in Engineering Software, vol. 79, pp. 70–80, 2015. View at: Publisher Site  Google Scholar
 B. Conan, A. Chaudhari, S. Aubrun, J. van Beeck, J. Hämäläinen, and A. Hellsten, “Experimental and numerical modelling of flow over complex terrain: the Bolund hill,” BoundaryLayer Meteorology, vol. 158, no. 2, pp. 183–208, 2016. View at: Publisher Site  Google Scholar
 M. Diebold, C. Higgins, J. Fang, A. Bechmann, and M. B. Parlange, “Flow over hills: a largeeddy simulation of the Bolund case,” BoundaryLayer Meteorology, vol. 148, no. 1, pp. 177–194, 2013. View at: Publisher Site  Google Scholar
 B. Conan, Wind resource accessment in complex terrain by wind tunnel modelling [Ph.D. thesis], Université d’Orléans, Orléans, France, 2012.
 T. S. Yeow, A. Cuerva, B. Conan, and J. PérezÁlvarez, “Wind tunnel analysis of the detachment bubble on Bolund Island,” Journal of Physics: Conference Series, vol. 555, no. 1, Article ID 012021, 2014. View at: Publisher Site  Google Scholar
 T. S. Yeow, A. CuervaTejero, and J. PérezÁlvarez, “Reproducing the Bolund experiment in wind tunnel,” Wind Energy, vol. 18, no. 1, pp. 153–169, 2015. View at: Publisher Site  Google Scholar
 A. Yoshizawa, “Bridging between eddyviscositytype and secondorder models using a twoscale DIA,” in Proceedings of the 9th International Symposium on Turbulent Shear Flow, vol. 3, pp. 23.1.1–23.1.6, Kyoto, Japan, 1993. View at: Google Scholar
 OpenCFD, OpenFOAM User and Programmer's Guide, OpenFOAM Foundation, OpenCFD Ltd., 2013.
 Pointwise, “Pointwise: mesh & grid generation software,” http://www.pointwise.com/ View at: Google Scholar
 S. Hanna and J. Chang, “Acceptance criteria for urban dispersion model evaluation,” Meteorology and Atmospheric Physics, vol. 116, no. 34, pp. 133–146, 2012. View at: Publisher Site  Google Scholar
 M. Schatzmann, H. R. Olesen, and J. Franke, COST 732 Model Evaluation Case Studies: Approach and Results, University of Hamburg, 2010.
 A. H. Maya Milliez, “Evaluation of dispersion models for improvement and guidance for the use of localscale emergency response tools; the ‘michelstadt’ modelling exercise—cost action es1006,” in Proceedings of the 6th International Symposium on Computational Wind Engineering (CWE '14), Hamburg, Germany, 2014. View at: Google Scholar
Copyright
Copyright © 2016 Ashvinkumar Chaudhari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.