Implicit Large Eddy Simulation of Wind Flow Over Rough Terrain Bofu Wang, Fei Cheng, Hui Xu, Xuerui Mao

Implicit Large eddy simulation of wind flow over rough terrain Bofu Wang, Fei Cheng, Hui Xu, Xuerui Mao

To cite this version:

Bofu Wang, Fei Cheng, Hui Xu, Xuerui Mao. Implicit Large eddy simulation of wind flow over rough terrain. 17th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery (ISROMAC2017), Dec 2017, Maui, United States. hal-02981011

HAL Id: hal-02981011 https://hal.archives-ouvertes.fr/hal-02981011 Submitted on 27 Oct 2020

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés.

Distributed under a Creative Commons Attribution| 4.0 International License Implicit Large eddy simulation of wind flow over rough terrain

Bofu Wang1, Fei Chen2, Hui Xu3, Xuerui Mao1*

Abstract TATIN RO G N M Atmospheric boundary layer (ABL) ﬂow over a rough terrain is studied numerically. A drag force O A A C I H S model is employed to represent surface roughness. e spectral vanishing viscosity (SVV) method I O N

R is adopted for implicit large eddy simulation (ILES). ese techniques enable us to capture the M

Y Y

S dominant features of turbulent flow over non-smooth terrain with inhomogeneous roughness at high Reynolds number without resolving the very fine scale turbulence near the wall. We then run a test case of wind resource assessment on rough terrain, namely ABL flow over a hill with ISROMAC 2017 small-scale surface roughness at Re = 106 (Re is defined based on hill height). e instantaneous flow, mean flow and Reynolds stress are analysed in order to interprete the obtained turbulent flow. International Symposium on e effect of roughness length and associated drag coefficient are also evaluated. Transport Phenomena Keywords and Large eddy simulation — Rough terrain — Wind energy — Spectral element method — Spectral Dynamics of Rotating vanishing viscosity. Machinery 1Faculty of Engineering, University of Noingham, Noingham, NG7 2RD, UK Maui, Hawaii 2Xinjiang Goldwind Science & Technology Co., Ltd., Beijing, 100176, China 3Department of Aeronautics, Imperial College London, London SW7 2AZ, UK December 16-21, 2017 *Corresponding author: xuerui.mao@noingham.ac.uk

