The Bolund Isolated Steep Hill (Bolund)

Status

June 2014

The Bolund benchmark [BSB+11] was revisited in Wakebench Phase 1. The results were presented at the Torque 2014 conference.

Highlights

The results are consistent with the original blind test with much better results at 5 m than at 2 m. Fig. 22 shows mean profiles along the B5 and A5 profiles where the \(k-\epsilon\) model shows superior performance than the \(k-\omega\) model due to a much finer grid around the escarpment zone.

../../_images/bolund_profiles.png

Fig. 24 Horizontal profiles of fractional speed-up ratio and turbulent kinetic energy along the A and AA transects. © Author(s) 2014. CC BY 3.0 License. Used with permission. [RGA+14]

December 2009

The Bolund blind comparison was organized by DTU in 2009. You can find input and validation data as well as results of the blind test in project website [BSB+11].

Highlights

In total 49 different simulations were submitted, composed of 3 physical models, 9 linearized numerical models and 37 CFD models (5 LES, 7 RANS 1-equation and 25 RANS 2-equation). The physical models predicted reasonably well the mean velocity profiles but under-predicted the turbulent kinetic energy. Linear models produced the worse results as they were not capable of reproducing the flow around the steep escarpment. RANS models provided the best results although the spread of the simulations was quite big, indicating user dependencies especially regarding mesh generation. LES-based models had problems but presented promising results with regard to turbulence modelling in the flow separation area just after the escarpment. The average error in the simulation of the mean wind speeds over all the measurement locations was 13-17% in the top 10 best models. This error dropped to 4-10% if only the measurements at 5 m were considered.

Scope and Objectives

The benchmark revisits the blind test of 2009 now allowing the participants to optimize their models to obtain the best match to the validation dataset.

The objectives are:

  • Test model fine-tuning strategies that will be applied in complex terrain sites.

  • Evaluate turbulence models in a test site with well defined boundary conditions.

Background

Bolund is a 12 m high, 130 m long and 75 m wide isolated hill situated to the North of RisøDTU in Roskilde Fjord, Denmark. It is surrounded by water in all directions except to the E, where a narrow isthmus leads to the mainland. The hill is characterized by a uniform roughness of 0.015 m and surrounded by water with a roughness length of 0.0003 m. An almost vertical escarpment in the prevailing W-SW sector ensures flow separation in the windward edge resulting in a complex flow field, quite challenging for flow models.

../../_images/bolund_topo.png

Fig. 25 Elevation map and instrument layout.

Measurement Campaign

The masts are positioned along two lines: A and B. Two additional masts (M0 and M9) were installed to measure the incoming undisturbed flow for westerly and easterly winds respectively. Mast M9 is placed in the coastline, where the roughness length is again 0.015 m. The masts are equipped with 23 sonic (Metek USA 1-Basic) and 12 cup anemometers (Risø Wind Sensor P2546) at heights between 2 and 15 m.

The Bolund experiment comprises a measurement campaign of three months between 2007 and 2008 [BBC+09] [BMB+11]. During the measurement campaign the absolute water level was monitored, which covered the isthmus most of the time. The campaign was designed for W-SW winds where the fetch ranges 4 to 7 km, ensuring undisturbed velocity profiles over water. The E sector is more difficult to characterize due to a more the heterogeneous land cover. The prevailing stability regime was neutral to slightly stable conditions (1/L < 0.04). The data was averaged over 30 min periods.

Previous Work

The original blind test after the experiment was conducted by DTU in 2009 [BSB+11]. Recent work on the Bolund hill include RANS simulations [PPC12], LES simulations [DHF+13] and wind tunnel experiments [YCTPA15] [CCA+16].

Input Data

