1 Introduction
The aim of this paper is to propose a continuation of the work initiated in Degond, Deluzet & Navoret (Reference Degond, Deluzet and Navoret2006), Degond et al. (Reference Degond, Deluzet, Navoret, Sun and Vignal2010) and Degond, Deluzet & Doyen (Reference Degond, Deluzet and Doyen2015) 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 those former achievements is aimed to be generalized here to more singular limits. In this series of prior 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 Degond et al. (Reference Degond, Deluzet and Doyen2015) can be interpreted as a kinetic description of the electron magnetohydrodynamics (MHD) (Gordeev, Kingsep & Rudakov Reference Gordeev, Kingsep and Rudakov1994; Swanekamp et al. Reference Swanekamp, Grossmann, Fruchtman, Oliver and Ottinger1996; Cho & Lazarian Reference Cho and Lazarian2004), 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 electron MHD with the aim of bridging the Vlasov–Maxwell system and 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 with 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 magnetized plasmas, or the Vlasov–Poisson system and the electron adiabatic response, also referred to as the Boltzmann relation (see Langmuir (Reference Langmuir1929) and Tonks & Langmuir (Reference Tonks and Langmuir1929a,Reference Tonks and Langmuirb) for seminal works and De Cecco et al. (Reference De Cecco, Deluzet, Negulescu and Possanner2017) 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 Crestetto, Crouseilles & Lemou (Reference Crestetto, Crouseilles and Lemou2012), Crouseilles & Lemou (Reference Crouseilles and Lemou2011), Dimarco, Mieussens & Rispoli (Reference Dimarco, Mieussens and Rispoli2014) and Lemou & Mieussens (Reference Lemou and Mieussens2008) for micro–macro methods, and to Degond, Dimarco & Pareschi (Reference Degond, Dimarco and Pareschi2011b) for the moment guided method, fluid-preconditioned fully implicit methods (Chen, Chacón & Barnes Reference Chen, Chacón and Barnes2011; Chen et al. Reference Chen, Chacón, Leibs, Knoll and Taitano2014; Chen & Chacón Reference Chen and Chacón2015) and asymptotic-preserving numerical methods (Jin Reference Jin1999; Degond & Deluzet Reference Degond and Deluzet2017). In these methods, macroscopic models are operated to either filter out the numerical noise of particle-in-cell methods (moment guided, micro–macro particle methods) or speed up the convergence of nonlinear implicit solvers (fully implicit particle-in-cell methods). The efficiency of these methods is closely related to that of the reduced models and adequacy of the kinetic plasma description. Another application can be envisioned with the hybrid coupling of particle-in-cell methods and MHD descriptions (Schumer et al. Reference Schumer, Swanekamp, Ottinger, Commisso, Weber, Smithe and Ludeking2001; Daldorff et al. Reference Daldorff, Töth, Gombosi, Lapenta, Amaya, Markidis and Brackbill2014) and more generally coupling strategies such as the current-coupling scheme and the pressure-coupling scheme (see Park et al. (Reference Park, Parker, Biglari, Chance, Chen, Cheng, Hahm, Lee, Kulsrud and Monticello1992), Tronci et al. (Reference Tronci, Tassi, Camporeale and Morrison2014) and 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 $\varepsilon _0$, the vacuum permittivity, go to zero (e.g. Jang & Masmoudi Reference Jang and Masmoudi2012; Tronci & Camporeale Reference Tronci and Camporeale2015) which is referred to as the full Maxwell or to the low-frequency pre-Maxwell asymptotic equations in Freidberg (Reference Freidberg2014, see § 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 (Freidberg Reference Freidberg2014; Klingenberg, Pirner & Puppo Reference Klingenberg, Pirner and Puppo2016) 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 a regime transition: the electron to ion mass ratio remains constant and the same property holds true for the vacuum permittivity.
The outline of the paper is the following. The plasma kinetic description is introduced in § 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 models is proposed in § 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 § 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 analyses is proposed in § 5 devoted to conclusions.
2 The Vlasov–Maxwell system in a dimensionless form
2.1 Objectives
In this section, the purpose is to unravel a series of asymptotic limits bridging the gap between the Vlasov–Maxwell system and a 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 achieve 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 Bhatnager–Gross–Krook (BGK) operator. This choice of collision operator is questionable from a strict modelling viewpoint; 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 intra- 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.
2.2 The Vlasov–BGK–Maxwell system
The most refined description of the plasma is constituted by two Vlasov equations, $f_\textrm {i}$ and $f_\textrm {e}$ being the ion and electron distribution functions:
In these equations, $q$ is the elementary charge and $m_\alpha$ is the mass of the species $\alpha$ ($\alpha =\textrm {e}$ for the electrons and $\textrm {i}$ for the ions). The BGK collision operator $\mathcal {Q}_\alpha$ is given by (Huba Reference Huba2011)
Here $\nu _{\alpha \alpha }$ and $\nu _{\alpha \beta }$ are the like-particle and inter-species collision frequencies which can be defined as (Degond Reference Degond2007; Spatschek Reference Spatschek2012)
where
$C$ denoting a constant with a magnitude equal to one, $\ln (\varLambda )$ the Coulomb logarithm and the ions being assumed mono-charged.
The Maxwellians $\mathcal {M}_{n_\alpha ,u_\alpha ,T_\alpha }$ and $\bar {\mathcal {M}}_{n_\alpha ,\bar {u}_\beta ,\bar {T}_\beta }$ are defined as
$D_v$ denoting the dimension of the velocity space and $k_\textrm {B}$ is the Boltzmann constant. The Maxwellian parameters $n_\alpha$, $u_\alpha$ and $T_\alpha$ are the density, mean velocity and temperature associated with the distribution function $f_\alpha$ and defined as
with $\gamma$ the specific heat ratio whose value depends on the dimensionality of the velocity space $D_v$ through
The collision operators verify the following conservation properties:
The temperature and the mean velocity $(\bar {u}_\beta ,\bar {T}_\beta )$ in the inter-species collision operator expression (2.5b) should be chosen with care in order to guarantee the total momentum and energy conservation. Indeed the identities
hold true for the operators defined by (2.5). The trivial choice $(\bar {u}_\beta ,\bar {T}_\beta )=({u}_\beta ,{T}_\beta )$ does ensure the plasma total momentum conservation, provided that $\nu _{\textrm {ei}} m_\textrm {e} n_\textrm {e} = \nu _{\textrm {ie}} m_\textrm {i} n_\textrm {i}$. However, in this case, the plasma total energy is not conserved. We refer to Greene (Reference Greene1973) for a seminal work, as well as to Klingenberg et al. (Reference Klingenberg, Pirner and Puppo2016) and references therein for recent advances, on the choice of these parameters compliant with the desired properties (2.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:
where $c$ is the speed of light, $\mu _0$ the vacuum permeability and $\epsilon _0$ the vacuum permittivity verifying $\mu _0\epsilon _0 c^2 = 1$. The Maxwell sources are the particle currents and densities
The definition of the collision frequencies as stated by (2.4) relates to different time scales. Indeed, because of their different masses, ions and electrons are not equally affected by collisions. These properties are more clearly emphasized working with dimensionless variables as proposed in the next section.
2.3 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:
We denote by $x_0$ and $t_0$ the characteristic space and time scales of the phenomena observed, which yields the velocity of interest $\vartheta _0 = x_0/t_0$. The magnitude of the thermal velocity for the species $\alpha$ is denoted $v_{0,\alpha }$ with $v_{0,\alpha }=k_\textrm {B} T_0/m_\alpha$. Due to the different masses, the thermal velocity of the electrons is not that of the ions. The reference thermal velocity $v_0$ will be defined by the ion one $v_{0}^2=k_\textrm {B} T_0/m_\textrm {i}$, hence $v_{0,\textrm {e}}=v_0/\varepsilon$ and $v_{0,\textrm {i}}=v_0$, where $\varepsilon ^2=m_\textrm {e}/m_\textrm {i}$. Finally, the scale of the electromagnetic field is denoted $(E_0,B_0)$, the particle current scale being defined as $J_0 = q n_0 u_0$. The dimensionless variables are defined according to
the collision frequencies verifying
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 $\varepsilon ^{-1}=\sqrt {m_\textrm {i}/m_\textrm {e}}$. Finally, the ions are almost unaffected by the collisions with the electrons. The relaxation of the ionic distribution function towards that of the electrons defines the largest time scale, by a factor $\varepsilon ^{-1}$ compared to the relaxation towards the thermodynamical equilibrium.
The dimensionless ionic and electronic Vlasov equations can be rewritten as (keeping the same notations for dimensionless variables)
together with the dimensionless Maxwell system written
This system is written using the following dimensionless parameters:
The dimensionless Maxwellians are defined by
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 electrons and ions. Accordingly, the relaxation of the electron mean velocity and temperature towards those 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 $\xi$ is intended to provide a measure of how the electronic and ionic dynamics are resolved. The choice $\xi =1$ means that the system is assumed to evolve at a speed comparable to the ionic thermal velocity $v_0$, while $\xi \varepsilon =1$ performs a rescaling of this typical velocity to the electron microscopic velocity. Setting $\xi =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 $\eta =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 $\beta M = \xi$ is also common in single-fluid plasma models. This amounts to assuming 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 latter assumption is in line with the Alfvén frozen theorem (Moreau Reference Moreau1990; Davidson Reference Davidson2001; Schnack Reference Schnack2009; Freidberg Reference Freidberg2014) characteristic of ideal MHD models: the magnetic field is frozen into the plasma and transported by its flow.
The derivation of reduced models consists of identifying small dimensionless parameters and letting 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 in explaining 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 vanishingly small parameters permits the derivation of 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
3.1.1 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 figure 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
for the evolution of the ions and electrons coupled to the dimensionless Maxwell system defined by (2.21).
From the scaling relations stated by (2.18a–c), discarding the inter-species collisions makes 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 intra- 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.
3.1.2 Handling the fluid limit
To identify easily a fluid regime, the distribution function is decomposed into a Maxwellian $\mathcal {M}_{n_\alpha ,u_\alpha ,T_\alpha }$ and a deviation from this Maxwellian ${\kappa } g_\alpha$ according to
the deviation verifying
With this decomposition, the Vlasov–Boltzmann equations (3.1) can be recast into a hydrodynamic set of equations with kinetic corrections, depending on the moment of the deviation $g_\alpha$, yielding
with
for the ions, and an equivalent system for the electrons,
with
These two systems are coupled to a set of equations (the Maxwell system (2.21)) driving the changes in the electromagnetic field, as well as an equation for the evolution of the deviations $g_\textrm {e}$ and $g_\textrm {i}$. The construction of these equations will be omitted in the present work; we refer for instance to Crouseilles & Lemou (Reference Crouseilles and Lemou2011), Crestetto et al. (Reference Crestetto, Crouseilles and Lemou2012) and De Cecco et al. (Reference De Cecco, Deluzet, Negulescu and Possanner2017) for examples of their derivation.
3.1.3 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:
This amounts to assuming 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 $\alpha \sim \lambda$ is essential to recover the low-frequency Ampère law, derived by neglecting the displacement current. This equation being common to MHD models, the quasi-neutral limit encompasses these two assumptions. With the vanishing of this generalized dimensionless Debye length $(\lambda ,\alpha )\to 0$, the Maxwell system degenerates into
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 either of these degenerate Gauss and Ampère equations. The remaining occurrence of the electric field is limited to the Faraday equation (3.6b). 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 (3.6): the electric field satisfying this system may be augmented by any gradient of a scalar potential (see also Degond & Deluzet (Reference Degond and Deluzet2017) for further details).
In the quasi-neutral limit, the electric field is provided by the particle current $J$ rather than the displacement current ($\partial E/\partial t$ originally present in Ampère'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.
3.2 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 dynamics of the electrons. The velocity of interest is the kinetic velocity of the electrons. This amounts to setting $\vartheta _0 = v_0/\varepsilon$ or equivalently $\xi \varepsilon = 1$, yielding
The collisions are assumed to be ineffective on the characteristic time scale $\kappa \varepsilon \gg {1}$, which amounts to neglecting the collision operator in (3.7), in particular for the ions, owing to $\varepsilon \ll 1$.
Performing the quasi-neutral limit $(\lambda =\alpha )\to 0$, the system at hand here is recast into (see Degond et al. Reference Degond, Deluzet and Doyen2015)
First, note that the formal time derivative of the Faraday equation (3.11) together with the curl of Ampère's law yields
which outlines that the electric field is known up to the gradient of a potential in this system. In Degond, Deluzet & Savelief (Reference Degond, Deluzet and Savelief2012a) and Degond et al. (Reference Degond, Deluzet and Doyen2015) 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 (3.3b) and (3.4b), yields
with, for any species $\alpha$ ($\alpha =\textrm {e}$ for the electrons and $\textrm {i}$ for the ions),
Inserting the identity (3.15a), together with $n=n_\textrm {i}=n_\textrm {e}$, into (3.14) gives
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 $\beta =1$ as well as $\eta =1$ which amounts to considering the thermal energy on a par with the electric one. This yields the following definition of the quasi-neutral model:
Note that the electric field provided by (3.17c) enforces a divergence-free particle current, or more precisely $\partial _t (\boldsymbol {\nabla } \boldsymbol {\cdot } J) = 0$. This yields, because of the continuity equation,
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 semi-discretization, the vanishing of the charge density, rather than its time double derivative, may be recovered (see Degond et al. Reference Degond, Deluzet and Doyen2015). Note also that the evolution of the ions only brings a marginal correction, proportional to $\varepsilon ^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 (Kingsep, Chukbar & Yank'kov Reference Kingsep, Chukbar and Yank'kov1990). In particular, this modelling accounts for the inertia of electrons. A noticeable difference from the electron MHD (see § 3.3.1) lies in the kinetic description of the plasma. An asymptotic-preserving method is proposed in Degond et al. (Reference Degond, Deluzet and Doyen2015) to bridge this quasi-neutral model and the Vlasov–Maxwell system. The properties of this quasi-neutral plasma description are investigated in Tronci & Camporeale (Reference Tronci and Camporeale2015) by means of a linear stability analysis.
3.2.2 A hybrid formulation of the Hall MHD
Hybrid modelling (Yin et al. Reference Yin, Winske, Gary and Birn2002; Winske et al. Reference Winske, Yin, Omidi, Karimabadi, Quest, Büchner, Scholer and Dum2003; Tronci et al. Reference Tronci, Tassi, Camporeale and Morrison2014) 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 (2.18a–c). 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 $\xi = 1$ resulting in the following system for the plasma:
The fluid limit for the electrons is selected assuming $(\varepsilon \kappa ) \ll 1$ meaning that the number of electron collisions during the typical time is large. The quasi-neutrality of the plasma amounts to setting $\lambda =\alpha \ll 1$. To overcome the degeneracy of the Maxwell system in the quasi-neutral limit, the electronic momentum is harnessed to provide the so-called generalized Ohm's law. The electronic system can be recast into
where the viscous stress tensor $\sigma _\textrm {e}$ as well as the thermal conductivity $\mu _\textrm {e}$ are defined in § 3 and, owing to the quasi-neutrality assumption, $n_\textrm {e} = n_\textrm {i} = n$.
The dynamics described by these equations is stiff. This is due to the smallness of $(\xi \varepsilon )$ 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\varepsilon )\ll 1$ gives rise to the following equilibria:
The classical massless approximation for the electrons is recovered with the generalized Ohm's law and a homogeneous electronic temperature. The definition of the mean velocity $u_\textrm {e}$ is derived from the particle current density $J= n (u_\textrm {i} - u_\textrm {e})$ together with Ampère's law (3.10), yielding
The hybrid plasma modelling is written (assuming $\eta =\beta =1$)
with
The derivation of a similar model is proposed in Acheritogaray et al. (Reference Acheritogaray, Degond, Frouvelle and Liu2011) with numerical investigations in Degond et al. (Reference Degond, Deluzet, Dimarco, Gallice, Santagati and Tessieras2011a).
3.3 A fluid hierarchy of quasi-neutral models
3.3.1 The electron MHD system
This model is obtained by letting $\kappa \to 0$ in (3.1). This yields the following set of equations for the electrons:
A similar system is derived for the ions, however with $\varepsilon =1$ and $\eta$ replaced by $-\eta$. These two sets of conservation laws are coupled to the Maxwell system (2.21). Performing the quasi-neutral limit in this system $(\lambda =\alpha )\to 0$ and focusing on the electronic dynamics with $\xi \varepsilon =1$ yields the quasi-neutral bi-fluid Euler–Maxwell system. This model is similar to that of § 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 Kingsep et al. (Reference Kingsep, Chukbar and Yank'kov1990) though different terminology may be used by other authors). It is implemented and investigated numerically in the framework of Asymptotic-Preserving methods in Degond et al. (Reference Degond, Deluzet and Savelief2012a).
3.3.2 The Hall MHD regime
The Hall MHD regime (see Lighthill Reference Lighthill1960; Witalis Reference Witalis1986; Schnack Reference Schnack2009) is recovered from the assumptions of the previous section but with a typical velocity equal to the plasma mean flow yielding $\xi =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_\textrm {i}$. The other parameters obey the classical scaling relations of MHD models: $\eta =1$ and $\beta =1$.
In the drift regime ($M\varepsilon \to 0)$, the electronic energy reduces to the internal energy
with the electronic momentum and energy verifying
The generalized Ohm's law (3.26b) is harnessed to compute the electric field. The electronic velocity $u_\textrm {e}$ is substituted by $u_\textrm {e}=u - J/n$.
The electric field is computed via the generalized Ohm's law, giving rise to
The first term on the right-hand side of equation (3.27) 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 (3.6b), the magnetic induction equation can be constructed, with
Finally the plasma mass density, momentum and total pressure $p$ and energy $W$ defined by
verify
The Hall velocity $\text {v}_\text {H}$, which can be interpreted as the electron velocity in the ion frame, is defined by
The ideal MHD equations are classically written under a conservative form using the system total pressure and energy:
writing the system as
This set of equations is supplemented with the electronic energy conservation (3.26c).
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 from departing 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 considered to derive the ideal MHD regime.
The drift approximation operated for the electrons amounts to vanishing of the electronic Mach numbers $(M\varepsilon )$. The scale separation introduced by the small electron to ion mass ratio $\varepsilon$ is not always sufficient to consider this limit independently of vanishing ionic Mach numbers $M \to 0$. For low ionic Mach number, a low-frequency filtering may be operated performing the limit of vanishing electronic Mach numbers jointly with ionic Mach numbers. This asymptotic defines the massless MHD regime (Besse et al. Reference Besse, Degond, Deluzet, Claudel, Gallice and Tessieras2004).
4 The electrostatic regime and the Boltzmann relation
4.1 Electrostatic limit of the Maxwell system
The electrostatic regime is recovered from the dimensionless Maxwell system (2.21) by letting $\alpha$ go to zero, independently of $\lambda$. This assumption shall be interpreted as a typical velocity negligible compared to the speed of light. From Ampère's law
the limit $\alpha \to 0$ provides $\boldsymbol {\nabla } \times B=0$. Nonetheless, the right-hand side of (4.1) remains an undetermined form. Therefore, Ampère'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 $\boldsymbol {\nabla }_x \times B = 0$ together with $\boldsymbol {\nabla }_x \boldsymbol {\cdot } B = 0$ yields $\partial _t B = 0$. Inserting this property into the Faraday equation (2.21b) provides an electrostatic electric field: $\boldsymbol {\nabla } \times E = 0$. Furthermore, the undetermined form in (4.1) is divergence-free. Therefore, computing the divergence of Ampère's law provides
which is a well-posed problem for the electric field under the condition $\boldsymbol {\nabla }_x \times E = 0$. Note that, owing to the continuity equation
originating from the conservation of the particle densities (3.3a) and (3.4a), the divergence of Ampère's law is equivalent to the time derivative of Gauss's law, with
Therefore, in the electrostatic regime, the Gauss equation is used to compute the electric field.
4.2 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 following is therefore
coupled to the Gauss equation
$\phi$ being the electrostatic potential, with $E = - \boldsymbol {\nabla } \phi$ and $n_\alpha = \int f_\alpha \,\textrm {d} v$.
The quasi-neutrality of the plasma is recovered for vanishing scaled Debye length $\lambda \to 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 using the equation of the electronic momentum conservation. The electric field is computed in order for this conservation to be satisfied. In Degond et al. (Reference Degond, Deluzet, Navoret, Sun and Vignal2010) an equivalent approach is proposed. It consists of using the time derivatives of Gauss's law to produce
In the quasi-neutral limit ($\lambda \to 0$), the electrostatic field is the Lagrangian multiplier of the quasi-neutrality constraint
From the system (3.4) (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:
Resuming the scaling relations of § 3.2.1, $\xi =1/\varepsilon$ and $\kappa \varepsilon > 1$, assuming that the ions are at rest, the evolution of the charge density $n_\textrm {i} -n_\textrm {e}$ is governed by (4.8) with
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 $\phi$ is obtained by inserting this relation into (4.6) and passing to the limit $\lambda \to 0$.
This yields the following quasi-neutral kinetic plasma description:
According to the values of $\kappa$, the collision term in (4.11a) may be disregarded, defining therefore a non-collisional kinetic description. Contrariwise, letting $\kappa \to 0$, a fluid description for the electrons may be derived, with
Note that a similar equation to (4.10), but with $P_\textrm {e} = p_\textrm {e} \mathbb {I}{\textrm {d}}$, may be worked out of the conservation of the electronic density (4.12a) and momentum (4.12b). This outlines that the electronic dynamics, in particular the electronic speed of sound, is resolved in this model. Comparable models are implemented and numerically explored in the context of asymptotic-preserving methods in Crispel, Degond & Vignal (Reference Crispel, Degond and Vignal2007) for fluid plasma description and Manfredi, Hirstoaga & Devaux (Reference Manfredi, Hirstoaga and Devaux2011) for kinetic equations.
4.3 Quasi-neutral models at the ionic scale
The typical velocity is chosen to be that of the ions with $\xi =1$ for the kinetic descriptions of the ions and $\xi =M$ for macroscopic models.
The hybrid modelling investigated in § 3.2.2 is defined by the scaling relations $\xi =1$, $(\varepsilon \kappa )\ll 1$ and $\eta =1$. The equilibria stated by (3.22) yield
with a homogeneous electronic temperature. This equation is integrated to provide the so-called Boltzmann relation
$n_0$ being a constant (independent of the space variable $x$) that should be determined from adequate conditions (Hagelaar Reference Hagelaar2007). Due to the Boltzmann relation, the quasi-neutral limit is no longer singular. Indeed plugging the Boltzmann relation (4.14) into the Gauss equation yields
This equation is not degenerate for the computation of the electric potential for vanishing $\lambda$. Indeed, the nonlinear part of the equation provides a means of computing $\phi$ in the quasi-neutral limit. This property is thoroughly investigated in Degond et al. (Reference Degond, Liu, Savelief and Vignal2012b).
The hybrid electrostatic model may be recast into
Letting $\kappa \to 0$ together with $\kappa /(M\varepsilon )\ll 1$ and $\xi =M$ yields the quasi-neutral fluid model:
In the models (4.16) and (4.17) 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\varepsilon ) \to 0$.
5 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 MHD model. The second framework is dedicated to electrostatic fields. In this context, the asymptotic analysis permits the derivation of 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 paper unravel different asymptotic parameters explaining the transition from one model to the other. The effort conducted in the present work consists of 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 scales, therefore a velocity scale as well, that the numerical method is aimed at capturing. This is related to the parameter $\xi$ used for the asymptotic analysis. Regarding the quasi-neutral modellings investigated, 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, in either 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 organization of § 4 is aimed at emphasizing this feature.
The second parameter, already established in prior works (see Degond et al. Reference Degond, Deluzet, Navoret, Sun and Vignal2010, Reference Degond, Deluzet and Savelief2012a, Reference Degond, Deluzet and Doyen2015), is the generalized scaled Debye length $\lambda$. It actually encompasses the scaled Debye length and the ratio of the typical velocity to the speed of light. Vanishing the generalized Debye length amounts to filtering 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. These 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\varepsilon )\ll 1$. In single-fluid plasma representation, the fast electron dynamics 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 (Chodura Reference Chodura, Post and Behrisch1986; Stangeby Reference Stangeby2000; Grismayer et al. Reference Grismayer, Mora, Adam and Héron2008; Manfredi et al. Reference Manfredi, Devaux, Shukla, Eliasson and Stenflo2008; Drouin et al. Reference Drouin, Gremillet, Adam and Héron2010). Accounting for this phenomenon is possible by selecting the appropriate typical velocity to resolve or filter the fast electron dynamics.
Finally the fluid assumption is classically related to a vanishing of the Knudsen number $\kappa \ll 1$ accounting for the relaxation of the distribution function towards the local thermodynamic equilibrium. The interplay of these four dimensionless parameters ($\xi$, $\lambda$, $M\varepsilon$, $\kappa$) 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 the study presented here draws the guidelines for the derivation of numerical methods implementing local up-scalings, therefore widening the use of numerical methods discretizing these reduced models.
Acknowledgements
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.C. 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.
Editor Giovanni Manfredi thanks the referees for their advice in evaluating this article.
Appendix A. Micro–macro decomposition, computation of the viscous terms
A.1 Introduction, definitions, elementary properties
The analyses carried out in this appendix are developed in the electrostatic framework and specified for electrons. The extensions for either electromagnetic fields or ions are straightforward and are therefore omitted for conciseness.
We first introduce the projector onto the Maxwellian. For any smooth function $\varphi$, the projector onto the Maxwellian denoted $\mathcal {M}_\textrm {e}$ is
For any smooth function $\varphi$, the projector onto $\mathcal {M}_\textrm {e}$, denoted $\varPi _{\mathcal {M}_\textrm {e}}$, is defined as (see Bennoune, Lemou & Mieussens Reference Bennoune, Lemou and Mieussens2008; Crestetto et al. Reference Crestetto, Crouseilles and Lemou2012)
where $\left \langle \varphi \right \rangle = \int \varphi \,\textrm {d} v$. For $k= 1,\ldots ,D_v$, we have the following properties:
Furthermore, if $g_\textrm {e}$ satisfies the micro–macro decomposition (3.2), the following identity holds true:
A.2 Computation of the deviation to the Maxwellian
The aim here is to characterize $g_\textrm {e}$ or more specifically an approximation to the first order in $\kappa$. Inserting the micro–macro decomposition into the Vlasov–Boltzmann equation (4.5b) yields
where
This provides, using properties (A 2),
It follows that
with
From properties (A 2), we can state the expression of the deviation to the Maxwellian:
A.3 Computation of the viscous terms
The viscous terms are defined by
Following the characterization (A 7) of $g_\textrm {e}$, we can write
Inserting in these definitions the following identities:
we obtain