24.02.2013 Views

Simulation numérique du contrôle actif par jets pulsés - Bibliothèque ...

Simulation numérique du contrôle actif par jets pulsés - Bibliothèque ...

Simulation numérique du contrôle actif par jets pulsés - Bibliothèque ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

ANNEXE A. PRÉSENTATION DU CODE DE CALCUL<br />

1 2<br />

-n<br />

+ + 1ou.Vk' - V.((1 + - .E_kn + pn<br />

-n n<br />

CE2 + CElS)E + un.v?+l P n<br />

- V.((i + )VE') -E +CeiPn<br />

Couplage avec les équations de Navier-Stokes<br />

Les équations discrétisées <strong>du</strong> champ moyen et les équations de turbulence sont intégrées<br />

de façon sé<strong>par</strong>ée. Au pas de temps courant, on résoud les équations de NAvIER-STOKES<br />

en utilisant les données de la turbulence évaluées au pas de temps précédent tandis que<br />

les équations de la turbulence sont résolues, <strong>par</strong> le même algorithme implicite itératif que<br />

les équations de NAVIER-STOKES, en utilisant les variables de l'écoulement calculé au pas<br />

de temps précédent.<br />

Modélisation de proche <strong>par</strong>oi<br />

Les équations de transport ci-dessus sont obtenues sous l'hypothèse d'un grand nombre<br />

de Reynolds. Par conséquent, dans une région proche de la <strong>par</strong>oi telle que la sous-couche<br />

visqueuse, où les effets moléculaires deviennent importants, ces équations ne sont plus<br />

valides. Deux approches ont été développées : l'une est basée sur une loi de <strong>par</strong>oi et la<br />

deuxième utilise un modèle bicouche.<br />

Afin d'éviter d'intégrer les équations de turbulence dans ces régions, nous pouvons<br />

utiliser une loi de <strong>par</strong>oi [157]. Ainsi, les équations pour le champ moyen et pour les variables<br />

turbulentes sont intégrées jusqu'à une distance y = Yw Pour y E [O; y], l'écoulement est<br />

supposé suivre la loi de <strong>par</strong>oi<br />

Ü = f (ôurY<br />

La fonction f sé<strong>par</strong>e la couche limite en trois <strong>par</strong>ties. A chaque pas de temps, une<br />

estimation de la vitesse 'û y = Yw peut être obtenue, puis la valeur de la vitesse de<br />

frottement Wr est calculée grâce à une loi de <strong>par</strong>oi (en utilisant une méthode de Newton)<br />

et cette valeur est utilisée pour obtenir des conditions de Dirchlet à la <strong>par</strong>oi pour k et E.<br />

En pratique, le point y = 11w est suffisament proche de la <strong>par</strong>oi pour être supposé faire<br />

<strong>par</strong>tie de celle-ci et les conditions limites peuvent donc être imposées sur la <strong>par</strong>oi. Pour<br />

les équations <strong>du</strong> champ moyen, cela correspond à des conditions de glissement.<br />

222<br />

Lt

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

Saved successfully!

Ooh no, something went wrong!