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

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

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

Saved successfully!

Ooh no, something went wrong!