1. INTRODUCTION include detachment and recirculation under gentle slope as- sumption, and it reproduces quite satisfactorily the measured Modelling atmospheric turbulent flow over a complex terrain speed-up on Askervein at the hill top [3] but performs less with roughness is of great interest for many engineering, successfully on the lee side. Another prominent measure- environmental, and geophysics applications, such as wind ment is the Bolund field campaign [4, 5, 6] over Bolund Hill energy assessment; weather predictions; wind effects on agri- in Denmark which provided new experimental datasets of culture and forest; and transport and dispersion of pollutants the mean flow and the turbulence properties. e hill has a in the atmospheric environment. Specially, for wind energy elatively smaller scale (12 m high) than Askervein Hill, but it application, accurate wind resource assessment is critical offers a challenging topography with steep slopes and a cliff in the early stage development of wind farms. Inland wind (90°). e selection of Bolund Hill is justified by the three turbines are usually sited on complex rough terrains with dimensional nature of the hill and by the flow separations grass, forests, hills and mountains which make the evalua- induced by the cliff that make it more challenging to model tion of wind resources extremely challenging as the height of than gentler terrains. Up to now, these two measurements roughness is commensurate with the turbines. Moreover, the have been widely used as complex-terrain benchmark cases Reynolds number of turbulent Atmospheric boundary layer 8 for validating experimental and numerical implementations (ABL) is very high (up to 10 ), imposing extra challenges to as well as evaluating turbulence models. numerical assessment of wind resources and wind turbine performances. During the past decades, serious aempts Regarding surface roughness, in early times, some model have been made to predict the flow over complex terrain experiments were performed in wind-tunnel to obtain tur- with roughness through conducting field experiments, em- bulent flow fields within and above a model plant canopy. ploying physical models in wind-tunnel experiments, and Finnigan and Brunet [7] carried out a ”Furry hill” experiment carrying out numerical simulations. e most widely studied with a two-dimensional ridge using and provided a first view case is flow over forested hills. of the main features of canopy flow over a hill. Neff and e area with forested hills is commonly difficult to ac- Meroney [8] conducted single-hot film anemometer mea- cess, which prevents field measurements by either masts or surements to evaluate how hill and vegetation affect wind LiDAR. A well-known hill experiment was conducted by an power availability. More recently, Poggi and Katul [9, 10, 11] international group on the Askervein Hill in the UK [1, 2]. conducted a series laser Doppler anemometry measurements e slope of this hill is relatively smooth and the mean wind over a train of gentle and narrow forested hills in a flume profiles are extensively used for research with particular in- facility. eir results confirmed the theoretically predicted terest in the so-called linearised model. is model does not intermient recirculation region downstream of a hill cov- LES of flow over rough terrain — 2/6 ered by a dense canopy. ey also observed sweep motions 2. METHODS dominating momentum transfer within and near the canopy e governing equations for the current problem are the top all along the hill, while ejections dominating further incompressible NS equations with a volume force term rep- above. Teunissen et al. used two wind tunnels and constant- resenting the force induced by the roughness, temperaturehot-wire anemometry to study the Askervein 1 2 hill at three scales: 1/800, 1/1200, and 1/2500. Turbulence ∂t u+ u ∇ u = ∇p+Re− u+F with ∇ u = 0, (1) statistics gave more scaer compared to field observations ( · ) − ∇ · with more than 50% error [12]. Discrepancies observed were where u, p and F are velocity, pressure and volume force, higher in the wake and weaker at the location of the maxi- respectively. Re is the Reynolds number defined based on mum speed-up. Recently, Yeow et al. [13, 14, 15] performed the free-stream velocity and the height of the hill. Re = 106 a series of wind-tunnel analyses with hot-wire anemometry is used throughout this work if not otherwise specified. and particle image velocimetry to study flow around the Bol- und island with a 1/115 scale model and gave details on the 2.1 Spectral element method unsteady flow behaviour. ere are also some other exper- A majority of numerical tools for wind flow simulation are iments with complex topographies inducing complex flow based on finite volume methods. ere are also some ex- phenomena, such as Chock and Cochran [16]. amples of successful applications of element-based Galerkin methods in numerical weather prediction [27]. In this study, we adopt the open-source soware Nektar++ [28] which is Measurement from meteorological towers and model designed to support the development of high-performance experiments are currently indispensable in industrial prac- scalable solvers for partial differential equations using the tices, but not sufficient to analyse the intricate wind pro- spectral/hp element method. As a high order finite element file on complex terrains. Numerical models have the poten- method, it can deal with arbitrary geometric complexity, tial to offer insight with high fidelity information in this and are capable of local mesh adaption by either increasing regard. Wind flow in complex terrains exhibits complex the number of elements or increasing the polynomial order phenomenons such as flow stagnation, rapid acceleration de- within elements. Moreover, it has excellent scalability for celeration, recirculation-reaachment, vortex shedding, etc. parallel simulation and is scalable to more than one million Such micro-scale events occur at higher frequencies as well processor cores. as at smaller spatial scales compared to meteorologic scale events, which imposes challenges to numerical simulation. 2.2 Roughness model Many efforts have been devoted to the behavior of the flow In LES applied in offshore wind energy, the roughness (waves) phenomena over terrain hills using Reynolds averaged mod- is considered to be sub-scale of the computational grid, and els [17, 18], which predict the mean flow fields using far less therefore can be modelled using wall turbulence models, computational resoures. However, such models are based on which are commonly in the form of the boundary condition a variety of assumptions and intrinsically ignore the flow un- for shear stress. However regarding the deployment of wind steadiness; they are also known to be inaccurate in the flow turbines on rough terrains, the scale of the roughness ele- separation region. In recent years, due to the development ments, e.g. buildings and forest, is commensurate with the of numerical techniques and high-performance computing turbine and these elements cannot be considered as sub-grid. facilities, more advanced CFD methods, e.g. large-eddy sim- In the present work, a drag force approach is introduced to ulation (LES) and direct numerical simulation (DNS) [19, 20] model roughness. e drag term can be expressed as: are employed in wind energy. It has been demonstrated that the LES techniques are able to reproduce many observed fea- Fi = Cd Af √ujujui, (2) tures of turbulent flow over flat terrain with homogeneous roughness [21] and downwind of forest leading edges [22]. where i=1, 2, 3 refers to streamwise, lateral and vertical com- ere are investigations applying LES to flow over a rough ponents respectively, Cd is the mean drag coefficient of the hill (see [23] and references therein). Recently, developing roughness, Af is the density of roughness, and ui is the wind- LES models in the open source soware OpenFOAM [24] velocity component. Furthermore, the dimensionless rough- becomes very popular for wind energy application [25, 26]. ness length is defined as z0. Both Cd and Af are assumed to be ree dimensional calculations for the Bolund hill case have independent on wind velocity. In our following simulations, been performed for validations. Improved comparison with we set Af = 1, Cd = 0.4 and z0 = 0.1 if not otherwise speci- measurements [4, 5] were obtained. fied. Without loss of generality, the roughness is assumed to be with a uniform height.

