A Crestetto
email: anais.crestetto@univ-nantes.fr

F Deluzet

D Doyen

Bridging kinetic plasma descriptions and single fluid models

Keywords: Plasma, Debye length, MHD, drift limit, fluid limit, Quasineutrality, Vlasov-Maxwell, Asymptotic-Preserving scheme

The purpose of this paper is to bridge kinetic plasma descriptions and low frequency single fluid models. More specifically, the asymptotics leading to Magneto-Hydro-Dynamic (MHD) regimes starting 1 from the Vlasov-Maxwell system are investigated. The analogy with the derivation, from the Vlasov-Poisson system, of a fluid representation for the ions coupled to the Boltzmann relation for electrons is also outlined. The aim is to identify asymptotic parameters explaining the transitions from one microscopic description to a macroscopic low frequency model. These investigations provide ground work for the derivation of multi-scale numerical methods, model coupling or physics based preconditioning.

Introduction

The aim of this paper is to propose a continuation of the work initiated in [START_REF] Degond | An asymptotically stable particle-in-cell (PIC) scheme for collisionless plasma simulations near quasineutrality[END_REF][START_REF] Degond | Asymptotic-preserving particle-in-cell method for the Vlasov-Poisson system near quasineutrality[END_REF][START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF] focusing on the derivation of asymptotic preserving schemes for kinetic plasma descriptions in the quasi-neutral limit. The purpose of these numerical methods is to provide a quasi-neutral description of the plasma with no constraints on the simulation parameters related to the Debye length but with the ability to perform local up-scalings with non neutral plasma descriptions. This brings a gain in the computational efficiency, since the discretization parameters can be set according to the physics of interest rather than the small scales (namely the Debye length) described by the model.

The methodology introduced in these former achievements is aimed to be generalized here to more singular limits. In this series of former works, the limit models remain kinetic and the scales of interest are related to the electron dynamics. For instance, the quasi-neutral limit of the Vlasov-Maxwell system investigated in [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF] can be interpreted as a kinetic description of the Electron-Magneto-Hydro-Dynamic (E-MHD) [START_REF] Gordeev | Electron magnetohydrodynamics[END_REF][START_REF] Swanekamp | Particleincell simulations of fast magnetic field penetration into plasmas due to the hall electric field[END_REF][START_REF] Cho | The Anisotropy of Electron MHD Turbulence[END_REF], accounting for the electron inertia, the massive ions being assumed at rest or slowly evolving. In the present paper, the objective is to go beyond the kinetic E-MHD with the aim to bridge the Vlasov-Maxwell system and Magneto-Hydro-Dynamic (MHD) models. In MHD systems, the scales of interest are defined by the overall plasma dynamic which is governed by the ions, the fast scales associated to the electron inertia being filtered out from the equations.

The present work is therefore devoted to the derivation of a model hierarchy bridging either the Vlasov-Maxwell system and MHD models for mag-netized plasmas or, the Vlasov-Poisson and the electron adiabatic response also referred to as Boltzmann relation (see [START_REF] Langmuir | The interaction of electron and positive ion space charges in cathode sheaths[END_REF][START_REF] Tonks | Oscillations in ionized gases[END_REF][START_REF] Tonks | A general theory of the plasma of an arc[END_REF] for seminal works and [START_REF] De Cecco | Asymptotic Transition from Kinetic to Adiabatic Electrons along Magnetic Field Lines[END_REF] for numerical investigations), for electrostatic frameworks. A wide range of applications of the present investigations can be named, specifically low variance Particle-In-Cell methods or more generally numerical discretization of kinetic models implementing a Micro-Macro decomposition of the distribution function. We refer for instance to [START_REF] Crestetto | Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles[END_REF][START_REF] Crouseilles | An asymptotic preserving scheme based on a micro-macro decomposition for Collisional Vlasov equations: diffusion and high-field scaling limits[END_REF][START_REF] Dimarco | An asymptotic preserving automatic domain decomposition method for the VlasovPoissonBGK system with applications to plasmas[END_REF][START_REF] Lemou | A New Asymptotic Preserving Scheme Based on Micro-Macro Formulation for Linear Kinetic Equations in the Diffusion Limit[END_REF]) for Micro-Macro methods, and to [START_REF] Degond | The moment-guided Monte Carlo method[END_REF] for the moment Guided method; fluid-preconditioned fully implicit methods [START_REF] Chen | An energy-and chargeconserving, implicit, electrostatic particle-in-cell algorithm[END_REF][START_REF] Chen | Fluid preconditioning for newtonkrylov-based, fully implicit, electrostatic particle-in-cell simulations[END_REF][START_REF] Chen | A multi-dimensional, energy-and chargeconserving, nonlinearly implicit, electromagnetic VlasovDarwin particlein-cell algorithm[END_REF]] and Asymptotic-Preserving numerical methods [START_REF] Jin | Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations[END_REF][START_REF] Degond | Asymptotic-Preserving methods and multiscale models for plasma physics[END_REF]. In these methods, macroscopic models are operated to either filter out the numerical noise of PIC methods (moment guided, Micro-Macro particle methods) or speed up the convergence of non linear implicit solvers (fully implicit PIC methods). The efficiency of these methods is closely related to that of the reduced models and its adequacy with the kinetic plasma description. Another application can be envisioned with the hybrid coupling of Particle-In-Cell methods and MHD descriptions [START_REF] Schumer | MHD-to-PIC transition for modeling of conduction and opening in a plasma opening switch[END_REF][START_REF] Daldorff | Two-way coupling of a global Hall magnetohydrodynamics model with a local implicit particle-in-cell model[END_REF] and more generally coupling strategies such as the Current-Coupling-Scheme (CCS) and the Pressure-Coupling-Scheme (PCS) (see [START_REF] Park | Three-dimensional hybrid gyrokineticmagnetohydrodynamics simulation[END_REF] [START_REF] Tronci | Hybrid Vlasov-MHD models: Hamiltonian vs. non-Hamiltonian[END_REF] and the references therein).

