The calibration coefficients of two commercial anemometers equipped with different rotors were studied. The rotor cups had the same conical shape, while the size and distance to the rotation axis varied. The analysis was based on the 2-cup positions analytical model, derived using perturbation methods to include second-order effects such as pressure distribution along the rotating cups and friction. The comparison with the experimental data indicates a nonuniform distribution of aerodynamic forces on the rotating cups, with higher forces closer to the rotating axis. The 2-cup analytical model is proven to be accurate enough to study the effect of complex forces on cup anemometer performance.

1. Introduction

1.1. Wind Speed Anemometry as an Important Tool in Wind Energy Generation

The importance of accuracy in wind speed measurements is emphasized as the wind energy sector is highly concerned with both wind turbine performance control and wind energy forecast on the field [1, 2]. The aforementioned accuracy of the measurements directly affects wind energy production, as this production is proportional to the third power of the wind speed [3].

On the other hand, it can be said that today the wind energy sector represents the larger demand of anemometers in the market, despite the increasing use of anemometers in other industries/applications. In addition, it seems that the demand of accurate anemometers will remain strong, although the investment in the wind energy industry has decreased in the traditionally leaders of the sector (Germany, Spain, and Denmark), new players are now very active (China, USA, India, and Brazil) [4].

Finally, the cup anemometer is at present the standardized instrument included in the most relevant code of practice regarding wind turbine power performance measurements (IEC 61400-12-1) [5].

1.2. Cup Anemometer Aerodynamics

A cup anemometer can be studied from two different perspectives: as a meteorological instrument or as a body in autorotation. As a meteorological instrument, the cup anemometer has been studied for a long time, using different techniques and mathematical models, under different climatic conditions and focusing on particular aspects of their performance and response. In addition to these, it should also be said that some important research projects regarding cup anemometers have been carried out based on public funds [69]. Table 1 summarizes some of these aspects of cup anemometers, along with the authors of the corresponding research contributions (an extensive review of the available literature has been included in the table).

In many cases, the purpose of the research conducted throughout the twentieth century involved studying certain characteristics of anemometer performance to obtain experimental data in order to develop mathematical models. It must be underlined that a validated mathematical model to predict anemometer performance under normal working conditions is a very useful and important tool in different fields, such as meteorology and the wind energy industry.

Mathematical models normally include Euler’s equation for describing the rotation of a rigid body (anemometer cup rotor), affected by both aerodynamic and friction torque [5]: where is the rotational speed of the anemometer rotor, is the moment of inertia, is the aerodynamic torque, and is the frictional torque that depends on the air temperature, , and the rotation speed, (from [10]: , where coefficients , , and are negative (the friction torque, , in expression (1) has a negative sign in the updated version of [5]. Therefore, coefficients , , and of the friction torque expression will be positive, if this is taken into account)). The frictional torque, , is generally neglected in all mathematical models [10], as it is very low compared to the aerodynamic torque within the normal wind speed range [11, 12]. The aerodynamic torque is modeled as a function that includes the wind speed, , the cup center rotation speed, ( is the cup center rotation radius, see Figure 1), and the vertical component of the wind speed, . This function is derived using nondimensional and perturbation analysis to the second-order Taylor polynomial, where the coefficients are measured by means of carefully planned testing [11, 1319]. This method was therefore developed to correlate a group of parameters with a specific individual anemometer, in order to obtain the highest possible accuracy in the predictions relating to dynamic behavior (i.e., accelerations and decelerations of the rotor from a steady situation, at constant rotational speed). This procedure requires no aerodynamic model of the rotating body or rotor.

However, a properly developed aerodynamic model of the rotor can be very interesting, as the effect of certain parameters can be identified to provide a better understanding of cup anemometer performance. Some efforts have been made to integrate rotational aerodynamics in the mathematical models with interesting results [14, 2024]. Nevertheless, it is fair to say that, due to the simplifications in the development of these models, some deviations remain when compared to the experimental results [21, 22].