is paper is oriented towards applying LES for ABL 2.3 Spectral vanishing viscosity method flows over complex terrains at very large Reynolds number. e spectral vanishing viscosity (SVV) method is adopted e boundary layer flow over an idealistic hilly terrain is stud- to account the small-scale turbulence. e SVV define a ied. is work paves the way for wind resource assessment viscosity in the spectral space that reaches a maximum at over an area with a mixed scale of roughnesses. high wave numbers, but vanishes for wave numbers below a LES of flow over rough terrain — 3/6 resolution-dependent threshold [29]. In this approach, one 50 adds an additional reaction term of the form to equations (1) Farfield ∂ ∂ ˆ u ǫSVV Q , (3) 1 ∂x ( ∗ ∂x ) High order outflow y where ǫSVV is a constant, * denotes the convolution operator, 0.5 and Qˆ is a kernel dictating which modes receive damping. U=1 0 is approach has been widely applied as a tool for implicit -1 0 1 U=(z/ )(1/7) large eddy simulation (ILES). e key point in SVV filtering Rough wall 0 -500 50 is that due to the shape of the kernel Qˆ k , x ( ) 2 N k Figure 2. Mesh in each x y plane. exp ( − ) , M < k N − Qˆ k =  (− M k 2 ) ≤ (4) ( ) ( − ) 0, k M   ≤ artificial viscosity for any mode number k is only applied 4. RESULTS AND DISCUSSION above a cutoff mode M. For these damped modes, the total 4.1 Instantaneous flow field viscosity can thus be expressed as 1 Re + ǫSVV . /

50 Farfield (a)

1 (b) High order Hill outflow y H 0.5 (c) W U=1 0 -1 0 1 (d) U=(z/ )(1/7) Rough wall 0 -50 0 50 x (e)

Figure 1. Schematic of the computational domain. Figure 3. Instantaneous contours of the streamwise velocity on horizontal cross sections (x-z planes) at (a) y=0.05, (b) y=0.1, (c) y=0.5, (d) y=1 and (e) y=2. 3. DISCRETIZATION e model of ABL flow over a small hill with uniform rough- We first evaluate the instantaneous flow fields. Figure ness near the wall is shown in Figure 1. e hill is quasi- 3 shows contours of streamwise velocity at various heights. three-dimensional (the spanwise direction is assumed to be Due to the roughness effect, upstream has already devel- homogeneous and discretised by Fourier transformation) and oped to turbulence. When this turbulent flow past the hill a only a two-dimensional (2d) plane is shown here. e shape turbulent wake region develops, the wind speed is reduced. of the hill is sinusoidal with height H = 1 and width W = 2. e upstream streamwise velocity is characterised by longi- At inflow, a power law velocity profile is applied. In order to tudinal paerns which can be regarded as the signature of trigger transition, a random perturbation with 5% turbulence coherent structures induced by the boom roughness. With intensity is applied at y < 2 in most of our calculations. A increasing height, these structures increase in size and look high order outflow condition is used at outlet [30]. At the like coherent streaks as usually observed in boundary layers. boom, no-slip wall condition is applied combined with a With further increasing height, fluctuations gradually de- uniform roughness with length z0. Far field condition is used crease and the laminar ambient flow is recovered. e effect at the top. e computational domain spans from -60 to 50, of the hill on instantaneous wind flow velocity is strong at 0 to 50 and 0 to π in the streamwise, vertical and spanwise the near wall region as shown in Figure 3 (a). When reaching directions, respectively. the hill, the longitudinal paern is almost eliminated in the Figure 2 shows the mesh generated by gmsh. ere are wake and then regenerated further downstream as shown in 2815 quadrilateral elements in each 2d plane, and at least one Figure 3 (b), (c), (d). e effect of the hill diminishes with in- element is set within the rough region near boom. In each creasing height which can be seen from the contrast between element, a spectral method is used to further decompose the upwind and downwind velocity. In Figure 3 (d) which shows element to a P+1 P+1 grid, where P represents a polyno- the contour at y=1 (hill height), stronger fluctuation can be mial order. P(= 6 is)×( adopted) in all our following simulations. observed manifesting higher turbulence intensity there. At In the homogeneous direction, a Fourier decomposition is y=2, the effect of the hill on instantaneous wind flow is no adopted and 32 Fourier modes are computed . longer observable (see Figure 3 e). LES of flow over rough terrain — 4/6 (a)