The conditions for simulating the Bolund flow field in neutral conditions are:

  • Digitized map of the Bolund hill with 25 cm resolution. Water level is set to 0.75 m.

  • Roughness digitized map: hill with \(z_0 = 0.015 m\), water with \(z_0 = 0.0003 m\), coastal (X > 325 m) with \(z_0 = 0.015 m\).

  • Inlet profiles: Measured at M0 for westerly winds and M9 for easterly winds.

  • Coordinates of met masts along lines A (239º) and B (270º).

  • No heat flux, gravity \(g = 9.81 m s^{-2}\), Coriolis parameter :math:` f_c = 1e-4 s^{-1}`.

  • Obukhov length: \(L = \infty\).

Use dry air with a density \(\rho = 1.225 kg m^3\) and dynamic viscosity \(\mu = 1.73e-5 kg m^{-1}s^{-1}\)

Validation Data

The validation dataset is composed of mean flow and turbulence data from cup and sonic anemometers at 10 met masts. Ensemble averages of 10 min averaged samples within ±8º wind direction sector, with wind speeds between 5 and 12 m/s at 5 m level and under neutral conditions (\(|1/L| < 0.004 m^{-1}\)) at the upstream masts, were used to derive the validation datasets which consists on:

  • Fractional-Speedup-Ratio (FSR) and normalized added turbulent kinetic energy (TKE) with respect to the reference inlet position, at 2 and 5 m above ground level along mast lines A and B

  • FSR and TKE vertical profiles at mast positions.

Velocity and TKE values will be normalized with the upstream friction velocity at the reference mast as in Bechmann et al. (2011). The validation dataset includes mean and standard deviation statistics from the ensemble profiles.

Model Runs

The inlet profile can be based on neutral M-O log-law, defined by the following input parameters:

  • Run 1: \(WD = 270º\), \(z_0 = 0.0003 m\), \(TKE/u_*^2 = 5.8\), \(u_* = 0.4 ms^{-1}\)

  • Run 2: \(WD = 255º\), \(z_0 = 0.0003 m\), \(TKE/u_*^2 = 5.8\), \(u_* = 0.4 ms^{-1}\)

  • Run 3: \(WD = 239º\), \(z_0 = 0.0003 m\), \(TKE/u_*^2 = 5.8\), \(u_* = 0.4 ms^{-1}\)

  • Run 4: \(WD = 90º\), \(z_0 = 0.015 m\), \(TKE/u_*^2 = 5.8\), \(u_* = 0.5 ms^{-1}\)

or by best fit to the measured inlet profiles (at M0 for runs 1,2 and 3 and M9 for run 4) if the participant considers that this can improve the results. The computational domain must extend at least to X = ±400 m in order to include the coastline to the East and make sure that the hill wake is completely covered. The origin of the coordinate system should be placed at M3 position with X pointing East, Y pointing North and Z pointing up.

Output Data

The simulated validation profiles consist on horizontal profiles along lines A and B at 2 and 5 m height and vertical profiles at mast positions of velocity components (U,V,W), turbulence kinetic energy (tke), dissipation rate (tdr), friction velocity (ust) and kinematic momentum fluxes (uu, vv, ww). The profiles should traverse the simulated domain from boundary to boundary. Hence, the required outputs are, in this order: X(m), Y(m), Z(m), U(m/s), V(m/s), W(m/s), tke(m2/s2), tdr(m2/s3), us(m/s), uu(m2/s2), vv(m2/s2), ww(m2/s2).

Use the file naming and format convention described in the Windbench user’s guide with profID = prof#, where # = [M0,M1,M2,M3,M5,M6,M7,M8,M9,A2,A5,B2,B5], i.e. 13 output files per user and model run. Additionally, for those users that participated in the blind test of 2009, please provide the output files that were obtained at that time. This will allow an assessment of the added value of onsite measurements for model tuning. Please follow the same format described before but with a BenchmarkID = Bolund_blind2009 to differentiate between the two sets of simulations.

Remarks

In order to evaluate the added value of model fine-tuning it is important that you describe how this is performed. Please report on the deviations with respect to default settings (those of the blind test). There are no guidelines on the definition of the computational mesh since this can have an important influence in the fine-tuning aspects of the model. Please describe how you integrate grid dependency in the evaluation process.

