Cette thèse porte sur le développement, l'analyse et la mise en oeuvre de nouvelles méthodes bi-grilles en éléments finis pour des équations de réaction-diffusion de type champs de phase (Allen-Cahn, Cahn-Hilliard) ainsi que leur couplage avec les équations de Navier-Stokes 2D incompressible. La présence d'un petit paramètre (largeur de l'interface) dans les modèles à interface diffuse demande à utiliser des schémas temporels implicites et coûteux tandis que ceux semi implicites sont rapides mais limités en stabilité. Les nouveaux schémas introduits ici reposent sur l'utilisation conjointe de deux espaces d'éléments finis, un grossier VH et un fin Vh, de plus grande dimension, permettant de décomposer la solution en partie principale portant les composantes bas modeset en une partie fluctuante portant les modes élevés. L'approche bi-grilles proposée consiste à appliquer les schémas stables (coûteux) sur VH (prédiction) et à effectuer une correction sur Vh à l'aide d'un schéma linéaire dont les composantes modes élevés sont stabilisées. Un gain important en temps CPU est obtenu au prix d'une faible perte de consistance. Dans ce contexte, de nouvelles méthodes numériques sont proposées pour les modèles de champs de phase et leur couplage fluide. Nous donnons des résultats de stabilité et validons l'approche sur des bancs d'essais / This thesis deals with the development, the analysis and the implementation of new bi-grid schemes in finite elements, when applied to phase-field models such as Allen-Cahn (AC) and Cahn-Hilliard (CH) equations but also their coupling with 2D incompressible Navier-Stokes equations. Due to the presence of a small parameter, namely the length of the diffuse interface, and in order to recover the intrinsic properties of the solution, (costly) implicit time schemes must be used; semi-implicit time schemes are fast but suffer from a hard time step limitation. The new schemes introduced in the present work are based on the use of two FEM spaces, one coarse VH and one fine Vh, of larger dimension. This allows to decompose the solution into a main part (containing only low mode components) and a fluctuant part capturing the high mode ones. The bi-grid approach consists then in applying as a prediction an unconditional stable scheme (costly) to VH and to update the solution in Vh by using a high mode stabilized linear scheme. A gain in CPU time is obtained while the consistency is not deteriorated. This approach is extended to NSE and to coupled models (AC/NSE) and (CH/NSE). Stability results are given, the numerical simulations are validated on reference benchmarks