Journal cover
Journal topic
**Wind Energy Science**
The interactive open-access journal of the European Academy of Wind Energy

Journal topic

**Brief communication**
29 Jan 2019

**Brief communication** | 29 Jan 2019

Brief communication: Wind inflow observation from load harmonics – wind tunnel validation of the rotationally symmetric formulation

^{1}Wind Energy Institute, Technische Universität München, Garching bei München 85748, Germany^{2}Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano, Milano 20156, Italy

^{1}Wind Energy Institute, Technische Universität München, Garching bei München 85748, Germany^{2}Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano, Milano 20156, Italy

**Correspondence**: Carlo L. Bottasso (carlo.bottasso@tum.de)

**Correspondence**: Carlo L. Bottasso (carlo.bottasso@tum.de)

Abstract

Back to toptop
The present paper further develops and experimentally validates the
previously published idea of estimating the wind inflow at a turbine rotor
disk from the machine response. A linear model is formulated that relates one
per revolution (1P) harmonics of the in- and out-of-plane blade root bending
moments to four wind parameters, representing vertical and horizontal shears
and misalignment angles. Improving on this concept, the present work exploits
the rotationally symmetric behavior of the rotor in the formulation of the
load-wind model. In a nutshell, this means that the effects on the loads of
the vertical shear and misalignment are the same as those of the horizontal
quantities, simply shifted by *π*∕2. This results in a simpler
identification of the model, which needs a reduced set of observations. The
performance of the proposed method is first tested in a simulation
environment and then validated with an experimental data set obtained with an
aeroelastically scaled turbine model in a boundary layer wind tunnel.

Download & links

How to cite

Back to top
top
How to cite.

Bertelè, M., Bottasso, C. L., and Cacciola, S.: Brief communication: Wind inflow observation from load harmonics – wind tunnel validation of the rotationally symmetric formulation, Wind Energ. Sci., 4, 89–97, https://doi.org/10.5194/wes-4-89-2019, 2019.

1 Introduction

Back to toptop
The ability to control a system is often intimately linked to the awareness of the surrounding environment. For a wind turbine, the environment is represented by the wind inflow, which is characterized by speed, direction, shears, veer, turbulence intensity, presence of impinging wakes, etc. Such parameters have a profound effect on the response of a single wind turbine as well as on clusters of interacting machines within a power plant. Better awareness of the wind environment can be translated into better turbine-level and plant-level operation and control.

The current standard equipment mounted on board wind turbines for the measurement of the wind inflow is composed of one or more anemometers and wind vanes, typically located at hub height, either on the nacelle or on the spinner. Even when properly calibrated, all such devices suffer from one inherent unavoidable limitation: they provide measurements at the single point in space where they are located. As such, they are necessarily blind to all wind characteristics that imply wind variations across the rotor disk. Alternative sensors are represented by lidars, which are, however, not yet routinely installed on board wind turbines because of cost, availability, reliability, effects due to weather conditions and lifetime issues. In this sense, current wind turbines have only a very limited awareness of the environment in which they operate.

The concept of the “rotor as a sensor” was developed to address the limitations of current wind measurement devices. The idea is conceptually very simple: changes in the wind inflow produce changes in the wind turbine response. If the wind-response map is known, one can then measure the response (for example, in the form of loads and/or accelerations) and estimate the inflow by inverting the map.

Various formulations have been proposed for this concept (Bottasso et al., 2010; Bottasso and Riboldi, 2014; Simley and Pao, 2014; Bottasso and Riboldi, 2015). In this paper we improve on the work described by Cacciola et al. (2016a) and Bertelè et al. (2017, 2018). The approach parameterizes the inflow in terms of four quantities: vertical and horizontal shears and misalignment angles. The wind-response map relates these four wind states to the 1P in- and out-of-plane blade root bending moments. Both linear and quadratic maps were considered in Bertelè et al. (2017), with a marginally better accuracy for the latter. System identification was used to find the model coefficients from simulations performed with an aeroservoelastic model in a variety of wind conditions, spanning the range of interest of the four wind states. Results indicate a better accuracy of the shears than the angles, although the latter are still well captured in their mean values.

Despite the more than promising results reported in Bertelè et al. (2017), the identification of the model relating wind states to load harmonics can be cumbersome. In fact, a data set is required that covers a desired range of the four wind states. While this is not a major issue in a simulation environment where one can generate all desired wind conditions, an identification based on field test data might not be easy or even possible. In fact, some wind parameters might not change much at a given site, e.g. upflow angle and horizontal shear. This would clearly be a major hurdle, as a model only knows what is in the data used for training it.

