27.02.2013 Views

Wind Energy

Wind Energy

Wind Energy

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

46<br />

Stability of the Tip Vortices in the Far Wake<br />

behind a <strong>Wind</strong> Turbine<br />

Valery L. Okulov and Jens N. Sørensen<br />

Summary. The work is a further development of a model developed by one of the<br />

authors Okulov (J. Fluid Mech. 521:319–342 2004) in which linear stability of N<br />

equally azimuthally spaced infinity long helical vortices with constant pitch were<br />

considered. A multiplicity of helical vortices approximates the tip vortices of the far<br />

wake behind a wind turbine. The present analysis is extended to include an assigned<br />

vorticity field due to root vortices and the hub of the wake. Thus the tip vortices are<br />

assumed to be embedded in an axisymmetric helical vortex field formed from the<br />

circulation of the inner part of the rotor blades and the hub. As examples of inner<br />

vortex fields we consider three generic vorticity distributions (Rankine, Gaussian<br />

and Scully vortices) at radial extents ranging from the core radius of a tip vortex to<br />

several rotor radii.<br />

46.1 Theory: Analysis of the Stability<br />

The influence of an assigned flow on the stability of multiple vortices has<br />

so far only been studied for circular arrays of point vortices or straight vortex<br />

filaments (the limiting case of a helical vortex with infinite pitch) [2, 3].<br />

A non-perturbed multiple of helical vortices (Fig. 46.1a rotates with constant<br />

angular velocity Ω = Ω0 + ΩInd + ΩSind and moves uniformly along<br />

its axis with velocity V = V0 + VInd + VSind. Note that in contrast to the<br />

plane case, where the total vortex motion consists of an assigned flow field<br />

(Ω0,V0) and the mutual induction of the vortices (ΩInd,VInd) with VInd =0,<br />

we also include self-induction of each of helical vortex (ΩSind,VSind) [1].<br />

From linear stability, introducing infinitesimal space displacements rk =<br />

α (t)exp(2πiks/N), χk = β (t)exp(2πiks/N), of the k-vortex from the its<br />

equilibrium position(a, 2πk/N + Ωt,V t) a correlation equation B =0was<br />

derived for dimensionless pitch τ = l/a and vortex radius ɛ = rcore/a

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!