Open access peer-reviewed chapter

# Effect of Turbulence on Fixed-Speed Wind Generators

By Hengameh Kojooyan Jafari

Submitted: March 7th 2012Reviewed: July 9th 2012Published: November 21st 2012

DOI: 10.5772/51407

## 1. Introduction

The influence of wind energy connection to the grid has increased greatly and turbulence or unreliable characteristics of wind energy are expected to produce frequency and voltage changes in power systems and protection system equipment. To prevent these changes, it is necessary to study the working point change due to turbulence. In other papers, the voltage and transient stability analysis have been studied during and after turbulence  and the impact of WTGs (wind turbine generators) on the system frequency, inertia response of different wind turbine technologies, and comparison between inertia response of single-fed and doubly-fed induction generators have been examined. Moreover study of the frequency change alone was conducted using Dig-SILENT simulator for FSWTs (fast-speed wind turbines) with one-mass shaft model .

In this chapter both frequency and grid voltage sag change are presented with MATLAB analytically and also by SIMULINK simulation in FSWTs with one- and two-mass shaft turbine models to compare both results and a new simulation of induction machine without limiter and switch blocks is presented as a new work. The first part of study is frequency change effect on wind station by SIMULINK that shows opposite direction of torque change in comparison with previous studies with Dig-SILENT. The second part of study is effect of frequency and voltage sag change on wind station torque due to turbulence in new simulation of induction generator that is new idea.

## 2. Wind turbine model

The equation of wind turbine power is

E1

The C p curve and equation are shown in Fig. 1 and given by equation (2) and (3)

E2
E3

whereθpitchis blade pitch angle, λis the tip speed ratio described by equation (4). The parameters are given in Table 1.

E4 Figure 1.Curve of C p for different tip speed ratios λ .

The curve of Fig.1 has positive slope before Cpmax and it has negative slope after Cpmax.

## 3. One-Mass Shaft Wind Station Model

Induction machine equation is

E5

Where, Tmis the mechanical torque, Teis the generator torque, Cis the system drag coefficient and Jis the total inertia.