To address this issue, the present work exploits the rotationally symmetric behavior of
the rotor. In fact, the effect caused by a horizontal shear on the rotor response is the
same as that caused by a vertical shear, only shifted by *π*∕2. Similarly, the effect
of a vertical upflow angle is the same of a horizontal yaw misalignment, again shifted by
*π*∕2. This means that one can collect data sets containing the desired changes in
vertical shears and yaw misalignments, and identify a model that is also capable of
representing the same range of horizontal shears and upflow angles.

The paper is organized as follows. Section 2 first introduces the wind parameterization and the wind-load map, and then uses the rotational symmetry of the rotor to eliminate some of the model coefficients from the identification problem unknowns. Section 3 compares the results of the new formulation to the original one first by simulations – conducted with an aeroservoelastic model – and then experimentally – using a scaled turbine in a wind tunnel. Finally, the work is closed by Sect. 4, where conclusions are drawn.

2 Formulation

Back to toptop
The wind inflow is characterized in terms of four so-called wind states, which are
defined as the vertical (upflow) and horizontal (yaw) misalignment angles *χ* and
*ϕ*, respectively, and the vertical and horizontal linear shears *κ*_{v}
and *κ*_{h}, respectively. These quantities should be regarded as
rotor-equivalent fits of the actual spatial distribution of the wind impinging on the
rotor disk at a certain instant of time.

