The UMIST-N Near-Wall Treatment Applied to Periodic Channel Flow
The UMIST-N Near-Wall Treatment Applied to Periodic Channel Flow
The UMIST-N Near-Wall Treatment Applied to Periodic Channel Flow
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 2. TURBULENCE MODELS 24<br />
ˆε is required in order <strong>to</strong> calculate ε for use in the k transport equation<br />
(Equation 2.16). ˆε is defined as<br />
ˆε = 2ν<br />
<br />
∂ √ 2 k<br />
∂y<br />
(2.30)<br />
Thus, ˆε does not require a transport equation, but may be solved from local<br />
quantities.<br />
νt is modified according <strong>to</strong><br />
where<br />
2.2.2 Yap Correction<br />
νt = Cµfµ<br />
fµ = exp<br />
⎡<br />
2 k<br />
˜ε<br />
⎣ −2.5<br />
<br />
1 + ˜ Ret<br />
50<br />
⎤<br />
(2.31)<br />
⎦<br />
(2.32)<br />
<strong>The</strong> ‘Yap correction’ refers <strong>to</strong> an additional source term in the ˜ε transport<br />
equation of the low-Reynolds-number k-ε model. <strong>The</strong> Yap correction was<br />
originally introduced by Yap [92] <strong>to</strong> improve the performance of the k-ε<br />
model in impinging and recirculating flows. Yap correction is often implicitly<br />
included in standard k-ε modelling, so it has been included in the present<br />
<strong>UMIST</strong>-N wall function for completeness.<br />
Using Yap correction, the ˜ε transport equation becomes<br />
2<br />
D˜ε ν + νt<br />
˜ε<br />
˜ε<br />
= ▽ ·<br />
▽ ˜ε + Cε1f1 Pk − Cε2f2 + E + Y (2.33)<br />
Dt σε<br />
k<br />
k<br />
with the additional Y source term defined as<br />
⎛⎡<br />
<br />
Y = max ⎝⎣0.83<br />
k 3<br />
2<br />
− 1<br />
2.5˜εy<br />
k 3<br />
2<br />
2.5˜εy<br />
2 2 ˜ε<br />
k<br />
⎤ ⎞<br />
⎦ , 0⎠<br />
(2.34)