# Moment of Inertia Dependence of Vertical Axis Wind Turbines in Pulsating Winds

Authors: Yutaka Hara; Koichi Hara; Tsutomu Hayashi
Journal: International Journal of Rotating Machinery
 (2012)

---

## Abstract

Vertical Axis Wind Turbines (VAWTs) are unaffected by changes in wind direction, and they have a simple structure and the potential for high efficiency due to their lift driving force. However, VAWTs are affected by changes in wind speed, owing to effects originating from the moment of inertia. In this study, changes in the rotational speed of a small VAWT in pulsating wind, generated by an unsteady wind tunnel, are investigated by varying the wind cycle and amplitude parameters. It is shown that the responses observed experimentally agree with simulations based on torque characteristics obtained under steady rotational conditions. Additionally, a simple equation expressing the relationship between the rotational change width and amplitude of the pulsating wind is presented. The energy efficiency in a pulsating wind remains constant with changes in both the moment of inertia and the wind cycle; however, the energy efficiency decreases when the wind amplitude is large.


---

## Body

## 1. Introduction

Vertical Axis Wind Turbines (VAWTs) [1] are promising energy devices [2–4]. They have a simple structure, their output power is unaffected by changes in wind direction, and based on their lift driving force, they have the potential for high efficiency. However, VAWTs are affected by wind speed variation, owing to their moment of inertia. As a result, rotor response to wind speed variation is important when considering the proper control of VAWTs.In general, for the design, analysis, and control of the more widely prevailing Horizontal Axis Wind Turbines (HAWTs), prediction of unsteady aerodynamic loads is challenging and crucial. For this reason, much effort has been put into modeling the unsteady aerodynamics of wind turbines [5]. In the case of VAWTs, even if the wind speed is constant, the angle of attack of the inflow to a blade periodically changes during rotation. Thus, dynamic analysis of a VAWT is essential for reliable design and reduction of cost [6]. In practice, the problems are even more complex, as the wind turbines are often installed under conditions of turbulence and abrupt wind-speed variations.Ordinarily, wind turbine performance predictions are compared with field or experimental data measured under constant wind speed. Recently, transient responses were measured for a straight-bladed VAWT subjected to steplike wind speed variation in an unsteady wind tunnel [7]. In the present study, to investigate the dependence of VAWT performance on the moment of inertia in unsteady wind, the response of a straight-bladed VAWT to a pulsating wind is measured, changing the rotor moment of inertia, the wind cycle, and the wind amplitude. The experimental results are compared with a numerical simulation, in which it is assumed that the experimental torque coefficients are obtained under conditions of constant wind speed. Furthermore, a simulation is carried out for a theoretical VAWT with a different rotor size from that used experimentally, and the dependence of the rotor response on the mean speed of the pulsating wind is investigated, as is the behavior of the rotor in the case of a long wind cycle.


## 2. Experimental Apparatus

A schematic of the experimental setup is shown in Figure1. The experimental VAWT has four straight blades with NACA 0012 airfoils (chord length: c=0.075 m). The rotor diameter is D=2R=0.6 m, and the height is H=0.47 m. The swept area of the rotor is  A=D×H=0.282 m2. The original moment of inertia of the rotor is I=0.1 kgm2, and this can be increased to up to 0.3 kgm2 by fixing weights on the axis under the rotor. Figure 2 displays a photograph of the VAWT rotor and the weights mounted upon the rotational axis. An Eiffel-type wind tunnel at Tottori University (square nozzle outlet: 1.5m×1.5m) can generate pulsating winds by means of an axial blower that can change its blade angle during constant speed rotation [8]. In the present experiment, a powder brake is manipulated to generate a constant load torque, TL. The distance between the wind tunnel outlet and the rotor axis is 1.5 m, and the wind speed, V, is measured with a hot wire probe placed at the nozzle outlet.Figure 1
Experimental VAWT setup.Figure 2
Experimental VAWT rotor with weights added to rotational axis to increase the moment of inertia.


## 3. Experimental Method