Table 1 shows the parameters of the one-mass shaft turbine model and induction generator.

 Generator Wind Turbine Rs= .011Ω c1=.44 Ls= .000054H c2 = 125 Lm= .00287H c3= 0 L′r= .000089H c4= 0 R′r= .0042 [Ω] c5= 0.1 Jm=.5 to 20.26 [kgm2] c6= 6.94 p(#pole   pairs)= 2 c7= 16.5 Pn= 2e6 [w] c6= 0.1 c9= -.002 R= 35 [m] A= πR2 [m2] ρ=1.2041 [kg/m3] vw= 6, 10, 13 [m/s] θpitch= 0 [º]

## 4. Two-Mass Shaft Induction Machine Model

This model is used to investigate the effect of the drive train or two-mass shaft, i.e., the masses of the machine and the shaft, according to the equation (8) , . In this equation,Jtis wind wheel inertia, JGis gear box inertia and generator’s rotor inertia connected through the elastic turbine shaft with a κas an angular stiffness coefficient and Cas an angular damping coefficient.

The angular shaft speed ωtcan be obtained from equations (6) and (7) , , .

The Parameters, defined above, are given in Table 2.

This model is described as equation (8).

E6
E7
E8
 υ 1/80 JG[kg.m2] .5 Jt [kg.m2] 1 C[Nm/rad2] 1e6 κ[Nm/rad] 6e7

### Table 2.

Parameters of two-mass shaft model.

## 5. Induction Machine and Kloss Theory

In a single-fed induction machine, the torque angular speed curve of equation (12)  is nonlinear, but by using the Kloss equation (13), equations (9), (10), and (11), this curve is linearly modified ,  as shown in Fig. 2. Therefore, the effect of frequency changes in wind power stations can be derived precisely by equation (12) and approximately using equation (13), as shown in Figs. 26.

E9
E10
E11
E12
E13

Equations (11) and (12) are given in per unit, but the associated resistances are in ohms. Figure 3.Mechanical and linear electrical torque versus slip. Figure 4.Mechanical and electrical torque versus frequency curves per unit with V sag = 10% Figure 5.Mechanical and electrical torque versus frequency per unit with V sag = 20%. Figure 6.Mechanical and electrical torque versus frequency per unit with V sag = 50%.

Figs. 3, 4, 5, and 6 illustrate that for lower wind speeds of 6 and 10 m/s, as the synchronous frequency f s and V sag change, the T e and TmMathType@MTEF@5@5@+=feaagCart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadsfadaWgaaWcbaGaamyBaaqabaaaaa@37E1@values of the rotor change in the same direction as the frequency of the network, as shown in Tables III, IV, V, and VI. These figures and tables give the results for V sag = 0% (i.e., only the frequency changes), 10%, 20%, and 50%. However, for a higher wind speed of 13 m/s, as f s and V sag change, the T e and TmMathType@MTEF@5@5@+=feaagCart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadsfadaWgaaWcbaGaamyBaaqabaaaaa@37E1@values of the rotor change in the opposite direction to the changes in the frequency of the network.

For small changes in the slip according to the Kloss approach in equation (13), the torque changes as follows :

E14

Then:

E15

and

E16

or

E17

Thus, the new angular operation speed is

E18
 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .96050 -.1157 1.0005 -.1064 1.0405 -.0974 10 .9621 -.5337 1.0021 -.491 1.0421 -.4493 13 .9631 -.7863 1.0035 -.8122 1.0439 -.8331

### Table 3.

Analytical MATLAB results for different frequencies.

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .9606 -.1156 1.0006 -.1064 1.0406 -.0974 10 .9625 -.5163 1.0027 -.5137 1.0429 -.5086 13 .9738 -.7868 1.0043 -.8127 1.0448 -.8335

### Table 4.

Analytical MATLAB results for Vsag= 10%.

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .9607 -.1156 1.0007 -.1064 1.0407 -.0974 10 .9632 -.5163 1.0034 -.5136 1.0437 -.5085 13 .9648 -.7875 1.0054 -.8134 1.0461 -.8341

### Table 5.

Analytical MATLAB results for Vsag= 20%

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .9618 -.1153 1.0018 -.1061 1.0418 -.0971 10 .9681 -.5161 1.0088 -.5131 1.0494 -.5076 13 .9724 -.7927 1.0139 -.8181 1.0555 -.8382

### Table 6.

Analytical MATLAB results for Vsag= 50%

## 6. Simulation of wind generator with frequency change

During turbulence and changes in the grid frequency, the torque speed (slip) curves change in such a way that as the frequency increases, the torque is increased at low wind speeds; 6 and 10 m/s, in contrast to Fig. 6 and decreases at a high speed of 13 m/s , as shown in Table 7 and Figs. 715.

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .9619 -.1148 1.0019 -.1057 1.0418 -.0969 10 .9684 -.5179 1.0091 -.5134 1.0494 -.5076 13 .9724 -.7945 1.0147 -.8177 1.0559 -.8373

### Table 7.

Simulink simulation results for one- and two-mass shaft models

Figs. 715 show the electrical torque and mechanical speed of the induction machine for the one- and two-mass shaft turbine models at wind speeds of 6, 10, and 13 m/s to validate Table 7.

## 7. Simulation of wind station with one-mass and two-mass shaft turbine models

The results of simulations of a simple grid, fixed-speed induction machine, and one-mass and two-mass shaft turbines are given in Tables 8 -10 and Figs. 1642. For an induction wind generator using the induction block in SIMULINK with high voltage sag i.e. 50% with frequencies 50 and 52 and  equal to 13, C p becomes negative, and the results are unrealistic. Then results of 50% voltage sag are realistic in new simulation of induction machine in Tables 8 -10.

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .9624 -.1152 1.0024 -.106 1.0423 -.097 10 .9703 -.516 1.0111 -.5128 1.0519 -.5071 13 .9757 -.795 1.0176 -.8201 1.0595 -.8399

### Table 8.

Simulation results by SIMULINK for one and two mass shaft model for Vsag= 10%

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .963 -.1151 1.003 -.1059 1.043 -.0969 10 .973 -.5159 1.014 -.5125 1.055 -.5066 13 .9799 -.7977 1.0223 -.8226 1.0648 -.842

### Table 9.

Simulation results by SIMULINK for one and two mass shaft model for Vsag= 20%

 υw fs= 48 fs= 50 fs= 52 ωm[pu] Te[pu] ωm[pu] Te[pu] ωm[pu] Te[pu] 6 .9674 -.114 1.0074 -.1048 1.0474 -.0959 10 .9933 -.5146 1.0364 -.5096 1.0796 -.502 13 1.0248 -.8239 1.0474 -.8347 1.0917 -.85

### Table 10.

Simulation results by SIMULINK for one and two mass shaft model for Vsag= 50%

## 8. New Simulation of Induction Machine

Figs. 33 and 42 show the results of new simulation of the induction machine model illustrated in Fig. 43 . The new simulation, which has no limiters and switches, is used because at high grid voltage drop-down or sag, the Simulink induction model does not yield realistic results.

The new simulation of induction machine is in dqosystem and synchronous reference frame simulation on the stator side; n (Transfer coefficient) is assumed to be 1. Circuit theory is used in this simulation, and it does not have saturation and switch blocks, unlike the MATLAB–SIMULINK Induction block. In Fig. 43, LMis the magnetic mutual inductance, and randLare the ohm resistance and self-inductance of the dqocircuits, respectively. The machine torque is given by equation (19). In this equation, id,qsandλd,qs, the current and flux parameters, respectively, are derived from linear equations (20)–(23); they are sinusoidal because the voltage sources are sinusoidal.

E19

Where Pis poles number, λdsand λqsare flux linkages and leakages, respectively, and iqsandidsare stator currents in qand dcircuits of dqosystem, respectively.

Then imatrix produced by the λmatrix is given by equation (20).

E20

where the inductance matrix parameters are given by (21), (22), (23).

E21
E22
E23

The linkage and leakage fluxes are given by (24) to (29).

E24
E25
E26
E27
E28
E29

To create the torque in equation (19), it is necessary to determine the currents in equations (30)–(33) from the stator and rotor currents by using current meters.

E30
E31
E32
E33

## 9. Conclusion

As frequency changes and voltage sag occurs because of turbulence in wind stations in ride-through faults, the system’s set point changes. The theoretical and simulation results results are similar for one mass shaft and two mass shaft turbine models. At lower wind speeds; 6 and 10 m/s, the directions of the changes in the new working point are the same as those of the frequency changes. At a higher wind speed; 13 m/s, the directions of these changes are opposite to the direction of the frequency changes. Simulation results of high grid voltage sag with SIMULINK induction block has error and new simulation of wind induction generator in synchronous reference frame is presented without error and in 50% voltage sag, new simulation of wind generator model has higher precision than that in 10% and 20% voltage sags; however, this model can simulate wind generator turbulence with voltage sags higher than 50%. Although results of new simulation of induction machine with wind turbine for 50% voltage sag and frequencies 50 and 52 have been presented in this chapter.

## 10. Nomenclature

P=Generator power

ρ=Air density

A=Turbine rotor area

Cp=Power Coefficient

υw=Wind speed

θpitch=Pitch angle

Te=Electrical torque

Tm=Mechanical torque

J = Inertia

ωm=Mechanical speed

C=Drag coefficient

ν=Gear box ration

Rs = Stator resistance

Ls = Stator inductance

 Lm= Mutual inductance

 Lr= Rotor inductance

Rr= Rotor resistance

p=Pole pairs

κ=Stiffness

λr,s=Rotor and stator flux

Kr,s=Rotor and stator park transformation in synchronous reference frame

ir,s=Rotor and stator current

vr,s=Rotor and stator voltage

## 11. Future Work

The new simulation of induction generator will be tested by new innovative rain turbine theory and model of the author.

## Acknowledgments

I appreciate Dr. Oriol Gomis Bellmunt for conceptualization, Discussions and new information and Dr. Andreas Sumper for discussions about first part of chapter, with special thanks to Dr. Joaquin Pedra for checking reference frame and starting point in new simulation of induction machine.

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2012 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Hengameh Kojooyan Jafari (November 21st 2012). Effect of Turbulence on Fixed-Speed Wind Generators, Advances in Wind Power, Rupp Carriveau, IntechOpen, DOI: 10.5772/51407. Available from:

### Related Content

#### This Book

Edited by Rupp Carriveau

Next chapter

#### Turbine Wake Dynamics

By Phillip McKay, Rupp Carriveau, David S-K Ting and Timothy Newson

#### Fundamental and Advanced Topics in Wind Power

Edited by Rupp Carriveau

First chapter

#### Aerodynamics of Wind Turbines

By Emrah Kulunk

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.