8 x=-25 x=-20 x=-15 x=-10 x=-5 x=3 y (z/ )(1/7) 4 (c)

0 0 0.4 0.8 1.2 U

Figure 4. Mean streamwise velocity proﬁles.

Figure 6. Averaged contour plots for (a) streamwise velocity and (b) spanwise velocity. 8 x=-25 x=-20 x=-15 x=-10 can be found in Figure 5. We can see that there is a small x=-5 zone with concentrated stream which forms low-level jet like x=3 y profile. In meteorological studies, low-level jets are caused 4 by a lot of mechanisms which is beyond the scope of current study. Figure 6 presents x-y slices of time- and z-averaged velocity components. From Figure 6 (a), the streamlines are distorted by the presence of the hill and a separated region on the lower lee side of the hill is generated as can be expected. 0 0 0.4 0.8 1.2 e flow reacts as if the hill and its separation region were a U single smooth obstacle. In Figure 6 (b), the average spanwise Figure 5. Mean streamwise velocity profiles with inflow velocity has a relative small magnitude and shows coherent perturbation applied at y < 4. structure similar in instantaneous field. is distribution suggests that the flow in spanwise direction is dominated by small scale structures. 4.2 Mean flow e average pressure filed is presented in Figure 7 (a). Figure 4 gives streamwise velocity profile sliced at different A large-scale horizontal pressure gradient appears due to streamwise positions x=-25, x=-20, x=-15, x=-10, x=-5 and topography, inducing the distortion of the mean flow. e x=3 as well as the 1/7 power law velocity profile. We can minimum pressure is located in wake region followed by an see the velocity profile does not have much variation with adverse pressure gradient. e surface pressure shown in upstream position (from x=-25 to x=-5), which is remarkably Figure 7 (b) also has a maximum before the hill. A minimum different from turbulent boundary layer flow. e extracted in surface pressure occurs as the flow passes the hill summit, profiles deviate from the power law profile in the near wall and large-magnitude oscillation is observed near hill top. e region where the velocity decreases a lile due to the surface variation of surface pressure is qualitative in agreement with roughness. e velocity profile indicates the boundary layer the wind tunnel measurement [32]. does not grow until near the hill side, a feature resulting from the uniform roughness. e profile at x=3 shows an inflection 4.3 Reynolds stress point near boundary which is caused by the inverse pressure Figure 8 shows the distribution of spanwise averaged stream- π T < u u >= 1 1 u u dzdt gradient in that region. From this profile, it is also observed wise normal Reynolds stress, ′ ′ π T ∫0 ∫0 ′ ′ , that the velocity above the hill is slightly larger than free where u′ is perturbation from the time averaged streamwise stream velocity due to the speed up. Note that one important velocity. It can be seen that the wake turbulence exhibit feature so called low-level jet [31] usually observed in ABL is plume-like structure [33]. e region of maximum produc- not captured in this case. We then increase the perturbation tion is some distance along the downstream of the hill. e region at inflow to y < 4. e streamwise velocity profiles structure extends downwind above the region of reversed LES of flow over rough terrain — 5/6 (a) 0.00012

z0=0.1, Cd=0.4 z =0.05, C =0.4 8E-05 0 d z0=0.1, Cd=0.2

4E-05 Wall shear stress shear Wall

0 0.2 -40 0 40 (b) x

Pressure Figure 9. Wall shear stress.

