aeolus header

Creating a Wind Field

The toolbox provides two methods for generating a wind field. One is interactive and the other is programatically from a Matlab script.

Both methods requires the set of data listed below

  • Average wind speed

  • Turbulence intensity

  • Grid spacing

  • Length of wind field

  • Width of wind field

  • Maximum simulation time

  • (for the new wind genetation method turbine positions are also needed)

The interactive method will ask for this data and generate the appropriate data structures etc. you only need to double-click the Generate Wind Field - Taylor or Generate Wind Field - No Taylor block in the Wind simulink library to start the procedure. Alternatively, you can also run the scripts create_windfield or create_windfield(farm), where farm is a farm structure.

The programmatic method requires that you place this data in appropriate data structures and invoke the gen_windfield or gen_windfield_no_taylor method.

This is illustrated with two examples for both the taylor and no_taylor versions. The output mat files are used when creating wind farm models (see creating wind farm models)

Wind field creation using Taylers frozen turbulence hypothesis method

%The average wind speed is 12 m/s
U0=12;

%The turbulence intensity is 0.1.
Ti=.1;

%The grid spacing is 15m
d=15;

%The length of the wind field is 1Km
Lx=1e3;

%The width of the wind field is 1.5Km
Ly=1.5e3;

%The maximum simulation time is 1000s
SimTime=1000;

%Generate the wind field
wind=gen_windfield(U0,Ti,d,Lx,Ly,SimTime);

%Save the wind field
save('benchmarkwind','wind');

Wind field creation using Taylers frozen turbulence hypothesis method

%The farm has three turbines, all of them the NREL5MW
farm.turbines={'NREL5MW', 'NREL5MW', 'NREL5MW'}

%The position of the turbines.
%(make sure that the turbines are within the grid of the windfield)
farm.pos=[0    100    200; %X values
	 200  300    400] %Y values

%The average wind speed is 12 m/s
U0=12;

%The turbulence intensity is 0.1.
Ti=.1;

%The grid spacing is 15m
d=15;

%The length of the wind field is 1Km
Lx=1e3;

%The width of the wind field is 1.5Km
Ly=1.5e3;

%The maximum simulation time is 1000s
SimTime=1000;

%Generate the wind field
wind=gen_windfield_no_taylor(U0,Ti,d,Lx,Ly,farm, SimTime);

%Save the wind field
save('benchmarkwind_no_taylor','wind')