Study of reduced kinetic models for plasma turbulence
Étude de modèles réduits cinétiques de la turbulence plasma
Résumé
Turbulent transport is one of the keys to improve the energy confinement time required for thermonuclear fusion reactors. The description of the kinetic turbulence of the plasma is a problem with 3 spatial coordinates and 3 velocity coordinates. Both theory and simulation of a problem of such high dimensionality are very difficult, and reduced models are helpfull to understand turbulence in Tokamaks. A widely used technique consists into averaging the cyclotron motion, which is much faster than the turbulence time scale. Such a reduction makes it possible to simplify the problem to three spatial coordinates of the particle guide centers, a parallel velocity or energy, and a perpendicular velocity appearing as the adiabatic invariant. Nonlinear gyrokinetic description requires massively parallel high performance numerical simulations. The difficulty lies in the non-linear terms (Poisson hooks) that describe multi-scale interactions, which is a challenge for both theory and simulation. Any reduced approach, based on well-controlled hypotheses, is therefore interesting to develop.On the basis of this ambition, this thesis concerns the turbulence of particles trapped in magnetized plasma. It is a 4D system, obtained after averaging the particle distribution function on cyclotron and bounce motions, which can be considered as a reduced form of standard gyrokinetic theory. We called it "bounce averaged gyrokinetics" during this work. Even if this description is greatly reduced compared to the gyrokinetic theory, nonlinear direct simulation remains a challenge.A description of the nonlinear polar coordinate terms is chosen, with a logarithmic grid along the norm of the wave vector, while the angles are discretized on a regular grid. The use of a logarithmic grid makes it possible to take into account a wide range of wave vectors, so physics on a very small scale. In a similar way to shell models for fluid turbulence, and in order to simplify the system, only the interactions between neighboring shells are considered.In a first step, the study of the linear system is presented, in particular the paraetric dependence of the instability thresholds and the linear growth rate, allowing to recover the strong anisotropy of the growth rates of the trapped ion modes (or TIM) and the modes of trapped electrons (or TEM). These studies also make it possible to validate the non-linear numerical codes with respect to an independently developer eigenvalue solver.In a second step, the isotropic hypothesis for nonlinear terms is used. Thus, there is no exact phase information for such 1D layer models, which leaves with a free parameter in the interaction coefficients. An original power law is evidenced, which is unaffected by the value of the free parameter, measuring the intensity of the nonlinear effects relative to the linear terms.From the simulation of the isotropic model, the phase information appears very important. Since the linear instability is anisotropic for the fusion, the simulation of the anisotropic model is thus carried out in a third time. The numerically resolved system is reduced to a kinetic species, assuming that the other species are adiabatic. Two different systems can thus be studied: kinetic ions + adiabatic electrons and kinetic electrons + adiabatic ions. Different spectra are observed in each of these two cases, and the validity of the adiabatic hypothesis is discussed for each species, based on a kinetic simulation with two species.
Le contrôle du transport turbulent est l'une des clés pour l'amélioration du temps de confinement nécessaire à la réalisation de la fusion thermonucléaire contrôlée. La description de la turbulence cinétique du plasma est un problème à 3 coordonnées spatiales et 3 coordonnées en vitesse. La théorie comme la simulation pour un problème de si haute dimensionnalité sont très difficiles, et des modèles réduits sont nécessaires pour comprendre la turbulence dans les Tokamaks. La technique largement utilisée est de faire moyenner le mouvement cyclotron, qui est beaucoup plus rapide que le phénomène de turbulence. Une telle réduction permet de simplifier le problème à trois coordonnées spatiales des centres-guides des particules, une vitesse parallèle ou énergie et une vitesse perpendiculaire apparaissant comme l'invariant adiabatique. La description gyrocinétique non linéaire requiert des simulations numériques de haute performance massivement parallèles. Toute la difficulté est due aux termes non linéaires (crochets de Poisson) qui décrivent les interactions multi-échelles, ce qui constitue un défi tant pour la théorie que pour la simulation. Toute approche réduite, basée sur des hypothèses bien contrôlées, est donc intéressante à développer.Sur la base de cette ambition, cette thèse concerne la turbulence des particules piégées dans le plasma magnétisé. C'est un système 4D, obtenu après avoir fait la moyenne de la fonction de distribution des particules sur les mouvements cyclotron et de rebond, ce qui peut être considéré comme une forme réduite de la théorie gyrocinétique standard. Nous l'avons appelé "bounce averaged gyrokinetics" pendant ce travail. Même si cette description est grandement réduite par rapport à la théorie gyrocinétique, la simulation directe non-linéaire reste un challenge.Une description des termes non linéaires en coordonnées polaires est choisie, avec une grille logarithmique en norme du vecteur d'onde, tandis que les angles sont discrétisés sur une grille régulière. L'utilisation d'une grille logarithmique permet de prendre en compte une large gamme de vecteurs d'ondes, donc la physique à très petite échelle. De manière analogue aux modèles en couches en turbulence fluide et afin de simplifier le système, seules les interactions entre couches voisines sont considérées.Dans un premier temps, l'étude du système linéaire est présentée, en particulier les seuils des paramètres et l'instabilité linéaire permettant de retrouver la forte anisotropie des taux de croissance des modes d'ions piégés (ou TIM) et des modes d'électrons piégés (ou TEM). Ces études permettent également de valider les codes numériques non-linéaires vis-à-vis d'un solveur aux valeurs propres développé indépendamment.Dans un second temps, l'hypothèse isotrope pour les termes non linéaires est utilisée. Ainsi il n'y a pas d'information de phase exacte pour de tels modèles en couches 1D, ce qui laisse un paramètre libre dans les coefficients d'interaction. Une loi de puissance originale est mise en évidence, qui n'est pas affectée par la valeur du paramètre libre, mesurant l'intensité des effets non-linéaires relativement aux termes linéaires.À partir de la simulation du modèle isotrope, l'information de phase apparaît très importante. Puisque l'instabilité linéaire est anisotrope pour la fusion, la simulation du modèle anisotrope est donc réalisée dans un troisième temps. Le système résolu numériquement est réduit à une espèce cinétique, en supposant que les autres espèces sont adiabatiques. Deux systèmes différents peuvent ainsi être étudiés: ions cinétiques + électrons adiabatiques et électrons cinétiques + ions adiabatiques. Des spectres différents sont observés dans chacun de ces deux cas, et la validité de l'hypothèse adiabatique est discutée pour chaque espèce, avec pour base de comparaison une simulation cinétique à deux espèces.
Origine : Version validée par le jury (STAR)
Loading...