Simulation des écoulements turbulents compressibles par une ...

Deuxième partie 41
U=Uxex+Urer (2.50)
On suppose de plus qu'elle ne comporte pas de composante tangentielle, les effets de giration ne sont
donc pas pris en compte. Les dérivées )4)/ar (composante radiale du gradient de 4)) sont décomposées
sous la forme:
òcbiò
(2.51)
Les contributions sont intégrés dans les termes de divergence, tandis que les quotients de la
forme 4)/r engendreront des termes sources supplémentaires. Les équations se mettent alors sous la
forme:
pu, + 1rpu =
ap a 2 ___ ia i a ia
a rarx
apur a
at
3Re r
ÒE a ia 'yae
=-[-(uxaxx+uraxr+)
(u arx+uraff+
Le vecteur des variables d'état s'exprime en fonction des composantes axiale et radiale de la vitesse:
p
pu,
pu
E
r
s
(2.52)
(2.53)
(2.54)
(2.55)
(2.56)