17.08.2013 Views

Publishers version - DTU Orbit

Publishers version - DTU Orbit

Publishers version - DTU Orbit

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.

where x is the center of the scanning circle and f(x) is any function of x. The corresponding<br />

spectral transfer function is given as,<br />

ˆTf(k1) = sinc 2<br />

<br />

k1Lf<br />

, (299)<br />

2<br />

where sinc(x) = sin(x)/x. The variances of uqq, vqq and wqq are given as,<br />

〈u ′2<br />

qq〉sin 2 <br />

φ = Φij(k)βi(k)β ∗ j(k) ˆ Tf(k1) dk, (300)<br />

〈v ′2<br />

qq〉sin 2 φ =<br />

〈w ′2<br />

qq〉cos 2 φ =<br />

<br />

<br />

Φij(k)γi(k)γ ∗ j(k) ˆ Tf(k1) dk, (301)<br />

Φij(k)αi(k)α ∗ j(k) ˆ Tf(k1) dk. (302)<br />

13.2.2 Systematic turbulence errors for the WindCube lidar<br />

The assumption made in section 13.2.1 that the mean wind direction comes from the North<br />

cannot be made for the WindCube, since it measures at four azimuth angles only (refer Fig.<br />

164), e.g. North, East, South and West. In this case the coordinate system is such that u is<br />

aligned in the mean wind direction. Thus,<br />

uwc = uNScosΘ+uEW sinΘ, (303)<br />

vwc = uNSsinΘ−uEW cosΘ, (304)<br />

whereuNS anduEW denotewindspeedsintheNorth-SouthandEast-Westdirectionsrespectively,<br />

Θ denotes the wind direction, and the subscript wc denotes the velocity components<br />

measured by the WindCube. From simple geometrical considerations (refer Fig. 164),<br />

uNS = ˜vrN − ˜vrS<br />

, (305)<br />

2sinφ<br />

uEW = ˜vrE − ˜vrW<br />

2sinφ<br />

, (306)<br />

where ˜vrN, ˜vrS, ˜vrE, ˜vrW are the weighted average radial velocities in the North, South,<br />

East and West directions respectively. For the w component,<br />

wwc = P(˜vrN + ˜vrS)+Q(˜vrE + ˜vrW)<br />

, (307)<br />

2cosφ<br />

where P and Q are the weights associated with the wind direction such that P + Q = 1.<br />

Leosphere uses P = cos 2 Θ and Q = sin 2 Θ, and hence, we use the same in our calculations.<br />

We proceed by deriving expressions for the uwc variance. The expressions for the (co-)<br />

variances of the remaining components of wind velocity can be derived in a similar manner.<br />

Substituting Eqs. (305), (306) into Eq. (303) we get,<br />

uwc = 1<br />

2sinφ [(˜vrN − ˜vrS)cosΘ+(˜vrE − ˜vrW)sinΘ]. (308)<br />

We define unit vectors in the four directions as,<br />

nN = n(−Θ),<br />

nS = n(π −Θ),<br />

nE = n( π<br />

2 −Θ),<br />

(309)<br />

nW = n( 3π<br />

2 −Θ),<br />

where nN, nS, nE and nW are the unit directional vectors in the North, South, East and<br />

West directions respectively. From Eq. (281), for the North direction,<br />

∞<br />

˜vrN =<br />

ϕ(s)nN ·v(snN +dfnN) ds. (310)<br />

−∞<br />

<strong>DTU</strong> Wind Energy-E-Report-0029(EN) 247

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

Saved successfully!

Ooh no, something went wrong!