Next Article in Journal
A Study of Factors Affecting National Energy Efficiency
Previous Article in Journal
Bibliometric Review and Technical Summary of PWR Small Modular Reactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Blade Design and Aerodynamic Performance Analysis of a 20 MW Wind Turbine for LCoE Reduction

1
Department of Aerospace Engineering, Graduate School, Jeonbuk National University, 567 Baekje-daero, Deokjin-gu, Jeonju 54896, Jeonbuk, Republic of Korea
2
Department of Aerospace Engineering, Jeonbuk National University, 567 Baekje-daero, Deokjin-gu, Jeonju 54896, Jeonbuk, Republic of Korea
*
Author to whom correspondence should be addressed.
Energies 2023, 16(13), 5169; https://doi.org/10.3390/en16135169
Submission received: 27 May 2023 / Revised: 25 June 2023 / Accepted: 3 July 2023 / Published: 5 July 2023
(This article belongs to the Topic Advances in Wind Energy Technology)

Abstract

:
The aim of this study is to develop a blade mass model that incorporates a low-induction rotor (LIR) and a low-specific power concept to reduce aerodynamic loads and lower the Levelized Cost of Energy (LCoE). This blade mass model replaces the traditional simple scaling rule and incorporates the concept of LCoE reduction, presenting not only the mass distribution in the blade span direction but also the stiffness distribution. In order to achieve the desired reduction in LCoE, we developed a mathematical model that expresses blade mass as a function of the axial induction factor, which influences the aerodynamic load on the blade. We used this model to determine geometries of various low-induction rotors for 20 MW class horizontal axis wind turbine, and to identify the axial induction factor that correlates with the lowest blade mass. The chord length and twist angle in the spanwise direction of the blade were determined using PROPID’s reverse design process, based on the specified axial induction factor. Since the low-induction concept is not aerodynamically optimal, a low-specific power design approach was also adopted. This involved increasing the blade length and shifting the power curve to the left. By doing so, the AEP is increased, directly contributing to a reduction in the LCoE. Mass per unit length of the blade was presented, reflecting the distribution of airfoil type, blade geometry, and shapes of internal structures such as spars and webs.

1. Introduction

The global need for new research and development of large-scale wind turbines is growing. New wind turbines will help lower the Levelized Cost of Energy (LCoE) [1], and efficiently harness the powerful and renewable energy source that is wind. Wind turbine manufacturers such as Vestas, Siemens-Gamesa, GE, and Mingyang have already commercialized or are preparing to commercialize 12 to 16 MW wind turbines, and national institutions are also actively researching the topic. The US National Renewable Energy Laboratory (NREL), for example, announced a 5 MW reference wind turbine (RWT) [2] in 2009. Although the NREL 5 MW wind turbine was not intended for manufacture, it was used to provide input data for FAST software version 7 [3] (an aeroelastic analysis tool for wind turbines) and served as a widely known reference among wind turbine manufacturers as well as researchers. The IEA presented a 15 MW class wind turbine [4], which was scaled-up based on a DTU 10 MW reference wind turbine [5], and also presented related aeroelastic analysis results using HAWC2 software version 12.8 [6] developed by DTU. The DTU 10 MW system adopted a medium speed multi-stage gearbox with a maximum tip speed of 90 m/s, whilst the IEA 15 MW system adopted a direct drive and had maximum tip speed of 95 m/s. The two systems were identical in their other parameters, as both belonged to the FFA-W3 airfoil series, IEC-1B turbine class, and possessed the same upwind-type rotor orientation and three-bladed rotor.
Reference [5] reports a 3-bladed 20 MW upwind rotor system based on the UpWind 5-MW reference turbine (version 8). The design employed the ECN’s PHATAS code for aeroelastic load calculation consistent with IEC61400-1 edition 2 Class II B [7]. Blade shape was determined using the classically up-scaled wind turbine geometry in BOT [8]. High Reynolds number aerodynamic performance data up to 20 million was investigated using RFOIL [9]. Other 20 MW wind turbines have been developed and presented through UpWind, INNWIND.EU, and other projects for reference wind turbines [10,11]. Furthermore, a phased study for the development of a 50 MW class wind turbine system is in progress in the United States of America [12,13,14,15].
The INNWIND.EU project of the European Union has developed innovative rotor designs, drivetrain components, and fixed and floating substructures that significantly reduce the LCoE (levelized cost of energy) while increasing the efficiency of 10 and 20 MW offshore wind turbines [16]. One of the project’s most important innovations is the low axial induction rotor (LIR), which captures more energy while constraining extreme loads at the blade root and large rotor diameters [16,17]. During the blade design stage, the amount of material used in the blade is the most significant cost determinant. Scaling up a 5 MW class wind turbine to a 20 MW class turbine results in a doubling of blade length, but an eight-fold increase in mass [1,18]. Unfortunately, the new generation of improved materials used for wind turbine manufacturing has not kept pace with the speed of load increase associated with larger wind turbines. To reinforce the stiffness of the blades, using more material may be necessary, which increases the LCoE. Therefore, an additional consideration in the design of very large blades, such as those in the 20 MW class, is how to reduce loads while meeting power curve requirements. To address these issues, the concept of low-induction rotor design has been introduced. The axial induction factor is the fractional decrease in axial wind speed between the far upstream and rotor plane, and the lower the axial induction factor, the less thrust is produced. The INNWIND.EU’s 20 MW LIR blade is 13% longer and 7.6% heavier than the blade of the 20 MW reference wind turbine developed through the same project [17]. Ultimately, the LCoE was reduced by 4% using this system instead of the 20 MW RWT [17].
One study attempted to obtain the lowest LCoE possible by changing the number of blades, tip speed ratio, chord length distribution, axial induction factor, and lift coefficient of 20 MW or higher-class reference wind turbines. Clearly, achieving a lower LCoE is a very important goal, even at the preliminary design stage. While estimating annual energy production during this stage is not difficult, estimating the production costs, a component of LCoE, requires a study. LCoE can be divided into two topics: wind turbines and wind farms. Shafiee et al. [19] developed a parametric whole life cycle cost model to identify the key cost drivers of offshore wind projects and parameters that significantly influence the LCOE. The proposed model was tested on a 500-MW offshore baseline wind farm project, and the results were compared with experimental findings reported in the literature. It is asserted that the proposed model can help evaluate project performance and reduce costs. Griffith et al. [20] studied blade manufacturing costs for the Sandia 100 m wind turbine blade using the Sandia Blade Manufacturing Cost Tool (version 1.0). They conducted sensitivity studies as examples to demonstrate the potential use of the tool for cost tradeoff analyses between materials, labor content, and equipment components involved in blade manufacturing. Ashuri [21] and Ashuri et al. [22] presented a method for multidisciplinary design analysis optimization (MDAO) of large-scale wind turbines. In their referenced work, they optimized the NREL 5 MW reference wind turbine and scaled it up to 10 MW and 20 MW to evaluate the effect on LCoE. Rotor and tower design parameters were optimized with the goal of reducing LCoE. In particular, the results obtained from the study of Reference [21] were used to identify the scaling rule tendency. Serafeim et al. [23] proposed an MDAO approach to reduce the LCoE of the DTU-10 MW Reference Wind Turbine. A cost model for the entire wind turbine was implemented by combining existing models from the literature with open data. The model considers the costs of composites, resin, adhesive, paint, bolts, lightning protection, as well as labor and other manufacturing process expenses. In addition, Bortolotti et al. [24] presented several other MDAO-based frameworks. In their study, they described a comprehensive blade cost model for wind turbine blades ranging from 30 to 100 m in length. The proposed cost model was applied to three specific blades: the WindPACT blade, the IEA 3.4 MW wind turbine blade, and the SNL-100-03 blade. Material costs accounted for 45–70% of the total cost of these three blades [20,24].
In this study, the LCoE calculation focuses specifically on wind turbines. The cost elements of a wind turbine, as identified in Fingersh et al.’s study [1], include the rotor, tower, drive train and nacelle, control and safety system, and balance of station. Additionally, for offshore wind turbines, there are supplementary costs associated with factors such as marinization, scour protection, and port and staging equipment. The most important elements of a wind turbine are the blade diameter and tower height. Cost of blade and tower account for 30% of the total cost [22].
The reasons for presenting the blade mass model in this study and the need for it are as follows. Firstly, the proposed blade mass model replaces the traditional simple scaling rule [1,21] and incorporates the concept of reducing LCoE. It can provide not only mass distribution along the blade span but also the stiffness distribution. A low-induction rotor and a low-specific power [25] concepts were adopted to reduce aerodynamic loads and lower the LCoE. The low-induction rotor design reduces the aerodynamic force exerted on the blade by minimizing the axial induction factor, resulting in a decreased amount of material required for blade fabrication. Additionally, the low-induction factor design compensates for lower power by increasing the blade length, achieving a low specific power and enabling the generation of rated power at lower wind speeds. This adjustment shifts the P-V curve to the left, ultimately increasing annual power generation. Moreover, by finding the axial induction factor that minimizes mass while satisfying the given P-V curve, the chord length and twist angle of the blade can be determined. The reduction in blade mass and the subsequent increase in annual energy production directly contribute to a lower LCoE, not only for the designed blade but also for the entire wind turbine system. Blade shapes, performance curves, and flapwise bending moments for various axial induction factors at the blade root are compared. Finally, the mass model of the blade proposed in this study can be used to obtain blade data for comprehensive aeroelastic analysis tools such as FAST [3], including mass and stiffness distributions along the span direction. It is important to note that these blade data already reflected the concept of reducing the LCoE of the wind turbine in an earlier step of the MDAO process.

