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.
Chapitre 4. <strong>Simu<strong>la</strong>tion</strong> <strong>de</strong> <strong>la</strong> Turbul<strong>en</strong>ce Homogène Isotrope (THI) 103<br />
Figure 4.3 – Visualisation du champ initial <strong>de</strong> vorticité (spectre VKP)<br />
4.1.3.3 Limites quant à l’utilisation du spectre VKP<br />
La définition d’un spectre VKP est plus contraignante que celle d’un spectre PP du fait <strong>de</strong><br />
<strong>la</strong> possibilité <strong>de</strong> paramétrer <strong>la</strong> taille <strong>de</strong> <strong>la</strong> zone inertielle. En effet pour un spectre PP, <strong>la</strong> taille<br />
<strong>de</strong> cette zone est fixée à :<br />
(√ )<br />
3<br />
2 − 1<br />
κ d − κ e =<br />
κ e (4.11)<br />
La situation est différ<strong>en</strong>te pour un spectre VKP pour lequel les propriétés turbul<strong>en</strong>tes sont<br />
définies à partir <strong>de</strong> <strong>la</strong> taille <strong>de</strong> cette zone. Cette procédure d’initialisation ne permet pas <strong>de</strong><br />
choisir, indép<strong>en</strong>damm<strong>en</strong>t <strong>de</strong>s autres paramètres, le montant énergétique <strong>de</strong> l’écoulem<strong>en</strong>t. En<br />
effet comme le montre <strong>la</strong> figure 4.4, le nombre <strong>de</strong> Reynolds turbul<strong>en</strong>t Re T admet, avant <strong>la</strong><br />
première itération, un minimum d’<strong>en</strong>viron 400 atteint pour K d /K e ≃ 11.6. Dans le paragraphe<br />
précé<strong>de</strong>nt, nous avons révélé l’exist<strong>en</strong>ce d’un écart, propre à l’usage du spectre VKP, <strong>en</strong>tre les<br />
valeurs initiales et celles calculées après une itération. Vu que nous conservons <strong>la</strong> <strong>de</strong>uxième<br />
catégorie <strong>de</strong> valeurs pour définir l’état initial <strong>de</strong> notre calcul, le minimum pour Re T ret<strong>en</strong>u est<br />
donc ram<strong>en</strong>é à 315. Cette valeur reste s<strong>en</strong>siblem<strong>en</strong>t élevée et le spectre VKP ne pourra donc pas<br />
être utilisé pour définir <strong>de</strong>s écoulem<strong>en</strong>ts à faible nombre <strong>de</strong> Reynolds turbul<strong>en</strong>t.<br />
Figure 4.4 – Nombre <strong>de</strong> Reynolds turbul<strong>en</strong>t <strong>en</strong> fonction <strong>de</strong> K d /K e<br />
Il est difficile <strong>de</strong> donner une signification physique au rapport K d /K e qui ne peut être interprété<br />
comme un indicateur <strong>de</strong> <strong>la</strong> taille <strong>de</strong> <strong>la</strong> zone inertielle. C’est pourquoi nous représ<strong>en</strong>tons