17.08.2013 Views

Publishers version - DTU Orbit

Publishers version - DTU Orbit

Publishers version - DTU Orbit

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

To further simplify the notation we define the translation operator Tδ acting on any scalar or<br />

vector field ξ(x),<br />

Tδξ(x) = ξ(x+δ). (311)<br />

We also define a convolution operator Cn acting on any scalar or vector field as,<br />

Cnv(x) =<br />

∞<br />

−∞<br />

For the North direction, Eq. (310) can be written as,<br />

ϕ(s)n·v(x+ns) ds. (312)<br />

˜vrN = CnNTdfnNv. (313)<br />

We get similar expressions for South, East and West directions. Eq. (308) can then be written<br />

as,<br />

uwc = 1<br />

2sinφ [(CnNTdfnN −CnSTdfnS)cosΘ+(CnETdfnE −CnWTdfnW)sinΘ]v (314)<br />

We also know that by definition,<br />

〈u ′2 <br />

〉 = 〈û(k)û ∗ (k)〉dk, (315)<br />

whereˆdenotes Fourier transform and ∗ denotes complex conjugation. In the Fourier space<br />

we have,<br />

Tδv(k) = e ik·δ ˆv(k), (316)<br />

Cnv(k) = ˆϕ(n·k)n· ˆv(k), (317)<br />

where ˆϕ(k) = sinc 2 (klp/2), considering that the weighting function for a pulsed lidar is<br />

commonly defined as, <br />

lp−|s|<br />

l ϕ(s) =<br />

2 p<br />

0<br />

for |s| < lp;<br />

elsewhere,<br />

(318)<br />

where lp is the half length of the ideally rectangular light pulse leaving the lidar assuming the<br />

matching time windowing (= 2lp/c, where c is the speed of light). Thus in Fourier space Eq.<br />

(314) can then be written as,<br />

ûwc(k) = 1<br />

2sinφ [(nNe idfk·nN sinc 2 (k·nNlp/2)−nSe idfk·nS sinc 2 (k·nSlp/2))cosΘ<br />

+(nEe idfk·nE sinc 2 (k·nElp/2)−nWe idfk·nW sinc 2 (k·nWlp/2))sinΘ]· ˆv(k)<br />

≡ b(k)· ˆv(k),<br />

and the variance (from Eq. 315),<br />

〈u ′2 wc 〉 =<br />

<br />

(319)<br />

Φij(k)bi(k)b ∗ j (k) dk, (320)<br />

where we have implicitly used the relation, Φij(k) = 〈ˆvi(k)ˆv ∗ j (k)〉. The (co-) variances of<br />

other components can be estimated in a similar manner by first estimating the corresponding<br />

weighting functions ci(k) and ai(k) for vwc and wwc components respectively.<br />

13.2.3 Definition of the systematic error<br />

For simplicity we define systematic error as the ratio of the lidar second-order moment to<br />

the true second-order moment. Thus a ratio equal to one would signify no systematic error,<br />

whereas deviations from unity signify systematic error. By definition, the true second-order<br />

moment of a velocity component is given as,<br />

〈v ′ iv ′ <br />

j〉 = Φij(k)dk. (321)<br />

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

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

Saved successfully!

Ooh no, something went wrong!