The aim of this paper is to continue with the postprocessing of the experimental results from the systematic cup anemometer testing campaign carried out at the IDR/UPM Institute during 2011 and to derive a simple mathematical model for studying some aspects of cup anemometer performance. The testing campaign consisted of 21 calibrations performed on 2 different anemometers (Climatronics 100075 and Ornytion 107A, see Figure 1), while varying the cup radius, , and the cup center rotation radius, , of the rotor (more information regarding the testing campaign can be found in [21]). Table 2 includes the calibration results from that campaign (details concerning the calibration facility are included in Appendix A). These results are the slope and offset of the transfer function (commonly known as calibration constants and [25, 26]): This transfer function, obtained experimentally in a wind tunnel, relates the wind speed to the output frequency of the anemometer, (see Appendix A). The slope of the transfer function given in terms of rotational frequency, , instead of output frequency, : is also included in Table 2. This new slope, , is the result of multiplying calibration constant by the number of pulses per turn, , given by the anemometer. (The Climatronics 100075 and Ornytion 107A anemometers give 30 and 2 pulses per turn, respectively [20].) Finally, the coefficients , , , and , (slope and offset) of the linear fittings with regard to the calibration coefficients, and , as a function of the cup center rotation radius, calculated for rotors with cups of the same radius, , are included in Table 2, along with the aforementioned results. The above equations were important for deriving the following experimentally fitted expressions [21]:

The most relevant conclusions reached were as follows.(i)The slope of the calibration transfer function, , depends on two different contributions, one related to the cup center rotation radius, , and the other related to the cups’ front area, , or cup radius,   ( ). The slope of expression (6), , seems to be directly related to the aerodynamic nondimensional coefficient of the cups, as very small differences in this coefficient were observed among the 42 calibrations performed on the 2 anemometers tested. That is, the fitting coefficient did not seem to depend on the anemometer, with the same value for both the Climatronics 100075 and the Ornytion 107A anemometers, whereas the other fitting coefficients, , , and , were different depending on the anemometer tested.(ii)The offset of the calibration transfer function, , also depends on the same shape parameters, although in this case each contribution is not totally independent from the next. In this case, all the fitting parameters, , , , , and , were different depending on the anemometer tested.

In the research described, cup anemometer performance was analyzed using the anemometer factor, , which relates the wind speed, , to the cup center rotation speed, ; that is, . In the calculations of this factor from the calibration results, the offset constant, , is usually neglected because it is small compared to the wind speed within the normal working range of the anemometer (average values of calibration constant for three Class-1 anemometers are 0.179 m s−1 (Risø P2546A), 0.248 m s−1 (Thies Clima 4.3350), and 0.184 m s−1 (Vector Instruments A100 L2) [20]. The calibration range of an anemometer, according to MEASNET is from 4 m s−1 to 16 m s−1 [25], although sometimes the upper limit of this calibration range is larger [20]). Nevertheless, the anemometer factor does indeed depend on the offset constant: Moreover, in a recent study at the IDR/UPM Institute, differences of up to 13.4% were observed at a wind speed of m/s regarding the anemometer factor, , if the offset constant is taken into account [4]. To avoid this possible source of misunderstandings, a simplified anemometer factor, , is proposed in this paper: Figure 2 includes the simplified anemometer factor (hereinafter, anemometer factor), , calculated with the results from Table 2. It must also be said that in previous studies [4, 22], the anemometer factor was regarded as displaying quadratic behavior as a function of the ratio of cup radius to cup center rotation radius, . This consideration was based on the assumption of asymptotic behavior for low values of the aforementioned parameter ( ). This seems reasonable because in this case the analytical models showed a lower dependence on [22]. However, the linear fittings to the data from Figure 2 show rather high determination coefficients, , with minimum and average values of and , respectively. This indicates linear rather than quadratic behavior within the normal range of parameter . Also, taking into account expressions (6) and (9), the following equation can be derived: which is in agreement with the linear behavior mentioned above. Obviously, in the expression above, and are derived from the fitting coefficients and in (6).

To develop a mathematical model, the aerodynamic torque should be considered proportional to the dynamic pressure relative to the cup. On the other hand, because a dynamic problem is considered, the aerodynamic torque should be described in terms of a nondimensional parameter, , formed by the rotational speed, , which is inversely proportional to the characteristic time of the movement, and the residence time, , comprised the wind speed, , and a characteristic length of the rotor, . This parameter, , is also called the cup center nondimensional velocity.

In this regard, Figure 3 shows the experimental results from Brevoort and Joyner [27] and Wyngaard et al. [13]. In these graphs, the aerodynamic torque, , is shown as a function of the dimensionless parameter . The nondimensional aerodynamic torque, (see expression (12)), is also shown. It can be observed that the nondimensional curves tend to collapse into a single curve, revealing the relationship between the aforementioned dimensionless parameters, and . This behavior was also analyzed by Pedersen [10], who found a second-order polynomial relationship between the nondimensional aerodynamic torque, , and the speed ratio (the speed ratio is defined in [5] as and , which are the wind speed at the calibration facility and the threshold velocity, respectively. The threshold velocity is derived by subtracting the friction effects from the calibration offset, ), .