The wind states are defined with respect to a nacelle-attached reference frame
(** x**,

$$\begin{array}{}\text{(1a)}& {\displaystyle}u(y,z)& {\displaystyle}=W(y,z)\mathrm{cos}\left(\mathit{\varphi}\right)\mathrm{cos}\left(\mathit{\chi}\right),\text{(1b)}& {\displaystyle}v(y,z)& {\displaystyle}=W(y,z)\mathrm{sin}\left(\mathit{\varphi}\right)\mathrm{cos}\left(\mathit{\chi}\right),\text{(1c)}& {\displaystyle}w(y,z)& {\displaystyle}=W(y,z)\mathrm{sin}\left(\mathit{\chi}\right),\end{array}$$

where *W*(*y*,*z*) is a linearly sheared wind field

$$\begin{array}{}\text{(2)}& W(y,z)={V}_{\mathrm{H}}\left(\mathrm{1}+{\displaystyle \frac{z}{R}}{\mathit{\kappa}}_{\mathrm{v}}+{\displaystyle \frac{y}{R}}{\mathit{\kappa}}_{\mathrm{h}}\right),\end{array}$$

*V*_{H} being the wind speed at hub height, and *R* the rotor radius.
According to this definition, the yaw misalignment and upflow angles are
positive when the wind blows from the left and the lower part of the rotor,
respectively, when looking upstream.

Notice that the formulation of Cacciola et al. (2016a) used a horizontal reference frame with respect to the terrain, while in the present case the frame is aligned with the rotor axis. Together with the assumed linearity of both shears, this is necessary in order to exploit the rotational symmetry of the rotor response. Hence, if the rotor is uptilted, one will have to transform the nacelle-frame wind components into a frame aligned with the ground if necessary.

Looking at Eq. (2), it appears that the
effect of the vertical shear *κ*_{v} on the velocity distribution
is the same of the one caused by the horizontal shear *κ*_{h},
when rotated by *π*∕2. On the other hand, looking at
Eq. (1a–c), the effect of
the angles *ϕ* and *χ* is more complex. To eliminate this problem, the
rotor-in-plane wind velocity components can be expressed in terms of the new
variables

$$\begin{array}{}\text{(3a)}& {\displaystyle}\stackrel{\mathrm{\u0303}}{v}& {\displaystyle}={\displaystyle \frac{v(\mathrm{0},\mathrm{0})}{{V}_{\mathrm{H}}}}=\mathrm{sin}\left(\mathit{\varphi}\right)\mathrm{cos}\left(\mathit{\chi}\right),\text{(3b)}& {\displaystyle}\stackrel{\mathrm{\u0303}}{w}& {\displaystyle}={\displaystyle \frac{w(\mathrm{0},\mathrm{0})}{{V}_{\mathrm{H}}}}=\mathrm{sin}\left(\mathit{\chi}\right),\end{array}$$

which, respectively, represent the nondimensional horizontal and vertical wind cross flows at the hub. This change of variables results in

$$\begin{array}{}\text{(4a)}& {\displaystyle}u(y,z)& {\displaystyle}=W(y,z)\sqrt{\mathrm{1}-{\stackrel{\mathrm{\u0303}}{v}}^{\mathrm{2}}-{\stackrel{\mathrm{\u0303}}{w}}^{\mathrm{2}}},\text{(4b)}& {\displaystyle}v(y,z)& {\displaystyle}=W(y,z)\stackrel{\mathrm{\u0303}}{v},\text{(4c)}& {\displaystyle}w(y,z)& {\displaystyle}=W(y,z)\stackrel{\mathrm{\u0303}}{w}.\end{array}$$

With this reformulation, the effect of $\stackrel{\mathrm{\u0303}}{v}$ on *v* is the same as the
effect of $\stackrel{\mathrm{\u0303}}{w}$ on *w*, when rotated by *π*∕2. Given $\stackrel{\mathrm{\u0303}}{v}$ and
$\stackrel{\mathrm{\u0303}}{w}$, the misalignment angle *ϕ* and upflow *χ* can be readily
recovered by inverting their respective
definitions (Eq. 3a–b):

$$\begin{array}{}\text{(5a)}& {\displaystyle}\mathit{\chi}=& {\displaystyle}\mathrm{arcsin}\left(\stackrel{\mathrm{\u0303}}{w}\right),\text{(5b)}& {\displaystyle}\mathit{\varphi}=& {\displaystyle}\mathrm{arcsin}(\stackrel{\mathrm{\u0303}}{v}/\mathrm{cos}\mathit{\chi}),\end{array}$$

although for small angles the difference between the two sets of variables will be negligible.

In this work, the linear model of Cacciola et al. (2016a) and Bertelè et al. (2017) is used to relate inflow conditions and machine response. The model writes

$$\begin{array}{ll}{\displaystyle}\mathit{m}& {\displaystyle}=\phantom{\rule{0.125em}{0ex}}\mathbf{F}(V,\mathit{\varrho})\mathit{\theta}+{\mathit{m}}_{\mathrm{0}}(V,\mathit{\varrho})\\ \text{(6)}& {\displaystyle}& {\displaystyle}=\left[\mathbf{F}(V,\mathit{\varrho})\phantom{\rule{0.33em}{0ex}}{\mathit{m}}_{\mathrm{0}}(V,\mathit{\varrho})\right]\left[\begin{array}{c}\mathit{\theta}\\ \mathrm{1}\end{array}\right]=\mathbf{T}\phantom{\rule{0.33em}{0ex}}\stackrel{\mathrm{\u203e}}{\mathit{\theta}},\end{array}$$

where ** m** is the load vector, $\mathit{\theta}={\left\{\stackrel{\mathrm{\u0303}}{v}\phantom{\rule{0.33em}{0ex}}{\mathit{\kappa}}_{\mathrm{v}}\phantom{\rule{0.33em}{0ex}}\stackrel{\mathrm{\u0303}}{w}\phantom{\rule{0.33em}{0ex}}{\mathit{\kappa}}_{\mathrm{h}}\right\}}^{\mathrm{T}}$ is the wind
state vector, while

$$\begin{array}{}\text{(7)}& \mathit{m}={\left\{{m}_{\mathrm{1}\mathrm{c}}^{\mathrm{OP}},\phantom{\rule{0.33em}{0ex}}{m}_{\mathrm{1}\mathrm{s}}^{\mathrm{OP}},\phantom{\rule{0.33em}{0ex}}{m}_{\mathrm{1}\mathrm{c}}^{\mathrm{IP}},\phantom{\rule{0.33em}{0ex}}{m}_{\mathrm{1}\mathrm{s}}^{\mathrm{IP}}\right\}}^{\mathrm{T}},\end{array}$$

where *m* indicates the blade bending moment, subscripts (⋅)_{1s} and
(⋅)_{1c}, respectively, indicate sine and cosine
harmonics, while superscripts (⋅)^{OP} and (⋅)^{IP},
respectively, out- and in-plane components. The load harmonics are readily computed via
the Coleman and Feingold transformation (Coleman and Feingold, 1958) once three
measured blade loads are available. For simplicity and brevity, the present paper only
considers a linear wind-response map. However, nonlinearities in the map can be readily
included, as shown by Bertelè et al. (2017).

To identify the model coefficients **T**, one should collect a rich
enough data set for which both wind states ** θ** and associated
blade loads

$$\begin{array}{}\text{(8)}& \mathbf{M}=\mathbf{T}\mathbf{\Theta}.\end{array}$$

Finally, the model coefficients are readily identified as

$$\begin{array}{}\text{(9)}& \mathbf{T}=\mathbf{M}{\mathbf{\Theta}}^{\mathrm{T}}(\mathbf{\Theta}{\mathbf{\Theta}}^{\mathrm{T}}{)}^{-\mathrm{1}}.\end{array}$$

The invertibility of the system is discussed in Bertelè et al. (2017).

Once the model expressed by Eq. (6) has been identified, it can be
used to express the dependency of given measured loads *m*_{M} on the wind
states,

$$\begin{array}{}\text{(10)}& {\mathit{m}}_{\mathrm{M}}=\mathbf{F}\mathit{\theta}+{\mathit{m}}_{\mathrm{0}}+\mathit{r},\end{array}$$

where ** r** is the measurement error, and the dependency on

$$\begin{array}{}\text{(11)}& {\mathit{\theta}}_{\mathrm{E}}={\left({\mathbf{F}}^{\mathrm{T}}{\mathbf{R}}^{-\mathrm{1}}\mathbf{F}\right)}^{-\mathrm{1}}{\mathbf{F}}^{\mathrm{T}}{\mathbf{R}}^{-\mathrm{1}}({\mathit{m}}_{\mathrm{M}}-{\mathit{m}}_{\mathrm{0}}),\end{array}$$

where **R**=**E**[*r**r*^{T}] is the covariance weighting
matrix. Given *θ*_{E}, the misalignment and upflow angles can be
recovered by using Eq. (5).

By considering the rotational symmetry of the rotor, the number of unknown coefficients
in **F** can be reduced. Indeed, a vertical shear will cause the same response of
an equivalent horizontal shear, simply shifted by an azimuthal delay of *π*∕2. The same
consideration holds for the vertical and horizontal cross flows. This rotational symmetry
is reflected in the derivatives of the loads with respect to the wind states, i.e., in
the coefficients of matrix **F**. By a rotation of *π*∕2, the load component
*m*_{1c} becomes *m*_{1s}, while the load component *m*_{1s}
becomes −*m*_{1c}. As a result, the following conditions apply between pairs of
model coefficients:

$$\begin{array}{}\text{(12a)}& {\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{c}}}{\partial \stackrel{\mathrm{\u0303}}{v}}}& {\displaystyle}={\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{s}}}{\partial \stackrel{\mathrm{\u0303}}{w}}},\text{(12b)}& {\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{s}}}{\partial \stackrel{\mathrm{\u0303}}{v}}}& {\displaystyle}=-{\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{c}}}{\partial \stackrel{\mathrm{\u0303}}{w}}},\text{(12c)}& {\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{c}}}{\partial {\mathit{\kappa}}_{\mathrm{h}}}}& {\displaystyle}={\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{s}}}{\partial {\mathit{\kappa}}_{\mathrm{v}}}},\text{(12d)}& {\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{s}}}{\partial {\mathit{\kappa}}_{\mathrm{h}}}}& {\displaystyle}=-{\displaystyle \frac{\partial {m}_{\mathrm{1}\mathrm{c}}}{\partial {\mathit{\kappa}}_{\mathrm{v}}}}.\end{array}$$