2. The Idea of the Low-Induction Rotor

The relationship between the mass model and the axial induction factor was determined using blade element momentum theorem [26]. In this theorem, a one-dimensional incompressible steady potential flow is assumed. The rotor and its surrounding flow were briefly modeled with a stream tube (Figure 1). The axial induction factor a = ( V 0 u ) / V 0 was defined using the free stream wind speed V 0 and wind speed u at the rotor plane. The power coefficient and thrust coefficient, which are the main coefficients of the wind turbine, were expressed as the axial induction factor. An optimal aerodynamic design aims to have an axial induction factor of a = 1 / 3 with a maximum power coefficient. As shown in Figure 2, as the axial induction factor was reduced from 1/3 to 1/4, thrust decreased by 16%, while power decreased by only 5%, i.e., thrust materially decreased whilst the power decreased only slightly. The low-induction rotor design, which reduces the aerodynamic loads acting on the blade, reflects this idea. Equations (1) and (2) are the power and thrust coefficients for one-dimensional axisymmetric flow, expressed as an axial induction factor using the 1D momentum theorem:
C P = 4 a ( 1 a ) 2
C T = 4 a ( 1 a )
This time, we want to determine the thrust acting on the blade using the blade element momentum theorem. Figure 3 illustrates the aerodynamic loads acting on a blade section located at a distance of r from the rotating axis and rotating at a speed of ω. Relative wind V r e l brought to the cross section of the wind turbine blade and the components of the decomposed forces. D and L represent drag and lift, while T and H represent the results of decomposition into the thrust in a direction perpendicular to the rotor plane and the tangential force applied to the rotor plane. H directly affects power by inducing torque, while T affects the bending moment of the blade. All forces were assessed in terms of [Newton/m]. According to the momentum theorem, power is a function of axial and tangential induction factors, but thrust is a function of only the axial induction factor. This is expressed as the following equation:
T = 4 π ρ V 0 2 a ( 1 a ) r / B
in which B denotes the number of blades. The reduced power, a function of the low-induction rotor design, was compensated for by increasing the length of the blade. All blade designs were compared and analyzed based on the same power curve assumptions.

3. Blade Mass Modeling to Reduce Turbine LCoE

3.1. LCoE of Wind Turbine

LCoE is expressed as an index that describes the economic feasibility of energy production and refers to the cost required to produce 1 kWh of energy. In wind turbine systems, LCoE is expressed as a function that incorporates initial installation and operating costs, as well as annual energy production, as shown in the following Equation [27]:
L C o E = C a p E x × F C R + O p E x A E P
where CapEx is the initial capital costs of fixed assets, FCR is the fixed charge rate, and OpEx is cost during their operation. FCR accounts for financial factors such as the economic growth rate and depreciation across the period of wind turbine operation. According to the traditional scaling rule, the denominator in Equation (4) is proportional to the square of the blade radius, while the numerator is proportional to the cube of the blade radius (blade mass). This discrepancy can pose a challenge when attempting to reduce LCoE. Therefore, it is crucial to minimize the mass as the blades increase in size.
In this study, we focused on the blade material costs in determining CapEx because material costs accounted for 50–75% of the total cost of blades [20,24], as shown in Table 1. All other costs except for blade material assumed to be the same. In Section 4.2.2, we will compare the LCoE between turbines with the same AEP but different blade lengths. Assuming that the remaining variables in Equation (4) are identical, except for the blade cost, we can express the LCoE ratio between the turbines as follows [27]:
β = C o s t b l a d e × F C R C a p E x × F C R + O p E x
L C o E L c o E r e f = 1 ( 1 C o s t b l a d e C o s t b l a d e , r e f ) β
in which β expresses the fraction of the blade cost over the total cost of a wind turbine system. The subscript ‘ref’ in Equation (6) denotes a selected blade among those with different lengths. In this study, we will choose the blade that minimizes the LCoE. C o s t b l a d e represents total cost of blade including materials, labor, tooling, and other related expenses [24]. The right-hand side of Equation (6) contains the relative proportion of the cost to the reference blade. While it is difficult to estimate the absolute cost of a blade, an increase or decrease in the relative cost to the reference blade is easy to express in a formula.

3.2. Blade Mass Modeling

