Approximation multi-échelles de l'équation de Vlasov
Langue Français
Ce document est protégé en vertu du Code de la Propriété Intellectuelle.
Auteur : Mouton, Alexandre
Date de soutenance : 16-09-2009
Directeur(s) de thèse : Frénod, Emmanuel - Sonnendrücker, Eric
Établissement de soutenance : Université de Strasbourg
École doctorale : École doctorale Mathématiques, sciences de l'information et de l'ingénieur (Strasbourg ; 1997-....)
Discipline : Mathématiques appliquées et sciences sociales
Classification : Mathématiques, Physique
Mots-clés libres : équation de Vlasov , rayon de Larmor
Mots-clés :
- Analyse numérique
- Thèses et écrits académiques
- Convergence (mathématiques)
- Thèses et écrits académiques
- Transport, Théorie du
- Thèses et écrits académiques
Résumé : Une des difficultés fondamentales de la simulation numérique de plasmas magnétisés est l'existence d'échelles de temps et d'espace mettant en jeu plusieurs ordres de grandeurs très différents. Afin de réaliser des simulations numériques efficaces de ces phénomènes physiques, il est essentiel de développer des modèles et des méthodes numériques adaptés à ces problèmes.A ce jour, la notion de convergence 2-échelles introduite par G. Allaire et G. Nguetseng est un des outils permettant de dériver rigoureusement des limites multi-échelles, ce qui nous permet d'obtenir des modèles limites qu'il est possible de discrétiser avec une méthode numérique usuelle : nous parlons alors d'une méthode numérique 2-échelles.L'objectif de cette thèse est de développer une méthode semi-lagrangienne 2 échelles sur un modèle de type Vlasov gyrocinétique afin de simuler un plasma soumis à un champ magnétique fort du même type que ceux utilisés pour le projet ITER. Cependant, comme les phénomènes physiques à simuler sont assez complexes et comme nous ne savons que peu de choses sur le comportement d'une méthode numérique 2-échelles sur un modèle non-linéaire, il convient de procéder par étapes avant de développer une telle méthode sur un modèle gyrocinétique.Dans une première partie, nous construisons une méthode de volumes finis 2-échelles sur les équations d'Euler 1D isentropiques faiblement compressibles. Bien que ce modèle soit assez différent d'un modèle de type Vlasov, il n'en est pas moins un cadre de travail relativement simple pour étudier le comportement d'une méthode numérique 2-échelles face à un modèle non-linéaire.Dans une seconde partie, nous nous basons sur le modèle limite développé par E. Frénod, F. Salvarani et E. Sonnendrücker afin de construire une méthode semi-lagrangienne 2-échelles pour simuler des faisceaux de particules en géométrie axisymétrique. Même si le modèle de Vlasov axisymétrique utilisé est différent d'un modèle gyrocinétique, il constitue un contexte idéal pour établir les bases d'une méthode semi-lagrangienne 2 échelles.Enfin, dans une troisième partie, nous utilisons la convergence 2-échelles afin d'améliorer les résultats de convergence faible-* établis par M. Bostan en 2007, et nous proposons une méthode semi-lagrangienne en avant permettant de valider numériquement ces résultats.
Entrepôt d'origine : 
Identifiant : ecrin-ori-298109
Type de ressource : Thèse