References

BSB+11(1,2,3,4)

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. Boundary-Layer Meteorol, 141(2):245, August 2011. URL: https://doi.org/10.1007/s10546-011-9637-x (visited on 2020-07-23), doi:10.1007/s10546-011-9637-x.

BBC+09(1,2)

Andreas Bechmann, Jacob Berg, Michael Courtney, Hans Ejsing Jørgensen, Jakob Mann, and Niels N. Sørensen. The Bolund Experiment: Overview and Background. Technical Report, Danmarks Tekniske Universitet, 2009. URL: https://orbit.dtu.dk/en/publications/the-bolund-experiment-overview-and-background (visited on 2020-07-23).

BMB+11(1,2)

J. Berg, J. Mann, A. Bechmann, M. S. Courtney, and H. E. Jørgensen. The Bolund Experiment, Part I: Flow Over a Steep, Three-Dimensional Hill. Boundary-Layer Meteorol, 141(2):219, July 2011. URL: https://doi.org/10.1007/s10546-011-9636-y (visited on 2020-07-23), doi:10.1007/s10546-011-9636-y.

CCA+16

Boris Conan, Ashvinkumar Chaudhari, Sandrine Aubrun, Jeroen van Beeck, Jari Hämäläinen, and Antti Hellsten. Experimental and Numerical Modelling of Flow over Complex Terrain: The Bolund Hill. Boundary-Layer Meteorol, 158(2):183–208, February 2016. URL: https://doi.org/10.1007/s10546-015-0082-0 (visited on 2020-07-23), doi:10.1007/s10546-015-0082-0.

DHF+13

Marc Diebold, Chad Higgins, Jiannong Fang, Andreas Bechmann, and Marc B. Parlange. Flow over Hills: A Large-Eddy Simulation of the Bolund Case. Boundary-Layer Meteorol, 148(1):177–194, July 2013. URL: https://doi.org/10.1007/s10546-013-9807-0 (visited on 2020-07-23), doi:10.1007/s10546-013-9807-0.

PPC12

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, 107-108:149–159, August 2012. URL: http://www.sciencedirect.com/science/article/pii/S0167610512001134 (visited on 2020-07-23), doi:10.1016/j.jweia.2012.04.011.

RGA+14(1,2)

Javier Sanz Rodrigo, Pawel Gancarski, Roberto Chavez Arroyo, Patrick Moriarty, Matthew Chuchfield, Jonathan W. Naughton, Kurt S. Hansen, Ewan Machefaux, Tilman Koblitz, Eoghan Maguire, Francesco Castellani, Ludovico Terzi, Simon-Philippe Breton, Yuko Ueda, John Prospathopoulos, Gregory S. Oxley, Carlos Peralta, Xiadong Zhang, and Björn Witha. IEA-Task 31 WAKEBENCH: Towards a protocol for wind farm flow model evaluation. Part 1: Flow-over-terrain models. J. Phys.: Conf. Ser., 524:012105, June 2014. Publisher: IOP Publishing. URL: https://doi.org/10.1088%2F1742-6596%2F524%2F1%2F012105 (visited on 2020-10-14), doi:10.1088/1742-6596/524/1/012105.

RM14

Javier Sanz Rodrigo and Patrick Moriarty. IEA-Task 31 WAKEBENCH: Towards a protocol for wind farm flow model evaluation. Part 1: Flow-over-terrain models. June 2014. URL: https://zenodo.org/record/4088287 (visited on 2020-10-14), doi:10.5281/zenodo.4088287.

YCTPA15

T.S. Yeow, A. Cuerva-Tejero, and J. Perez-Alvarez. Reproducing the Bolund experiment in wind tunnel. Wind Energy, 18(1):153–169, 2015. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/we.1688. URL: https://onlinelibrary.wiley.com/doi/abs/10.1002/we.1688 (visited on 2020-07-23), doi:10.1002/we.1688.