Comparison of Power Coefficients in Wind Turbines Considering the Tip Speed Ratio and Blade Pitch Angle

: This paper presents a review of the power and torque coefﬁcients of various wind generation systems, which involve the real characteristics of the wind turbine as a function of the generated power. The coefﬁcients are described by mathematical functions that depend on the trip speed ratio and blade pitch angle of the wind turbines. These mathematical functions are based on polynomial, sinusoidal, and exponential equations. Once the mathematical functions have been described, an analysis of the grouped coefﬁcients according to their function is performed with the purpose of considering the variations in the trip speed ratio for all the coefﬁcients based on sinusoidal and exponential functions, and with the variations in the blade pitch angle. This analysis allows us to determine the different coefﬁcients of power and torque used in wind generation systems, with the objective of developing algorithms for searching for the point of maximum power generated and for the active control of wind turbines with variations in the blade pitch angle.


Introduction
The progressive growth in global electricity demand results in a deterioration of the environment, which manifests in higher emissions of carbon dioxide (CO 2 ) and an increase in the greenhouse effect.This is mainly due to the fact that this energy demand has mostly been met by fossil fuels (oil, coal, etc.) [1,2].The high dependence on fossil fuels, their involvement in "climate change", and the progressive increase in their costs has motivated interest in other alternative forms of electricity generation that do not depend on fossil fuels [3].These alternative forms are known as renewable energy sources, and they have led to a diversification of energy generation [4,5]; these new forms include wind energy and solar energy.The wind is a clean and inexhaustible resource that is available in most parts of the world; however, there are specific areas in which this natural resource is abundant [1,[6][7][8].Different statistical methods have been developed to calculate the average wind speed, the wind energy density, and its load factor for a specific geographic area [9]; this allows for the identification of suitable zones for the installation of wind farms of different MWs and low power systems in the order of kWs.Several projects have helped to increase the production of clean electrical energy.At present, there are numerous projects aiming to start up several wind farms, increasing the production of clean electric energy.Meanwhile, the administrations of several countries have established incentives and subsidies for these projects, to comply with the Kyoto Protocol and achieve a reduction in CO 2 emissions.
The recent advances in wind generation systems (WGS) have led to a reduction in costs, bringing them to highly competitive levels compared to the systems of electric generation using conventional energies (oil, coal, etc.) [10].
The process of converting wind energy to electric power originates from the wind gusts colliding with the blades of a wind turbine, which turns an electric generator because it is interconnected with a gearbox.The electric generator provides electric power that is variable in both amplitude and frequency, which is not suitable for common working loads; in addition, the wind is random in nature, so the amplitude and frequency of voltages provided by the generator vary constantly.In order to condition the variable voltages of the generator and to supply loads with amplitude and frequency values according to the regulations of a given country, electronic power systems are used to carry out this conversion [11,12].
One of the most important parameters for evaluating the performance of a wind turbine is the power coefficient, as established by the International Electrotechnical Commission (IEC) in its recommendation IEC 61400-12-1 [13].The power coefficient (Cp) refers to the relationship between actual power production and the power available in the ambient wind flowing towards the turbine blades [14].The power coefficient is one of the parameters that is required to calculate the productive efficiency of wind turbines or wind farms, where parameters such as turbulence, temperature, and air density are also involved; in addition, data science is used to analyze all these parameters and make them more efficient [15,16].
The present paper focuses on the behavior of the different wind turbines that are presented in the literature.The wind turbines analyzed in this work are the horizontal axis and three-blade types, because they are the most commonly found.This paper is structured as follows.It begins by showing the characteristics of a wind turbine in order to establish which are the most important parameters that establish its behavior and cause each turbine to have particular characteristics in the process of converting wind energy.With this, it is established that one of the most significant parameters that differentiates one wind turbine from another is the Power Coefficient.The power coefficient determines the power variation that can be obtained from a wind turbine; it should be noted that each model of wind turbine has its own power coefficient.This paper shows the mathematical functions that are used to represent the behavior of the power coefficient in relation to the specific speed and the angle of attack of the blade, establishing general formulas that allow for a better analysis and comparison of the power coefficients.Finally, a discussion of the review based on the practical cases presented in the literature is carried out; this enables us to establish the limitations and the advantages of each of the power coefficients considered in the analysis.The power coefficient can be used for other analyses of wind generation systems, such as the production efficiency of wind turbines.

Power and Torque in a Wind Turbine
The generated power of a wind turbine is stable based on the kinetic energy in the air of a mass object (m), which moves at a speed (v).Therefore, the wind kinetic energy is given by: The power of the moving air is obtained from the wind energy, assuming that the wind speed (v ω ) is constant and can be expressed by: where .
m is the mass airflow per second.If the air passes through an area (A) and the density of the air is (ρ), the mass airflow is expressed by: .m = ρAv ω (3) Energies 2023, 16, 2774 3 of 23 substituting (3) into (2), we obtain: where A is the area covered by the blades of the turbine, and can be expressed by: substituting ( 5) into (4), we obtain: However, the wind kinetic energy is not converted in its entirety into mechanical energy.This is because the principle of continuity of the air flow passing through the turbine must be met, where the wind passing through the turbine is braked, so it leaves at a lower speed, but is never equal to zero.In an ideal wind turbine, the amount of kinetic energy of the wind that turns into mechanical energy is 59.25%, which is established in Betz's law.This maximum conversion value is known as the Betz limit.Therefore, a term is added that relates the power of the wind and the power that can be obtained in the turbine, this relationship is known as the power coefficient (C p ).Therefore, the power in the turbine (P t ) is expressed by: The power coefficient depends on several factors, among which are the type of turbine, that is, if it is a horizontal or vertical axis type; in the case of a horizontal axis, this also depends on the number of blades, the specific speed or trip speed ratio (λ), and the blade pitch angle of the turbine (β).
In this way, for the same type of turbine, the power coefficient depends on λ and the blade pitch angle, C p ( λ, β).
The symbol λ relates to the tangential velocity of the end of the turbine blade between the wind speed.This is expressed by: where ω m is the rotational speed of the turbine.The power coefficient C p depends on both λ and the type of turbine, as shown in Figure 1.The variation of C p with respect to β and λ is shown in Figure 2.
Energies 2023, 16, x FOR PEER REVIEW 3 of 23 substituting (3) into (2), we obtain: where  is the area covered by the blades of the turbine, and can be expressed by: substituting ( 5) into (4) , we obtain: However, the wind kinetic energy is not converted in its entirety into mechanical energy.This is because the principle of continuity of the air flow passing through the turbine must be met, where the wind passing through the turbine is braked, so it leaves at a lower speed, but is never equal to zero.In an ideal wind turbine, the amount of kinetic energy of the wind that turns into mechanical energy is 59.25%, which is established in Betz's law.This maximum conversion value is known as the Betz limit.Therefore, a term is added that relates the power of the wind and the power that can be obtained in the turbine, this relationship is known as the power coefficient (  ).Therefore, the power in the turbine (  ) is expressed by: The power coefficient depends on several factors, among which are the type of turbine, that is, if it is a horizontal or vertical axis type; in the case of a horizontal axis, this also depends on the number of blades, the specific speed or trip speed ratio (), and the blade pitch angle of the turbine ().
In this way, for the same type of turbine, the power coefficient depends on λ and the blade pitch angle,   ( , ).
The symbol  relates to the tangential velocity of the end of the turbine blade between the wind speed.This is expressed by: where   is the rotational speed of the turbine.The power coefficient   depends on both  and the type of turbine, as shown in Figure 1.The variation of   with respect to β and λ is shown in Figure 2.  The other important parameter involved in the process of converting the energy in the turbine is the torque of the turbine shaft, which is related to the power in the turbine.This relationship is expressed by:   =     (9) substituting ( 7) into ( 9) and solving   , we obtain: substituting ( 8) into (10) gives: From (11), the torque coefficient (  ) is established by: substituting ( 12) into (11) gives: The variation of   with respect to λ and the type of turbine is shown in Figure 3, and the same variation in relation to β and λ is shown in Figure 4.The other important parameter involved in the process of converting the energy in the turbine is the torque of the turbine shaft, which is related to the power in the turbine.This relationship is expressed by: substituting ( 7) into ( 9) and solving T t , we obtain: substituting ( 8) into (10) gives: From (11), the torque coefficient (C t ) is established by: substituting ( 12) into (11) gives: The variation of C t with respect to λ and the type of turbine is shown in Figure 3, and the same variation in relation to β and λ is shown in Figure 4.The other important parameter involved in the process of converting the energy in the turbine is the torque of the turbine shaft, which is related to the power in the turbine.This relationship is expressed by:   =     (9) substituting ( 7) into ( 9) and solving   , we obtain: substituting ( 8) into (10) gives: From (11), the torque coefficient (  ) is established by: substituting ( 12) into (11) gives: The variation of   with respect to λ and the type of turbine is shown in Figure 3, and the same variation in relation to β and λ is shown in Figure 4.

Power Coefficient and Torque Coefficient
In the literature, there are different graphical representations of   ; in some ca these differ considerably.In the case of commercial wind turbines, manufacturers prov documentation where the   behavior graphs are provided.These graphs are used to tain a mathematical approximation that allows us to understand their behavior in an alytical way.These approximations are carried out by means of optimal numerical m ods [18].
From conducting an exhaustive bibliographical review, it was found that the ma matical approximations of the variation in the power coefficient are expressed thro three mathematical functions: The following section shows all of the potency coefficients found in the literat grouped based on their mathematical function, as well as the parameters of the wind g eration system where it is used, in order to establish their main characteristics.

Polynomial Function
Four power coefficients were found from a review of the literature, and their ma matical expression has the form of a polynomial function.Each of the power coefficie represented by a polynomial equation is detailed below; these are identified by the or of the polynomial.From the bibliographic review, four power coefficients were fou where their mathematical expression has the form of an exponential function.These identified according to the order of the polynomial function presented by the author the paper.

Third-Order Polynomial Function
In [19,20], the authors show the hardware implementation of control algorithms constant wind speed turbine emulator based on a DC-Machine.This system includ wind speed profile, a mathematical model of a horizontal axis wind turbine, a separa excited DC-Machine, and a fourth quadrant chopper.The power coefficient only depe on λ and is expressed by:   () = −3.7 × 10 −5  3 − 0.004834 2 + 0.1063 − 0.02086 ( This power coefficient is used for a wind generation system with a blade radiu 1.086 m, and the system has a power of 3 kW.This paper presents the simulation res

Power Coefficient and Torque Coefficient
In the literature, there are different graphical representations of C p ; in some cases, these differ considerably.In the case of commercial wind turbines, manufacturers provide documentation where the C p behavior graphs are provided.These graphs are used to obtain a mathematical approximation that allows us to understand their behavior in an analytical way.These approximations are carried out by means of optimal numerical methods [18].
From conducting an exhaustive bibliographical review, it was found that the mathematical approximations of the variation in the power coefficient are expressed through three mathematical functions: The following section shows all of the potency coefficients found in the literature grouped based on their mathematical function, as well as the parameters of the wind generation system where it is used, in order to establish their main characteristics.

Polynomial Function
Four power coefficients were found from a review of the literature, and their mathematical expression has the form of a polynomial function.Each of the power coefficients represented by a polynomial equation is detailed below; these are identified by the order of the polynomial.From the bibliographic review, four power coefficients were found, where their mathematical expression has the form of an exponential function.These are identified according to the order of the polynomial function presented by the authors in the paper.

Third-Order Polynomial Function
In [19,20], the authors show the hardware implementation of control algorithms for constant wind speed turbine emulator based on a DC-Machine.This system includes a wind speed profile, a mathematical model of a horizontal axis wind turbine, a separately excited DC-Machine, and a fourth quadrant chopper.The power coefficient only depends on λ and is expressed by: C p (λ) = −3.7 × 10 −5 λ 3 − 0.004834λ 2 + 0.1063λ − 0.02086 (14) This power coefficient is used for a wind generation system with a blade radius of 1.086 m, and the system has a power of 3 kW.This paper presents the simulation results in Matlab/Simulink.In [21,22], the authors present an isolated small wind turbine emulator based on a separately excited DC motor.It was developed to emulate and evaluate the performance of a small wind turbine using different control strategies.The power coefficient only depends on λ and is expressed by: C p (λ) = 0.00044λ 4 − 0.012λ 3 + 0.097λ 2 − 0.2λ + 0.11 (15) This paper presents experimental results in Matlab/Simulink.

Fifth-Order Polynomial Function
In [23,24], an improved maximum-power-point tracking algorithm for wind-energy conversion systems is presented.The proposed method significantly reduces the turbine mechanical stress related to conventional techniques, so that both the maintenance needs and the medium time between failures are expected to be improved.The power coefficient only depends on λ and is expressed by: C p (λ) = 0.043 − 0.108λ + 0.146λ 2 − 0.0605λ 3 + 0.0104λ 4 − 0.0006λ 5 (16) This power coefficient is used for a wind generation system with a blade radius of 1.525 m and the system has a power of 2 kW.This paper presents simulation and experimental results.

Sixth-Order Polynomial Function
In [25], the authors present a wind turbine emulator (WTE) that is designed and implemented considering different requirements from the development and testing of a control strategy in a doubly fed wind-power-generating system.The power coefficient only depends on λ and is expressed by: C p (λ) = 0.051 − 0.0022λ + 0.0052λ 2 − 5.1425 × 10 −4 λ 3 − 2.795 × 10 −5 λ 4 + 4.6313 × 10 −6 λ 5 −1.331 × 10 −7 λ 6 (17) This power coefficient is used for a wind generation system with a blade radius of 2.5 m and the system has a power of 4.2 kW.This paper presents experimental results.

General Exponential Function
Based on the general polynomial function shown in [25,26], this is expressed by: All the power coefficients based on a polynomial function are grouped according to (18).Considering that the maximum order found in the functions is the sixth, the equation expressed by (19) is obtained.The constants of each polynomial function are shown in Table 1.
The behavior of the four power coefficients based on a polynomial function in relation to λ are shown in Figure 5a.The graph indicates the maximum value of the power coefficient (C p max ), as well as the value of λ for that point (λ max ), in each of the functions found.None of the coefficients exceed the Betz limit.In the established interval of λ between 0 to 18, the C p of the third, fifth and sixth orders have adequate behavior; however, the fourth-order function must be limited to an λ of 12.16 because, for higher values of λ, the value of C p begins to increase again, which does not correspond to a real behavior of a wind turbine.The torque coefficient is obtained using (12).Therefore, in Figure 5b, the behavior of the power coefficient in relation to λ is shown, based on polynomial functions of the power coefficient.The behavior of the four power coefficients based on a polynomial function in relation to λ are shown in Figure 5a.The graph indicates the maximum value of the power coefficient (   ), as well as the value of λ for that point (  ), in each of the functions found.None of the coefficients exceed the Betz limit.In the established interval of  between 0 to 18, the   of the third, fifth and sixth orders have adequate behavior; however, the fourth-order function must be limited to an  of 12.16 because, for higher values of , the value of   begins to increase again, which does not correspond to a real behavior of a wind turbine.The torque coefficient is obtained using (12).Therefore, in Figure 5b, the behavior of the power coefficient in relation to λ is shown, based on polynomial functions of the power coefficient.

Sinusoidal Function
From the bibliographic review, five power coefficients were found where their mathematical expressions have the form of a sinusoidal function.Then, each of them is detailed; these are identified by the authors of the paper where the power coefficient was found.[19] In [19], the authors show the hardware implementation of control algorithms for a constant wind speed turbine emulator based on a DC-Machine.This system includes the wind speed profile, mathematical model of a horizontal axis wind turbine, a separately excited DC-Machine, and a fourth quadrant chopper.The power coefficient used depends on both λ and β in degrees, and is expressed by:

Moussa, Bouallegue, and Kehedher
This power coefficient is used for a wind generation system with a blade radius of 1.086 m and the system has a power of 3 kW.This paper presents simulation results in Matlab/Simulink.3.2.2.Coto, García, Díaz, and Gómez [27] In [27], the authors explain the design, implementation, and simulation of a model that represents the possible conditions of wind that can occur in a given wind installation.The power coefficient used depends on both λ and β in degrees, and is expressed by:

Sinusoidal Function
From the bibliographic review, five power coefficients were found where their mathematical expressions have the form of a sinusoidal function.Then, each of them is detailed; these are identified by the authors of the paper where the power coefficient was found.[19] In [19], the authors show the hardware implementation of control algorithms for a constant wind speed turbine emulator based on a DC-Machine.This system includes the wind speed profile, mathematical model of a horizontal axis wind turbine, a separately excited DC-Machine, and a fourth quadrant chopper.The power coefficient used depends on both λ and β in degrees, and is expressed by:

Moussa, Bouallegue, and Kehedher
This power coefficient is used for a wind generation system with a blade radius of 1.086 m and the system has a power of 3 kW.This paper presents simulation results in Matlab/Simulink.
Energies 2023, 16, 2774 8 of 23 3.2.2.Coto, García, Díaz, and Gómez [27] In [27], the authors explain the design, implementation, and simulation of a model that represents the possible conditions of wind that can occur in a given wind installation.The power coefficient used depends on both λ and β in degrees, and is expressed by: This power coefficient is used for a wind generation system with a blade radius of 48.2 m and a power of 750 kW.This paper presents simulation results in Matlab/Simulink.3.2.3.Xin, Wanli, Bin, and Pengcheng [28] In [28], the authors propose a sliding-mode variable-structure controller based on the analysis of the feature model of the variable pitch and the speed wind turbine generator system, in order to improve the dynamic performance in the operational areas of constant power output in wind generation system.The power coefficient used depends on both λ and β in degrees and it is expressed by: This power coefficient is used for a wind generation system with a blade radius of 35 m and a power of 2 MW.This paper presents simulation results.[29] In [29], the authors present the DC-voltage regulation of a five-level neutral-point clamp in a closed loop.It consists of regulating the average value of the DC-voltage by using one loop instead of four loops.The modeling and the control of the different components of the wind energy conversion system are presented.The wind turbine is controlled using the maximum power point tracking algorithm based on the wind speed estimation.The power coefficient used depends on both λ and β in degrees, and is expressed by:

Merahi, Mekhilef, and Madjid
This power coefficient is used for a wind generation system with a blade radius of 70.5 m and a power system of 1.5 MW.This paper presents simulation results.[18] In [18], the authors present the implementation of an experimental wind energy board using a DC-machine as a wind emulator.The power coefficient used depends on both λ and β in degrees, and is expressed by:

Nouira and Khedher
This power coefficient is used for a wind generation system with a blade radius of 20 m and power of 10 kW.This paper presents experimental results.

General Sinusoidal Function
All the power coefficients registered in a sinusoidal function are grouped for general equation purposes based on a sinusoidal function.The general equation is expressed by (25), and the constants of each sinusoidal function are shown in the Table 2.
Energies 2023, 16, 2774 The behavior of the five power coefficients based on a sinusoidal function in relation to λ and β equal to zero are shown in Figure 6a.β is equal to zero for the purpose of comparing its behavior with the power coefficients based on a polynomial function, since this does not include β.These coefficients are identified according to the first author of the paper where the power coefficient is presented.The graph indicates the maximum value of the power coefficient (C p max ), as well as the value of λ for that point (λ max ), in each of the functions found.None of the coefficients exceed the Betz limit.In the established interval of λ between 0 to 18, the C p values established by Moussa [19], Coto [27], Xin [28] and Nouira [18] have adequate behavior.However, the power coefficient established by Merahi [29] must be limited to a λ of 10.05, because, for higher values of λ, the values of C p are negative, which does not correspond to a real behavior of a wind turbine.The torque coefficient is obtained using (12).Therefore, in Figure 6b, the behavior of the power coefficient in relation to λ is shown, based on sinusoidal functions of the power coefficient.
The behavior of the five power coefficients based on a sinusoidal function in relation to λ and β equal to zero are shown in Figure 6a.β is equal to zero for the purpose of comparing its behavior with the power coefficients based on a polynomial function, since this does not include β.These coefficients are identified according to the first author of the paper where the power coefficient is presented.The graph indicates the maximum value of the power coefficient (   ), as well as the value of λ for that point (  ), in each of the functions found.None of the coefficients exceed the Betz limit.In the established interval of  between 0 to 18, the   values established by Moussa [19], Coto [27], Xin [28] and Nouira [18] have adequate behavior.However, the power coefficient established by Merahi [29] must be limited to a  of 10.05, because, for higher values of , the values of   are negative, which does not correspond to a real behavior of a wind turbine.The torque coefficient is obtained using (12).Therefore, in Figure 6b, the behavior of the power coefficient in relation to λ is shown, based on sinusoidal functions of the power coefficient.Moussa [19], Coto [27], Xin [28], Merahi [29] and Nouira [18].
To more thoroughly evaluate the power coefficients based on sinusoidal functions, an analysis is conducted with variations of the specific velocity and β.
The behavior of the   and   proposed by Moussa, Bouallegue, and Kehedher [19], with variations of  and , is shown by Figure 7.It is observed that, with high β values, abnormal behavior occurs, because   increases again after having decreased.On the other hand, the   increases considerably when  increases, so this coefficient is not recommended for use when implementing control techniques through β.Moussa [19], Coto [27], Xin [28], Merahi [29] and Nouira [18].
To more thoroughly evaluate the power coefficients based on sinusoidal functions, an analysis is conducted with variations of the specific velocity and β.
The behavior of the C p and C t proposed by Moussa, Bouallegue, and Kehedher [19], with variations of λ and β, is shown by Figure 7.It is observed that, with high β values, abnormal behavior occurs, because C p increases again after having decreased.On the other hand, the C t increases considerably when β increases, so this coefficient is not recommended for use when implementing control techniques through β.The behavior of the   and   proposed by Xin, Wanli, Bin, and Pengcheng [28] with variations of  and  is shown in Figure 8.It is observed that the behavior of   and   is adequate, because, with high β values, the coefficients decrease, so a control through β can be implemented with this proposed coefficient.The behavior of the   and   proposed by Merahi, Mekhilef, and Madjid [29] with variations of  and  is shown in Figure 9.The behavior of   is adequate because, as β increases, the coefficient decreases; however, the same does not happen with   , because, with low β and high β, the value of   increases considerably, which does not allow for control techniques to be applied through β.The behavior of the C p and C t proposed by Xin, Wanli, Bin, and Pengcheng [28] with variations of λ and β is shown in Figure 8.It is observed that the behavior of C p and C t is adequate, because, with high β values, the coefficients decrease, so a control through β can be implemented with this proposed coefficient.The behavior of the   and   proposed by Xin, Wanli, Bin, and Pengcheng [28] with variations of  and  is shown in Figure 8.It is observed that the behavior of   and   is adequate, because, with high β values, the coefficients decrease, so a control through β can be implemented with this proposed coefficient.The behavior of the   and   proposed by Merahi, Mekhilef, and Madjid [29] with variations of  and  is shown in Figure 9.The behavior of   is adequate because, as β increases, the coefficient decreases; however, the same does not happen with   , because, with low β and high β, the value of   increases considerably, which does not allow for control techniques to be applied through β.The behavior of the C p and C t proposed by Merahi, Mekhilef, and Madjid [29] with variations of λ and β is shown in Figure 9.The behavior of C p is adequate because, as β increases, the coefficient decreases; however, the same does not happen with C t , because, with low β and high β, the value of C t increases considerably, which does not allow for control techniques to be applied through β.The behavior of the   and   proposed by Nouira and Khedher [18] with variations of  and  is shown in Figure 10.The behavior of   does not vary considerably with the variation of β; in addition,   increases considerably with low β and high blade pitch angles, which does not allow for control techniques to be applied through β.

Exponencial Function
From the bibliographic review, eight power coefficients were found where their mathematical expression has the form of an exponential function; then, each of them is detailed.These coefficients are identified by the authors of the paper where the power coefficient was found.
3.3.1.Kotti, Janakiraman, and Shireen [30,31] In [30,31], the authors present a novel sensorless adaptive control for maximum power point tracking in wind energy conversion systems.The proposed control allows the generator to track the optimal operating point of the wind energy conversion systems The behavior of the C p and C t proposed by Nouira and Khedher [18] with variations of λ and β is shown in Figure 10.The behavior of C p does not vary considerably with the variation of β; in addition, C t increases considerably with low β and high blade pitch angles, which does not allow for control techniques to be applied through β.The behavior of the   and   proposed by Nouira and Khedher [18] with variations of  and  is shown in Figure 10.The behavior of   does not vary considerably with the variation of β; in addition,   increases considerably with low β and high blade pitch angles, which does not allow for control techniques to be applied through β.

Exponencial Function
From the bibliographic review, eight power coefficients were found where their mathematical expression has the form of an exponential function; then, each of them is detailed.These coefficients are identified by the authors of the paper where the power coefficient was found.
3.3.1.Kotti, Janakiraman, and Shireen [30,31] In [30,31], the authors present a novel sensorless adaptive control for maximum power point tracking in wind energy conversion systems.The proposed control allows the generator to track the optimal operating point of the wind energy conversion systems

Exponencial Function
From the bibliographic review, eight power coefficients were found where their mathematical expression has the form of an exponential function; then, each of them is detailed.These coefficients are identified by the authors of the paper where the power coefficient was found.
3.3.1.Kotti, Janakiraman, and Shireen [30,31] In [30,31], the authors present a novel sensorless adaptive control for maximum power point tracking in wind energy conversion systems.The proposed control allows the generator to track the optimal operating point of the wind energy conversion systems under fluctuating wind conditions by an efficient two-step tracking process.The power coefficient used depends on both λ and β in degrees and is expressed by: where: here, β is defined as the angle between the plane of rotation and the blade cross section chord.This power coefficient is used for a wind generation system of 2.25 MW with simulation results and 320 W for experimental results.In both cases, the results shown only report variations in wind speed and β is constant.[32] In [32], the authors show how a variable-speed wind turbine can be used to generate a fixed value of voltage at the output with the help of a PI controller; this is achieved by varying the pitch angle of the blades.Pitch angle control is the most common means of adjusting the aerodynamic torque of the wind turbine when the wind speed is above the rated speed and various controlling variables may be chosen, such as wind speed, generator speed, and generator power.The power coefficient used depends on both λ and β in degrees and is expressed by:

Khajuria and Kaur
where: 1 In the results shown, there are variations in the wind speed and in β.

Ovando, Aguayo and Cotorogea [33-41]
In [33][34][35][36][37][38][39][40][41], the authors show the experimental implementation of a wind turbine emulator for wind energy conversion systems using a separately excited DC motor.The model of the wind turbine (WT) and the control of the DC motor are implemented in MATLAB/SIMULINK.The power coefficient used depends on both λ and β in degrees and is expressed by: where: 1 This power coefficient is used for a wind generation system with a blade radius of 1 m and the system has a power of 300 W. The results shown only report variations in wind speed and β is constant 3.3.4.Feng Gao, Da-Ping Xui, and Yue-Gang Lv [42][43][44][45][46][47] In [42][43][44][45][46][47], the authors show the typical features of a hybrid system by analyzing its characteristics.Therefore, the hybrid dynamic model of a wind generation system was built based on theory of hybrid automata; at the same time, a hybrid control system was designed to solve the problem of global automation according to the wind generation system control strategy.The power coefficient used depends on both λ and β in degrees and is expressed by: C p (λ, β) = 0.22 116 where: 1 This power coefficient is used for a wind generation system of 300 kW.The results are presented only in simulation, with variation in the wind speed.[48] In [48], the authors present the modelling and commissioning of a wind turbine emulator test rig and a performance comparison of four advanced control techniques.The power coefficient used depends on both λ and β in degrees and is expressed by: C p (λ, β) = 0.5 72.5

Llano, Tatlow, and McMahon
where: 1 This power coefficient is used for a wind generation system with a blade radius of 1.5 m and a power of 1.8 kW.The results are shown for simulations and experiments.
3.3.6.Shi, Zhu, Cai, Wang, and Yao [49,50] In [49,50], the authors present a generalized average model of the DC wind turbine, taking into consideration electromechanical transients.In the proposed wind turbine model, the power conversion system is simplified based on the average model and it can be adapted to any DC wind turbine with an active front end and two step-up levels.The power coefficient used depends on both λ and β in degrees and it is expressed by: C p (λ, β) = 0.73 151 where: 1 This power coefficient is used for a wind generation system with a blade diameter of 126 m and a power of 5 MW.The results shown are for a simulation.
3.3.7.Bustos, Vargas, Milla, Saez, Zareipour, and Nuñez [51] In [51], the authors present a model comparison of a fixed-speed wind turbine operating on a real wind farm.By relying on real data obtained from a wind farm operating in the Chilean interconnected system, three different models are identified and analyzed.The power coefficient used depends on both λ and β in degrees and is expressed by: C p (λ, β) = 0.44 124.99 where: 1 This power coefficient is used for a wind generation system with a blade diameter of 126 m and a power of 2 MW.This paper presents experimental results.
Energies 2023, 16, 2774 14 of 23 3.3.8.Ahmed, Karim, and Ahmad [52] In [52], the authors propose an efficient and low-speed gearless wind-based microgeneration system.This system exploits wind energy in urban and rural parts of the country, which is converted to electrical energy to meet daily domestic energy requirements.The power coefficient used depends on both λ and β and can be expressed as: where: 1 This power coefficient is used for a wind generation system with a blade diameter of 4 m and a power of 1.5 kW.This paper presents simulation results.

General Exponential Function
All the power coefficients registered in an exponential function are grouped for general equation purposes based on an exponential function.This general equation is expressed by ( 42) and ( 43) and the constants of each exponential function are shown in Table 3.
where: 1 Figure 11a shows the response of an exponential function with eight power coefficients and with β equal to zero, for the purpose of comparing its response with the response of the power coefficients of a polynomial function.These coefficients are identified according to the first author of the paper where the power coefficient is presented.The graph indicates the maximum value of the power coefficient (C p max ) and λ (λ max ) for each function.It can be seen that none of the coefficients exceed the Betz limit.All the power coefficients have adequate compression because it starts from zero and increases its value until reaching a maximum value and then decreases to zero as λ increases.The torque coefficient is obtained by applying (12).Therefore in Figure 11b, the response of the power coefficient in relation to λ is shown, based on exponential functions of the power coefficient.
To further evaluate the power coefficients based on exponential functions, an analysis is carried out with variations of the specific velocity and β.
To further evaluate the power coefficients based on exponential functions, an analysis is carried out with variations of the specific velocity and β.
The response of the C p and C t presented by Kotti, Janakiraman, and Shireen [30,31] with variations of λ and β is shown in Figure 12, where it can be observed that the response of C p and C t is adequate because, as β increases, the coefficients decrease, so the control can be implemented through β.
To further evaluate the power coefficients based on exponential functions, an analysis is carried out with variations of the specific velocity and β.
The response of the   and   presented by Kotti, Janakiraman, and Shireen [30,31] with variations of  and  is shown in Figure 12, where it can be observed that the response of   and   is adequate because, as β increases, the coefficients decrease, so the control can be implemented through β.The response of the   and   proposed by Khajuria and Kaur [32] with variations of  and  is shown in Figure 13.Here, the response of   and   is adequate because, for high β values, the coefficients decrease; therefore, a control through β can be implemented with this proposed coefficient.The response of the C p and C t proposed by Khajuria and Kaur [32] with variations of λ and β is shown in Figure 13.Here, the response of C p and C t is adequate because, for high β values, the coefficients decrease; therefore, a control through β can be implemented with this proposed coefficient.The behavior of the curves of   and   proposed by Ovando, Aguayo, and Cotorogea [33][34][35][36][37][38][39][40][41] with variations of  and  is shown in Figure 14.Here, it can be observed that the response of   is adequate because, as β increases, the coefficient decreases; however, this is not the case with   , where, at low and high β values, an undesired behavior is observed, which does not allow for control techniques to be applied through β.The behavior of curves of the   and   proposed by Feng Gao, Da-Ping Xui, and Yue-Gang Lv [42][43][44][45][46][47] with variations of  and  is shown in Figure 15.Here, it is observed that the response of   is adequate because, as β increases, the coefficient decreases.However, this does not happen with   curves, where, for low β and high β values, an undesired behavior is observed, which does not allow for control techniques to be applied through β.The behavior of the curves of C p and C t proposed by Ovando, Aguayo, and Cotorogea [33][34][35][36][37][38][39][40][41] with variations of λ and β is shown in Figure 14.Here, it can be observed that the response of C p is adequate because, as β increases, the coefficient decreases; however, this is not the case with C t , where, at low and high β values, an undesired behavior is observed, which does not allow for control techniques to be applied through β.The behavior of the curves of   and   proposed by Ovando, Aguayo, and Cotorogea [33][34][35][36][37][38][39][40][41] with variations of  and  is shown in Figure 14.Here, it can be observed that the response of   is adequate because, as β increases, the coefficient decreases; however, this is not the case with   , where, at low and high β values, an undesired behavior is observed, which does not allow for control techniques to be applied through β.The behavior of curves of the   and   proposed by Feng Gao, Da-Ping Xui, and Yue-Gang Lv [42][43][44][45][46][47] with variations of  and  is shown in Figure 15.Here, it is observed that the response of   is adequate because, as β increases, the coefficient decreases.However, this does not happen with   curves, where, for low β and high β values, an undesired behavior is observed, which does not allow for control techniques to be applied through β.The behavior of curves of the C p and C t proposed by Feng Gao, Da-Ping Xui, and Yue-Gang Lv [42][43][44][45][46][47] with variations of λ and β is shown in Figure 15.Here, it is observed that the response of C p is adequate because, as β increases, the coefficient decreases.However, this does not happen with C t curves, where, for low β and high β values, an undesired behavior is observed, which does not allow for control techniques to be applied through β.The responses of   and   proposed by Llano, Tatlow, and McMahon [48] with variations of  and  are shown in Figure 16.It is observed that the behavior of the   curve is adequate because, as β increases, the coefficient decreases.However, with   this does not occur because, at low β and high β values, an undesired behavior appears that prevents the application of control techniques through β.The responses of   and   proposed by Shi, Zhu, Cai, Wang, and Yao [49,50] with variations of  and  are shown in Figure 17.Here, the behavior of the curves   and   is adequate because, at high values of β, the coefficients decrease, so that a control through β can be implemented.The responses of C p and C t proposed by Llano, Tatlow, and McMahon [48] with variations of λ and β are shown in Figure 16.It is observed that the behavior of the C p curve is adequate because, as β increases, the coefficient decreases.However, with C t this does not occur because, at low β and high β values, an undesired behavior appears that prevents the application of control techniques through β.The responses of   and   proposed by Llano, Tatlow, and McMahon [48] with variations of  and  are shown in Figure 16.It is observed that the behavior of the   curve is adequate because, as β increases, the coefficient decreases.However, with   this does not occur because, at low β and high β values, an undesired behavior appears that prevents the application of control techniques through β.The responses of   and   proposed by Shi, Zhu, Cai, Wang, and Yao [49,50] with variations of  and  are shown in Figure 17.Here, the behavior of the curves   and   is adequate because, at high values of β, the coefficients decrease, so that a control through β can be implemented.The responses of C p and C t proposed by Shi, Zhu, Cai, Wang, and Yao [49,50] with variations of λ and β are shown in Figure 17.Here, the behavior of the curves C p and C t is adequate because, at high values of β, the coefficients decrease, so that a control through can be implemented.The responses of   and   proposed by Bustos, Vargas, Milla, Saez, Zareipour, and Nuñez [51] with variations of  and  are shown in Figure 18.Here, it is observed that the behavior of   is adequate because, as β increases, the coefficient decreases.However, this does not occur with   because, for low β and high β values, an undesired behavior is observed, which does not allow for the application of control techniques through β.Finally, the responses of   and   proposed by Ahmed, Karim, and Ahmad [52] with variations of  and  are shown in Figure 19.Here, the behavior of curves   and   is adequate because, for high β values, the coefficients decrease, and a control through β can be implemented.The responses of C p and C t proposed by Bustos, Vargas, Milla, Saez, Zareipour, and Nuñez [51] with variations of λ and β are shown in Figure 18.Here, it is observed that the behavior of C p is adequate because, as β increases, the coefficient decreases.However, this does not occur with C t because, for low β and high β values, an undesired behavior is observed, which does not allow for the application of control techniques through β.The responses of   and   proposed by Bustos, Vargas, Milla, Saez, Zareipour, and Nuñez [51] with variations of  and  are shown in Figure 18.Here, it is observed that the behavior of   is adequate because, as β increases, the coefficient decreases.However, this does not occur with   because, for low β and high β values, an undesired behavior is observed, which does not allow for the application of control techniques through β.Finally, the responses of   and   proposed by Ahmed, Karim, and Ahmad [52] with variations of  and  are shown in Figure 19.Here, the behavior of curves   and   is adequate because, for high β values, the coefficients decrease, and a control through β can be implemented.Finally, the responses of C p and C t proposed by Ahmed, Karim, and Ahmad [52] with variations of λ and β are shown in Figure 19.Here, the behavior of curves C p and C t is adequate because, for high β values, the coefficients decrease, and a control through β can be implemented.

Discussion
The coefficients based on polynomial functions depend only on λ, since they consider a constant β.This is possible since these models are normally used in low-power wind turbines, where control of β is not applied.This is unlike the coefficients based on sinusoidal and exponential functions, in which the coefficients depend on the blade pitch angle (β) and the specific velocity (λ).
The analysis was carried out considering two cases.In the first case, only λ is varied, in order to evaluate the power coefficients based on the three mathematical functions.In the second case, both λ and β are varied, but they are only performed for the power coefficients based on sinusoidal and exponential mathematical functions.

Discussion Considering Variations of λ and β = 0
Regarding the coefficients based on polynomial functions, it is observed that they do not follow a specific pattern.Instead, each one varies in a different way, due to which it is necessary in some cases to limit λ so that the coefficient is viable and valid.This is the case for the power coefficient proposed by Arifujjaman, Iqbal, and Quaicoe [21]; here, for values of λ greater than 12.16, the value should be limited, which does not occur in a real turbine.With regard to the torque coefficient, for the case of the coefficient based on a function of the fourth and fifth orders, the value tends to infinity for small values of λ.From the 4 analyzed coefficients based on a polynomial function, the one with the highest   value is the coefficient based on a third-order function proposed by Moussa, Bouallegue, and Kehedher [19], where    is 0.5218 with β equal to 9.892.
For the power coefficients based on sinusoidal functions, it is observed that all the power coefficients have the same behavior and only vary in amplitude and when they start and end.Another disadvantage is that the power coefficient graph does not describe the real behavior of the wind turbine, so it is necessary to limit the value of λ for each power coefficient.This need to limit λ is observed more when the torque coefficient is analyzed, because, for some cases, the coefficient tends to infinity for small values of λ.From the 5 analyzed coefficients, which were based on a sinusoidal function, the one with the highest value of   is that proposed by de Merahi, Mekhilef, and Madjid [29], since a    of 0.5397 is obtained at a β of 4.789; this is the power coefficient with highest maximum value (   ) of the 3 types of power coefficient analyzed.
Regarding the coefficients based on exponential functions, significant variety is observed, but all start at zero for a β of zero, which does not occur in most cases of the

Discussion
The coefficients based on polynomial functions depend only on λ, since they consider a constant β.This is possible since these models are normally used in low-power wind turbines, where control of β is not applied.This is unlike the coefficients based on sinusoidal and exponential functions, in which the coefficients depend on the blade pitch angle (β) and the specific velocity (λ).
The analysis was carried out considering two cases.In the first case, only λ is varied, in order to evaluate the power coefficients based on the three mathematical functions.In the second case, both λ and β are varied, but they are only performed for the power coefficients based on sinusoidal and exponential mathematical functions.

Discussion Considering Variations of λ and β = 0
Regarding the coefficients based on polynomial functions, it is observed that they do not follow a specific pattern.Instead, each one varies in a different way, due to which it is necessary in some cases to limit λ so that the coefficient is viable and valid.This is the case for the power coefficient proposed by Arifujjaman, Iqbal, and Quaicoe [21]; here, for values of λ greater than 12.16, the value should be limited, which does not occur in a real turbine.With regard to the torque coefficient, for the case of the coefficient based on a function of the fourth and fifth orders, the value tends to infinity for small values of λ.From the 4 analyzed coefficients based on a polynomial function, the one with the highest C p value is the coefficient based on a third-order function proposed by Moussa, Bouallegue, and Kehedher [19], where C p max is 0.5218 with β equal to 9.892.
For the power coefficients based on sinusoidal functions, it is observed that all the power coefficients have the same behavior and only vary in amplitude and when they start and end.Another disadvantage is that the power coefficient graph does not describe the real behavior of the wind turbine, so it is necessary to limit the value of λ for each power coefficient.This need to limit λ is observed more when the torque coefficient is analyzed, because, for some cases, the coefficient tends to infinity for small values of λ.From the 5 analyzed coefficients, which were based on a sinusoidal function, the one with the highest value of C p is that proposed by de Merahi, Mekhilef, and Madjid [29], since a C p max of 0.5397 is obtained at a β of 4.789; this is the power coefficient with highest maximum value (C p max ) of the 3 types of power coefficient analyzed.
Regarding the coefficients based on exponential functions, significant variety is observed, but all start at zero for a β of zero, which does not occur in most cases of the coefficients based on polynomial and sinusoidal functions.Moreover, they all reach a maximum value and then descend until reaching zero at different values of β, thus achieving the real behavior of a turbine.In addition, the behavior observed for the coefficient of torque is also very close to reality because almost all start at zero, except for the one proposed by Ovando, Aguayo and Cotorogea [33][34][35][36][37][38][39][40][41], which starts at a value of 0.008, but none of them tends towards infinity for any value of λ.Of the 4 analyzed coefficients based on an exponential function, the one with the highest C p value is the coefficient based on a third-order function proposed by Ovando, Aguayo, and Cotorogea [33][34][35][36][37][38][39][40][41], where C p max is 0.48 at a β of 8.144.

Discussion Considering Both Variations of λ and β
For the cases where variations in λ and β are analyzed, it is observed that the point where the value of the coefficients is higher occurs when β is equal to zero, except for the power coefficients proposed by Moussa, Bouallegue, and Kehedher [19] and Nouira and Khedher [13], which are based on a sinusoidal function, where the maximum value of C p is observed at other values different from β = 0. Additionally, the coefficient proposed by Moussa, Bouallegue, and Kehedher [19] has a problem for high values of both λ and β, because, once the value of C p decreases, it increases again, which does not happen with real wind turbines.
The coefficients based on a single sinusoidal function that presents an adequate behavior for both the C p and the C t is that proposed by Xin, Wanli, Bin, and Pengcheng [28], because, as the value of β increases, the values of C p and C t decrease; it is therefore suitable to apply a control based on β.For the coefficient proposed by Merahi, Mekhilef, and Madjid [29], the behavior of the power coefficient is adequate; however, the power coefficient tends to infinity through variations of λ and β.
In relation to the coefficients based on exponential functions, those that present real behaviors are those proposed by Kotti, Janakiraman, and Shireen [30], Khajuria and Kaur [32], Shi, Zhu, Cai, Wang, and Yao [49,50], and Ahmed, Karim, and Ahmad [42].This is because they are the ones where the values of beta C p and C t increase, meaning that these systems can apply control techniques based on β.The remaining coefficients based on a sine function show the real behavior of C p ; however, the behavior of C t does not correspond to a real turbine, because the value of C t begins to decrease as β begins to increase, but, as β continues to increase, the value of C t also increases, which is not real and therefore cannot use control techniques based on β.

Conclusions
From the literature review carried out with the objective of analyzing the behavior of the power and torque coefficient proposed by different authors, it was found that they are expressed according to mathematical functions based on polynomial, sinusoidal, and exponential equations.These are obtained according to the physical characteristics of the turbines, such as the number of blades, construction materials, the shape of the blades, and β.Therefore, from the mathematical analysis, these coefficients can be seen to essentially depend on three mathematical functions based on each of them in a specific function (polynomial, sinusoidal, or exponential).In this paper, they were classified according to these functions; in the case of coefficients based on sinusoidal and exponential functions, a general equation was obtained.In relation to the general equations obtained, Tables 2 and 3 were established; in these tables, the constants for each analyzed coefficient are established and grouped according to their mathematical model.The coefficients based on polynomial functions only depend on λ, and its operating range must be adequately limited in order to avoid generating non-real turbine behaviors.For example, for low values of the specific velocity, coefficient values are obtained of power that tends to infinity.The coefficients based on sinusoidal and exponential functions depend on λ and β; in this case, it is observed that the maximum value of the power coefficient is obtained when β is zero.In addition, it should be noted that, as β increases, the values of the power and torque coefficient must decrease in order to apply active control techniques of the turbines based on β.However, not all of the coefficients that have been proposed have this characteristic, so it is desirable that their behavior be adequately evaluated.This review offers an overview of all the power coefficients presented in the literature that serve to develop emulators of wind generation systems, develop and evaluate search algorithms of the maximum power point, and have active control techniques by varying β, among other applications.

Figure 1 .
Figure 1.A power coefficient that depends on  and the type of wind turbine[17].

Figure 1 . 23 Figure 2 .
Figure1.A power coefficient that depends on λ and the type of wind turbine[17].

Figure 3 .
Figure3.Torque coefficient that depends on λ and the type of wind turbine[17].

Figure 2 .
Figure 2. A power coefficient for a specific type of turbine, with variations of λ and β [17].

Energies 2023 , 23 Figure 2 .
Figure 2. A power coefficient for a specific type of turbine, with variations of  and β [17].

Figure 3 .
Figure3.Torque coefficient that depends on λ and the type of wind turbine[17].

Figure 4 .
Figure 4. Torque coefficient for a specific type of turbine, with variations of λ and β.

Table 1 .
Constants of the power coefficients found in relation to the order of the general polynomial function.

Table 2 .
Constants of the power coefficients found in relation to the general sinusoidal function.

Table 3 .
Constants of the power coefficients found in relation to the general exponential function.