These conditions apply to both the out- and the in-plane components.

The term *m*_{0} in Eq. (6) represents the effects of gravity on
the loads (Bertelè et al., 2017). Since this term is nonsymmetric, no reduction of these
coefficients is possible in this case.

The advantage of this approach is not only in the reduced number of unknown model coefficients, but, most importantly, in the reduced datapoints necessary for identification. In fact, by eliminating the coefficients of horizontal shear and upflow angle, one can use tests in which only yaw misalignment angle and vertical shear are changing. Therefore, since the model is linear and depends on two parameters, a minimum of only three operating conditions is required for identification.

3 Results

Back to toptop
The proposed method was first tested by numerical simulations, using the model of a
horizontal-axis three-bladed 3 MW wind turbine. The machine has a rotor diameter of
93 m; a hub height of 80 m; 4.5^{∘} of nacelle uptilt;
and cut-in, rated and cut-out speeds equal to 3, 12.5 and 25 m s^{−1}, respectively.
A transition region *I**I* 1∕2 connects the partial- and
full-load regimes, extending between 9 and 12.5 m s^{−1}. The machine response was
simulated by the aeroservoelastic multibody software `Cp-Lambda`

(Bauchau et al., 2003; Bottasso and Croce, 2006), which is based on a
geometrically exact finite element formulation. The model includes flexible blades, tower
and drive train, and compliant foundations. The collective pitch and torque controller is
implemented according to Riboldi (2012) and
Bottasso et al. (2012), while generator and pitch actuators are
modeled as first- and second-order dynamical systems, respectively. The aerodynamic rotor
model is based on blade element momentum theory (BEM), augmented by classical tip and
root losses, unsteady aerodynamics and dynamic stall models. Turbulent wind time
histories were generated with the `TurbSim`