The classical model for aerodynamic torque is represented by the following expression: where is the air density, the number of cups, the front area of the cups, the wind speed, the rotational speed, (concave side) and (convex side) are the aerodynamic drag force coefficients of the cups, and and are coefficients that take into account the effectiveness of the aerodynamic simplification. The average torque produced by each cup is then expressed as a function of the forces at only two positions. See Figure 4 for the normal-to-the-cup aerodynamic coefficient, , with regard to a nonrotating cup expressed as a function of the wind direction, . The dashed line in the figure represents the approximation of aerodynamic torque produced by one cup shown in expression (11). That expression can be rewritten in nondimensional form: where the dependence on the nondimensional parameter is reflected. Also, a new parameter is included in the equation, the drag coefficient ratio, , With this model, the equilibrium point, , is obtained for

Theoretical equation (12) has been fitted to the graphs for the nondimensional torque included in Figure 3. The results are in good correspondence with the experimental results, indicating the validity of this analytical approach. The simplification made with (11) is also known as the 2-cup analytical model, which models the anemometer’s behavior with a rotor consisting of cups whose aerodynamic moment along one turn is characterized by the two more relevant positions with respect to the wind, and , see Figure 4. For the equilibrium point, , the rotor can be idealized as a 2-cup rotor with cups, respectively, positioned at the aforementioned angles (wind pointing to the concave-side of one cup) and (wind pointing to the convex-side of the other cup). This situation is sketched in Figure 4.

1.3. Aim of the Present Work

As explained in the previous subsection, the aim of the present paper is to derive a mathematical model with a strong aerodynamic basis, which could help to explain its behavior as a function of certain parameters. To do so the paper has been organized as follows. Section 2 describes a mathematical (analytical) model for analyzing cup anemometer behavior. This model is based on the classical 2-cup positions model explained earlier. Although the results from the classical model show some limitations for drag coefficient ratios outside the normal ranges [21] (Brevoort and Joyner [27] measured drag coefficients ratio, , for different cup shapes in the bracket , whereas Schrenk [28] measured for hemispherical cups), it has proven strong enough to provide well-founded explanations of cup anemometer performance [14, 28] and also correlated well with data resulting from the analysis of several commercial cup anemometers [20]. In Section 3, the model is correlated with the results from Table 2. Finally, conclusions are summarized in Section 4.

2. Mathematical Model

As stated above, the aim of this paper is to derive a model simple enough to analytically study cup anemometer behavior but also complex enough to include the parameters most relevant to the problem. In this sense, simple classical analytical approximations are unable to reflect the influence of cup size (i.e., cup radius, ), in the solution for the equilibrium state ( ; see expressions (12), (13), and (14)) [21, 22]. To overcome this limitation, a nonconstant force distribution on the cup is considered.

The normal-to-the-cup force, (see Figure 4), is generally expressed as where is the relative-to-the-cup wind speed, is the cup’s front area, and is the normal-to-the-cup force coefficient (which depends on the angle of with respect to the cup). This force coefficient is normally measured by static wind tunnel testing of an isolated cup, which does not take into account any rotating flow around the cup [21, 22, 27]. Nevertheless, expression (15) is a good approximation for low values of the parameter (i.e., large cup center rotation radius, , in comparison to the cup radius, ). Some efforts have already been made to include the aforementioned rotating flow effect, considering the dynamic pressure relative to the cup as a function of the distance to the rotation axis [22]. However, the results were not relevant because the normal-to-the-cup aerodynamic coefficient, , was still considered constant all over the front surface of the cup, and was based on static wind tunnel measurements. As a result, this paper considers a non-constant pressure distribution on the front area of the cups. The aerodynamic forces on the cup arm are also considered in the mathematical approximation.

