Wind field simulations at Askervein hill - WindSim
Wind field simulations at Askervein hill - WindSim
Wind field simulations at Askervein hill - WindSim
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 -