27.02.2013 Views

Wind Energy

Wind Energy

Wind Energy

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

28 J. Tambke et al.<br />

where u∗ is the friction velocity of the air flow and w∗ the one of the water<br />

flow while ρair and ρwater are the respective densities. This is equivalent<br />

to assuming that the shear stress is constant across the interface between air<br />

and water. The constant stress defines the wave boundary layer, extending<br />

from a height zB above to −zB below the water level. This wave boundary<br />

layer is similar to the surface layer onshore where the logarithmic wind profile<br />

is valid.<br />

Second, the inertial coupling relation [1] represents the shear stress in the<br />

wave boundary layer in the form of a drag law with regard to the inertially<br />

weighted fluid speeds at the limits of this layer, where KI is a specific drag<br />

coefficient: τwave = KI[ √ ρair u(zB) − √ ρwater u(−zB)] 2 .<br />

Third, analogous to the first similarity relation the turbulent viscosities<br />

νair and νwater of the two Ekman layers are also assumed to be weighted<br />

according to the inverse density ratio of the two fluids: νwater<br />

νair<br />

= ρair<br />

ρwater .<br />

These conditions allow the determination of the velocities of the two fluids<br />

at certain heights as an implicit function of a given geostrophic wind. The<br />

full profiles can be derived in the following way: Due to the constant shear<br />

stress the wind profile in the wave boundary layer has a logarithmic shape for<br />

neutral thermal stratification. Non-neutral conditions can be modelled with<br />

an additional stability function Ψm and the Monin–Obukhov (MO) length<br />

L. Expressed in a coordinate system where the horizontal component of the<br />

stress tensor τ h is parallel to the x-axis the profile can be written as<br />

�<br />

u(z) = uL + u∗<br />

κ ln<br />

� �<br />

z<br />

+ u∗<br />

κ Ψm<br />

�<br />

z<br />

� �<br />

,vL , (5.1)<br />

L<br />

zR<br />

and is valid for zR ≤ z ≤ zB; κ is the von Karman constant and zR is the<br />

theoretical height where the momentum transfer from the air to the wave<br />

field is centred. The above mentioned drift velocities are uair = u(zR) =uL<br />

and uwater = u(−zR) ≈ 1<br />

29uL. As a result of the above assumptions and the<br />

choice of the coordinate system the drift (uL,vL) at the very small height<br />

zR is collinear to the geostrophic wind: (uL,vL) = 1<br />

2<br />

G, where G =(ug,vg).<br />

This theoretical, though curious relation has no relevance in-situ, since neither<br />

(uL,vL) close to the waves nor G can be measured directly. Note that in the<br />

chosen coordinate system ug > 0andvg < 0.<br />

Concluding from [2], the relation that connects the friction velocity u∗ at<br />

the water surface to the geostrophic wind is given by<br />

|G| =<br />

√ r 2 +1<br />

|r +1|<br />

u∗<br />

√ . (5.2)<br />

KI<br />

The angle of rotation of the surface stress tensor to the left hand side<br />

of the geostrophic velocity (in the northern hemisphere) is atan(-1/r) with<br />

−1

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

Saved successfully!

Ooh no, something went wrong!