The torque on one cup at position (see Figure 4) is then defined by the following equation: where is the distance to the rotation axis, is the width of the cup at position ,    is the drag coefficient of each section of the cup, and ,   ,   , and are the wind speed, rotational speed, cup radius, and cup center rotation radius, respectively. Finally, is the diameter of the arm that attaches the cup to the rotor’s head and is the drag coefficient of the cross section of that arm (obviously, referred to the diameter ). In addition, at position , the torque on the cup is defined by the following equation: where is the drag coefficient of each section of the cup in the new position (note that torques and are defined in different rotational directions). For convenience, the new variable is considered hereinafter. If the classical model is considered, the average aerodynamic torque on a cup in a rotor is characterized by the two positions mentioned ( and ). So, as in the equilibrium (autorotation) state, this aerodynamic torque on the rotor must counterbalance the friction torque, . That condition can be expressed as where is now a proportional constant that depends on the number of cups. In the classical 2-cup model, expression (11), stands for the exact number of cups, nevertheless, if friction is taken into account and the model is fitted to experimental results without introducing any coefficient to take into account the aerodynamic simplification ( and in expression (11)), it seems more reasonable to consider the ratio between aerodynamic and friction torques related to the number of cups, and proportional to it in a first approximation. In addition, a linear law is considered for the nonconstant drag coefficients: where and are the average drag coefficients of the cups ( and in (11) and (12)), is the distance measured from the center of the cup (see Figure 1), and and are slope constants of the mentioned linear laws. Finally, a friction torque coefficient, , is defined as Taking into account these last expressions, (18) can be rewritten as If the same linear law is considered for both and ; that is, and coefficient is now defined as then the following expression can be derived from (21): where As previously indicated, is the ratio of cup radius to cup center rotation radius in the equations above.

Finally, a simple expression can be obtained: where It must also be said that expression (26) is dependent on coefficients ( ,   ,   ,   , and   ), which are low when compared to 1. This makes it possible to derive solutions using Taylor series expansions. Also, if friction is considered negligible ( ), the drag coefficient distributions along the cups are considered constant ( ), and no aerodynamic drag is considered on the cup arms ( ); expression (26) turns into and the solutions are Obviously, the only logical solution is the one with a negative sign ( ): As stated in Section 1, this is the solution for the simpler 2-cup modeling of cup anemometer behavior.

Going back to (26), if small perturbations are considered (drag distribution on cups, friction, etc.), it is possible to find an approximate asymptotic solution: where . In addition, considering as stated above, the following expression can be derived (see Appendix B): Finally, (31) can be expressed as

In Figure 5, the exact and approximate solutions for the proposed method ((24) and (33), resp.) are fitted to the experimental results (Table 2), without taking into account the offset of the transfer function; that is, (see expression (9)). These solutions, and , were calculated without considering the effect of the rotor arm or friction ( and ). As shown in the figure, there is good agreement between the solutions. Several conclusions can be derived from (32). As expected, the contribution of the rotor arm and friction, and , tends to reduce the steady rotational speed, . On the other hand, the fitting to the experimental results reveals a negative value of , indicating a higher average aerodynamic load on the closest-to-the-rotating-axis area of the cup. This is a rather surprising result, which could be produced by both unsteady and rotational aerodynamic effects.

In order to analyze anemometer transfer function (3) using the model developed to introduce the effects of the rotor arms, friction, and the aforementioned “inverse” distribution of the aerodynamic load on the cup, the anemometer factor, , must be expressed as a function of : and then The first term in the above equation directly affects the first term in transfer function (3), , and can be divided into three other terms. The first, only depends on the cup aerodynamics and not on the slight effects considered. The second one, implies a correction of the slope, linear with , that would be included in the second term in expression (10), respecting the negative sign as . Finally, the third term, can also be separated into three different contributions to the slope of the transfer function:(1)one correction not dependent on (2)one correction linear with , which also respects the negative sign in the second term in expression (10) (3)and, finally, a nonlinear term

With regard to the second term in (35), it represents a friction velocity term, , which can be expressed as If the second-order polynomial approximation for the friction torque mentioned in Section 1 is considered: then (42) turns into where the first term, indicates the deviation from the ideal velocity (36) at low rotational speeds. This effect is mentioned in Section 1, where the non-linear behavior of anemometer constant at low wind speeds is described. The second term, involves a contribution that does not depend on , so it represents the offset constant, , in transfer function (3). Finally, the third term, represents a contribution to transfer function slope .

From expression (42), it can also be assumed that the friction term, , is inversely proportional to the cup surface, , meaning that bigger cups would be translated into a lower friction term. Nevertheless, it must be taken into account that bigger cups also experience greater lateral forces on the anemometer rotor shaft, increasing the friction.

