16.08.2013 Views

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

SHOW MORE
SHOW LESS

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)

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

Saved successfully!

Ooh no, something went wrong!