The mean speed of the pulsating wind was fixed atV0=10 m/s, and the amplitude was fixed at ΔV=1 or 2 m/s. The wind cycle was changed from T=2 to 8 s. The load torque was held at around TL=0.3 Nm, consistent with the maximum power condition of the experimental VAWT at 10 m/s constant wind speed. The sampling frequency of the rotational and wind speeds was 8 Hz. The recording time for the torque, rotational speed, and wind speed depended on the pulsating wind cycle and was set to between approximately 200 and 250 s. This corresponded to about 50 cycles for T=4 s, 42 cycles for T=6 s, and 31 cycles for T=8 s. All of the measured data were phase-averaged, based on the measured wind speed data, and an example of the averaged data for the wind and rotational speeds is given in Figure 3. Several time series datasets were obtained for each experimental condition, and the output for specific characteristics, such as time delay, was averaged such that a representative value could be found for each condition.Figure 3
Response of rotational speed to pulsating wind (I=0.1 kgm2, T=3 s, ΔV = 1 m/s).


## 4. Simulation

The behavior of angular velocity,ω, in response to the pulsating wind was simulated by using a Runge-Kutta method to numerically integrate the motion equation of the wind turbine:(1)Idωdt=Q-TL,
where(2)Q=TW-TB,(3)TW=12ρCq(λ)V2AR,(4)λ=RωV.
In (1), I is the moment of inertia for the wind turbine rotor only, and the moments of inertia for the torque detector and powder brake are regarded as insignificant. Q is the effective torque and is defined as the difference between the wind turbine torque, TW, and the braking torque, TB, coming from the bearings (2). TL is the load torque, which is generated by the powder brake in this case. The wind turbine torque, TW, is defined via (3), where ρ is air density.In the simulation, the torque coefficientCq(λ), which is dependent on the tip speed ratio,λ, and wind speed, V, was given by interpolating tabulated data of this coefficient measured under a constant wind speed. Figure 4 presents the torque coefficients of the VAWT used in this study. These torque characteristics were measured without any weights placed upon the rotor axis; however, it was shown that the addition of weights did not result in a significant change to Cq(λ). Finally, it was assumed that the flow speed at the rotor was the same as the upstream main flow when the latter was set as a sinusoidal wave with an average wind speed of 10 m/s.Figure 4
Torque coefficients,Cq(λ), for a VAWT with NACA 0012 airfoils.


## 5.  Results and Discussion

### 5.1. Time Delay and Rotational Speed Width

An example of the experimental results is shown in Figure3, where the phase-averaged wind and rotational speeds for about 50 periods are graphically plotted for a case where ΔV=1 m/s, T=3 s, and I=0.1 kgm2. It can be seen that the rotational speed varies periodically with a constant time-delay,τ, to the wind speed variation. For the case shown in Figure 3, the time-averaged rotational speed is Nave=620.5 rpm, and the rotational speed width (p-p) is ΔN=6.3 rpm.The data in Figure5(a) are the time delays for the case where ΔV=1 m/s. These delays were corrected by uniformly subtracting 0.15 s to account for the wind speed being measured 1.5 m upstream from the rotor center. Additionally, the delay in the wind variation resulting from the axial distance was measured for the unsteady wind generated by the wind tunnel described in Section 2 [9]. Figure 5(a) shows that the time delay increases with the wind cycle, but stays constant for increases in the moment of inertia. The time delays show a similar dependence on the wind cycle, T, and the moment of inertia, I, for the case when ΔV = 2 m/s in Figure 5(b).Time delay,τ, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)The experimental error for the time delays,τ(s), measured under each experimental condition, is shown in Figures 6(a) and 6(b). The experimental error is defined as the difference between the maximum and the minimum values found for the 3–5 measurements of τ expressed as a percentage of the average value of the measurements. The experimental errors were 15.8% when ΔV= 1 m/s and 12.9% when ΔV = 2 m/s.Experimental error of time delay for wind amplitude (a)ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)The dependence of the rotational speed width,ΔN, on both the moment of inertia and wind cycle, for ΔV = 1 m/s, is given in Figure 7(a) and shows that the speed width, ΔN, increases with wind cycle, T, and decreases with the moment of inertia, I. When ΔV = 2 m/s, the dependence of ΔN on T and I is nearly identical to that when ΔV = 1 m/s, although the values of ΔN are almost double (Figure 7(b)). As above, the experimental error for the rotational speed width, ΔN (rpm), measured for each experimental condition, are shown in Figures 8(a) and 8(b), and the experimental errors for ΔN are 12.6% when ΔV = 1 m/s and 11.2% when ΔV = 2 m/s.Rotational speed width,ΔN, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)Experimental error of rotational speed width for wind amplitude (a)ΔV = 1 and (b) ΔV= 2 m/s.
(a)(b)