The aim here is to clarify how the asymptotic parameters interact with each other and define reduced models, but also, to relate these parameters to meaningful physical quantities. The MHD regime is sometimes derived by letting ε 0 , the vacuum permittivity, go to zero (see for instance [START_REF] Jang | Derivation of Ohm's Law from the Kinetic Equations[END_REF][START_REF] Tronci | Neutral Vlasov kinetic theory of magnetized plasmas[END_REF]) which is referred to as the full Maxwell to the low frequency pre-Maxwell's equations asymptotic in [26, see section 2. 3.3]. It is also common to let the electron to ion mass ratio go to zero to explain the vanishing of the electron inertia [START_REF] Freidberg | Ideal MHD[END_REF][START_REF] Klingenberg | A consistent Kinetic Model for a Two-Component Mixture with an application to plasma[END_REF] in deriving either MHD modelling or the Boltzmann relation. Although the right asymptotic models are recovered by this means, these assumptions do not account for changes in the system characteristics that may explain for a regime transition: the electron to ion mass ratio remains constant and the same property holds true for the vacuum permittivity.

The outlines of the paper are the following. The plasma kinetic description is introduced in Sec. 2 together with the Maxwell system driving the evolution of the electromagnetic field. A dimensionless form of the system is stated in order to develop an asymptotic analysis and the derivation of reduced models. A hierarchy of quasi-neutral model is proposed in Sec. 3 for the Vlasov-Maxwell system. It encompasses fully kinetic, hybrid as well as single fluid (MHD) plasma descriptions. The electrostatic framework is investigated in Sec. 4. The electrostatic limit of the Maxwell system is performed.

A hierarchy of models, similar to that of the electromagnetic framework is derived. Finally, a synthesis of these asymptotic analysis is proposed in Sec. 5 devoted to conclusions.

2 The Vlasov-Maxwell system in a dimensionless form

Objectives

In this section, the purpose is to unravel a series of asymptotic limits bridging the gap between the Vlasov-Maxwell system and a Magneto-Hydro-Dynamic (MHD) model. The difficulty is therefore to identify parameters explaining the transition from one description to the other and to relate these parameters to specific characteristics of the system. The tools mobilized to meet this aim are based on the asymptotic analysis of the Vlasov-Maxwell system. Since the low frequency plasma modelling is related to a fluid plasma description, the kinetic model is upgraded with collision operators. Therefore, the most refined modelling consists of a Vlasov equation for the electrons and the ions, augmented with a collision operator and coupled to the Maxwell system. Even if the physical model is non collisional or weakly collisional, the transition towards a fluid limit is accounted for by a collisional process, thanks to a BGK operator. This choice of collision operator is questionable from a strict modelling view point, nonetheless, the purpose here is to easily derive the fluid limit at a limited computational cost. In this respect the BGK collision operator is a good candidate. First, the whole collisional processes are considered, including both inner and inter species collisions. Nonetheless, only the minimal collisional process will be accounted for to derive a MHD regime from the kinetic model. This point will be outlined in the following sections. The introduction of non dimensional quantities will naturally reveal dimensionless parameters in the equations. Letting some of these parameters go to zero shapes the hierarchy of models derived for the Vlasov-Maxwell system and bridging the gap with MHD models.

The Vlasov-BGK-Maxwell system

The most refined description of the plasma is constituted by two Vlasov equations, f i and f e being the ion and electron distribution functions

∂ t f i + v • ∇ x f i + q m i (E + v × B) • ∇ v f i = Q i , (1)
∂ t f e + v • ∇ x f e - q m e (E + v × B) • ∇ v f e = Q e . (2)
In these equations, q is the elementary charge, m α is the mass of the species α (α = e for the electrons and i for the ions). The BGK collision operator Q α are given by [31]

Q α = Q αα + Q αβ , Q αα = ν αα (M nα,uα,Tα -f α ) , Q αβ = ν αβ M nα,u β ,T β -f α , (3)
ν αα and ν αβ being the like-particle and inter species collision frequencies which can be defined as [START_REF] Degond | Chapter One -Asymptotic continuum models for plasmas and disparate mass gaseous binary mixtures[END_REF][START_REF] Spatschek | High temperature plasmas: theory and mathematical tools for laser and fusion plasmas[END_REF]

ν ii = K 0 n i (k B T i ) 3 2 √ 2 √ m i , ν ie = K 0 n e (k B T i ) 3 2 √ m e m i , (4a)
ν ee = K 0 n e (k B T e ) 3 2 √ 2 √ m e , ν ei = K 0 n i (k B T e ) 3 2 1 √ m e , (4b)
where

K 0 = C q 2 4πǫ 0 2 ln(Λ) , (4c)
C denoting a constant with a magnitude equal to one, ln(Λ) the Coulomb logarithm and the ions being assumed mono-charged. The Maxwellians M nα,uα,Tα and M nα,u β ,T β are defined as

M nα,uα,Tα = n α (x, t) m α 2πk B T α (x, t) Dv 2 exp - m α |u α (x, t) -v| 2 2k B T α (x, t) , (5a)
M nα,u β ,T β = n α (x, t) m α 2πk B T β (x, t) Dv 2 exp - m α |u β (x, t) -v| 2 2k B T β (x, t) , (5b)
D v denoting the dimension of the velocity space, k B is the Boltzmann constant. The Maxwellian parameters are n α , u α and T α the density, mean velocity and temperature associated to the distribution function f α and defined as

n α = Ωv f α dv , n α u α = Ωv v f α dv , 1 γ -1 n α k B T α = m α 2 Ωv |v-u| 2 f α dv , (6)
with γ the specific heat ratio whose value depends on the dimensionality of the velocity space

D v through γ -1 = 2 D v . (7)
The collision operators verify the following conservation properties

Q αα m α   1 v |v| 2   dv = 0 , (8a)
Q αβ m α   1 v |v| 2   dv + Q βα m β   1 v |v| 2   dv = 0 . (8b)
The temperature and the mean velocity (u β , T β ) in the inter-species collision operator expression (5b) should be chosen with care in order to guarantee the total momentum and energy conservation. Indeed the following identities

Q αβ m α v dv = ν αβ m α n α (u β -u α ) , (9a)
Q αβ m α |v| 2 2 dv = ν αβ D v 2 n α k B T β -T α + 1 2 m α n α |u β | 2 -|u α | 2 (9b)
hold true for the operators defined by [START_REF] Chen | Fluid preconditioning for newtonkrylov-based, fully implicit, electrostatic particle-in-cell simulations[END_REF]. The trivial choice (u β , T β ) = (u β , T β ) does ensure the plasma total momentum conservation, provided that ν ei m e n e = ν ie m i n i . However in this case, the plasma total energy is not conserved. We refer to [START_REF] Greene | Improved Bhatnagar-Gross-Krook model of electron-ion collisions[END_REF] for a seminal work, as well as [START_REF] Klingenberg | A consistent Kinetic Model for a Two-Component Mixture with an application to plasma[END_REF] and the references therein for recent advances, on the choice of these parameters compliant with the desired properties (8b) of the inter-species collision operators therefore providing the conservation of the plasma momentum and energy.

The electromagnetic field (E, B) evolution is driven by the Maxwell system:

1 c 2 ∂ t E -∇ x × B = -µ 0 J , (10)
∂ t B + ∇ x × E = 0 , (11)
∇ x • E = ρ ǫ 0 , (12)
∇ x • B = 0 , ( 13
)
where c is the speed of light, µ 0 the vacuum permeability and ǫ 0 the vacuum permittivity verifying µ 0 ǫ 0 c 2 = 1. The Maxwell sources are the particle currents and densities

ρ = q(n i -n e ) , (14a)
J = q(n i u i -n e u e ) . (14b)
The definition of the collision frequencies as stated by Eqs.(4) relates different time scales. Indeed, because of their different masses, ions and electrons are not equally affected by collisions. This properties are more clearly emphasized working with dimensionless variables as proposed in the next section.

Scaling of the Vlasov-Maxwell system

The equations are written with dimensionless quantities in order to easily identify different regimes. The scaling is introduced under a priori assumptions that the electronic and ionic temperatures, densities and mean velocities are comparable with a magnitude denoted T 0 , n 0 and u 0 . These scales define the typical Debye length as well as the electron plasma period

λ D = ǫ 0 k B T 0 q 2 n 0 , τ pe = m e ǫ 0 q 2 n 0 .
We denote by x 0 and t 0 the characteristic space and time scales of the phenomena observed, which yields to the velocity of interest ϑ 0 = x 0 /t 0 . The magnitude of the thermal velocity for the species α is denoted v 0,α with v 0,α = k B T 0 /m α . Due to the different masses, the thermal velocity of the electron is not that of the ions. The reference thermal velocity v 0 will be defined by the ion one v 2 0 = k B T 0 /m i , hence v 0,e = v 0 /ε and v 0,i = v 0 , where ε 2 = m e /m i . Finally, the scale of the electromagnetic field is denoted (E 0 , B 0 ), the particle current scale being defined as J 0 = qn 0 u 0 . The dimensionless variables are defined according to

x * = x x 0 , t * = t t 0 , v * = v v 0,α , f * = f n 0 /(v 0,α ) Dv , n * = n n 0 , J * = J qn 0 u 0 , E * = E E 0 , B * = B B 0 ,
the collision frequencies verifying

ν ee,0 = ν ei,0 = 1 ε ν ii,0 , ν ie,0 = εν ii,0 , ε = m e m i . (15)
On the fastest time scales, the electron distribution function relaxes towards a Maxwellian. On the same time scale, the electron mean velocity and temperature relax towards those of the ions. The relaxation of the ionic distribution function towards the local equilibrium is slower, by a factor ε -1 = m i /m e . Finally, the ions are almost unaffected by the collisions against the electrons. The relaxation of the ionic distribution function towards that of the electrons define the largest time scale, by a factor ε -1 compared to the the relaxation towards the thermodynamical equilibrium.

The dimensionless ionic and electronic Vlasov equations can be rewritten as (keeping the same notations for dimensionless variables):

ξ∂ t f i + v • ∇ x f i + η(E + β ξ v × B) • ∇ v f i = ξ κ ν ii (M n i ,u i ,T i -f i ) + εν ie M n i ,ue,Te -f i , (16)
ξε∂ t f e + v • ∇ x f e -η(E + β εξ v × B) • ∇ v f e = ξ κ ν ee (M ne,ue,Te -f e ) + ν ei M ne,u i ,T i -f e , (17)
together with the dimensionless Maxwell system writting

λ 2 η(α 2 ∂E ∂t -β∇ x × B) = -α 2 M ξ J , (18a)
β∂ t B + ∇ x × E = 0 , (18b) λ 2 η∇ x • E = n i -n e , (18c)
∇ x • B = 0 , (18d) J = n i u i -n e u e . ( 18e
)
This system is written thanks to the following dimensionless parameters

                                                   ε 2 = m e m i
the ratio of the electronic and ionic masses ,

λ = λ D x 0 the scaled Debye length , M = u 0 v 0 the ionic Mach number, with v 0 = k B T 0 m i the ionic speed of sound , ξ = ϑ 0 v 0
the ratio of the typical velocity to the ionic speed of sound , α = ϑ 0 c the ratio of the typical velocity to the speed of light , η = qx 0 E 0 k B T 0 the ratio of the electric and plasma internal energies , β = ϑ 0 B 0 E 0 the induced electric field relative to the total electric field , κ -1 = ν ii,0 t 0 the number of ion-ion collisions during the typical time .

The dimensionless Maxwellians are defined by

M ne,uα,Tα = M ne,uα,Tα = n e (x, t) 1 2πT α (x, t) Dv 2 exp - |Mεu α (x, t) -v| 2 2T α (x, t) , (20a)
M n i ,uα,Tα = M n i ,uα,Tα = n i (x, t) 1 2πT α (x, t) Dv 2 exp - |Mu α (x, t) -v| 2 2T α (x, t) . (20b)
Some comments can be stated regarding the meaning of these parameters and the scaling relations.

The typical mean velocity and temperature are assumed to be the same for the electrons and the ions. Accordingly, the relaxation of the electron mean velocity and temperature towards that of the ions may be assumed to marginally contribute to the evolution of the system. This assumption is therefore consistent with the investigation of resistive-less plasma modellings and the neglect of the inter-species collisions.

The parameter ξ is intended to provide a measure of how the electronic and ionic dynamics are resolved. The choice ξ = 1 means that the system is assumed to evolve at a speed comparable to the ionic thermal velocity v 0 , while ξε = 1 performs a rescaling of this typical velocity to the electron microscopic velocity. Setting ξ = M relates the typical speed of the system to the ionic mean velocity u 0 . Actually, the Mach number measures the gap between the microscopic (thermal) and macroscopic velocity scales.

The scaling relation η = 1 is generally assumed in single fluid plasma representation. The plasma internal energy is then on a par with the electric energy. This equilibrium is fundamental in the derivation of the Boltzmann relation. The identity βM = ξ is also common in single-fluid plasma models. This amounts to assume that the induced electric field scales as the product of the plasma mean velocity and the typical magnetic field: E 0 = u 0 B 0 . In other words, the magnetic field is essentially transported with the plasma flow. This later assumption is in line with the Alfven's frozen theorem [START_REF] Moreau | Magnetohydrodynamics[END_REF][START_REF] Freidberg | Ideal MHD[END_REF][START_REF] Davidson | An introduction to magnetohydrodynamics[END_REF][START_REF] Schnack | Lectures in Magnetohydrodynamics[END_REF] characteristic of ideal MHD models: the magnetic field is frozen into the plasma and transported by its flow.

The derivation of reduced models consists in identifying small dimensionless parameters and let them go to zero. The smallness of the scaled Debye length refers to a typical space scale much larger than the physical Debye length. This means that the charge separations, occurring on space scales comparable to the Debye length, are assumed unimportant to explain the evolution of the system. Sending the scaled Debye length to zero performs a low frequency filtering into the equation deriving thus a quasi-neutral model. In the context of the derivation of numerical methods, the typical length relates to the mesh size. This outlines the advantage of reduced models: the low frequency filtering operated by vanishing small parameters permits to derive numerical methods with discretization parameters (mesh size and time step) unconstrained by the small scales filtered out from the original equations.

3 A hierarchy of quasi-neutral models bridging the Vlasov-Maxwell system and the Hall-MHD regime

3.1 Handling the fluid and quasi-neutral limits

A hierarchy of fluid and kinetic quasi-neutral models

The aim here is to reduce the number of free dimensionless parameters, deriving by this means different reduced models well suited for the description of low frequency phenomena. As depicted in Fig. 1, the starting point of this hierarchy of models implements the minimal upgrades of the Vlasov-Maxwell system to recover a MHD regime. Precisely, only the inter-species collisions are taken into account in the initial model in order for the distribution function to relax towards the local equilibrium. This yields

ξ∂ t f i + v • ∇ x f i + η(E + β ξ v × B) • ∇ v f i = ξ κ ν ii (M n i ,u i ,T i -f i ) , (21a)
ξε∂ t f e + v • ∇ x f e -η(E + β εξ v × B) • ∇ v f e = ξ κ ν ee (M ne,ue,Te -f e ) , (21b)
for the evolution of the ions and electrons coupled to the dimensionless Maxwell system defined by Eqs. [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF].

From the scaling relations stated by Eq. ( 15), discarding the inter-species collisions make sense for the ions. Due to their large mass, the ions are almost unaffected by encounters with electrons. For the electrons, this assumption is not in line with the scaling of the like and inter-species collision frequencies. However, the purpose here is to propose a physically meaningful framework to clarify the foundation of a numerical method bridging the gap between a kinetic description of a weakly (or non) collisional magnetized plasma with a MHD regime. The interspecies collisions give rise to the resistivity in the macroscopic system which is not the targeted class of modelling for this work.

Handling the fluid limit

To identify easily a fluid regime, the distribution function is decomposed into a Maxwellian M nα,uα,Tα and a deviation from this Maxwellian κ g α according to

f α = M nα,uα,Tα + κ g α , (22a)
  1 v |v| 2   g α = Ωv   1 v |v| 2   g α dv = 0 . ( 22b
)
With this decomposition, the Vlasov-Boltzmann equations ( 21) can be recast into a hydrodynamic set of equations with kinetic corrections, depending on the moment of the deviation g α , yielding

ξ M ∂ t n i + ∇ x • (n i u i ) = 0 (23a) M 2 ξ M ∂ t (n i u i ) + ∇ x • (n i u i ⊗ u i ) + ∇ x p i -η n i E + βM ξ u i × B = -κ∇ x • v ⊗ vg i , (23b) ξ M ∂ t W i + ∇ x • (W i + p i )u i -ηn i E • u i = - κ M ∇ x • |v| 2 2 vg i , (23c)
with

W i = (M) 2 1 2 n i |u i | 2 + p i γ -1 , p i = n i T i , (23d)
for the ions, and an equivalent system for the electrons,

ξ M ∂ t n e + ∇ x • (n e u e ) = 0 (24a) (Mε) 2 ξ M ∂ t (n e u e ) + ∇ x • (n e u e ⊗ u e ) + ∇ x p e + η n e E + βM ξ u e × B = -κ∇ x • v ⊗ vg e , (24b) ξ M ∂ t W e + ∇ x • (W e + p e )u e + ηn e E • u e = - κ Mε ∇ x • |v| 2 2 vg e , (24c)
with

W e = (Mε) 2 1 2 n e |u e | 2 + p e γ -1 , p e = n e T e . (24d)
These two systems are coupled to a set of equations (the Maxwell system ( 18)) driving the changes in the electromagnetic field, as well as an equation for the evolution of the deviations g e and g i . The construction of these equations will be omitted in the present work, we refer for instance to [START_REF] Crouseilles | An asymptotic preserving scheme based on a micro-macro decomposition for Collisional Vlasov equations: diffusion and high-field scaling limits[END_REF][START_REF] Crestetto | Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles[END_REF][START_REF] De Cecco | Asymptotic Transition from Kinetic to Adiabatic Electrons along Magnetic Field Lines[END_REF] for examples of their derivation.

On the quasi-neutral limit

Omitting the collisions, the fastest velocity in this system is the propagation of waves at the speed of light described by the Maxwell system. The Debye length as well as the plasma period also define small space and time scales for large plasma densities. The quasi-neutral limit is defined by the following scaling relations:

(α, λ) → 0 , α ∼ λ . (25)
This amounts to assume that the scaled Debye length is small compared to the typical length and that the system evolves at a speed lower than the speed of light. By this means, the small scales related to these parameters are filtered out of the equations. The last hypothesis α ∼ λ is essential to recover the low frequency Ampere's law, derived by neglecting the displacement current. This equation being common to Magneto-Hydro-Dynamic models, the quasi-neutral limit encompasses these two assumptions. With the vanishing of this generalised dimensionless Debye length (λ, α) → 0, the Maxwell system degenerates into

β∇ x × B = M ξ J , ( 26a
)
β∂ t B + ∇ x × E = 0 , (26b)
n i = n e , (26c)
∇ x • B = 0 . ( 26d
)
From Gauss's law, the property of the electronic density to match that of the ions is recovered, which genuinely enforces quasi-neutrality of the plasma.

The electric field has no contribution in both these degenerate Gauss and the Ampere equations. The remaining occurrence of the electric field is limited to the Faraday equation (26b). Therefore, this set of equations is not well suited for the computation of the electric field. Indeed, the electrostatic component of the electric field can be arbitrarily chosen in Eqs. ( 26): the electric field satisfying this system may be augmented by any gradient of a scalar potential (see also [START_REF] Degond | Asymptotic-Preserving methods and multiscale models for plasma physics[END_REF] for further details). In the quasi-neutral limit, the electric field is provided by the particle current J rather than the displacement current (∂E/∂t originally present in Ampere's law). To close the system, the dependence of J with respect to E shall be explained to restore uniqueness of the electric field. This is related to the model describing the plasma.

A hierarchy of kinetic models for quasi-neutral plasmas 3.2.1 A kinetic formulation of the Electron-MHD

The aim here is to follow the microscopic dynamic of the electrons. The velocity of interest is the kinetic velocity of the electrons. This amounts to set ϑ 0 = v 0 /ε or equivalently ξε = 1, yielding

∂ t f i + ε (v • ∇ x f i + η(E + εβv × B) • ∇ v f i ) = 1 κ ν ii (M n i ,u i ,T i -f i ) , ( 27a
)
∂ t f e + v • ∇ x f e -η(E + βv × B) • ∇ v f e = 1 εκ ν ee (M ne,ue,Te -f e ) , (27b)
The collisions are assumed to be ineffective on the characteristic time scale: κε ≫ 1 which amounts to neglect the collision operator in Eqs. [START_REF] Gordeev | Electron magnetohydrodynamics[END_REF], in particular for the ions, owing to ε ≪ 1.

Performing the quasi-neutral limit (λ = α) → 0, the system at hand here is recast into (see [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF])

∂ t f i + ε (v • ∇ x f i + η(E + εβv × B) • ∇ v f i ) = 0 , ( 28
) ∂ t f e + v • ∇ x f e -η(E + βv × B) • ∇ v f e = 0 , (29)
β∇ x × B = (Mε)J , ( 30
) β∂ t B + ∇ x × E = 0 , (31)
n e = n i = n , (32)
∇ x • B = 0 .
First, note that the formal time derivative of the Faraday equation (31) together with the curl of Ampere's law yields

∇ x × ∇ x × E = -(Mε)∂ t J , (33)
which outlines that the electric field is known up to the gradient of a potential in this system. In [START_REF] Degond | Numerical approximation of the Euler-Maxwell model in the quasineutral limit[END_REF][START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF] the ill-posed nature of this equation is corrected by explaining the relation between the current density and the electric field. The conservation of the ionic and electronic momentum, as stated by Eqs. (23b) and (24b), yields

(Mε) ∂J ∂t = -(Mε) 2 ∇ x • (S i -S e ) -∇ x • ε 2 P i -P e + η(ε 2 n i + n e )E + ηβ(Mε) ε 2 n i u i + n e u e × B , (34a)
with, for any species α (α = e for electrons and i for the ions),

S α = n α u α ⊗ u α , P α = p α Id + κ v ⊗ vg α . (34b)
Inserting the identity (34a), together with n = n i = n e , into Eq. [START_REF] Jin | Efficient asymptotic-preserving (AP) schemes for some multiscale kinetic equations[END_REF] gives

n(1 + ε 2 )E + ∇ x × ∇ x × E = -(Mε)n ε 2 u i + u e × B + (Mε) 2 ∇ x • (S i -S e ) + ∇ x • ε 2 P i -P e . ( 35
)
This equation is well posed in the quasi-neutral limit (n > 0) and can be used for the computation of the electric field. It is written under the assumption β = 1 as well as η = 1 which amounts to consider the thermal energy on a par with the electric one. This yields the following definition of the quasi-neutral model

∂ t f i + ε (v • ∇ x f i + (E + εv × B) • ∇ v f i ) = 0 , (36a) ∂ t f e + v • ∇ x f e -(E + v × B) • ∇ v f e = 0 , (36b)
n(1 + ε 2 )E + ∇ x × ∇ x × E = -(Mε)n ε 2 u i + u e × B + (Mε) 2 ∇ x • (S i -S e ) + ∇ x • ε 2 P i -P e , (36c)
∂ t B + ∇ x × E = 0 , (36d) ∇ x • B = 0 .
Note that the electric field provided by Eq. (36c) enforces a divergence free particle current, or more precisely ∂ t (∇ • J) = 0. This yields, thanks to the continuity equation:

∂ 2 ρ ∂t 2 = 0 . This proves the consistency of this model with the quasi-neutrality assumption (matching of the electronic and ionic densities) as soon as the initial data are compliant with this regime. Note that, working a time semidiscretization, the vanishing of the charge density, rather than its time double derivative, may be recovered (see [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF]). Note also that the evolution of the ions only brings a marginal correction, proportional to ε 2 , to the electron dynamics. On this time scale, the ions may be considered at rest.

The characteristics of this model are similar to the so-called Electron MHD: the time scale of interest is that of the electrons, the ions merely creating a motionless background for the fast electron flows [START_REF] Kingsep | Electron Magnetohydrodynamics[END_REF]. In particular, this modelling accounts for the inertia of electrons. A noticeable difference with the Electron-MHD (see Sec. 3.4) lies in the kinetic description of the plasma. An Asymptotic-Preserving method is proposed in [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF] to bridge this quasi-neutral model and the Vlasov-Maxwell system. The properties of this quasi-neutral plasma description are investigated in [START_REF] Tronci | Neutral Vlasov kinetic theory of magnetized plasmas[END_REF] by means of a linear stability analysis.

A hybrid formulation of the Hall-MHD

Hybrid modelling [START_REF] Winske | Hybrid Simulation Codes: Past, Present and FutureA Tutorial[END_REF][START_REF] Yin | Hybrid and Hallmagnetohydrodynamics simulations of collisionless reconnection: Effect of the ion pressure tensor and particle Hall-magnetohydrodynamics[END_REF][START_REF] Tronci | Hybrid Vlasov-MHD models: Hamiltonian vs. non-Hamiltonian[END_REF] refers to a class of plasma models where the ions are described by a kinetic equation while the fluid limit is assumed for the electrons. This is in line with the scaling relations of the collision frequencies stated by Eqs. [START_REF] Degond | Chapter One -Asymptotic continuum models for plasmas and disparate mass gaseous binary mixtures[END_REF]. The relaxation of the electronic distribution function towards the local equilibrium is indeed faster than for the ions. The aim of this modelling is to filter out of the equations the fast scales carried by the electron dynamics. Therefore, a zero inertia regime is also assumed for the electrons together with the fluid limit and the quasi-neutrality of the plasma.

The typical velocity selected here is the microscopic (thermal) velocity of the ions. This translates into the identity ξ = 1 resulting in the following system for the plasma:

∂ t f i + v • ∇ x f i + η(E + β ξ v × B) • ∇ v f i = 1 κ ν ii (M n i ,u i ,T i -f i ) , (37a)
∂ t f e + 1 ε v • ∇ x f e -η(E + β εξ v × B) • ∇ v f e = 1 εκ ν ee (M ne,ue,Te -f e ) . (37b)
The fluid limit for the electrons is selected assuming (εκ) ≪ 1 meaning that the number of electron collisions during the typical time is large. The quasi-neutrality of the plasma amounts to set λ = α ≪ 1. To overcome the degeneracy of the Maxwell system in the quasi-neutral limit, the electronic momentum is harnessed to provide the so-called generalised Ohm's law. The electronic system can be recast into (Mε) 2 where the viscous stress tensor σ e as well as the thermal conductivity µ e are defined in Sec. A.3 and, owing to the quasi-neutrality assumption, n e = n i = n.

The dynamic described by these equations is stiff, this is due to the smallness of (ξε) in this regime: the thermal velocity of the ions (defined as the typical velocity) is small compared to that of the electrons. Therefore, the electrons are in a low Mach regime. Assuming (Mε) ≪ 1 gives rise to the following equilibria

∇ x (nT e ) + η n E + (βM)u e × B = 0 , (38a)
∇ x T e = 0 . (38b)
The classical massless approximation for the electrons is recovered with the generalised Ohm's law and a homogeneous electronic temperature. The definition of the mean velocity u e is derived from the particle current density J = n(u i -u e ) together with Ampère's law [START_REF] Hagelaar | How to normalize MaxwellBoltzmann electrons in transient plasma models[END_REF], yielding

u e = u i - β M ∇ x × B n . ( 39
)
The hybrid plasma modelling writes (assuming η = β = 1)

∂ t f i + v • ∇ x f i + (E + v × B) • ∇ v f i = 0 , (40a)
E = -Mu × B + M ∇ x × B n × B -T e ∇ x n n , (40b)
∂B ∂t + ∇ x × E = 0 , ∇ x • B = 0 , (40c)
with

n = f i dv , nu = vf i dv . (40d)
The derivation of a similar model is proposed in [START_REF] Acheritogaray | Kinetic formulation and global existence for the hall-magneto-hydrodynamics system[END_REF] with numerical investigations in [START_REF] Degond | Simulation of non equilibrium plasmas with a numerical noise reduced particle in cell method[END_REF].

A fluid hierarchy of quasi-neutral models

The Electron MHD-system

This model is obtain by letting κ → 0 in Eqs. [START_REF] Degond | Numerical approximation of the Euler-Maxwell model in the quasineutral limit[END_REF]. This yields the following set of equations for the electrons

ξ M ∂ t n e + ∇ x • (n e u e ) = 0 (41a) ξ M ∂ t (n e u e ) + ∇ x • (n e u e ⊗ u e ) + 1 (Mε) 2 ∇ x p e + η n e (E + βM ξ u e × B) = 0 , (41b) ξ M ∂ t W e + ∇ x • (W e + p e )u e + ηn e E • u e = 0 . (41c)
A similar system is derived for the ions however with ε = 1 and η replaced by -η. These two sets of conservation laws are coupled to the Maxwell system [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF]. Performing the quasi-neutral limit in this system (λ = α) → 0 and focusing on the electronic dynamics with ξε = 1 yields the quasi-neutral bi-fluid Euler-Maxwell system. This model is similar to that of Sec. 3.2.1 but with a fluid description for the plasma. Following the same terminology, this model is referred to as Electron-MHD system (also used in [START_REF] Kingsep | Electron Magnetohydrodynamics[END_REF] though different terminology may be used by other authors). It is implemented and investigated numerically in the framework of Asymptotic-Preserving methods in [START_REF] Degond | Numerical approximation of the Euler-Maxwell model in the quasineutral limit[END_REF].

The Hall-MHD regime

The Hall-MHD regime (see [START_REF] Lighthill | Studies on Magneto-Hydrodynamic Waves and other Anisotropic Wave Motions[END_REF][START_REF] Witalis | Hall Magnetohydrodynamics and Its Applications to Laboratory and Cosmic Plasma[END_REF][START_REF] Schnack | Lectures in Magnetohydrodynamics[END_REF]) is recovered from the assumptions of the precedent section but with a typical velocity equal to the plasma mean flow yielding ξ = M. The fast electronic dynamics is filtered out from the equations to provide a low frequency modelling for the plasma driven by the evolution of the massive ions. The plasma velocity, denoted u is defined as that of the heavy species u = u i . The other parameters obey the classical scaling relations of MHD models: η = 1 and β = 1.

In the drift regime (Mε → 0), the electronic energy reduces to the internal energy E e = p e /(γ -1) ,

with the electronic momentum and energy verifying

∇ x p e = -n (E + u e × B) , (42b)
∂ t E e + ∇ x • (E e + p e )u e = -nE • u e . ( 42c
)
The generalised Ohm's law (42b) is harnessed to compute the electric field. The electronic velocity u e is substituted by u e = u -J/n. The electric field is computed thanks to the generalized Ohm's law, giving rise to

E = -u × B + ∇ x × B n × B - ∇ x p e n . ( 43
)
The first term of the right hand side of the equation ( 43) is the classical frozen field term, explaining the convection of the magnetic field together with the plasma . The second and third contributions are the so-called Hall and diamagnetic terms. Inserting this definition in the Faraday law (26b), the magnetic induction equation can be constructed, with

∂ t B + ∇ x • (u ⊗ B -B ⊗ u) = -∇ x × ∇ x × B n × B + ∇ x × ∇ x p e n . (44)
Finally the plasma mass density, momentum and total pressure p and energy W defined by

p = p i + p e = n(T i + T e ) , (45a)
W = W i + E e = (M) 2 1 2 n|u| 2 + p γ -1 , (45b)
verify

∂ t n + ∇ x • (nu) = 0 (45c) ∂ t (nu) + ∇ x • (nu ⊗ u) + 1 M 2 ∇ x p = 1 M 2 J × B (45d)
∂ t W + ∇ x • (W + p)u + ∇ x • (E e + p e )v H = E • J , (45e)
the Hall velocity v H , which can be interpreted as the electron velocity in the ion frame, is defined by

v H = - J n . (46)
The ideal MHD equations are classically written under a conservative form using the system total pressure and energy

p TOT = p + B 2 2 , W TOT = W + B 2 2 , (47a)
writing the system as

∂ t n + ∇ x • (nu) = 0 (47b) ∂ t (nu) + ∇ x • nu ⊗ u - 1 M 2 B ⊗ B + 1 M 2 ∇ x p TOT = 0 (47c) ∂ t W TOT + ∇ x • (W TOT + p TOT )u -(B • u)B = -∇ x • (E e + p e )v H , (47d)
∂ t B -∇ x • B ⊗ (u + v H ) -(u + v H ) ⊗ B = ∇ x × ∇ x p e n . (47e)
This set of equations is supplemented with the electronic energy conservation (42c).

The ideal MHD equations are recovered from this system assuming an ideal Ohm's law where the current density is assumed small compared to the ion mean velocity and therefore neglected. However, in this simplified framework (omitting the unlike particle collisions), there are no mechanisms preventing the electron mean velocity to depart from that of the ions. Consequently, the generalized Ohm's law incorporates the Hall velocity in complement to the so called ideal Ohm's law. The effect of the resistivity should be consider to derive the ideal MHD regime.

The drift approximation operated for the electrons amounts to vanishing the electronic Mach numbers (Mε). The scale separation introduced by the small electron to ion mass ratio ε is not always sufficient to consider this limit independently to vanishing ionic Mach numbers M → 0. For low ionic Mach number a low frequency filtering may be operated performing the limit of vanishing electronic Mach numbers jointly with the ionic Mach numbers. This asymptotic defines the Massless MHD regime [START_REF] Besse | A model hierarchy for ionospheric plasma modeling[END_REF]. [START_REF] Chen | An energy-and chargeconserving, implicit, electrostatic particle-in-cell algorithm[END_REF] The electrostatic regime and the Boltzmann relation

Electrostatic limit of the Maxwell system

The electrostatic regime is recovered from the dimensionless Maxwell system (18) by letting α go to zero, independently of λ. This assumption shall be interpreted as a typical velocity negligible compared to the speed of light. From Ampère's law

λ 2 η ∂E ∂t + M ξ J = βλ 2 η α 2 ∇ x × B , (48)
the limit α → 0 provides ∇ × B = 0. Nonetheless, the right hand side of Eq. ( 48) remains an undetermined form. Therefore the Ampere's law is not well suited for the computation of the electric field in the electrostatic limit. However, subjected to convenient boundary conditions, the property

∇ x × B = 0 together with ∇ x • B = 0 yields ∂ t B = 0.
Inserting this property into the Faraday equation (18b) provides an electrostatic electric field: ∇ × E = 0. Furthermore, the undetermined form in Eq. ( 48) is divergence free. Therefore, computing the divergence of Ampere's law provides

λ 2 η ∂∇ x • E ∂t = - M ξ ∇ x • J ,
which is a well posed problem for the electric field under the condition ∇ x × E = 0. Note that, owing to the continuity equation

∂ρ ∂t + M ξ ∇ x • J = 0 ,
originating from the conservation of the particle densities (23a) and (24a), the divergence of Ampère's law is equivalent to the time derivative of Gauss law, with

∂ ∂t λ 2 η∇ x • E -ρ = 0 .
Therefore, in the electrostatic regime, the Gauss equation is used to compute the electric field.

Quasi-neutral models at the electronic scale

This analysis is carried out under the assumption of a vanishing magnetic field (B = 0). The plasma description considered in the sequel is therefore

ξ∂ t f i + v • ∇ x f i + ηE • ∇ v f i = ξ κ ν ii (M n i ,u i ,T i -f i ) , (49a)
ξε∂ t f e + v • ∇ x f e -ηE • ∇ v f e = ξ κ ν ee (M ne,ue,Te -f e ) (49b)
coupled to the Gauss equation

-λ 2 η∆φ = n i -n e , (49c)
φ being the electrostatic potential, with E = -∇φ and n α = f α dv. The quasi-neutrality of the plasma is recovered for vanishing scaled Debye length λ → 0. In this regime, similarly to the electromagnetic case, an equation needs to be manufactured from the motion of the particles, to compute the electric field. This is classically obtained thanks to the equation of the electronic momentum conservation. The electric field is computed in order for this conservation to be satisfied. In [START_REF] Degond | Asymptotic-preserving particle-in-cell method for the Vlasov-Poisson system near quasineutrality[END_REF] an equivalent approach is proposed. It consists in using the time derivatives of Gauss law to produce

-λ 2 η ∂ 2 ∂t 2 ∆φ = ∂ 2 ∂t 2 (n i -n e ) . (50)
In the quasi-neutral limit (λ → 0), the electrostatic field is the Lagrangian multiplier of the quasi-neutrality constraint ∂ 2 ∂t 2 (n i -n e ) = 0 .

From the system (24) (with B = 0), formally computing the time derivative of the density conservation together with the divergence of the momentum equation, the following identity is recovered

ξ M 2 ∂ 2 n e ∂t 2 = -∇ x • ∇ x • (n e u e ⊗ u e ) - 1 (Mε) 2 ∇ x • ∇ x P e + ηn e E ; (51)
P e = p e Id + κ v ⊗ vg e , with v ⊗ vg e = (v ⊗ v)g e dv . (52)
Resuming the scaling relations of Sec. 3.2.1: ξ = 1/ε and κε > 1, assuming that the ions are at rest, the evolution of the charge density n i -n e is governed by Eq. ( 51) with

∂ 2 n e ∂t 2 = -(Mε) 2 ∇ x • ∇ x • (n e u e ⊗ u e ) -∇ x • ∇ x P e + ηn e E . ( 53
)
The evolution of the density is barely independent of the mean flow velocity but relies on the balance between the pressure and the electric forces.

The equation providing the electric potential φ is obtained by inserting this relation into Eq. ( 50) and passing to the limit λ → 0. This yields the following quasi-neutral kinetic plasma description

∂ t f e + v • ∇ x f e + η∇ x φ • ∇ v f e = 1 εκ ν ee (M ne,ue,Te -f e ) , (54a)
-η∇ x • (n e ∇ x φ) = -∇ x • ∇ x • P e -(Mε) 2 ∇ x • (n e u e ⊗ u e ) . (54b)
According the values of κ, the collision term in Eq. (54a) may be disregarded, defining therefore a non collisional kinetic description. Contrariwise, letting κ → 0, a fluid description for the electrons may be derived, with

Note that a similar equation to ( 53), but with P e = p e Id, may be worked out of the conservation of the electronic density (55a) and momentum (55b). This outlines that the electronic dynamics, in particular the electronic speed of sound, is resolved in this model. Comparable models are implemented and numerically experienced in the context of Asymptotic-Preserving methods in [START_REF] Crispel | An asymptotic preserving scheme for the two-fluid Euler-Poisson model in the quasineutral limit[END_REF] for fluid plasma description and [START_REF] Manfredi | Vlasov modelling of parallel transport in a tokamak scrape-off layer[END_REF] for kinetic equations.

Quasi-neutral models at the ionic scale

The typical velocity is chosen to be that of the ions with either ξ = 1 for the kinetic descriptions of the ions and ξ = M for macroscopic models. The hybrid modelling investigated in Sec. 3.2.2 is defined by the scaling relations ξ = 1, (εκ) ≪ 1 and η = 1. The equilibria stated by Eqs. [START_REF] Lighthill | Studies on Magneto-Hydrodynamic Waves and other Anisotropic Wave Motions[END_REF] yields:

T e ∇ x n e = n e ∇ x φ ,

with an homogeneous electronic temperature. This equation is integrated to provide the so-called Boltzmann relation

n e = n 0 exp - φ T e , (57)
n 0 being a constant (independent of the space variable x) that should be determined from adequate conditions [START_REF] Hagelaar | How to normalize MaxwellBoltzmann electrons in transient plasma models[END_REF]. Due to the Boltzmann relation, the quasi-neutral limit is not singular any more. Indeed plugging the Boltzmann relation (57) into the Gauss equation yields

-λ 2 ∆φ = n i -n 0 exp - φ T e . ( 58
)
This equation is not degenerate for the computation of the electric potential for vanishing λ. Indeed, the non linear part of the equation provides a means of computing φ in the quasi-neutral limit. This property is thoroughly investigated in [START_REF] Degond | Numerical Approximation of the Euler-Poisson-Boltzmann Model in the Quasineutral Limit[END_REF].

The hybrid electrostatic model may be recast into

∂ t f i + v • ∇ x f i -∇ x φ • ∇ v f i = 1 κ ν ii (M n i ,u i ,T i -f i ) , (59a)
φ = -T e ln n i n 0 . ( 59b
)
Letting κ → 0 together with κ/(Mε) ≪ 1 and ξ = M yields the quasineutral fluid model

∂ t n i + ∇ x • (n i u i ) = 0 (60a) ∂ t (n i u i ) + ∇ x • (n i u i ⊗ u i ) + 1 M 2 (∇ x p i + n i ∇ x φ) = 0 , (60b)
∂ t W i + ∇ x • (W i + p i )u i + n i ∇ x φ • u i = 0 , (60c) φ
= -T e ln n i n 0 . (60d)
In the models (59) and (60) following the evolution of the plasma at the ionic scale, the fast electronic dynamics introduced by the electron inertia is filtered out of the equations by performing the low frequency limit (Mε) → 0.

Conclusions

In this paper, we propose an asymptotic analysis bridging kinetic plasma descriptions coupled to the Maxwell system and single plasma modelling. Two frameworks are investigated. The first one is devoted to electromagnetic fields. The plasma is represented by a hierarchy of models starting with the bi-kinetic Vlasov-Maxwell system while ending with the single fluid Hall-Magneto-Hydro-Dynamic model. The second framework is dedicated to electrostatic fields. In this context, the asymptotic analysis permits to derive a hierarchy of models bridging the bi-kinetic Vlasov-Poisson system to a single fluid representation consisting of a fluid system for the ions coupled to the Boltzmann relation for the electrons. The investigations proposed within this document unravel different asymptotic parameters explaining the transition from one model to the other. The effort conducted in the present work consists in relating these asymptotic parameters to characteristics of the system. This means that the transition from one model to the other may be explained by a change in the plasma characteristics or the typical scales at which the plasma is observed. This last notion is important in the perspective of designing a numerical method. Indeed, the discretization of these equations requires the use of a mesh interval as well as a time step. These two numerical parameters define the typical space and time scale, therefore a velocity scale as well, the numerical method is aimed at capturing. This is related to the parameter ξ used for the asymptotic analysis. Regarding the quasi-neutral modellings investigated within this document, different choices are operated for this parameter. The fastest scales are related to the electron thermal velocity when the fast electron dynamics is intended to be captured by the model. This is for instance the value selected for Electron-MHD regimes, either in the fluid or kinetic frameworks. For hybrid or single fluid plasma representations, the velocity scale is reduced to that of the mean flow of the plasma defined by the massive ions. The organisation of Sec. 4 is aimed at emphasizing this feature.

The second parameter, already established in precedent works (see [START_REF] Degond | Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit[END_REF][START_REF] Degond | Numerical approximation of the Euler-Maxwell model in the quasineutral limit[END_REF][START_REF] Degond | Asymptotic-preserving particle-in-cell method for the Vlasov-Poisson system near quasineutrality[END_REF]), is the generalized scaled Debye length λ. It actually encompasses the scaled Debye length and the ratio of the typical velocity to the speed of light. Vanishing the generalised Debye length amounts to filter out from the equations the small scales attached to the charge separation as well as those related to the propagation of electromagnetic waves at the speed of light. The quasi-neutral limit is therefore a low frequency limit. Quasi-neutrality breakdowns may be explained by a refinement of the typical length scale or, for instance, a decrease of the plasma density. This changes are well accounted for by the asymptotic parameters set up to perform the analysis.

The vanishing of the electron inertia is related to a low Mach regime (Mε) ≪ 1. In single fluid plasma representation, the fast electron dynamic is dropped out of the equations to perform a low frequency filtering, the system being assumed to evolve at a lower speed attached to the massive ions. Nonetheless, the nature of the flow may be subjected to significant changes explaining that the particle inertia becomes significant again to account for the system evolution. This is illustrated in studies of plasma flows in sheaths, with supersonic particles, while the mean plasma velocity is small compared to the speed of sound in the plasma bulk [START_REF] Chodura | Plasma Flow in the Sheath and the Presheath of a Scrape-Off Layer[END_REF][START_REF] Stangeby | The plasma boundary of magnetic fusion devices[END_REF][START_REF] Grismayer | Electron kinetic effects in plasma expansion and ion acceleration[END_REF][START_REF] Drouin | Particle-in-cell modeling of relativistic laserplasma interaction with the adjustable-damping, direct implicit method[END_REF][START_REF] Manfredi | Plasma-wall transition in weakly collisional plasmas[END_REF]. Accounting for this phenomena is possible selecting the appropriate typical velocity to resolve or filter the fast electron dynamic.

Finally the fluid assumption is classically related to a vanishing of the Knudsen number κ ≪ 1 accounting for the relaxation of the distribution function towards the local thermodynamic equilibrium. The interplay of these four dimensionless parameters (ξ, λ, Mε, κ) defines a hierarchy of reduced models bridging kinetic plasma descriptions coupled to the Maxwell system to quasi-neutral plasma representations including kinetic, hybrid and single fluid modellings. Other reduced models may also be investigated considering different closures, to give access to a class of models with non scalar pressure, or fluid models evolving higher order moments. These low frequency models are widely used to design efficient numerical methods. The benefit of the asymptotic analysis conducted within this document draws the guidelines for the derivation of numerical methods implementing local up-scalings, therefore widening the use of numerical methods discretizing these reduced models.

where ϕ = ϕ dv. For k = 1, . . . , D v , we have the following properties

(Id -Π Me )(∂ t M e ) = (Id -Π Me )(E • ∇ v M e ) = 0 , (62a)

Figure 1 :

 1 Figure 1: Fluid and kinetic (quasi-neutral) model hierarchies derived from the Vlasov-BGK-Maxwell system.

  ∂ t n e + ∇ x • (n e u e ) ∂ t (n e u e ) + ∇ x • (n e u e ⊗ u e ) + 1 (Mε) 2 (∇ x p e -ηn e ∇ x φ) = 0 , (55b) 1 (Mε) ∂ t W e + ∇ x • (W e + p e)u e -ηn e ∇ x φ • u e = 0 .

( 2 2T 2 e- D v + 2 2 2T 2 e-A. 3 2 2T 2 e+

 22222322 Id -Π Me ) 1 n e M e v k ∂ x k n e = 0 ; (62b) (Id -Π Me ) v -(Mε)u e • ∂ x k u e v k M e = -v -(Mε)u e 2 D v ∂ x k u e,k+ v k -(Mε)u e,k v -(Mε)u k • ∂ x k u e M e ;(62c)(Id -Π Me ) |v -(Mε)u e | 2 2T 2 e -D v 2T e M e v k ∂ x k T e = M e v -(Mε)u e 2T e v k -(Mε)u e,k ∂ x k T e .(62d)Furthermore, if g e satisfies the Micro-Macro decomposition[START_REF] Degond | The moment-guided Monte Carlo method[END_REF], the following identities holds true Π Me (g e ) = Π Me (∂ t g e ) = 0 (62e)A.2 Computation of the deviation to the MaxwellianThe aim here is to characterize g e or more specifically an approximation to the first order in κ. Inserting the micro-macro decomposition into the Vlasov-Boltzmann equation (49b) yieldsξε∂ t M e +v•∇ x M e -ηE •∇ v M e +κ (ξε∂ t g e + v • ∇ x g e -ηE • ∇ v g e ) = L Me g e ,where L Me g e = -ξν ee g e .This provides, using properties (62),L Me g e = (Id-Π Me ) (v • ∇ x M e )+κ ξε∂ t g e +(Id-Π Me ) v•∇ x g e -ηE•∇ v g e .It followsg e = (L Me ) -1 (Id -Π Me ) (v • ∇ x M e ) + O(κ) ,(63a)withv k ∂ x k M e = ∂ x k n e n e + Mε T e v-(Mε)u e •∂ x k u e + (v -(Mε)u e ) 2 2T 2 e -D v 2T e ∂ x k T e v k M e .(63b) From properties (62), we can state the expression of the deviation to the Maxwelliang e =x k u e,k + v k -(Mε)u e,k v -(Mε)u e • ∂ x k u e + v -(Mε)u e D v + 2 2T e v k -(Mε)u e,k ∂ x k T e + O κ ξ .Computation of the viscous termsThe viscous terms are defined byv ⊗ vg e = -2 σ e • u e + q e + O κ ξ .Following the characterization (64) of g e we can writeσ e =v)M e v l -(Mε)u e,l Dv k=1 v k -(Mε)u e,k ∂ x l u e,k , Mε)u e • ∇ x T e ,Inserting in these definitions, the following identities|v| 2 2 (v ⊗ v)M e v -(Mε)u e D v + 2 2T e v -(Mε)u e • ∇ x T e = 0 , v)M e v l -(Mε)u e,l Dv k=1 v k -(Mε)u e,k ∂ x l u e,k = (Mε)σ e • u e ; we obtain σ e = -1 ν ee (n e T e ) ∇ x u e + ∇ x u T e -2 D v (∇ x • u e ) Id , (65a) q e = -D v + 2 2ν ee (n e T e ) ∇ x T e . (65b)

  function ϕ, the projector onto the Maxwellian denoted M e with M e = n e (x, t)For any smooth function ϕ, the projector onto M e , denoted Π Me , is defined as (see[START_REF] Bennoune | Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier-Stokes asymptotics[END_REF][START_REF] Crestetto | Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles[END_REF])

			1 2πT e (x, t)	Dv/2 exp	|v -(Mε)u e (x, t)| 2 2T e (x, t)	.		(61a)
	Π Me (ϕ) = ϕ +	v -(Mε)u e T e	• v -(Mε)u e ϕ
	+	2 D v	|v -(Mε)u e | 2 2T e	-	D v 2	|v -(Mε)u e | 2 2T e	-	D v 2	ϕ	M e n e	.

Acknowledgments

This work has been supported by "Fédération de recherche pour la Fusion par Confinement Magnétique" (FrFCM) in the frame of the project "BRIDIPIC: BRIDging Particle-In-Cell methods and low frequency numerical models of plasmas". A. Crestetto acknowledges support from the french "Agence Nationale pour la Recherche (ANR)" in the frame of the projects MoHyCon ANR-17-CE40-0027-01 and MUFFIN ANR-19-CE46-0004.

A Micro-Macro decomposition, computation of the viscous terms.

A.1 Introduction, definitions, elementary properties

The analysis carried out in this section are developed in the electrostatic framework and specified for the electrons. The extensions for either electromagnetic fields or the ions are straightforward and are therefore omitted for conciseness.

We first introduce the projector onto the Maxwellian For any smooth