Auteur / Autrice : | Martin Vohralík |
Direction : | Danielle Hilhorst, Jiř̌í Maryš̌ka |
Type : | Thèse de doctorat |
Discipline(s) : | Mathématiques |
Date : | Soutenance en 2004 |
Etablissement(s) : | Paris 11 en cotutelle avec Université technique tchèque, Prague |
Partenaire(s) de recherche : | autre partenaire : Université de Paris-Sud. Faculté des sciences d'Orsay (Essonne) - Université technique tchèque. Faculté de Génie nucléaire et de Sciences de l'ingénieur (Prague) |
Résumé
Nous étudions des méthodes numériques pour la simulation de l'écoulement et du transport de contaminants en milieux poreux et fracturés. Au chapitre 1, nous proposons un schéma permettant une discrétisation efficace, robuste, conservative et stable des équations de convection–réaction–diffusion paraboliques dégénérées sur des maillages non structurés en dimensions 2 ou 3 d'espace. Nous discrétisons le terme de diffusion, en général anisotrope, par la méthode des éléments finis non conformes et les autres termes par celle des volumes finis et démontrons l'existence et l'unicité d'une solution discrète et sa convergence vers une solution faible. Nous proposons finalement une variante de ce schéma pour des maillages qui ne se raccordent pas et nous l'appliquons aux simulations réelles. Au chapitre 2, nous présentons une démonstration constructive des inégalités de Poincaré–Friedrichs discrètes et indiquons les valeurs optimales des constantes dans ces inégalités. Ces résultats sont importants dans l'analyse de méthodes numériques non conformes. Au chapitre 3, nous montrons que la méthode des éléments finis mixtes de Raviart–Thomas de plus bas degré est équivalente à un schéma de volumes finis à plusieurs points. Cette approche permet de réduire le temps de calcul de la méthode des éléments finis mixtes, tout en conservant sa grande précision, ce qui est confirmé par les tests numériques. Enfin, au chapitre 4, nous proposons une version de la méthode des éléments finis mixtes pour la simulation de l'écoulement dans un réseau de fractures perturbant un massif rocheux, démontrons qu'elle est bien posée et étudions sa relation avec la méthode des éléments finis non conformes.