Material costs constitute the largest component of blade costs, followed by labor costs. Other expenses encompass overhead, capital, tooling, utility, and maintenance costs [24]. As the blade cost represents approximately 3–28% of the total turbine cost [1,11,21,28] as shown in Table 2, with material costs accounting for approximately 50–75% of the blade cost [20,24] as shown in Table 1, the contribution of blade materials to the overall turbine cost is estimated to be around 1.5–21%. The cost of blades typically accounts for around 10% of the total turbine cost. Additionally, implementing load reduction via the low-induction rotor concept and enhancing AEP through the low specific power design concept can decrease the LCoE associated with the turbine. We developed a blade mass ratio model for horizontal axis wind turbine to compare the LCoE of wind turbines with different blade numbers or axial induction factors. The cost of manufacturing, processing, and transporting each blade was not considered, only the cost of the material used to create each blade. To simplify the cost calculation process, first, each blade is divided into several blade elements of a certain length (Figure 4). Any blade element i and baseline blade element i b a s e that were compared were composed of the same materials at the same ratio. As blade elements were divided by the same length, the cost of each blade was calculated as the sum of the cost of each blade element. The ratio expressing the cost of a designed blade compared to the cost of the baseline blade was expressed as:
C o s t b l a d e C o s t b l a d e , b a s e = 1 + α ( C o s t b l a d e ,   m a t e r i a l C o s t b l a d e ,   m a t e r i a l , b a s e 1 )
M a s s b l a d e M a s s b l a d e , b a s e = ( c o s t   p e r   k g ) ( c o s t   p e r   k g ) ( d e n s i t y ) ( d e n s i t y ) R R b a s e 1 i A r e a i 1 i A r e a i , b a s e ( = C o s t b l a d e ,   m a t C o s t b l a d e ,   m a t , b a s e )
in which Area describe the blade sectional area of the blade element, and R denotes the tip radius of the rotor blade, respectively. The subscripts ‘base’ and ‘mat’ in Equations (7) and (8) denote a baseline turbine (NREL 5 MW reference wind turbine in this study) and material, respectively. Each blade mass ratio reflects blade length and the cross-sectional area of the material used at each blade element. Each blade element can be divided into skins, spars, and webs. In Equation (7), α represents the percentage of material cost in the total cost of the blade and specifically ranges from 50% to 75% [20,24]. To focus on the blade mass only, we convert the blade cost ratio presented in Equation (7) to the blade mass ratio, as shown in Equation (8).
As shown in Figure 5, skin was modeled simply as a hollow elliptical tube, while spars and webs were similarly modeled as a hollow rectangular tube. The hollow elliptical tube can be described in terms of chord length and airfoil thickness, which must be determined in order to perform an aerodynamic analysis of each blade, while skin thickness t s k i n is an unknown variable. In our model blade section (Figure 5), the thickness of the spar and the thickness of the web were assumed to be identical ( t s p a r = t w e b ). Web height, distance, and thickness were calculated as airfoil thickness ratio τ multiplied by 0.9 [29], chord length multiplied by w (with w = 0.3 ) [30], and skin thickness multiplied by γ ( γ = 1.5 ), respectively. Through these dimensional determinations, the cross-sectional area of the blade structure was expressed using Equation (9). The parameter f A was then derived as a function of w ,   γ ,   and   τ as shown in Equation (10). Ratios between airfoil thickness and chord length, i.e., airfoil thickness ratios, for all airfoils used for the wind turbine blade design are shown in the third column of Table 3. P e / C in Table 3 expresses the perimeter length divided by the chord length.
The mass ratio between the two blades can be expressed as a function of chord length and skin thickness, as shown in Equation (11). The right columns of Table 3 present the perimeter of an ellipse with the same thickness ratio as the corresponding airfoil perimeter, along with the relative error between the two values. The perimeter plays a crucial role in determining the mass and moment of inertia per unit length of the blade cross section. According to the results in Table 3, the relative error in the perimeter for both geometries is no more than 9%. Therefore, substituting the airfoil with an ellipse does not introduce a significant error. On the other hand, utilizing the ellipse simplification in the blade mass model offers the advantage of a highly concise expression, as shown in Equation (11).
A r e a i = C i t s k i n , i f A
f A = π 2 ( 1 + τ i ) + 2 γ ( w i + 0.9 τ i )
M a s s b l a d e M a s s b l a d e , b a s e = R R r e f B B r e f 1 i C i t s k i n , i f A 1 i C i , r e f t s k i n , i , b a s e f A , b a s e
in which sectional parameters are presented in Figure 5.
Of the various loads acting on the rotor, we considered the flapwise bending moment at the blade root [31], which is directly affected by thrust acting on the rotor plane. Thus, we modelled the blade as a simple cantilever beam (Figure 6). We assumed that the bending stresses at the root of both the designed and the baseline blades were identical. Chord length and airfoil thickness at each blade section are known values, and the moment of inertia can be obtained using shape modeling (Figure 5). We assumed that each structural component—spar, skin, and web—resisted bending stress, similar to the wings of high-speed aircraft and very large wind turbine blades. Equation (16) was then obtained by rearranging Equations (12)–(15) into the equation of thickness:
σ = σ b a s e
σ = M τ i C i / 2 I
I i = 1 8 τ i 2 C i 3 t s k i n , i f I
f I = π 4 ( 3 + τ i ) + 27 25 γ ( 3 w i + 0.9 τ i )
t s k i n , i = t s k i n , i , b a s e M M b a s e τ i , b a s e C i , b a s e 2 f I , b a s e τ i C i 2 f I
in which I and M represent the area moment of inertia and flapwise bending moment at the root of the blade, while f I denotes the coefficient of the area moment of inertia.
The flapwise bending moment at the blade root is caused by aerodynamic thrust, as confirmed through Equation (3). The flapwise bending moment along the blade span can be expressed through Equation (17), which incorporates our simplified cantilever beam assumption. The force acting on the cantilever beam can be expressed as a triangular distributed load along the spanwise direction (Figure 6). F n is the force per unit length acting on the end of the cantilever beam. Equation (19) can be obtained from Equations (11), (13), (17) and (18). Equation (19) shows that the ratio of the mass of the designed blade to the mass of the baseline blade (NREL 5 MW blade) is expressed as a function of airfoil thickness, chord length, axial induction, and blade length. According to the numerical simulation, the effects of the γ and w i on the mass ratio determined through Equation (19) are negligible.
M = F n 6 R ( 2 R 3 + 3 R 2 r r 3 )
F n = 4 ρ π V 0 2 a ( 1 a ) r / B
M a s s b l a d e M a s s b l a d e , b a s e = a ( 1 a ) a r e f ( 1 a b a s e ) R 4 R b a s e 4 1 i C i , r e f t s k i n , i , b a s e f A τ i , b a s e C i , b a s e f I , b a s e τ i C i f I 1 i C i , b a s e t s k i n , i , b a s e f A , b a s e
in which r shows the spanwise location of the blade element from the center of the rotor hub.

4. Results and Discussion

4.1. Considerations for Aerodynamic Rotor Design

4.1.1. Aerodynamic Data of Airfoils

One of differences between the NREL 5 MW wind turbine and our 20 MW wind turbine is the distribution of airfoils along the blade. In particular, the NACA64-618 airfoil was not applied to the outboard section of our blades, as the longer blade length of the low-induction rotor required increased stiffness at the tip to prevent excessive twist or deflection. Instead, we used a DU93_W_210LM airfoil at the outboard section of the blade, as shown in Table 3 [32,33].
To design a 20 MW class wind turbine rotor, the NREL 5 MW wind turbine was used as a baseline turbine, and the need for an increased Reynolds number consistent with an increase in blade chord length and rotational speed was considered for aerodynamic data correction. The Reynolds number of the 20 MW class wind turbine blade is approximately double (around Re   = 20 × 10 6 ) that of the NREL 5 MW wind turbine. This value is high enough to influence aerodynamic performance, and lift and drag coefficients were revised to account for the airfoil’s angle of attack at a high Reynolds number accordingly [33]. The results were used as the airfoil input data for the design of a 20 MW class low-induction rotor.

4.1.2. Low-Specific Power Design

NREL 5 MW and DTU 10 MW reference wind turbines have some known shortcomings, including a too high specific power of 401 W/m2 (Table 4). In 2016, the IEA performed a survey to identify the use cases and need for reference wind turbines for research and development applications [4]. Nearly 58% of the respondents suggested a specific power lower than 350 W/m2. The specific power of the IEA 15 MW reference wind turbine developed by NREL, DTU, and the University of Maine [4] was 332 W/m2. Table 4 shows the specific powers of various reference wind turbines, including the NREL 5 MW, LEANWIND 8 MW [34], DTU 10 MW, AVATAR 10 MW [35], IEA 15 MW, ECN 20 MW, and INNWIND.EU 20 MW RWTs. Table 4 also presents relevant information concerning three commercial turbines: the GE Heliade-X 12 MW, and two Korean manufacturer’s turbines. Notably, the two Korean system manufacturers plan to release wind turbines that reflect low specific power designs in accommodation of the wind conditions around the Korean Peninsula. Based on a wind resource assessment of measured wind data around the Korean Peninsula, the wind class around the sea areas of the Korean Peninsula at which fixed sub-structure type offshore wind turbines can be appropriately installed is 3 or more [36]. In follows that turbines with a lowered specific power design are suitable for the West Sea area of the Korean Peninsula. Fundamentally, the low-induction rotor is attractive in low wind speed areas because it has the property of producing a lower specific power, which makes it more efficient in such conditions.

4.1.3. Basic Design Parameters for 20 MW Wind Turbine

PROPID [37] is a wind turbine aerodynamic design tool that allows for an inverse design process in which desired variables, such as axial induction factor or lift coefficient along the blade span, are inputted first. Therefore, using this software, unlike an optimal aerodynamic design in which an axial induction factor is set to a = 1 / 3 , a user can set the desired axial induction factor to obtain an appropriate blade geometry. We set axial induction factor values at a = 0.2   ~   0.3 for our low-induction rotor, and the value of the axial induction factor was set to be constant at each blade section (except for the inboard cylinder part of the blade). The basic parameters of our low-induction rotor for use in a 20 MW class wind turbine are shown in Table 5.
The aerodynamic design process for our 20 MW RWT proceeded with a rated wind speed of 10 m/s. The tip speed ratio of the NREL 5 MW wind turbine is 7, while that of our 20 MW wind turbine is as high as 9. As the 20 MW wind turbine is sufficiently far from the seashore, concerns over noise impacts are not acute, and aerodynamic performance can be improved by increasing the rotational speed (or Reynolds number).
Additional aerodynamic analysis of a rotor for a 20 MW wind turbine was performed using FAST software to examine the aerodynamic design of the blade [3]. Prior to this study, the results of the aerodynamic analyses of two tools, PROPID and FAST, had been cross-validated [38], and we likewise obtained the same numerical results of the C P λ curves (Figure 7) of an aerodynamically optimized blade for a 20 MW wind turbine rotor.