code (Jonkman and Kilcher, 2012) in
accordance with the Kaimal model, at the nodes of a square grid overlapping the rotor
disk. “Ground truth” values of the wind states – to be used for assessing the quality
of observed quantities – were obtained by fitting the instantaneous wind field at the
grid nodes to the rotor swept area.

Turbulent simulations were run for a duration of 10 min, according to standard practice. The 1P harmonics were computed by the Coleman and Feingold transformation (Coleman and Feingold, 1958), using in- and out-of-plane bending moment components measured by strain gauges placed at the root of each blade. The resulting signal was finally cleaned with a low-pass filter; on-line adaption of the filter parameters was used to account for changes in rotational speed due to turbulent wind fluctuations.

Two observation models were identified. The first is the linear formulation of Bertelè et al. (2017), which does not exploit the rotational symmetry of the rotor, while the second is the linear rotationally symmetric formulation of the present paper. In the first case, the model was identified from nonturbulent wind cases corresponding to all combinations of the following wind parameters:

$$\begin{array}{}\text{(13a)}& {\displaystyle}& {\displaystyle}\mathit{\varphi}=\left[\mathrm{0}\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}\mathrm{16}\right]{}^{\circ},\text{(13b)}& {\displaystyle}& {\displaystyle}{\mathit{\kappa}}_{\mathrm{v}}=\left[\mathrm{0.06}\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}\mathrm{0.18}\right],\text{(13c)}& {\displaystyle}& {\displaystyle}\mathit{\chi}=\left[\mathrm{4.5}\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}\mathrm{16.5}\right]{}^{\circ},\text{(13d)}& {\displaystyle}& {\displaystyle}{\mathit{\kappa}}_{\mathrm{h}}=[\mathrm{0}\phantom{\rule{0.33em}{0ex}}-\mathrm{0.1}].\end{array}$$

A separate identification was performed for each
wind speed, considering the values *V*=[3 4 5 6 7 8 9 11 15 19] m s^{−1}. A second model was obtained by exploiting symmetry and
linearity. Accordingly, the identification set was reduced to the following wind
parameter combinations:

$$\begin{array}{}\text{(14a)}& {\displaystyle}& {\displaystyle}\mathit{\varphi}=\left[\mathrm{0}\phantom{\rule{0.33em}{0ex}}\mathrm{16}\right]{}^{\circ},\text{(14b)}& {\displaystyle}& {\displaystyle}{\mathit{\kappa}}_{\mathrm{v}}=\left[\mathrm{0.06}\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}\mathrm{0.18}\right],\text{(14c)}& {\displaystyle}& {\displaystyle}\mathit{\chi}=\mathrm{4.5}{}^{\circ},\text{(14d)}& {\displaystyle}& {\displaystyle}{\mathit{\kappa}}_{\mathrm{h}}=\mathrm{0},\end{array}$$

therefore assuming both upflow *χ* and horizontal shear
*κ*_{h} to be constant. Notice that the upflow angle is set to
4.5^{∘}, which corresponds to the rotor uptilt.

