Research article 03 May 2021
Research article  03 May 2021
On the scaling of wind turbine rotors
 ^{1}Wind Energy Institute, Technische Universität München, 85748 Garching b. München, Germany
 ^{2}National Renewable Energy Laboratory, Golden, CO 80401, USA
 ^{1}Wind Energy Institute, Technische Universität München, 85748 Garching b. München, Germany
 ^{2}National Renewable Energy Laboratory, Golden, CO 80401, USA
Correspondence: Carlo L. Bottasso (carlo.bottasso@tum.de)
Hide author detailsCorrespondence: Carlo L. Bottasso (carlo.bottasso@tum.de)
This paper formulates laws for scaling wind turbine rotors. Although the analysis is general, the article primarily focuses on the subscaling problem, i.e., on the design of a smallersized model that mimics a fullscale machine. The present study considers both the steadystate and transient response cases, including the effects of aerodynamic, elastic, inertial, and gravitational forces. The analysis reveals the changes to physical characteristics induced by a generic change of scale, indicates which characteristics can be matched faithfully by a subscaled model, and states the conditions that must be fulfilled for desired matchings to hold.
Based on the scaling laws formulated here, the article continues by considering the problem of designing scaled rotors that match desired indicators of a fullscale reference. To better illustrate the challenges implicit in scaling and the necessary tradeoffs and approximations, two different approaches are contrasted. The first consists in a straightforward geometric zooming. An analysis of the consequences of zooming reveals that, although apparently simple, this method is often not applicable in practice, because of physical and manufacturing limitations. This motivates the formulation of scaling as a constrained optimal aerodynamic and structural matching problem of wide applicability.
Practical illustrations are given considering the scaling of a large reference 10 MW wind turbine of about 180 m in diameter down to three different sizes of 54, 27, and 2.8 m. Results indicate that, with the proper choices, even models characterized by very significant scaling factors can accurately match several key performance indicators. Additionally, when an exact match is not possible, relevant trends can at least be captured.
This article is concerned with the aeroservoelastic scaling of wind turbine rotors. The general scaling problem includes both up and subscaling (or downscaling). This work primarily focuses on the latter aspect – i.e., on the design of subscaled models – but also briefly touches upon the former. Specifically, this work tries to answer the following scientific questions:

What are the effects of a change of scale (i.e., both in the case of up and subscaling) on the steady and transient response of a wind turbine?

What steady and transient characteristics of the response of a fullscale wind turbine can be matched by a subscaled model?