### 5.2. Phase Delay

Phase delays,ϕ(=τ/T×2π), between the rotational speed variation and wind variation are displayed in Figures 9(a) and 9(b), for the cases where ΔV = 1 m/s and ΔV = 2 m/s, respectively. For the former, the phase delay is increased for short wind cycles, T. Conversely, when ΔV= 2 m/s, the phase delay remains nearly constant throughout. The average simulated phase delay, represented by the broken line in Figure 9, is found to be approximately π/2in both cases.Phase delay,ϕ, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)


### 5.3.  Rate of Change of Rotational Speed

The rate of change of the rotational speed divided by the wind cycle,(ΔN/Nave)/T, is shown for both wind amplitudes in Figures 10(a) and 10(b). For each ΔV case, the values of (ΔN/Nave)/T are expressed by a curve inversely proportional to the moment of inertia, I, and independent of the wind cycle, T. For ΔV = 2 m/s, a difference exists between the experimental results and those from the simulation depicted by the broken curve. Nevertheless, for both cases of ΔV, the experimental and simulated outputs largely conform to each other.Rate of change of rotational speed divided by wind cycle,(ΔN/Nave)/T, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)Comparing coefficients of the two simulated curves, the coefficient whenΔV = 2 m/s is around twice as large as that for ΔV = 1 m/s. Therefore, the rotational speed width may be expressed by(5)ΔN≈C(NaveI)TΔV,
where C is a constant and takes a value of around 6.8×10-4 for the wind turbine used here.


### 5.4. Energy Efficiency

The experimental energy efficiency,η(%), of the wind turbine is shown in Figures 11(a) and 11(b), together with simulated results depicted by broken lines. Figure 12 presents the experimental error of η for each experimental condition. For the case when ΔV = 1 m/s, the experimental energy efficiency error is 4.0%, and when ΔV = 2 m/s, the error is 5.1%. In this study, η is defined as the ratio between the mechanical power, PT, of the turbine and the input wind power, PW, expressed as a percentage:(6)η=100(PTPW),
where(7)PT=1T∑0TωTLΔt,PW=1T∑0T12ρAV3Δt.
In (7), Δt is the data sampling time interval and  ∑ indicates summation over a wind cycle.Energy efficiency,η, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)Experimental error of energy efficiency for wind amplitude (a)ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)From Figures11(a) and 11(b), the energy efficiency tends to decrease with increasing moment of inertia, I. However, this tendency was caused by an increase in the braking torque of the bearings, which affected by the weights added to change the moment of inertia. For both cases of ΔV, the simulated energy efficiency is not dependent on I or T, but the energy efficiency does decrease when the wind amplitude is large. This decrease in efficiency results from the fact that, although raising the wind amplitude leads to increased input energy, the average power of the wind turbine remains unchanged, or even decreases slightly, when the wind amplitude is altered. Figure 13 shows the loci of the simulated mechanical power, PT, for the wind turbine, in the case T = 4 s. The average powers (20.6 W for ΔV = 1 m/s, 19.8 W for ΔV = 2 m/s) are almost the same, regardless of ΔV, whereas the input wind power, PW, increases by 4.7% when ΔV is doubled (171 W for ΔV = 1 m/s, 179 W for ΔV = 2 m/s).Figure 13
Loci of simulated mechanical power of the wind turbine in the caseT = 4 s. (Note: power curves of the wind turbine measured at five wind speeds from 8–12 m/s are also shown.)


