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 />

13. FEJEZET. A NAGY ÖRVÉNY SZIMULÁCIÓ 87<br />

ahogy ezt korábban láttuk.<br />

13.2. A nagy örvény szimuláció alapgondolata<br />

A nagy örvény szimuláció a közvetlen numerikus szimuláció és Reynolds átlagolt<br />

modellezés előnyeit próbálja ötvözni, emellett a hátrányokat is sikerül ötvöznie.<br />

A megközelítés alapgondolata, hogy az energia kaszkád koncepció helyes és minden<br />

esetben alkalmazható az a megközelítés, hogy egy bizonyos mérettartomány<br />

alatt a turbulencia közel univerzális. Ezek alapján az javasoljuk, hogy a nagy örvényeket<br />

(Large-Eddy) szimuláljuk és csak a kisebbeket modellezzük, mivel azokat<br />

jóval könnyebb. Így a nehezen modellezhető nagy léptékeket szimuláljuk, tehát az<br />

eredmény pontosabb lesz, kicsiket modellezzük, tehát nem kell a Kolmogorov létékek<br />

felbontására is képes finom hálón számolnunk, így a szimulációnk olcsóbb<br />

lesz. A hátrányok is természetesen megjelennek, az eredmény nem lesz tökéletes,<br />

de jóval lassabb, drágább lesz mint egy RANS. Ezek után nézzük meg, hogyan<br />

vezethető le a megoldandó egyenlet.<br />

13.3. A LES egyenlet levezetése<br />

13.3.1. A szűrés<br />

A nagy örvény szimuláció alapegyenletét hagyományosan a Navier-Stokes egyenlet<br />

térbeli szűrésével állítják elő, mi is ezt az utat követjük itt, lássuk hogyan simítható<br />

térben egy áramlási mennyiség. Definiáljuk az átlagolást a következőképpen:<br />

〈ϕ〉 (xj, t) def<br />

�<br />

= G∆(ri; xj) ϕ(xj − ri, t)dri (13.1)<br />

V<br />

ami egy G∆, egy ∆ tipikus szélességű magfüggvénnyel képzett konvolúciós integrál,<br />

melyben V a teljes vizsgált tartomány, ahol ϕ(xj, t) értelmezve van. Itt G∆-ra<br />

igaz, hogy az első változóban kompakt tartójú (matematikailag kevésbé szabatosan:<br />

a nem nulla értékkészletű tartománya zárt). Ezen felül, hogy egy konstans<br />

függvény szűrtje önmaga legyen, a következőnek is igaznak kell lennie.<br />

�<br />

G∆(ri; xj)dri = 1 (13.2)<br />

Vázlat verzió<br />

V<br />

Ha G∆(ri; xj) homogén (a második változóban) és izotrop (az első változóban)<br />

azaz csak ri abszolút értéke számít a függvényben, azaz G∆(|ri|) egyváltozós<br />

függvény. Pár tipikusan használható ilyen magfüggvény látható a 13.1 ábrán.<br />

Saját használatra

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

Saved successfully!

Ooh no, something went wrong!