3D Turbulent Wind Field Simulation by means of Sandia Method for wind energy applications

3D turbulent wind field generator with Kaimal model

Description

This code is meant for all those people interested or actively involved within the wind energy business.

It tries to couple the best features of already existing tools, such as TurbSim and Bladed.

The code has been quite tested and reviewed since the previous version, in order to overcome issues dealing with the Coherence function as well as the Gaussian statistics.

Easy to use and free-source.

Inputs:

1) Hub Height [m]
2) Grid Width [m]
3) Grid Height[m]
4) Number of points along Y
5) Number of points along Z
6) Number of FFT points (2^12 recommended)
7) Analysis Time (generally 600s)
8) Hub Average Wind Speed [m/s]
9) Coherence Model (TurbSim or Bladed)

Outputs:

1) U: Turbulent time series for v-component [m/s]
2) V: Turbulent time series for v-component [m/s]
3) W: Turbulent time series for v-component [m/s]

• Required Products: Statistics Toolbox
• MATLAB release: MATLAB 8.0 (R2012b)

[Jean]

Thank you for this script.
I think there is a mistake in the IEC 61400-1 implementation, the power spectrum for Kaimal should be with a 5/3 power and not a 5/6:

PSD(i,IVec)=SigmaLU(IVec)/(1.0 + LU(IVec)*Freq(i))^Exp3

with Exp3=5/3
Moreover the calculated time series is artifically scaled to the desired mean & standard deviation, but the calculation output seems to be very far from the desired value. I therefore suspect an error in the code.
A lot of lines could be vectorized.

[Hamid Sarlak]

 

[Francesco Perrone]

You should download the lapack package here

http://www.mathworks.com/matlabcentral/fileexchange/16777

as said within "Other Requirements!

[Allan]

Thank you for your sharing work.
I met the following problem. It seems that function lapack can't be found.
Is this subroutine needs additional toolbox?

Undefined function 'lapack' for input arguments of type 'char'.

Error in Sandia_3D_method (line 287)
U = lapack('SPPTRF','L',Ntot,H,0);