Résumé : On présente une méthode d'éléments finis mixte pour la résolution des équations bidimensionnelles de Quasi-Stokes. Cette méthode doit servir de solveur de base pour la résolution numérique des équations de Navier-Stoke- s instationnaires.Ces dernières sont linéarisées en utilisant une méthode des caractéristiques du premier ordre pour traiter le terme de dérivée particulaire.La formulation mixte de Quasi-Stokes a pour inconnues la vorticité et la fonction courant discrétisées par des fonctions continues et affines localement. Une approche directe entraînant la perte d'un ordre d'erreur dans ce type de méthode,on utilise une technique de régularisation-st- abilisation par adjonction d'une nouvelle forme discrète.Cette nouvelle formulation reste consistante et est inconditionnellement convergente.En cas de régularité suffisante, une estimation d'erreur optimale en $O(h)$ est obtenue.Des essais numériques sont présentés permettant de confirmer les résultats annoncés dans l'étude théorique.