4.2. Low LCoE Design for 20 MW Wind Turbine Rotor

4.2.1. Validation of Mass Model

To compare the blade mass using Equation (19), blade skin thickness had to be determined. However, since the skin thickness of the 20 MW class low-induction rotor is unknown, Equation (19) was used to express the ratio of blade mass. In this manner, the mass of each low-induction rotor could be simultaneously compared. In other words, although the absolute mass of each blade was unknown, the relative mass of each was expressed by Equation (20). We assumed that the web distance w and the ratio γ of the skin of the corresponding blade element of the 20 MW class low-induction blade had the same value as the NREL 5 MW wind turbine blade. It is assumed that the normalized positions of the airfoils used in the blades targeted in this study are the same.
M a s s b l a d e M a s s b l a d e , c a s e 1 = a ( 1 a ) R 4 C i , 5 M W t s k i n , i , 5 M W f A τ i , 5 M W C i , 5 M W f I , i , 5 M W τ i C i f I , i a c a s e 1 ( 1 a c a s e 1 ) R c a s e 1 4 C i , 5 M W t s k i n , i , 5 M W f A , c a s e 1 τ i , 5 M W C i , 5 M W f I , i , 5 M W τ i , c a s e 1 C i , c a s e 1 f I , i , c a s e 1
In Equation (20), the subscripts “Case 1” and “5 MW” denote the values of the blade parameters for Case 1 and the NREL 5 MW blade [2], respectively. These values are used to estimate the mass distribution of the wind turbine blade using the mass model presented in this study.
By utilizing Equation (9), the mass at each node ‘i’ of the wind turbine blade can be expressed as Equation (21). Furthermore, Equation (21) can be transformed into Equation (22) through the use of Equations (16)–(18). This allows for the estimation of the mass at each node ‘i’ of an arbitrary wind turbine blade, leveraging the known information from the NREL 5 MW of wind turbine blades (baseline blades).
M a s s b l a d e , i = C i t s k i n , i f A , i Δ R i ( d e n s i t y i )
M a s s b l a d e , i = C i , 5 M W t s k i n , i , 5 M W f A , i Δ R i ( d e n s i t y i , 5 M W ) a ( 1 a ) R 4 a 5 M W ( 1 a 5 M W ) R 5 M W 4 τ i , 5 M W C i , 5 M W f I , i , 5 M W τ i C i f I , i
To validate Equation (22), we used two wind turbine rotors proposed by different sources: a 3.4 MW rotor proposed by the IEA [39] and a 10 MW reference rotor proposed by DTU [5]. We estimated the mass distributions of both blades using the mass model, scaling rule and compared the results, as shown in Figure 8 and Figure 9. The simple scaling rule shows the largest error with the blade mass distribution. Based on the comparison results, we can conclude that the mass distribution modeling used in this study is reliable. The blade masses for several reference wind turbines and the masses calculated using Equation (22) are presented in Table 6. It is evident that the presence or absence of MDAO has an impact on the difference between the two blade masses. Most of the blades that underwent MDAO are lighter compared with the masses obtained using Equation (22) for those that did not. In particular, the SNL 100-03 blade, which is derived from the SNL 100-00 blade and utilizes lighter materials through stepwise optimization, exhibits a significant improvement in mass reduction.
Figure 10 illustrates the non-dimensional mass distribution of the NREL 5 MW, IEA 3.4 MW, and DTU 10 MW wind turbine blades. The length and mass of each blade element are non-dimensionalized by the blade radius and blade mass, respectively. The outboard section of the blades exhibits good agreement, while discrepancies arise in the inboard section. These differences suggest variations in material density within the inboard blade elements or differences in the non-dimensional dimensions of structures such as the skin, spar, and web (see the dimensions of the structure of Figure 5 in Section 3.2). These disparities can manifest during the MDAO process, leading to observable variations in the mass distribution, particularly in the inboard part of the blade, as shown in Figure 8 and Figure 9.

4.2.2. Axial Induction Factor for Lowest LCoE

Two low-induction rotors were selected as comparators to confirm the difference according to the reference for comparison. The aerodynamically optimal 20 MW rotor with a = 0.3 is referred to as Case 1, which is close to a = 1 / 3 . The first low-induction rotor with a = 0.2 is referred to as Case 2. The blade corresponding to the lowest mass is named Case 3. The specifications for these blades are shown in Table 7. If the value of the axial induction factor corresponding to the lowest value of the blade mass ratio changes depending on the value of the denominator in the left-hand-side term of Equation (20), it would suggest that there is a problem in determining the axial induction factor with the lowest blade mass model. To confirm this, the variation of the blade mass ratios according to the axial induction factor was determined as the denominator of Equation (20) was changed to corresponding blade mass of Case 1 and Case 2 rotors. The results confirmed that there were no problems in determining the lowest mass axial induction factor as shown in Figure 10. Notably, the blade mass was lowest at a = 0.26 .
We then checked whether the change of β in Equation (5) affected the selection of the axial induction factor corresponding to the lowest LCoE. The range of the β generally being between 3% and 27% [1,11,21,28], the LCoE ratio of the low-induction rotor was compared for β = 0.05, 0.1, 0.15 and 0.2 as shown Figure 11. As we already saw in Figure 12, the axial induction coefficient for the lowest LCoE was still 0.26, as shown in Figure 11, suggesting that variations in β also did not affect the determination of which low-induction rotor resulted in the lowest LCoE. Accordingly, Case 3 includes the rotor with the lowest LCoE, as designed by an axial induction factor of 0.26. The mass of the Case 3 rotor is 8% less than the Case 1 rotor based on Figure 12 and Table 7.
Table 8 presents a comparison of the LCoE results for the 20 MW wind turbine studied by Ashuri et al. [11] and the Case 1, Case 2, and Case 3 rotors in this study. The mass and cost models employed in this study utilize the scaling rule to estimate the mass and cost of individual wind turbine components. Specifically, the blade mass and cost model can be expressed as Equations (23) and (24), while the tower mass and cost model can be described by Equations (25) and (26) [1,21].
M a s s b l a d e = 0.1452 × R 2.9158
C o s t b l a d e = ( 13.084 × M a s s b l a d e 4452.2 ) × ( 1 + PPI )
M a s s t o w e r = 0.697 × ( s w e p t   a r e a ) × ( h u b   h e i g h t ) 1414
C o s t t o w e r = 1.5 × M a s s t o w e r × ( 1 + PPI )
The PPI in Equations (24) and (26) is an index that is updated on a monthly basis to track the changes of costs of products and materials over a wide range of industries and industrial products [21].
To calculate the LCoE in Table 8, we use the blade masses of Case 1, Case 2, and Case 3 calculated by the mass model in this study. We can calculate the cost of pitch system, hub cone, main shaft, shaft bearing, yaw system, nacelle frame, nacelle railing as a function of blade length, and the cost of Gearbox, Electronics, Nacelle cover, turbine connection, cooling and the hydraulics system, Monitoring and safety system as a function of wind turbine rotor rated power. With the aforementioned assumptions, it is assumed that BOS (Balance of Station), LRC (Levelized Replacement Cost), M&O (Maintenance and Operation), IR (Interest Rate, FCR) are equal [11]. For AEP estimation in Case 1, Case 2, and Case 3, we employ the Weibull distribution with coefficients c = 9.47 and k = 2. Comparing the LCoE of the wind turbine system with the Case 1, Case 2, and Case 3 rotors, we observed that the Case 3 rotor shows a decrease in LCoE of approximately 1% compared to the Case 1 rotor, 6.8% compared to the Case 2, and 7% compared to the Ashuri et al. 20 MW rotor. In addition, the LCoE value in parentheses bottom of Table 8 is the LCoE when considering a wind turbine without a gearbox, and the Case 3 rotor has an LCoE reduction of about 1.2% compared with the Case 1 rotor, about 7.8% compared with the Case 2 rotor, and about 6.8% compared with the Ahuri et al. 20 MW rotor. Based on these results, it is confirmed that the Case 3 rotor, obtained by applying the mass model proposed in this study, achieves the lowest LCoE.