The two models were then tested and compared in turbulent wind conditions. Three
different combinations of inflow angles and shears (not included in the identification
set) were considered, each using four different turbulent realizations, for a total of 12
tests performed at each given wind speed and turbulence intensity (TI).
Figures 2 and 3 show, respectively, the mean (over 10 min and
over all turbulent seeds) absolute error *ϵ* and standard deviation *σ* as
functions of wind speed, for two different levels of TI, equal to 5 % and 12 %.
The results of the reference full model are shown using solid lines, while the ones of
the rotationally symmetric formulation using dashed lines. The two formulations appear to
be characterized by a very similar performance. Actually, notwithstanding its reduced
identification set, the symmetric method obtains marginally better results. As expected,
TI has a negative effect on the quality of the estimates. In addition, as already noticed
in Bertelè et al. (2017), angle estimates appear to be less precise than shear
estimates. Nonetheless, for 12 % TI at 15 m s^{−1}, the yaw misalignment mean
error is about 2.5^{∘}. This appears to be a good result when compared to the
typical accuracy of nacelle-mounted anemometers.

Next, the proposed formulation was tested using an aeroelastically scaled wind turbine
operated in a boundary layer wind tunnel. The scaled model represents a three-bladed
horizontal-axis wind turbine with a hub height of about 1.8 m, a rotor diameter of 2 m
and a rated wind speed of 6 m s^{−1} (Bottasso et al., 2014). The turbine design
preserves the tip speed ratio, Lock number, and
placement of the lowest tower and rotor nondimensional frequencies of the reference
machine, resulting in a scaled model of realistic aeroelastic behavior
(Bottasso et al., 2014). Each of the flexible scaled blades is equipped with strain
gauges at the blade roots, which measure the flapwise and edgewise bending moments, while
an optical incremental encoder is used to measure the blade azimuthal position.

Tests were performed in the boundary layer test section of the wind tunnel of Politecnico di Milano (Bottasso et al., 2014). Two different boundary layer conditions, characterized by different mean vertical shears and TI levels, were obtained by the use of suitable turbulence generators at the chamber inlet and roughness elements placed on the floor. Such inflow conditions were then accurately mapped over the rotor swept area with triple hot-wire probes, providing a reference mean inflow that can be considered the “ground truth”. The lower turbulence condition was characterized by a TI of 3.8 % and a linear vertical shear of 0.03, while the higher turbulence case by a TI of 8.5 % and a linear vertical shear of 0.12.

For various wind speeds, several tests were performed for different combinations of yaw
misalignment, vertical shear and upflow angle as reported in Table 1.
Changes in mean vertical shear were obtained by changing the wind tunnel boundary layer
conditions. Changes in mean misalignment angle were realized by yawing the turbine model
with respect to the wind. To create different upflow angles, the wind turbine tower foot
was installed on a tiltable ramp. By changing the ramp angle, the turbine can be pitched
by ±6^{∘}. As the rotor has an uptilt angle of 6^{∘} with respect to the
tower, the use of the ramp allows one to obtain upflow angles between 0 and 12^{∘}.
Finally, the horizontal shear for all tests can be considered null, as the flow in the
wind tunnel is essentially uniform in the lateral direction.

A total number of 174 different conditions were tested. The entire set of experiments was
then divided into two subsets. The first one was used for identifying the observer model,
and it contains two combinations of vertical shear and misalignment angle per wind speed,
with an upflow of 6^{∘}; these test points are indicated with “×” symbols
in Table 1. The second subset was instead used for validating the
observer performance. This second subset contains all the other experiments, indicated
with “∘” symbols in Table 1. Notice that the second set of
experiments correspond to upflow angles of 0 and 12^{∘}, values that are not
contained in the identification set. This is possible thanks to the symmetry of the
rotor: the information contained in the identification set on the effect of the
misalignment angle is used to infer the effect of the upflow, although no operating
points at different upflows are used during training.

To validate the performance of the observer, the machine response during each test was
averaged over a time window of 180 s in order to estimate the corresponding mean inflow
parameters. The length of the time window is dictated in this case not only by the need
to average out turbulent fluctuations, but also by the dynamic characteristics of this
particular closed-return wind tunnel. Figure 4 shows an
excerpt of the results obtained at a wind speed of 7 m s^{−1}
(Fig. 4a), which corresponds to the beginning of the
full load region, and a speed of 5.5 m s^{−1}
(Fig. 4b), which corresponds to the end of the partial
load region. In each panel, the reference (true) wind parameter is shown on the *x* axis,
while the corresponding observed quantity is given on the *y* axis. It follows that an
ideal match would be represented by the bisector of the quadrant. The yaw misalignment
estimation (Fig. 4a) appears to be quite accurate and
has a maximum error of less than 1.3^{∘}. Better accuracy can be achieved for high
positive yaw angles; this is to be expected, since such conditions are included in the
identification set (cf. Table 1). Even the upflow estimation
(Fig. 4b) appears to be quite accurate, with a maximum
error of about 1.5^{∘}. Note that the accuracy in the upflow estimation validates
the assumption of rotational symmetry of the parameters, as no upflow changes were
present in the data set used for identifying the load-wind model (again, cf.
Table 1). Indeed, the model coefficients related to this parameter
were obtained using the symmetry conditions given by
Eq. (12a–d).

