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