Simulation of directly gridconnected wind turbines for voltage fluctuation evaluation.Abstract Power and voltage generated by a wind turbine are more variable than that produced by conventional generators. Therefore, with the increasing penetration of wind energy into the grid, it is very important to simulate precisely the wind turbine/grid interaction to evaluate voltage fluctuation Fluctuation A price or interest rate change. and other power quality aspects. In this paper, the wind turbine instantaneously in·stan·ta·ne·ous adj. 1. Occurring or completed without perceptible delay: Relief was instantaneous. 2. generated power and voltage at the point of common connection (PCC PCC prothrombin complex concentrate. ) with grid are simulated by considering all the aerodynamical and mechanical effects, which could affect them. The inherent effect of the wind speed on the entire blade swept area is simulated in the model of the wind speed. This model takes into account the effects of tower shadow and the rotational sampling of turbulence turbulence, state of violent or agitated behavior in a fluid. Turbulent behavior is characteristic of systems of large numbers of particles, and its unpredictability and randomness has long thwarted attempts to fully understand it, even with such powerful tools as . The generated power is obtained by the simulation of the wind speed time series into a wind turbine model. The proposed wind turbine model includes detailed submodels of aerodynamic rotor rotor: see generator; motor, electric. , drive train and ac generator. A dq arbitrary reference frame is used for modeling of asynchronous Refers to events that are not synchronized, or coordinated, in time. The following are considered asynchronous operations. The interval between transmitting A and B is not the same as between B and C. The ability to initiate a transmission at either end. generator. The flux flux In metallurgy, any substance introduced in the smelting of ores to promote fluidity and to remove objectionable impurities in the form of slag. Limestone is commonly used for this purpose in smelting iron ores. variations in the stator stator: see generator; motor, electric. circuit are involved. The flickermeter model which expresses voltage fluctuations is simulated according to according to prep. 1. As stated or indicated by; on the authority of: according to historians. 2. In keeping with: according to instructions. 3. the IEC (International Electrotechnical Commission, Geneva, Switzerland, www.iec.ch) An organization that sets international electrical and electronics standards founded in 1906. It is made up of national committees from over 60 countries. IEC  International Electrotechnical Commission standard 61000415. The wind speed, wind turbine and flickermeter models are simulated in Matlab/Simulink software. Both of grid and site parameters, which affect voltage fluctuation, are investigated. These parameters have a wide influence on voltage fluctuation and flicker flicker: see woodpecker. flicker Any of six species of New World woodpeckers (genus Colaptes) noted for spending much time on the ground eating ants. emission levels. Introduction Wind turbine generators are increasingly becoming among the prominent components of power systems. Due to the stochastic By guesswork; by chance; using or containing random values. stochastic  probabilistic nature of wind, electrical power delivered by this type of generation possesses similar features. Furthermore, wind turbines are usually integrated with the grid at remote terminals, far from central loads or conventional generation. This is raising certain reluctance on the part of utility companies to inject in·ject v. 1. To introduce a substance, such as a drug or vaccine, into a body part. 2. To treat by means of injection. that source of power with unknown behavior and to evaluate the accompanied power quality considerations. Due to the importance of windturbine power quality considerations, standard IEC 6140021 [1] provides procedures for determining the power quality characteristics of wind turbines. The sources of fluctuations in the generated power are due to stochastic aspects that determine wind speed at different times and heights, and to deterministic 1. (probability) deterministic  Describes a system whose time evolution can be predicted exactly. Contrast probabilistic. 2. (algorithm) deterministic  Describes an algorithm in which the correct next step depends only on the current state. or periodic effects. The largest periodic effect known as the tower shadow is at the frequency at which rotor blades pass by the tower. In the common three bladed horizontal axis wind turbines, this frequency is known as a 3p frequency which is three times the rotational frequency [2]. The reactive power reactive power: see power, electric. consumption of the wind turbine asynchronous generator depends on the active generated power. The drawn reactive power increases with the increase of the generated active power. Therefore, the reactive power consumption of the generator fluctuates as the wind speed fluctuates. Due to the fluctuations in the active and reactive power, the voltage at PCC fluctuates. Voltage fluctuation is a serious issue particularly for directconnected wind turbines because these turbines produce power dependent on the variations of the wind speed and inject it without conditioning into the grid. Voltage fluctuation disturbs the sensitive electric and electronic equipment. This may lead to a great reduction in the life span of most equipment [3]. The lighting flicker level is generally used to measure voltage fluctuation. A case studying the voltage level profile when the power system is integrated with wind generation is given [4]. The influence of wind turbines on consumer voltage quality is studied [5]. A frequency domain approach to wind turbines for flicker analysis is presented [6]. The need to accurately simulate the wind turbines and investigate their interaction with the grid is becoming so important, since the penetration of wind in certain areas reaches significant levels. A suggested nonlinear A system in which the output is not a uniform relationship to the input. nonlinear  (Scientific computation) A property of a system whose output is not proportional to its input. simulation depending on the collected measured data is presented [7]. This simulation employs the neural network neural network or neural computing, computer architecture modeled upon the human brain's interconnected system of neurons. Neural networks imitate the brain's ability to sort out patterns and learn from trial and error, discerning and extracting technique to predict the wind turbine output power. The modeling of wind turbines for power system studies is presented [8]. The aerodynamic loads of the wind turbines are represented and simulated in frequency domain [9]. A comprehensive model of mechanical part consists of a number of lumped inertias, elastically coupled to each other is presented [10], [11]. The problem in this model that it needs more manufacturer's design data about turbine elements which is not available in most cases. Soft tools and techniques used for modeling and design simulation of wind turbines are reviewed [12]. A simple approach to aggregated wind farm equivalent for the analysis of power system operation is presented [13]. This paper presents the comprehensive timedomain modeling of wind speed, wind turbine and flickermeter and investigates the influence different factors on voltage fluctuation caused by wind turbines. The wind speed model considers the effects of tower shadow and the interaction of blade rotation with wind turbulence. The wind turbine model is composed of three sub models connected together, aerodynamic rotor, drive train and induction generator In`duc´tion gen´er`a`tor 1. A machine built as an induction motor and driven above synchronous speed, thus acting as an alternatingcurrent generator;  called also asynchronous generator ltname>. models. A soft shaft along with the mass inertia inertia (ĭnûr`shə), in physics, the resistance of a body to any alteration in its state of motion, i.e., the resistance of a body at rest to being set in motion or of a body in motion to any change of speed or change in direction of of both the generator and wind turbine rotors is simulated in a 2mass drive train model. The induction generator model, in the dq0 arbitrary reference frame, involves the transients in the stator circuit. The flickermeter is implemented according to the IEC standard 61000415 [14]. The wind speed, wind turbine and flickermeter models are implemented in Matlab/Simulink. The wind speed produced from wind speed model is applied to the aerodynamic model to extract the aerodynamic torque. This torque is fed to the drive train/generator models to simulate the electrical power from wind turbine. Modeling The wind speed, wind turbine and flickermeter models are explained in the following subsections. Wind Speed Model The main idea in the wind speed model is to represent the interaction of turbine blades with wind speed distribution over the rotor swept area by an equivalent time series data at hub level. The equivalent time series data simulates the wind and when it is applied to a specific aerodynamic model, it reproduces the aerodynamic torque from a wind turbine. Equivalent wind speed ([W.sub.eq]), produced from wind speed model, takes into account the deterministic and stochastic wind speed actions on the rotor area. The [W.sub.eq] is given as follows [15]; [W.sub.eq] (t, [psi PSI  Portable Scheme Interpreter ])  [U.sub.eqdet] ([psi]) + [U.sub.eqsto] (t, [psi]) (2) where: [U.sub.eqdet] is the deterministic part of the wind speed; [U.sub.eqsto] is the stochastic part of the wind speed; [psi] is swept angle of the rotor blades; and t is the time. Deterministic Part The deterministic part of the wind is time independent and it depends only on the rotor position, [psi]. It consists of the mean wind speed, tower shadow and wind shear wind shear, a sudden, drastic change in wind direction or speed over a comparatively short distance. Most winds travel horizontally, as does most wind shear, but under certain conditions, including thunderstorms and strong frontal systems, wind shear will travel in a . The wind shear is very important to the analyses of blade loads but it is not transmitted to the electrical power so it is neglected. The deterministic part components are explained in the following items. a) Mean wind speed The mean wind speed is estimated from the site wind data or wind Atlas A wind atlas contains data on the wind speed and wind direction in a region. These data include maps, but also time series or frequency distributions. A climatological wind atlas covers hourly averages at a standard height (10 meter) over even longer periods (30 years). . It is converted from its measuring level to the wind turbine hub level according to the following formula [16]; [U.sub.o] = ([h.sub.o]/[h.sub.1])[g.sup.*] x [W.sub.1] (2) where: [W.sub.1] is the average wind speed measured at height [h.sub.1]; [U.sub.o] is the mean wind speed at hub level, [h.sub.o]; [g.sup.*] is an exponent exponent, in mathematics, a number, letter, or algebraic expression written above and to the right of another number, letter, or expression called the base. In the expressions x^{2} and x^{n}, the number 2 and the letter n varying from (0.16 to 0.40). [g.sup.*] = 0.16 for open land with a few obstacles and [g.sup.*] = 0.4 for land with big obstacles). b) Tower shadow Towers, in horizontal axis wind turbines, are obstacles to the free wind that modifies the wind flow. The upstream flow, which is interested in this paper, is reduced in front the tower and increased laterally. The model is limited to the effects on horizontal, upwind, and threebladed wind turbines. When the rotor blade crosses the tower, a drop in the aerodynamic torque will be occurred. Consequently, the Fourier transformation of the aerodynamic torque influenced by tower shadow shows the 3np (i.e., 3p, 6p, 9p, ....) components included in it. The wind turbine dynamics remove the higher 3np frequency components from the electrical output power and the dominant frequency of tower shadow is 3p [15]. Hence, the simulation of deterministic part of wind speed model is illustrated as block diagrams A chart that contains squares and rectangles connected with arrows to depict hardware and software interconnections. For program flow charts, information system flow charts, circuit diagrams and communications networks, more elaborate graphical representations are usually used. in Fig. 1 [15]. Where "[U.sub.det3]" is the contribution amplitude amplitude (ăm`plĭtd'), in physics, maximum displacement from a zero value or rest position. of the tower shadow, "[U.sub.o]" is the mean wind speed and "[U.sub.eq_det]" is the output deterministic wind speed. [FIGURE 1 OMITTED] Stochastic part The stochastic part depends on both the time and the rotor position, [psi]. Low frequency wind variations have large amplitudes and vice versa VICE VERSA. On the contrary; on opposite sides. . To model this characteristic, the Kaimal spectral spectral /spec·tral/ (spek´tral) pertaining to a spectrum; performed by means of a spectrum. spec·tral adj. Of, relating to, or produced by a spectrum. description of turbulence signal, K(f), is commonly used [15], [17]. Equation (3) defines the Kaimal turbulence filter turbulence filter the effect of the structure of the turbinates, trachea and bronchi in directing air flow against mucous membranes, serving to remove particulate matter from inspired air. model. f.K(f) /[[sigma].sup.2] = f.[X.sub.L]/[U.sub.o]/[(1+[1.5.sup.*] (f.[X.sub.L]/[U.sub.o])).sup.5/3] (3) Where: f is the frequency of the turbulence; [x.sub.L] is the turbulence length scale; and The interaction between the rotation of wind turbine blades and the variable wind speed due to wind turbulences is called the rotational sampling of turbulence. The rotational sampling of turbulence has been reported as the main cause of flicker from wind turbines. Therefore, the contribution of rotational sampling of turbulence on wind speed model is very necessary. The stochastic part, which includes the rotational sampling turbulence, is simulated by same way of deterministic component and depends mainly on (0p and 3p harmonic harmonic. 1 Physical term describing the vibration in segments of a soundproducing body (see sound). A string vibrates simultaneously in its whole length and in segments of halves, thirds, fourths, etc. ). However, the amplitudes are time dependent as shown in Fig. 2. The turbulence is simulated by random numbers filtered by Kaimal filter to eliminate high frequency components. Two secondorder filters are employed as admittance Admittance The ratio of the current to the voltage in an alternatingcurrent circuit. In terms of complex current I and voltage V, the admittance of a circuit is given by Eq. (1), and is related to the impedance of the circuit Z by Eq. (2). functions for the smoothing of the amplitudes of 0p and 3p harmonics har·mon·ic adj. 1. a. Of or relating to harmony. b. Pleasing to the ear: harmonic orchestral effects. c. after applying Kaimal filter [17]. The [W.sub.eq] is the sum of the stochastic and deterministic parts as presented in Eqn. 1.The rotor position, [psi], is obtained from the integral of the angular angular /an·gu·lar/ (ang´gular) sharply bent; having corners or angles. rotor speed. [FIGURE 2 OMITTED] Wind Turbine Model The wind turbine structure is generally composed of the rotor blades, a mean of drive train, and the electrical generator This article is about machines that produce electricity. For other uses, see Generator. “Dynamo” redirects here. For other uses, see Dynamo (disambiguation). . In this paper, the rotor blades have been connected to an induction generator through a flexible coupling followed by gear box. The induction generator is directly connected to the grid. The dynamic wind turbine parameters of a fixedspeed, stallregulated wind turbine are investigated. Figure 3 shows a layout of the proposed wind turbine model. Each submodel presented in Fig. 3 will be given in the following subsections. [FIGURE 3 OMITTED] Aerodynamic Rotor Model The aerodynamic rotor converts the wind into mechanical torque. The aerodynamic model uses the [W.sub.eq] and the rotor speed as input data to compute To perform mathematical operations or general computer processing. For an explanation of "The 3 C's," or how the computer processes data, see computer. the aerodynamic torque on the wind turbine shaft. The aerodynamic torque, [T.sub.aer], can be expressed in (4) [11]. [T.sub.aer] = [rho] x [pi]* [R.sup.2] x [W.sup.3.sub.eg] x [C.sub.p]([lambda])/ 2 [(omega).sub.wtr] (4) where: [[rho].sub.x] is the air density; R is the rotor radius; [C.sub.p]([lambda]) is the power coefficient of stall regulated wind turbine; [[omega].sub.wtr] is the rotational rotor speed. [lambda] is the tip speed ratio The tip speed ratio λ (lambda) or TSR for wind turbines is the difference between the rotational speed of the tip of a blade and the actual velocity of the wind. If the velocity of the tip is exactly the same as the wind speed the tip speed ratio is 1. (R x [[omega].sub.wtr] / [W.sub.eq]). Drive Train Model The drive train model simulates the elastic elastic Of or relating to the demand for a good or service when the quantity purchased varies significantly in response to price changes in the good or service. connection between the lowspeed wind turbine shaft and the highspeed generator shaft. It consists of the inertia of both the turbine and generator (2mass system). The elastic coupling Dr. Ing. Leonhard Geislinger invented the Elastic coupling in 1958. Originally it was a side product of the development for thermopneumatic locomotives, but unlike the trains, it was a big success. between the two masses is modeled by spring, k, and damper damp·er n. 1. One that deadens, restrains, or depresses: Rain put a damper on our picnic plans. 2. An adjustable plate, as in the flue of a furnace or stove, for controlling the draft. , B. The equivalent model for the drive train of the wind turbine referred to the generator side is shown in Fig. 4. [FIGURE 4 OMITTED] The dynamics of the drive train on the generator side can be written in (5) and (6) [8]: [T.sub.aer] = [J.sub.T] d[[omega].sub.wtr]/dt + k x [DELTA]v + B x [DELTA][omega]) (5) k x [DELTA]v + B x [DELTA][omega]) = [J.sub.G] d[[omega].sub.g]/dt + [T.sub.e] (6) Where: [DELTA][omega] = [[omega].sub.wtr]  [[omega].sub.g] [DELTA]v = [v.sub.wtr]  [v.sub.g] [[omega].sub.wtr] = d[v.sub.wtr]/dt and [[omega].sub.g] = d[v.sub.g]/dt [J.sub.T] (kg x [m.sup.2]) is the turbine moment of inertia, [J.sub.G] (kg x [m.sup.2]) is the generator moment of inertia, k (N x m/rad) is the stiffness of the shaft, B (N x m x sec/rad) is the absorption of the shaft, [T.sub.aer] (N x m) is the wind turbine torque, [T.sub.e] (N.m) is the electromagnetic generator torque, [[omega].sub.wtr], [[omega].sub.g] (rad [s.sup.1]) are the angular speed of the turbine and of the generator, respectively, and [v.sub.wtr] [v.sub.g] (rad) are the positions of the turbine and generator shafts, respectively. Induction Generator Model The asynchronous machine has been used as generator for many years in wind turbines. This type of machine has high reliability, low price, and low maintenance. The speed flexibility (slip), when compared to the synchronous Refers to events that are synchronized, or coordinated, in time. For example, the interval between transmitting A and B is the same as between B and C, and completing the current operation before the next one is started are considered synchronous operations. Contrast with asynchronous. machines reduces the current spikes due to wind gusts. The main disadvantage is that asynchronous machine consumes reactive power, which influences the wind turbine/grid integration in normal and transient conditions. A dq arbitrary reference frame is extensively used for modeling of several types of induction generators [18], [19], [20]. In this model the dynamic equations of induction machine in the dq0 arbitrary reference frame are employed. The variations in the stator circuit are involved. Flickermeter Model Flicker is defined as "impression of unsteadiness of visual sensation induced by a light stimulus whose luminance The amount of brightness, measured in lumens, that is given off by a pixel or area on a screen. For example, dark red and bright red would have the same chrominance, but a different luminance. or spectral distribution fluctuates with time" [21]. The flicker level depends on the amplitude, shape and repetition frequency of voltage fluctuation. Flicker level is measured by the use of a flickermeter described in IEC 61000415 [14]. The flickermeter takes voltage as input and gives the shortterm flicker index, [P.sub.st], as output. The flickermeter sequence consists of five steps; Step1: input voltage adaptor An alternate spelling of "adapter." See adapter. (tool) Adaptor  (Automatic DAta Parallelism TranslatOR) A source to source transformation tool that transforms data parallel programs written in Fortran 77 with array extensions, parallel loops, and layout directives to Step2: square law demodulator See demodulate. Demodulator A device used to recover the original modulating signal from a modulated wave. A demodulator is also known as a detector. Step3: weighting filters Step4: squaring and smoothing Step5: statistical manipulation The steps 14 are simulated in Matlab/simulink. The output of step 4 is the instantaneous in·stan·ta·ne·ous adj. 1. Occurring or completed without perceptible delay: Relief was instantaneous. 2. flicker level (IFL IFL International Fight League (mixed martial arts) IFL Integrated Facility for Linux (IBM) IFL Intense Football League IFL Indoor Football League IFL I'm Feeling Lucky (Google) ). Using the statistical manipulation over IFL for an observation period of 10 min, the [P.sub.st], can be estimated. The normalized response of the flickermeter (Fig. 5) shows that a quite small voltage fluctuation at certain frequency (8.8 Hz) can be irritable irritable /ir·ri·ta·ble/ (ir´itahb'l) 1. capable of reacting to a stimulus. 2. abnormally sensitive to stimuli. 3. prone to excessive anger, annoyance, or impatience. . The integration of wind turbines with the distribution networks may affect different issues of power quality. The flicker level ([P.sub.st][less than or equal to]1) can be a limiting condition for connecting wind turbines to low voltage Low voltage is an electrical engineering term that broadly identifies safety considerations of an electricity supply system based on the voltage used. While different definitions exist for the exact voltage range covered by "low voltage", the most commonly used ones include "mains networks [21]. [FIGURE 5 OMITTED] Voltage Fluctuation Caused by Wind Turbine The wind turbine induction generator injects fluctuated active power into the grid and correspondingly it draws fluctuated reactive power from the grid. Figure 6 shows a simplified diagram of wind turbine generator connected to grid. If the variation of the active power injected in·ject·ed adj. 1. Of or relating to a substance introduced into the body. 2. Of or relating to a blood vessel that is visibly distended with blood. injected 1. introduced by injection. 2. congested. to the grid is [DELTA]P and the corresponding variation of the reactive power absorbed from the grid is [DELTA]Q, then voltage fluctuation, at PCC, [DELTA]V/V v/v volume (of solute) per volume (of solvent). v/v volume (of solute) per volume (of solvent). , is given in (7). The nominal voltage is 1 pu. [DELTA]V/V [approximately equal to] [DELTA]P x R  [DELTA]Q x X (7) which can be rewritten as follows; [DELTA]V/V [approximately equal to][DELTA]S x Z x cos([theta Theta A measure of the rate of decline in the value of an option due to the passage of time. Theta can also be referred to as the time decay on the value of an option. If everything is held constant, then the option will lose value as time moves closer to the maturity of the option. ] + ([phi]) (8) where: R is the resistance of the grid impedance impedance, in electricity, measure in ohms of the degree to which an electric circuit resists the flow of electric current when a voltage is impressed across its terminals. (pu); X is the reactance of the grid impedance (pu); [DELTA]S is the apparent power variation, [DELTA]S = [square root of [P.sup.2] + [DELTA][Q.sup.2] (pu); Z is the grid impedance amplitude (pu); [theta] is the grid impedance angle; and [phi] is [tan.sup.1]([DELTA]Q/[DELTA]P). [FIGURE 6 OMITTED] Simulation Results Samples of Results of Wind Speed and Wind Turbine Modeling Using the parameters of the wind turbine generator, which are given in the appendix, the wind speed and wind turbine models have been implemented under Matlab/Simulink. The wind turbine ratings are taken as base quantities. The grid fault level and X/R ratio are assumed 50 pu and 10 respectively. The noload reactive power demand of the induction generator under study is compensated by capacitor capacitor or condenser, device for the storage of electric charge. Simple capacitors consist of two plates made of an electrically conducting material (e.g., a metal) and separated by a nonconducting material or dielectric (e.g. banks installed at PCC. Figure 7 presents a sample of selected model outputs; equivalent wind speed, active/reactive power and voltage at PCC. [FIGURE 7 OMITTED] The wind speed variation, assuming turbulence intensity of 10%, and [U.sub.o] = 12 m/sec., is shown in Fig. 7(a). The corresponding aerodynamic power and electrical output power are shown in Fig. 7(b). Not all variations in the aerodynamic power are transmitted to the electrical power. It means that the soft shaft coupling and the inertia of wind turbine rotor and generator rotor masses damp and smooth the variations of the aerodynamic power. Figure 7(c) shows the variation of the absorbed reactive power from grid. It is clear that the voltage fluctuation, in Fig. 7(d), is very small due the high fault level of the grid. Figure 8 traces the variation of the absorbed reactive power with the injected active power for the integrated wind energy system and the induction generator under study. The drawn reactive power increases with the increase of the generated active power. This relation is known as QP characteristic of the induction generator. [FIGURE 8 OMITTED] Influences of Grid Parameters on Flicker Caused by Wind Turbines The grid parameters affecting the wind turbine flicker emission are the fault level, and X/R ratio of the gird impedance. Fault Level Figure 9 shows the variation of the shortterm flicker index with different grid fault levels. This simulation test is carried out with two cases of grid impedance angles; 45[degrees] and 70[degrees]. The mean wind speed at hub level is maintained at 12 m/sec. and site turbulence is 10%. From Fig. 9, it can be seen that the flicker level decreases with the increase of fault levels. [FIGURE 9 OMITTED] XIR XIR Externally Initiated Reset XIR Xtreme Instant Rumble Clan (gaming group) Ratio of grid impedance The X/R ratio of the gird impedance is studied in terms of the impedance angle, [tan]=[tan.sup.1](X/R). Figure 10 shows the variation of the shortterm flicker index with the impedance angle. The fault level is maintained at 10 pu. The mean wind speed at hub level is 12 m/sec. and site turbulence is 10%. From Fig. 6 it can be seen that the flicker decreases with the increase of [theta] untill the minimum point, then the slope is reversed and flicker increases with increase of [theta]. The 'V' curve trend for the relation of [P.sub.st] with the grid impedance angle agrees with those given in [22], [23]. However, in [24], there is no indication to the minimum point and states that flicker increases with low X/R ratio. The minimum point of voltage fluctuation is occurred when [theta] + [phi] = 90[degrees]. The angle [phi] can be obtained according to the incremental Additional or increased growth, bulk, quantity, number, or value; enlarged. Incremental cost is additional or increased cost of an item or service apart from its actual cost. variation, [DELTA]Q/[DELTA]P, of QP characteristic of the induction generator (Fig. 8). The QP characteristic of the generator determines the [phi]P relation. Then the mean generated power (operating point) determines the value of [phi] and consequently [theta]=90[degrees]  [psi]. at the point of minimum flicker emission. Figure 11 shows the variation of angle [phi] with the active generated power. The mean active generated power is estimated accordance to the wind turbine power curve at mean wind speed. In this case, it is approximately 0.6 pu. as shown in Fig. 7(b). Then [phi] [approximately equal to] 8.6[degrees] which gives [theta] [approximately equal to] 81.4[degrees] at minimum flicker emission as shown in Fig. 10 The minimum point of flicker emission cannot be zero because the power swings up and down the mean value. [FIGURE 10 OMITTED] [FIGURE 11 OMITTED] Effects of Site Parameters on Wind Turbine Flicker Emission The effects of mean wind speed and the site turbulence on flicker emission will be discussed. Mean Wind Speed The flicker severity is calculated for the same site turbulence ([I.sub.u] = 10%) and different mean wind speed, [U.sub.o]. The fault level of the grid is 10 pu. The [P.sub.st] values versus the mean wind speed are estimated for two grid impedance angles, 45[degrees]and 70[degrees]. [FIGURE 12 OMITTED] [FIGURE 13 OMITTED] The [P.sub.st] variation with the mean wind speed is illustrated in Fig. 12. It can be explained by the wind turbine power curve, shown in Fig. 13. From Fig. 13, it can concluded that the output power and its fluctuations in the low wind region are low and therefore the induced voltage fluctuation is small. As the wind speed increases from cut in speed to 13 m/sec., the output power fluctuations and [P.sub.st] increase, approximately in linear relation with the mean wind speed. However, in the stall region (greater than 13 m/sec.), the rate of change of the aerodynamic power curve is reduced, resulting in a corresponding reduction in the output electrical power variability. Then [P.sub.st] values vary in small range with the increase of the wind speed as shown in Fig. 12. The trend of [P.sub.st] versus the mean wind speed, shown in Fig. 12, agrees with this given in literature [22]. However, in [23], [25], the flicker emission increases with increasing wind speed. Turbulence Intensity In this case, the mean wind speed at hub level is maintained at 12 m/sec., the fault level is 10 pu., and the grid impedance angle is 45[degrees]. Figure 14 shows the variation of [P.sub.st] with the turbulence intensity. The increase in the wind speed turbulence increases the power variability then the flicker emission increases with the increase of the turbulence intensity as shown in Fig. 14. [FIGURE 14 OMITTED] Conclusion The comprehensive wind speed and wind turbine models can be applied for voltage fluctuation and power quality studies. The wind speed model involves the tower shadow and the rotational sampling of turbulence. A 2mass drive train model, including soft shaft coupling, is used for damping damping In physics, the restraint of vibratory motion, such as mechanical oscillations, noise, and alternating electric currents, by dissipating energy. Unless a child keeps pumping a swing, the backandforth motion decreases; damping by the air's friction opposes the and smoothing the output electrical power variations. A dq arbitrary reference frame is used for modeling of the asynchronous generator. The variations in the stator circuit are involved. The flicker level caused by voltage fluctuation is evaluated by the flickermeter, described in IEC 61000415. From simulation results, voltage fluctuations are widely affected by the grid strength and X/R ratio of grid internal impedance. The flicker emission is decreased with higher fault levels. The risk of voltage fluctuation increases in the resistive resistive /re·sis·tive/ (rezis´tiv) pertaining to or characterized by resistance. grids. The wind turbine operating point and the QP characteristic of the generator determine the point of minimum flicker emission. The trend of flicker variation with the mean wind speed depends mainly on the wind turbine power curve. The power variability and consequently flicker emission increases with turbulence increase. A wide look on the results indicates that grid parameters have more effect on flicker emission than site parameters. Appendix The Wind Turbine Data, Drive Train Data and Generator Data are stated as follows: rated power 180 kW; hub height 30 m; rotor diameter 23.2 m; number of blades three; rotor speed 42 r/min; blade profile NACA63 200; gearbox ratio 23.75. stall regulated. turbine inertia 102.8 [kg.m.sup.2]; generator inertia 4.5 [kg.m.sup.2]; stiffness of the shaft 2700 N.m/rad; (all data referred to the high speed shaft). nominal voltage 400 V; number of polepairs three; stator resistance 0.0092 Ohm; rotor resistance (referred to the stator) 0.0061 Ohm; stator leakage inductance Leakage inductance is that property of an electrical transformer that causes a winding to appear to have some selfinductance in series with the mutuallycoupled transformer windings. 0.186 mH; rotor leakage inductance (referred to the stator) 0.427 mH; magnetizing inductance inductance, quantity that measures the electromagnetic induction of an electric circuit component; it is a property of the component itself rather than of the circuit as a whole. 6.7 mH. References [1] IEC 6140021, "Wind turbine generator systemsPart 21: Measurement and assessment of power quality characteristics of grid connected wind turbines" first edition, 2001. [2] D. S. L. Dolan and P. W. Lehn, "Realtime wind turbine emulator Hardware, software or a combination of the two that enables a computer to act like another computer and run applications written for that computer. In the past, it was often a hardware addon that actually contained an instruction execution module for the emulated computer. suitable for power quality and dynamic control studies" proceedings of the International Conference on Power Systems Transients (IPST'05), Montreal, Canada, June 1923, 2005. [3] M.I. Marei, E.F. ElSaadany and M.M.A. Salama, "Estimation techniques for voltage flicker envelope tracking", Electric Power System Research, Vol. 70, pp. 3037, 2004. [4] M. A. H. ElSayed, "Voltage level evaluation in distribution systems with integrated wind generation", Mepcon' 8th Conf., University of Helwan, Cairo, Egypt, 2931 December 2001. [5] J. O. Tande, "Applying power quality characteristics of wind turbines for assessing impact on voltage quality", Wind Energy, Vol. 5, Issue 1, pp. 3752, 2002. [6] C. Vilar, J. Usaola and H. Amaris, "A frequency domain approach to wind turbines for flicker analysis", IEEE (Institute of Electrical and Electronics Engineers, New York, www.ieee.org) A membership organization that includes engineers, scientists and students in electronics and allied fields. Trans. on Energy Conversion, Vol. 18, No. 2, pp. 335341, June 2003. [7] S. Kelouwani and K. Agbossou, " Nonlinear model identification of wind turbine with a neural network", IEEE Trans. on Energy Conversion, Vol. 19, No. 3, September 2004. [8] T. Petru and T. Thiringer, "Modeling of wind turbines for power system studies", IEEE Trans. on Power systems, Vol. 17, No. 4, November 2002. [9] P. Soerensen, "Frequency domain modeling of wind turbines structures", RisoR749 (EN), 1994. [10] S.A. Papathanassiou and M.P. Papadopoulos, "Mechanical stresses in fixed speed wind turbine due to network disturbances", IEEE Trans. on Energy Conversion, Vol. 16, No. 4, pp. 361367, December 2001. [11] P. Ledesma, J. Usaola, and J. L. Rodriguez, "Models of WECS WECS Wind Energy Conversion System WECS Workshop on Education in Computer Security WECS Wärtsilä Engine Control System WECS Wymondham & Edmondthorpe Civic Society (UK) WECS Web Enabled Collaboration System for power system dynamic studies," Proceedings of the UPEC'98, Edinburgh, 1998. [12] F. Iov, A. D. Hansen, C. Jauch, P. Soerensen, and F. Blaabjerg, "Advanced tools for modeling, design and optimization of wind turbine systems", Nordic Wind Power Conf., Chalmers University of Technology (body, education) Chalmers University of Technology  A Swedish university founded in 1829 offering master of science and doctoral degrees. Research is carried out in the main engineering sciences as well as in technology related mathematical and natural sciences. , Sweden, March 2004. [13] Z. Lubosny, K. Dobrzynski, J. Klucznik and R. Zajczyk, "Wind farm equivalenting for electric power system operation analysis", Nordic Wind Power Conf., Chalmers University of Technology, Sweden, March 2004. [14] IEC 61000415, "Electromagnetic Compatibility (hardware, testing) Electromagnetic Compatibility  (EMC) The extent to which a piece of hardware will tolerate electrical interference from other equipment, and will interfere with other equipment. (EMC (1) (EMC Corporation, Hopkinton, MA, www.emc.com) The leading supplier of storage products for midrange computers and mainframes. Founded in 1979 by Richard J. Egan and Roger Marino, EMC has developed advanced storage and retrieval technologies for the world's largest companies. )Part 4: Testing and measurements techniques Section 15: Flickermeter, Functional and design specifications", first edition, 1997. [15] P. Rosas, "Dynamic influences of wind power on the power system", Riso R1408 (EN), March 2003. [16] N. K. Bansal, M. Kleemann and M Meliss, "Renewable energy Renewable energy utilizes natural resources such as sunlight, wind, tides and geothermal heat, which are naturally replenished. Renewable energy technologies range from solar power, wind power, and hydroelectricity to biomass and biofuels for transportation. sources and conversion technology", book, Tata McGrawHill, 1990. [17] P. Soerensen A. D. Hansen and P. Rosas, " Wind models for simulation of power fluctuations from wind farms", Journal of Wind Engineering, Vol. 90, pp. 13811402, 2002. [18] M. P. Papadopoulos, S.A. Papathanassiou, S.T. Tentzerakis, "Modelling of induction machine main flux saturation saturation, of an organic compound saturation, of an organic compound, condition occurring when its molecules contain no double or triple bonds and thus cannot undergo addition reactions. in the arbitrary reference frameA case study on wind turbine selfexcitation", Proceedings of ICEM'96, Vigo, Spain, September 1996. [19] A. Feijoo, J. Cidras and C. Carrillo, "A third order model for the doublyfed induction machine", Electric Power Systems Research, Vol. 56, pp. 121127, 2000. [20] J. B. Ekanayake, L. Holdsworth and N. Jenkins, "Comparison of 5th order and 3rd order machine models for doubly fed induction generator (DFIG DFIG DoublyFed Induction Generator ) wind turbines", Electric Power Systems Research, Vol. 67, pp. 207 215, 2003. [21] IEC 6100033, "Electromagnetic Compatibility (EMC)Part 3: LimitsSection 3: Limitation of voltage fluctuations and flicker in lowvoltage supply systems for equipment with rated current [less than or equal to] 16 A." l" edition 1994, amendment 2001. [22] S. A. Papathanassiou, S.J. Kiartzis, M.P. Papadopoulos, A.G. Kladas, "Wind turbine flicker calculation using neural networks", Wind Engineering, Vol. 24, No. 5, pp. 317335, September 2000. [23] P. Sorensen, J. O. Tande, L. M. Sondergaard and J. D. Kledal, "Flicker emission levels from wind turbines" Wind Engineering, Vol. 20, No. 1, pp. 3946, 1996. [24] T. Thiringer, "Power quality measurements performed on a lowvoltage grid equipped with two wind turbines", IEEE Trans. on Energy Conversion, Vol. 11, No. 3, pp. 601606, September 1996. [25] P. Sorensen, G. Gerdes, F. Santjer, N. Robertson, W. Davy, M. Koulouvari, E. Morfiadakis, A. Larsson, "Standards for measurements and testing of wind turbine power quality" European Wind Energy Conf., Nice, France, pp. 721724, 15 March 1999. Hassan H. ElTamaly *, Mohamed A. A. Wahab and Ali H. Kasem Faculty of Engineering, ELMinia University, ELMinia 61111, EGYPT. * Corresponding author. Address: Faculty of Engineering ELMinia University, ELMinia 61111, EGYPT Email: dr_h_tamaly@yahoo.com 

Reader Opinion