Dynamic Modeling of Wind Turbines Based on Estimated Wind Speed under Turbulent Conditions

: Large-scale wind turbines with a large blade radius rotates under ﬂuctuating conditions depending on the blade position. The wind speed is maximum in the highest point when the blade in the upward position and minimum in the lowest point when the blade in the downward position. The spatial distribution of wind speed, which is known as the wind shear, leads to periodic ﬂuctuations in the turbine rotor, which causes ﬂuctuations in the generator output voltage and power. In addition, the turbine torque is a ﬀ ected by other factors such as tower shadow and turbine inertia. The space between the blade and tower, the tower diameter, and the blade diameter are very critical design factors that should be considered to reduce the output power ﬂuctuations of a wind turbine generator. To model realistic characteristics while considering the critical factors of a wind turbine system, a wind turbine model is implemented using a squirrel-cage induction motor. Since the wind speed is the most important factor in modeling the aerodynamics of wind turbine, an accurate measurement or estimation is essential to have a valid model. This paper estimates the average wind speed, instead of measuring, from the generator power and rotating speed and models the turbine’s aerodynamics, including tower shadow and wind shear components, without having to measure the wind speed at any height. The proposed algorithm overcomes the errors of measuring wind speed in single or multiple locations by estimating the wind speed with estimation error less than 2%.


Introduction
Currently, electrical power is the main energy source in homes, industries, and workplaces. Industrial development and population growth have led to increased power consumption in the last three decades. Natural resources such as gas, petroleum, and coal are quickly depleting because of their growing use in power plants, vehicles, and industries. This serious problem motivates countries to find alternative energy sources through inexhaustible natural resources. One solution is found in clean, sustainable, inexhaustible, environmentally friendly, and renewable sources such as solar energy, wind energy, and tidal ocean currents. Wind energy has attracted more interest for several reasons and has increasingly become the most harvested renewable energy source [1]. Increased interest in wind turbines to generate electrical power entails studying and modeling the steady state and dynamic behavior of the wind turbine in laboratory conditions to prevent possible problems during installation and later use. In addition, over the last few decades, new control algorithms have been established to improve the performances of wind-generation systems. Before applying a new algorithm, it must be

Wind Turbine Model
The extracted power from the wind at any wind speed and turbine rotational speed is given by [18]  (1) Equation (1) indicates that the power coefficient is influenced by the TSR and the blade pitch angle, which is always below the rated wind speed. In wind energy systems, when the pitch angle is kept constant, the power captured will be a function only in the turbine rotational speed, maximum at a certain speed. Hence, to maximize the captured power, λ should be kept at its optimum value opt λ , and then the captured power is technically computed from the wind speed and the blade radius. When the wind speed is above the rated speed, the pitch angle adjusts to reduce the impact of excessive wind speed on the generated power. Here, β is kept constant for simplicity. Figure 2 presents the turbine blade's power variation with wind speed and rotational speed. The maximum power output occurs at a particular rotational speed.

Wind Shear and Tower Shadow Model
Different wind speeds at different altitudes and the passing of the blade in front of the tower result in periodic torque pulsations in wind turbines.

Wind Turbine Model
The extracted power from the wind at any wind speed and turbine rotational speed is given by [18] Equation (1) indicates that the power coefficient is influenced by the TSR and the blade pitch angle, which is always below the rated wind speed. In wind energy systems, when the pitch angle is kept constant, the power captured will be a function only in the turbine rotational speed, maximum at a certain speed. Hence, to maximize the captured power, λ should be kept at its optimum value λ opt , and then the captured power is technically computed from the wind speed and the blade radius. When the wind speed is above the rated speed, the pitch angle adjusts to reduce the impact of excessive wind speed on the generated power. Here, β is kept constant for simplicity. Figure 2 presents the turbine blade's power variation with wind speed and rotational speed. The maximum power output occurs at a particular rotational speed.

Wind Turbine Model
The extracted power from the wind at any wind speed and turbine rotational speed is given by [18]  (1) Equation (1) indicates that the power coefficient is influenced by the TSR and the blade pitch angle, which is always below the rated wind speed. In wind energy systems, when the pitch angle is kept constant, the power captured will be a function only in the turbine rotational speed, maximum at a certain speed. Hence, to maximize the captured power, λ should be kept at its optimum value opt λ , and then the captured power is technically computed from the wind speed and the blade radius. When the wind speed is above the rated speed, the pitch angle adjusts to reduce the impact of excessive wind speed on the generated power. Here, β is kept constant for simplicity. Figure 2 presents the turbine blade's power variation with wind speed and rotational speed. The maximum power output occurs at a particular rotational speed.

Wind Shear and Tower Shadow Model
Different wind speeds at different altitudes and the passing of the blade in front of the tower result in periodic torque pulsations in wind turbines.

Wind Shear and Tower Shadow Model
Different wind speeds at different altitudes and the passing of the blade in front of the tower result in periodic torque pulsations in wind turbines. Based on the height, the wind speed changes as follows [19,20]: (2) Table 1 summarizes the geographical features and represents a general wind shear index [1]. As a result of height-dependent wind speeds, the turbine blades oscillate in response to the wind shear. Additionally, the amplitude of the wind speed increases with the height for different rotating angles from 0 • to 360 • . The wind speed in Equation (2) can be rewritten with respect to the speed at the turbine hub as follows [21]: Considering the blade's rotation and the unique wind shear at each angle, instead of expressing the wind shear in only the z direction, the wind shear in Equation (6) can be expressed as a function of the blade rotational angle ψ b and radial distance r from the rotor axis. The wind shear torque can be modeled as [22]: which can be written as Equation (5) can be represented as a function of the wind shear disturbance F(r, ψ b ): Due to wind shear, the wind speed can be expanded using Fourier transform as shown in the following Equation [22]: In low altitude, wind shear is influenced by the large amount of change in wind speed due to the friction of the surface. In addition, the longer the blade length of the wind turbine, the wider the turning radius of the rotor, and the larger the wind shear effect is. Figure 3 shows the wind shear effect of the turbine blade length, r hub height, and H empirical wind shear index as in Equation (7). In low altitude, wind shear is influenced by the large amount of change in wind speed due to the friction of the surface. In addition, the longer the blade length of the wind turbine, the wider the turning radius of the rotor, and the larger the wind shear effect is. Figure 3 shows the wind shear effect of the turbine blade length, r hub height, and H empirical wind shear index as in Equation (7). In the case of a large-diameter wind turbine, the blades experience a wide range of wind speeds in each revolution. With three-blade wind turbines, the frequent rotation through this wide range of speeds results in torque and power oscillations at three times the rotor rotational speed, called 3p frequency. In the case of a large-diameter wind turbine, the blades experience a wide range of wind speeds in each revolution. With three-blade wind turbines, the frequent rotation through this wide range of speeds results in torque and power oscillations at three times the rotor rotational speed, called 3p frequency.
The tower radius (a), the distance of the blade origin from the tower midline (x), and the lateral distance of the blade from the tower midline (y) affect the tower shadow amplitude, as shown in Figure 4. The tower radius (a), the distance of the blade origin from the tower midline ( ), and the lateral distance of the blade from the tower midline ( ) affect the tower shadow amplitude, as shown in Figure 4. In front of the tower, the disturbance is nonnegligible and affects the wind speed and turbine output power. The wind speed reduces by a value of Δ in front of the tower and within the region W of 3-6 times of the tower diameter, as shown in Figure 5. The tower shadow as a function of x and y is given by [23] ( ) Expressing the tower shadow as a function of r and b ψ , In front of the tower, the disturbance is nonnegligible and affects the wind speed and turbine output power. The wind speed reduces by a value of ∆ in front of the tower and within the region W of 3-6 times of the tower diameter, as shown in Figure 5. The tower radius (a), the distance of the blade origin from the tower midline ( ), and the lateral distance of the blade from the tower midline ( ) affect the tower shadow amplitude, as shown in Figure 4. In front of the tower, the disturbance is nonnegligible and affects the wind speed and turbine output power. The wind speed reduces by a value of Δ in front of the tower and within the region W of 3-6 times of the tower diameter, as shown in Figure 5. The tower shadow as a function of x and y is given by [23] ( ) Expressing the tower shadow as a function of r and b ψ , The tower shadow as a function of x and y is given by [23] V(x, y) Expressing the tower shadow as a function of r and ψ b , where m = 1 + α(α−1)R 2 8H 2 . Figure 6 shows wind shear and tower shadow effects at the tip of each blade at an average wind speed of 8 m/s. The tower shadow shape was calculated for different radial circulating points, from 8 m to 20 m in diameter, measured from the hub center, with an incremental distance of 8 m, as shown in Figure 6. When the diameter increases, the points on this circle perimeter experience short tower shadow intervals. The inner curve in Figure 6 indicates the circle with a diameter of 20 m, which has the lowest tower shadow interval. In contrast, the outer and widest curve in Figure 7 illustrates the shadow where 2 8 H . Figure 6 shows wind shear and tower shadow effects at the tip of each blade at an average wind speed of 8 m/s. The tower shadow shape was calculated for different radial circulating points, from 8 m to 20 m in diameter, measured from the hub center, with an incremental distance of 8 m, as shown in Figure 6. When the diameter increases, the points on this circle perimeter experience short tower shadow intervals. The inner curve in Figure 6 indicates the circle with a diameter of 20 m, which has the lowest tower shadow interval. In contrast, the outer and widest curve in Figure   If the radial distance r is kept constant and the distance between the blade origin to the tower midline x varies from 3 m to 5 m, as Figure 8 illustrates, the tower shadow effect increases when the distance x is small. Therefore, the blades should be installed far from the tower to reduce the tower shadow effect. The tower shadow is given by   Figure 6 shows wind shear and tower shadow effects at the tip of each blade at an average wind speed of 8 m/s. The tower shadow shape was calculated for different radial circulating points, from 8 m to 20 m in diameter, measured from the hub center, with an incremental distance of 8 m, as shown in Figure 6. When the diameter increases, the points on this circle perimeter experience short tower shadow intervals. The inner curve in Figure 6 indicates the circle with a diameter of 20 m, which has the lowest tower shadow interval. In contrast, the outer and widest curve in Figure 7   If the radial distance r is kept constant and the distance between the blade origin to the tower midline x varies from 3 m to 5 m, as Figure 8 illustrates, the tower shadow effect increases when the distance x is small. Therefore, the blades should be installed far from the tower to reduce the tower shadow effect. The tower shadow is given by If the radial distance r is kept constant and the distance between the blade origin to the tower midline x varies from 3 m to 5 m, as Figure 8 illustrates, the tower shadow effect increases when the distance x is small. Therefore, the blades should be installed far from the tower to reduce the tower shadow effect. The tower shadow is given by For any wind speed value, the turbine torque T t can be expressed as: The wind shear module implements the disturbance in torque due to wind shear given as Rv ws (13) and the tower shadow module implements the disturbance in torque due to the tower shadow expressed as Due to turbine torque, wind shear, and tower shadow, the resultant aerodynamic torque is defined as Equation (15a). For any wind speed value, the turbine torque Tt can be expressed as: The wind shear module implements the disturbance in torque due to wind shear given as (13) and the tower shadow module implements the disturbance in torque due to the tower shadow expressed as Due to turbine torque, wind shear, and tower shadow, the resultant aerodynamic torque is defined as Equation (15a).
However, the wind speed at the hub level is not accurate when modeling the turbine torques. Hence, a wind speed estimation method can be used to estimate the average wind speed to resemble real characteristics. The resultant torque in Equation (15a) can be expressed as Equation (15b):

Dynamic of the Rotating Masses in the Wind Turbine
The wind turbine rotor can be coupled directly to the generator shaft, or through a gearbox. The model traditionally used to represent the dynamic behavior of the wind group is the one that uses a system composed of two masses and their respective moments of inertia, with one mass However, the wind speed at the hub level is not accurate when modeling the turbine torques. Hence, a wind speed estimation method can be used to estimate the average wind speed to resemble real characteristics. The resultant torque in Equation (15a) can be expressed as Equation (15b):

Dynamic of the Rotating Masses in the Wind Turbine
The wind turbine rotor can be coupled directly to the generator shaft, or through a gearbox. The model traditionally used to represent the dynamic behavior of the wind group is the one that uses a system composed of two masses and their respective moments of inertia, with one mass representing the complete wind turbine and the other the generator rotor. These masses are connected by an axis that has its stiffness and damping modeled, as shown in Figure 9. A gearbox with an n ratio is used to couple the low-speed side of the turbine to the high-speed side of the generator. representing the complete wind turbine and the other the generator rotor. These masses are connected by an axis that has its stiffness and damping modeled, as shown in Figure 9. A gearbox with an n ratio is used to couple the low-speed side of the turbine to the high-speed side of the generator. Figure 9. Dynamic model of a wind turbine.

The Wind Speed Simulator
To model the torque fluctuation of a running wind turbine, a stochastic model for wind speed is required. This model is a sum of harmonics and is expressed as [24]    For each selected harmonic component, the amplitude can be written as a function of Dryden

The Wind Speed Simulator
To model the torque fluctuation of a running wind turbine, a stochastic model for wind speed is required. This model is a sum of harmonics and is expressed as [24]  For each selected harmonic component, the amplitude can be written as a function of Dryden spectra φ u as [25] The longitudinal Dryden power spectrum is defined in terms of turbulence intensity σ, turbulence length, and harmonic frequency as follows [26]:

Induction Motor Control
The induction machine can be modeled by writing the dynamic equations in the d-q synchronous reference frame form as [27] v ds = i ds R s − ω e λ qs + pλ ds (19) v qs = i qs R s + ω e λ ds + pλ qs (20) The electromagnetic torque in the induction machine is fundamentally defined as The motor torque and flux connected to the power source are controlled through the back-to-back converter. Whereas the motor flux is kept constant, the torque is controlled to provide the wind turbine characteristics. To generate a continuous variable torque, the torque current i qs is continuously controlled as defined in Equation (23). In experiment, the turbine reference torque and wind speed pattern are implemented by Equations (15b) and (16), respectively. The torque current reference i * qs is then calculated from the torque equation by dividing it by the torque constant K t , as shown in Figure 10a.
The electromagnetic torque in the induction machine is fundamentally defined as The motor torque and flux connected to the power source are controlled through the back-to-back converter. Whereas the motor flux is kept constant, the torque is controlled to provide the wind turbine characteristics. To generate a continuous variable torque, the torque current qs i is continuously controlled as defined in Equation (23). In experiment, the turbine reference torque and wind speed pattern are implemented by Equations (15b) and (16), respectively. The torque current reference * qs i is then calculated from the torque equation by dividing it by the torque constant t K , as shown in Figure 10a.
In fact, the motor reference d-axis current is kept constant at the value that maintains the minimum iron loss. In contrast, in the grid-side converter, the DC-link voltage is kept constant, and the output power factor is adjusted to unity by controlling the grid-side converter q-axis current and the d-axis current. The dc-link voltage is controlled to the desired value by using a proportional-integral (PI) controller and the change in the dc-link voltage represents a change in the q-axis current. A current feed-forward control loop is also used here to improve the dc-link voltage response to load disturbance. For unity power factor, the demand for the d-axis current is zero. Figure 10b shows a control block diagram of the grid-side converter.

Control Scheme of DFIG
To control the active and reactive power in the stator side, a d-q reference frame synchronization with the stator flux is selected. The flux vector in the stator is aligned with the d-axis. The flux linkage of the rotor and stator is defined as [28] Rotor voltages in the d-q reference frame can be written in terms of rotor and stator magnetizing currents.
The stator flux angle is expressed as In fact, the motor reference d-axis current is kept constant at the value that maintains the minimum iron loss. In contrast, in the grid-side converter, the DC-link voltage is kept constant, and the output power factor is adjusted to unity by controlling the grid-side converter q-axis current and the d-axis current. The dc-link voltage is controlled to the desired value by using a proportional-integral (PI) controller and the change in the dc-link voltage represents a change in the q-axis current. A current feed-forward control loop is also used here to improve the dc-link voltage response to load disturbance. For unity power factor, the demand for the d-axis current is zero. Figure 10b shows a control block diagram of the grid-side converter.

Control Scheme of DFIG
To control the active and reactive power in the stator side, a d-q reference frame synchronization with the stator flux is selected. The flux vector in the stator is aligned with the d-axis. The flux linkage of the rotor and stator is defined as [28] λ s = λ ds = L m i ms = L s i ds + L m i dr (24) λ dr = L 2 m L s i ms + σL r i dr (25) λ qr = σL r i dr (26) Rotor voltages in the d-q reference frame can be written in terms of rotor and stator magnetizing currents.
The stator flux angle is expressed as λ s qs = (v s qs − R s i s qs )dt (31) θ e = tan −1 λ s qs λ s ds (32) where the superscript s represents the parameters in the stationary reference frame.
Regulating the q-axis component of the rotor current controls either the active power of the DFIG stator side or the developed torque in the generator.
To directly control the stator-side reactive power, the rotor d-axis current component should be regulated. This reactive power can be expressed as The stator active and reactive powers are proportional to the current components i qr , and i dr , respectively.
The DFIG wind turbine system configuration consists of a DFIG connected to the utility grid through its stator. The back-to-back Pulse-Width Modulation (PWM) power converter enables a bidirectional power-flow control. As a result, the DFIG can operate either in sub-synchronous (ω r < ω e ) or in super-synchronous modes (ω r > ω e ). By calculating the stator active power, the control block diagram of the stator power is presented in the figure. The PI controller gains for the active and reactive power controllers are K pp , K ip , K pq , and K iq , respectively. The rotor reference q-axis current i * qr is then calculated as the output of the active power controller and input for the inner current control loop. The rotor instantaneous q-axis current is then calculated from the sensed three-phase rotor currents and controlled to produce a reference q-axis rotor voltage. Normally, the output reactive power of the wind power conversion system is controlled as zero to keep unity power factor of the stator voltage and current. The stator reactive power is controlled to the desired value to produce the reference d-axis rotor voltage as shown in Figure 11b. One can see how the outer stator power feedback loop produces the rotor reference d-axis current i * dr for the inner current feedback control loop. The reference q-axis rotor voltage is then produced by controlling the rotor d-axis component. The schematic diagram of the simplified control scheme is shown in Figure 11c.
The DFIG optimum power P * is the reference power value determined by the wind speed and rotor angular speed of the power controller loop. A simple proportional-integral (PI) controller can be utilized to regulate the d-q component of the rotor current only if the reference frames for both the measured and reference current vectors match [29].
To achieve the full control of the grid-side current, the dc-link voltage must be boosted to a level higher than the amplitude of the line-line voltage. The power flow of the grid-side converter is controlled so as to keep the dc-link voltage constant. To maintain the dc-link voltage constant and to ensure the reactive power flowing into the grid at null, the grid-side converter currents are controlled using the d-q vector control approach. The dc-link voltage is controlled to the desired value by using an IP controller and the change in the dc-link voltage represents a change in the q-axis current. A simple voltage control scheme is shown in Figure 12a. A current feed-forward control loop is also used here to improve the dc-link voltage response to load disturbance. For unity power factor, the demand for the d-axis current is zero. Figure 12b   To achieve the full control of the grid-side current, the dc-link voltage must be boosted to a level higher than the amplitude of the line-line voltage. The power flow of the grid-side converter is controlled so as to keep the dc-link voltage constant. To maintain the dc-link voltage constant and to ensure the reactive power flowing into the grid at null, the grid-side converter currents are controlled using the d-q vector control approach. The dc-link voltage is controlled to the desired value by using an IP controller and the change in the dc-link voltage represents a change in the q-axis current. A simple voltage control scheme is shown in Figure 12a. A current feed-forward control loop is also used here to improve the dc-link voltage response to load disturbance. For unity  power factor, the demand for the d-axis current is zero. Figure 12b shows a control block diagram of the grid-side converter.

Support Vector Regression
The SVR is an algorithm that estimates a function that predicts unknown mapping between a system's l-dimensional inputs and real output using a series of training data. When the map is accurately constructed, the relation between the inputs and output is then used to predict the target output. The regression algorithm approximates the unknown desired function in the following form [30]: These values are extracted from the available training data to determine the target function. The first step is to obtain both w and b to minimize the real data differences outside the ε -insensitive region, which are determined by slack variables ξ and * ξ . To minimize the empirical risk, the slack variables are applied as [31] ( ) ( )

Support Vector Regression
The SVR is an algorithm that estimates a function that predicts unknown mapping between a system's l-dimensional inputs and real output using a series of training data. When the map is accurately constructed, the relation between the inputs and output is then used to predict the target output. The regression algorithm approximates the unknown desired function in the following form [30]: These values are extracted from the available training data to determine the target function. The first step is to obtain both w and b to minimize the real data differences outside the ε-insensitive region, which are determined by slack variables ξ and ξ * . To minimize the empirical risk, the slack variables are applied as [31] R reg ( f ) = 1 2 subject to The penalty factor C determines the compromise value between reduction to minimize both empirical and global errors and the model complexity term w 2 . SVR reduces the error to zero if the penalty factor is very high. This option leads to a complex model. However, a low C leads to high estimation errors.
The regression problem in Equation (24) is then solved using dual multipliers: is the kernel function that enables a dot product in high-dimensional feature space using low-dimensional space data input without calculating the function. The radial base function that is used in the framework is expressed as shown in the following equation:

Particle Swarm Optimization
The PSO method is a population-based search technique that starts with a population of random solutions called particles [32]. Every solution within the swarm is called a particle [33]. The swarm is arbitrarily initialized and updated in each reiteration t to fully adapt the fitness function. The i th d-dimensional particle is characterized by its location vector P t i = P t i1 , P t i2 , . . . , P t id and its speed vector Each particle knows its best personal location P t i,Best and the entire population's best global solution G t Best . The position of the population at iteration t is P t = P t 1 , P t 2 , , . . . , P t N , where N is the population size. Each particle updates its location according to Equations (40) and (42): The constriction coefficients introduced in [33,34] are used to set c 1 and c 2 : where φ = φ 1 + φ 2 > 4. In this paper, φ 1 = φ 2 = 2.05; thus,c 1 = c 2 = 1.4962. Location and speed boundaries are adjusted such that if a particle tries to pass the allowable boundaries, a limiting procedure brings it back within the permissible limit. If the best personal solution P i,t Best has a higher fitness than the current best global solution G t−1 Best , then G t Best is equal to P i,t Best [35]. The end points are determined when the best global solution results in a permissible fitness or when a predetermined number of iterations is achieved.
The wind speed estimation using PSO-SVR method can be described in steps as follows: • Enter original data for estimation, known as the data preparation step; • Initialize the particles with random velocities and positions, known as particles optimization; • Perform an offline training process for SVR with training samples assessing each particle fitness value of the PSO for the SVR; • Update each particle velocity and position until the termination condition is satisfied; • Construct and retrain the SVR estimation model based on the optimal parameters.
Training data are obtained from the turbine power equation with pre-specified rotor speed and wind velocity samples. During on-line operation, the turbine power is calculated as For each sample, the rotor speed and the corresponding turbine power ( ω i , p i,j i = 1, . . . 50, j = 1, . . . 50 ) are combined as a pair for the inputs of SVR estimator as shown in Figure 13a. For any operating condition, if the turbine power and the rotor speed are known, the wind speed can be calculated as described in the framework of a PSO-SVR method in Figure 13b [36].
• Initialize the particles with random velocities and positions, known as particles optimization; • Perform an offline training process for SVR with training samples assessing each particle fitness value of the PSO for the SVR; • Update each particle velocity and position until the termination condition is satisfied; • Construct and retrain the SVR estimation model based on the optimal parameters.
Training data are obtained from the turbine power equation with pre-specified rotor speed and wind velocity samples. During on-line operation, the turbine power is calculated as   The detailed implementation process of the proposed method is shown in Figure 14.
Initialization and parameter PSO fitness calculation Online Implementation The detailed implementation process of the proposed method is shown in Figure 14.
The detailed implementation process of the proposed method is shown in Figure 14. Online Implementation Figure 14. Implementation process diagram.

Experimental Results
It is desirable that the wind turbine simulator is implemented in hardware (M-G set) to test the WECS operation in the laboratory. The hardware setup consists of the DFIG-driven by the SCIM, PWM converters, and digital signal processor (DSP) control boards as shown in Figure 15. The specifications of the wind turbine blade modeled for the simulator, the SCIM motor ratings and parameters, and the experimental setup are listed in Tables A1 and A2 in the Appendix. Figure 16 shows the turbine performance for turbulent wind speed, of which the spectrum is similar to the real wind speed pattern. High estimation accuracy and very high PSO-SVR dynamic performance is illustrated in Figure 16a. It is noticeable that the estimated wind speed has a slight difference from the real value due to the trade-off between the minimizing error and the model complexity. The Figure 14. Implementation process diagram.

Experimental Results
It is desirable that the wind turbine simulator is implemented in hardware (M-G set) to test the WECS operation in the laboratory. The hardware setup consists of the DFIG-driven by the SCIM, PWM converters, and digital signal processor (DSP) control boards as shown in Figure 15. The specifications of the wind turbine blade modeled for the simulator, the SCIM motor ratings and parameters, and the experimental setup are listed in Tables A1 and A2 in the Appendix A. Figure 16 shows the turbine performance for turbulent wind speed, of which the spectrum is similar to the real wind speed pattern. High estimation accuracy and very high PSO-SVR dynamic performance is illustrated in Figure 16a. It is noticeable that the estimated wind speed has a slight difference from the real value due to the trade-off between the minimizing error and the model complexity. The generator rotational speed follows the wind speed pattern in all wind speed values, as shown in Figure 16b. Figure 16c plots the generator output power versus variable wind speed. By adding the effect of wind shear and tower shadow, the total turbine torque oscillates in the same frequency. However, the inertia of the turbine leads to the reduction of the effect of both disturbances as presented in Figure 17a,b. By adding the effect of wind shear and tower shadow, the total turbine torque oscillates in the same frequency. However, the inertia of the turbine leads to the reduction of the effect of both disturbances as presented in Figure 17a,b. It can be observed that periodic components are concentrated on frequency multiples of the blade passing frequency, mainly 3p, 6p … etc., as it is shown in Figure 18. It should be expected that the amplitude decreases when they increase in order.  It can be observed that periodic components are concentrated on frequency multiples of the blade passing frequency, mainly 3p, 6p . . . etc., as it is shown in Figure 18. It should be expected that the amplitude decreases when they increase in order.
Equation (16) is then used to generate a stochastic model of the wind speed to model the real nature of instantaneous wind speed, as shown in Figure 19a. The turbulent wind speed spectrum is generated by the wind speed simulator. The parameters are chosen as follows: V o = 9 [m/s], N = 15 in (11), L = 90 [m] and σ = 0.2 in Equation (13). It is obvious that the tower shadow and the wind shear are quite significant for the variations of the turbine torque and power. The range of harmonic content in the wind speed profile is from 0.1 Hz to 10 Hz. From Figure 19, it can be observed that the turbine performance for turbulent wind speed continuously follows the wind speed pattern, as shown in Figure 19b-e. During this low wind speed, the pitch angle controller keeps the pitch angle position in the minimum value, which is zero. However, when the average wind speed increases to 12 m/s, the generator rotational speed and generator output power increase to the maximum limits, which activate the pitch angle controller to increase the pitch angle to the predetermined position based on the wind speed value. In this case, the generator output power and the generator rotating speed are constant and equal to the rated values. Figure 20 depicts the transition between the maximum power point tracking (MPPT) mode and constant power mode. The generator output power is constant when the speed goes higher than 13 m/s in Figure 20a, as shown in Figure 20c. The generator speed, rotor q-axis current, torque and have the same trend as in Figure 20b, d and e, respectively. It can be observed that periodic components are concentrated on frequency multiples of the blade passing frequency, mainly 3p, 6p … etc., as it is shown in Figure 18. It should be expected that the amplitude decreases when they increase in order.  Figure 19, it can be observed that the turbine performance for turbulent wind speed continuously follows the wind speed pattern, as shown in Figure 19b-e. During this low wind speed, the pitch angle controller keeps the pitch angle position in the minimum value, which is zero. However, when the average wind speed increases to 12 m/s, the generator rotational speed and generator output power increase to the maximum limits, which activate the pitch angle controller to increase the pitch angle to the predetermined position based on the wind speed value. In this case, the generator output power and the generator rotating speed are constant and equal to the rated values. Figure 20 depicts the transition between the maximum power point tracking (MPPT) mode and constant power mode. The generator output power is constant when the speed goes higher than 13 m/s in Figure 20a, as shown in Figure 20c. The generator speed, rotor q-axis current, torque and have the same trend as in Figure 20b, d and e, respectively.     The wind variation is modeled as a sum of wide range harmonics with frequencies 0.1-10 Hz. The spectrum of the turbulent wind speed in Figure 21 shows the low frequency components of the wind speed. The measured stator and rotor currents are shown in Figure 22a,b. From the stator current, it is obvious that the currents are almost sinusoidal.     The wind variation is modeled as a sum of wide range harmonics with frequencies 0.1-10 Hz. The spectrum of the turbulent wind speed in Figure 21 shows the low frequency components of the wind speed. The measured stator and rotor currents are shown in Figure 22a,b. From the stator current, it is obvious that the currents are almost sinusoidal.   The operating modes of a wind turbine are obtained, as in Figure 23. This curve is divided into four regions, as follows: The operating modes of a wind turbine are obtained, as in Figure 23. This curve is divided into four regions, as follows: -Region AB, when the wind speed is less than the cut-in speed when the rotational speed is less than the minimum rotational speed for optimum operation. -Region BC, when the wind speed is higher than the cut-in speed and less than the rated value. The output power is given by P opt = K opt υ 3 . -Region CD, when the rotational speed approaches to its rated value. -Region DE, when the wind speed is beyond the limits and the generator output power is controlled to its rated value. The blade pitch controller is activated in this region.  The operating modes of a wind turbine are obtained, as in Figure 23. This curve is divided into four regions, as follows: -Region AB, when the wind speed is less than the cut-in speed when the rotational speed is less than the minimum rotational speed for optimum operation. -Region BC, when the wind speed is higher than the cut-in speed and less than the rated value. The output power is given by -Region CD, when the rotational speed approaches to its rated value.
-Region DE, when the wind speed is beyond the limits and the generator output power is controlled to its rated value. The blade pitch controller is activated in this region.        Figure 24a,b show the accuracy of the proposed wind estimation algorithm to estimate the wind speed. The estimation error in continuous change of wind speed is in the range of 0.2 m/s, which is considered a small value with respect to using an anemometer or measuring the average wind speed in different locations.

Conclusion
In this paper, a novel sensorless-based modeling for wind energy system is presented. It is developed using a torque-controlled SCIM as the wind turbine simulator. Stochastic wind speed profiles, wind shear, and torque shadow effects are simulated: 1. The turbine steady-state characteristics are modeled using the turbine torque equation, which is a function of the wind speed, tip-speed ratio, and blade pitch angle.

Conclusions
In this paper, a novel sensorless-based modeling for wind energy system is presented. It is developed using a torque-controlled SCIM as the wind turbine simulator. Stochastic wind speed profiles, wind shear, and torque shadow effects are simulated: 1.
The turbine steady-state characteristics are modeled using the turbine torque equation, which is a function of the wind speed, tip-speed ratio, and blade pitch angle.

2.
A derivation for the wind shear and tower shadow torque components are obtained as a function of the turbine dimensions and the estimated effective wind speed.
The main cause of system oscillations is tower shadow, wind shear shows significant oscillations on the total output torque, and the spectrum analysis shows the oscillation frequencies, which are from the third-order components. The electromagnetic transient models of DFIG and control strategies were presented and analyzed.
All the characteristics are simulated based on the estimated wind speed using the PSO-SVR method, which gives an accurate simulation of the dynamic characteristics of wind turbine.

Conflicts of Interest:
The authors declare no conflict of interest. wind speed at hub height ψ b blade rotational angle r radial distance from the rotor axis F(r, ψ b ) the wind shear disturbance a the tower radius x the distance from the blade origin to the tower midline y the lateral distance of the blade from the tower midline V ts the tower shadow T t turbine torque T ts disturbance in torque due to tower shadow T ws disturbance in torque due to wind shear J T moment of inertia of the wind turbine J G moment of inertia of the generator k stiffness of coupling shaft D coupling shaft damping n transmission ratio of the gearbox P g the generator power ω m the turbine speed J the system moment of inertia υ(t) instantaneous wind speed V o effective average value of wind speed N harmonic samples order ω i harmonic frequency A i harmonic amplitude σ turbulence intensity L u turbulence length v ds , v qs the d-q voltages in the stator i ds , i qs the d-q currents in the stator ω e the frequency in the rotor λ ds , λ qs the d-q flux linkages in the stator side L m equivalent magnetizing inductance L s equivalent stator self-inductance L r equivalent rotor self-inductance λ ds , λ qs equivalent d-q flux linkage in the stator λ dr , λ qr equivalent rotor d -q flux linkage i ms , i ds , i dr magnetizing d-axis currents in the stator and rotor K t the torque constant P s , Q s Stator active and reactive power v dr , v qr the d-q voltages in the rotor R r rotor equivalent resistance ω sl slip angular frequency R s the stator resistance θ e the synchronous frame angle w the SVR weight vector b the bias ε the allowable error C the penalty factor α i , α * i Lagrange multipliers. c 1 and w t is the velocity weight at iteration t. c 2 social learning factor r 1 , r 2 cognition factor are random numbers in the range [0, 1]

Appendix A
The specification of the induction machine used for the test is three-phase, four poles, 230 V, 50 Hz, 3 kW.