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 59<br />
P ij = −<br />
(<br />
)<br />
∂U j ∂U i<br />
u i u k + u j u k<br />
∂x k ∂x k<br />
est le terme <strong>de</strong> production par le mouvem<strong>en</strong>t moy<strong>en</strong> ;<br />
Dij u = − ∂u iu j u k<br />
est <strong>la</strong> corré<strong>la</strong>tion triple <strong>de</strong>s vitesse (diffusion) ;<br />
∂x k<br />
D p ij = − 1 (<br />
∂pui<br />
+ ∂pu )<br />
j<br />
est <strong>la</strong> diffusion par fluctuation <strong>de</strong> pression ;<br />
ρ ∂x j ∂x i<br />
Dij ν = ν ∂2 u i u j<br />
∂x k ∂x k<br />
est <strong>la</strong> diffusion molécu<strong>la</strong>ire ;<br />
Π ij = p (<br />
∂ui<br />
+ ∂u )<br />
j<br />
ρ ∂x j ∂x i<br />
est le terme <strong>de</strong> corré<strong>la</strong>tion pression-déformation ;<br />
ε ij =<br />
( )<br />
∂ui ∂u j<br />
2ν<br />
∂x j ∂x i<br />
est <strong>la</strong> pseudo-dissipation.<br />
Ces équations mett<strong>en</strong>t <strong>en</strong> jeu <strong>de</strong>s termes <strong>de</strong> transport (advection et diffusion) qui disparaiss<strong>en</strong>t<br />
dans le cadre d’une THI. En revanche, dans le cadre d’écoulem<strong>en</strong>ts confinés, ces termes<br />
<strong>de</strong>vi<strong>en</strong>n<strong>en</strong>t un <strong>en</strong>semble <strong>de</strong> sources et <strong>de</strong> puits locaux qui provoqu<strong>en</strong>t, soit une variation <strong>de</strong><br />
l’énergie cinétique turbul<strong>en</strong>te, soi<strong>en</strong>t <strong>de</strong>s transferts intercomposantes. Lors <strong>de</strong>s simu<strong>la</strong>tions d’écoulem<strong>en</strong>ts<br />
<strong>de</strong> proche paroi, il s’agira alors <strong>de</strong> modéliser les termes <strong>de</strong> production, <strong>de</strong> corré<strong>la</strong>tion<br />
pression-déformation et <strong>de</strong> dissipation afin <strong>de</strong> vali<strong>de</strong>r <strong>la</strong> simu<strong>la</strong>tion <strong>de</strong> <strong>la</strong> <strong>turbul<strong>en</strong>ce</strong>.<br />
En abandonnant l’hypothèse <strong>de</strong> viscosité <strong>de</strong> <strong>turbul<strong>en</strong>ce</strong>, les modèles RSTE permett<strong>en</strong>t <strong>de</strong><br />
s’affranchir d’une re<strong>la</strong>tion locale <strong>en</strong>tre les t<strong>en</strong>sions <strong>de</strong> Reynolds et l’écoulem<strong>en</strong>t moy<strong>en</strong>. Ils permett<strong>en</strong>t<br />
donc <strong>de</strong> mieux pr<strong>en</strong>dre <strong>en</strong> compte les effets d’histoire ou <strong>en</strong>core d’anisotropie <strong>de</strong> <strong>la</strong><br />
<strong>turbul<strong>en</strong>ce</strong>. Des modèles d’ordre supérieur, basés sur <strong>la</strong> fermeture d’équations <strong>de</strong> transport d’ordre<br />
égal ou supérieur à trois, sont difficilem<strong>en</strong>t <strong>en</strong>visageables, d’une part par le manque <strong>de</strong><br />
données expérim<strong>en</strong>tales sur ce type <strong>de</strong> corré<strong>la</strong>tions d’ordre élevé, et d’autre part par <strong>la</strong> lour<strong>de</strong>ur<br />
<strong>de</strong> <strong>la</strong> démarche <strong>en</strong>gagée. Les modèles RSTE constitu<strong>en</strong>t donc un bon intermédiaire <strong>en</strong>tre les<br />
modèles à viscosité <strong>de</strong> <strong>turbul<strong>en</strong>ce</strong> et les métho<strong>de</strong>s <strong>de</strong> simu<strong>la</strong>tion <strong>numérique</strong> <strong>directe</strong>.<br />
2.2.4 Analyse corré<strong>la</strong>toire<br />
C’est Taylor [60] qui étudia le premier les corré<strong>la</strong>tions statistiques <strong>de</strong>s fluctuations <strong>de</strong> vitesse<br />
<strong>en</strong> <strong>de</strong>ux points distincts d’un champ turbul<strong>en</strong>t. Leur emploi est un prolongem<strong>en</strong>t naturel <strong>de</strong><br />
l’analyse corré<strong>la</strong>toire <strong>en</strong> un point, dont on a vu l’importance, notamm<strong>en</strong>t avec le t<strong>en</strong>seur <strong>de</strong>s<br />
contraintes <strong>de</strong> Reynolds −ρu i u j .<br />
Bi<strong>en</strong> que l’utilisation <strong>de</strong> ces corré<strong>la</strong>tions doubles ne soit pas aisée dans le cas général, leur<br />
étu<strong>de</strong> <strong>en</strong> situation <strong>de</strong> <strong>turbul<strong>en</strong>ce</strong> homogène et isotrope pr<strong>en</strong>d tout son s<strong>en</strong>s. En effet, <strong>de</strong>s propriétés<br />
remarquables sont alors i<strong>de</strong>ntifiables et nous permett<strong>en</strong>t <strong>de</strong> nous doter d’élém<strong>en</strong>ts <strong>de</strong><br />
validation supplém<strong>en</strong>taires.<br />
2.2.4.1 Définition <strong>de</strong>s outils mathématiques<br />
Soi<strong>en</strong>t <strong>de</strong>ux points P 1 ( ⃗x 1 , t) et P 2 ( ⃗x 2 , t) du domaine spatial muni d’un repère orthonormé<br />
(⃗i,⃗j, ⃗ k). La corré<strong>la</strong>tion double spatio-temporelle <strong>de</strong> vitesse, s’exprime par :<br />
Q i,j (P 1 , P 2 , t 1 , t 2 ) = u i (P 1 , t 1 )u j (P 2 , t 2 ) (2.68)