Analyse et estimation d’erreur en volumes finis. Application aux écoulements en milieu poreux et à l’adaptation de maillage
Auteur / Autrice : | Amadou Mahamane |
Direction : | Fayssal Benkhaldoun |
Type : | Thèse de doctorat |
Discipline(s) : | Mathématiques appliquées |
Date : | Soutenance en 2009 |
Etablissement(s) : | Paris 13 |
Mots clés
Mots clés contrôlés
Résumé
Dans cette thèse nous nous intéressons dans une première partie à la simulation numérique d'un écoulement diphasique en milieu poreux par une méthode de volumes finis adaptatifs. La modélisation proposée par G. CHAVENT, conduisant au couplage d'une famille d'équations elliptiques en pression, et d'une équation de convection-diffusion non linéaire en saturation a été retenue. Nous discrétisons la convection par décentrement et la diffusion par le schéma VFdiamant. Nous montrons la L puissance infinie - stabilité du schéma de discrétisation obtenu dans des cas simples. Nous implémentons ce schéma à l'aide des cas tests académiques sur des maillages non structurés adaptatifs; les résultats obtenus reproduisent ceux de la littérature. La seconde partie de ce travail est consacré à l'étude de certains schémas volumes finis dédiés à l'approximation des opérateurs de diffusion. En particulier, nous considérons les schémas : VFdiamant, DDFV développé par P. Omnes et K. Domelevo, VFmixte dû à J. Droniou et à R. Eymard et CVFE développé par B. Amaziane et M. Afif. Ainsi, l'analyse du schéma VFmixte dans le cadre de l'approximation d'une équation de convection-diffusion-réaction a montré la convergence forte de la solution numériques dans L puissance Q (oméga) pour tout q < 2d /(d-2) et la convergence faible du gradient discret dans L puissance carré (oméga)puissance d ; oméga étant un ouvert de IRexposant d, d=2,3. Une analyse d'erreur a posteriori a également été menée, aussi bien pour le schéma DDFV que pour le schéma VFmixte, dans le cas d'une équation de diffusion. L'implémentation des indicateurs d'erreur pour DDFV a montré leur pertinence en termes de localisation de l'erreur. Nous menons enfin une étude numérique comparative des schémas CVFE, DDFV et VFdiamant dans le cadre de l'approximation de l'équation de la chaleur.