WIN1 Part 01-1
WIN1 Part 01-1
2023
Contents
1
12.10.2023
Literature
2
12.10.2023
Principle
• Sailing against the wind using an airfoil
• Even faster than the wind
Flettner Rotor
1920ies
Rotating rotors
Magnus-Effekt
Auxiliary drive
Fuel saving concept
3
12.10.2023
Sky-Sails
Tractive power support for cargo ships by
automated kite
Exploitation of wind at heights
> 100 m
Kite is steered in an 8-shape
4
12.10.2023
European Bock wind mill (12th century) Tower wind mill – widespread in the Mediterranean
5
12.10.2023
Source: https://en.wikipedia.org/wiki/Windcatcher
Halladay windmill – invented around 1850 by D. Halladay Dutch windmill (16th century)
6
12.10.2023
7
12.10.2023
8
12.10.2023
9
12.10.2023
Enercon E-126
‘20
10
12.10.2023
‘20
11
12.10.2023
12
12.10.2023
13
12.10.2023
14
12.10.2023
Orography
Roughness of the ground
• Temperature
• Roughness
• Obstacles
• Orography
Up wind
Land wind
Sea wind
Sea Beach Grass Trees Hills Forests Houses Cities Valleys Mountains
15
12.10.2023
Characteristic Properties
Wind Statistics
1
𝑣 = 𝑣
𝑁
Standard deviation
1
𝜎= 𝑣 −𝑣 )²
𝑁−1
16
12.10.2023
12 5000 40
Beitrag zur
contribution Leistungsdichte
to power = xHäufigkeit
density = frequency power density * Leistungsdichte
Häufigkeit
frequency Leistungsdichte
power density 4500 35
10
4000
power densityin(W/m²)
power densityin(W/m²)
W/m²
W/m²
30
3500
8
in %
Häufigkeit(%)
3000 25
Leistungsdichte
Leistungsdichte
frequency
6 2500 20
2000 15
4 1500
10
1000
2
500 5
0 0 0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
>19
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
>19
Windgeschwindigkeitsklasse
wind speed class Windgeschwindigkeitsklasse
wind speed class
17
12.10.2023
Weibull distribution
𝑘 𝑣
𝑓 𝑣 = 𝑒
𝐴 𝐴
𝑣 = 𝐴 (−𝑙𝑛(1 − 𝐹 𝑣 )
18
12.10.2023
The average wind speed 𝑣̅ and the variance 𝜎 can be calculated from the
Weibull distribution by:
1
𝑣̅ = A Γ 1 +
𝑘
𝜎 =𝐴 Γ 1+ -𝐴 Γ 1+
Wind Rose
19
12.10.2023
Wind Rose
https://www.weather.gov/ama/dalhartwindroseinformation
Geostrophic
wind
Surface
20
12.10.2023
200
ℎ D v = 3.2m/s
𝑣 𝑙𝑛
𝑧
= 150
𝑣 ℎ
𝑙𝑛
𝑧
Height (m)
100
50
h1 height 1 [m] z0 = 0,40 z0 = 0,20 z0 = 0,05
h2 height 2 [m]
v1 wind speed in h1 [m/s] 0 D v = 0.6m/s
v2 wind speed in h2 [m/s]
0 2 4 6 8 10
z0 roughness length [m]
Average wind speed (m/s)
21
12.10.2023
Height
free atmosphere
1000 m
Ekman layer
turbulent boundary layer
4,500
kW
2,500
100 m kW
1,500
kW
h1 Höhe 1 [m]
600
h2 Höhe 2 [m] 500 kW
Prandtl
v1 layer
Windgeschwindigkeit 300
in kW
h1 [m/s]
kW
v2 Windgeschwindigkeit
50 in h2 [m/s]
z0 5 mm kW
Rauhigkeitslänge [m]
laminar boundary layer
Turbulence Intensity Iv
𝜎
𝐼 =
𝑣̅
Gusts:
Deviation of the average wind speed in
a range of about 10 sec
22
12.10.2023
cup anemometer
propeller anemometer
23
12.10.2023
ultrasound anemometer
2D 3D
24
12.10.2023
Wind Data
European Windatlas
Austrian Wind Potential Analysis
AuWiPot
www.windatlas.at
Weather stations
25
12.10.2023
Assumption:
𝑣(𝐴) = const.
𝜌 = const.
𝑣
• Kinetic Energy:
1
𝐸 = m𝑣
2
• Power:
d𝐸
P=
𝑑𝑡
• v = const.
1 d(m𝑣 ) 1 𝑑𝑚 1
P= = v² = 𝑚̇v²
2 𝑑𝑡 2 𝑑𝑡 2
• Mass flow 𝑚̇
𝑚̇ = ϱ 𝑉̇ = ϱ A v
26
12.10.2023
𝑃 = 𝜚𝐴𝑣³ = 𝜚𝑣³
v [m/s] E [kWh/m²·a]
4 343
5 671
6 1.159
7 1.840
8 2.747
27
12.10.2023
Air Density
Fluid Mechanics
Bernoulli law:
1
𝑝 + 𝜌𝑣² = const
2
28
12.10.2023
Fluid Mechanics
Incompressible Fluid
𝐴 ,𝑣 𝐴 ,𝑣
1 1 1
𝑃 = 𝑃 − 𝑃 = 𝜚𝐴 𝑣 − 𝜚𝐴 𝑣 = 𝜚 𝐴 𝑣 − 𝐴 𝑣
2 2 2
29
12.10.2023
𝐴
𝜌𝐴 𝑣 = 𝜌𝐴 𝑣 → 𝑣 = 𝑣
𝐴
1 1
→𝑃 = 𝜚𝑣 𝐴 𝑣 − 𝑣 = 𝑚̇ 𝑣 − 𝑣
2 2
→𝑃 → maximum
𝐹 = 𝑚̇ ∆𝑣 = 𝑚̇ 𝑣 − 𝑣
𝑃 = 𝐹𝑣 = 𝑚̇ 𝑣 − 𝑣 𝑣
1
𝑚̇ 𝑣 − 𝑣 𝑣 = 𝑚̇ 𝑣 − 𝑣
2
1
𝑣 = 𝑣 +𝑣 (Froude-Rankine Theorem)
2
1
→ 𝑚̇ = 𝜌𝐴 𝑣 + 𝑣
2
1
→𝑃 = 𝜌𝐴 𝑣 + 𝑣 𝑣 −𝑣
4
For the derivation of the Froude-Rankine Theorem see e.g. Gasch/Twele.
30
12.10.2023
1
𝑃 𝜌𝐴 𝑣 + 𝑣 𝑣 −𝑣 1 𝑣 +𝑣 𝑣 −𝑣
𝑐 = =4 =
𝑃 1 2 𝑣
𝜚𝐴𝑣
2
1 𝑣 𝑣
= 1+ 1−
2 𝑣 𝑣
Ratio →𝜉
1
𝑐 = 1−𝜉 1+𝜉
2
Betz coefficient
16
𝑐 , = = 0.593
27
Specific Power
1000 1
Ideal
900 0,9
Betz
800 0,8
Pitch
Coefficient of power
700 0,7
densityin(W/m²)
Stall
W/m²
500 0,5
400 0,4
Power
300 0,3
200 0,2
100 0,1
0 0
0 5 10 15 20 25
Windgeschwindigkeit in m/s
Wind speed (m/s)
31
12.10.2023
Coefficient of power cP
16
𝑐 , = = 0.59
27
ratio
32
12.10.2023
Two Principles:
Aerodynamic Drag
Drag Force 𝐹 : A
cw = 1.11
vw
1
𝐹 =𝑐 𝜌𝐴𝑣
2
1 cw = 1.33
vw
∆𝑝 = 𝜌𝑣
2
33
12.10.2023
1
𝐹 =𝑐 𝜌𝐴 𝑣 − 𝑢 ²
2
1
𝑃=𝐹 𝑢=𝑐 𝜌𝐴 𝑣 − 𝑢 ² 𝑢
2
𝑐 coefficient of drag
𝑣 wind speed
𝑢 speed of the object
𝐹 drag force
𝐴 area approached by incident flow
34
12.10.2023
𝑢
𝜆= <1
𝑣
𝑢 circumferential speed
𝑣 wind speed
Coefficient of power 𝑐
4
𝑐 , = 𝑐
27
*dt. Schnelllaufzahl
35
12.10.2023
Aerodynamic Lift
𝒗𝒘
36
12.10.2023
Lift force: 𝐹 𝐹
1
𝐹 = 𝜌𝑣 ²𝑐 𝐴
2
Drag force:
1
𝐹 = 𝜌𝑣 ²𝑐 𝐴
2
𝐹
Resulting force: 𝑣
𝐹 = 𝐹 ²+𝐹 ²
37
12.10.2023
Goal of a lift driven rotor: max. lift force at lowest possible drag force!
Lift/drag ratio* 𝜀:
𝐹 𝑐 (𝛼 )
𝜀(𝛼 ) = =
𝐹 𝑐 (𝛼 )
𝑐 lift coefficient
𝑐 drag coefficient
𝛼 angle of attack
*dt. Gleitzahl
𝜶
𝐹
𝑣
38
12.10.2023
𝜶
𝑣
𝐹
Lift and drag force get smaller with increasing angle of attack.
Angle of attack 𝛼 = 0
𝑣
𝐹
39
12.10.2023
𝐹 𝐹
𝑣 𝐹
At too high angle of attack 𝛼 the flow stalls at the trailing edge.
𝐹
𝑣
𝐹 𝐹
40
12.10.2023
Angle of attack 𝛼
𝑐 = 𝑐 (𝛼)
𝑐 = 𝑐 (𝛼)
Angle of attack 𝜶
𝑢
𝜆= >1
𝑣
𝑢 circumferential speed
𝑣 wind speed
* dt. Langsamläufer
** dt. Schnellläufer
41
12.10.2023
Depends on:
wind speed
circumferential speed
Rotor Blade
42
12.10.2023
Theoretical optimum
near optimum
trapez
rectangular
Rotor Blade
43
12.10.2023
Rotor Blade
ideal
Rotor Blade
44
12.10.2023
Blade tip-losses
Another source of losses is the flow around the tip of the blade from the
pressure side (lower side of the profile) to the suction side (upper side).
This causes the blade lift to decreases toward the tip.
45
12.10.2023
Stall Control
Rotor blade is fixed
Pitch angle fixed
Angle of attack depends on the rotor revolution speed and wind speed
Power < 1 MW
46
12.10.2023
Active Stall
Ratio of the Forces for a Lift Device E-66, gearless drive train (ENERCON)
47
12.10.2023
Pitch Control
• Rotor blade is rotatable around its axis (profile nose towards incoming flow)
Angle of attack depends on the rotor revolution speed and wind speed
• Power > 500 kW
48
12.10.2023
49
12.10.2023
Power Curves
Power
Power
P rated P rated
Rotor
The rotor area determines the energy yield of
the wind turbine.
The largest wind turbines (e.g. Enercon E-
127) has a rotor diameter of 127 m and a
swept area of 12.700 m²
Rotor speed:
5 – 20 rpm (> 1,5 MW)
15 – 50 rpm (< 1,5 MW)
~ 40 m 500 kW
50 – 60 m 1 MW
70 – 90 m 2 MW
90 – 115 m 3 MW
115 – 130 m 5 MW
50
12.10.2023
51