Dans cette thèse, nous combinons des méthodes itératives basées sur l'élément fini P₂ hiérarchique et une stratégie d'adaptation de maillage anisotrope fondée sur l'intersection de métriques pour résoudre les écoulements turbulents avec ou sans surface libre. Les métriques sont obtenues en considérant la matrice hessienne des différentes composantes du gradient de chaque variable. Le modèle k-e , en formulation logarithmique, est utilisé pour modéliser la turbulence, et la surface libre est calculée par la méthode des surfaces de niveau («level sets»). Après avoir résolu quelques problèmes classiques qui mettent en évidence l'efficacité du solveur itératif et de la stratégie d'adaptation de maillage anisotrope, nous résolvons un problème d'écoulement tourbillonnaire entrant dans un diffuseur conique. Cette simulation a beaucoup d'intérêt pour les écoulements dans un aspirateur hydraulique. L'ensemble de notre approche est ensuite appliquée pour simuler un écoulement turbulent autour d'un cylindre perçant la surface libre. Finalement, les résultats positifs obtenus dans cette étude nous permettent d'aborder le phénomène d'aquaplanage en simulant l'écoulement turbulent à surface libre autour d'un pneu.