Terrain Effects on Wind Flow: Simulations With an Immersed Boundary Method

The modelling of the wind resource over arbitrary topography is required to optimize the micrositing of wind turbines. Most solvers use classical body-fitted grid for simulations. In such an approach, to cover the wind rose using a rectangular domain, a dedicated mesh must be generated for each direction. Moreover, over complex terrain, additional numerical errors are introduced in the solver due to coordinate transformations. To overcome these challenges and to facilitate the grid generation process, an immersed boundary method is developed in connection with a RANS solver in order to simulate turbulent atmospheric flows over arbitrary topography. In this method, a Cartesian grid is used and the boundary condition on the terrain surface is modelled within the solver using a “direct forcing” approach. With the immersed boundary method a rectangular grid can be used to simulate the flow field for all wind directions and only a rotation of the digital elevation map is required. Ghost cells are used to enforce the desired boundary condition at the immersed surface. The immersed boundary method developed in this work is used to simulate the flow in connection with both Baldwin-Lomax and kω turbulence models. The performance of the method is examined for the flow over a two-dimensional hill. Results are compared with experimental data as well as a classical body-fitted grid to isolate the effect of the boundary conditions. The comparisons show good agreement among all the results. The results for the three-dimensional wind flow simulation over the Askervein Hill test case are also presented, and show the capability of the immersed boundary method in a full-scale scenario.Copyright © 2011 by ASME

[1]  N. Mortensen,et al.  Exploring the limits of WAsP the wind atlas analysis and application program , 1996 .

[2]  G. Iaccarino,et al.  RANS solvers with adaptive structured boundary non-conforming grids , 2002 .

[3]  Yannick Meillier,et al.  Measurements of Boundary Layer Profiles with In Situ Sensors and Doppler Lidar , 2008 .

[4]  A. S. Grove,et al.  An experimental investigation of the steady separated flow past a circular cylinder , 1964, Journal of Fluid Mechanics.

[5]  E. F. Bradley,et al.  Boundary-layer flow over low hills , 1987 .

[6]  H. Lomax,et al.  Thin-layer approximation and algebraic model for separated turbulent flows , 1978 .

[7]  P. Taylor,et al.  The Kettles Hill Project: Field observations, wind-tunnel simulations and numerical model predictions for flow over a low hill , 1988 .

[8]  H. I. Andersson,et al.  Non-Linear, Microscale Modelling of the Flow Over Askervein Hill , 2006 .

[9]  Reza S. Abhari,et al.  Seven-Sensor Fast-Response Probe for Full-Scale Wind Turbine Flowfield Measurements , 2011 .

[10]  J. Palma,et al.  Simulation of the Askervein flow. Part 2: Large-eddy simulations , 2007 .

[11]  V. C. Patel,et al.  Numerical simulation of wind flow over hilly terrain , 2000 .

[12]  Jacob Berg,et al.  The Bolund experiment: Overview and background , 2009 .

[13]  Peter A. Taylor,et al.  Boundary-layer flow over topography : Impacts of the Askervein study , 1996 .

[14]  José M. L. M. Palma,et al.  Numerical simulation of the atmospheric flow in a mountainous region of the North of Portugal , 1998 .

[15]  Pollutant dispersion over two-dimensional hilly terrain , 1998 .

[16]  F. Porté-Agel,et al.  A scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric boundary layer , 2000, Journal of Fluid Mechanics.

[17]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[18]  Reza S. Abhari,et al.  Three-Dimensional Flow Prediction and Improvement of Holes Arrangement of a Film-Cooled Turbine Blade Using a Feature-Based Jet Model , 2007 .

[19]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[20]  Almerindo D. Ferreira,et al.  Experimental and numerical simulation of flow around two-dimensional hills , 1995 .

[21]  R. Bouard,et al.  Experimental determination of the main features of the viscous flow in the wake of a circular cylinder in uniform translation. Part 1. Steady flow , 1977, Journal of Fluid Mechanics.

[22]  A. R. Brown,et al.  Large-Eddy Simulation Of Turbulent Separated Flow Over Rough Hills , 2002 .

[23]  J. Palma,et al.  Simulation of the Askervein Flow. Part 1: Reynolds Averaged Navier–Stokes Equations (k∈ Turbulence Model) , 2003 .

[24]  N. Ron-Ho,et al.  A Multiple-Grid Scheme for Solving the Euler Equations , 1982 .

[25]  S. Biringen,et al.  Numerical Simulation of a Cylinder in Uniform Flow , 1996 .

[26]  J. Lundquist,et al.  An Immersed Boundary Method for the Weather Research and Forecasting Model , 2014 .

[27]  P. A. Taylor,et al.  Some numerical studies of surface boundary-layer flow above gentle topography , 1977 .

[28]  G. D. Raithby,et al.  The Askervein hill project: A finite control volume prediction of three-dimensional flows over the hill , 1987 .