Adaptive Space-Time Finite Element Method in High Temperature Superconductivity

Cette thèse porte sur le développement d'une méthode d'éléments finis adaptative pour discrétiser un modèle électromagnétique issu du domaine de la supraconductivité à haute température critique. Dans les faits, ce modèle consiste en une version non-linéaire du problème de courants de Foucault classique en électromagnétisme, dans lequel la non-linéarité engendre un état mixte de régions supraconductrices et de régions normales. La méthode développée dans cette thèse peut être appliquée directement à la conception et à l'optimisation de dispositifs supraconducteurs à haute température critique, sans s'y limiter. Le problème mathématique correspondant à ce modèle, que l'on appelle p-rotationnel (p-curl en anglais), est une équation différentielle aux dérivées partielles de type évolutionnaire monotone dont les solutions de la forme faible appartiennent à un espace fonctionnel de Sobolev de type L^p. Étant donnée la présence de singularités dans les solutions, les méthodes numériques développées à ce jour pour résoudre ce probléme se généralisent mal aux domaines 2D et 3D. La principale difficulté provient de l'absence d'estimateur d'erreur qui permettrait un contrôle adaptatif de la finesse du maillage et du pas de temps. Cette thèse présente deux contributions principales. Premièrement, nous avons développé et implémenté une méthode d'éléments finis adaptative basée sur une formulation espace-temps. Afin de faciliter l'adaptivité, nous avons conçu une structure arborescente espace-temps, qui se construit de façon récursive en fonction du raffinement tout en préservant l'irrégularité de premier niveau (1-irregularity) du maillage espace-temps. De plus, nous avons développé un opérateur d'interpolation qui permet de préserver la continuité des degrés de liberté sur les arêtes “sans voisins” (hanging edges). La seconde contribution principale de la thèse est le développement d'un estimateur d'erreur a posteriori basé sur le résidu, et nous avons prouvé mathématiquement sa fiabilité dans le cas semi-discret. Un élément clé de cette preuve fut d'utiliser une nouvelle version de la décomposition d'Helmholtz pour l'espace 〖W_0〗^p (curl;Ω), requis pour démonter une variante de l'orthogonalité de Galerkin. La fiabilité d'une grandeur physique d'intérêt, appelée pertes AC, a aussi été démontrée. Des résultats numériques en 1D et 2D sont aussi présentés dans les cas uniformes et adaptatifs. En bref, cette recherche se distingue des travaux précédents sur le p-rotationnel parce qu'elle se base sur une analyse mathématique théorique rigoureuse pour guider le développement et l'analyse des nouvelles techniques proposées.