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) 97<br />

4.1 Initialisation du champ turbul<strong>en</strong>t<br />

Pour initialiser <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong>, nous utiliserons <strong>de</strong>ux types <strong>de</strong> spectres : le spectre <strong>de</strong> Passot-<br />

Pouquet (PP) et le spectre <strong>de</strong> Von-Kármán modifié par Pao (VKP). Avant <strong>de</strong> définir les propriétés<br />

<strong>de</strong> chacun d’eux, nous prés<strong>en</strong>tons les caractéristiques générales <strong>de</strong>s simu<strong>la</strong>tions <strong>de</strong> <strong>la</strong><br />

THI.<br />

4.1.1 Configuration d’écoulem<strong>en</strong>t étudiée<br />

4.1.1.1 État <strong>de</strong> référ<strong>en</strong>ce<br />

L’adim<strong>en</strong>sionnem<strong>en</strong>t <strong>de</strong>s équations simulées nécessite <strong>la</strong> définition <strong>de</strong>s états <strong>de</strong> référ<strong>en</strong>ces<br />

(2.24). Considérant un flui<strong>de</strong> à une température T ∞ <strong>de</strong> 400 K et une pression atmosphérique<br />

P atm égale à 1.01325 × 10 5 P a, nous obt<strong>en</strong>ons, selon <strong>la</strong> nature du spectre utilisée, les états <strong>de</strong><br />

référ<strong>en</strong>ce suivants lors <strong>de</strong> l’initialisation :<br />

Variables Spectre PP Spectre VKP Unité<br />

Re ac 500 1000 -<br />

T ref 160 160 K<br />

a ref 43.16 43.16 m.s −1<br />

t ref 1.91×10 −6 3.82×10 −6 s<br />

C p,ref 903.3 903.3 J.kg −1 .K −1<br />

P ref 1.412×10 5 1.412×10 5 P a<br />

ρ ref 2.405 2.405 kg.m −3<br />

µ ref 1.711×10 −5 1.711×10 −5 P a.s (ou P l)<br />

ν ref 7.113×10 −6 7.113×10 −6 m 2 .s −1<br />

L ref 8.241×10 −5 1.648×10 −5 m<br />

Table 4.1 – Paramètres <strong>de</strong> référ<strong>en</strong>ce associés aux spectres initiaux<br />

Dans <strong>la</strong> suite, pour redonner une « réalité dim<strong>en</strong>sionnelle » aux résultats prés<strong>en</strong>tés, il suffira<br />

<strong>de</strong> multiplier les valeurs calculées par les paramètres <strong>de</strong> référ<strong>en</strong>ce du tableau 4.1 afin d’obt<strong>en</strong>ir<br />

les valeurs physiques réelles associées.<br />

4.1.1.2 Simplifications <strong>de</strong>s écoulem<strong>en</strong>ts simulés<br />

Le flui<strong>de</strong> considéré pour cette étu<strong>de</strong> est un flui<strong>de</strong> satisfaisant <strong>la</strong> loi <strong>de</strong>s gaz parfaits. Nous<br />

admettons que ce flui<strong>de</strong> est assimi<strong>la</strong>ble au dioxygène dont <strong>la</strong> viscosité et <strong>la</strong> <strong>de</strong>nsité volumique<br />

sont supposées constantes. On se p<strong>la</strong>ce donc dans le cas d’un écoulem<strong>en</strong>t incompressible, monoespèce<br />

afin <strong>de</strong> s’affranchir <strong>de</strong>s effets <strong>de</strong> <strong>la</strong> cinétique chimique. Les nombres <strong>de</strong> Lewis, Prandtl et<br />

Schmidt sont égaux à :<br />

Le = 1 (4.1)<br />

P r = 0.7 (4.2)<br />

Sc = 1 (4.3)

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

Saved successfully!

Ooh no, something went wrong!