Eléments finis courbes et accélération pour le transport de neutrons

• Main Author: Moller, Jean-Yves
• Language: FRE
• Published: Université Henri Poincaré - Nancy I 2012
• Subjects: [MATH:MATH_MP] Mathematics/Mathematical Physics; [MATH:MATH_NA] Mathematics/Numerical Analysis; transport neutronique; éléments finis; mathématiques appliquées; accélération synthétique
• Online Access: http://tel.archives-ouvertes.fr/tel-00659200; http://tel.archives-ouvertes.fr/docs/00/65/92/00/PDF/these.pdf

La modélisation des réacteurs nucléaires repose sur la résolution de l'équation de Boltzmann linéaire. Nous nous sommes intéressés à la résolution spatiale de la forme stationnaire de cette équation. Après discrétisation en énergie et en angle, l'équation hyperbolique est résolue numériquement par la méthode des éléments finis discontinus. Le solveur MINARET utilise cette méthode sur un maillage triangulaire non structuré afin de pouvoir traiter des géométries complexes (comprenant entre autres des arcs de cercle). Cependant, l'utilisation d'arêtes droites introduit une approximation de la géométrie. Autoriser l'existence d'arêtes courbes permet de coller parfaitement à la géométrie, et dans certains cas de diminuer le nombre de triangles du maillage. L'objectif principal de cette thèse est l'étude d'éléments finis sur des triangles possédant un ou plusieurs bords courbes. Le choix des fonctions de base est un des points importants pour ce type d'éléments finis. Un résultat de convergence a été obtenu sous réserve que les triangles courbes ne soient pas trop éloignés des triangles droits associés. D'autre part, un solveur courbe a été développé pour traiter des triangles avec un, deux ou trois bords courbes. Une autre partie de ce travail porte sur l'accélération de la convergence des calculs. En effet, la résolution du problème est itérative et peut, dans certains cas, converger très lentement. Une méthode d'accélération dite DSA (Diffusion Synthetic Acceleration) permet de diminuer le nombre d'itérations et le temps de calcul : un calcul de diffusion est ajouté à chaque itération. L'opérateur de diffusion est un préconditionneur de l'opérateur de transport. La DSA a été mise en oeuvre en utilisant une technique issue des méthodes de pénalisation intérieure. Une analyse de Fourier en 1D et 2D permet d'évaluer l'accélération dans le cas de milieux infinis périodiques et de vérifier la stabilité du schéma lorsque de fortes hétérogénéités existent.