10.01.2015 Views

Simulation numérique directe de la turbulence en présence d ... - ISAE

Simulation numérique directe de la turbulence en présence d ... - ISAE

Simulation numérique directe de la turbulence en présence d ... - ISAE

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.

Chapitre 4. <strong>Simu<strong>la</strong>tion</strong> <strong>de</strong> <strong>la</strong> Turbul<strong>en</strong>ce Homogène Isotrope (THI) 121<br />

<strong>la</strong> dissipation par effets visqueux et <strong>de</strong> réinjecter cette dissipation sous <strong>la</strong> forme d’une force<br />

proportionnelle à <strong>la</strong> fluctuation <strong>de</strong> vitesse. Il a montré que ce forçage permet <strong>de</strong> recouvrir une<br />

gamme <strong>de</strong> nombre d’on<strong>de</strong> al<strong>la</strong>nt <strong>de</strong> l’échelle intégrale à <strong>la</strong> zone inertielle. Autrem<strong>en</strong>t dit, le<br />

forçage agit principalem<strong>en</strong>t sur les tourbillons porteurs d’énergie. Son principal avantage est qu’il<br />

s’applique <strong>directe</strong>m<strong>en</strong>t dans le domaine physique avec <strong>de</strong>s conditions aux limites non périodiques.<br />

Le terme <strong>de</strong> forçage s’ajoute aux équations <strong>de</strong> conservation <strong>de</strong> <strong>la</strong> quantité <strong>de</strong> mouvem<strong>en</strong>t<br />

(2.4) avec :<br />

( )<br />

∂ui<br />

ρ<br />

∂t + u ∂u i<br />

j = ρF i − ∂P + µ ∂2 u i<br />

+ ρBu i (4.43)<br />

∂x j ∂x i ∂x j ∂x j<br />

où le paramètre B est donnée par <strong>la</strong> re<strong>la</strong>tion :<br />

B = ɛ + u i.∇p<br />

ρu i .u i<br />

= ɛ + u i.∇p<br />

3ρu ′2 (4.44)<br />

La figure 4.24 représ<strong>en</strong>te <strong>la</strong> variation <strong>de</strong> l’énergie cinétique turbul<strong>en</strong>te moy<strong>en</strong>ne k et du taux<br />

moy<strong>en</strong> <strong>de</strong> dissipation turbul<strong>en</strong>te ε dans le cas sans et avec forçage. Ce forçage linéaire est<br />

décl<strong>en</strong>ché une fois <strong>la</strong> THI établie, dans le cas prés<strong>en</strong>t à t = 5. Il apparaît que <strong>la</strong> valeur <strong>de</strong> k reste<br />

constante pour t > 5 tandis que celle <strong>de</strong> ε ne cesse <strong>de</strong> décroître.<br />

(a) Évolution <strong>de</strong> k<br />

(b) Évolution <strong>de</strong> ε<br />

Figure 4.24 – Conséqu<strong>en</strong>ces du forçage sur les évolution <strong>de</strong> k (a) et ε (b)<br />

Ceci a pour conséqu<strong>en</strong>ce <strong>de</strong> modifier les profils <strong>de</strong> vorticité au cours du forçage comme<br />

l’indique <strong>la</strong> figure 4.25. Nous remarquons sur cette figure que les petites structures tourbillonnaires<br />

sont bi<strong>en</strong> simulées à l’instant initial, mais celles-ci sont dissipées par effets visqueux et<br />

ne sont pas alim<strong>en</strong>tées par les grosses structures. Cette évolution est due au terme <strong>de</strong> forçage<br />

<strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> qui ne <strong>la</strong>isse pas <strong>de</strong> temps aux grosses structures tourbillonnaires <strong>de</strong> recréer <strong>la</strong><br />

casca<strong>de</strong> énergétique du fait <strong>de</strong> <strong>la</strong> ré-injection <strong>de</strong> <strong>la</strong> dissipation turbul<strong>en</strong>te sur <strong>de</strong>s nombres d’on<strong>de</strong><br />

trop petits.<br />

4.4.1.3 Une taille <strong>de</strong> rugosités à atteindre<br />

La principale motivation <strong>de</strong> l’utilisation d’un forçage est l’obt<strong>en</strong>tion <strong>de</strong> profils ab<strong>la</strong>tés possédant<br />

<strong>de</strong>s rugosités <strong>de</strong> taille significative (<strong>de</strong> l’ordre du µm). En l’occurr<strong>en</strong>ce, nous v<strong>en</strong>ons <strong>de</strong><br />

montrer les difficultés r<strong>en</strong>contrées avec l’utilisation d’un forçage dans le domaine physique. En

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

Saved successfully!

Ooh no, something went wrong!