Raffinement local adaptatif et méthodes multiniveaux pour la simulation d'écoulements multipĥasiques.

This thesis is devoted to the study of some numerical and mathematical aspects of incompressible multiphase flows simulations with a diffuse interface Cahn-Hillliard/Navier-Stokes model. The space discretisation is performed thanks to the finite elements method. The presence of different scales in the system suggests the use of a local adaptive refinement method. The algorithm, that we introduced, allows to implicitly handle the non conformities of the generated meshes to produce conformal finite elements approximation spaces. Moreover, we show that this method can be exploited to build multigrid preconditioners. Concerning the time discretization, we begin by the study of the Cahn-Hilliard system. A semi-implicit scheme, which ensure the decrease of the discrete energy, is proposed to remedy to convergence failures of the Newton method used to solve this (non linear) system. We show existence and convergence of discrete solutions towards the weak solution of the system. We then continue this study by providing an inconditionnaly stable time discretization of the complete Cahn-Hilliard/Navier-Stokes model. An important point is that this discretization does not strongly couple the Cahn-Hilliard and Navier-Stokes systems allowing to independently solve the two systems in each time step. We show the existence of discrete solutions and, in the case where the three fluids have the same densities, their convergence towards weak solutions. We study, to finish this part, different issues linked to the use of the incremental projection method. Finally, the last part presents several examples of numerical simulations, diphasic and triphasic, in two and three dimensions.Cette thèse est consacrée à l'étude de certains aspects numériques et mathématiques liés à la simulation d'écoulements incompressibles triphasiques à l'aide d'un modèle à interfaces diffuses de type Cahn-Hilliard/Navier-Stokes. La discrétisation spatiale est effectuée par éléments finis. La présence d'échelles très différentes dans le système suggère l'utilisation d'une méthode de raffinement local adaptatif. La procédure mise en place permet de tenir compte implicitement des non conformités des maillages générés, pour produire in fine des espaces d'approximation conformes. Nous montrons, en outre, qu'il est possible d'exploiter cette méthode pour construire des préconditionneurs multigrilles. Concernant la discrétisation en temps, notre étude a commencé par celle du système de Cahn-Hilliard. Pour remédier aux problèmes de convergence de la méthode de Newton utilisée pour résoudre ce système (non linéaire), nous proposons un schéma semi-implicite permettant de garantir la décroissance de l'énergie. Nous montrons l'existence et la convergence des solutions discrètes. Nous poursuivons ensuite cette étude en donnant une discrétisation en temps inconditionnellement stable du modèle complet Cahn-Hilliard/Navier-Stokes ne couplant pas fortement les deux systèmes. Nous montrons l'existence des solutions discrètes et, dans le cas où les trois fluides ont la même densité, nous montrons leur convergence. Nous étudions, pour terminer cette partie, diverses problématiques liées à l'utilisation de la méthode de projection incrémentale. Enfin, la dernière partie présente plusieurs exemples de simulations numériques, diphasiques et triphasiques, en deux et trois dimensions