3.6 Illustrations numériques et exemples de codes de calcul
87
Fig. 3.13: Champ pression calculés par FreeFEM++ (figure de gauche) et
Matlab (figure de droite) à t = 0.4.
schéma WENO en sa version WENO-Z [29, 3] pour la partie advection (ce
schéma est décrit dans l’annexe). Le code est donc basé sur des méthodes
tout à fait différentes de celles du code éléments finis. On obtient des résultats
numériques remarquablement proches, compte tenu du fait que l’interpolation
P 1 de la fonction convect provoque un étalement numérique plus important
de l’interface, visible sur les courbes de pression. On voit dans la version
différences finies que le gradient de pression est capturé essentiellement sur
une maille, ce qui correspond à la valeur de ε choisie dans ce cas. L’étalement
de la force inhérent à la formulation Level Set n’induit donc pas de dissipation
dans l’écoulement au cours du temps.
0
0.1
0.2
0.3
0.4
0.5
0.6
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.1
0.2
0.3
0.4
0.5
0.6
Fig. 3.14: Profil de pression calculés par FreeFEM++ (figure de gauche) et
Matlab (figure de droite) à t = 0.4.