As mentioned, the results included in Table 2 summarize the results from [21]. Bearing in mind these results (see expressions (3) to (7)), it is possible to derive the following equation: where is the contribution to the anemometer factor that depends on the cup center rotation radius, , and: is the contribution to the anemometer factor that does not depend on the cup center rotation radius but on the cup radius, . Both coefficients, and , are dimensionless and can be related to the anemometer factor defined in expression (9) through the following equation: If we compare expression (48) to the same expression obtained from the theoretical model (35), neglecting the effect of rotor arms or friction ( ), that is, the two terms of the aforementioned (48) can be explained from a theoretical basis. First, the slope of the anemometer calibration curve (a transfer function based on the rotation frequency; that is, expression (3) instead of expression (2)), has one term, , proportional to the cup center rotation radius, , with the proportionality constant . Second, the last term of the aforementioned anemometer transfer function slope, , which does not depend on , depends on both the load (pressure) distribution along the cups, , and the cup radius, . It should be underlined that this second term in the equation of the anemometer factor only appears in the theoretical model if the “inverse” force distribution along the cup is considered.

Comparing expressions (51) and (52) leads to Finally, combining expressions (49) and (53) with the above equation it is possible to rewrite expression (51) in terms of the mentioned ratio of cup radius to cup center rotation radius, , as which, taking into account the negative value of coefficient , results in an equation very similar to the one derived from the experimental results (10). Nevertheless, it must also be noted that both terms in expression (10) include the slight effects of rotor arms and friction, which are not considered in the above equation.

3. Results and Discussion

In Figure 6, anemometer factor , obtained from the calibrations performed on the anemometers used in the testing campaign, is shown as a function of the parameter , for rotors with the same cup center rotation radius, . The figure also compares these experimental results to those from the proposed analytical model, without taking into account the effect of rotor arms or friction ( ). In each case, the complete expression from the analytical model (24) is fitted by adjusting both and . Table 3 includes the values of these parameters for each case, along with the slope, , and offset, , of the linear fittings and the corresponding determination coefficient, . Obviously, and , from expression (51).

Good agreement between the proposed model and the testing results is observed in Figure 6; in other words, the effect of parameter is correctly reflected by the model. However, it is also fair to say that greater differences between the model and the experimental results are observed in the figure for higher values of . This dispersion of the testing results could be explained by the fact that a longer cup center rotation radius, , lowers the rotation speed, making the rotation movement less constant and more influenced by the third harmonic term [4, 136] (the rotation speed of a cup anemometer can be decomposed into different harmonic terms using the Fourier expansion: . Leaving aside the average term, , the third harmonic term, , is the most important term of the expansion in normal circumstances due to the 3-cup rotor shape).

To return to the linear fittings in Table 3, taking into account the high values of determination coefficient , it seems that the linear approximation previously derived (51) reflects the anemometer behavior quite well. It is also possible to observe the variation of both coefficients, and , as a function of , which is plotted in Figure 7. For both anemometers, Climatronics 100075 and Ornytion 107A, these coefficients show approximate linear behavior for  mm: For lower cup center rotation radius values,  mm, the tendency of both coefficients, and , to decrease the anemometer factor as a function of is alleviated, probably due to cup-wake interaction. Table 4 includes the coefficients of the aforementioned linear fitting to and for  mm (expressions (55)). The following conclusions can be reached from these fittings.(i)The behavior of coefficient is very similar for both anemometers, an average expression being: . ( expressed in mm). This parameter seems to depend only on the rotor shape (cups size and rotor diameter) and not on the anemometer shape (anemometer body).(ii)The linear fittings to parameter show similar values for the slope, , but different values for the offset, . This difference could be explained by different friction on the shaft/bearings system or a different anemometer shape (“neck” thickness).From expressions (51) and (55) it is possible to obtain, in combination with the anemometer factor equation, a new expression for wind speed, that is, an expression indicating a term related to the anemometer transfer function slope, , which does not depend on . The effect of this term was experimentally observed [21]; see expression (6).

To return to (53), it is possible to estimate the gradient along the cup of the normal-to-the-cup aerodynamic force, (see expressions (19) and (22)), Figure 8 shows coefficient as a function of the anemometer cup radius, . Considering an average value between (Climatronics 100075) and (Ornytion 107A), the above equation can be rewritten as The above equation clearly indicates the aforementioned “inverse” aerodynamic load distribution on the anemometer cups. In other words, the aerodynamic load on the rotating cup is higher in the area closest to the rotating axis. This effect, which could be attributed to local changes in wind direction along the cup, to the effect of rotating flow on the pressure distribution along the cup, or to wake interaction, should not be neglected in the development of new models for studying cup anemometers.