4.2.3. Case Study for Different Rotor Design

Figure 13 shows the power coefficient curves according to the tip speed ratio for various low-induction factors in cases in which a = 0.2 (Case 2 rotor in Table 7), 0.26 (Case 3 rotor in Table 7), and 0.3 (Case 1 rotor in Table 7), and confirms that the design axial induction factor affects the power coefficient curve. The axial induction factor a = 0.2 was chosen to maximizes the objective function of C P / C M 2 / 3 without losses, based on the study by Chaviaropoulos and Sieros [41]. The maximum value of the power coefficient of the Case 3 rotor is 4.1% smaller than the Case 1 rotor shown in Table 7. The reduction has a direct impact on annual energy production (AEP). Recalling from Equation (27) that the square of the rotor length and the power coefficient are inversely proportional, to equalize all AEPs for the low-induction rotors, the blade length of each low-induction rotor was adjusted to match Figure 14 or the power curve. The adjusted blade lengths were also applied to the blade mass ratio model shown in Equation (19).
P = C p 1 2 ρ A V 0 3
All blades were designed based on identical power curve assumptions. The power curves for Case 1, Case 2, and Case 3 rotors, as determined through PROPID, are shown in Figure 15 and all curves were confirmed to match well with each other. In this study, the 20 MW reference wind turbine model is a direct-drive type without a gear box, and the rated power, accounting for both the efficiencies of the generator and the power converter, is 21.4 MW at a rated wind speed of 10 m/s. This reflects a mechanical transmission efficiency of 99%, a generator efficiency of 96%, and a power conversion efficiency of 98.5% [32].
The shapes and specifications of the Case 1, Case 2, and Case 3 blades are shown in Figure 16 and Figure 17, and Table 7, respectively. Comparing each Case, we note that a decrease in the axial induction factor correlates with an increase in the twist angle (Figure 17). However, the Case 3 blade has the longest chord length among the low-induction rotors (Figure 16). Figure 18 highlights that the low-induction rotors, including Case 2 and Case 3, had a lower flapwise bending moment at the blade root than the Case 1 rotor. In particular, variation of the bending moment according to the wind speed at the blade roots for the Case 3 rotor which has both a = 0.26 and the lowest LCoE are not significantly different compared to the loads generated by the Case 1 rotor with a = 0.2. This means that the low-induction rotor design with a slightly higher axial induction factor does not compromise the structural integrity of the blade and can still meet the load requirements. This is beneficial because it allows for a reduction in the amount of material needed to reinforce the blade, which can lead to lower costs and a more efficient design. Finally, the mass distributions along the blade span of the Case 1, Case 2, and Case 3 blades are shown in Figure 19.
Figure 20 shows that the normalized LCoE ratio of the aerodynamically optimized two-bladed rotor is lower than that of the aerodynamically optimized three-bladed rotor (Case 1), at least when the rated wind speed and tip speed ratio for the rotors are the same (Table 5). Figure 21 highlights that the total mass of the two-bladed wind turbine rotor with the lowest LCoE is 27% less than that of the three-bladed Case 1 rotor. Considering that the ratio of blades to total wind turbine price is 0.1, a 27% decrease in blade cost corresponds to a 2.7% decrease in LCoE. However, the maximum value of the power coefficient for two-bladed rotor with the lowest LCoE is 0.4563 as shown in Figure 21; thus, the maximum value of the power coefficient for the two-bladed rotor is 2.6% less than that of the three-bladed rotor with the lowest LCoE (Case 3). As a general matter, the two bladed rotor, when compared to its three-bladed counterpart, does not have enough margin to avoid resonance between the tower frequency and the blade passage frequency [42]. The two-bladed rotor also experiences a larger tip loss effect, and would require a teetering hinge at the hub to reduce the loads due to the noted aerodynamic unbalance [42]. Moreover, the mass moment of inertia of the wind turbine rotor for the yawing axis is a function of the azimuth angle of the blade. We note that the operation of the two-bladed wind turbine rotor generates moment fluctuation during yawing motion. For these reasons, this study adopted the three-bladed rotor as the default RWT. Blade cost can be reduced to 8% less than that of the aerodynamically optimized blade.

5. Conclusions

We developed a blade mass model and implemented it in the design of a 20 MW offshore rotor, utilizing a low-induction and low-specific power approach. The primary objectives of this design effort are to decrease the aerodynamic load on the blades through the low-induction concept, enhance the annual energy production (AEP) through the low-specific power concept, reduce blade costs, and ultimately achieve a lower levelized cost of energy (LCoE) for the wind turbine. A mathematical formula to express the mass ratio, which is the function of the axial induction factor, a, and blade shape, was also presented. All of airfoils’ aerodynamic data were calibrated to fit the Reynolds number of the 20 MW wind turbine blade.
The mass ratio formula showed that the lowest LCoE for the 20 MW wind turbine was achievable at a = 0.26, and this result was consistent when the blade’s share of the total system price was changed. Both the blade geometry along the spanwise direction and the power coefficient Cp according to the tip-speed-ratio λ of three selected axial induction factors (a = 0.2, 0.26, 0.3) were compared using PROPID’s reverse design process. We confirmed that the bending moment of the low-induction rotor (a = 0.26) at the blade root was less than that of the aerodynamically optimized rotor (a = 0.3). However, the maximum power coefficient value of the low-induction rotor was 4.1% less than that of the optimized rotor. This reduced performance was compensated for by increasing blade length to satisfy the same power curve assumptions. Overall, the low-induction rotor at a = 0.26 achieved an 8% reduction in mass compared to the aerodynamically optimized rotor at a = 0.3. By evaluating the LCoE parameters of each turbine component, the LCoE of the three 20 MW turbines developed in this study was calculated and compared to the LCoE of another wind turbine of the same capacity. The results showed that the LCoE of the low-induction rotor turbine (Case 3) not only had the lowest value, but it was also 7% smaller than the LCoE of the fourth wind turbine used for comparison.
Although a two-bladed rotor reduces costs by 27% more than a three-bladed rotor, we selected a three-bladed rotor for the 20 MW reference wind turbine. This decision was made because the two-bladed rotor experiences a reduced resonance avoidance margin with the tower, moment fluctuations during yawing motion, a diminished maximum power coefficient, and requires a teetering hub.
Finally, distribution of mass per unit length along the blade span was presented for the three-bladed upwind type rotor, reflecting the airfoil types, blade geometry, and shapes of internal structures such as spars and webs. The devised blade mass ratio model has shown itself to be a suitable formula for the preliminary design of blade shapes and estimation of mass distribution along the blade span for use in ultra-large scale reference wind turbines. The mass model of the blade proposed in this study can be utilized to obtain blade data, including mass and stiffness distribution along the blade span. This blade data is valuable for conducting comprehensive analyses using aeroelastic tools. Unlike the mass distribution result obtained through a simple scaling rule, this blade data reflects the concept of reducing the LCoE of the wind turbine.

Author Contributions

Conceptualization, K.-H.J. and K.-W.R.; methodology, K.-H.J. and K.-W.R.; software, K.-H.J.; validation, K.-W.R.; formal analysis K.-H.J. and K.-W.R.; investigation, K.-H.J. and K.-W.R.; resources, K.-H.J. and K.-W.R.; data curation, K.-H.J. and K.-W.R.; writing—original draft preparation, K.-H.J.; writing—review and editing, K.-W.R.; visualization, K.-H.J. and K.-W.R.; supervision, K.-W.R.; project administration, K.-W.R.; funding acquisition, K.-W.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Korea Electric Power Corporation, grant number R21XO01-6 and supported by Research Base Construction Fund Support Program funded by Jeonbuk National University in 2021.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

