Abraham_9781606504840_Chapter 06
Abraham_9781606504840_Chapter 06
Numerical Simulations of
Small Wind Turbines—
HAWT Style
Jimmy C. K. Tong
Here, the author introduces the important topic of product certification and
code compliance. Next, he focuses on cases that are considered for numerical
investigation of horizontal-axis wind turbine (HAWT) performance. Among
the cases are different design environments with varying wind gusts, wind
direction, and turbulence. Creation of the numerical mesh and application of
the flow boundary conditions are discussed in detail. So too are the various
turbulence models that modern computational simulators employ.
An emerging field in this area is the coupled interaction of the fluid
and the solid object. This analysis approach, termed Fluid Structural
Interaction, is finding use among research scientists as computer capa-
bilities increase. Finally, the author discusses the types of information
that are obtained from numerical simulation and how that information
guides the design process.
6.1 INTRODUCTION
In order to introduce a new wind turbine model into a market, most coun-
tries require or encourage product certification that is in compliance with
local or international standards. To qualify for the certification, the wind
turbine will need to pass the safety and standard performance testing and
design verifications specified by the designated authority. This rigorous
process is employed to ensure that the wind turbine will meet a certain
quality standard and is safe to use, and the advertised performance is com-
parable to other products available in the market.
Design verification is part of the certification process, and numeri-
cal simulation can be helpful to demonstrate that the structural and safety
requirements are met under the operating conditions. In general, numerical
simulation uses the code-specified incoming wind condition to determine
the wind loading acting on the blades and the turbine structure. Hence, the
wind turbine can be designed to meet required material safety factors. This
will be discussed in more detail later in this chapter.
Beyond design verification, the numerical simulation is also bene-
ficial for designers to improve the performance of the wind turbine. In
particular, a few computer codes are developed for blade design in order
to capture kinetic energy from moving air. As the entire wind turbine com-
prises more components, the general CFD software is applicable for the
complete system analyses. Details about the approach of numerical simu-
lation, the analysis process with inputs and outputs, and the design itera-
tion are discussed later in this chapter.
IEC 61400-2: 2006 Ed 2.0 Design requirements for small wind tur-
bines [1].
IEC 61400-11: 2012 Ed 3.0 Acoustic noise measurement techniques [2].
IEC 61400-12-1: 2005 Ed 1.0 Power performance measurements of
electricity-producing wind turbines [3].
IEC 61400-14: 2005 Ed 1.0 Declaration of apparent sound level and
tonality values [4].
IEC 61400-21: 2008 Ed 2.0 Measurement and assessment of power
quality characteristics of grid-connected wind turbines [5].
IEC 61400-22 2010 Ed 1.0 Conformity testing and certification [6].
IEC 61400-23: 2001 Ed 1.0 Full scale structural testing of rotor
blades [7].
While there are local code variations, most of the codes are very similar
to the IEC 61400-2 [1]. Taking IEC 61400-2 as the basis for wind turbine
design is commonly accepted by developers who intend their products for
the global market. The standard specifies the design methodology, which
includes simplified load equations, aeroelastic modeling, and mechanical
134 • SMALL-SCALE WIND POWER
Table 6.1. Set of Design Load Cases for aeroelastic models from IEC
61400-2.
Design DLC Wind condition Other Type of
situation conditions analysis
1) Power 1.1 NTM Vin < Vhub < F, U
production Vout or 3Vave
1.2 ECD Vhub < Vdesign U
(Continued)
NUMERICAL SIMULATIONS OF SMALL WIND TURBINES • 135
With the code and standard, which cover the certification process and
define the safety and performance requirements of a wind turbine,
numerical simulation is a common tool for the developers to properly
assess the design and to effectively make design improvement through
the iterative product development cycle. Out of the various important
components of a wind turbine, like the blades, generator, brakes, con-
trols of torque, pitch and yaw, and the tower, the numerical simula-
tion discussed in this chapter only covers the mechanical aspect of the
design. Other aspects, such as electrical and electromagnetic aspects, are
not discussed.
As the blades and the rest of the turbine system (rotor, nacelle, and
tower) are components that directly interact with the wind, numerical sim-
ulation for mechanical/structural loading and power generation (through
transforming the kinetic energy from moving air into rotating mechanical
energy for the generator to turn into electricity) can be determined through
aeroelastic modeling and CFD modeling. For aeroelastic modeling, the
Blade Element Momentum (BEM) theory has been commonly employed
for blade design and airfoil selection. BEM allows simple calculations of
steady loads, thrust, and power for different settings of wind speed, rota-
tional speed, and pitch angle. Beyond steady calculations, and with the
advancement of computation power, CFD is becoming the standard tool
applied to wind turbines subjected to complex flow situations. Designers
and engineers are still using both tools, independently or in combination,
NUMERICAL SIMULATIONS OF SMALL WIND TURBINES • 137
for the turbine design. Figure 6.1 shows streamline patterns from a typical
CFD simulation for a conventional three-bladed HAWT.
Figure 6.3. Meshing for wind turbine blade and rotor CFD simulation.
Source: Reproduced from Bazilevs et al. [11].
140 • SMALL-SCALE WIND POWER
blade and rotor surfaces to resolve the boundary layer of the fluid on the
leading and trailing edges and tip of the blade. With only one-third of the
domain modeled, it is possible to obtain a solution using a numerical tech-
nique of periodicity. The periodic boundary conditions can be used when
the flows going out through one boundary reappear as the flows going in
through the opposite boundary. If the flow pattern is too complex for the
periodicity to be appropriate, the same mesh can be repeated two more
times every 60°. When the complete face of the blade and rotor is modeled,
then the usual CFD analysis without using periodicity can be conducted.
The other type of analysis covers the entire wind turbine system,
which includes all the blades, rotor, nacelle, and tower structure. This
analysis can be useful to simulate a real site condition or a wind tunnel
test situation. A sample mesh for this simulation is shown in Figure 6.4.
As seen in Figure 6.4, the solution domain is chosen for a wind tunnel.
In both the wind tunnel and real site location, the recommended distances to
form the solution domain are 2–5 times the rotor diameter as the upstream
distance, 10–20 times as the downstream distance, and 1–3 times from the
edge as the distance from both sides. Mesh refinement is done near all
surfaces to allow finer resolution of the flow pattern to be developed within
the boundary layer. In particular, special attention on mesh refinement is
required near the leading and trailing edges and the tip of the blades.
Since wind turbines are a type of rotating machinery, a numerical
simulation technique using a rotating frame of reference is appropriate
Interface zone 1
Interface zone 2
Figure 6.4. Meshing for entire wind turbine system CFD simulation.
Source: Reprinted from Moa et al. [12].
NUMERICAL SIMULATIONS OF SMALL WIND TURBINES • 141
to simulate the moving blades and rotor. The fluid inside this rotating
domain is considered to be moving with the frame, and this simplifies the
computation and allows faster completion of the calculation. As shown in
Figure 6.4, the surfaces of blades and rotor are meshed in the rectangle
part and the rectangle part is inside the cylinder part, which causes the full
rotation of the blades and rotor. Both the rectangle and cylinder parts are
considered moving domains, whereas the wind tunnel part is treated as a
stationary domain. The sliding mesh technique is employed and the fluid
interfaces are defined between the aforementioned parts. Hence, interface
zone 1 is located between the rectangle and cylinder parts, and interface
zone 2 is located between the cylinder and wind tunnel parts. Although
this technique may be computationally demanding, it is the known method
for achieving high accuracy unsteady solutions for moving and stationary
domains that comprise multiple frames of reference.
The boundary conditions are specified by the standard, and the inputs
described in the previous section are to be translated as the boundary
conditions. Beyond the conditions from product certification and code
compliance, conditions from the actual site location can also be used in
the simulation.
Focus is now turned back to CFD simulation to further determine the next
step for completing the analysis. The wind blows in open air; therefore,
most flow situations are turbulent flows. Due to the complexity of the
different flow patterns, numerous models are developed to balance the
accuracy and speed of the computation. In addition to flow modeling, sim-
ulation with multi-physics is becoming more common for achieving bet-
ter designs. One combined model worth mentioning is the fluid–structure
interaction (FSI) model, which is discussed later.
flow patterns to more generalized codes that cover wide applications. The
guiding principles of formulating the CFD code include the law of con-
servation of mass and Newton’s second law of motion. In the specialized
literature for fluid mechanics, Newton’s second law for flowing fluid is
sometimes called momentum conservation or, alternatively, the Navier–
Stokes equations. Conservation of mass states that mass cannot be created
or destroyed, while Newton’s second law balances forces with changes
of momentum. The following list of turbulent flow models is available
through most commonly known software:
6.3.8 SOFTWARE
0.5
-0.5
-1
0 0.2 0.4 0.6 0.8 1
Y/c
Figure 6.5. Pressure coefficient (Cp) on the blade surface along the entire blade
length.
Source: Reproduced from Gómez-Iradi, Barakos, and Munduate [13].
12
10
Power generation (kW)
8
Model #1
6
Model #2
4 Model #3
0
0 5 10 15 20 25 30
Wind speed at hub height (m/s)
Table 6.2. Suggested format to present load results for extreme (top) and
fatigue (bottom) loads.
Results of the extreme load evaluation at a component
Design Safety Fx Fy Fz Fres Mx My Mz Mres
Load factor [kN] [kN] [kN] [kN] [kNm] [kNm] [kNm] [kNm]
case
DLC with
Fx Min
DLC with
Fx Max
DLC with
Fy Min
DLC with
Fy Max
DLC with
Fz Min
DLC with
Fz Max
DLC with
Fres
Min
DLC with
Fres
Max
DLC with
Mx Min
DLC with
Mx Max
DLC with
My Min
DLC with
My Max
DLC with
Mz Min
DLC with
Mz Max
DLC with
Mres
Min
DLC with
Mres
Max
NUMERICAL SIMULATIONS OF SMALL WIND TURBINES • 147
for this process. Therefore, careful planning for each of the aforemen-
tioned steps is recommended.
6.6 ACKNOWLEDGMENTS
The author would like to express appreciation to Janet Pau and Wai Heng
Ken Cheng for their support for this chapter. In addition, the author would
like to thank professor Ephraim Sparrow for his friendship and inspiration
for professional excellence.
REFERENCES