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 2. Les équations du problème physique 53<br />
• Approximation <strong>de</strong> <strong>la</strong> conductivité thermique λ<br />
Afin d’estimer <strong>la</strong> conductivité thermique d’un gaz supposé parfait, <strong>la</strong> première approximation<br />
utilisée est l’hypothèse du nombre <strong>de</strong> Prandtl constant avec :<br />
P r = µC p<br />
λ<br />
En ce qui concerne les gaz multi-espèces, nous appliquerons, là <strong>en</strong>core, <strong>la</strong> loi <strong>de</strong> Wilke avec :<br />
λ =<br />
n e ∑<br />
α=1<br />
(2.42)<br />
X α λ α<br />
∑ ne<br />
β=1 X (2.43)<br />
βψ λ,αβ<br />
où λ α est <strong>la</strong> conductivité thermique du gaz α et ψ λ,αβ correspond à l’interaction <strong>de</strong> conductivité<br />
thermique <strong>en</strong>tre les espèces α et β par :<br />
ψ λ,αβ =<br />
( ( ) 1 ( ) 1<br />
) 2<br />
1 + λα 2 M β 4<br />
λβ M α<br />
[ ( )] 1<br />
8 1 + Mα 2<br />
M β<br />
(2.44)<br />
Maint<strong>en</strong>ant que nous avons rappelé l’<strong>en</strong>semble <strong>de</strong>s équations utilisées pour simuler les écoulem<strong>en</strong>ts<br />
turbul<strong>en</strong>ts composés d’espèces réactives, nous allons prés<strong>en</strong>ter les outils, mis <strong>en</strong> p<strong>la</strong>ce<br />
durant cette thèse, nécessaires à <strong>la</strong> validation <strong>de</strong> <strong>la</strong> simu<strong>la</strong>tion <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong>.<br />
2.2 L’approche statistique <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong><br />
La mise <strong>en</strong> p<strong>la</strong>ce d’une théorie statistique <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> s’explique par le caractère aléatoire<br />
<strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong> (cf. 1.2.1.1). Le fait <strong>de</strong> considérer un <strong>en</strong>semble <strong>de</strong> réalisations, plutôt que chaque<br />
réalisation prise individuellem<strong>en</strong>t, va nous permettre d’i<strong>de</strong>ntifier les propriétés et les phénomènes<br />
reproductibles <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong>. Nous définissons ici les outils statistiques adéquats <strong>de</strong>stinés à<br />
vali<strong>de</strong>r le phénomène <strong>de</strong> <strong>turbul<strong>en</strong>ce</strong> modélisé par les équations du paragraphe 2.1.<br />
2.2.1 Motivations<br />
Que ce soit à l’époque <strong>de</strong>s premières observations <strong>de</strong> Léonard <strong>de</strong> Vinci ou <strong>de</strong>s t<strong>en</strong>tatives <strong>de</strong><br />
caractérisation d’Osborne Reynolds, il n’était pas possible d’accé<strong>de</strong>r, par le calcul, à tous les<br />
détails <strong>de</strong>s fluctuations <strong>de</strong> l’agitation turbul<strong>en</strong>te. La situation a changé au début <strong>de</strong>s années<br />
1970 avec les premiers travaux d’Orszag [41] qui a introduit <strong>la</strong> résolution <strong>numérique</strong> <strong>directe</strong><br />
<strong>de</strong>s équations <strong>de</strong> Navier-Stokes. Ces calculs, ambitieux, nécessit<strong>en</strong>t <strong>en</strong>core <strong>de</strong>s temps <strong>de</strong> calcul<br />
très longs et pour <strong>la</strong> plupart <strong>de</strong>s problèmes intéressant l’ingénieur, <strong>la</strong> résolution d’équations<br />
moy<strong>en</strong>nées suffit. Dans les <strong>de</strong>ux cas, l’approche statistique est considérée car <strong>la</strong> détermination <strong>de</strong><br />
données statistiques <strong>en</strong> <strong>turbul<strong>en</strong>ce</strong> peut s’<strong>en</strong>visager selon <strong>de</strong>ux types <strong>de</strong> démarches : <strong>la</strong> statistique<br />
avant et après résolution.<br />
2.2.1.1 Statistique avant résolution<br />
Il s’agit <strong>de</strong> l’approche probabiliste standard où <strong>la</strong> vitesse, <strong>la</strong> pression, <strong>la</strong> température, <strong>en</strong>tre<br />
autres, sont considérées comme <strong>de</strong>s fonctions aléatoires, décomposables <strong>en</strong> une somme d’une<br />
partie moy<strong>en</strong>ne et d’une partie fluctuante (cf. 2.2.2.3). Sur le p<strong>la</strong>n théorique, l’évolution d’une<br />
mesure <strong>de</strong> probabilité <strong>de</strong>s équations du mouvem<strong>en</strong>t n’a pas été résolue dans le cas général.