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
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
126 Chapitre 4. <strong>Simu<strong>la</strong>tion</strong> <strong>de</strong> <strong>la</strong> Turbul<strong>en</strong>ce Homogène Isotrope (THI)<br />
(a) Échelles <strong>de</strong> longueurs<br />
(b) Nombres <strong>de</strong> Reynolds turbul<strong>en</strong>ts<br />
Figure 4.29 – Comparaison <strong>de</strong>s forçages linéaire et spectral variable S f 4<br />
atteinte à l’instant t et celle visée k ⋆ . Cette courbe révèle d’abord que le nombre d’on<strong>de</strong> <strong>de</strong><br />
(a) Évolution <strong>de</strong> κ f<br />
(b) Évolution du rapport k/k ⋆<br />
Figure 4.30 – Évolution <strong>de</strong> κ f et du rapport k/k ⋆ pour <strong>de</strong>s forçages fixe et variable<br />
forçage théorique <strong>en</strong> THI est bi<strong>en</strong> décroissant au cours du calcul, ceci est <strong>en</strong> adéquation avec les<br />
remarques du paragraphe 4.4.1. Les autres nombres d’on<strong>de</strong> possè<strong>de</strong>nt bi<strong>en</strong> une valeur fixée tout<br />
au long <strong>de</strong>s simu<strong>la</strong>tions. La figure 4.30 permet aussi <strong>de</strong> constater que le montant énergétique<br />
k ⋆ est efficacem<strong>en</strong>t maint<strong>en</strong>u. En d’autres termes, il n’y a pas <strong>de</strong> phase transitoire du forçage<br />
spectral quelles que soi<strong>en</strong>t les configurations <strong>en</strong>visagées. Cet avantage est dû à <strong>la</strong> distribution<br />
aléatoire <strong>de</strong>s structures turbul<strong>en</strong>tes à chaque génération du champ <strong>de</strong> forçage Ψ 2 .<br />
L’étu<strong>de</strong> <strong>de</strong> <strong>la</strong> figure 4.26 illustrant l’évolution <strong>de</strong> k et ε pour les cas <strong>de</strong> forçages spectraux<br />
étudiés permet <strong>de</strong> juger l’efficacité <strong>de</strong> chacune <strong>de</strong>s configurations quant à <strong>la</strong> manière d’<strong>en</strong>tret<strong>en</strong>ir<br />
<strong>la</strong> <strong>turbul<strong>en</strong>ce</strong>. Ainsi, comme dans le cas du forçage linéaire, on constate que le taux <strong>de</strong> dissipation<br />
turbul<strong>en</strong>te ε est décroissant pour le cas S f 4 suggérant que les échelles <strong>de</strong> longueurs ne sont pas<br />
constantes. Cette propriété est vérifiée sur <strong>la</strong> figure 4.31 qui illustre l’augm<strong>en</strong>tation <strong>de</strong>s échelles<br />
η, λ T et L ii . Ce comportem<strong>en</strong>t est bi<strong>en</strong> <strong>en</strong> adéquation avec l’évolution théorique <strong>de</strong>s ces échelles<br />
dans le cas d’une THI <strong>en</strong> décroissance. Quant aux cas S f 1 , Sf 2 et Sf 3 , ils permett<strong>en</strong>t <strong>de</strong> fixer<br />
correctem<strong>en</strong>t <strong>la</strong> valeur <strong>de</strong>s échelles. La comparaison <strong>de</strong>s champs <strong>de</strong> vorticité obt<strong>en</strong>us pour un