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.
100 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 />
4.1.2.2 Valeurs initiales <strong>de</strong>s paramètres <strong>de</strong> l’écoulem<strong>en</strong>t<br />
Les champs spectraux caractérisés par le spectre PP sont <strong>en</strong>suite transportés dans l’espace<br />
physique (cf. 3.3.2). À l’issue <strong>de</strong> <strong>la</strong> 1 ère itération durant <strong>la</strong>quelle les équations <strong>de</strong> Navier-Stokes<br />
sont résolues, on extrait <strong>la</strong> valeur <strong>de</strong>s paramètres initiaux (Tab. 4.4).<br />
Paramètres Référ<strong>en</strong>ce Ensemble A Ensemble B<br />
initiaux S1 A = SB 2 S2 A S3 A S1 B S3<br />
B<br />
Couleur rouge vert bleu cyan orange<br />
u ′ 0 0.163 0.229 0.325 0.109 0.215<br />
k 0 3.99 10 −2 7.87 10 −2 0.159 1.79 10 −2 6.95 10 −2<br />
ε 0 7.76 10 −3 1.53 10 −3 3.09 10 −2 1.59 10 −3 2.37 10 −2<br />
ν 0 2.03 10 −3 2.03 10 −3 2.03 10 −3 2.03 10 −3 2.03 10 −3<br />
Re T 0 100 200 400 100 100<br />
Re λ0 25.8 36.5 51.6 25.8 25.8<br />
∆t 6.66 10 −3 5.85 10 −3 4.67 10 −3 7.97 10 −3 5.85 10 −3<br />
L T 0 1.03 1.44 2.05 1.51 0.77<br />
L ii0 0.40 0.40 0.40 0.60 0.31<br />
λ T 0 0.14 0.14 0.14 0.21 0.11<br />
η 0 3.21 10 −2 2.71 10 −2 2.28 10 −2 4.78 10 −2 2.43 10 −2<br />
κ max 140 140 140 140 140<br />
κ e 6 6 6 4 8<br />
κ d 7.35 7.35 7.35 4.90 9.80<br />
Table 4.4 – Valeurs initiales <strong>de</strong>s paramètres <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> initialisée par un spectre PP<br />
4.1.2.3 Visualisation <strong>de</strong>s champs initiaux <strong>de</strong> l’écoulem<strong>en</strong>t<br />
La <strong>de</strong>scription précé<strong>de</strong>nte <strong>de</strong> l’état initial <strong>de</strong>s écoulem<strong>en</strong>ts (Tab. 4.4) est complétée par une<br />
visualisation <strong>de</strong>s champs <strong>de</strong> vorticité correspondants (Fig. 4.2). Ces représ<strong>en</strong>tations permett<strong>en</strong>t<br />
<strong>de</strong> constater l’effet du nombres <strong>de</strong> Reynolds Re T et du nombre d’on<strong>de</strong> κ e sur l’écoulem<strong>en</strong>t généré.<br />
Ainsi, les spectres S A montre qu’une augm<strong>en</strong>tation <strong>de</strong> Re T provoque une hausse <strong>de</strong> l’amplitu<strong>de</strong><br />
<strong>de</strong> l’agitation turbul<strong>en</strong>te alors que les spectres S B illustr<strong>en</strong>t l’influ<strong>en</strong>ce du choix <strong>de</strong> κ e sur <strong>la</strong><br />
taille <strong>de</strong>s structures porteuses d’énergie.<br />
4.1.3 Initialisation par un spectre <strong>de</strong> Von-Kármán corrigé par Pao (VKP)<br />
Les efforts cons<strong>en</strong>tis pour implém<strong>en</strong>ter un spectre VKP afin d’initialiser <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> sont<br />
motivés par <strong>la</strong> capacité <strong>de</strong> ce spectre à modéliser <strong>de</strong>s zones inertielles plus gran<strong>de</strong>s que dans le<br />
cas d’un spectre PP. Ainsi, il jouit d’une meilleure représ<strong>en</strong>tativité <strong>de</strong>s écoulem<strong>en</strong>ts turbul<strong>en</strong>ts.