On s'intéresse dans ce travail à l'étude algorithmique parallèle et à l'implémentation sous MPI (Message Passing Interface) d'un solveur numérique simulant un écoulement diphasique eau-air modélisant l'aération mécanique dans un lac eutrophe. Ce modèle est gouverné par les équations de Navier-Stokes bidimentionnelles écrites en formulation fonction courant-tourbillon. La méthode de résolution combine la méthode des caractéristiques pour la discrétisation temporelle et une méthode d'éléments finis C^0 optimal stabilisé pour l'approximation spatiale. Le système matriciel discret obtenu est résolu par la méthode du Bigradient conjugué stabilisé avec préconditionneur. L'analyse de performance en temps de calcul du code séquentiel implémenté montre que 90 Après une phase d'analyse algorithmique parallèle de cette partie identifiant l'interdépendance des différentes tâches et données, une implémentation du code sous MPI a été mise en place sur une machine virtuelle à mémoire distribuée. On présente la version parallélisé du code qui f t testée pour l'aspect performance. Les résultats de performance parallèle, tant sur le plan temps CPU que Elapsed, confirment l'efficacité de l'algorithme implémenté.