23.06.2015 Views

Wind field simulations at Askervein hill - WindSim

Wind field simulations at Askervein hill - WindSim

Wind field simulations at Askervein hill - WindSim

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

ε<br />

t<br />

= ε p<br />

Here<br />

by:<br />

K<br />

⎛<br />

⎜<br />

U<br />

=<br />

⎝ C<br />

β =<br />

τ<br />

µ<br />

C<br />

2<br />

⎞<br />

⎟ , Y is the normal distance of the first grid point from the wall, and β is given<br />

⎠<br />

µ<br />

( ε ) ( )<br />

−1<br />

p<br />

C<br />

p3<br />

K<br />

2<br />

C<br />

- C<br />

2<br />

p1<br />

- C<br />

p2<br />

which yields a value of β=0.25.<br />

Non equilibrium log-law wall functions<br />

This law uses the root of the kinetic energy<br />

characteristic length scale.<br />

Thus, equ<strong>at</strong>ion (1) is replaced by :<br />

k r<strong>at</strong>her than the friction velocity as the<br />

U k<br />

U<br />

2<br />

r<br />

2<br />

τ<br />

2<br />

Et<br />

k y<br />

= ⎛ ⎝ ⎜ ⎞<br />

ln ⎟<br />

κ ⎠<br />

t<br />

where<br />

κ<br />

t<br />

= κ C<br />

µ<br />

E<br />

EC<br />

0.<br />

25<br />

t<br />

=<br />

µ<br />

0 . 2 5<br />

The equ<strong>at</strong>ions (2) to (5) remain unchanged. The value of k <strong>at</strong> the near wall point is calcul<strong>at</strong>ed<br />

from its transport equ<strong>at</strong>ion, considering th<strong>at</strong> energy diffusion is set to zero to the wall.<br />

For the2-scale k-ε model, the tre<strong>at</strong>ment is similar to th<strong>at</strong> employed for the standard k-ε model,<br />

with k p and k t determined from their respective balance equ<strong>at</strong>ions and for simplicity it is<br />

assumed th<strong>at</strong> ε p =ε t .<br />

2.3.2.3. Other boundaries<br />

On the l<strong>at</strong>eral sides and on the top of the domain are imposed zero normal deriv<strong>at</strong>ives and<br />

zero values for velocities, k and ε. At the outlet, zero gradient boundary conditions are<br />

imposed, meaning th<strong>at</strong> plane geometry downwind the domain is assumed. Pressure is set to be<br />

equal to zero.<br />

- 14 -

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

Saved successfully!

Ooh no, something went wrong!