## 5.1. Time Delay and Rotational Speed Width

An example of the experimental results is shown in Figure3, where the phase-averaged wind and rotational speeds for about 50 periods are graphically plotted for a case where ΔV=1 m/s, T=3 s, and I=0.1 kgm2. It can be seen that the rotational speed varies periodically with a constant time-delay,τ, to the wind speed variation. For the case shown in Figure 3, the time-averaged rotational speed is Nave=620.5 rpm, and the rotational speed width (p-p) is ΔN=6.3 rpm.The data in Figure5(a) are the time delays for the case where ΔV=1 m/s. These delays were corrected by uniformly subtracting 0.15 s to account for the wind speed being measured 1.5 m upstream from the rotor center. Additionally, the delay in the wind variation resulting from the axial distance was measured for the unsteady wind generated by the wind tunnel described in Section 2 [9]. Figure 5(a) shows that the time delay increases with the wind cycle, but stays constant for increases in the moment of inertia. The time delays show a similar dependence on the wind cycle, T, and the moment of inertia, I, for the case when ΔV = 2 m/s in Figure 5(b).Time delay,τ, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)The experimental error for the time delays,τ(s), measured under each experimental condition, is shown in Figures 6(a) and 6(b). The experimental error is defined as the difference between the maximum and the minimum values found for the 3–5 measurements of τ expressed as a percentage of the average value of the measurements. The experimental errors were 15.8% when ΔV= 1 m/s and 12.9% when ΔV = 2 m/s.Experimental error of time delay for wind amplitude (a)ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)The dependence of the rotational speed width,ΔN, on both the moment of inertia and wind cycle, for ΔV = 1 m/s, is given in Figure 7(a) and shows that the speed width, ΔN, increases with wind cycle, T, and decreases with the moment of inertia, I. When ΔV = 2 m/s, the dependence of ΔN on T and I is nearly identical to that when ΔV = 1 m/s, although the values of ΔN are almost double (Figure 7(b)). As above, the experimental error for the rotational speed width, ΔN (rpm), measured for each experimental condition, are shown in Figures 8(a) and 8(b), and the experimental errors for ΔN are 12.6% when ΔV = 1 m/s and 11.2% when ΔV = 2 m/s.Rotational speed width,ΔN, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)Experimental error of rotational speed width for wind amplitude (a)ΔV = 1 and (b) ΔV= 2 m/s.
(a)(b)


## 5.2. Phase Delay

Phase delays,ϕ(=τ/T×2π), between the rotational speed variation and wind variation are displayed in Figures 9(a) and 9(b), for the cases where ΔV = 1 m/s and ΔV = 2 m/s, respectively. For the former, the phase delay is increased for short wind cycles, T. Conversely, when ΔV= 2 m/s, the phase delay remains nearly constant throughout. The average simulated phase delay, represented by the broken line in Figure 9, is found to be approximately π/2in both cases.Phase delay,ϕ, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)


## 5.3.  Rate of Change of Rotational Speed

The rate of change of the rotational speed divided by the wind cycle,(ΔN/Nave)/T, is shown for both wind amplitudes in Figures 10(a) and 10(b). For each ΔV case, the values of (ΔN/Nave)/T are expressed by a curve inversely proportional to the moment of inertia, I, and independent of the wind cycle, T. For ΔV = 2 m/s, a difference exists between the experimental results and those from the simulation depicted by the broken curve. Nevertheless, for both cases of ΔV, the experimental and simulated outputs largely conform to each other.Rate of change of rotational speed divided by wind cycle,(ΔN/Nave)/T, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)Comparing coefficients of the two simulated curves, the coefficient whenΔV = 2 m/s is around twice as large as that for ΔV = 1 m/s. Therefore, the rotational speed width may be expressed by(5)ΔN≈C(NaveI)TΔV,
where C is a constant and takes a value of around 6.8×10-4 for the wind turbine used here.


## 5.4. Energy Efficiency

