Simulation des écoulements turbulents compressibles par une ...

Deuxième partie 44
fq r dr dx 22Y + + 7Yk f N. r dr dx (2.63)
ac1r'T '8Yi9Yj9Yk T
où j et k représentent les deux autres noeuds du triangle. Pour que la formulation reste cohérente, les
termes visqueux sont multipliés par ce coefficient correcteur et la formulation faible mixte volumes
fïnis/éléments finis s'écrit:
f(+divF(w) + si(W))p1rdrdx
= y 22Yi7Yj7Yk
l8y. + 9YJ + 9Yk ¡
T
T,ieT '
(2.64)
(divK(W,gradW) + SradW)) N1r dr dx i=1,N
2.2.3 Calcul des intégrales volumiques et des termes sources
Dans l'intégration du terme de dérivation en temps, les intégrales de surface sur la cellules sont
remplaçées par des intégrales sur les éléments toriques. A 2ir près, cela revient à calculer les
intégrales:
r dr dx
(2.65)
Soit K(i) l'ensemble des noeuds voisins de i , on transforme l'intégrale de surface en une intégrale
sur le contour à l'aide de la formule de Stokes:
Ici rdrdx
):df &
jEK(i) ac1
(2.66)