Une méthode de calcul de couches limites tridimensionnelles se développant autour d'obstacles profiles quelconques a été mise au point. Le domaine d'application de celle-ci va, en laminaire, de l'incompressible au supersonique. En écoulement turbulent, il est limité par la validité du modèle de longueur de mélange utilisé. Le schéma numérique adopté est de type semi-implicite. Les équations de quantité de mouvement et d'énergie sont discrétisées suivant les lignes de courant locales, cela afin de respecter rigoureusement les contraintes liées à la nature parabolique du système d'équations. Contrairement aux méthodes habituelles qui utilisent les coordonnées généralisées, on se sert comme intermédiaire de calcul du plan tangent osculateur en chaque nœud du maillage de surface lié à l'obstacle, évitant ainsi d'imposer une hypothèse de régularité supplémentaire pour les lignes de maillage de données. Il est ainsi possible de limiter, voire de supprimer, la phase de prétraitement des données venant du calcul fluide parfait.