What are the most suitable ways to design the aerodynamic and structural characteristics of a subscaled model?
The understanding of both up and subscaling is relevant to contemporary wind energy technology.
Regarding upscaling, wind turbines have experienced a continuous growth in size in the past decades. This trend has been mostly driven by a desire for increased capacity factors, which can be obtained essentially through two main design parameters: by lowering the specific power – which, for a given power rating, means a larger rotorswept area – and by designing taller towers, which reach higher above ground, where wind blows faster. In turn, improved capacity factors have a positive effect on the cost of energy, which has helped propel the penetration of wind in the energy mix. The design of the nextgeneration wind turbines, especially for offshore applications, is expected to follow this same path, with rotor diameters of present and future products already exceeding 200 m (IRENA, 2019; GE, 2019; Siemens Gamesa, 2020). Unfortunately, larger blades cannot be obtained by simply scaling up existing smaller blades but must be designed to beat the cubic law of growth. In fact, weight (and hence cost) grows with volume – i.e., with the cube of size – whereas power capture only grows with rotor swept area, i.e., with the square of size (Sieros et al., 2012). Against this background, it is clearly useful to understand the changes that can be expected in a turbine response as the result of an increase in size.
Subscaling, on the other hand, is useful as a research tool: by designing and testing smallerscale versions of fullscale references, one can validate simulation tools, explore ideas, compare alternative solutions, and deepen the knowledge and understanding of complex physical phenomena. Among other advantages, scaled testing is usually much cheaper and less risky than fullscale testing. In addition, fullscale testing is typically performed on prototypes or even commercial products, which raises often unsurmountable issues because of intellectual property rights and trade secrecy. In turn, this limits opportunities for publication, data sharing, and full exploitation of the results from the scientific community. With commercial turbine sizes expected to grow even further in the future, it is becoming more important than ever to fully understand how to best employ subscaling as a research tool.
Two subscaled testing activities are possible: wind tunnel testing with smallscale models and field testing with small turbines. In both cases, the goal is to match at least some of the characteristics of the original fullscale problem. Clearly, this requires a complete understanding of the effects of a change (in this case, a reduction) of scale on the response of a wind turbine.
Wind tunnel testing of subscaled wind turbine models offers some unique opportunities. First, the operating conditions in a wind tunnel are to a large extent controllable and typically highly repeatable. Second, measurements – especially of flow quantities – that are possible in the lab environment are generally more difficult, are less precise, and have a lower resolution at full scale. Third, costs and risks are much more limited than in the case of field testing, and the time for the conduction of the experiments is shorter (not only because of the reduced challenges but also because of time acceleration, as explained later). Fourth, since a smallscale model cannot exactly match a fullscale product, property right issues are typically much less of a constraint.
The first wind tunnel experiments on wind turbine aerodynamics were conducted in the last decades of the 20th century, as summarized in Vermeer et al. (2003). Studies carried out during the Unsteady Aerodynamics Experiment (Simms et al., 2001) with a stallregulated 10 m diameter, 20 kW turbine were, among others, key to uncovering the importance of specific flow phenomena, such as dynamic stall, 3D rotational effects, and tower–wake interactions. Later, the 4.5 m diameter scaled models designed for the Model rotor EXperiments In controlled COnditions (MEXICO) project enabled the validation of multiple aerodynamic models, ranging from blade element momentum (BEM) to computational fluid dynamics (CFD) (Snel et al., 2009). These wind turbine models were designed following a set of scaling laws aimed at replicating as accurately as possible the aerodynamic behavior of fullscale machines. More recently, the inclusion of closedloop controls and aeroservoelastic considerations in the scaling process expanded the scope of wind tunnel testing beyond aerodynamics (Campagnolo et al., 2014). Nowadays, wind tunnel tests are extensively used to gain a better understanding of wake effects, to validate simulation tools, and to help develop novel control strategies (Bottasso and Campagnolo, 2020). The recent study of Wang et al. (2020) tries to quantify the level of realism of wakes generated by smallscale models tested in a boundary layer wind tunnel.
Unfortunately, the exact matching of all relevant physical processes between fullscale and subscale models is typically not possible. This mismatch increases with the scale ratio, and it becomes especially problematic when large wind turbines (with rotor sizes on the order of 10^{2} m and power ratings on the order of 10^{6}–10^{7} W) are scaled to very small size wind tunnel models (characterized by rotors on the order of 10^{−1}–10^{0} m and power ratings on the order of 10^{0}–10^{2} W). To limit the scale factor, instead of using very small models in a wind tunnel, testing can be conducted in the field with smallsize wind turbines (with a rotor on the order of 10^{1} m and power ratings on the order of 10^{5} W).
Examples of stateoftheart experimental test sites realized with smallsize wind turbines are the Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas (Berg et al., 2014), which uses three 27 m diameter Vestas V27 225 kW turbines, or the soontobeready WINSENT complexterrain facility in the German Swabian Alps (ZSW, 2016), which uses two 54 m diameter S&G 750 kW turbines.
Reducing the scaling ratios and moving to the field offers the opportunity to overcome some of the constraints typically present in wind tunnel testing, although some of the advantages of wind tunnels are clearly lost. Indeed, the range of testing conditions cannot be controlled at will, measurements are more difficult, and costs are higher. Here research has so far mainly focused on steadystate aerodynamics and wake metrics. For example, within the National Rotor Testbed project (Resor and Maniaci, 2014), teams at the University of Virginia, Sandia National Laboratories, and National Renewable Energy Laboratory have designed a blade for the SWiFT experimental facility, replacing the original Vestas V27 blade. The scaling laws were specifically chosen to replicate the wake of a commercial 1.5 MW rotor at the subscale size of the V27 turbine. To capture the dynamic behavior of very large wind turbines, additional effects must, however, be considered in the scaling laws. For example, Loth et al. (2017) have recently proposed a methodology to include gravity in the scaling process, and they have demonstrated their approach to scale a 100 m blade down to a 25 m size. Gravity is also crucially important in floating offshore applications (Azcona et al., 2016) to balance buoyancy and correctly represent flotation dynamics, with its effects on loads, stability, and performance and with implications for control design.
This paper considers the general problem of scaling a wind turbine rotor to a different size, including the effects caused by aerodynamic, elastic, inertial, and gravitational forces. The study is structured in two main parts.
Initially, an analysis of the problem of scaling is presented. The main steady and transient characteristics of a rotor in terms of performance, aeroservoelasticity, and wake shedding are considered, and the effects caused by a generic change of scale are determined. The analysis reveals that, in principle, most of the response features can be faithfully represented by a subscaled model. However, an exact matching of all features is typically impossible because of chordbased Reynolds effects, which lead to changes in the aerodynamic behavior of the system. Another limit comes from wind conditions: the wind field is not scaled when using utilitysize models in the field, and wind tunnel flows can only partially match the characteristics of the atmospheric boundary layer. The analysis also shows that scaling is essentially governed by two parameters: the geometric (length) scaling factor and the time scaling factor. Based on these two parameters, all matched and unmatched quantities can be fully characterized.
In the second part, the paper continues by looking at the problem of designing a subscaled model. Two different approaches are considered. The first is a straightforward zooming down of all blade characteristics based on a pure geometrical scaling (Loth et al., 2017), which is appealing for its apparent simplicity. The second is based on a complete aerostructural redesign, which is formulated here in terms of two constrained optimizations: the aerodynamic one defines the external shape of the blade, whereas the structural optimization sizes the structural components. Both strategies aim at replicating the dynamic behavior (including gravitational effects) of a fullscale wind turbine at a smaller scale, and they are therefore based on the same scaling laws. Clearly, the complete redesign is a more complicated process than the pure geometric zoomingdown approach. However, the main goal of scaling is that of designing a rotor that matches at scale the behavior of a target fullscale machine as closely as possible. From this point of view, the simplicity of design – which is a oneoff activity – is less of a concern, especially today, when sophisticated automated rotor design tools are available (Bortolotti et al., 2016). Apart from simplicity, zooming is very often simply not possible for large scale factors because of unrealistically small sizes (especially the thickness of shell structures), nonachievable material characteristics, or impossibletoduplicate manufacturing processes (Wan and Cesnik, 2014; Ricciardi et al., 2016). In all those cases, a different aerodynamic shape, a different structural configuration, and different materials are used to obtain the desired behavior, as shown, for example, in the design of a smallsize aeroelastically scaled rotor by Bottasso et al. (2014), or as customarily done in the design of scaled flutter models for aeronautical applications (Busan, 1998).
Although the intrinsic limits of the straightforward zoomingdown approach are probably well understood, these two alternative methodologies are compared here in order to give a better appreciation of the complexities that one has to face in the design of scaled models. To give practical and concrete examples, a very large rotor is scaled down to three different model sizes, including two different utility wind turbines and a smallscale wind tunnel model. For each model, the zoomingdown approach is adopted when possible for its simplicity and then replaced by the redesign method when fidelity or physical limits make it impractical or impossible.
Furthermore, the paper analyzes the accuracy with which the subscale models successfully mirror relevant key characteristics of the fullscale reference, in terms both of absolute values and of trends. This is indeed an important aspect of scaling: even if the exact matching of certain quantities is sometimes not possible, scaled models can still be highly valuable if they are able to at least capture trends. As an example of such a trend analysis, the subscale models are used here to explore changes in loading between unwaked and waked inflow conditions, which are then validated against the corresponding loading changes of the fullscale machine. Results indicate that even the smallest model is capable of capturing complex details of wake interaction, including an interesting lack of symmetry for left/right wake impingements caused by rotor uptilt.
A final section completes the paper, listing the main conclusions that can be drawn from the results and highlighting their limits.
Buckingham's Π theorem (Buckingham, 1914) states that a scaled model (labeled (⋅)_{M}) has the same behavior as a fullscale physical system (labeled (⋅)_{P}) if all the m relevant nondimensional variables π_{i} are matched between the two systems. In other words, when the governing equations are written as
the two systems are similar if
Depending on the scaled testing conditions, not all dimensional quantities can usually be matched. In the present case, we consider testing performed in air, either in a wind tunnel or in the field, neglecting hydrodynamics.
The length (geometric) scale factor between scaled and fullscale systems is defined as
where l is a characteristic length (for example the rotor radius R), whereas the scale factor for time t is defined as
As a consequence of these two definitions, one can determine the angular velocity and wind speed scaling factors, which are respectively written as ${n}_{\mathrm{\Omega}}={\mathrm{\Omega}}_{\mathrm{M}}/{\mathrm{\Omega}}_{\mathrm{P}}=\mathrm{1}/{n}_{\mathrm{t}}$ and ${n}_{\mathrm{v}}={V}_{\mathrm{M}}/{V}_{\mathrm{P}}={n}_{\mathrm{l}}/{n}_{\mathrm{t}}$. A nondimensional time can be defined as τ=tΩ_{r}, where Ω_{r} is a reference rotor speed, for example the rated one. It is readily verified that, by the previous expressions, nondimensional time is matched between the model and physical system; i.e., τ_{M}=τ_{P}. The two factors n_{l} and n_{t} condition, to a large extent, the characteristics of a scaled model.
The following Sects. 2.1 and 2.2 analyze the main steady and transient characteristics of a rotor in terms of performance, aeroservoelasticity, and wake shedding. The analysis discusses which of these characteristics can be matched by a scaled model and which conditions are required for the matchings to hold. Next, Sect. 2.3 offers an overview on the main scaling relationships and discusses the choice of scaling parameters.
2.1 Steady state
2.1.1 Rotor aerodynamics
The power coefficient characterizes the steadystate performance of a rotor, and it is defined as ${C}_{\mathrm{P}}=P/(\mathrm{1}/\mathrm{2}\mathit{\rho}A{V}^{\mathrm{3}})$, where P is the aerodynamic power, ρ the density of air, A=πR^{2} the rotor disk area, and V the ambient wind speed. The thrust coefficient characterizes the wake deficit and the rotor loading and is defined as ${C}_{\mathrm{T}}=T/(\mathrm{1}/\mathrm{2}\mathit{\rho}A{V}^{\mathrm{2}})$, where T is the thrust force. For a given rotor, the power and thrust coefficients depend on tipspeed ratio (TSR) $\mathit{\lambda}=\mathrm{\Omega}R/V$ and blade pitch β, i.e., ${C}_{\mathrm{P}}={C}_{\mathrm{P}}(\mathit{\lambda},\mathit{\beta})$ and ${C}_{\mathrm{T}}={C}_{\mathrm{T}}(\mathit{\lambda},\mathit{\beta})$.
It is readily verified that λ_{M}=λ_{P} for any n_{l} and n_{t}, which means that it is always possible to match the scaled and fullscale TSR. This ensures the same velocity triangle at the blade sections and the same wake helix pitch.
Ideally, a scaled model should match the C_{P} and C_{T} coefficients of a given fullscale target; it is clearly desirable for the match not to hold at a single operating point but over a range of conditions. BEM theory (Manwell et al., 2002) shows that both rotor coefficients depend on the steadystate aerodynamic characteristics of the airfoils. In turn, the lift C_{L} and drag C_{D} coefficients of the aerodynamic profiles depend on the angle of attack and on the Mach and Reynolds numbers.
The local Mach number accounts for compressibility effects and is defined as $\mathit{\text{Ma}}=W/{a}_{\mathrm{s}}$, where W is the flow speed relative to a blade section, and a_{s} is the speed of sound. Using the previous expressions, the Mach number of the scaled model is ${\mathit{\text{Ma}}}_{\mathrm{M}}={\mathit{\text{Ma}}}_{\mathrm{P}}{n}_{\mathrm{l}}/{n}_{\mathrm{t}}^{\mathrm{2}}$. Because of typical tip speeds, compressibility does not play a significant role in wind turbines. Hence, the matching of the Mach number can be usually neglected for current wind turbines. The situation might change for future offshore applications, where, without the constraints imposed by noise emissions, higher tipspeed and TSR rotors may have interesting advantages.
The Reynolds number represents the ratio of inertial to viscous forces and is defined as $\mathit{\text{Re}}=\mathit{\rho}lu/\mathit{\mu}$, where l is a characteristic length, u a characteristic speed, and μ the dynamic viscosity. In the present context, the most relevant definition of the Reynolds number is the one based on the blade sections, where l=c is the chord length, and u=W is the flow speed relative to the blade section. In fact, the Reynolds number has a strong effect on the characteristics and behavior of the boundary layer that develops over the blade surface, which in turn, through the airfoil polars, affects the performance and loading of the rotor. Testing in air in a wind tunnel or in the field (hence with similar ρ and μ but with a reduced chord c) leads to a mismatch between the scaled and fullscale chordbased Reynolds numbers, as ${\mathit{\text{Re}}}_{\mathrm{M}}={\mathit{\text{Re}}}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{2}}/{n}_{\mathrm{t}}$.
The effects due to a chordbased Reynolds mismatch can be mitigated by replacing the airfoils of the fullscale system with others better suited for the typical Reynolds conditions of the scaled model (Bottasso et al., 2014). A second approach is to increase the chord of the scaled model. This, however, has the effect of increasing the rotor solidity – defined as $\mathrm{\Sigma}=B{A}_{\mathrm{b}}/A$, where B is the number of blades and A_{b} the blade planform area – which may have additional consequences. In fact, the TSR of the maximum power coefficient is directly related to rotor solidity. This can be shown by using classical BEM theory with wake swirl, which gives the optimal blade design conditions by maximizing power at a given design TSR λ_{d}. By neglecting drag, the optimal design problem can be solved analytically to give the chord distribution of the optimal blade along the spanwise coordinate r (Manwell et al., 2002):
Although based on a simplified model that neglects some effects, this expression shows that chord distribution and design TSR are linked. This means that, if one increases solidity (and hence chord) to contrast the Reynolds mismatch while keeping C_{L} fixed, the resulting rotor will have a lower TSR for the maximum power coefficient. Therefore, this technique of correcting the Reynolds number moves the optimal TSR away from the one of the fullscale reference, which may or may not be acceptable, depending on the goals of the model. For example, if one wants to match the behavior of the C_{P}–λ curves over a range of TSRs, such an approach would not be suitable. As shown by Eq. (5), this effect can be eliminated or mitigated by changing the design C_{L} accordingly; however, if this moves the operating condition of the airfoil away from its point of maximum efficiency, a lower maximum power coefficient will be obtained.
In addition, chord c and lift C_{L} are further constrained by the circulation $\mathrm{\Gamma}=\mathrm{1}/\mathrm{2}c{C}_{\mathrm{L}}W$ (Burton et al., 2001), which plays an important role in the aerodynamics of the rotor and its wake.
Considering first the rotor, the lift and drag generated by the airfoils located close to the blade root are modified by the combined effects of centrifugal and Coriolis forces. In fact, the former causes a radial pumping of the flow that, as a result, moves outboard in the spanwise direction. This radial motion over a rotating body generates chordwise Coriolis forces that alleviate the adverse pressure gradient on the airfoils and, in turn, delay stall. As shown by the dimensional analysis developed by Dowler and Schmitz (2015), rotational augmentation causes multiplicative corrections, denoted ${g}_{{C}_{\mathrm{L}}}$ and ${g}_{{C}_{\mathrm{D}}}$, to the nonrotating lift and drag coefficients, which can be written, respectively, as
where Δθ is the total blade twist from root to tip. Equations (6a) and (6b) show that, in order to match the effects of rotational augmentation, the model and fullscale system should have the same blade nondimensional chord and twist distributions; the same nondimensional circulation $\mathrm{\Gamma}/\left(RW\right)$; and the same Rossby number $\mathit{\text{Ro}}=\mathrm{\Omega}r/\left(\mathrm{2}\phantom{\rule{0.125em}{0ex}}W\right)$, which represents the ratio of inertia to Coriolis forces. Matching nondimensional circulation between the two systems implies matching either both the planform shape $c/R$ and the lift coefficient C_{L} or the product of the two. As previously noted, some of these options may lead to a different TSR of optimal C_{P}. On the other hand, it is readily verified that the Rossby number is always matched for any choice of n_{l} and n_{t}.
2.1.2 Wake aerodynamics
The circulation is relevant not only for rotational augmentation but also for wake behavior. In fact, each blade sheds trailing vorticity that is proportional to the spanwise gradient $\mathrm{d}\mathrm{\Gamma}/\mathrm{d}r$ (Schmitz, 2020). Therefore, designing a blade that matches the spanwise distribution of Γ (and, hence, also its spanwise gradient) ensures that the scaled rotor sheds the same trailed vorticity. Additionally, a matched circulation ensures also a matched thrust, which is largely responsible for the speed deficit in the wake and for its deflection in misaligned conditions (Jiménez et al., 2010).
The Reynolds mismatch derived earlier applies also to its rotorbased definition, which is relevant to wake behavior and is obtained by using l=2R and u=V. However, Chamorro et al. (2012) showed that the wake is largely unaffected by this parameter as long as Re>10^{5}, which is typically the case unless extremely small model turbines are used. The same is true for the terrainheightbased Reynolds number definition that applies to flows over complex terrains, where Reynoldsindependent results are obtained when Re>10^{4} (McAuliffe and Larose, 2012).
The detailed characterization of the behavior of scaled wakes is considered as out of the scope of the present investigation, and the interested reader is referred to Wang et al. (2020) for a specific study on this important topic.
2.1.3 Gravity
The Froude number represents the ratio of aerodynamic to gravitational forces and is written as $\mathit{\text{Fr}}={V}^{\mathrm{2}}/gR$, where g is the acceleration of gravity. The Froude number of the scaled model is readily found to be ${\mathit{\text{Fr}}}_{\mathrm{M}}={\mathit{\text{Fr}}}_{\mathrm{P}}{n}_{\mathrm{l}}/{n}_{\mathrm{t}}^{\mathrm{2}}$. Enforcing Froude (Fr_{M}=Fr_{P}) results in the time scaling factor being set to ${n}_{\mathrm{t}}=\sqrt{{n}_{\mathrm{l}}}$. This condition determines the only remaining unknown in the scaling laws, so that the scalings of all nondimensional parameters can now be expressed in terms of the sole geometric scaling factor n_{l}. Froude scaling is used when gravity plays an important role, for example in the loading of very large rotors or for floating offshore applications where weight and buoyancy forces should be in equilibrium.
2.1.4 Elasticity
The steady deflections due to aerodynamic loading of the scaled and fullscale wind turbines can be matched by adjusting the stiffness of the scaled model. In fact, consider the very simplified model of a blade represented by a clamped beam of length R under a uniformly distributed aerodynamic load per unit span, denoted $q=\mathrm{1}/\mathrm{2}\mathit{\rho}{W}^{\mathrm{2}}c{C}_{\mathrm{L}}$. The beam nondimensional tip deflection is $s/R=q{R}^{\mathrm{3}}/\left(\mathrm{8}EJ\right)$, where EJ is the bending stiffness, E is Young's modulus, and J is the crosssectional moment of inertia. By the previous definitions of length scale and timescale, one gets $(s/R{)}_{\mathrm{M}}=(s/R{)}_{\mathrm{P}}$ if $(EJ{)}_{\mathrm{M}}=(EJ{)}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{6}}/{n}_{\mathrm{t}}^{\mathrm{2}}$. Hence, nondimensional deflections can be matched, provided that the stiffness can be adjusted as shown. Matching this requirement may imply changing the material and/or the configuration of the structure, because of technological, manufacturing, and material property constraints (Busan, 1998; Ricciardi et al., 2016), as discussed more in detail later on.
2.2 Transient response
A scaled model should obey some additional conditions in order for the transient response of the fullscale system to be matched.
2.2.1 Rotor aerodynamics and inflow
As mentioned earlier, any aerodynamically scaled model can always be designed to enforce the TSR without additional conditions. To extend the similitude to dynamics, the nondimensional time derivative of the TSR should also be matched, i.e., ${\mathit{\lambda}}_{\mathrm{M}}^{\prime}={\mathit{\lambda}}_{\mathrm{P}}^{\prime}$, where a nondimensional time derivative is denoted as $(\cdot {)}^{\prime}=\mathrm{d}\cdot /\mathrm{d}\mathit{\tau}$. By using the definition of λ, one gets
The rotor dynamic torque balance equilibrium is written as $I\dot{\mathrm{\Omega}}=Q$. In this expression, I is the rotor polar moment of inertia, $\dot{(\cdot )}=\mathrm{d}\cdot /\mathrm{d}t$ indicates a derivative with respect to time, and $Q={Q}_{\mathrm{a}}({Q}_{\mathrm{e}}+{Q}_{\mathrm{m}})$ is the shaft torque. The aerodynamic torque is denoted as ${Q}_{\mathrm{a}}=\mathrm{1}/\mathrm{2}\mathit{\rho}AR{C}_{\mathrm{P}}/\mathit{\lambda}$, while Q_{e} is the electrical torque provided by the generator and Q_{m} the mechanical losses. The aerodynamic torque scales as ${Q}_{{a}_{\mathrm{M}}}={Q}_{{a}_{\mathrm{P}}}{n}_{\mathrm{l}}^{\mathrm{5}}/{n}_{\mathrm{t}}^{\mathrm{2}}$, and clearly Q_{e}+Q_{m} must scale accordingly. Since the mechanical losses depend on friction, it might be difficult to always match Q_{m}, especially in a smallscale model. This problem, however, can be eliminated by simply providing the necessary electrical torque to generate the correct term, Q_{e}+Q_{m}. Considering that the dimensions of I are [I]=[ρ_{m}][l]^{5}, where ρ_{m} is the material density and l a characteristic length, the first term ${\mathrm{\Omega}}^{\prime}R/V$ in Eq. (7) is matched between the two models if the material density is matched, i.e., if ${\mathit{\rho}}_{{m}_{\mathrm{M}}}={\mathit{\rho}}_{{m}_{\mathrm{P}}}$.
The second term, $\mathit{\lambda}{V}^{\prime}/V$, in Eq. (7) is matched if the two systems operate at the same TSR and if the wind speed has the same spectrum of the wind in the field. The matching of wind fluctuations (clearly, only in a statistical sense) induces not only the same variations in the TSR, and hence in the rotor response, but also the same recovery of the wake, which is primarily dictated by the ambient turbulence intensity (Vermeer et al., 2003).
Matching of the wind spectrum is in principle possible in a boundary layer wind tunnel if a flow of the desired characteristics can be generated. Turbulent flows can be obtained by active (Hideharu, 1991; Mydlarski, 2017) or passive means (Armitt and Counihan, 1968; Counihan, 1969). Active solutions are more complex and expensive but also more flexible and capable of generating a wider range of conditions. When testing in the field, the flow is invariably not scaled. This will have various effects on the scaled model response, which might be beneficial or not depending on the goals of scaled testing. In fact, the acceleration of time (t_{M}=t_{P}n_{t}) implies a shift in the wind frequency spectrum. Among other effects, this means that lowprobability (extreme) events happen more frequently than at full scale. Similarly, the scaling of speed (${V}_{\mathrm{M}}={V}_{\mathrm{P}}{n}_{\mathrm{l}}/{n}_{\mathrm{t}}$) implies higher amplitudes of turbulent fluctuations and gusts than at full scale.
Magnitude and phase of the aerodynamic response of an airfoil (as for example modeled by Theodorsen's theory (Bisplinghoff and Ashley, 2002)) are governed by the reduced frequency $\mathit{\kappa}={\mathit{\omega}}_{\mathrm{m}}c/\left(\mathrm{2}W\right)$, where ω_{m} is the circular frequency of motion. Harmonic changes in angle of attack take place at various frequencies ${\mathit{\omega}}_{{\mathrm{m}}_{j}}$ and are caused by the inhomogeneities of the flow (shears, misalignment between rotor axis and wind vector), blade pitching, and structural vibrations in bending and twisting. The reduced frequency can be written as ${\mathit{\kappa}}_{j}={\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{{\mathrm{m}}_{j}}\mathrm{\Omega}c/\left(\mathrm{2}W\right)$, where ${\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{{\mathrm{m}}_{j}}={\mathit{\omega}}_{{\mathrm{m}}_{j}}/\mathrm{\Omega}$ indicates a nondimensional frequency. This expressions shows that, once the nondimensional frequencies ${\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{{\mathrm{m}}_{j}}$ (due to inflow, pitch, and vibrations) are matched, the corresponding reduced frequencies are also matched, as the term $\mathrm{\Omega}c/\left(\mathrm{2}W\right)$ is always automatically preserved between scaled and fullscale systems for any n_{l} and n_{t}.
Dynamic stall effects depend on reduced frequency κ and chordbased Reynolds number. Typical dynamic stall models depend on the lift, drag, and moment static characteristics of an airfoil and various time constants that describe its unsteady inviscid and viscous response (Hansen et al., 2004). As previously argued, κ can be matched, and all time constants are also automatically matched by the matching of nondimensional time. However, a mismatch of the chordbased Reynolds number is typically unavoidable and will imply differences in the dynamic stall behavior of the scaled and fullscale models, which will have to be quantified on a casebycase basis.
2.2.2 Wake aerodynamics
The Strouhal number is associated with vortex shedding, which has relevance in tower and rotor wake behavior; the Strouhal number has also been recently used to describe the enhanced wake recovery obtained by dynamic induction control (Frederik et al., 2019). A rotor–wakerelevant definition of this nondimensional parameter is $\mathit{\text{St}}=f\mathrm{2}R/V$, where f is a characteristic frequency. Using the previous relationships, it is readily shown that ${\mathit{\text{St}}}_{\mathrm{M}}={\mathit{\text{St}}}_{\mathrm{P}}{n}_{\mathrm{l}}/\left({n}_{\mathrm{t}}{n}_{\mathrm{v}}\right)=\mathrm{1}$; i.e., the Strouhal number is always exactly matched between scaled and fullscale models for any n_{l} and n_{t} when TSR is matched.
During transients, spanwise vorticity is shed that is proportional to its temporal gradient. Using BEM theory (Manwell et al., 2002, p. 175), the nondimensional spanwise circulation distribution is computed as
In this expression, C_{L,α} is the slope of the lift curve, θ is the sectional pitch angle, and U_{P} and U_{T} are the flow velocity components at the blade section, respectively perpendicular and tangent to the rotor disk plane, such that ${W}^{\mathrm{2}}={U}_{\mathrm{P}}^{\mathrm{2}}+{U}_{\mathrm{T}}^{\mathrm{2}}$. The flow speed component tangential to the rotor disk is ${U}_{\mathrm{T}}=\mathrm{\Omega}r+{u}_{\mathrm{T}}$, where u_{T} contains terms due to wake swirl and yaw misalignment. The flow speed component perpendicular to the rotor disk is ${U}_{\mathrm{P}}=(\mathrm{1}a)V+\dot{d}+{u}_{\mathrm{P}}$, where a is the axial induction factor, $\dot{d}$ the outofplane blade section flapping speed, and u_{P} the contribution due to yaw misalignment and vertical shear. Neglecting u_{P} and u_{T} and using Eq. (8), the nondimensional time rate of change of the circulation becomes
For a correct similitude between scaled and fullscale systems, the nondimensional derivatives λ^{′}, a^{′}, θ^{′}, and $(\dot{d}/V{)}^{\prime}$ should be matched.
The matching of λ^{′} has already been addressed. The term a^{′} accounts for dynamic changes in the induction, which are due to the speed of actuation (of torque and blade pitch) and the intrinsic dynamics of the wake. The speed of actuation is matched if the actuators of the scaled model are capable of realizing the same rates of change of the fullscale system, i.e., if θ^{′} is matched. The intrinsic dynamics of the wake are typically modeled by a firstorder differential equation (Pitt and Peters, 1981):
where a represents inflow states and A is a matrix of coefficients proportional to $V/R$. It is readily verified that the matching of nondimensional time results in the matching of a^{′}. Finally, the term $(\dot{d}/V{)}^{\prime}$ is due to the elastic deformation of the blade, which is addressed next.
2.2.3 Elasticity
Considering blade flapping, the Lock number Lo is defined as
where I_{b} is the blade flapping inertia. Matching the Lock number ensures the same ratio of aerodynamic to inertial forces. Considering that the flapping inertia is dimensionally proportional to [ρ_{m}][l]^{5}, where ρ_{m} is the material density and l a characteristic length, matching the Lock number can be obtained by simply matching the material density of the blade, i.e., ρ_{mM}=ρ_{mP}. A similar definition of the Lock number can be developed for the fore–aft motion of the rotor due to the flexibility of the tower, leading to the same conclusion.
The system ith nondimensional natural frequency is defined as ${\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{i}={\mathit{\omega}}_{i}/\mathrm{\Omega}$, where ω_{i} is the ith dimensional natural frequency. Matching the lowest N nondimensional frequencies means that the corresponding eigenfrequencies in the scaled and fullscale system have the same relative placement among themselves and with respect to the harmonic excitations at the multiple of the rotor harmonics. In other words, the two systems have the same Campbell diagram (Eggleston and Stoddard, 1987). In addition, by matching nondimensional frequencies, the ratio of elastic to inertial forces is correctly scaled. Considering that the bending natural frequency of a blade is dimensionally proportional to $\sqrt{EJ/{\mathit{\rho}}_{\mathrm{m}}{l}^{\mathrm{6}}}$, the matching of nondimensional natural frequencies implies $(EJ{)}_{\mathrm{M}}=(EJ{)}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{6}}/{n}_{\mathrm{t}}^{\mathrm{2}}$, which is the same result obtained in the steady case for the matching of static deflections under aerodynamic loading. The same conclusions are obtained when considering deformation modes other than bending, so that in general one can write ${K}_{\mathrm{M}}={K}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{6}}/{n}_{\mathrm{t}}^{\mathrm{2}}$, where K is stiffness. Here again, it can be concluded that, for each given n_{l} and n_{t}, one can match the frequencies by adjusting the stiffness of the scaled model.
It should be remarked that this condition only defines the stiffnesses that should be realized in the scaled model, not how these are actually obtained. As noted earlier, it is typically difficult if not impossible to simply zoom down a complex realistic structure, and the model design may require a different configuration and choice of materials (Busan, 1998). An optimizationbased approach to the structural matching problem is described later in this work.
It is worth noting that matching both the Lock number and the placement of nondimensional natural frequencies implies that structural deflections caused by aerodynamic loads are correctly scaled. In fact, the Lock number is the ratio of aerodynamic to inertial forces, while ${\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{i}^{\mathrm{2}}$ is proportional to the ratio of elastic to inertial forces. Therefore, if both ratios are preserved, then $\mathit{\text{Lo}}/{\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{i}^{\mathrm{2}}$, being the ratio of aerodynamic to elastic forces, is also preserved. In symbols, this ratio is written as
where the righthand side is indeed proportional to the nondimensional tip deflection $\stackrel{\mathrm{\u0303}}{s}=s/R$ of a clamped beam subjected to a distributed load $q={C}_{\mathrm{L},\mathit{\alpha}}\mathit{\rho}c(R\mathrm{\Omega}{)}^{\mathrm{2}}$.
The matching of frequencies is also relevant to the matching of transient vorticity shedding in the wake, as mentioned earlier. In fact, assume that the blade flapping motion can be expressed as the single mode $d={d}_{\mathrm{0}}{e}^{{\mathit{\omega}}_{\mathrm{f}}t}$, where d is the flapping displacement and ω_{f} the flapping eigenfrequency. Then, the term $(\dot{d}/V{)}^{\prime}$ of Eq. (9) becomes
where ${\stackrel{\mathrm{\u0303}}{\mathit{\omega}}}_{\mathrm{f}}={\mathit{\omega}}_{\mathrm{f}}/\mathrm{\Omega}$ is the nondimensional flapping frequency. This term is matched between the scaled and fullscale models if the nondimensional flapping frequency is matched.
2.3 Subscaling criteria
As shown earlier, scaling is essentially governed by two parameters: the geometric scaling factor n_{l} and the time scaling factor n_{t}. No matter what choice is made for these parameters, the exact matching of some nondimensional parameters can always be guaranteed; these include nondimensional time, TSR, and Strouhal and Rossby numbers. In addition, the matching of other nondimensional quantities can be obtained by properly scaling some model parameters, again independently from the choice of n_{l} and n_{t}. For example, selecting the material density as ρ_{mM}=ρ_{mP} enforces the matching of the Lock number, whereas scaling the stiffness as ${K}_{\mathrm{M}}={K}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{6}}/{n}_{\mathrm{t}}^{\mathrm{2}}$ ensures the proper scaling of the system nondimensional natural frequencies. This way, several steady and unsteady characteristics of the fullscale system can be replicated by the scaled system. Other quantities, however, cannot be simultaneously matched, and one has to make a choice.
Table 1 summarizes the main scaling relationships described earlier. The reader is referred to the text for a more comprehensive overview of all relevant scalings.
The choice of the scaling parameters n_{l} and n_{t} is highly problem dependent. Indeed, given a fullscale reference, n_{l} is set by the size of its scaled replica, which is usually predefined to a large extent. For instance, the choice of the subscale size for a wind tunnel model depends on the characteristics of the target tunnel, to limit blockage (Barlow et al., 1999). When scaling down to a utility size, one typically chooses to reblade an existing turbine (Berg et al., 2014; Resor and Maniaci, 2014), thereby setting the scaling factor. The choice of n_{t} is often not straightforward and typically implies tradeoffs among quantities that cannot all be simultaneously matched.
For example, when the effects of gravity have to be correctly represented by the scaled model, the matching of the Froude number must be enforced. By setting Fr_{M}=Fr_{P}, one obtains the condition for the time scaling factor ${n}_{\mathrm{t}}=\sqrt{{n}_{\mathrm{l}}}$. Having set n_{t}, the scalings of all nondimensional parameters can now be expressed in terms of the sole geometric scaling factor n_{l}.
Another example is given by the design of smallscale wind turbine models for wind tunnel testing, which typically leads to small geometric scaling factors n_{l}. Bottasso et al. (2014) defined an optimal scaling by minimizing the error in the Reynolds number and the acceleration of scaled time. The latter criterion was selected to relax the requirements on closedloop control sampling time: since ${\mathit{\text{Re}}}_{\mathrm{M}}={\mathit{\text{Re}}}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{2}}/{n}_{\mathrm{t}}$, small geometric scaling factors might require very fast scaled times and hence high sampling rates, which could be difficult to achieve in practice for closedloop control models. Bottasso and Campagnolo (2020) used a different criterion, where a best compromise between the Reynolds mismatch and power density is sought. In fact, power density (defined as power P over volume or, in symbols, ${\mathit{\rho}}_{\mathrm{P}}=P/{R}^{\mathrm{3}}$) scales as ${\mathit{\rho}}_{{\mathrm{P}}_{\mathrm{M}}}/{\mathit{\rho}}_{{\mathrm{P}}_{\mathrm{P}}}={n}_{\mathrm{l}}^{\mathrm{2}}/{n}_{\mathrm{t}}^{\mathrm{3}}$ and, hence, increases rapidly for small n_{t}. For small n_{l} it becomes increasingly difficult, if not altogether impossible, to equip the scaled models with functional components (e.g., drivetrain, generator, actuation systems, sensors) that fit in the dimensions prescribed by the scaling factors. The adoption of larger components can be acceptable or not, depending on the nonphysical effects that are generated by their bigger dimensions and the goals of the model.
Yet another example of how delicate these choices can be is found in the experiments described by Kress et al. (2015). In this work, a scaled rotor was designed for experiments in a water tank, with the goal of comparing upwind and downwind turbine configurations. The rotor of the model was scaled geometrically from a fullscale reference; however, the same scaling ratio could not be used for the nacelle because of the need to house the necessary mechanical components. As a result, the model was equipped with an unrealistically large nacelle that, combined with the lower Reynolds number (which causes a thicker boundary layer), likely increased the redirection of the flow towards the outerblade portions in the downwind configuration. In turn, this led to the conclusion that nacelle blockage improves power production in downwind rotors. Although this may be true for the scaled experiment, there is little evidence that the same conclusion holds for a fullscale machine (Anderson et al., 2020). Because of miniaturization constraints, a larger nacelle is also used in the TUM G1 scaled turbine (Bottasso and Campagnolo, 2020), a machine designed to support wake studies and wind farm control research. The effects of the outofscale nacelle on the wake have, however, been verified and appear in this case to be very modest (Wang et al., 2020).
Additionally, particular combinations of n_{l} and n_{t} can make it difficult to find suitable designs. A clear example is found in the structural redesign of an aeroelastically subscaled blade. Indeed, as previously discussed, the scaled blade should have stiffnesses that scale as ${K}_{\mathrm{M}}={K}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{6}}/{n}_{\mathrm{t}}^{\mathrm{2}}$ and a mass density that scales as ρ_{mH}=ρ_{mP} to ensure the same nondimensional frequencies and Lock number. Depending on the values of the scaling parameters chosen, these scaling relationships might lead to unconventional combinations of stiffness and mass properties, which can be challenging to fulfill as shown in the next section.
Upscaling is a design effort driven by different criteria, including, among others, annual energy production (AEP), cost of material and manufacturing, logistics, and transportation. The situation is different for subscaling. In fact, the previous section has clarified the scaling relationships that exist between a fullscale system and its scaled model. The analysis has revealed that in general several steady and unsteady characteristics of the original system can be preserved in the scaled one. The question is now how to design such a scaled model in order to satisfy the desired matching conditions. This problem is discussed in this section.
3.1 Straightforward zooming down
This approach is based on the exact geometric zooming of the blade, including both its external and internal shape, and it has been advocated by Loth et al. (2017).
Regarding the external blade shape, geometric zooming implies that the same airfoils are used for both the scaled and the fullscale models. The mismatch of the Reynolds number (which is ${\mathit{\text{Re}}}_{\mathrm{M}}={\mathit{\text{Re}}}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{3}/\mathrm{2}}$ for Froude scaling) may imply a different behavior of the polars, especially for large values of n_{l}. On the other hand, as shown earlier, a geometric scaling ensures the near matching (up to the effects due to changes in the polars) of various characteristics, such as optimum TSR, nondimensional circulation, rotational augmentation, and vorticity shedding.
Regarding the internal blade shape, the skin, shear webs, and spar caps are also geometrically scaled down when using straightforward zooming. It should be noted that, for large geometric scaling factors n_{l}, the thickness of elements such as the skin or the shear webs may become very thin, possibly less than typical composite plies.
The zoomed scaling has to satisfy two constraints on the properties of the materials used for its realization.
The first constraint is represented by the matching of material density (ρ_{mM}=ρ_{mP}), which is necessary to ensure the same Lock number. It should be remarked that the overall material density of the blade includes not only the density of the main structural elements but also contributions from coatings, adhesive, and lightning protection. These components of the blade may not be simply scaled down, so this problem may deserve some attention.
The second constraint is represented by the scaling of the stiffness, which is necessary for ensuring the matching of nondimensional natural frequencies. For Froude scaling, stiffness changes as ${K}_{\mathrm{M}}={K}_{\mathrm{P}}{n}_{\mathrm{l}}^{\mathrm{5}}$. Considering bending, the stiffness is K=EJ. For a blade made of layered composite materials, the bending stiffness is more complicated than the simple expression EJ, and it will typically need to be computed with an ad hoc methodology, for example using the anisotropic beam theory of Giavotto et al. (1983). However, the simple expression EJ is sufficient for the dimensional analysis required to understand the effects of scaling. Since the sectional moment of inertia J is dimensionally proportional to l^{4}, with l being a characteristic length of the blade cross section, this constraint requires Young's modulus to change according to E_{M}=E_{P}n_{l}. This implies that all materials used for the scaled blade, including the core, should have a lower level of stiffness than (and the same density as) the materials used at full scale; as shown later, this constraint is not easily met.
Since strain ϵ is defined as the ratio of a displacement and a reference length, it follows that ϵ_{M}=ϵ_{P}. Therefore, given that E_{M}=E_{P}n_{l}, it follows that σ_{M}=σ_{P}n_{l}, and the stresses in the scaled model are reduced compared to the ones in the fullscale model. Still, one would have to verify that the admissible stresses and strains of the material chosen for the scaled blade are sufficient to ensure integrity.
The critical buckling stress of a curved rectangular plate is
where k_{c} is a coefficient that depends on the aspect ratio of the panel, its curvature, and its boundary conditions; ν is Poisson's ratio; d is the panel thickness; and b is the length of the loaded edges of the plate (Jones, 2006). Here again, the expression of the critical stress of a layered anisotropic composite plate would be more complex than the one reported in Eq. (14), but this is enough for the present dimensional analysis. By using the scaling relationships for length and for E, Eq. (14) readily leads to ${\mathit{\sigma}}_{{\text{cr}}_{\mathrm{M}}}={\mathit{\sigma}}_{{\text{cr}}_{\mathrm{P}}}{n}_{\mathrm{l}}$. This means that if the fullscale blade is buckling free, so is the scaled one, as both the critical buckling stress and the stresses themselves scale in the same manner.
3.2 Aerostructural redesign
An alternative approach to the design of a subscale model is to identify an external shape and an internal structure that match, as closely as possible, the aeroelastic behavior of the fullscale blade. This approach offers more degrees of freedom, at the cost of an increased design complexity; indeed, one designs a new blade that, although completely different from the fullscale one, matches some of its characteristics.
In this second approach, the first step consists of defining a blade shape that can mimic the aerodynamic behavior of the fullscale system. As previously discussed, this can be obtained according to different criteria. Here, the following three conditions are considered. First, a new set of airfoils is selected to match as closely as possible, despite the different Reynolds number of operation, the polar coefficients of the airfoils of the fullscale blade; this is relevant for the matching of the performance and loading of the rotor. Second, the two rotors should have similarly shaped power coefficient curves, which is relevant for performance on and off the design point. Finally, the blades should have the same spanwise circulation distribution, which is relevant for a similar aerodynamic loading of the blade and wake behavior. The resulting scaled blade shape (both in terms of cross sections, because of the changed airfoils, and in terms of chord and twist distributions) will be different from the fullscale rotor. However, this is clearly irrelevant, as the goal is to match some quantities of interest between the two rotors, not their shape.
The aerodynamic design problem can be formally expressed as
Vector p_{a} indicates the aerodynamic design variables, which include the chord and twist distributions c(η) and θ(η), appropriately discretized in the spanwise direction, while J_{a} is a design figure of merit, m_{a} are matching constraints, and finally c_{a} are additional design conditions. This formulation of the aerodynamic design problem is very general, and different choices of the figure of merit and of the constraints are possible, depending on the goals of the scaled model.
In the present work, the aerodynamic optimization cost function is formulated as
This cost drives the design towards the power coefficient of the target fullscale model ${\widehat{C}}_{\mathrm{P}}$ at ${N}_{{C}_{\mathrm{P}}}$ control stations. This cost function ensures that the subscale model – whose airfoils generally present a reduced efficiency due to the lower chordbased Reynolds number – has a C_{P} that is as close as possible to the fullscale model. Using ${N}_{{C}_{\mathrm{P}}}=\mathrm{1}$ leads to a design with the best C_{P} at the TSR λ_{1}.
Within the vector of matching equality constraints, one set of conditions enforces the matching of the spanwise distribution of the circulation $\widehat{\mathrm{\Gamma}}$ at N_{Γ} control stations:
where $\widehat{(\cdot )}$ indicates in general a tobematched scaled quantity of the target fullscale model. Another constraint may be added to prescribe the maximum power coefficient to take place at the same design TSR, i.e., ${\mathit{\lambda}}_{max\left({C}_{\mathrm{P}}\right)}={\mathit{\lambda}}_{max\left({\widehat{C}}_{\mathrm{P}}\right)}$. Finally, vector c_{a} specifies additional design inequality constraints, which may include a margin to stall, maximum chord, and others, depending on the application.
Once the new aerodynamic shape is identified, the second step consists in the design of an internal blade structure that can mimic the fullscale aeroelastic behavior while ensuring integrity and satisfying manufacturing and realizability constraints. This approach allows for more freedom than the zoomingdown approach; for example, one can use different materials than the ones used for the fullscale design, and nonstructural masses can be added without affecting the matching characteristics of the scaled blade.
The structural design problem can be formally expressed as
Vector p_{s} indicates the structural design variables, which include the size of the various blade structural elements (skin, spar caps, shear webs, and leading and trailingedge reinforcements), discretized span and chordwise. Here again, this formulation is very general, and specific goals will lead to different choices of the merit function and of the constraints.
For example, assuming the blade to be modeled as a beam, the structural optimization cost can be formulated as
where w_{s} is a tuning weight, M_{p} and K_{q} are elements of the mass and stiffness matrices, and the sets 𝒮_{M} and 𝒮_{K} identify the elements that should be considered within the generally fully populated symmetric mass and stiffness matrices. The first term in the cost aims at the matching of the scaled target mass distribution, while the second aims at the stiffness distribution. Vector m_{s} indicates the matching equality constraints. These may include the matching of a desired number of natural frequencies ${\mathit{\omega}}_{i}={\widehat{\mathit{\omega}}}_{i}$ and the matching of a desired number of mode shapes and/or static deflections ${\mathit{u}}_{j}\left({\mathit{\eta}}_{i}\right)={\widehat{\mathit{u}}}_{j}\left({\mathit{\eta}}_{i}\right)$ at a given number of spanwise stations η_{i}. Finally, vector c_{s} specifies the additional design inequality constraints. These constraints express all other necessary and desired conditions that must be satisfied in order for the structural design to be viable and in general include maximum stresses and strains for integrity, maximum tip deflection for safety, buckling, and manufacturing and technological conditions.
It should be noted that the matching of the scaled beam stiffness and mass distributions – if it can be achieved – is an extremely powerful condition. In fact, a geometrically exact nonlinear beam model is fully characterized entirely in terms of its reference curve, stiffness, and mass matrices (Bottasso and Borri, 1998). This means that exactly matching all of these quantities would ensure the same nonlinear structural dynamic behavior of the fullscale target. As shown later, this is not always possible because of limits due to technological processes, material characteristics, chosen configuration of the scaled model, etc. In this case, there is a partial match between the fullscale and scaled beam models, and the sets 𝒮_{M} and 𝒮_{K} include only some elements of the mass and stiffness matrices. When this happens, additional matching constraints can help in ensuring as similar a behavior as possible between the scaled and fullscale structures, for example by including static deflections and/or modal shapes, as shown later.
The two strategies of straightforward zooming and aerostructural redesign are applied here to the subscaling of a 10 MW machine, developed in Bottasso et al. (2016) as an evolution of the original Danmarks Tekniske Universitet (DTU) 10 MW reference wind turbine (Bak et al., 2013). The main characteristics of the turbine are reported in Table 2. Some of the principal blade characteristics are given in Table 3, which reports the position of the airfoils, whereas Table 4 details the blade structural configuration and Table 5 summarizes the material properties.
Three different subscalings are considered here. The first subscale model, denominated the W model, is based on the German WINSENT test site (ZSW, 2016), which is equipped with two 750 kW turbines with a rotor diameter of 54 m (ZSW, 2017). The reference rotor blades are scaled down to match the span of the WINSENT blades; reblading one of the WINSENT turbines yields a subscale model of the fullscale 10 MW turbine suitable for field testing. The second model, denominated the S model, is based on the SWiFT test site, which is equipped with Vestas V27 turbines. Here, the fullscale rotor is scaled down to a diameter of 27 m. Finally, the T model is a wind tunnel model with a rotor diameter of 2.8 m, which is similar to the scaled floating turbine tested in the Nantes wave tank in the INNWIND.EU project (Azcona et al., 2016).
Table 6 reports the different geometric scaling factors and a few additional key quantities of the three subscale models. For all, Froude scaling is used, which sets the timescale factor as previously explained. The application of the scaling laws to the fullscale turbine results in the characteristics listed in Table 7. Independently of the approach chosen to define the internal and external shape, the scaled models must fulfill these conditions to correctly mirror the dynamic behavior of the fullscale wind turbine.
The gravoaeroservoelastic scaling laws lead to very light and flexible subscale blades. For instance, the standard blades of the V27 weigh 600 kg (Vestas, 1994), which is 4 times more than the gravoaeroservoelastically scaled blades of the S model. It should, however, be remarked that this ratio would be smaller for a modern blade, since the V27 was designed more than 25 years ago and its blades are heavier than the ones based on contemporary technology.
The following sections detail the design of the external and internal shape of the three subscale blades. Section 4.1 describes the aeroservoelastic and design tools used to this end. Then, Sects. 4.2 and 4.3 discuss, respectively, the strengths and limitations of each design strategy for each subscale model.
4.1 Aeroservoelastic and design tools
The aeroservoelastic models are implemented in CpLambda
(Bottasso et al., 2012). The code is based on a multibody formulation for flexible systems with general topologies described in Cartesian coordinates. A complete library of elements – including rigid bodies, nonlinear flexible elements, joints, actuators, and aerodynamic models – is available, as well as sensor and control elements.
The aerodynamic characteristics of the blade are described through lifting lines, including spanwise chord and twist distribution and aerodynamic coefficients. The code is coupled with aerodynamic models based on the BEM model, formulated according to streamtube theory with annular and azimuthally variable axial and swirl inductions, unsteady corrections, root and blade tip losses, and a dynamic stall model.
The tower and rotor blades are modeled by nonlinear, geometrically exact beams of arbitrary initially undeformed shapes, which are bending, shear, axial, and torsion deformable. The structural and inertial characteristics of each beam section are computed with ANBA
(Giavotto et al., 1983), a 2D finiteelement crosssectional model. Finally, fullfield turbulent wind grids are computed with TurbSim
(Jonkman et al., 2009) and used as input flow conditions for the aeroservoelastic simulations.
CpMax
(Bortolotti et al., 2016) is a design framework wrapped around CpLambda
, which implements optimization algorithms to perform the coupled aerostructural design optimization of the blades and, optionally, of the tower. For the present work, the code was modified to implement also the scaled design matching optimizations defined by Eqs. (15) and (18). All optimization procedures are solved with a sequential quadratic programming algorithm, in which gradients are computed by means of finite differences.
4.2 External shape design
For all three models, the design of the subscale external blade shape aims at replicating the aerodynamic characteristics of the fullscale rotor, including its wake. As long as the chordbased Reynolds numbers are sufficiently large, a zoomingdown approach is clearly the simplest strategy for designing the external shape of a scaled blade.
Airfoil FFAW3241 equips the outermost part of the fullscale blade (see Table 3). Its performance at the three typical Reynolds numbers of the fullscale, W, and S models was computed with ANSYS Fluent
(ANSYS, Inc., 2019). The results are reported in Fig. 1. The performance of the airfoil is clearly affected by the Reynolds number, with a particularly significant drop in efficiency for the lowest Reynolds case. Notwithstanding these Reynolds effects, the zoomingdown approach is selected for the W and S models, since the airfoils are still performing well at their corresponding typical subscale Reynolds number. A redesign approach with alternative airfoils was not attempted here, and would probably lead only to marginal improvements of the aerodynamic performance.
On the other hand, for the small geometric scaling factor of the T model, the aerodynamic redesign approach is necessary. In general, smooth airfoils present a large reduction in aerodynamic efficiency below a critical Reynolds number of about 70 000 (Selig et al., 1995). Efficient profiles specifically developed for lowReynoldsnumber applications are generally necessary in order to get a good matching of the fullscale aerodynamic performance. As an alternative to the original airfoil, the 14 % thick airfoil RG14 (Selig et al., 1995) is selected, because its aerodynamic characteristics at the scaled Reynolds number are in reasonable agreement with the ones of the original airfoil at its fullscale Reynolds number (Fig. 1). The blade is then completely redesigned, using the RG14 airfoil along its full span.
The blade shape is parameterized by means of chord and twist spanwise distributions. The design problem is formulated as the maximization of the power coefficient at the design TSR λ_{d} of the fullscale rotor, solving Eq. (15) with the cost given by Eq. (16) for ${N}_{{C}_{\mathrm{P}}}=\mathrm{1}$ and λ_{1}=λ_{d}. The nonlinear constraints expressed by Eq. (17) enforce the same spanwise nondimensional circulation distribution of the fullscale blade.
Figure 2 shows the external shapes of the fullscale blade and the three subscale models in terms of chord, relative thickness, twist, and Reynolds number. Clearly, the shape curves for the W and S models overlap with the fullscale ones, because zooming is used in these two cases, as previously explained.
The three subscale models have the same TSR in region II as the fullscale machine and the correspondingly subscaled rated rotor speeds. The rated wind speeds do not exactly match the subscale ones, on account of the differences in the C_{P}TSR curves caused by the Reynolds effect.
4.3 Design of the internal structure
The definition of the internal structure has to achieve the following goals: the matching of the fullscale aeroelastic behavior, the integrity of the blade under loading, and the feasibility of the manufacturing process. In the next two sections, the zoomingdown and the redesign approaches are applied to the structure of the three subscale blades.
4.3.1 Limits of the zoomingdown approach
The straightforward zoomingdown approach can be applied to the internal structure of the W and Smodel blades, as their external geometrical shape has also been defined following this approach. The resulting structures satisfy all scaling constraints but present some critical challenges.
First, the thicknesses of some of the components are unrealistically low. The blade root of the W model is, for example, only 20 mm thick and is therefore unable to accommodate the rootbolted connections. Furthermore, the scaling of the outer shell skin leads to a laminate thickness of less than one ply. The third web of the Smodel blade is also extremely thin (less than 1 mm) and very close to the trailing edge.
Additionally, the scaled structure requires materials characterized by very peculiar mechanical properties. Indeed, as previously shown, the scaling laws require the modulus of elasticity to obey the relationship E_{M}=E_{P}n_{l} and the material density to be ρ_{mM}=ρ_{mP}. For example, the outer shell of the Wmodel blade requires an elasticity modulus of 6.6 GPa and a density of 1845 kg m^{−3}, which are not typical values of conventional materials (see Fig. 3). Finally, nonstructural masses – such as glue, paint, and lightning protection – cannot be exactly zoomed down by geometric scaling and need to be treated separately.
One may try to relax some of these hurdles by increasing the necessary component thicknesses and choosing materials with mechanical properties that compensate for this increase. For example, a 3fold increase of the skin thickness in the W model would be able to accommodate the rootbolted connection and would satisfy manufacturing tolerances. To meet the mass and inertia constraints, a material should be used that has a lower density, ${\mathit{\rho}}_{\text{mM}}={\mathit{\rho}}_{\text{mP}}/\mathrm{3}$, and a lowerelasticity modulus, ${E}_{\mathrm{M}}={E}_{\mathrm{P}}{n}_{\mathrm{l}}/\mathrm{3}$. Figure 3 reports Ashby's diagram of Young's modulus vs. density (Cambridge University Engineering Department, 2003). In this plot, the values corresponding to the outer shell skin materials have been marked with × symbols. A red symbol indicates the fullscale blade, a yellow symbol is used for the W model considering the exact zoomingdown approach, and a green symbol indicates the solution with a 3fold thickness increase. It should be noted that, although the properties of the scaled models do correspond to existing materials, these are typically not employed for the manufacturing of blades. Therefore, their actual use for the present application might indeed pose some challenges.
Overall, the zoomingdown approach for the structural design is not really straightforward and is significantly more complicated than in the case of the aerodynamic design. An alternative is offered by a complete redesign of the internal structure, which is illustrated in the next section.
4.3.2 Redesign of the W and S models
An alternative to the zoomingdown approach is the redesign of the internal structure. This consists of a typical blade design process, subjected not only to additional constraints that enforce the desired scaling relationships but, crucially, also to all other conditions that are necessary to make the design viable. For example, here a lower bound to the thickness of all structural components is set to 1 mm, while a minimum thickness of 60 mm is assumed at the root to accommodate the bolted connection of the W and S models.
Additionally, one has greater freedom in the choice of materials. For the present applications, the glassfiberreinforced plastic (GFRP) composites of the fullscale blade appear to be suitable choices also for the W model. On the other hand, these materials are too stiff for the S model, due to its smaller geometric scaling. An alternative was found within the family of thermoplastic materials that have typical stiffness values between 1–3 GPa and densities between 900 and 1400 kg m^{−3} (Brondsted et al., 2005). Although not strictly of interest here, thermoplastics also have interesting advantages over thermosets, such as reduced cycle times, lower capital costs of tooling and equipment, smaller energy consumption during manufacturing, and enhanced recyclability at the end of their life (Murray et al., 2018).
During the design phase of the subscale models, more careful attention can also be paid to the distributions of nonstructural masses. Specifically, masses from shell and sandwich cores must be recomputed for the new scaled structure in order to prevent buckling of the sandwich panels. Additional masses from surface finishing and painting are also recomputed according to the surface of the external shell. In fact, if a zoomingdown strategy is chosen for the design of the external geometry, these masses will scale with the length scale factor. Masses from resin uptake in the outer shell and shear webs are recomputed for the scaled structure assuming a constant area density. Indeed, this value does not change from the full scale to the subscale, since it depends on the material and manufacturing process. A different assumption is taken for the masses of bonding plies and adhesive along the shear webs and leading and trailing edge. Since these masses are chordwise dependent, the linear density of these materials in the subscale size must be corrected by the length scale factor. Finally, the linear density of the lightning protection system is assumed to be constant for all sizes.
The structural design is formulated as the matching optimization problem expressed by Eq. (18). The cost function given by Eq. (19) considers the sole spanwise matching of the mass distribution, i.e., it neglects inertia terms in 𝒮_{M} and uses w_{s}=0. The matching constraints m_{s} include the lowest three natural frequencies, and the static deflected shape of the outboard 40 % section of the blade. This static condition was chosen to represent the maximum tip displacement resulting from turbulent simulations in power production for the fullscale machine (design load case (DLC) 1.1; see IEC (2005)). Finally, the additional design constraints c_{s} include stresses, strains, fatigue and technological constraints in the form of bounds on thickness and thickness rate of change of the laminates.
The structural design for the W and S models is based on a typical thinwalled composite configuration, where the design variables are defined as the spanwise thicknesses of the skin, shear webs, spar caps, and leading and trailingedge reinforcements. Given the smaller size of the scaled blades, one single shear web is used instead of the three used in the fullscale 10 MW model. Table 8 describes the mechanical properties of the materials used for these two blades, while Table 9 associates the various structural elements with the materials.
For the S model, the thermoplastic materials polymethyl methacrylate (PMMA) and polyoxymethylene (POM) are chosen because of their lower level of stiffness. The use of polymer materials reduces the nonstructural masses, as the adhesive is no longer necessary. Due to the reduced fatigue characteristics of these materials, the blade lifetime is limited to 5 years. This is assumed to be acceptable in the present case, given the research nature of these blades. Constraints on maximum stresses and strains are satisfied with an ample margin for these blades. However, the inclusion of a larger set of DLCs (including extreme events and parked conditions) might create more challenging situations, which could increase the requirements regarding material strength, possibly eventually leading to the selection of different materials.
Figure 4 reports the internal structure of the W and S models, as well as the overall mass distributions, including realistic nonstructural masses. The scaled mass distribution follows quite closely the reference one along the blade span, with the exception of the root because of the additional thickness that must be ensured to accommodate the bolted connection. The blade satisfies the scaling inertial and elastic constraints within a tolerance of less than 5 %.
4.3.3 Redesign of the T model
The very small size of the wind tunnel model blade prevents the use of a typical thinwalled solution. Following Bottasso et al. (2014) and Campagnolo et al. (2014), this scaled blade is not hollow but presents a full cross section obtained by machining a foamy material. Two unidirectional spar caps provide the required flapwise stiffness distribution. The surface smoothness is obtained by a very thin layer of skin made of glue. Although Bottasso et al. (2014) and Campagnolo et al. (2014) considered different scaling laws, their blade design configuration was found to be a suitable choice even in the present gravoaeroservoelastic scaling exercise. The selection of appropriate materials represents a critical aspect of the problem, and the mechanical properties listed in the Cambridge University Materials Data Book (Cambridge University Engineering Department, 2003) were used to guide the material selection for the spar caps and core. A rigid polymer foam is chosen as filler, because of its relatively high level of stiffness and lightness. For the spar caps, thermoplastic polymers are again found to be the most suitable solution even though their stiffnesstodensity ratio is much lower than materials traditionally used for spar caps. Moreover, the use of thermoplastics allows for alternative and simpler manufacturing processes, leading to a higher flexibility in the spar cap design. From this family of materials, polypropylene is chosen because of its low stiffness modulus. Finally, the external shell is covered by a very thin layer of the epoxy structural adhesive Scotch Weld AF 32 (3M Adhesives Division, 2000).
The design variables are represented by the spanwise thickness and width of the two spars. The design problem is formulated according to the constrained matching optimization expressed by Eq. (18). The cost function of Eq. (19) considers the spanwise mass distribution in 𝒮_{M} and the flapwise stiffness distribution in 𝒮_{K}. The matching constraints m_{s} include the lowest three natural frequencies and the flapwise static extreme tip deflection. Both the cost and the constraints only consider the flapwise characteristics of the blade, because the structural configuration consisting of a solid core and two spar caps allows for limited control of the edgewise characteristics. As a result, the scaled blade presents a higher level of edgewise stiffness than the fullscale reference.
Figure 5 reports the results of the design optimization. The desired matching of mass and flapwise stiffness is achieved, except at blade root. Even though the placement of the first flapwise natural frequency with respect to the rotor speed is ensured, the constraint on the lowest edgewise natural frequency could not be exactly matched due to the large chord. Small disparities in mass distribution introduce a difference of about 1 % in the blade flapping inertia.
In this section, the behavior of the scaled models is compared to the fullscale machine. The main goal here is to assess to what extent the subscale models are capable of successfully mirroring relevant key characteristics and load trends of the fullscale reference.
The same collectivepitch/torque controller governs all machines. The controller uses a lookup table for torque to operate at rated TSR in region II and a proportional–integral–derivative (PID) pitch loop to maintain constant rated power in region III. The PID gains used for the scaled models are obtained by transforming the ones of the fullscale machine using the scaling laws, and the regulation trajectory is adapted to each model to account for differences in the C_{P}TSR curves. Notice that the scaling of gains is a conservative approach: in the case of an exact matching at scale of all aeroelastic characteristics of the turbines, the use of a scaled controller will also ensure an identical closedloop response. However, if the scaled models do not exactly represent the fullscale reference – which is invariably the case in practice – an ad hoc retuned controller (i.e., a controller specifically optimized for the scaled model) will in general have better performance than the one obtained by the scaling of the gains. The choice of gain scaling instead of retuning was made here to consider a worstcase scenario.
5.1 Relevant key indicators
The models are simulated in a power production state at five different wind speeds from cutin to cutout. The winds of the scaled simulations are obtained by velocity scaling the turbulent winds used for the fullscale machine (i.e., the integral space and timescales are both correctly scaled). The matching between the scaled and fullscale turbines is assessed with the help of 10 different indicators: AEP; maximum flapwise tip displacement (MFTD); maximum thrust at main shaft (ThS); maximum combined blade root moment (CBRM); maximum flapwise bending root moment (FBRM); maximum edgewise bending root moment (EBRM); and the Weibullaveraged damage equivalent load (DEL) for ThS, CBRM, FBRM, and EBRM.
5.1.1 Utilityscale models
As previously discussed, the design both of the external shape and of the internal structure may induce differences in the behavior of a scaled model with respect to its fullscale reference. To better understand the effects of these differences and their origins, three different sets of results are presented in Fig. 6.
The first plot (a) compares the indicators of the fullscale turbine with the upscaled ones of the W and S models. Both the internal structure and the external shape are obtained by zooming down, and Reynolds effects are accounted for by CFDcomputed polars. Although a zoomeddown structure cannot really be a practical solution – as discussed earlier – because of excessively thin structural elements or the need for peculiar material properties, this solution is shown here because it highlights the sole effects of the Reynolds mismatch. In other words, since this is a purely numerical study, the thicknesses and mechanical properties were used exactly as produced by scaling, resulting in a nearly exact satisfaction of the matching of all structural characteristics. Therefore, the differences of the indicators between the fullscale and scaled models shown in this plot can be entirely attributed to Reynolds effects. The fullscale and utilitysize models are equipped with airfoil polars at different Reynolds numbers computed with the CFD code ANSYS Fluent
(ANSYS, Inc., 2019).
The second plot (b) compares the indicators for the W and S models featuring a zoomeddown external shape (which neglects Reynolds effects) and a redesigned internal structure. Although Reynolds effects would, in reality, be present, by neglecting them here – which is again possible because this is a purely numerical study – one can assess from this solution the sole effects of the structural redesign on the matching of the indicators.
Finally, the third and last plot (c) considers the solution obtained by zooming down the aerodynamic shape, considering Reynolds effects and a redesigned internal structure. As argued earlier, this is indeed the solution that is practically realizable, and, therefore, these are the more realistic results of the set considered here. Hence, differences between the fullscale and scaled models are due to mismatches caused by both the Reynolds number and the redesign procedure.
As expected by the size difference, results shown in the first plot suggest a larger effect of the Reynolds number mismatch for the S model than for the W model. This results in a drop in all indicators because of the decreased airfoil efficiency.
The second plot shows a similar matching for both models. Indeed, most of the key loads are matched within 5 % for both the W and the S model. A larger difference between the two models is found for EBRM and DEL EBRM, which are only poorly matched by the W model, whereas they are quite accurate for the S model. The mismatch is due to a slightly higher sectional mass in the last 20 % of the blade of the W model, as shown in Fig. 4. A significant difference with respect to full scale is also observed for the maximum flapwise tip displacement of both the W and S models. This difference is caused by a slightly different dynamic behavior induced by mismatches in the flapwise and torsional stiffness distributions. Even though FBRM matches very well for both the W and S model at the root, these differences lead to a poorer match at sections toward the blade tip, which in the end impacts MFTD.
Overall, both models are capable of matching the key indicators of the fullscale target reasonably well, considering both Reynolds effects and a redesigned structure.
5.1.2 Wind tunnel model
The behavior of the T model is compared with the 10 MW baseline in Fig. 7. The additional indicator, maximum edgewise tip deflection (METD), is considered in this case. The polars for the T model are computed with Xfoil (Drela, 2013).
The comparison shows satisfactory behavior of the wind tunnel model for most key indicators, notwithstanding the very different Reynolds numbers (about 10^{7} for the fullscale reference and about 2×10^{4} for the T model). As expected, the largest mismatch is found for the maximum edgewise tip displacement. This can be justified by the inability of the structural design variables (limited to the two caps) in controlling the edgewise stiffness.
5.2 Load trends in waked conditions
Scaled models can also be used to capture trends, instead of absolute values. Indeed, the goal of scaled testing is often to understand the trends generated on some metric by, for example, a control technology or by a particular operating condition or other factors, whereas the exact quantitative assessment of the induced effects must be left to a final fullscale verification.
As an example of the analysis of trends, the scaled models designed here are used to explore changes in loading between unwaked and waked inflow conditions. To this end, the fullscale turbine is simulated with an average inflow velocity of 7 m s^{−1}, considering a shear exponent of 0.2 and a turbulence intensity of 8 %. The wake deficit generated by an upstream 10 MW machine is then added to this inflow, in order to simulate a waked condition. The wake is modeled by the superposition of a turbulent wind grid generated with TurbSim
(Jonkman et al., 2009) and the firstorder solution of the deficit of the Larsen model (EWTSII model) (Bottasso et al., 2017). The downstream turbine is located at a longitudinal downstream distance of 4D from the upstream machine, and its lateral distance from the wake center is varied from −1.25D (right, looking downwind) to 1.25D (left), realizing different degrees of wake–rotor overlap. The scaled models are simulated by velocityscaling the fullscale inflows. The key indicators considered are AEP; ThS; FBRM; and DELs for CBRM, FBRM, and EBRM.
Figure 8 reports changes in key indicators at several degrees of wake overlap with respect to unwaked inflow conditions. The fullscale machine presents the largest reduction in AEP and ThS in full wake overlap. An asymmetrical load trend of the DELs for FBRM, EBRM, and CBRM is visible when the rotor is operating in partial wake. This behavior is mostly due to the rotor uptilt angle, which introduces an additional velocity component in the rotor plane. In fact, for a clockwiserotating (when looking downstream) rotor, this extra velocity component increases the inplane velocity at the blade sections when the blade is on the right side of the rotor (i.e., during the downstroke; here left and right are defined for an observer looking downstream). Additionally, when a wake impinges on the right side of the rotor, the outofplane velocity component decreases, because of the wake deficit. Both of these effects tend to decrease the angle of attack at the blade sections. On the other hand, when a wake impinges on the left portion of the rotor, the effect of the decreased outofplane component is in part balanced by the also decreased inplane component. Because of this different behavior, larger load fluctuations (and hence higher fatigue loads) are observed for right wake impingements than for left ones. A similar effect is caused by the elasticity of the tower: under the push of the thrust, the tower bends backwards, which in turn tilts the rotor upward, adding to the previously described phenomenon. Other minor effects are also due to the elastic deformations caused by gravity, which again contribute to breaking the symmetry of the problem.
Overall, the largest scaled models follow the trends very well, with the S model performing slightly better than the W model. Indeed, the W model is better than the S model when looking at Weibullaveraged quantities (Fig. 6), but the S model presents a slightly superior matching of blade loads at the specific speed at which the load trend study is performed. The trends are also reasonably captured by the smallerscale T model, but with significant differences in DEL FBRM. Specifically, there is an overestimation of this quantity around the −0.5D lateral wake center position. A detailed analysis of the results revealed this behavior to be caused by the blade operating at angles of attack close to the stalling point. This indicates another possible limit of models with largescale factors, whose airfoils may have very different stall and poststall behavior than their fullscale counterparts.
This paper has analyzed the scaling conditions that should be met by a subscale model to match a fullscale reference in terms of its full aeroservoelastic response. The analysis has shown that many relevant key aspects of the steady and unsteady response of a machine, considered as flexible, can indeed be matched. Part of this analysis can also be used to understand expected changes due to upscaling, which can be useful in the design of larger rotors. To the authors' knowledge, this is one of the most comprehensive analyses of the problems of scaling wind turbines presented thus far.
Within this framework, this paper has considered two alternative ways of designing a scaled rotor. The first is based on the idea of exactly zooming down the fullscale reference to obtain the subscale model. An alternative strategy is to completely redesign the rotor, from both an aerodynamic and a structural point of view. This produces a scaled blade that, although possibly very different from the fullscale one, matches some of its key characteristics as closely as possible.
These two alternative strategies have been tested on the gravoaeroservoelastic scaling of a conceptual 10 MW blade to three different subscale models: two utilityscale ones to be used for the reblading of small existing turbines and one for equipping a very small model turbine to conduct experiments in the controlled environment of a wind tunnel.
The following conclusions can be drawn from the application of the two strategies to these three different scaling problems.
The simplest strategy to design the external shape of utilityscale blades is the straightforward zoomingdown approach, as long as the subscale Reynolds number is sufficiently high. This strategy benefits from a simple implementation and leads to an acceptable match of the blade aerodynamic performance. However, when the blade aerodynamic performance is compromised by the Reynolds mismatch – which is the typical case of wind tunnel models – the alternative but more complex strategy of redesigning the aerodynamic shape becomes preferable if not altogether indispensable. Special lowReynoldsnumber airfoils may be used to mitigate the effects caused by the reduced Reynolds regime. However, different behavior at and around stall might lead to different loads when operating at large angles of attack.
The straightforward zooming down of the blade internal structure is instead typically very difficult for all scaling ratios. In fact, the need for materials of unusual characteristics and the nonscalability of nonstructural masses unfortunately hinder the applicability of this simple approach. An alternative is found in the structural redesign strategy, which offers more flexibility at the price of increased complexity. Even here, however, the problem is nontrivial. For example, materials may play a critical role, due to the very flexible nature of some of these scaled blades.
The aeroservoelastic analyses conducted herein have shown that, in general, it is not possible to exactly match all the characteristics of a fullscale machine with a subscale model. However, with the proper choices, some key indicators are nicely captured. In addition, changes in operating conditions are represented quite well even at the smaller scale. For example, it was shown that changes in loading from an unwaked to a waked condition are accurately represented by all scaled models, which successfully capture intricate and possibly unexpected couplings with design aspects such as nacelle uptilt and tower deflection. The good performance of the models in capturing such complex effects opens up a range of applications and use cases. For example, with the right design choices, scaled models can be employed to better understand rotor–wake interactions or test sophisticated control strategies at the turbine and/or plant levels.
Further improvements in the performance of the subscale models are certainly possible. Indeed, while some of the limitations result from the choice of quantities to be matched, others can be overcome by technological advances. For instance, improvements in measurement technology can relax the requirements on the scaling of time, allowing for a better match of other quantities. Additionally, advances in material and manufacturing may ease the application of unconventional materials; relax sizing constraints; and lead to more accurate, simpler, fastertodevelop, and cheaper models.
This work has exclusively focused on the wind turbine itself, and the effects of scaling have been quantified for the aerodynamic performance and loading of the rotor. The recent study of Wang et al. (2020) expands this analysis by considering the effects of scaling on wake behavior. Even in that case the conclusion is that properly scaled models can produce very realistic wakes.
Further work should focus on expanding the scope of the scaling analysis, introducing the effect of hydrodynamics. Indeed, as floating wind energy is expected to significantly grow in the coming years, it is becoming increasingly important to better understand which aspects of the aerohydroservoelastic response of these machines can be matched and how to best design subscale models. This is, however, only part of the problem. Research efforts are also necessary to better understand how to replicate the inflow conditions that fullscale machines face in various types of atmospheric and terrain conditions. This is a challenging task, since it requires a deep understanding of atmospheric flows, their interaction with the terrain orography and the vegetation, and technology to replicate these flows at scale.
It is the hope of the authors that the results shown in this paper will increase the confidence in scaled testing, in the belief that scaled model have a significant role to play in the advancement of wind energy science.
a  Axial induction factor 
a_{s}  Speed of sound 
c  Chord length 
d  Outofplane blade section flapping displacement 
f  Characteristic frequency 
g  Acceleration of gravity 
l  Characteristic length 
n_{l}  Geometric scaling factor, i.e., ${l}_{\mathrm{M}}/{l}_{\mathrm{P}}$ 
n_{t}  Time scaling factor, i.e., ${t}_{\mathrm{M}}/{t}_{\mathrm{P}}$ 
n_{Ω}  Angular velocity scaling factor, i.e., ${\mathrm{\Omega}}_{\mathrm{M}}/{\mathrm{\Omega}}_{\mathrm{P}}$ 
n_{v}  Wind speed scaling factor, i.e., ${V}_{\mathrm{M}}/{V}_{\mathrm{P}}$ 
p  Vector of design parameters 
r  Spanwise coordinate 
s  Tip deflection 
t  Time 
u  Characteristic speed 
A  Rotor disk area 
A_{b}  Blade planform area 
B  Number of blades 
C_{D}  Drag coefficient 
C_{L}  Lift coefficient 
C_{L,α}  Slope of the lift curve 
C_{P}  Power coefficient 
C_{T}  Thrust coefficient 
E  Young's modulus or airfoil efficiency, i.e., ${C}_{\mathrm{L}}/{C}_{\mathrm{D}}$ 
EJ  Bending stiffness 
Fr  Froude number 
I  Rotor polar moment of inertia 
I_{b}  Blade flapping inertia 
J  Cost function 
K  Stiffness 
Lo  Lock number 
M  Mass 
Ma  Mach number 
P  Aerodynamic power 
Q  Torque 
R  Rotor radius 
Re  Reynolds number 
Ro  Rossby number 
St  Strouhal number 
T  Thrust force 
U_{P}  Flow velocity component perpendicular to the rotor disk plane 
U_{T}  Flow velocity tangent to the rotor disk plane 
V  Wind speed 
W  Flow speed relative to a blade section 
β  Blade pitch 
ϵ  Strain 
θ  Sectional pitch angle 
κ  Reduced frequency 
λ  Tipspeed ratio 
λ_{d}  Design TSR 
μ  Fluid dynamic viscosity 
ν  Poisson coefficient 
ρ  Air density 
ρ_{m}  Material density 
ρ_{P}  Power density 
σ  Stress 
τ  Nondimensional time 
ω  Natural frequency 
Γ  Circulation 
Δθ  Total blade twist from root to tip 
Σ  Rotor solidity 
Φ  Rotor uptilt angle 
Ξ  Rotor cone angle 
Ω  Rotor angular velocity 
(⋅)_{a}  Pertaining to the aerodynamic design 
(⋅)_{s}  Pertaining to the structural design 
(⋅)_{M}  Scaled system 
(⋅)_{P}  Fullscale physical system 
$\dot{(\cdot )}$  Derivative with respect to time, i.e., $\mathrm{d}\cdot /\mathrm{d}t$ 
$(\cdot {)}^{\prime}$  Derivative with respect to nondimensional time, i.e., $\mathrm{d}\cdot /\mathrm{d}\mathit{\tau}$ 
$\stackrel{\mathrm{\u0303}}{(\cdot )}$  Nondimensional quantity 
$\widehat{(\cdot )}$  Tobematched scaled quantity 
AEP  Annual energy production 
BEM  Blade element momentum 
Bx  Biaxial 
CBRM  Combined bending root moment 
CFD  Computational fluid dynamics 
CFRP  Carbonfiberreinforced plastic 
DEL  Damage equivalent load 
DLC  Design load case 
EBRM  Edgewise bending root moment 
FBRM  Flapwise bending root moment 
GFRP  Glassfiberreinforced plastic 
LD  Low density 
LE  Leading edge 
MFTD  Maximum flapwise tip displacement 
METD  Maximum edgewise tip displacement 
PID  Proportional integral derivative 
PMMA  Polymethyl methacrylate 
POM  Polyoxymethylene 
PP  Polypropylene 
SQP  Sequential quadratic programming 
ThS  Thrust at main shaft 
TSR  Tipspeed ratio 
TE  Trailing edge 
Tx  Triaxial 
Ux  Uniaxial 
The data used for the present analysis can be obtained by contacting the authors.
HC modified the CpMax
code to support the scaled matching optimization, designed the subscale models, performed the simulations, and analyzed the results; CLB devised the original idea of this research, performed the theoretical scaling analysis, formulated the matching optimization problem, and supervised the work; and PB collaborated in the modification of the software, the design of the subscale models, and the conduction of the numerical simulations. HC and CLB wrote the manuscript. All authors provided important input to this research work through discussions and feedback and by improving the manuscript.
The authors declare that they have no conflict of interest.
The authors would like to thank Chengyu Wang and Daniel J. Barreiro of the Technical University of Munich for the computation of the airfoil polars using CFD for multiple Reynolds numbers. Additionally, credit goes to Eric Loth of the University of Virginia for having introduced the authors to the zooming approach and to Filippo Campagnolo of the Technical University of Munich for fruitful discussions and support. This work was authored in part by the National Renewable Energy Laboratory, operated by the Alliance for Sustainable Energy, LLC, for the US Department of Energy (DOE) under contract no. DEAC3608GO28308. Funding was provided by the US Department of Energy Office of Energy Efficiency and Renewable Energy Wind Energy Technologies Office. The views expressed in the article do not necessarily represent the views of the DOE or the US Government. The US Government retains and the publisher, by accepting the article for publication, acknowledges that the US Government retains a nonexclusive, paidup, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for US Government purposes.
This research has been supported by the BMWi through the WINSENT project (grant no. 0324129F) and the DOE (grant no. DEAC3608GO28308).
This paper was edited by Katherine Dykes and reviewed by two anonymous referees.
Anderson, B., Branlard, E., Vijayakumar, G., and Johnson, N.: Investigation of the nacelle blockage effect for downwind wind turbines, J. Phys. Conf. Ser., 1618, 062062, https://doi.org/10.1088/17426596/1618/6/062062, 2020. a
ANSYS Fluent: https://www.ansys.com/products/fluids/ansysfluent (last access: 18 December 2019), 2019. a, b
Armitt, J. and Counihan, J.: The simulation of the atmospheric boundary layer in a wind tunnel, J. Atmos. Environ., 2, 49–61, https://doi.org/10.1016/00046981(68)90019X, 1968. a
Azcona, J., Lemmer, F., Matha, D., Amann, F., Bottasso, C. L., Montinari, P., Chassapoyannis, P., Diakakis, K., Spyros, V., Pereira, R., Bredmose, H., Mikkelsen, R., Laugesen, R., and Hansen, A. M.: INNWIND. EU Deliverable D4.24: Results of wave tank tests, http://www.innwind.eu/publications/deliverablereports (last access: 18 December 2019), 2016. a, b
Bak, C., Zahle, F., Bitsche, R., Kim, T., Yde, A., Natarajan, A., and Hansen, M. H.: INNWIND. EU Deliverable D1.21: Reference Wind Turbine Report, http://www.innwind.eu/publications/deliverablereports (last access: 18 December 2019), 2013. a
Barlow, J. B., Rae, W. H., and Pope, A.: Lowspeed wind tunnel testing, 3rd Edn., Wiley, Hoboken, New Jersey, USA, 1999. a
Berg, J., Bryant, J., LeBlanc, B., Maniaci, D., Naughton, B., Paquette, J., Resor, B., and White, J.: Scaled Wind Farm Technology Facility Overview, in: 32nd ASME Wind Energy Symposium, AIAA SchiTech Forum, 13–17 January 2014, National Harbor, Maryland, https://doi.org/10.2514/6.20141088, 2014. a, b
Bisplinghoff, R. L. and Ashley, H., Principles of Aeroelasticity, Dover Publications, Mineola, New York, USA, 2002. a
Bortolotti, P., Bottasso, C. L., and Croce, A.: Combined preliminarydetailed design of wind turbines, Wind Energ. Sci., 1, 71–88, https://doi.org/10.5194/wes1712016, 2016. a, b
Bottasso, C. L. and Borri, M.: Integrating finite rotations, Comput. Method. Appl. Mech., 164, 307–331, https://doi.org/10.1016/S00457825(98)000310, 1998. a
Bottasso, C. L. and Campagnolo, F.: Wind tunnel testing of wind turbines and farms, Handbook of Wind Energy Aerodynamics, edited by: Stoevesandt, B., Schepers, G., Fuglsang, P., Sun, Y., Springer Nature, Cham, https://doi.org/10.1007/9783030054557_541, 2021. a, b, c
Bottasso, C. L., Campagnolo, F., and Croce, A.: Multidisciplinary constrained optimization of wind turbines, Multibody Syst. Dyn., 27, 21–53, https://doi.org/10.1007/s110440119271x, 2012. a
Bottasso, C. L., Campagnolo, F., and Petrovic, V.: Wind tunnel testing of scaled wind turbine models: Beyond aerodynamics, J. Wind Eng. Ind. Aerodyn., 127, 11–28, https://doi.org/10.1016/j.jweia.2014.01.009, 2014. a, b, c, d, e
Bottasso, C. L., Bortolotti, P., Croce, A., and Gualdoni, F.: Integrated aerostructural optimization of wind turbines, Multibody Syst. Dyn., 4, 317–344, https://doi.org/10.1007/s1104401594881, 2016. a, b
Bottasso, C. L., Cacciola, S., and Schreiber, J.: Local wind speed estimation, with application to wake impingement detection, Renew. Energ., 116, 155168, https://doi.org/10.1016/j.renene.2017.09.044, 2017. a
Brondsted, P., Lilholt, H., and Lystrup, A.: Composite Materials For Wind Power Turbine Blades, Annu. Rev. Mater. Res., 35, 505–538, https://doi.org/10.1146/annurev.matsci.35.100303.110641, 2005. a
Buckingham, E.: On Physically Similar Systems, Illustrations of the Use of Dimensional Equations, Phys. Rev., 4, 345–376, https://doi.org/10.1103/PhysRev.4.345, 1914. a
Burton, T., Jenkins, N., Sharpe, D., and Bossanyi, E.: Wind energy handbook, John Wiley & Sons, West Sussex, UK, 2001. a
Busan, R.: Flutter Model Technology, WLTR973074, WrightPatterson Air Force Base, OH, USA, 1998. a, b, c
Cambridge University Engineering Department: Materials Data Book, http://wwwmdp.eng.cam.ac.uk/web/library/enginfo/cueddatabooks/materials.pdf (last access: 18 December 2019), 2003. a, b, c
Campagnolo, F., Bottasso, C. L., and Bettini, P.: Design, manufacturing and characterization of aeroelastically scaled wind turbine blades for testing active and passive load alleviation techniques within a ABL wind tunnel, J. Phys. Conf. Ser., 524, 012061, https://doi.org/10.1088/17426596/524/1/012061, 2014. a, b, c
Chamorro, L. P., Arndt, R. E. A., and Sotiropoulos, F.: Reynolds number dependence of turbulence statistics in the wake of wind turbines, Wind Energy, 15, 733–742, https://doi.org/10.1002/we.501, 2012. a
Counihan, J.: An improved method of simulating an atmospheric boundary layer in a wind tunnel, Atmos. Environ., s3, 197–200, https://doi.org/10.1016/00046981(69)900080, 1969. a
Drela, M.: Xfoil 6.99 Documentation, http://web.mit.edu/drela/Public/web/xfoil/ (last access: 18 December 2019), 2017. a
Dowler, J. L. and Schmitz, S.: A solutionbased stall delay model for horizontalaxis wind turbines, Wind Energy, 18, 17931813, https://doi.org/10.1002/we.1791, 2015. a
Eggleston, D. M. and Stoddard, F. S.: Wind Turbine Engineering Design, Van Nostrand Reinhold, New York, NY, USA, 1987. a
Frederik, J. A., Weber, R., Cacciola, S., Campagnolo, F., Croce, A., Bottasso, C., and van Wingerden, J.W.: Periodic dynamic induction control of wind farms: proving the potential in simulations and wind tunnel experiments, Wind Energ. Sci., 5, 245–257, https://doi.org/10.5194/wes52452020, 2020. a
General Electric: GE Renewable Energy unveils the first HaliadeX 12 MW, the world's most powerful offshore wind turbine, Press release, available at: https://www.ge.com/news/pressreleases/gerenewableenergyunveilsfirsthaliadex12mwworlds (last access: 26 April 2021), 2019. a
Giavotto, V., Borri, M., Mantegazza, P. and Ghiringhelli, G.: Anisotropic beam theory and applications, Comput. Struct., 16, 403–413, https://doi.org/10.1016/00457949(83)901797, 1983. a, b
Hansen, M. H., Gaunaa, M., and Madsen, H. A.: A BeddoesLeishman type dynamic stall model in statespace and indicial formulations, Technical University of Denmark, Riso, Denmark, https://backend.orbit.dtu.dk/ws/portalfiles/portal/7711084/ris_r_1354.pdf (last access: 18 December 2019), 2004. a
Hideharu, M.: Realization of a largescale turbulence field in a small wind tunnel, Fluid Dyn. Res., 8, 1–4, https://doi.org/10.1016/01695983(91)90030M, 1991. a
International Electrotechnical Commission: International Electrotechnical Commission, IEC 614001 Edn. 3: Wind turbines – Part 1: Design requirements, IEC, Geneva, Switzerland, 2005. a
International Renewable Energy Agency: Future of wind: Deployment, investment, technology, grid integration and socioeconomic aspects (A Global Energy Transformation paper), Abu Dhabi, 2019. a
Jiménez, Á., Crespo, A., and Migoya E.: Application of a LES technique to characterize the wake deflection of a wind turbine in yaw, Wind Energy, 13, 559–572, https://doi.org/10.1002/we.380, 2010. a
Jones, R. M.: Buckling of Bars, Plates, and Shells, Bull Ridge Publishing, Virginia, 2006. a
Jonkman, J.: TurbSim User's Guide, NREL Report TP50036970, NREL, Golden, CO, USA, https://doi.org/10.2172/15020326, 2009. a, b
Kress, C., Chokani, N., and Abhari, R. S.: Downwind wind turbine yaw stability and performance, Renew. Energ., 83, 1157–1165, https://doi.org/10.1016/j.renene.2015.05.040, 2015. a
Loth, E., Kaminski, M., Qin, C., Fingersh, L. J., and Griffith, D. T.: GravoAeroelastic Scaling for ExtremeScale Wind Turbines, in: 35th AIAA Applied Aerodynamics Conference, AIAA AVIATION Forum, Denver, CO, https://doi.org/10.2514/6.20174215, 2017. a, b, c
McAuliffe, B., Larose, G.: Reynoldsnumber and surfacemodeling sensitivities for experimental simulation of flow over complex topography, J. Wind Eng. Ind. Aerod., 104–106, 603–613, https://doi.org/10.1016/j.jweia.2012.03.016, 2012. a
Manwell, J. F., McGowan, J. G., and Rogers, A. L.: Wind energy explained: theory, design and application, Second Edition, John Wiley & Sons Publication, West Sussex, United Kingdom, 2009. a, b, c
Murray, R. E., Jenne, S., Snowberg, D., Berry, D., and Cousins, D.: TechnoEconomic Analysis of a MegawattScale Thermoplastic Resin Wind Turbine Blade, Renew. Energ., 131, 111–119, https://doi.org/10.1016/j.renene.2018.07.032, 2018. a
Mydlarski, L.: A turbulent quarter century of active grids: from Makita (1991) to the present, Fluid Dyn. Res., 49, 061401, https://doi.org/10.1088/18737005/aa7786, 2017. a
Pitt, D. M. and Peters, D. A.: Theoretical prediction of dynamicinflow derivatives, Vertica, 5, 21–34, 1981. a
Resor, B. R. and Maniaci, D. C.: Definition of the National Rotor Testbed: An Aeroelastically Relevant ResearchScale Wind Turbine Rotor, in: 32nd ASME Wind Energy Symposium, AIAA SciTech Forum, 13–17 January 2014, National Harbor, MD, https://doi.org/10.2514/6.20140357, 2014. a, b
Ricciardi A. P., Canfield, R., Patil, M. J., and Lindsley, N.: Nonlinear aeroelastic scaledmodel design, J. Aircraft, 53, 20–32, https://doi.org/10.2514/1.C033171, 2016. a, b
Schmitz, S.: Aerodynamics of Wind Turbines, A Physical Basis for Analysis and Design, John Wiley & Sons Ltd, Hoboken, NJ, USA, 2020. a
Selig, M., Guglielmo, J., Broeren, A., and Giguère, P.: Summary of LowSpeed Airfoil Data, SoarTech Publications, Virginia, 1995. a, b
Siemens Gamesa: Powered by change: Siemens Gamesa launches 14 MW offshore direct drive turbine with 222meter rotor, Press release, available at: https://www.siemensgamesa.com/newsroom/2020/05/200519siemensgamesaturbine14222dd (last access: 26 April 2021), 2020. a
Sieros, G., Chaviaropoulos, P., Sorensen, J. D., Bulder, B. H., and Jamieson, P.: Upscaling wind turbines: theoretical and practical aspects and their impact on the cost of energy, Wind Energy, 15, 3–17, https://doi.org/10.1002/we.527, 2012. a
Simms, D., Schreck, S., Hand, M., and Fingersh, L. J.: NREL Unsteady Aerodynamics Experiment in the NASAAmes Wind Tunnel: A Comparison of Predictions to Measurements, available at: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.452.974&rep=rep1&type=pdf (last access: 26 April 2021), 2001. a
Snel, H., Schepers, J. G., and Siccama, N. B.: Mexico Project: The Database and Results of Data Processing and Interpretation, in: 47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition, Orlando, FL, https://doi.org/10.2514/6.20091217, 2009. a
3M Adhesives Division: ScotchWeld Structural Adhesive Film AF 32, Technical Data Issue No. 3, available at: https://multimedia.3m.com/mws/media/241415O/3mtmscotchweldtmstructuraladhesivefilmaf32.pdf (last access: 18 December 2019), 2000. a
Vermeer, L. J., Sorensen, J. N., and Crespo, A.: Wind Turbine Wake Aerodynamic, Prog. Aerosp. Sci., 39, 467–510, https://doi.org/10.1016/S03760421(03)000782, 2003. a, b
Vestas General Specification: Vestas V27–225 kW, 50 Hz Windturbine with Tubular/Lattice Tower, Version 1.2.0.24, 1994. a
Wan, Z. and Cesnik, C. E. S.: Geometrically nonlinear aeroelastic scaling for very flexible aircraft, AIAA J., 52, 2251–2260, https://doi.org/10.2514/1.J052855, 2014. a
Wang, C., Campagnolo, F., Canet, H., Barreiro, D. J., and Bottasso, C. L.: How realistic are the wakes of scaled wind turbine models?, Wind Energ. Sci. Discuss. [preprint], https://doi.org/10.5194/wes2020115, in review, 2020. a, b, c, d
ZSW – Zentrum für Solarenergie und WasserstoffForschung BadenWürttemberg: New WindForS project: Wind Energy Research in the Swabian Alps, available at: https://www.zswbw.de/en/newsroom/news/newsdetail/news/detail/News/newwindforsprojectwindenergyresearchintheswabianalps.html (last access: 18 December 2019), 2016. a, b
ZSW – Zentrum für Solarenergie und WasserstoffForschung BadenWürttemberg: ZSW and S&G Engineering Join Forces to Set Up Wind Power FieldTest Site, available at: https://www.zswbw.de/en/newsroom/news/newsdetail/news/detail/News/zswandsgengineeringjoinforcestosetupwindpowerfieldtestsite.html (last access: 18 December 2019), 2017. a