Résumé (FR)

Cette thèse porte sur la modélisation numérique pour l'écoulement et le transport dans les milieux poreux. Nous présentons une nouvelle méthode de couplage entre les réations chimiques et le transport en utilisant une méthode de Newton-Krylov, et nous étudions un modèle d'écoulement qui traite l'intersection des fractures par une méthode de décomposition de domaine.Ce travail est divisé en trois parties: la première partie est une analyse de différents schémas numériques pour la discrétisation du problème d'advection-diffusion, notament par une technique de séparation d'opérateurs, ainsi que leur mise en oeuvre informatique dans un code industriel.La deuxième partie, qui est la contribution majeure de cette thèse, est consacrée à la modélisation et l'implémentation d'une méthode de couplage globale pour le transport réactif. Le système couplé transport-chimie est décrit, après discrétisation en temps, par un système d'équations non linéaires. La taille du système sous-jacent, à savoir le nombre de points de grille multiplié par le nombre d'espèces chimiques, interdit la résolution du système linéaire par une méthode directe. Pour remédier à cette difficulté, nous utilisons une méthode de Newton-Krylov qui évite de former et de factoriser la matrice Jacobienne.Dans la dernière partie, nous présentons un modèle d'écoulement dans un milieu fracturé tridimentionnel, basée sur une méhode de décomposition de domaine, et qui traite l'intersection des fracures. Nous démontrons l'existence et l'unicité de la solution, et nous validons le modèle par des tests numériques.