Finally, to better understand the performance of the observer, mean inflow parameters
were estimated and compared to the respective ground truth for each test not included in
the identification set. For each wind speed, such mean errors were averaged over the
number of tests and reported in Fig. 5. Here again, results
appear to be significantly accurate; in fact, for both turbulence levels, a maximum mean
error smaller than 1^{∘} is observed in the angle estimates, while the error in the
shear estimates is less than $\mathrm{6}\times {\mathrm{10}}^{-\mathrm{3}}$.

Comparing the experimental results with the numerical ones in the low TI cases (equal to 3.8 % and 5 %, respectively), one should notice that the mean estimation errors present the same range of accuracy, as one can appreciate by comparing Fig. 5 with Fig. 2. This can be considered an additional proof of the general applicability of the method, since these results were obtained with two different models applied to two very different machines, using numerical and experimental data sets.

4 Conclusions

Back to toptop
Following the work presented in Cacciola et al. (2016a) and Bertelè et al. (2017), this paper has further developed and experimentally validated a method to estimate the inflow at the rotor disk. Specifically, a linear model was formulated to estimate four wind parameters: the vertical and horizontal shears, and the vertical and horizontal wind misalignments. Improving on the previous publications, the rotationally symmetric behavior of the rotor was exploited in order to simplify the model identification procedure, by reducing the number of necessary measured operating conditions.

The performance of the proposed rotationally symmetric model was tested both in
simulation and with an aeroelastically scaled wind turbine model in a boundary layer wind
tunnel. Results indicate no significant difference in the accuracy of the new
rotationally symmetric formulation with respect to the original one, even if the number
of tests required for identification is significantly decreased. The expected mean error
in the angle estimation is less than 1 and 2.5^{∘} for low and high TI levels,
respectively. An even higher accuracy can be obtained for the estimation of shears.
Moreover, the experimental results are well in line with the ones obtained by numerical
simulations.

Data availability

Back to toptop
Data availability.

Data can be provided upon request. Please con-

tact the
corresponding author Carlo L. Bottasso

(carlo.bottasso@tum.de).

Appendix A: Nomenclature

Back to toptop
m |
Generic blade moment |

m |
Vector of moment harmonics |

R |
Rotor radius |

V |
Wind speed |

V |
Wind vector |

$\stackrel{\mathrm{\u0303}}{v}$ | Non-dimensional horizontal cross flow at the hub |

$\stackrel{\mathrm{\u0303}}{w}$ | Non-dimensional vertical cross flow at the hub |

ϱ |
Air density |

ϕ |
Yaw misalignment angle |

χ |
Upflow angle |

κ_{v} |
Vertical shear |

κ_{h} |
Horizontal shear |

ϵ |
Mean error |

σ |
Standard deviation |

θ |
Wind state vector |

(⋅)^{T} |
Transpose |

(⋅)_{E} |
Estimated quantity |

(⋅)^{OP} |
Out-of-plane quantity |

(⋅)^{IP} |
In-plane quantity |

(⋅)_{1c} |
1P cosine amplitude |

(⋅)_{1s} |
1P sine amplitude |

BEM | Blade element momentum |

Lidar | Light detection and ranging |

TI | Turbulence intensity |

1P | Once per revolution |

Author contributions

Back to toptop
Author contributions.

All authors equally contributed to this work.

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements

Back to toptop
Acknowledgements.

This work has been partially supported by the CL-Windcon project,
which receives funding from the European Union Horizon 2020 research and
innovation program under grant agreement No. 727477.

This work was supported by the German Research

Foundation (DFG) and the Technical
University of Munich (TUM) in the framework of the Open Access Publishing Program.

Edited by: Sandrine Aubrun

Reviewed by: two anonymous referees

References

