Développement et comparaison de bases d'éléments finis pour la résolution de l'équation de transport sur des maillages hexagonaux

Cette proposition de thèse s'inscrit dans le cadre mathématique de la résolution numérique d'équations aux dérivées partielles par le biais d'une discrétisation des variables dépendantes. Elle s'intéresse, dans un formalisme d'éléments finis, à la construction de schémas d'ordre élevé (au sens de l'espace discrétisé) sur des maillages polygonaux. Le cadre industriel applicatif est la résolution de l'équation de Boltzmann appliquée au transport des neutrons dans le cœur d'un réacteur nucléaire. Dans ce contexte, beaucoup de codes modernes de calcul s'appuient sur une discrétisation par éléments finis (plus précisément, un schéma Galerkin discontinu décentré amont) pour des maillages triangulaires ou quadrilatères du domaine spatial. Cette méthode performante et polyvalente est propice au développement de méthodes de raffinement de maillage adaptatives. L'objectif de cette thèse est l'analyse et le développement d'un tel schéma optimisé pour des maillages en “nid d'abeille” (pavage hexagonal régulier) caractéristiques des réseaux d'assemblages des réacteurs rapides refroidis au sodium. Du point de vue applicatif, cette optimisation est requise pour l'amélioration des simulations multiphysiques coûteuses des transitoires accidentels de ces réacteurs. Pour atteindre cet objectif, on cherchera à mettre au point des bases d'éléments finis d'ordre élevé écrites directement sur l'hexagone. En particulier, on s'intéressera pour ce nouveau schéma au raffinement adaptatif en ordre qui requiert une estimation de l'erreur de discrétisation a posteriori pour permettre d'augmenter localement l'ordre du schéma afin de réduire l'erreur de manière optimale vis-à-vis du coût de calcul. In the realm of applied mathematics, this Ph.D. proposal is related to the numerical solution of partial differential equations through the discretization of the associated dependent variables. In the frame of finite element analysis, it is focused on the construction of high-order schemes on polygonal meshes. The industrial application of such developments is the solution to the Boltzmann equation that describes the transport of neutrons within the core of a nuclear reactor. In this context, many modern codes use a finite-element-based discretization (more precisely, an upwind discontinuous Galerkin scheme) on triangular or quadrilateral meshes spanning the spatial domain. This efficient and versatile method is well suited for the development of adaptive mesh refinement techniques. The object