08.11.2012 Views

Jegyzet kezdemény (BMEGEÁTMO10 ... - BME Áramlástan Tanszék

Jegyzet kezdemény (BMEGEÁTMO10 ... - BME Áramlástan Tanszék

Jegyzet kezdemény (BMEGEÁTMO10 ... - BME Áramlástan Tanszék

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

c○Copyright 2010 - Lohász Máté Márton, Régert Tamás<br />

12. FEJEZET. A RANS MODELLEZÉS 82<br />

közelíteni. Így a disszipációra pusztán analógia alapján a k egyenlethez hasonlóan<br />

írjuk föl:<br />

∂tε + uj ∂jε = C1εP ε ε<br />

� �<br />

νt<br />

− C2εε − ∂j ∂jε (12.17)<br />

k k σε<br />

A baloldalra írtunk lokális és konvektív változást, a jobb oldalon a produkciót<br />

és disszipációt ε<br />

k segítségével a helyes dimenzióra váltottuk át és C1ε-t és C2ε-t<br />

használjuk további korrekcióra. A transzport tagot szintén gradiens diffúziós hipotézissel<br />

írtuk σε Schmidt szám felhasználásával.<br />

A standard modell konstansai<br />

A modell egyenletek konstansaira az eredeti verzióban a következő értékeket javasolták:<br />

12.2.3. A k-epszilon modell tulajdonságai<br />

Cν = 0,09 (12.18)<br />

C1ε = 1,44 (12.19)<br />

C2ε = 1,92 (12.20)<br />

σk = 1 (12.21)<br />

σε = 1,3 (12.22)<br />

Ebben a fejezetben néhány speciális esetben nézzük meg milyen megoldásokat ad<br />

az egyenletrendszer. Ezen elemzésben vizsgáljuk először a homogén turbulencia<br />

esetét, ekkor az egyenletrendszer a következő alakot ölti:<br />

Csillapodó turbulencia<br />

dtk = P − ε (12.23)<br />

dtε = C1εP ε ε<br />

− C2εε (12.24)<br />

k k<br />

Vázlat verzió<br />

Csillapodó turbulencia esetén P = 0 így az egyenletrendszer könnyen megoldható:<br />

� �−n t<br />

k(t) = k0<br />

(12.25)<br />

Saját használatra<br />

t0<br />

� �−n−1 t<br />

ε(t) = ε0<br />

t0<br />

(12.26)

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

Saved successfully!

Ooh no, something went wrong!