Optimisation de méthodes numériques pour la physique des plasmas : application aux faisceaux de particules chargées

par Anaïs Crestetto

Thèse de doctorat en Mathématiques appliquées et sciences sociales

Sous la direction de Philippe Helluy.

Le président du jury était Francis Filbet.

Le jury était composé de Eric Sonnendrucker, Nicolas Besse, Nicolas Crouseilles.

Les rapporteurs étaient Xavier Ferrieres.


  • Résumé

    Cette thèse propose différentes méthodes numériques pour simuler les plasmas ou les faisceaux de particules chargées à coût réduit. Le mouvement de particules chargées soumises à un champ électromagnétique est régi par l'équation de Vlasov, couplée aux équations de Maxwell ou de Poisson. Dans la première partie, une méthode multi-fluides est utilisée pour la résolution du système de Vlasov-Poisson 1D. Elle est basée sur la connaissance a priori de la forme prise par la fonction de distribution f. Ce type de méthodes est plutôt adapté aux systèmes restant proches de l'état d'équilibre. La deuxième partie propose de décomposer f en une partie d'équilibre et une perturbation. L'équilibre est résolu par une méthode fluide, la perturbation par une méthode cinétique plus précise. On construit un schéma préservant l'asymptotique pour le système de Vlasov-Poisson-BGK basé sur une telle décomposition. On étudie dans la troisième partie la méthode PIC en géométrie 2D axisymétrique. Un travail basé sur l'analyse isogéométrique est présenté ainsi qu'un code PIC - Galerkin Discontinu parallélisé sur carte graphique.

  • Titre traduit

    Optimisation of numerical methods for plasma physics : application to charged particle beams


  • Résumé

    This thesis presents different numerical methods for the simulation of plasmas or charged particles beams with reduced cost. Movement of charged particles in an electromagnetic field is given by the Vlasov equation, coupled to the Maxwell equations for the electromagnetic field, or to the Poisson equation. In the first part, a multi-fluid method is used for solving the 1D Vlasov-Poisson system. It is based on the a priori knowledge of the shape of f. This kind of methods is rather adapted to systems staying close to the equilibrium. The second part presents the decomposition of f between an equilibrium part and a perturbation. The equilibrium part is solved by a fluid method whereas we use a kinetic method for the perturbation. We construct an asymptotic preserving scheme for the Vlasov-Poisson-BGK system using such a decomposition. The third part deals with the PIC method in 2D axisymmetric geometry. A work based on isogeometric analysis is presented, and then a PIC - Discontinuous Galerkin program computed on graphic card.


Il est disponible au sein de la bibliothèque de l'établissement de soutenance.

Consulter en bibliothèque

La version de soutenance existe

Où se trouve cette thèse ?

  • Bibliothèque : Université de Strasbourg. Service commun de la documentation. Bibliothèque électronique 063.
Voir dans le Sudoc, catalogue collectif des bibliothèques de l'enseignement supérieur et de la recherche.