In a previous paper [21], the offset of the transfer function (coefficient in expressions (2) and (3)) was experimentally fitted to the expression, depending on the cup center rotation radius, , and the front area of the cups, (see (7)). Taking that expression into account, offset depends linearly on . Figure 9 shows coefficient for the calibrations performed (see Table 2) as a function of the cup center rotation radius, , for each cup radius, . Linear fittings to the data are also shown in the graphs (the coefficients for these fittings, and , are included in Table 2). The linear trend mentioned earlier is shown in the graphs, although it is much clearer in the case of the Ornytion 107A. This could be explained by the different effects of aging or wear and tear, on both anemometers. In the testing campaign considered in this paper, the Ornytion 107A was new, whereas the Climatronics 100075 had been used for internal procedures at the IDR for several years (some degree of degradation regarding this anemometer was previously illustrated in [82]). Another effect observed in the curves corresponding to and for the Climatronics 100075 is a deviation from the linear trend for lower values of . This suggests an interaction between the cup and the wake generated at the anemometer’s “neck.” This effect is not observed in the graph for the Ornytion 107A, probably because this anemometer has a thinner “neck,” which also produces a thinner wake. Finally, a decrease in the slope, , along with the cup radius, , can be observed in Figure 9.

4. Conclusions

In this study, cup anemometer response was analyzed using the 2-cup analytical model. The model was fitted to experimental data from 42 calibrations performed on two different cup anemometers (Climatronics 10075 and Ornytion 107A), equipped with 21 different cup rotors (conical cups, with varying cup sizes and distances to the rotation axis). The major conclusions resulting from this study are as follows.(i)The anemometer factor, , for cup anemometers equipped with conical cups depends linearly on the shape parameter (the ratio of the cup radius to the cup center rotation radius, ), , within the range studied (from to ). The slope of this linear behavior, , seems to depend only on the rotor shape but not on the anemometer body (for the anemometers tested).(ii)The results indicate that the aerodynamic force distribution along the rotating cups is not uniform, as the calculated aerodynamic force is higher closer to the rotating axis, which explains the linear correction term to the anemometer factor.(iii)The classical 2-cup analytical model modified with the nonconstant force distribution along the rotating cup of the anemometer seems to be accurate enough for studying the complex aerodynamics effects involved in the rotor performance.



Anemometer calibrations are performed in the S4 wind tunnel at the IDR/UPM Institute. This is an open-circuit wind tunnel with a closed test section measuring 0.9 by 0.9 m. It is served by four 7.5 kW fans with a flow uniformity under 0.2% in the testing area. The wind speed in the testing chamber is measured by Airflow 0.48 Pitot tube connected to a Druck LPM 9481 high-precision pressure transducer, with the electrical signal from the pressure transducer measured by a Keithley 2000 digital multimeter. Temperature and humidity sensors (Vaisala PTU 200 and Vaisala HMP45D) are used to determine the air density value. The rotation frequency of the anemometer is measured with an Agilent 53131A universal counter. Another digital multimeter is used to measure the voltage or current output from the anemometer when required.

An example of the calibrations performed on this wind tunnel is included in Figure 10. The wind speed measured by the wind tunnel instruments, , is plotted in comparison to the anemometer’s frequency output, . The calibration curve (i.e., the linear transfer function) corresponding to the AC calibration is plotted in the aforementioned figure (its mathematical expression is also included). The uncertainty levels of the calibrations performed at the S4 calibration wind tunnel are specified following the ISO/IEC 17025 standard [137], these levels being 0.1 m s−1 for wind speeds from 4 m s−1 to 10 m s−1 and 0.01   m s−1 for wind speeds, , from 10 m s−1 to 23 m s−1.


The aim of this Appendix is to find the roots of equation: in the case .

If the following perturbation solution is assumed: then The solution is as follows, if : where If (B.4) is subtracted from (B.3), taking into account , the following expression is obtained: and then where Once the problem has been linearized, correction depends on perturbations and .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors are indebted to Enrique Vega, Alejandro Martínez, Encarnación Meseguer, Luis García, and Eduardo Cortés for the help regarding the testing campaign and, more generally, all the experimental work carried out. The authors are grateful to Tania Tate for her kind help on improving the style of the text. Prof. Pindado and the other coauthors are truly indebted to Luis Gómez, Eva Díaz, Javier Medel, Clarisa Buendía, and the rest of the staff of the Library at the Aeronautics and Space Engineering School (Escuela de Ingeniería Aeronáutica y del Espacio) of the Polytechnic University of Madrid (Universidad Politécnica de Madrid), for their constant support to the research carried out regarding cup anemometer performances.