aAxial induction factor
a’Tangential induction factor
BNumber of blades
CiChord length at the ith blade element
CPPower coefficient
CTThrust coefficient
F`nNormal aerodynamic force at a rotor plane
IArea moment of inertia
LCoELevelized cost of energy
MBending moment at radial position r
rRadial position from rotating axis
RRotor radius
tThickness of airfoil skin, spar, and web
uWind speed at rotor plane
V0Free stream wind speed
VrelEffective wind speed
wDistance between webs
αRatio of a blade material cost to the total blade cost of a wind turbine system
βRatio of a blade cost to the total cost of a wind turbine system
γCoefficient of airfoil web thickness to airfoil skin thickness
λTip speed ratio
ρAir density
σStress
σsolSolidity
τAirfoil thickness
φInflow angle

References

  1. Fingersh, L.; Hand, M.; Laxson, A. Wind Turbine Design Cost and Scaling Model; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2006. [Google Scholar]
  2. Jonkman, J.; Butterfield, S.; Musial, W.; Scott, G. Definition of a 5-MW Reference Wind Turbine for Offshore System Development; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2009. [Google Scholar]
  3. Jonkman, J.M.; Buhl, M.L., Jr. FAST User‘s Guide; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2005. [Google Scholar]
  4. Gaertner, E.; Rinker, J.; Sethuraman, L.; Zahle, F.; Anderson, B.; Barter, G.; Abbas, N.; Meng, F.; Bortolotti, P.; Skrzypinski, W.; et al. Definition of the IEA Wind 15-Megawatt Offshore Reference Wind Turbine; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2020. [Google Scholar]
  5. Bak, C.; Zahle, F.; Bitsche, R.; Kim, T.S.; Yde, A.; Henriksen, L.C.; Natarajan, A.; Hansen, M. Description of the DTU 10 MW Reference Wind Turbine; DTU Wind Energy Report-I-0092; Technical University of Denmark (DTU): Lyngby, Denmark, 2013. [Google Scholar]
  6. Larsen, T.J.; Hansen, A.M. How 2 HAWC2, the User’s Manual; Risø DTU National Laboratory: Roskilde, Denmark, 2019. [Google Scholar]
  7. IEC 61400-1; Wind Energy Generation Systems—Part 1: Design Requirements. IEC: Geneva, Switzerland, 2019.
  8. Bot, E.T.G.; Ceyhan, O. Blade Optimisation Tool User’s Manual; The Energy Research Centre of the Nethelands (ECN): Petten, The Netherlands, 2009. [Google Scholar]
  9. Montgomerie Jensen, B.O.G.; Brand, A.J.; Bosschers, J.; van Rooij, R.P.J.O.M. Three-Dimensional Effects in Stall; The Energy Research Centre of the Nethelands (ECN): Petten, The Netherlands, 1997. [Google Scholar]
  10. Peeringa, J.; Brood, R.; Ceyhan, O.; Engels, W.; de Winkel, G. UpWind 20MW Wind Turbine Pre-Design Blade Design and Control; The Energy Research Centre of the Nethelands (ECN): Petten, The Netherlands, 2011. [Google Scholar]
  11. Ashuri, T.; Martins, J.R.R.A.; Zaaijer, M.B.; van Kuik, G.A.M.; van Bussel, G.J.W. Aeroservoelastic Design Definition of a 20 MW Common Research and Turbine Model. Wind. Energy 2016, 19, 2071–2087. [Google Scholar] [CrossRef]
  12. Ananda, G.K.; Bansal, S.; Selig, M.S. Aerodynamic Design of the 13.2 MW SUMR-13i Wind Turbine Rotor. In Proceedings of the 2018 Wind Energy Symposium, Kissimmee, FL, USA, 8–12 January 2018. [Google Scholar] [CrossRef]
  13. Pao, L.Y.; Zalkind, D.S.; Griffith, D.T.; Chetan, M.; Selig, M.S.; Ananda, G.K.; Bay, C.J.; Stehly, T.; Loth, E. Control co-design of 13 MW Downwind Two-bladed Rotors to Achieve 25% Reduction in Levelized Cost of Energy. Annu. Rev. Control 2021, 51, 331–343. [Google Scholar] [CrossRef]
  14. Qin, C.; Loth, E.; Zalkind, D.S.; Pao, L.Y.; Yao, S.; Griffith, D.T.; Selig, M.S.; Damiani, R. Downwind Coning Concept Rotor for a 25 MW Offshore Wind Turbine. Renew. Energy 2020, 156, 314–327. [Google Scholar] [CrossRef]
  15. Yao, S.; Chetan, M.; Griffith, D.T.; Mendoza, A.S.E.; Selig, M.S.; Martin, D.; Kianbakht, S.; Johnson, K.; Loth, E. Aero-structural Design and Optimization of 50 MW Wind Turbine with over 250-m Blades. Wind. Eng. 2021, 46, 273–295. [Google Scholar] [CrossRef]
  16. Innovative Wind Conversion Systems (10–20 MW) for Offshore Applications; Final Report Summary—INNWIND.EU, 7th Framework Project, CORDIS; Eucopean Commission: Luxembourg, 2018; Available online: https://cordis.europa,eu/project/id/308974/reporting (accessed on 11 December 2018).
  17. Jensen, P.H.; Chaviaropoulos, T.; Natarajan, A. LCoE Reduction for the Next Generation Offshore Wind Turbine; Outcomes from the INNWIND.EU Project; European Commission: Luxembourg, 2017. [Google Scholar]
  18. Kesikbaş, O. Investigations of Upscaling Effects for Aerodynamic Design of Large Wind Turbine Rotors by Using Bem Theory and Optimization. Master’s Thesis, Middle East Technical University, Ankara, Turkiye, 2019. [Google Scholar]
  19. Shafiee, M.; Brennan, F.; Espinosa, I.A. A Parametric Whole Life Cost Model for Offshore Wind Farm. Int. J. Life Cycle Assess. 2016, 21, 961–975. [Google Scholar] [CrossRef] [Green Version]
  20. Griffith, D.T.; Johanns, W. Large Blade Manufacturing Cost Studies Using the Sandia Blade Manufacturing Cost Tool and Sandia 100-Meter Blades; SANDIA REPORT; Sandia National Laboratories: Albuquerque, NM, USA, 2013. [Google Scholar]
  21. Ashuri, T. Beyond Classical Upscaling: Integrated Aeroservoelastic Design and Optimization of Large Offshore Wind Turbines. Ph.D. Thesis, Delft University of Technology, Zuid-Holland, The Netherlands, 2012. [Google Scholar]
  22. Ashuri, T.; Zaaijer, M.B.; Martins, J.R.R.A.; van Bussel, G.J.W.; van Kuik, G.A.M. Multidisciplinary Design Optimization of Offshore Wind Turbines for Minimum Levelized Cost of Energy. Renew. Energy 2014, 68, 893–905. [Google Scholar] [CrossRef] [Green Version]
  23. Serafeim, G.; Manoas, D.; Riziotis, V.; Chaviaropoulos, P. Multidisciplinary Aeroelastic Optimization of a 10 MW-scale Wind Turbine Rotor Targeting to Reduced LCoE. J. Phys. Conf. Ser. 2022, 2265, 042051. [Google Scholar] [CrossRef]
  24. Bortolotti, P.; Berry, D.; Murray, R.; Gaertner, E.; Jenne, D.; Damiani, R.; Barter, G.; Dykes, K. A Detailed Wind Turbine Blade Cost Model; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2019. [Google Scholar]
  25. Kannan, B.; Jeyyaraj, B. Low Specific Power Wind Turbines for Reduced Levelized Cost of Energy; IntechOpen: London, UK, 2022. [Google Scholar]
  26. Hansen, M.O.L. Aerodynamics of Wind Turbine, 2nd ed.; EARTHSCAN: London, UK, 2008. [Google Scholar]
  27. Jang, K.H. 20 MW Class Low-Induction Rotor Shape Design and Analysis for LCoE Reduction. Master’s Thesis, Jeonbuk National University, Jeonju, Republic of Korea, 2023. [Google Scholar]
  28. Chaviaropoulos, P.K.; Chortis, D.; Lekou, D. Definition of the Reference Wind Turbine—Analysis of Rotor Design Paramete; INNWIND.EU Project, Deliverable 1.2.1; European Commission: Luxembourg, 2013. [Google Scholar]
  29. Giguère, P.; Selig, M.S.; Tangler, J.L. Blade Design Trade-Offs Using Low-Lift Airfoils for Stall-Regulated HAWTs. J. Sol. Energy Eng. 1999, 121, 217–223. [Google Scholar] [CrossRef] [Green Version]
  30. Giguère, P.; Selig, M.S. Blade Geometry Optimization for the Design of Wind Turbine rotors. In Proceedings of the 2000 ASME Wind Energy Symposium, AIAA-2000-0045, 2000 ASME Wind Energy Symposium, Reno, NV, USA, 10–13 January 2000. [Google Scholar] [CrossRef] [Green Version]
  31. Hibbeler, R.C. Mechanics of Materials, 9th ed.; PEARSON: London, UK, 2013. [Google Scholar]
  32. Ryu, K.W.; Seo, Y.S.; Jeon, H.M.; Cho, H.S.; Park, S.Y.; Jang, K.H.; Choi, J.H.; Song, P.M.; Jung, I.H.; Cheon, S.W.; et al. Conceptual Design for 20 MW Class Offshore Wind Power System around Korean Peninsula; R21XO01-6; Korea Electric Power Corporation: Naju-si, Republic of Korea, 2022. [Google Scholar]
  33. Ryu, K.W.; Jang, K.H. Aerodynamic Design and Loads Analysis for 20 MW Offshore Wind Turbine System. In Proceedings of the 2022 Spring Conference of Korea Wind Energy Association, Jeju, Republic of Korea, 20–22 June 2022. [Google Scholar]
  34. Desmond, C.; Murphy, J.; Blonk, L.; Haans, W. Description of an 8 MW Reference Wind Turbine. J. Phys. Conf. Ser. 2016, 753, 92013. [Google Scholar] [CrossRef]
  35. Schepers, J.G.; Ceyhan, O.; Savenije, F.J.; Stettner, M.; Kooijman, H.J.; Chaviarapoulos, P.; Sieros, G.; Simao Ferreira, C.; Sørensen, N.; Wächter, M.; et al. AVATAR: AdVanced Aerodynamic Tools for Large Rotors. In Proceedings of the 33rd Wind Energy Symposium, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar] [CrossRef] [Green Version]
  36. Oh, K.Y.; Kim, J.Y.; Lee, J.S.; Ryu, K.W. Wind resource assessment around Korean Peninsula for feasibility study on 100 MW offshore wind farm. Renew. Energy 2012, 42, 217–226. [Google Scholar] [CrossRef]
  37. Selig, M.S. PROPID User Manual; University of Illinois at Urbana-Champaign: Champaign, IL, USA, 2012. [Google Scholar]
  38. Jang, K.H.; Ryu, K.W. Blade Shape Design with Minimization of LCoE. In Proceedings of the 2022 Spring Conference of Korea Wind Energy Association, Jeju, Republic of Korea, 20–22 June 2022. [Google Scholar]
  39. Bortolotti, P.; Tarrés, H.C.; Dykes, K.; Merz, K.; Sethuraman, L.; Verelst, D.; Zahle, F. Systems Engineering in Wind Energy WP2.1 Reference Wind Turbines; Technical Report; IEA Wind; National Renewable Energy Laboratory: Oak Ridge, TN, USA, 2019. [Google Scholar]
  40. Griffith, D.T.; Ashwill, T.D. The Sandia 100-Meter All-Glass Baseline Wind Turbine Blade: SNL 100-00; SANSDIA REPORT; Sandia National Laboratories: Albuquerque, NM, USA, 2011. [Google Scholar]
  41. Chaviaropoulos, P.; Sieros, G. Design of low induction rotors for use in large offshore wind farms. In Proceedings of the Scientific Track, European Wind Energy Association (EWEA), Barcelona, Spain, 10–13 March 2014. [Google Scholar]
  42. Bergami, L.; Madsen, H.A.; Rasmussen, F. A two-bladed teetering hub configuration for the DTU 10 MW RWT: Loads considerations. In Proceedings of the EWEA 2014, Barcelona, Spain, 10–13 March 2014. [Google Scholar]
Figure 1. Control volume around a horizontal axis wind turbine.
Figure 1. Control volume around a horizontal axis wind turbine.
Energies 16 05169 g001
Figure 2. C T and C P vs. axial induction factor a.
Figure 2. C T and C P vs. axial induction factor a.
Energies 16 05169 g002
Figure 3. Aerodynamic loads in a blade section.
Figure 3. Aerodynamic loads in a blade section.
Energies 16 05169 g003
Figure 4. A blade element for the ith node.
Figure 4. A blade element for the ith node.
Energies 16 05169 g004
Figure 5. Structural modeling of a blade section for LCoE comparison.
Figure 5. Structural modeling of a blade section for LCoE comparison.
Energies 16 05169 g005
Figure 6. Load distribution along the blade span for LCoE comparison.
Figure 6. Load distribution along the blade span for LCoE comparison.
Energies 16 05169 g006
Figure 7. Cp   λ curve of a test rotor compared using PROPID and FAST.
Figure 7. Cp   λ curve of a test rotor compared using PROPID and FAST.
Energies 16 05169 g007
Figure 8. Comparison of spanwise mass distribution for the IEA 3.4 MW [39] wind turbine.
Figure 8. Comparison of spanwise mass distribution for the IEA 3.4 MW [39] wind turbine.
Energies 16 05169 g008
Figure 9. Comparison of spanwise mass distribution for the DTU 10 MW [5] wind turbine.
Figure 9. Comparison of spanwise mass distribution for the DTU 10 MW [5] wind turbine.
Energies 16 05169 g009
Figure 10. Comparison of spanwise normalized mass distribution for three selected reference wind turbines.
Figure 10. Comparison of spanwise normalized mass distribution for three selected reference wind turbines.
Energies 16 05169 g010
Figure 11. Normalized LCoE for various β values.
Figure 11. Normalized LCoE for various β values.
Energies 16 05169 g011
Figure 12. Mass ratio of Case 1 and Case 2 according to axial induction factor.
Figure 12. Mass ratio of Case 1 and Case 2 according to axial induction factor.
Energies 16 05169 g012
Figure 13. Cp λ curves for various axial induction factors.
Figure 13. Cp λ curves for various axial induction factors.
Energies 16 05169 g013
Figure 14. Rotor radius to compensate for Cp,max reduction.
Figure 14. Rotor radius to compensate for Cp,max reduction.
Energies 16 05169 g014
Figure 15. P-V curve for Case 1, Case 2, and Case 3 rotors.
Figure 15. P-V curve for Case 1, Case 2, and Case 3 rotors.
Energies 16 05169 g015
Figure 16. Chord length along the blade span for Case 1, Case 2, and Case 3 rotors.
Figure 16. Chord length along the blade span for Case 1, Case 2, and Case 3 rotors.
Energies 16 05169 g016
Figure 17. Twist angle along the blade span for Case 1, Case 2, and Case 3 rotors.
Figure 17. Twist angle along the blade span for Case 1, Case 2, and Case 3 rotors.
Energies 16 05169 g017
Figure 18. Load-V curve for Case 1, Case 2, and Case 3 rotors.
Figure 18. Load-V curve for Case 1, Case 2, and Case 3 rotors.
Energies 16 05169 g018
Figure 19. Mass per unit length along the blade span for Case 1, Case 2, and Case 3 rotors using Equation (22).
Figure 19. Mass per unit length along the blade span for Case 1, Case 2, and Case 3 rotors using Equation (22).
Energies 16 05169 g019
Figure 20. Normalized LCoE ratio of the two-bladed and three-bladed rotor.
Figure 20. Normalized LCoE ratio of the two-bladed and three-bladed rotor.
Energies 16 05169 g020
Figure 21. Cp λ curves for Case 3 and a two-bladed rotor.
Figure 21. Cp λ curves for Case 3 and a two-bladed rotor.
Energies 16 05169 g021
Table 1. Cost break-down for wind turbine blades.
Table 1. Cost break-down for wind turbine blades.
Turbine/ProjectRated Power [MW]Material [%]Labor [%]Etc. [%]
WindPACT [24]1.5472627
IEA 3.4 MW [24]3.4482527
SNL 100-00 [20]13.2721414
SNL 100-01 [20]13.2751213
SNL 100-03 [24]13.2601624
Table 2. Shares of the overall costs of various wind turbines.
Table 2. Shares of the overall costs of various wind turbines.
Turbine/ProjectRated Power [MW]Blades [%]ICC [%]O&M [%]Etc [%]
Baseline of [1]33.672.920.66.5
Ashur et al. [11]2028.487.93.98.4
Ashuri [21]105.870.026.13.9
Ashuri [21]208.273.523.23.3
INNWIND.EU [28]103.7362.537.5-
INNWIND.EU [28]103.7861.638.4-
Table 3. Airfoil distribution along each blade span.
Table 3. Airfoil distribution along each blade span.
r/RAirfoilThickness Ratio Perimeter   Ratio   of   Airfoil   ( P e / C ) Perimeter   Ratio   of   Ellipse   ( P e / C )   Relative Error [%]
0.046Cylinder1.006.286.280
0.089Cylinder1.006.286.280
0.132Cylinder1.006.286.280
0.187DU99_W_405LM0.402.282.394.8
0.252DU99_3500.352.222.355.9
0.317DU99_3500.352.222.355.9
0.382DU97_W_3000.302.162.327.4
0.447DU97_W_3000.302.162.327.4
0.512DU97_W_3000.302.162.327.4
0.577DU91_W2_2500.252.122.298.0
0.642DU91_W2_2500.252.122.298.0
0.707DU91_W2_2500.252.122.298.0
0.772DU93_W_210LM0.212.092.278.6
0.837DU93_W_210LM0.212.092.278.6
0.892DU93_W_210LM0.212.092.278.6
0.935DU93_W_210LM0.212.092.278.6
0.978DU93_W_210LM0.212.092.278.6
Table 4. Powers for reference turbines and selected commercial turbines.
Table 4. Powers for reference turbines and selected commercial turbines.
Turbine/ProjectRated Power [MW]Rated Wind Speed [m/s]Rotor Diameter [m]Specific Power [W/m2]
NREL 5 MW [2]511.4126401
DTU 10 MW [5]1011.4178.3401
IEA 15 MW [4]1510.59240332
ECN 20 MW [10]2010.0252401
INNWIND.EU 20 MW [17]2011.4252401
LEANWIND [34]812.5164379
AVATAR [35]1010.4205.8300
GE 12 MW12 220316
Unison 10 MW1012.0209292
Doosan 8 MW810.5205242
Table 5. Basic design parameters for low-induction rotors.
Table 5. Basic design parameters for low-induction rotors.
DescriptionsValues
Rated power [MW]20
Number of blades3
Blade length [m]150.81~162.55
Rated wind speed [m/s]10
Tip speed ratio [-]9
Rated rotational speed [rpm]5.54~5.7
Cp,max [-]0.421~0.489
Axial induction factor a 0.2~0.3
Table 6. Comparison of the mass predicted by the mass model for the reference wind turbines.
Table 6. Comparison of the mass predicted by the mass model for the reference wind turbines.
TurbineMDAORated Power [MW]Blade Radius
[m]
Blade Mass
[kg]
Predicted Mass [kg]
NREL 5 MW [2]X56317,74017,740
IEA 15 MW [4]O1512065,250189,169
DTU 10 MW [5]O1089.1541,72247,172
ECN 20 MW [10]X20123161,000135,383
Ashuri 20 MW [11]O20135259,000215,072
SNL 100-03 [13]X13.210049,519101,232
IEA 3.4 MW [39]O3.46516,44122,147
SNL 100-00 [40]X13.2102.5114,17280,581
Table 7. Specifications of designed 20 MW rotors.
Table 7. Specifications of designed 20 MW rotors.
DescriptionsCase 1 RotorCase 2 RotorCase 3 Rotor
Number of blades3
Blade length [m]150.81162.55154.02
Blade length/Blade length_Case 1 [%]100107.8102.1
Cut-in wind speed [m/s]3
Rated wind speed [m/s]10
Cut-out wind speed [m/s]25
Tip speed ratio [-]9
Rated rotational speed [rpm]5.75.295.58
Cp,max [-]0.48880.42070.4686
Cp,max/Cp,max_Case 1 [%]10086.195.9
Designed axial induction factor a [-]0.30.20.26
Blade mass/Blade mass_Case 1 [%]100102.292.0
Table 8. Comparison of LCoE calculation results for a 20 MW wind turbine [11].
Table 8. Comparison of LCoE calculation results for a 20 MW wind turbine [11].
Ashuri 20 MW [11]Case 1Case 2Case 3
EquipmentCost [k$]Mass [ton]Cost [k$]Mass [ton]Cost [k$]Mass [ton]Cost [k$]Mass [ton]
Blade4051.72593735.8228.43831.4234.23430209.7
Hub1456.9252.81288.4223.61320.5229.11185.8205.8
Pitch system1945.32362468.6208.33005.6213.52604.3191.4
Hub cone34.64.638.25.141.55.539.15.2
Main shaft1605.3159.12073.6205.62575255.32203.7218.5
Shaft bearing1013.442.51380.257.81795.475.2148662.3
Gearbox4955.5161.94981.9 4981.9 4981.9
Drive train brake44.4444.6444.6444.64
Generator1592.259.8159960159960159960
Electronics1572.8 1572.8 1572.8 1572.8
Yaw system1495176.81945237.32429.1304.22070.3254.5
Nacelle frame752.6280.8895.1348.61036.3403.6932.7363.2
Nacelle railing414.235.1514.143.6595.150.4535.645.4
Nacelle cover279.623.4291.823.5291.823.5291.823.5
Turbine connection1235.5 1247.2 1247.2 1247.2
Cooling and hydraulic system3091.6310.31.6310.31.6310.31.6
Monitoring and safety system65.4 65.4 65.4 65.4
Tower39711588.33971.01588.33971.01588.33971.01588.3
Turbine Capital Costs34,897.8 35,888.5 38,376.6
Foundation290.7 290.7 290.7 290.7
Installation363.1 363.1 363.1 363.1
Farm connection838.2 838.2 838.2 838.2
Site assessment and permit934.5 934.5 934.5 934.5
Balance of Station2426.5 2426.5 2426.5 2426.5
Initial Capital Cost37,324.3 38,315 40,803.1 37,858
Levelized Replacement Cost249.3 249.3 249.3 249.3
Maintenance and
Operation
108.7 118.1 118.1 118.1
Interest Rate (FCR)0.07 0.07 0.07 0.07
AEP [GWh]86 93.4 93.4 93.4
LCoE [$/kWh]0.03454(0.03051)0.03263(0.02890)0.03450(0.03077)0.03229(0.02856)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jang, K.-H.; Ryu, K.-W. Blade Design and Aerodynamic Performance Analysis of a 20 MW Wind Turbine for LCoE Reduction. Energies 2023, 16, 5169. https://doi.org/10.3390/en16135169

AMA Style

Jang K-H, Ryu K-W. Blade Design and Aerodynamic Performance Analysis of a 20 MW Wind Turbine for LCoE Reduction. Energies. 2023; 16(13):5169. https://doi.org/10.3390/en16135169

Chicago/Turabian Style

Jang, Kang-Ho, and Ki-Wahn Ryu. 2023. "Blade Design and Aerodynamic Performance Analysis of a 20 MW Wind Turbine for LCoE Reduction" Energies 16, no. 13: 5169. https://doi.org/10.3390/en16135169

APA Style

Jang, K. -H., & Ryu, K. -W. (2023). Blade Design and Aerodynamic Performance Analysis of a 20 MW Wind Turbine for LCoE Reduction. Energies, 16(13), 5169. https://doi.org/10.3390/en16135169

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop