Time-DomainModeling of Tower Shadow andWind Shear in Wind Turbines

Tower shadow and wind shear contribute to periodic fluctuations in electrical power output of a wind turbine generator. The frequency of the periodic fluctuations is n times the blade rotational frequency p, where n is the number of blades. For three-bladed wind turbines, this inherent characteristic is known as the 3p effect. In a weak-power system, it results in voltage fluctuation or flicker at the point of common coupling of the wind turbine to the grid. The phenomenon is important to model so as to evaluate the flicker magnitude at the design level. Hence, the paper aims to develop a detailed time-domain upwind fixed speed wind turbine model which includes the turbine’s aerodynamic, mechanical, electrical, as well as tower shadow and wind shear components. The model allows users to input factors such as terrain, tower height, and tower diameter to calculate the 3p oscillations. The model can be expanded to suit studies involving variable speed wind turbines. Six case studies demonstrate how the model can be used for studying wind turbine interconnection and voltage flicker analysis. Results indicate that the model performs as expected.


Introduction
In addition to the stochastic pulsations caused by wind, electrical power output from wind turbines connected to the power system grid experience short-term fluctuations at n times the rotational frequency of the blades (p) or np frequency, where n is the number of blades.Wind turbines generally have three blades and hence this is known as the 3p effect.The 3p fluctuations are due to two effects known as wind shear and tower shadow [1,2].Wind shear refers to the variation of wind with height while tower shadow is used to describe the redirection in wind due to the obstruction presented by the tower [1,3].These two effects are an inherent characteristic of any wind turbine and will result in periodic fluctuations in output power.
For wind turbines connected to a weak grid, the 3p frequency in the range of 0.5-2 Hz results in fluctuations in system rms voltage at the point of common coupling (PCC) and are known as voltage flicker or fluctuations [4,5].In a system with significant wind generation wherein wind farms supply majority of the power to the grid, synchronization of the blades in the wind farm contribute to larger power and voltage fluctuations [2].From a consumer perspective, voltage flicker is a serious power quality issue [6] and hence wind turbine models have become indispensable for wind power interconnection study.Using such a model, the designer of the wind turbine must be capable of calculating at the design stage, the expected flicker magnitude and whether the flicker is below the threshold of objection [6].The designer must also be able to assess the dependence of voltage flicker on wind speed, grid characteristics, as well as the loads connected at the PCC [4].Many papers on wind power integration have been reported in the literature [4,7,8].Due to their focus on voltage flicker and other power quality issues, these models often tend to oversimplify the torque considerations due to tower shadow and wind shear.They do not factor in wind turbine specifications such as hub height, tower diameter, and terrain on which the wind turbine is located which affect the magnitude of tower shadow and wind shear.These parameters should be an input to the wind turbine simulation model, offering insight into the wind turbine parameters which can be optimized to reduce 3p oscillations at the design level.Other papers oversimplify wind turbine operation and rotor performance characteristics [9].
Based on the aforementioned background, the objective of this paper is to develop a detailed time-domain wind turbine model upon which more sophisticated models can be developed.The paper focuses on modeling a fixedspeed wind turbine with an upwind rotor since they are most common due to reduced tower shading [1,10].The paper builds upon the most basic model of a fixed-speed wind turbine developed in [11] and includes the turbine's aerodynamic, mechanical, electrical, as well as tower shadow and wind shear components.The complete model can be easily implemented in popular simulation packages such as PSCAD/EMTDC and MATLAB.Since the model has been broken down into blocks/modules and each module is separately quantified, the model is flexible, allowing users to separately study the impact of wind shear and tower shadow on output power.The model can also be used to study the wind turbine integration to the grid, voltage flicker analysis, mitigation techniques of 3p effects, and expanded to suit studies involving variable speed wind turbines.
The paper is organized as follows: Sections 2 and 3 study the modeling and identifies the factors affecting wind shear and tower shadow.Section 4 discusses the comprehensive modeling of the wind turbine.Section 5 presents six case studies to evaluate the validity of the model along with the dataset used in simulating the wind turbine so that it can be easily reproduced by the users.The results are consistent with those reported in the literature and hence the model performs as expected.

Wind Shear
Wind turbine blades experience maximum wind speed when facing directly upwards and minimum when the blades are facing downwards.This variation in wind with height is known as wind shear.For three-bladed wind turbines, each of the three blades will experience minimum wind speed in one complete rotation and hence wind shear contributes to the 3p frequency [1].This section studies the wind shear model and identifies the factors which contribute to the wind shear effect.

Modeling Wind Shear.
A simple wind shear model is given by the power law as [1,3] where z, H, V h , and α are the elevation above ground, hub height, wind speed at hub height, and empirical wind shear exponent, respectively.Writing (1) as a function of r (radial distance from rotor axis) and θ (blade azimuthal angle), the wind shear model becomes where W s (r, θ) is the disturbance in wind speed due to wind shear.It can be approximated by a third order truncated Taylor series as, (3)

Factors Affecting
Wind Shear.From ( 2) and ( 3), disturbance in wind due to wind shear is dependent on V h , H, α, and r. Figure 1(a) shows the variation of wind shear with V h .For V h = 13.5 m/s, the maximum wind speed experienced by the blades of a wind turbine is 14.75 m/s while the minimum wind speed is 13.04 m/s.Hence, the variation in wind speed, that is, the difference between maximum and minimum wind speed for one complete blade rotation is 1.71 m/s.Similarly for V h = 14 m/s, the variation in wind speed is 1.77 m/s.Thus wind shear increases with the increase in V h .Also note that the minimum wind speed occurs at θ = 180 • , that is, when the blades are pointing downwards and maximum wind speed at θ = 0 • , that is, blades pointing upwards.Figure 1(b) shows the impact of α on wind shear.The values of α have been chosen depending on terrain, as shown in Table 1 [12].As α increases, wind shear increases.Figure 1(c) shows that wind shear decreases with increase in hub height.For example, at H = 120 m, the variation in wind is 0.72 m/s while at H = 80 m, the wind variation is 1.77 m/s. Figure 1(d) shows the wind shear experienced by different blade elements.Minimum wind shear is experienced by the blade element closest to the hub (r = 5 m).

Tower Shadow
Tower shadow effect is the alteration in uniform flow of wind due to the presence of the tower.For an upwind turbine, when the blade is directly in front of the tower, it experiences minimum wind.For a three-bladed wind turbine, each of the blades experience minimum wind in one complete rotation and hence tower shadow contributes to the 3p effect [1,2].This section models tower shadow effect for an upwind turbine and identifies the factors affecting tower shadow.

Modeling the Disturbance in Wind due to Tower Shadow.
Figure 2 shows an aerial view of the tower immersed in uniform wind flow.As seen from the figure, when the uniform wind is incident on the tower, it tries to bend around the tower.At the point marked as "source," the wind streamlines radiate outwards.At the point marked as "sink,"  the wind streamlines converge.A source and sink at very close spacing (separated only by the tower diameter) form a doublet [13].
Hence, from potential flow theory, stream function of wind around the tower can be derived by superposing a doublet on uniform flow as [13] Differentiating stream function ψ with respect to y yields the total wind flow velocity V (y, x) in the x-direction as ( The second term in (5), that is, V tower (y, x) represents the disturbance in wind due to the presence of the tower [1,10].Here a is the tower radius, x is the distance of blade origin from the tower midline and y is the lateral distance of the blade from the tower midline, as shown in Figure 3. Writing V tower (y, x) as a function of r and θ yields

ISRN Renewable Energy
It should be noted that ( 6) is valid only when the blades are in the tower shadow zone, that is, 90 • ≤ θ ≤ 270 • , as shown in Figure 3.As justified in [1], the hub height wind speed in ( 5) is assumed to be equal to the spatial average wind speed.

Factors Affecting Tower Shadow.
From ( 6), the disturbance in wind due to tower-shadow effect depends on three factors, r, θ, and x.Hence, in the first case, r is varied from 5 m to 20 m in (6) while θ is varied from 90 • to 270 • , that is, the tower-shadow zone, as shown in Figure 4(a).It is observed that the blade elements closest to the hub (r = 5 m) experience a longer duration of the tower shadow.However at θ = 180 • , that is, blade directly in front of the tower, the same wind deficit is seen by all the blade elements.
Next, the radial distance of blade element from the hub, that is, r is kept constant at 20 m and the distance of blade origin from the tower midline, that is, x is varied from 2 m to 5 m.As seen from Figure 4(b), the tower shadow effect is most pronounced when the blades are closer to the tower (x = 2 m).Hence, mounting the blades away from the tower minimizes tower-shadow effect.Also note that there is an acceleration in wind when the blades are at either side of the tower, as seen in Figure 4.This is true since when the streamlines bend around the tower, they become concentrated on the sides of the tower, as shown in Figure 2, resulting in an acceleration of wind.

Time-Domain Modeling of Wind Turbines
From a modeling standpoint, the primary components of a wind turbine consist of turbine rotor (prime mover), shaft and gearbox unit (drive-train and speed changer), induction generator, and a control system [11].Wind turbines extract the kinetic energy available in winds and convert it to electrical energy and hence require interaction between aerodynamic, mechanical, and electrical components, as illustrated in Figure 5.The modeling of these four building blocks are discussed in this section.The modeling of the first three modules have been discussed in details in [11].This subsection discusses the modeling of the aerodynamic torque module and the two additional blocks on torque due to wind shear and tower shadow.Note that the input to the aerodynamic block is the hub height wind speed (V h ) which can be set to any value during simulation run-time or can also be in the form of time-series data of an actual wind speed measurement [11].

Torque Computation.
The wind turbine experiences a relatively constant hub height wind speed V h .However, due to tower shadow and wind shear, the blades of a wind turbine experience a varying wind speed as it rotates through a complete cycle, as seen in Sections 2 and 3. Hence the approach to formulating a pragmatic torque equation incorporating the 3p effect is to first derive the total spatially varying wind field, which includes the contribution from tower shadow and wind shear, and then compute the resultant aerodynamic torque [1].The approach is briefly summarized below.
The total wind experienced by the wind turbine blades due to tower shadow and wind shear is where W s (r, θ) and V tower (r, θ, x) are given by ( 3) and (6).
Next, the total aerodynamic torque produced by a threebladed wind turbine when immersed in a wind field of v(t, r, θ) is,

Mechanical block
Electrical block where T aero (t, θ) is the aerodynamic torque, r 0 the radius at which the blade profile begins, R the radius of the rotor disk, A is the cross-sectional area of rotor, ρ is the air density, ω rotor is the rotor angular velocity, and s is defined to be equal to 1 − (r 0 /R) 2 .This equation has been obtained by linearizing the classical torque equation about the hub height wind speed [1,10].To remove the radial dependence in v(t, r, θ), [10] defines an equivalent wind speed v eq (t, θ), which would give rise to the same aerodynamic torque T aero (t, θ), such that, ISRN Renewable Energy Equating ( 8) and ( 9) and simplifying, the expression for equivalent wind speed is obtained as Equation ( 10) has been split into three components: v eq 0 due to the hub height wind speed, v eq ws due to wind shear, and v eq ts due to the tower shadow, where Note that ( 12) and ( 13) have been simplified by assuming r 0 = 0.The resultant aerodynamic torque due to v eq (t, θ) can be then obtained by substituting v eq (t, θ) in ( 9) and simplifying Equation ( 14) is the complete torque equation representing aerodynamic torque due to kinetic energy in wind as well as torque due to tower shadow and wind shear and has been implemented for modeling.The aerodynamic torque module, shown in Figure 5 implements the classical torque equation, that is, first term in (14) as The wind shear module implements the disturbance in torque due to wind shear, that is, the second term in ( 14) as The tower shadow module implements the disturbance in torque due to tower shadow, that is, the third term in ( 14) as Modeling ( 14) into three separate modules gives the wind turbine model greater flexibility allowing tower shadow and wind shear to be studied separately, as will be demonstrated in Section 5. Also note that the tower shadow and wind shear modules input the wind turbine specifications (H, x, a, α) for calculating the torque oscillations due to 3p effect, as shown in Figure 5. Few important points when modeling the above three modules are (i) C p (λ 0 ) and λ 0 should be computed at the hub height wind speed V h ; (ii) computation of v eq ws and v eq ts in ( 16) and ( 17) require instantaneous values of the three blade angles.
The azimuthal angle of the first blade (θ b1 ) can be obtained by integrating rotor angular velocity.Integrator has to be reset once the first blade angle reaches 360 • .The other two blade angles can be obtained by adding 120 • and 240 • to θ b1 , respectively;.
(iii) v eq ts should be computed only when the individual blade is in the tower-shadow zone, as discussed in Section 3.For example, the first blade (reference blade) will be in the tower-shadow zone only when θ b1 is between 90 • and 270 • .A logical operator should be used so that when the first blade is outside the tower-shadow zone, contribution to v eq ts from the first blade should be zero.
The total aerodynamic torque can be then obtained by summing up T classical aero , T wind shear aero , and T tower shadow aero , scaled by the gear ratio and used as an input to the drive-train model, as shown in Figure 5.

Mechanical Block.
The mechanical block consisting of the wind turbine shaft, generator shaft, and gearbox has been modeled using a two-mass inertia representation, since the two-mass model is simpler to model and shows the dynamic response of the wind turbine generator.State-space equations governing the two-mass model have been detailed in [11].

Control
Block.This block can be developed for adjusting the pitch angle in variable speed wind turbines or controlling reactive compensation devices to meet wind turbine reactive power demand.For a fixed-speed wind turbine, no control is required [11].

Case Studies: Performance Evaluation of Wind Turbine Model
The aforementioned wind turbine components are integrated to form a 1.5 MW fixed speed wind turbine model.The numerical data set used for simulation of wind turbine model is listed in Table 2.The wind turbine is then connected to a weak grid via a 2-MVA 0.6/20 kV delta-wye-grounded transformer with a leakage inductance of 5% and a 6 mile medium voltage cable [14].The grid has been intentionally chosen to be weak so as to study the impact of tower shadow and wind shear on power quality.The grid short-circuit capacity (SCC) was selected to be 25 MVA, since the grid is considered electrically weak when [15] SCC where P n is the rated power of the wind turbine, that is, 1.5 MW in this case.X/R ratio of the grid is assumed to be 6.At the point of common coupling (PCC), a 1-MVA, 0.98 power factor load is connected.Figure 6 shows the complete system arrangement.Case studies presented in this section aim at evaluating the wind turbine model performance and demonstrate the analyses that can be done using this model.The simulation package used is PSCAD/EMTDC.

Performance of the Wind Turbine without 3p
Effect.The performance of the wind turbine without 3p effect (tower shadow and wind shear modules disconnected) is studied to determine the rotor angular velocity at the rated wind speed.The rated wind speed can be obtained from the power curve shown in Figure 7.The wind turbine model is subjected to hub height wind speed (V h ) ranging from 0 to 20 m/s.Since this is a fixed speed wind turbine, pitch angle remains constant at −6.15 • .The rated wind speed is 15 m/s and rated power is 1.5 MW.The rotor angular velocity ω rotor is 1.8 rad/s or 0.286 Hz. require the values of C p (λ 0 ) and λ 0 at the hub height wind speed V h .Sample/hold circuits are used in conjunction with switch SW such that when SW is switched "ON," the S/H circuits hold the values of C p (λ 0 ) and λ 0 , as illustrated in Figure 8. Depending on the simulation time step, the initialization step requires a few seconds, that is, 10 seconds in PSCAD/EMTDC.

Output Power from Wind Turbine.
The wind turbine is subjected to rated V h .After initialization, the wind shear and tower shadow modules are activated by toggling the switch (SW = 1).The output power from the wind turbine starts to oscillate as soon as the tower shadow and wind shear modules are activated, as shown in Figure 9.An FFT block with fundamental frequency equal to rotational frequency of the blades (ω rotor = 0.286 Hz) is used to determine the frequency components present in the output power.As shown in Figure 10, FFT indicates the presence of a significant third harmonic, that is, 0.86 Hz (which is three times ω rotor ), whose magnitude is 0.0516 MW.

Increase in Output
Power due to 3p Effect.The output power swings beyond the rated power (1.5 MW) when tower shadow and wind shear modules are switched "ON," as shown in Figure 9.This is counterintuitive since output power should decrease as the 3p effect leads to a net decrease in the wind speed.Figure 11 shows the equivalent wind speed and aerodynamic torque scaled down to bring them on the same scale as the output power.As seen from the figure, a decrease in equivalent wind due to 3p effect results in a corresponding decrease in torque.The output power, however, seems to peak at that time instant.The reason for the output power swinging beyond the rated power was determined to be the transient response of the generator.The oscillations in torque continually disturb the system, that is, the generator is not allowed to reach steady state, which is why the output power is swinging above the rated power.To prove the above theory, the inertia of the turbine J rot was increased from 1000 kgmm to 2000 kgmm, since inertia plays a significant role in controlling the transient response.Now, the 3p oscillations in the output power does not swing beyond the rated power, as shown in Figure 11.

Impact of Wind Shear and Tower Shadow Individually on
Output Power.This section demonstrates how the model can be used to study the impact of tower shadow and wind shear individually on the magnitude of output power oscillations.The wind turbine model is initialized by subjecting it to rated V h and then switch SW is used to activate the tower shadow module only.The output power oscillates between 1.528 MW and 1.429 MW, as shown in Figure 12.When switch SW activates the wind shear module only, the output power oscillates between 1.469 MW and 1.461 MW.The magnitude of the oscillations are very small, indicating that tower shadow effect has a more dominant impact on output power that wind shear.This is consistent with the analysis reported in [1].during the initialization of the model.Once the tower shadow and wind shear modules are switched "ON," the rms voltage shows the presence of the 3p oscillations.These fluctuations in system rms voltage may cause consumer annoyance and complaint and are known as voltage flicker.The usual method of expressing flicker is similar to that of voltage modulation, that is, change in voltage with respect to the average voltage over a certain time period as [6] % Here V max and V min are the maximum and minimum values of modulated signal while V 0 is the average value of normal operating voltage.Using (19), the voltage flicker (expressed in percentage of voltage modulation) was calculated to be 0.186%.This is below the threshold of objection.However, in a wind farm, if many wind turbine blades synchronize, the magnitude of voltage flicker would be severe [2].

Factors Affecting Voltage Flicker.
Voltage flicker due to tower shadow and wind shear depends on [4,5] (i) hub height wind speed V h ; (ii) grid short-circuit capacity; (iii) grid impedance angle or X/R ratio.
To investigate how these parameters influence voltage flicker, they are varied one at a time.The parameters used in Figure 6 will be used as the base case, that is, V h , SCC, and X/R ratio of grid are 15 m/s, 25 MVA, and 6, respectively.Turbulence intensity and loads connected at the PCC also affect voltage flicker, however, they are not studied in this section.

Impact of V h on Voltage
Flicker.The wind turbine was subjected to wind speeds ranging from cut-in to cutout speed and the voltage flicker at each wind speed was calculated using (19).As seen from Figure 14, voltage flicker increases almost linearly with wind speed till the rated wind speed is reached.For V h greater than the rated speed, voltage flicker starts decreasing due to the stall regulation of the rotor which significantly reduces the turbulence-induced fluctuations.Hence, voltage flicker is maximum at the rated wind speed or, in other words, voltage flicker is maximum when wind power output is maximum.

Impact of Grid Short-Circuit Capacity (SCC) on Voltage
Flicker.The SCC or strength of the grid is varied at each simulation.Figure 15 illustrates the approximately inversely proportionally relationship with short-circuit capacity (SCC) of grid.Higher the SCC, stronger the grid is, and hence voltage flicker is minimum.

Impact of Network Impedance Angle (X/R Ratio) on
Voltage Flicker.In this subsection, X/R ratio of the grid is varied from 1 to 7 while the grid SCC is kept constant at 25 MVA.As shown in Figure 16, minimum voltage flicker is obtained for an X/R ratio of 2.
The V-shaped curve can be explained by considering voltage change (ΔU) across the grid impedance (R and X) as [5] ΔU = PR cos θ − θ pf Here P is the real power produced by the wind turbine, θ is the network impedance angle, θ pf is the power factor at which the wind turbine is operating while U g represent the voltage at the PCC.From (20), the determining factor for minimizing the voltage drop is cos(θ − θ pf ).When (θ − θ pf ) approaches 90 • , voltage flicker is minimum, as shown in Figure 17.Since θ pf is approximately 26 • in the simulation case, minimum voltage flicker occurs for θ = 64 • or X/R = 2.

Conclusion
This paper presents a detailed model of a fixed speed wind turbine with an upwind rotor.The modeling includes a detailed representation of tower shadow and wind shear modules which are often approximated or neglected in existing models.The model will allow the user to input wind turbine specifications and observe how those parameters will directly affect the magnitude of 3p oscillations.Another advantage is that the model is flexible, allowing the user to study tower shadow and wind shear individually.Case studies have been carried out to use this model for studying the integration of wind turbine to the grid and analyzing voltage flicker at the PCC.The results are consistent with those reported in the literature and hence the model performs as expected.In addition, the model can also be used to study the mitigation of 3p effect and expanded to include power and pitch controllers to suit other studies involving variable speed wind turbines.

5 H
13.5 m/s V h = 14 m/s (a) Effect of wind speed (V h ) on wind shear = 80 m H = 90 m H = 100 m H = 110 m H = 120 m (c) Effect of hub height (H) on wind shear Effect of radial distance (r) on wind shear

Figure 3 :
Figure 3: Dimensions used in tower shadow formula.

Figure 4 :
Figure 4: Comparison of tower shadow at different (a) radial distances r and (b) lateral distances from tower midline x.

Figure 5 :
Figure 5: Block diagram of a wind turbine.

Figure 6 :
Figure 6: Fixed-speed wind turbine connected to the grid.

Figure 7 :
Figure 7: Power curve of a 1.5 MW fixed speed wind turbine.

Figure 12 :
Figure 12: Magnitude of 3p oscillations due to tower shadow and wind shear individually.

Table 1 :
Variation of α with terrain.

Table 2 :
Specifications of wind turbine model.
5.5.Voltage at the PCC.Rms voltage was measured at PCC as shown in Figure13.The voltage remains constant at 11.28 kV