0 and drag coeﬃcient can strongly increase the wall shear stress.

-0.2-10 0 10 x REFERENCES Figure 7. (a) Average pressure distribution, (b) Surface [1] Taylor P. and Teunissen H. W. Askervein’82: report pressure. on the september/october 1982 experiment to study boundary-layer flow over askervein, south uist. Ontario technical report MSRB-83-8. Atmospheric Environment Service, Downsview, 1983. [2] Taylor P. and Teunissen H. W. e askervein hill project: report on the sept./oct. 1983, main field experiment. On- tario technical report MSRB-84-6. Atmospheric Environ- ment Service, Downsview, 1984. [3] Hunt J. C. R., Leibovich S., and Richards K. J. Turbu- lent shear flows over low hills. Q. J. Roy. Meteorol. Soc., 114:1435, 1988. < > Figure 8. Reynolds stress u′u′ . [4] Berg J., Mann J., Bechmann A., Courtney M. S., and Jørgensen H. E. e bolund experiment, part i: flow mean flow. is structure indicate high turbulence level in over a steep, three-dimensional hill. Boundary-Layer this region and it is associated with strong shear layer for- Meteorol., 141:219–243, 2011. mulated on the top of the recirculation zone. [5] Bechmann A., Sørensen N. N., Berg J., Mann J., and Figure 9 shows wall shear stress at boom surface with Rethor¨ e¨ P. E. e bolund experiment, part ii: blind varies roughness length z0 and drag coefficient Cd. Take comparison of microscale flow model. Boundary-Layer z0 = 0.1, Cd = 0.4 as a reference case, either decreasing z0 or Meteorol., 141:245–271, 2011. C decreasing d will lead to an increase of wall shear stress. A [6] Bechmann A., Berg J., Courtney M., Jørgensen H., Mann high peak in wall shear stress is observed at hill top. J., and Sørensen N. N. e bolund experiment: Overview and background. Danmarks Tekniske Universitet Risø 5. SUMMARY National Laboratory for Sustainable Energy, 2009. In this paper, we have carried out ILES to investigate the [7] Finnigan J. J. and Brunet Y. Turbulent airflow in forests 6 ABL flows over a small hill at Re = 10 . It is found that on flat and hilly terrain. in Wind and trees, Cous MP, under the effect of surface roughness, the flow upstream of Grace J (eds). Cambridge University Press: Cambridge, the hill behaves like a fully developed turbulent boundary UK., pages 3–40, 1995. layer. Streak-like structures can clearly be observed from [8] instantaneous flow fields. Low-level-jet like velocity profile Neff D. E. and Meroney R. N. Wind-tunnel modeling of can be captured due to roughness effect. As the turbulent hill and vegetation influence on wind-power availability. flow past the hill, a wake region is developed on the lee side of J. Wind Eng. Ind. Aerodyn., 74-76:335–343, 1998. the hill. e wake region is characterised by a recirculation [9] Poggi D. and Katul G. G. An experimental investigation zone, a strong elevated shear layer downstream from the of the mean momentum budget inside dense canopies summit and high turbulence levels. e wake turbulence is on narrow gentle hilly terrain. Agric. Forest Meteorol., characterised by plume-like structure. e roughness length 144:1–13, 2007. LES of flow over rough terrain — 6/6

