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 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)

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

Saved successfully!

Ooh no, something went wrong!