The experimental energy efficiency,η(%), of the wind turbine is shown in Figures 11(a) and 11(b), together with simulated results depicted by broken lines. Figure 12 presents the experimental error of η for each experimental condition. For the case when ΔV = 1 m/s, the experimental energy efficiency error is 4.0%, and when ΔV = 2 m/s, the error is 5.1%. In this study, η is defined as the ratio between the mechanical power, PT, of the turbine and the input wind power, PW, expressed as a percentage:(6)η=100(PTPW),
where(7)PT=1T∑0TωTLΔt,PW=1T∑0T12ρAV3Δt.
In (7), Δt is the data sampling time interval and  ∑ indicates summation over a wind cycle.Energy efficiency,η, for wind amplitude (a) ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)Experimental error of energy efficiency for wind amplitude (a)ΔV = 1 and (b) ΔV = 2 m/s.
(a)(b)From Figures11(a) and 11(b), the energy efficiency tends to decrease with increasing moment of inertia, I. However, this tendency was caused by an increase in the braking torque of the bearings, which affected by the weights added to change the moment of inertia. For both cases of ΔV, the simulated energy efficiency is not dependent on I or T, but the energy efficiency does decrease when the wind amplitude is large. This decrease in efficiency results from the fact that, although raising the wind amplitude leads to increased input energy, the average power of the wind turbine remains unchanged, or even decreases slightly, when the wind amplitude is altered. Figure 13 shows the loci of the simulated mechanical power, PT, for the wind turbine, in the case T = 4 s. The average powers (20.6 W for ΔV = 1 m/s, 19.8 W for ΔV = 2 m/s) are almost the same, regardless of ΔV, whereas the input wind power, PW, increases by 4.7% when ΔV is doubled (171 W for ΔV = 1 m/s, 179 W for ΔV = 2 m/s).Figure 13
Loci of simulated mechanical power of the wind turbine in the caseT = 4 s. (Note: power curves of the wind turbine measured at five wind speeds from 8–12 m/s are also shown.)


## 6. Numerical Analysis of a Theoretical VAWT