[10] Poggi D. and Katul G. G. Turbulent flows on forested [24] Opencfd ltd. user and programmer’s guide. hilly terrain: e recirculation region. Q. J. R. Meteorol. hp://www.openfoam, 2012. Soc., 133:1027–1039, 2007. [25] Churchfield M. J., Lee S., Michalakes J., and Moriarty P. [11] Poggi D., Porporato A., Ridolfi L., Albertson J. D., J. A numerical study of the effects of atmospheric and and Katul G. G. e effect of vegetation density on wake turbulence on wind turbine dynamics. J. Turbul., canopy sub-layer turbulence. Boundary-Layer Meteorol., 13:N14, 2012. 111:565–578, 2004. [26] Vuorinen V., Chaudhari A., and Keskinen J. P. Large- [12] Teunissen H. W., Shokr M., Bowen A., Wood C., and eddy simulation in a complex hill terrain enabled by a Green D. e askervein hill project: wind-tunnel simu- compact fractional step openFOAM solver. Adv. Eng. lations at three length scales. Boundary-Layer Meteorol., So., 79:70–80, 2015. 40:1–29, 1987. [27] Simone Marras, James F. Kelly, Margarida Moragues, An- [13] Yeow T. S., Cuerva A., Conan B, and Perez´ Alvarez´ J. dreas Muller,¨ Michal A. Kopera, Mariano Vazquez,´ Fran- Wind tunnel analysis of the detachment bubble on bol- cis X. Giraldo, Guillaume Houzeaux, and Oriol Jorba. A und island. In: Proceedings of the science of making torque review of element-based galerkin methods for numerical from wind, Oldenburg, Germany, 2012. weather prediction: Finite elements, spectral elements, [14] Yeow T. S., Cuerva A., Conan B, and Perez´ Alvarez´ and discontinuous galerkin. Arch. Comput. Method E., J. Wind tunnel analysis of the detachment bubble on 23:673–722, 2016. bolund island. Journal of Physics: Conference Series, [28] Cantwell C.D., Moxey D., Comerford A., Bolis A., Rocco 555:012021, 2014. G., Mengaldo G., Grazia D. De, Yakovlev S., Lombard [15] Yeow T. S., Cuerva A., and Perez´ Alvarez´ J. Reproducing J.-E., Ekelschot D., Jordi B., Xu H., Mohamied Y., Eskils- the bolund experiment in wind tunnel. Wind Energy, son C., Nelson B., Vos P., Bioo C., Kirby R. M., and 18:153–169, Wind Energy. Sherwin S. J. Nektar++: An open-source spectral/hp element framework. Comput. Phys. Comm., 192:205–219, [16] Chock G. and Cochran L. Modeling of topographic 2015. wind speed effects in hawaii. J. Wind Eng. Ind. Aerodyn., 93:623–638, 2005. [29] Tadmor E. Convergence of spectral methods for nonlin- ear conservation laws. SIAM J. Numer. Anal., 26:30–44, [17] Finardi S., Brusasca G., Morselli M., Trombei F., and 1989. Tampieri F. Boundarylayer flow over analytical two- dimensional hills: A systematic comparison of different [30] Dong S., Karniadakis G. E., and Chryssostomidis C. A models with wind tunnel data. Boundary Layer Meteorol., robust and accurate outflow boundary condition for in- 63:259–291, 1993. compressible flow simulations on severely-truncated unbounded domains. J. Comput. Phys., 261:95–136, 2014. [18] Ying R. and Canuto V. Numerical simulation of flow over two-dimensional hills using a second-order turbulence [31] Wagner A. Uber die tageswinde in der freien atmosphire. closure model. Boundary Layer Meteorol., 85:447–474, Grrl. Britr. Grophys., 1939. 1997. [32] Brunet Y., Finnigan J. J., and Raupach M. R. A wind [19] E. Kleusberg, R. F. Mikkelsen, P. Schlaer, S. Ivanell, and tunnel study of air flow in waving wheat: Single-point D. S. Henningson. High-order numerical simulations of velocity statistics. Boundary Layer Meteorol., 70:95–132, wind turbine wakes. J. Phys. Conf. Ser., 854:012025, 2017. 1994. [20] A. Ducoin, M. S. Shadloo, and S. Roy. Direct numerical [33] Keith W. Ayoe, Peter P. Sullivan, and Edward G. Paon. simulation of flow instabilities over savonius style wind Les and wind tunnel modelling over hills varying steep- turbine blades. Renew. Energ., 105:374–385, 2017. ness and roughess. e Fih International Symposium [21] Su H. B., Shaw R. H., Paw U. K. T., Moeng C. H., and Sul- on Computational Wind Engineering, 2010. livan P. P. Turbulent statistics of neutrally stratified flow within and above a sparse forest from large-eddy simulation and field observations. Boundary-Layer Meteorol., 88:363–397, 1998. [22] Yang B., Raupach M., Shaw R. H., Paw U. K. T., and Morse A. P. Large eddy simulation of turbulent flow across a forest edge. part i: Flow statistics. Boundary- Layer Meteorol., 120:377–412, 2006. [23] Dupont S., Bruneta Y., and Finnigan J. J. Large-eddy simulation of turbulent flow over a forested hill: Validation and coherent structure identification. Q. J. R. Meteorol. Soc., 134:1911–1929, 2008.