Back to toptop
Bauchau, O. A., Bottasso, C. L., and Trainelli, L.: Robust integration schemes for flexible multibody systems, Comput. Method. Appl. M., 192, 395–420, https://doi.org/10.1016/S0045-7825(02)00519-4, 2003. a

Bertelè, M., Bottasso, C. L., Cacciola, S., Daher Adegas, F., and Delport, S.: Wind inflow observation from load harmonics, Wind Energ. Sci., 2, 615–640, https://doi.org/10.5194/wes-2-615-2017, 2017. a, b, c, d, e, f, g, h, i, j

Bertelè, M., Bottasso, C. L., and Cacciola, S.: Simultaneous estimation of wind shears and misalignments from rotor loads: formulation for IPC-controlled wind turbines, J. Phys. Conf. Ser., 1037, 032007, https://doi.org/10.1088/1742-6596/1037/3/032007, 2018. a

Bottasso, C. L. and Croce, A.: Cp-Lambda: user's manual, Technical Report, Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano, Milano, Italy, 2006. a

Bottasso, C. L. and Riboldi, C. E. D.: Estimation of wind misalignment and vertical shear from blade loads, Renew. Energ., 62, 293–302, https://doi.org/10.1016/j.renene.2013.07.021, 2014. a

Bottasso, C. L. and Riboldi, C. E. D.: Validation of a wind misalignment observer using field test data, Renew. Energ., 74, 298–306, https://doi.org/10.1016/j.renene.2014.07.048, 2015. a

Bottasso, C. L., Croce, A., and Riboldi, C. E. D.: Spatial estimation of wind states from the aeroelastic response of a wind turbine, The Science of Making Torque from Wind (TORQUE 2010) 28–30 June 2010, Heraklion, Crete, Greece, 2010. a

Bottasso, C. L., Croce, A., Nam, Y., and Riboldi, C. E. D.: Power curve tracking in the presence of a tip speed constraint, Renew. Energ., 40, 1–12, https://doi.org/10.1016/j.renene.2011.07.045, 2012. a

Bottasso, C. L., Campagnolo, F., and Petrovic̀: Wind tunnel testing of scaled wind turbine models: beyond aerodynamics, J. Wind. Eng. Ind. Aerod., 127, 11–28, https://doi.org/10.1016/j.jweia.2014.01.009, 2014. a, b, c

Bottasso, C. L., Cacciola, S., and Schreiber, J.: A wake detector for wind farm control, J. Phys. Conf. Ser., 625, 012007, https://doi.org/10.1088/1742-6596/625/1/012007, 2015.

Bottasso, C. L., Cacciola, S., and Schreiber, J.: Local wind speed estimation, with application to wake impingement detection, Renew. Energ., 116, 155–168, 2018.

Cacciola, S., Bertelè, M., and Bottasso, C. L.: Simultaneous observation of wind shears and misalignments from rotor loads, J. Phys. Conf. Ser., 753, 052002, https://doi.org/10.1088/1742-6596/753/5/052002, 2016a. a, b, c, d

Cacciola, S., Bertelè, M., Schreiber, J., and Bottasso, C. L.: Wake center position tracking using downstream wind turbine hub loads, J. Phys. Conf. Ser., 753, 032036-1–6, https://doi.org/10.1088/1742-6596/753/3/032036, 2016b.

Coleman, R. P. and Feingold, A. M.: Theory of self-excited mechanical oscillations of helicopter rotors with hinged blades, Technical Report, NACA TN 1351, 1958. a, b

Jonkman, B. J. and Kilcher, L.: TurbSim user's guide: version 1.06.00, NREL Technical report, 2012. a

Riboldi, C. E. D.: Advanced control laws for variable-speed wind turbines and supporting enabling technologies, PhD thesis, Politecnico di Milano, Milano, Italy, 2012. a

Simley, E. and Pao, L. Y.: Evaluation of a wind speed estimator for effective hub-height and shear components, Wind Energy, 19, 167–184, https://doi.org/10.1002/we.1817, 2014. a

Short summary

This paper describes a new formulation for estimating the wind
inflow at the rotor disk, based on measurements of the blade loads.
The new method improves on previous formulations by exploiting the
rotational symmetry of the problem. Experimental results obtained
with an aeroelastically scaled model in a boundary layer wind
tunnel are used for validating the proposed approach.

This paper describes a new formulation for estimating the wind
inflow at the rotor disk, based...

Wind Energy Science

The interactive open-access journal of the European Academy of Wind Energy