Analyse mathématique et numérique de modèles gyrocinétiques

Mathematical and numerical analysis of gyro-kinetic models

• Auteur / Autrice : Céline Caldini-Queiros
• Direction : Mihaï Bostan
• Type : Thèse de doctorat
• Discipline(s) : Mathématiques et applications
• Date : Soutenance le 15/11/2013
• Etablissement(s) : Besançon
• Ecole(s) doctorale(s) : École doctorale Sciences pour l'ingénieur et microtechniques (Besançon ; 1991-....)
• Partenaire(s) de recherche : Laboratoire : Laboratoire de Mathématiques de Besançon (Besançon) - Laboratoire de Mathématiques de Besançon
• Jury : Président / Présidente : Alexei Lozinski
• Jury : Examinateurs / Examinatrices : Mihaï Bostan, Alexei Lozinski, Maxime Hauray, Nicolas Besse, Nicolas Crouseilles, Virginie Grandgirard
• Jury : Rapporteurs / Rapporteuses : Eric Sonnendrücker, Emmanuel Frénod

Mots clés

Mots clés contrôlés

Fokker-Planck, Équation de

Mots clés libres

Equations gyro-cinétiques

Equation de Vlasov

Rayon de Larmor fini

Opérateur de Fokker-Planck-Landau

Collisions

Modèles micro-macro

Schéma asymptotiquement préservatif

Gyro-kinetic equations

Vlasov equations

Finite Larmor radius

Guiding-center model

Fokker-Planck-Landau equations

Collisions

Asymptotic preserving schemes

Micro-macro decomposition

Résumé

Cette thèse porte sur les équations gyro-cinétiques et traite un développement rigoureux deslimites de l'équation de Vlasov avec différents opérateurs de collision dans un champ magnétiquefort, ainsi que du développement de méthodes numériques.On commence par une étude de l'opérateur de moyenne. L'opérateur de moyenne a été développé parM. Bostan dans le cadre général d'une équation pour laquelle une partie du transport estfortement pénalisée. Puis, on applique ces résultats généraux aux deux régimes limites que nousétudions : le régime du rayon de Larmor fini et le régime centre-guide.On s'intéresse au calcul précis et explicite de la moyenne de l'opérateur de Fokker-Planck-Landau. On se place pour cela dans le cas du régime du rayon de Larmor fini. Avant de réaliserles calculs sur l'opérateur de Fokker-Planck-Landau, qui contient des convolutions et des termesde diffusion, il semble raisonnable de calculer la moyenne de l'opérateur de relaxation deBoltzmann, dont l'expression est plus simple.On se place ensuite dans le cas du régime centre-guide et on présente un schéma numérique basésur une décomposition micro-macro de la fonction de distribution des particules qui provientd'un travail en collaboration avec N. Crouseilles et M. Lemou. On obtient un schéma uniformémentconsistant avec le modèle continu, pour tout ordre du champ magnétique. Des simulationsnumériques, basées sur cette approche, ont été réalisées à l'aide d'un code de calcul 2D quel'on a développé durant cette thèse.On présente ensuite un projet réalisé dans le cadre du Cemracs 2012, consacré à la modélisationdes écoulements sanguins dans le réseau veineux cérébral.

The main subject of this thesis is the gyro-kinetic equation. We present a rigourousdeveloppement of the Vlasov equation limits with different collision operator in a strongmagnetic field and numerical methods.We start with a study of the gyro-average operator. The average operator has been introduced byM. Bostan in the case of an equation where part of the transport is highly penalised. Then weapply our results at the two approximation we study : the finite Larmor radius approximation andthe guiding-center approximation.We first focus on the precise and explicit computation of the Fokker-Planck-Landau operatoraverage in the finite Larmor radius approximation. The Fokker-Planck-Landau operator containsconvolution and diffusion terms, it is then reasonable to first compute the average of theBoltzmann relaxation operator.We then focus on the guiding-center approximation and present a numerical scheme based on amicro-macro decomposition of the particles distribution fonction which comes from a joint workwith N. Crouseilles and M. Lemou. We obtain a scheme which is uniformly consistant with thecontinuous model for any order of the magnetic field. Numerical simulation based on thisapproach are presented.The last chapter of this thesis presents a project which was realised during the Cemracs 2012concerning the modelisation of blood flow in cerebral veins.