In the previous section, numerical simulations showed good agreement with experimental results. In this section, similar simulations are reported for a theoretical VAWT of different size with different airfoils. This new VAWT has four straight blades with NACA 0018 airfoils (c=0.15m),  D=2m, and H=1.5 m. The swept area (A=3 m2) is about 10-fold that of the experiment rotor (0.282 m2). I of the new VAWT is estimated to be about 5 kgm2 if the blades are made of a lightweight material with a density of 300 kg/m3. Power and torque were calculated by using blade element momentum theory [1, 10–12]. Simulated power (Cp) and torque (Cq) coefficients are presented in Figure 14 for nine wind speeds from 6 to 14 m/s. In these performance calculations, the double-multiple streamtube model was used for flow field modeling and the modified Gormont model [1] was applied to consider the effect of dynamic stall. The parameter AM in the dynamic stall model was set to 1000 in this study. Aerodynamic data (lift coefficient, CL, and drag coefficient, CD) needed for the performance calculation were taken from [13] for Re=1.0×104-1.6×105 and from [14] for Re>1.6×105. Two peaks appear in the Cq-λ curves at each constant wind speed in Figure 14(b). The right-hand peak corresponds to the maximum torque based on the static aerodynamic data, and the left-hand peak is attributed to the dynamic stall effects.Simulated power (Cp) and torque (Cq) coefficients for a theoretical VAWT with NACA 0018 airfoil.
(a)(b)To simulate the response of the theoretical VAWT, the amplitude of the pulsating wind was fixed atΔV = 1 m/s only, whereas the mean speed of the pulsating wind was fixed at  V0=7,10, and 13 m/s, and the wind cycle at T=4,8,16, and 32 s. Simulations were performed following the method described in Section 4, and  (1) was numerically integrated from an initial state, which was taken to be the maximum power state of the rotor at the mean wind speed. The calculation time was 300 s for a time interval of 0.01 s, and the steady-state data between 200 and 300 s were used in the analysis. Three torque curves were input as data into the simulation. For example, when V0 = 7 m/s, the torque coefficients,Cq(λ), for V=6,7, and 8 m/s shown in Figure 14 were used. Furthermore, for the simulations in this section, the braking torque, TB, was ignored (TB = 0; i.e., Q=TW).Phase delays,ϕ=(τ/T×2π), at each mean wind speed are displayed in Figure 15. In cases where V0 = 7 or 10 m/s and T=4,8, or 16 s all phase delays are approximately π/2. However, when V0 = 7 or 10 m/s and T=32 s, the phase delay is less thanπ/2 when the moment of inertia is small (I = 5 or 7 kgm2). When the mean wind speed is large (V0 = 13 m/s), this tendency becomes more remarkable.Phase delay,ϕ, for mean wind speed (a) V0 = 7, (b) V0 = 10, and (c) V0 = 13 m/s.
(a)(b)(c)Dependence of(ΔN/Nave)/Ton I is shown for each mean wind speed in Figure 16. At all mean wind speeds, when T=4,8, or 16 s, the values of (ΔN/Nave)/T are independent of T and match well with the approximation in (5). The curve constants C of the theoretical VAWT are 0.0508, 0.0499, and 0.0700 for V0 = 7, 10, and 13 m/s, respectively. Roughly speaking, the C increases in accordance with the order of the moment of inertia (C=6.8×10-4 for I=0.1 kgm2 of the experimental rotor; C≈5.0×10-2  for I≈5 kgm2 of the theoretical rotor). However, for a fixed rotor size and shape, the value of C shows a weak dependence on the mean speed of the pulsating wind. For T = 32 s in Figure 16, the values of  (ΔN/Nave)/T  tend to diverge from approximations given by (5) when I<10 kgm2. However, the divergence of the large wind speed case is different than those for small and medium wind speeds.Rate of change of rotational speed divided by wind cycle,(ΔN/Nave)/T, for mean wind speed (a) V0 = 7, (b) V0 = 10, and (c) V0 = 13 m/s.
(a)(b)(c)Energy efficiencies,η (%), are shown in Figure 17. The efficiency reflects the effect of the Reynolds number on the blade aerodynamic performance and η≈28% when V0 = 7 m/s, and η≈ 36% when V0 = 10 or 13 m/s. At all mean wind speeds, when T = 4, 8, or 16 s, η is almost independent of both T and I. However, in the case T= 32 s, efficiency varies when I is small. Moreover, for V0 = 10 m/s, η decreases with decreasing moment of inertia. In contrast, when V0 = 13 m/s, η increases abruptly when I<7 kgm2.Energy efficiency,η, for mean wind speed (a) V0 = 7, (b) V0 = 10, and (c) V0 = 13 m/s.
(a)(b)(c)To elucidate the above-mentioned behavior of the energy efficiency whenT is large and I is small, time profiles of the rotational speed, N, and the loci of the effective torque, Q, are next presented. The initial state of each profile and locus is the maximum power state of the mean wind speed.For the case whereV0 = 13 m/s, T= 32 s, and I = 7 kgm2, the time profile of N is shown in Figure 18 and the corresponding locus of Q is presented in Figure 19. Similarly, for the case where V0 = 13 m/s, T = 32 s, and I= 3 kgm2, the time profile and locus are shown in Figures 20 and 21, respectively. Comparing the cases where I = 7 kgm2(Figures 18 and 19) and I = 3 kgm2(Figures 20 and 21), a rapid change in the rotational speed due to the small moment of inertia (I= 3 kgm2) causes increases in the rotational speed width, ΔN, the average rotational speed, Nave (Nave = 353.9 and 377.1 rpm when I= 7 and 3 kgm2, resp.), and the energy efficiency, η.Figure 18
Time profile of the rotational speed of pulsating wind in the case whereV0 = 13 m/s, T = 32 s, and I= 7 kgm2.Figure 19
Locus of torque,Q, of pulsating wind in the case where V0 = 13 m/s, T = 32 s, and I = 7 kgm2. Constant load torque is TL=38.3 Nm. Final average rotational speed is Nave=353.9 rpm.Figure 20
Time profile of rotational speed of pulsating wind in the case whereV0 = 13 m/s, T = 32 s, and I = 3 kgm2.Figure 21
Locus of torque,Q, of pulsating wind in the case where V0 = 13 m/s, T = 32 s, and I = 3 kgm2. Constant load torque is TL = 38.3 Nm. Final average rotational speed is Nave = 377.1 rpm.For the case whereV0 = 10 m/s, I = 5 kgm2, and T = 16 s, the time profile of N is shown in Figure 22 and the corresponding locus of Q is presented in Figure 23, and when V0 = 10 m/s, I= 5 kgm2, and T = 32 s, the time profile and locus are shown in Figures 24 and 25. Comparing between these cases, a slow change in the wind speed due to the long cycle (T = 32 s) causes an increase in ΔN and a decrease in Nave (Nave = 299.4 and 296.2 rpm when T = 16 and 32 s, resp.). As a result, η is somewhat decreased.Figure 22
Time profile of rotational speed of pulsating wind in the case whereV0 = 10 m/s, T = 16 s, and I=5 kgm2.Figure 23
Locus of torque,Q, of pulsating wind in the case where V0 = 10 m/s, T = 16 s, and I = 5 kgm2. Constant load torque is TL = 20.8 Nm. Final average rotational speed is Nave=299.4 rpm.Figure 24
Time profile of rotational speed of pulsating wind in the case whereV0 = 10 m/s, T = 32 s, and I = 5 kgm2.Figure 25
Locus of torque,Q, of pulsating wind in the case where V0 = 10 m/s, T = 32 s, I = 5 kgm2. Constant load torque is TL = 20.8 Nm. Final average rotational speed is Nave=296.2 rpm.As shown in Figures21 and 25, for cases with large wind cycle and small moment of inertia, the torque locus is enlarged in the direction of the axis of the rotational speed and the energy efficiency can vary. However, the direction of the change in the energy efficiency appears to be dependent on the shape of torque curve or on the relation between T and I.Finally, as an example of an extreme case, the time profile ofN and the corresponding locus of Q when V0 = 10 m/s, T = 32 s, and I = 3 kgm2 are shown in Figures 26 and 27, respectively. As shown in Figure 27, with a long wind cycle and a small moment of inertia, the right-hand side of (1)—that is, the value Q-TL—takes negative values after the first peak in wind speed. Therefore, the rotational speed of the wind turbine decreases (dω/dt<0) and finally becomes zero. This unsustainable rotation condition is encountered even at larger values of I if T or ΔV becomes large. However, because the load torque of a wind power generator generally varies according to the rotational speed, the behavior of an actual wind turbine might be somewhat different from the present result.Figure 26
Time profile of rotational speed of pulsating wind in the case whereV0 = 10 m/s, T = 32 s, and I = 3 kgm2.Figure 27
Locus of torque,Q, of pulsating wind in the case where V0 = 10 m/s, T = 32 s, and I = 3 kgm2. Constant load torque is TL = 20.8 Nm.


## 7.  Conclusions

The response of a straight-bladed VAWT to pulsating wind was investigated by changing the moment of inertia, the wind cycle, and the wind amplitude. The phase delay between the rotational speed and wind variation held a constant value, becoming approximatelyπ/2when the wind amplitude was large. The rate of change in the rotational speed divided by the wind cycle, (ΔN/Nave)/T, was inversely proportional to the moment of inertia and independent of the wind cycle. This result suggests a relation between the rotational speed width and pulsating wind amplitude. The energy efficiency of the VAWT in a pulsating wind with unvarying amplitude stayed almost constant under both changes in the moment of inertia and in the wind cycle, but the energy efficiency decreased when the wind amplitude was large.Similarly, numerical simulation of a theoretical VAWT of a larger size than the experimental rotor showed that the energy efficiency was almost independent of the wind cycle and the moment of inertia under the condition of constant load torque. However, in the case of a long wind cycle and a small moment of inertia, the energy efficiency is expected to vary and is dependent on the torque curves.



---

*Source: 910940-2012-04-18.xml*