27.02.2013 Views

Wind Energy

Wind Energy

Wind Energy

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

7<br />

Models for Computer Simulation of <strong>Wind</strong> Flow<br />

over Sparsely Forested Regions<br />

J.C. Lopes da Costa, F.A. Castro and J.M. L.M. Palma<br />

7.1 Introduction<br />

Because of the flow complexity and wind power reduction, the regions near or<br />

within forests have not been considered for wind park installation. However,<br />

because of the shortage of good clear sites, this situation was reverted, adding<br />

new motivation to the study of the wind flow over forested regions.<br />

The so-called k − ε model, the standard in flow modelling of engineering<br />

applications, has also appeared in atmospheric flow studies, circumventing the<br />

limitations of simpler approaches [1]. However, in the case of canopy flows,<br />

there are open questions on the actual performance of the k − ε compared<br />

with lower order models, and the mathematical form of the additional terms<br />

in the transport equations for the turbulence kinetic energy k and its rate of<br />

dissipation ε [2].<br />

7.2 Mathematical Models<br />

The canopy models as in [4, 6] and a version of k − ε − v2 − f [3] are studied<br />

over model forests. For sake of compactness, the model description will be<br />

restricted to their differences from the standard versions, i.e. without the terms<br />

accounting for the canopy. The reader may refer to the original publications<br />

for a complete description of each of the models being used here.<br />

The presence of the canopy requires an additional term in the momentum<br />

equation accounting for the extra drag, and re-definition of the production<br />

terms of k and ε in the transport equations, where additional terms Sk and<br />

Sε will appear (Pk − ρε + Sk) and (C1ε/kPk − C2ρε2 /k + Sε). In case of<br />

k −ε−v 2 −f turbulence model, the equation for f contains also an additional<br />

Sk, Cf2(Pk + Sk)/k.<br />

The terms Sk and Sε will read<br />

Sk = 1<br />

2 ραCD<br />

�<br />

βp<br />

�<br />

�U � �3 �<br />

− βd �U � �<br />

� k , (7.1)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!