Section: Scientific Foundations

Kinetic models for plasma and beam physics

plasma physics, beam physics, kinetic models, reduced models, Vlasov equation, modeling, mathematical analysis, asymptotic analysis, existence, uniqueness

Plasmas and particle beams can be described by a hierarchy of models including N-body interaction, kinetic models and fluid models. Kinetic models in particular are posed in phase-space and involve specific difficulties. We perform a mathematical analysis of such models and try to find and justify approximate models using asymptotic analysis.

Models for plasma and beam physics

The plasma state can be considered as the fourth state of matter, obtained for example by bringing a gas to a very high temperature (10 4 K or more). The thermal energy of the molecules and atoms constituting the gas is then sufficient to start ionization when particles collide. A globally neutral gas of neutral and charged particles, called plasma, is then obtained. Intense charged particle beams, called nonneutral plasmas by some authors, obey similar physical laws.

The hierarchy of models describing the evolution of charged particles within a plasma or a particle beam includes N-body models where each particle interacts directly with all the others, kinetic models based on a statistical description of the particles and fluid models valid when the particles are at a thermodynamical equilibrium.

In a so-called kinetic model, each particle species s in a plasma or a particle beam is described by a distribution function f s (𝐱,𝐯,t) corresponding to the statistical average of the particle distribution in phase-space corresponding to many realisations of the physical system under investigation. The product f s d𝐱d𝐯 is the average number of particles of the considered species, the position and velocity of which are located in a bin of volume d𝐱d𝐯 centered around (𝐱,𝐯). The distribution function contains a lot more information than what can be obtained from a fluid description, as it also includes information about the velocity distribution of the particles.

A kinetic description is necessary in collective plasmas where the distribution function is very different from the Maxwell-Boltzmann (or Maxwellian) distribution which corresponds to the thermodynamical equilibrium, otherwise a fluid description is generally sufficient. In the limit when collective effects are dominant with respect to binary collisions, the corresponding kinetic equation is the Vlasov equation

βˆ‚f s βˆ‚t+π―Β·βˆ‚f s βˆ‚π±+q m(𝐄+𝐯×𝐁)Β·βˆ‚f s βˆ‚π―=0,

which expresses that the distribution function f is conserved along the particle trajectories which are determined by their motion in their mean electromagnetic field. The Vlasov equation which involves a self-consistent electromagnetic field needs to be coupled to the Maxwell equations in order to compute this field

-1 c 2 βˆ‚π„ βˆ‚t+βˆ‡Γ—π=ΞΌ 0 𝐉,βˆ‚π βˆ‚t+βˆ‡Γ—π„=0,÷𝐄=ρ Ξ΅ 0 ,÷𝐁=0,

which describes the evolution of the electromagnetic field generated by the charge density

ρ(𝐱,t)=βˆ‘ s q s ∫f s (𝐱,𝐯,t)d𝐯,

and current density

𝐉(𝐱,t)=βˆ‘ s q s ∫f s (𝐱,𝐯,t)𝐯d𝐯,

associated to the charged particles.

When binary particle-particle interactions are dominant with respect to the mean-field effects then the distribution function f obeys the Boltzmann equation

βˆ‚f βˆ‚t+π―Β·βˆ‚f βˆ‚π±=Q(f,f),

where Q is the nonlinear Boltzmann collision operator. In some intermediate cases, a collision operator needs to be added to the Vlasov equation.

The numerical solution of the three-dimensional Vlasov-Maxwell system represents a considerable challenge due to the huge size of the problem. Indeed, the Vlasov-Maxwell system is nonlinear and posed in phase space. It thus depends on seven variables: three configuration space variables, three velocity space variables and time, for each species of particles. This feature makes it essential to use every possible option to find a reduced model wherever possible, in particular when there are geometrical symmetries or small terms which can be neglected.

Mathematical and asymptotic analysis of kinetic models

The mathematical analysis of the Vlasov equation is essential for a thorough understanding of the model as well for physical as for numerical purposes. It has attracted many researchers since the end of the 1970s. Among the most important results which have been obtained, we can cite the existence of strong and weak solutions of the Vlasov-Poisson system by Horst and Hunze [74] , see also Bardos and Degond [51] . The existence of a weak solution for the Vlasov-Maxwell system has been proved by Di Perna and Lions [59] . An overview of the theory is presented in a book by Glassey [71] .

Many questions concerning for example uniqueness or existence of strong solutions for the three-dimensional Vlasov-Maxwell system are still open. Moreover, their is a realm of approached models that need to be investigated. In particular, the Vlasov-Darwin model for which we could recently prove the existence of global solutions for small initial data [52] .

On the other hand, the asymptotic study of the Vlasov equation in different physical situations is important in order to find or justify reduced models. One situation of major importance in tokamaks, used for magnetic fusion as well as in atmospheric plasmas, is the case of a large external magnetic field used for confining the particles. The magnetic field tends to incurve the particle trajectories which eventually, when the magnetic field is large, are confined along the magnetic field lines. Moreover, when an electric field is present, the particles drift in a direction perpendicular to the magnetic and to the electric field. The new time scale linked to the cyclotron frequency, which is the frequency of rotation around the magnetic field lines, comes in addition to the other time scales present in the system like the plasma frequencies of the different particle species. Thus, many different time scales as well as length scales linked in particular to the different Debye length are present in the system. Depending on the effects that need to be studied, asymptotic techniques allow to find reduced models. In this spirit, in the case of large magnetic fields, recent results have been obtained by Golse and Saint-Raymond [72] , [80] as well as by Brenier [57] . Our group has also contributed to this problem using homogenization techniques to justify the guiding center model and the finite Larmor radius model which are used by physicist in this setting [67] , [65] , [66] .

Another important asymptotic problem yielding reduced models for the Vlasov-Maxwell system is the fluid limit of collisionless plasmas. In some specific physical situations, the infinite system of velocity moments of the Vlasov equations can be closed after a few of those, thus yielding fluid models.