Généralisation de la méthode Nitsche XFEM pour la discrétisation de problèmess d'interface elliptiques

NXFEM generalization for elliptic interface problems discretization

• Auteur / Autrice : Nelly Barrau
• Direction : Robert Luce
• Type : Thèse de doctorat
• Discipline(s) : Mathématiques appliquées
• Date : Soutenance le 10/10/2013
• Etablissement(s) : Pau
• Ecole(s) doctorale(s) : École doctorale sciences exactes et leurs applications (Pau, Pyrénées Atlantiques ; 1995-)

Mots clés

Mots clés contrôlés

Éléments finis, Méthode des

Calcul d'erreur

Simulation par ordinateur

Analyse numérique

Équations aux dérivées partielles

Mots clés libres

NXFEM

Généralisation à l’ordre supérieur

Maillages 2D-3D

Méthodes adaptatives

NXFEM

High-order generalization

2D-3D meshes

Adaptive methods

Résumé

Cette thèse porte sur la généralisation de la méthode NXFEM proposée par A. et P. Hansbo pour le problème d’interface elliptique. La modélisation et simulation numérique d’écoulements dans des domaines fracturés sont au coeur de nombreuses applications, telles que le milieu pétrolier (modélisation de réservoirs, présence de failles, propagation d’un signal, repérage de couches), l’aérospatiale (problème de chocs, de rupture), en génie civil (fissuration du béton), mais également dans la biologie cellulaire (déformation des globules rouges). En outre, de nombreux projets de recherche nécessitent le développement des méthodes robustes pour la prise en compte de singularités, ce qui fait partie des motivations et des objectifs de l'équipe Concha, ainsi que de cette thèse. Une modification de cette méthode a tout d’abord été proposée afin d’obtenir la robustesse à la fois par rapport à la géométrie du maillage coupé par l’interface et par rapport aux paramètres de diffusion. Nous nous sommes ensuite intéressés à sa généralisation à tout type de maillages 2D-3D (triangles, quadrilatères, tétraèdres, hexaèdres), et pour tout type d’éléments finis (conformes, non conformes, Galerkin discontinus) pour des interfaces planes et courbes. Les applications ont été orientées vers des problèmes d’écoulements en milieux poreux fracturés : adaptation de la méthode NXFEM à la résolution d’un modèle asymptotique de failles, à des problèmes instationnaires, de transports, ou encore à des domaines multi-fracturés.

This thesis focuses on the generalization of the NXFEM method proposed by A. and P. Hansbo for elliptic interface problem. Numerical modeling and simulation of flow in fractured media are at the heart of many applications, such as petroleum and porous media (reservoir modeling, presence of faults, signal propagation, identification of layers ...), aerospace (problems of shock, rupture), civil engineering (concrete cracking), but also in cell biology (deformation of red blood cells). In addition, many research projects require the development of robust methods for the consideration of singularities, which is one of the motivations and objectives of the Concha team and of this thesis. First a modification of this method was proposed to obtain a robust method not only with respect to the mesh-interface geometry, but also with respect to the diffusion parameters. We then looked to its generalization to any type of 2D-3D meshes (triangles, quadrilaterals, tetrahedra, hexahedra), and for any type of finites elements (conforming, nonconforming, Galerkin discontinuous) for plane and curved interfaces. The applications have been referred to the flow problems in fractured porous media : adaptation of NXFEM method to solve an asymptotic model of faults, to unsteady problems, transport problems, or to multi-fractured domains.