Issue 
A&A
Volume 635, March 2020



Article Number  A87  
Number of page(s)  15  
Section  The Sun and the Heliosphere  
DOI  https://doi.org/10.1051/00046361/201834686  
Published online  12 March 2020 
Consistent transport properties in multicomponent twotemperature magnetized plasmas
Application to the Sun atmosphere
^{1}
Centre de Mathématiques Appliquées, Ecole Polytechnique, CNRS, Palaiseau, France
email: quentin.wargnier@polytechnique.edu
^{2}
Laboratoire de Physique des Plasmas, Ecole polytechnique, CNRS, Palaiseau, France
^{3}
Aeronautics and Aerospace Department, von Karman Institute for Fluid Dynamics, RhodeSaintGenèse, Belgium
^{4}
NASA Ames Research Center, Moffett Field, CA, USA
Received:
19
November
2018
Accepted:
25
December
2019
Aims. We present a fluid model that has been developed for multicomponent twotemperature magnetized plasmas in chemical nonequilibrium for the partially to fully ionized collisional regimes. We focus on transport phenomena with the aim of representing the atmosphere of the Sun.
Methods. This study is based on an asymptotic fluid model for multicomponent plasmas derived from kinetic theory, yielding a rigorous description of the dissipative effects. The governing equations and consistent transport properties are obtained using a multiscale ChapmanEnskog perturbative solution to the Boltzmann equation based on a dimensional analysis. The mass disparity between free electrons and heavy particles is accounted for, as well as the influence of the electromagnetic field. We couple this model to the Maxwell equations for the electromagnetic field and derive the generalized Ohm’s law for multicomponent plasmas. The model inherits a wellidentified mathematical structure leading to an extended range of validity for the Sun’s atmospheric conditions. We compute consistent transport properties by means of a spectral Galerkin method using the LaguerreSonine polynomial approximation. Two nonvanishing polynomial terms are used when deriving the transport systems for electrons, whereas only one term is retained for heavy particles.
Results. In a simplified framework where the plasma is fully ionized, we compare the transport properties for the lower solar atmosphere to conventional expressions for magnetized plasmas attributed to Braginskii, showing a good agreement between both results. For more general partially ionized conditions, representative of the lower solar atmosphere, we compute the muticomponent transport properties corresponding to the species diffusion velocities, heavyparticle and electron heat fluxes, and viscous stress tensor of the model for a heliumhydrogen mixture in local thermodynamic equilibrium. The model is assessed for the 3D radiative magnetohydrodynamic simulation of a pore at the Sun photosphere. The resistive term is found to dominate mainly the dynamics of the electric field at the pore location. The battery term for heavy particles appears to be higher at the pore location and at some intergranulation boundaries.
Key words: plasmas / Sun: chromosphere / Sun: atmosphere
© Q. Wargnier et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The lower atmosphere of the Sun is a complex and dynamic layer where the plasma is found in a wide range of different regimes – from weakly ionized and nonmagnetized at the bottom of the photosphere to fully ionized and magnetized at the top of the transition region. In the Sun’s chromosphere, the pressure varies from a thousand pascals just above the photosphere to a few pascals in the high chromosphere (see Vernazza et al. 1981). Similarly, the magnitude of the magnetic field is large in active regions, around thousands of gauss in sunspots, whereas it is just a few gauss in quietSun regions, as shown in Wiegelmann et al. (2014). It is still a great challenge nowadays to develop a unified model that can be used for both partially and fully ionized regimes under the large disparity of plasma parameters in the lower atmosphere.
The study of partially ionized plasmas in the presence of a magnetic field, such as in prominences and the lower atmosphere of the Sun, demands models that are beyond the ideal singlefluid magnetohydrodynamic (MHD) description. Phenomena necessary to fully understand the behavior of plasmas in the Sun chromosphere, such as Cowling’s resistivity, thermal conduction, heating due to ionneutral friction, heat transfer due to collisions, charge exchange collisions, and ionization energy losses, are usually disregarded in ideal MHD models, or only described by means of adhoc terms.
Two types of fluid models for studying the lower Sun atmosphere are found in the literature. First, the singlefluid MHD description considers the plasma as a conducting fluid in the presence of a magnetic field. It has the main drawback of assuming thermal equilibrium conditions, where all the species are considered to be at the same temperature. This model is assumed to be valid at the photosphere, in a highly collisional framework, allowing us to study the formation of magnetic field concentrations at the solar surface in sunspots, magnetic pores, and the largescale flow patterns associated with them (see Hartlep et al. 2012). It is also used for simulating the lower part of the atmosphere of the Sun, for example, incorporating subgridscale turbulence models for the transport of heat and electrical resistivity, as presented by Kitiashvili et al. (2015). The full MHD equations are solved in Martinez Sykora et al. (2015) accounting for nongrey radiative transfer and thermal conduction outside local thermodynamic equilibrium in order to study the effects of the partial ionization of the Sun chromosphere.
More recently, multifluid MHD models have been used to represent the nonequilibrium conditions of the chromosphere, based on continuity, momentum, and energy conservation equations for each species considered in the mixture (see Leake et al. 2013; Khomenko et al. 2014; Khomenko 2017; Alvarez Laguna et al. 2016, 2018; Ni et al. 2018). These models lead to very stiff systems that are difficult to solve numerically as they exhibit characteristic times that range from the convective and diffusive times of each fluid down to the collisional and chemical kinetics time scales, Alvarez Laguna et al. (2016). Leake et al. (2013) performed a multifluid simulation of magnetic reconnection for a weakly ionized reacting plasma, with a particular focus on the solar chromosphere, by considering collisional transport, chemical reactions between species, as well as radiative losses. Braginskii (1965) has derived rigorous expressions for the transport properties of fullyionized plasmas starting from the Boltzmann equation. Khomenko et al. (2014) proposed a model for the description of a multicomponent partially ionized solar plasma. Deriving rigorous transport properties for such multifluid model is complex and, so far, the theory has not yet been developed to the same level of accuracy as Braginskii’s.
In this paper, we propose a novel approach for studying the lower atmosphere of the Sun that is neither a singlefluid nor a multifluid MHD model, but a multicomponent driftdiffusion model derived by Graille et al. (2009) from kinetic theory. This approach is able to capture most of the multifluid phenomena, that is different velocities between species, collisional exchange of mass, momentum and energy, chemical reactions, thermal nonequilibrium, and magnetized transport. This system of equations is less stiff as it solves for only one momentum equation, like in the singlefluid MHD approach. The electrons are assumed to move at non relativistic speeds in a collisional fluid regime. The transport properties are retrieved through a generalized ChapmanEnskog solution to the Boltzmann equation using a multiscale perturbation method. These developments lead to a model with an extended range of validity from partially to fully ionized plasmas, with or without the presence of a magnetic field. We couple this model to the Maxwell equations for the electromagnetic field and derive the corresponding generalized Ohm’s law for multicomponent plasmas. As in Braginskii’s theory, our model includes anisotropy in the transport properties of electrons that is created by the magnetic field. These properties are computed by solving for the integrodifferential systems presented by Scoggins et al. (2016). We use a spectral Galerkin method based on the LaguerreSonine polynomial approximation previously studied in depth for various applications (see Devoto 1969; Ferziger & Kaper 1973; Woods 1995; Zhdanov 2002; Magin & Degrez 2004a; Bruno et al. 2011; Capitelli et al. 2013). The transport systems are implemented in the Mutation++ library that compiles stateoftheart transport collision integral data for the different pairs of species in the mixture, see Scoggins & Magin (2014). In a simplified framework where the plasma is fully ionized, we compare the transport properties for the lower atmosphere of the Sun to the conventional expressions for magnetized plasmas attributed to Braginskii. For more general partially ionized conditions representative of the Sun’s lower atmosphere, we compute the muticomponent transport properties corresponding to the species diffusion velocities, heavyparticle and electron heat fluxes, and viscous stress tensor for a heliumhydrogen mixture in local thermodynamic equilibrium. Finally, the model is assessed for the 3D radiative MHD simulation of a pore in the highly turbulent upper layer of the solar convective zone. We compute the thermal conductivity, electrical conductivity, species diffusion coefficients, and components of the generalized Ohm’s law based on an existing fluid simulation. We draw conclusions about the importance of the contribution of its components, in particular, of the resistive and battery terms. The coupling of the model developed for transport with a flow solver is beyond the scope of this paper.
The structure of the paper is as follows. In Sect. 2, the nondimensonal analysis used for the generalized ChapmanEnskog expansion is presented, together with the multicomponent driftdiffusion model for twotemperature magnetized plasmas, the transport fluxes, and the generalized Ohm law. In Sect. 3, we describe the mixture considered, the conditions representative of the lower Sun atmosphere and the method used for computing the transport properties. In Sect. 4, we verify the model proposed on a fullyionized case by comparing the results with those obtained by means of Braginskii’s theory. Finally, in Sect. 5, we discuss all the transport properties for a partially ionized case. Additionally, we compute the transport properties and the components of the generalized Ohm law for 3D radiative MHD simulations of a pore in the low Sun atmosphere.
2. Driftdiffusion model for multicomponent plasmas
In this section, we present the multicomponent driftdiffusion model for twotemperature magnetized plasmas. It was derived from kinetic theory by Graille et al. (2009) as a generalized ChapmanEnskog solution to the Boltzmann equation, using a multiscale perturbation method based on a dimensional analysis. Additionally, we compare this model to the multifluid description widely used for the Sun’s lower atmosphere.
2.1. Multiscale analysis of the Boltzmann equation
We consider a multicomponent plasma composed of electrons, denoted here by the index 𝔢, and heavy particles (atoms and molecules, neutral or ionized), denoted by the subscript 𝔥. The species are assumed to be point particles, neglecting their internal energy. We combine the equations derived by Graille et al. (2009), Magin et al. (2009) for the fully magnetized case and the Maxwellian regime for reactive collisions.
In order to apply the ChapmanEnskog method, Graille et al. (2009) perform a multiscale analysis on the nondimensional Boltzmann equations for electrons and heavy particles. The order of magnitude of the different terms in the Boltzmann equation is studied through carefully chosen reference quantities. In the asymptotic fluid limit, the Knudsen number is assumed to be as the same order of magnitude as the square root of the mass ratio between electrons and heavy particles, defined as , where m_{𝔢} and m_{𝔥} are the mass of electrons and reference heavy particles respectively. This small parameter drives thermal nonequilibrium between the electron and heavyparticle baths. In the strongly magnetized regime, the Hall parameter is assumed to scale as ε^{0} = 1. The species distribution functions are expanded in the multiscale perturbation parameter ε following Enskog’s approach. As opposed to Braginskii, no assumption is made a priori on the zeroorder distribution function. The asymptotic analysis of the Boltzmann equation is performed at successive orders of ε. The main results occurring at different timescales are summarized in Table 1.
Time scales hierarchy and macroscopic equations derived using the ChapmanEnskog method, Graille et al. (2009).
In time scales of order , the electron population thermalize at the temperature T_{𝔢}. The electron distribution function is a MaxwellBoltzmann distribution obtained by solving the electron Boltzmann equation at the order ε^{−2}. At order ε^{−1} that corresponds to the time scale , heavy particles thermalize at temperature T_{𝔥}. At the zeroth order ε^{0} that corresponds to the convective time scales, Euler equations for heavy particles and zeroorder driftdiffusion equations for electrons are obtained. Finally, at order ε, corresponding to the diffusive time scale, we obtain NavierStokes equations for heavy particles and firstorder driftdiffusion equations for electrons. It is important to mention that in Braginskii’s approach, the macroscopic equations are retrieved by taking moments of the Boltzmann equation by considering adhoc perturbations to the MaxwellBoltzmann distributions. Only a correct scaling deduced from dimensional analysis can yield a sound multicomponent treatment.
2.2. Multicomponent equations
We denote by symbol H the set of indices for the heavy particles of the mixture considered. First, the mass conservation equations for electrons and heavy particles are, as follows,
Here, ρ_{𝔢} is the mass density of electrons, ρ_{i} is the mass density of heavy particle i ∈ H, v_{𝔥} is the heavyparticle hydrodynamic velocity that has been chosen as the velocity reference frame, V_{𝔢} is the electron diffusion velocity, V_{i}, i ∈ H is the diffusion velocity of heavy particle i in the reference frame, such that ∑_{i ∈ H}ρ_{i}V_{i} = 0. Quantities ω_{𝔢} and ω_{i}, with i ∈ H, are respectively the chemical production rates of electrons and heavy particles.
Second, a momentum equation for all the particles within the plasma is found to be
where, p = p_{𝔢} + p_{𝔥} is the mixture pressure, which is the sum of the partial pressure of electrons p_{𝔢} and heavy particles p_{𝔥}, Π_{𝔥} is the viscous stress tensor. The quantity ρ_{𝔥} is the mass density of heavy particles, 𝔫 = 𝔫_{𝔢} + 𝔫_{𝔥} is the mixture number density (𝔫_{𝔢} the number density of electrons and 𝔫_{𝔥} the number density of heavy particles), 𝔫𝔮 is the mxiture charge defined by 𝔫𝔮 = 𝔫_{𝔢}𝔮_{𝔢} + ∑_{i ∈ H}𝔫_{i}𝔮_{i}, where 𝔮_{i} is the charge of species i, E the electric field, and I the total current density defined as
where J_{𝔥} is the heavyparticle conduction current density, J_{𝔢} is the electron conduction current density, and B is the magnetic field. Note that electrons only participate in the momentum balance through the total pressure gradient and Lorentz force due to the large mass disparity between electrons and heavy particles.
Conservation of energy for free electrons and heavy particles is read as
where ρ_{𝔢}e_{𝔢} and ρ_{𝔥}e_{𝔥} are the internal energies of electrons and heavy particles respectively, and q_{𝔢} and q_{𝔥}, the electron and heavyparticle heat fluxes, respectively. Quantity E′ = E + v_{𝔥∧}B is the electric field in the heavyparticle reference frame, and are the energy relaxation terms at order ε^{0} and ε, respectively. Quantities J_{𝔥} ⋅ E′ and J_{𝔢} ⋅ E′ are the power that is developed by the heavyparticle and electron current densities, and Ω_{𝔢} and Ω_{𝔥} are the energy production rate for electrons and heavy particles.
Summing Eqs. (5) and (6), with an equation for kinetic energy, the total energy equation is obtained as
where the total energy is , and I ⋅ E the power developed by the electromagnetic field.
The system of Eqs. (1)–(7) is coupled to the set of Maxwell’s equations (8) :
where ε_{0} is the vacuum permittivity and μ_{0} the vacuum permeability.
The electron transport fluxes such as the electron diffusion velocity V_{𝔢}, the electron heat flux q_{𝔢}, the electron current density J_{𝔢} are composed of two terms: 1 a term at the convective time scale, at order ε^{0}, corresponding to the Euler equations for heavy particles, and 2 a term which is a first order correction, at order ε, at the dissipative time scale corresponding to the NavierStokes equations for heavy particles. Similarly, the heavyparticle transport fluxes such as the heavyparticle diffusion velocity V_{i}, i ∈ H, heat flux q_{𝔥}, viscous stress tensor Π_{𝔥}, and current density J_{𝔥} are defined at the dissipative time scale, at the order ε of the generalized ChapmanEnskog expansion.
The governing Eqs. (1)–(8) differ from the multifluid models used for partially ionized plasmas. Whereas multifluid models consider one hydrodynamic velocity distinct for each species, here only one common hydrodynamic velocity is used for the heavy species while each species diffuse in this reference frame. In addition, the structure of the governing equations is symmetrizable hyperbolic (Graille et al. 2009), which can be regarded as an important property for the numerical discretization of the system. Nevertheless, it is necessary to close the model by computing the transport properties. This computation is presented in Sect. 5 for a heliumhydrogen mixture. Additionally, by using the definition of the total current density I and the Maxwell equations Eq. (8), a generalized Ohm’s law for this particular model is derived in Sect. 2.5.
2.3. Transport fluxes for heavy particles
With the same formalism that is used by Graille et al. (2009), we introduce some extra notations in order to express the anisotropic transport properties in the presence of a magnetic field. First, a unit vector for the magnetic field ℬ = B/B is defined as well as the three direction matrices
such that we have for any vector x in three dimensions
In the (x, ℬ) plane, the vector x^{‖} is the component of x that is parallel to the magnetic field and x^{⊥} is the perpendicular component. Therefore, we have x = x^{‖} + x^{⊥}. The vector x^{⊙} lies in the direction transverse to the (x, ℬ) plane. The three vectors x^{‖}, x^{⊥}, and x^{⊙} are then mutually orthogonal. The anisotropic transport coefficients are expressed by means of the matrix notation
In the weakly or unmagnetized plasma regime the parallel and perpendicular component of the transport coefficients are identical, leading to μ^{‖} = μ^{⊥}. In such framework, the transverse component of the transport coefficient (or in other words, perpendicular to the (x, ℬ) plane), is vanishing, leading to μ^{⊙} = 0.
We recall that the heavyparticle transport fluxes presented in Sect. (2.2) are:
The viscous stress tensor Π_{𝔥} is defined as
where η_{𝔥} is the viscosity of heavy particles. Then, the heavyparticle diffusion velocity V_{i} is given as
where D_{ij} is the multicomponent diffusion coefficient of heavy particles, d_{j} is the diffusion driving force of the particle j that is interacting with the heavy particle i, and χ_{𝔥j} is the heavyparticle thermal diffusion ratio of particle j. The diffusion driving force d_{j} is defined as
which is composed of three forces: 1 the force due to the gradient of the partial pressure ∂_{x}p_{j}, 2 the Lorentz force, and 3 F_{j𝔢} that is an average electron force acting on the heavy particle j. Quantity F_{j𝔢} belongs to the category of diffusion driving forces and allows for a coupling between the heavy particles and electrons (Graille et al. 2009). This average force is defined as
where and are anisotropic transport coefficients. Finally, the heavyparticle heat flux reads
where λ_{𝔥} is the heavyparticle thermal conductivity and ρ_{j}h_{j} is the enthalpy of heavy particle j. The second term of Eq. (14) corresponds to thermal diffusion.
In the previous transport fluxes, some of the usual terms can be identified. The viscous stress tensor for heavy particles Eq. (10) is proportional to the strain tensor. Similarly, the first term of V_{i} in Eq. (11) is a generalized Fick’s law where the flux is proportional to the diffusion driving force d_{j}. Also, the first term of the heavy particle heat flux q_{𝔥} in Eq. (14) is the usual Fourier law. In addition, V_{i} includes a term that is proportional to ∂_{x}lnT_{𝔥}, known as the Soret effect, described in Magin & Degrez (2004b), Giovangigli & Graille (2003, 2009).
In summary, the transport coefficients for heavy particles to be computed in the following sections are
In addition, the anisotropic transport coefficients associated to the coupling terms between electron and heavy particles are
2.4. Transport fluxes for electrons
The electron transport fluxes are
The electron diffusion velocity is defined as
where is the tensor for the diffusion coefficient of electrons, d_{𝔢} is the electron diffusion driving force and is the electron thermal diffusion ratio. The electron diffusion velocity V_{𝔢} is split into two terms : 1 the terms proportional to d_{𝔢} and to ∂_{x}lnT_{𝔢} are at order ε^{0}, at the heavy particle convective timescale and 2 the terms proportional to V_{i} are at order ε at the heavy particle dissipative timescale. The electron diffusion driving force is defined as
which is composed of two forces: 1 the force due to the gradient of the partial pressure of electron ∂_{x}p_{𝔢} and 2 the Lorentz force. The electron heat flux reads:
where is the electron thermal conductivity tensor and ρ_{𝔢}h_{𝔢} is the enthalpy of electrons. Quantity q_{𝔢} is split into two terms: 1 the terms proportional to ∂_{x}T_{𝔢} and V_{𝔢} are at the heavyparticle convective timescale, and 2 the terms proportional to V_{i} that are at the heavyparticle dissipative timescale.
As in the heavy species transport properties, some usual terms can be identified that are Fick’s and Fourier’s laws. Additionally, terms that are coupled to the heavyparticle diffusion are present at the first order of the generalized ChapmanEnskog expansion. In summary, the anisotropic transport coefficients associated to the transport fluxes for electrons are
In this section, a list of the transport fluxes and the corresponding transport coefficients has been presented. The method used for computing the latter will be presented in the next Sect. 3.
2.5. Generalized Ohm’s law
In the following, we derive a general expression for Ohm’s law in the previous set of governing Eqs. (1)–(8). In order to do so, we rewrite the expression of the electric current by grouping the terms in each of the driving forces. By doing this, we obtain a general algebraic expression for the electric field E as a function of the transport coefficients and the corresponding driving forces. Note that, in this derivation, the displacement current has been neglected.
With the electron diffusion velocity from Eq. (16) and the heavyparticle diffusion velocity from Eq. (11) we find the total current I as follows
where the multicomponent electromagnetic matrices are defined as:
and the tensor is defined as
After some algebra, using Eq. (18) and neglecting the displacement current, the general expression of the electric field is obtained as
The expression of the multicomponent electromagnetic matrices can be subdivided into two terms: (1) a term which depends on the coupled heavyparticle and electron transport properties, such as , which scales at the dissipative timescale for the heavy particles at order ε, and (2) a term which depends only on the electron transport properties , which scales as the convective timescale for the heavy particles at order ε^{0}.
Some usual terms can be identified in the general expression of the electric field Eq. (25). The first term of Eq. (25) is the resistive term, where the expression of the resistivity tensor is defined as
The second term and third term of Eq. (25), can be identified as a general expression of the battery term for a multicomponent plasma due to the pressure gradients of electrons and heavy particles. The fourth and last term of Eq. (25) are additional terms due to the presence of Soret and Dufour terms in the equations of the diffusion velocities Eqs. (11) and (16).
In Appendix B, a simplified fully ionized plasma case has been considered which leads to a simplified expression for the electric field. In this case, the multicomponent electromagnetic matrices can be simplified, and the usual expression of the electric field and magnetic induction equation are retrieved (see in Appendix B).
3. Methodology
A heliumhydrogen mixture, composed of 92% Hydrogen and 8% of Helium in mole fractions, which is typical of the Sun’s atmosphere, Asplund et al. (2009), is considered. The set of species considered in this mixture is denoted as
The heavy species such as carbon, oxygen or metals are not taken into account. We assume that they do not impact the transport properties as they are trace elements, i.e., the mole fractions are very small.
We study the transport coefficients for the previous mixture within a range of temperature, pressure, and magnetic field that is largely representative of the lower Sun atmosphere, as shown by Vernazza et al. (1981), Carlsson & Stein (1995): the temperature varies from 1000 K to 30 000 K, the pressure from 1 Pa to 10^{4} Pa, and the magnetic field from a few Gauss to thousands of Gauss, as shown by Wiegelmann et al. (2014). In the following, the plasma beta parameter is defined as β_{p} = 2μ_{0}p/B^{2}, where p is the total pressure of the plasma in Pascal and B is the magnetic field in Tesla.
Consequently, for a range of temperature between 1000 K and 30 000 K, two cases have been considered. The case A, where the total pressure is p = 10^{4} Pa, and β_{p} = 10. The case B, where the total pressure is p = 1 Pa, and β_{p} = 0.1. Finally, the case A is chosen as a thermally pressure dominated case whereas, the case B is a magnetically pressure plasma case. In this framework, in the case A, the transport coefficients are shown to be isotropic and in the case B the latter are shown to be anisotropic. The latter example may be representative for conditions that we can find 1 in a cool sunspot, where the average temperature is generally around 4000 K and the magnitude of the magnetic field is about 1000 G, as shown in Fröhlich & Lean (2004), Solanki (2003), and 2 in the quiet sun photosphere, in the lower and upper part of the chromosphere, where the temperature is varying from 5000 K to 10 000 K, and the pressure from 10^{4} Pa to 0.1 Pa, as shown in Vernazza et al. (1981), Russell (1929), Asplund et al. (2009) (in orders of magnitude).
Based on the chosen conditions, we compute the thermochemical equilibrium composition. The mole fraction and the ionization degree of the heliumhydrogen mixture S_{1} for the case A and case B are shown in Figs. 1–3. These results are obtained with a method that is based on the minimization of the Gibbs free energy with suitable mass balance constraints, developed by Scoggins & Magin (2015) in thermal equilibrium. The compositions that are shown in Figs. 1 and 2 will be used to study the transport properties in the following sections.
Fig. 1. Mole fraction of the heliumhydrogen mixture S_{1}, for P = 10^{4} Pa (case A), as a function of temperature. 

Open with DEXTER 
Fig. 2. Mole fraction of the heliumhydrogen mixture S_{1}, for P = 1 Pa (case A), as a function of temperature. 

Open with DEXTER 
Fig. 3. Ionization degree the heliumhydrogen mixture of S_{1}, for case A (full line) and case B (dashed line) as function of the temperature. 

Open with DEXTER 
Similarly as Scoggins et al. (2016), Magin & Degrez (2004a), the calculation of the transport coefficients is based on the solution of integrodifferential equations. In order to solve these equations, a spectral Galerkin method is applied. This method expands the coefficients in a series of orthogonal LaguerreSonine polynomials that are truncated at a given order of approximation. The calculation is thus reduced to a linear algebraic system of equations. As a result, the transport coefficients can be obtained by the resolution of linear systems with matrices that are known functions of the macroscopic parameters that are the field variables and the collision integrals between particles. The solution of these systems allows for the transport coefficients to be written as linear combinations of the collision integrals, which take into account the interaction potential for a collision between two particles. These linear combinations are derived by extending the definition and the calculation of bracket integrals introduced by Ferziger & Kaper (1973) or in Woods (1995), Zhdanov (2002), Balescu (1988) to the thermal nonequilibrium case, studied in depth by Kolesnikov (2003). According to Magin & Degrez (2004a), Kolesnikov (2003), Tirsky (1993), the transport coefficients involving collisions between heavy particles and electrons converge for expansions in second order non vanishing LaguerreSonine polynomials and higher. In this work, we use the third order Laguerre Sonine polynomials approximation in order to compute the transport properties. This method has been widely used in the literature, for example Bruno et al. (2011), Woods (1995), Magin & Degrez (2004b), Zhdanov (2002), Balescu (1988), Capitelli et al. (2013).
The explicit relations for and in terms of the solutions to the transport systems can be found in Scoggins et al. (2016) and in Appendix C. The heavy particle transport systems for η_{𝔥}, D_{ij}, χ_{𝔥j}, λ_{𝔥}, are found in Magin & Degrez (2004b) with the difference that the mole fractions are given in terms of heavy species only, excluding electrons. These are also presented in Appendix C. The transport collision integrals for the heliumhydrogen mixture are taken from Bruno et al. (2010).
4. Verification of the method in a fully ionized plasma case S_{2}
In order to verify the presented method, we perform a comparison with Braginskii’s transport properties. In the method of Braginskii (1965), the computation of the transport properties as well as the derivation of the governing equations are valid only for fully ionized plasmas. As it can be seen in Fig. 3, the heliumhydrogen mixture S_{1} can be considered to be fullyionized, mainly composed of S_{2} = {H^{+}, e^{−}}, when the temperature is higher than 15 000 K. The comparison will be thus performed in conditions where the mixture is S_{2} in a range of temperatures from T = 15 000 K to T = 30 000 K for the case A and case B. To illustrate the comparison, we focus on the properties , , η_{𝔥} and λ_{𝔥}, although the rest of them show similar behaviour.
On the one hand, in Braginskii (1965), the derivation of the governing equations can be summarized in three main steps: 1 A fully ionized ionelectron plasma is considered in a constant magnetic field, 2 The Landau collision operators are used, simplified by the Lorentz process, and 3 an adhoc ChapmanEnskog method is used based on the square root of the mass ratio between electron and ions, Balescu (1988). On the other hand, in Graille et al. (2009), a general multicomponent plasma that can be partially or fully ionized is considered in a constant magnetic field, the Chapman and Cowling collision operators highly studied in Ferziger & Kaper (1973), Woods (1995) are used and the ChapmanEnskog expansion is performed after a dimensional analysis of the Boltzmann equation. Finally, the two methods lead to distinct governing equations.
Although the governing equations between the two models are different, the integrodifferential systems for computing the transport properties are similar or even identical in the case of a fully ionized plasma. In both models, the anisotropic electron transport properties have the same integrodifferential systems. However, only the systems related to the parallel component of the heavy particle transport properties are identical to those from the model derived by Graille et al. (2009). Consequently, only the parallel component of the heavy particle transport properties can be compared with those from the model of Graille et al. (2009). This is due to the fact that both models are based on the ChapmanEnskog expansion. However, the differences result from the scale analysis from the Boltzmann equation that is carried out by Graille et al. (2009) before applying the expansion.
In both models, the transport coefficients are expanded in a series of orthogonal LaguerreSonine polynomials. The latter are written as linear combinations of collision integrals that are simplified by potential interactions, based on the usual Coulomb interaction screened by the Debyelength. This approximation assumes collisions with large impact parameters and small scattering angles. However, in Braginskii (1965), the series are truncated at the secondorder approximation whereas a thirdorder approximation has been performed in the case of the model presented. The expression of the transport coefficients depends on the mean collision times and defined as
where log(Λ) is the Coulomb logarithm defined by Spitzer (1963), k_{B}, Boltzmann’s constant, and Z, the charge number. The mean collision times as defined in Eq. (28), can be seen as a firstorder Chapman Cowling approximation of the collision time for electronion, and, ionion collisions, as shown by Woods (1995). Correction terms depending on Z are used for the computation of the transport coefficients. This method leads to simplified expressions of the transport coefficients that depend only on the mean collision times and the charge number of the fully ionized plasma considered, as presented by Balescu (1988), Woods (1995).
In Braginskii (1965) (see Eq. (4.37)), the parallel and perpendicular components of the electron thermal conductivity tensor are defined as
where Br denotes the computation of the transport coefficient as derived by Braginskii (1965). and ω_{𝔢} = 𝔮_{𝔢}B/m_{𝔢} and the values in brackets correspond to Braginskii’s coefficients for a charge number Z = 1.
Figures 4 and 5 show the parallel and perpendicular component of the electron thermal conductivity tensor , as function of the temperature, for the case A and the case B, for the fully ionized plasma S_{2}. Here, we compare the expressions from Braginskii Eq. (29) and Eq. (30) with those that are given by Scoggins et al. (2016) that are based on a thirdorder LaguerreSonine polynomials approximation. Strong similarities are obtained in all the considered cases. In Braginskii (1965), the components of the electron thermal conductivity tensor are underestimated leading to differences that are less than 20%. These differences are increasing at high temperatures. Similar results have been obtained for all the other electron transport properties.
Fig. 4. Parallel component of the electron thermal conductivity tensor for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER 
Fig. 5. Perpendicular component of the electron thermal conductivity tensor for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER 
Similarly, the parallel component of the heavy thermal conductivity and of the heavy particle viscosity of the model of Braginskii (1965, see Eq. (4.44)), have been compared with the expression from Scoggins et al. (2016) and Magin & Degrez (2004a). In Braginskii (1965), the heavy thermal conductivity and heavy viscosity are defined as
Figure 6 and 7 show the heavy thermal conductivity λ_{𝔥} and heavy viscosity η_{𝔥} respectively, in the same conditions as in Figs. 4 and 5. As before, strong similarities have been obtained in all the considered cases for the chosen conditions, which leads to differences that are smaller than 20%. In addition, it can be shown that the heavy transport properties from Braginskii (1965) are isotropic at the chosen conditions.
Fig. 6. Heavy thermal conductivity λ_{𝔥} for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER 
Fig. 7. Heavy particle viscosity η_{𝔥} for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER 
In summary, we can conclude that the proposed method is verified for the fully ionized case. The main differences that are obtained between the two models are due to 1the order of Laguerre Sonine polynomials that was used that is second order in Braginskii’s model, Balescu (1988), and third order in the proposed method, and 2 the nature of the collision operators used, Landau collision operators in the model of Braginskii, as opposed to Chapman and Cowling collision operators in the model of Graille et al. (2009). Additionally, the formulation of the transport properties that are considered in this paper are generalized for any type of partially ionized mixture.
5. Transport properties for a partially ionized heliumhydrogen plasma
5.1. Transport fluxes in thermochemical equilibrium
In order to simplify the analysis of the presented transport systems, we consider thermochemical equilibrium T_{𝔢} = T_{𝔥} = T, isobaric mixtures at rest. The total heat flux is entirely a function of the temperature gradient and magnetic field and may be written as
where the Soret and reactive thermal conductivities may be written as
where θ_{𝔢} and θ_{i}, i ∈ H are defined as
Here, θ_{𝔢} and θ_{i} correspond to diffusion velocities for a temperature gradient of 1, that are V_{𝔢} = θ_{𝔢}∂_{x}T and V_{i} = θ_{i}∂_{x}T. We compute all the transport properties for the heliumhydrogen mixture S_{1} for case A and case B.
Figures 8 and 9 present the parallel, perpendicular and transverse components of the electron thermal conductivity tensor as a function of the temperature, for both cases. According to Fig. 8 (case A), the perpendicular component is equal to the parallel component for the entire range of temperatures, thus, the electron thermal conductivity is isotropic. Indeed, the pressure forces are dominating the magnetic pressure forces, so the plasma is unmagnetized. On the other hand, in Fig. 9 (case B), for temperatures higher than T = 5000 K, the electron thermal conductivity is anisotropic since the magnitude of magnetic field is higher. This results in a transverse component that is higher than the perpendicular component of . Similar results have been obtained for the other anisotropic electron transport properties such as and .
Fig. 8. Components of the electron thermal conductivity tensor for the isotropic case A, at the thirdorder Laguerre Sonine polynomials, for the heliumhydrogen mixture S_{1} as a function of temperature. 

Open with DEXTER 
Fig. 9. Components of the electron thermal conductivity tensor for the anisotropic case B, at the thirdorder Laguerre Sonine polynomials, for the heliumhydrogen mixture S_{1} as a function of the temperature. 

Open with DEXTER 
Figure 10 shows the heavyparticle thermal conductivity λ_{𝔥}, as a function of the temperature, for the case A and case B. In Fig. 10, strong differences between the two cases for a temperature higher than 6000 K can be seen. In the case A, λ_{𝔥} increases from 1000 K to 9000 K, which is expected since λ_{𝔥} is an increasing function of the temperature. However, in the case A after 9000 K, λ_{𝔥} decreases. This decrease is due to the ionization of hydrogen. Indeed, the heavy particle thermal conductivity is related to a combination of the cross sections variations of all the heavy species in the mixture, which are proportional to the mole fractions of each heavy particles. This result is coherent with Fig. 1, which shows that the mole fraction of H is decreasing after 9000 K. Similar behavior as the case A have been observed for the case B, except that the ionization of H starts at 6000 K for this pressure. In Fig. 10, the second modulation observed around 12 000 K is due to the ionization of Helium as shown in Fig. 2.
Fig. 10. Heavy thermal conductivity λ_{𝔥}, at the secondorder of Laguerre Sonine approximation, for case A and case B, for the heliumhydrogen mixture S_{1} as a function of temperature. 

Open with DEXTER 
Figure 11 shows the components of the total heat flux (Eq. (33)) as a function of the temperature, for the isotropic case A. It is clear that the reactive thermal conductivity is higher than the other components for certain ranges of temperature between 2200 K and 4300 K, where hydrogen dissociation occurs, and for temperature higher than 10 000 K, where hydrogen ionization takes place. The heavy thermal conductivity λ_{𝔥} is the second term which dominates the total heat flux, and is higher than for a range of temperature between 4200 K and 10 000 K. The results here obtained are consistent with those of Scoggins et al. (2016).
Fig. 11. Component of the total heat flux (33) as a function of the temperature for the isotropic case A for the heliumhydrogen mixture S_{1}. 

Open with DEXTER 
Figure 12 and 13 show the parallel component of each term of the electron heavyparticle transport coefficients , as a function of the temperature, for the case A and case B. As before, each term of the electron heavyparticle transport tensor , is proportional to the mole fraction.
Fig. 12. Parallel component of the electron/heavyparticle transport coefficient , for the isotropic case A, for the heliumhydrogen mixture S_{1}. 

Open with DEXTER 
Fig. 13. Parallel component of the electron/heavyparticle transport coefficient , for the anisotropic case B, for the heliumhydrogen mixture S_{1}. 

Open with DEXTER 
5.2. Transport properties in a pore at the photosphere of the Sun
As done in the previous section, the transport coefficients of the previous heliumhydrogen mixture are computed for the conditions found in the upper layer of the solar convective zone from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). The simulation results are obtained for the computational domain of 6.4 × 6.4 × 5.5 Mm^{3} with the grid sizes: 50 × 50 × 43 km^{3}, 25 × 25 × 21.7 km^{3} and 12.5 × 12.5 × 11 km^{3} (128^{2} × 127, 256^{2} × 253 and 512^{2} × 505 mesh points). The domain includes an upper 5 Mmdeep layer of the convective zone and the solar chromosphere. In this section, the results and quantities are obtained from simulations of Kitiashvili et al. (2010), Wray et al. (2015, 2018) via a singlefluid model. This is a postprocessed calculation for quantities that belong to the multicomponent model presented. Since the electron heavyparticle collision frequency is high in the conditions chosen, a thermal equilibrium case T_{𝔢} = T_{𝔥} has been considered for the computation of the transport coefficients. For the sake of clarity, only results from a slice at a constant geometrical height z = −0.5 Mm, in the lower photosphere has been presented.
Figures 14–16 show snapshots of the distribution of the plasma beta parameter β_{p}, temperature T, and total mass density, respectively. As it can be seen, the temperature is varying from 4000 K to 6500 K, the plasma beta parameter is varying on a large range of magnitude, from weakly to stronglymagnetized. In the snapshot of the simulation, a characteristic granulation pattern with the relatively hot (T > 5500 K) and less dense upflowing weaklymagnetized plasma in the middle of the granular cells can be observed. In addition, the lower temperature (T < 4500 K) and higher density downflowing strongly magnetized plasma at the intergranulation boundaries can be perceived (red lines of granulation). A strongly magnetized cold plasma can be seen in the middle of the snapshot.
Fig. 14. Plasma beta coefficient β_{p} distribution from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 15. Temperature (K) distribution from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 16. Total mass density (kg m^{−3}) distribution from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Figures 17–20 present the distribution of the heavyparticle heat flux λ_{𝔥}∂_{x}T and the ratios , /λ_{𝔥}, and , respectively. Figure 18 shows that the electron thermal conductivity tensor is almost isotropic everywhere, except in the middle of the snapshot where . In Fig. 19, the peak of reactive thermal conductivity is attributed to the dissociation reaction of hydrogen. Figure 20 shows that the electron thermal conductivity is small compared to the heavy thermal conductivity. These results are related to the results from Figs. 1 and 2 that show that the mole fraction of electrons is very small compared to the mole fraction of heavy particles in that range of temperature between 4000 K to 6500 K.
Fig. 17. Distribution of the heavy particle heat flux λ_{𝔥}∂_{x}T, for the heliumhydrogen mixture S_{1} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 18. Ratio distribution, computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 19. Ratio /λ_{𝔥} distribution, computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 20. Ratio distribution, computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{2} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
5.3. Components of the generalized Ohm’s law in a pore at the Sun photosphere
As in Sect. 5.2, we compute the components of the generalized Ohm’s law from Eq. (25) using a heliumhydrogen mixture, from the simulation by Kitiashvili et al. (2010). According to the result found in Fig. 18, we assume an isotropic distribution of the transport properties. Figures 21–25 show the distribution of the resistive term, the electron battery term, the heavyparticle battery term and the Soret terms for electron and heavy particles respectively.
Fig. 21. Distribution of the resistive term (first term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 22. Distribution of the electron battery term (second term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 23. Distribution of the heavy battery term (third term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 24. Distribution of the Soret/Dufour term (fourth term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
Fig. 25. Distribution of the Soret/Dufour term (last term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER 
In the results obtained, the resistive term appears to be the higher term in the generalized Ohm’s law in both inside (magnitude 10^{2}) and outside the pore (magnitude 1 − 10, not shown in Fig. 21). However, outside the pore, the battery term for heavy particles appears to be the second higher term. All the other terms are negligible inside and outside the pore (magnitude 10^{−7}). In addition, the Soret and battery terms for electrons are negligible compared to the other terms of the generalized Ohm’s law. Indeed, this is due to the mole fraction of electrons which is very small compared to heavy particles under these conditions. These results are coherent with the mole fraction distribution presented in Figs. 1–3.
6. Conclusion
The present model, derived from the kinetic theory of Graille et al. (2009), is neither a singlefluid MHD nor a multifluid model. It is an intermediary model between the two that is a driftdiffusion multicomponent model. Although the multifluid model and the multicomponent model are catching the same physics, their main difference is the scaling that has been used in the generalized ChapmanEnskog expansion. Their scaling leads to a thermal nonequilibrium multicomponent model with one momentum equation, where the electrons and each heavy particle diffuse in the hydrodynamic heavyparticle reference frame. These developments lead to an extended range of validity for partially and fully ionized plasma, as well as weakly, strongly and nonmagnetized plasmas, and for a general multicomponent mixtures, which can be applied to conditions of the Sun’s atmophere. From the numerical point of view, the multicomponent model is prone to less stiffness compared to the conventional multifluid models, where source terms and coupling terms are considered in the momentum equations and energy equations of each particle. On the other hand, the singlefluid MHD model is shown to be a thermal equilibrium model with no relaxation term and decoupling between the particles. However, in this framework, this model shows less numerical stiffness in absence of source or relaxation terms. Fluidmodels are valid only in cases where the plasma is in collisional regimes leading to small values of the Knudsen number. Thus, only perturbations of Maxwellian distribution functions are considered here.
From the set of governing equations, a generalized Ohm’s law has been derived. A general expression of the resistive term, as well as the battery term, has been obtained for a general multicomponent plasma. This general expression of the electric field can be simplified in a fully ionized plasma case (See Appendix B).
General conditions largely representative of the lower solar atmosphere have been chosen in order to compute all the transport properties for a heliumhydrogen mixture S_{1}. The latter are obtained by solving transport systems which are presented in Appendix C. A spectral Galerkin method based on a thirdorder LaguerreSonine polynomials approximation has been used. For convenience, we have implemented the transport model described in the paper in Mutation++, an opensource library^{1}. In order to validate our model, a comparison with the model of Braginskii (1965) has been performed in the case of a fully ionized plasma S_{2}. Both models are derived from the kinetic theory based on a ChapmanEnskog method. Differences are observed in 1the structure of the governing equations and 2 the nature of the collision operators. While the heavy transport properties are anisotropic in Braginskii (1965), in our model, they remain isotropic. Nevertheless, under the chosen range of conditions, both behave as isotropic. In Braginskii (1965), the corresponding series of LaguerreSonine polynomials are truncated at the second order, whereas a thirdorder accuracy is reached in our model. Good agreement has been obtained for the considered fully ionized S_{2} mixture in the chosen conditions.
Finally, using the Mutation++ library, the method allows us to compute all the transport properties for a partially ionized plasma for a given mixture. The obtained results strongly depend on the mole fraction between the species of the mixture. We have been able to identify the behavior of the transport coefficients related to the chemistry of the species in the partially ionized mixture S_{1}. Considering the complexity of the multicomponent model and the new terms that appear, postprocessing calculations were performed based on the results of a pore simulation for a partially ionized plasma of a mixture of Hydrogen and Helium S_{1} in thermochemical equilibrium in the highly turbulent upper layer of the solar convective zone by Kitiashvili et al. (2010). These results allow us to both understand the effect and magnitude of these terms under the conditions we chose, especially the terms of the generalized Ohm’s law. Generally, the plasma being weakly ionized under these conditions, all the transport fluxes related to the electrons are negligible compared to those from heavy particles. In particular, we show a new formulation of the resistivity which depends on the multicomponent transport coefficients taking into account all possible interactions in the mixture. It appears to be much more general than standard resistivity such as ambipolar resistivity or Spitzer resistivity. Under conditions where the plasma is dominated by electromagnetic forces, which is the case in the center of the pore in Sect. 5, we find qualitatively that this term dominates the dynamics of the weakly ionized heliumhydrogen plasma. On the other hand, at the interboundary granulations (or more generally in the case of quiet sun photosphere conditions) when the plasma is weakly ionized and thermal pressure dominated, the new terms such as the Soret/Dufour effects (related to the temperature or pressure gradients) are dominating the plasma dynamics. In such conditions, these new terms allow magnetic reconnection process to occur. Indeed, according to Eq. (25), the parallel component of these terms may produce an electric field parallel to the magnetic field which may result in reconnection. These terms do not appear in the classical multifluid or singlefluid model (see Leake et al. 2013) since only a resistive term is generally considered, which allows magnetic reconnection to occur only in the magnetic pressure dominated plasma regime. Preliminary results have been obtained by Wargnier et al. (2018) for a fully ionized plasma in thermal nonequilibrium in the framework of a magnetic reconnection event.
The presented transport systems have been implemented into MUTATION++ Library. The GitHub link repository is https://github.com/mutationpp/Mutationpp
Acknowledgments
The research of Q. Wargnier is funded by an Idex ParisSaclay interdisciplinary IDI PhD grant, and relies on the the support of NASA Ames Research Center (ARC), Advanced Supercomputing Division, von Karman Institute for Fluid Dynamics, CMAP – Initiative HPC@Maths from Ecole Polytechnique and of Ecole doctorale de Mathématiques Hadamard. T. Magin is supported by a Jean d’Alembert chair of University ParisSaclay. Part of this work was conducted during the 2018 NASA Summer Program at ARC. We would like to thank I.N. Kitiashvili for providing us with the data in order to evaluate the transport properties in Sect. 5.2. A. Alvarez Laguna and J. B. Scoggins were each funded by a postdoctoral fellowship from the Fondation Mathematique Jacques Hadamard (FMJH), LabEx Mathématique Hadamard  ANR11LABX0056.
References
 Alvarez Laguna, A., Lani, A., Deconinck, H., Mansour, N. N., & Poedts, S. 2016, J. Comput. Phys., 318, 252 [NASA ADS] [CrossRef] [Google Scholar]
 Alvarez Laguna, A., Ozak, N., Lani, A., et al. 2018, J. Phys. Conf. Ser., 1031, 012015 [CrossRef] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Balescu, R. 1988, Transport Processes in Plasmas, Vol. 1: Classical Transport; Vol. 2: Neoclassical Transport (NorthHolland Press) [Google Scholar]
 Braginskii, S. I. 1965, Rev. Plasma Phys., 1, 205 [Google Scholar]
 Bruno, D., Capitelli, M., Catalfamo, C., & Giordano, D. 2011, Phys. Plasmas, 18, 012308 [NASA ADS] [CrossRef] [Google Scholar]
 Bruno, D., Catalfamo, C., Capitelli, M., et al. 2010, Phys. Plasmas, 17, 112315 [NASA ADS] [CrossRef] [Google Scholar]
 Capitelli, M., Bruno, D., & Laricchiuta, A. 2013, Fundamental Aspects of Plasma Chemical Physics I (Springer) [CrossRef] [Google Scholar]
 Carlsson, M., & Stein, R. F. 1995, ApJ, 440, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Devoto, R. S. 1969, J. Plasma Phys., 2 [Google Scholar]
 Ferziger, J. H., & Kaper, H. G. 1973, Am. J. Phys., 41, 601 [NASA ADS] [CrossRef] [Google Scholar]
 Fröhlich, C., & Lean, J. 2004, A&ARv., 12, 273 [NASA ADS] [CrossRef] [Google Scholar]
 Giovangigli, V., & Graille, B. 2003, Phys. A, 327, 313 [CrossRef] [Google Scholar]
 Giovangigli, V., & Graille, B. 2009, Phys. A, 42, 025503 [NASA ADS] [CrossRef] [Google Scholar]
 Graille, B., Magin, T. E., & Massot, M. 2009, Math. Models Methods Appl. Sci., 19, 527 [CrossRef] [Google Scholar]
 Hartlep, T., Busse, F. H., Hurlburt, N. E., & Kosovichev, A. G. 2012, MNRAS, 419, 2325 [NASA ADS] [CrossRef] [Google Scholar]
 Khomenko, E. 2017, Plasma Phys. Controlled Fusion, 59, 014038 [NASA ADS] [CrossRef] [Google Scholar]
 Khomenko, E., Collados, M., Diaz, A., & Vitas, N. 2014, Phys. Plasmas, 21, 092901 [NASA ADS] [CrossRef] [Google Scholar]
 Kitiashvili, I. N., Couvidat, S., & Lagg, A. 2015, ApJ, 808, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Kitiashvili, I. N., Kosovichev, A. G., Wray, A. A., & Mansour, N. N. 2010, ApJ, 719, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Kolesnikov, A. 2003, 41st Aerospace Sciences Meeting and Exhibit [Google Scholar]
 Leake, J. E., Lukin, V. S., & Linton, M. G. 2013, Phys. Plasmas, 20, 061202 [NASA ADS] [CrossRef] [Google Scholar]
 Magin, T. E., & Degrez, G. 2004a, J. Comput. Phys., 198, 424 [NASA ADS] [CrossRef] [Google Scholar]
 Magin, T., & Degrez, G. 2004b, Phys. Rev. E, 70, 046412 [NASA ADS] [CrossRef] [Google Scholar]
 Magin, T. E., Graille, B., & Massot, M. 2009, Center Turbul. Res. Annu. Res. Briefs, 71 [Google Scholar]
 Martinez Sykora, J., De Pontieu, B., Hansteen, V. H., & Carlsson, M. 2015, SH, 31B [Google Scholar]
 Ni, L., Lukin, V. S., Murphy, N. A., & Lin, J. 2018, Phys. Plasmas, 25, 042903 [NASA ADS] [CrossRef] [Google Scholar]
 Russell, H. N. 1929, ApJ, 70, 11 [NASA ADS] [CrossRef] [Google Scholar]
 Scoggins, J. B., Knisely, C. P., & Magin, T. E. 2016, AIP Conf. Proc., 1786, 130002 [CrossRef] [Google Scholar]
 Scoggins, J. B., & Magin, T. E. 2014, 11th AIAA/ASME Joint Thermophysics and Heat Transfer Conference [Google Scholar]
 Scoggins, J. B., & Magin, T. E. 2015, Combust. Flame, 162, 4514 [CrossRef] [Google Scholar]
 Solanki, S. 2003, A&ARv., 11, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Spitzer, L. 1963, Am. J. Phys., 31, 890 [NASA ADS] [CrossRef] [Google Scholar]
 Tirsky, G. A. 1993, Annu. Rev. Fluid Mech., 25, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 [Google Scholar]
 Wargnier, Q., Laguna, A. A., Kestener, P., et al. 2018, J. Phys.: Conf. Ser., 1125, 012021 [CrossRef] [Google Scholar]
 Wiegelmann, T., Thalmann, J. K., & Solanki, S. K. 2014, A&ARv., 22, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Woods, L. C. 1995, An Introduction to the Kinetic Theory of Gases and Magnetoplasmas (Cambridge University Press), 286, 406 [Google Scholar]
 Wray, A. A., Bensassi, K., Kitiashvili, I. N., Mansour, N., & Kosovichev, A. G. 2015, ArXiv eprints [arXiv:1507.07999] [Google Scholar]
 Wray, A. A., Bensassi, K., Kitiashvili, I. N., Mansour, N., & Kosovichev, A. G. 2018, in Realistic simulations of Stellar Radiative MHD, ed. E. B. J.P. Rozelot (EDP Sciences), 39 [Google Scholar]
 Zhdanov, V. 2002, Transport Processes in Multicomponent Plasma (Taylor & Francis Group), 44 [Google Scholar]
Appendix A: Nondimensional Boltzmann equations
In Graille et al. (2009), the nondimensional Boltzmann equations for electron and heavy particles is
and for each heavy species,
Where the collision operators are defined as
Where C_{i}, i ∈ H and C_{𝔢} are the peculiar velocities for heavy species and electron respectively.
Appendix B: Ohm’s law for a fully ionized plasma S_{1}
In a fully ionized plasma, the definition of the current density is
where the electron current density is defined by Eq. (4). It should be noted that in the particular case of a fully ionized plasma, no diffusion velocities of heavy species are considered, because the presented model is in the reference frame of heavy particles.
Then, similarly as the general case, the total current density can be expressed in function of fluxes, as follows,
It should be noted that the multicomponent electromagnetic matrices in a fully ionized plasma are
Using the MaxwellAmpere’s law in a nonrelativistic context, the expression of the electric field as a function of fluxes can be obtained
where the electron resistivity tensor is .
Appendix C: Transport systems
C.1. Anisotropic transport systems
Anisotropic transport coefficients are computed in terms of solutions to complex linear systems. For μ ∈ {D_{𝔢}, (D_{i})_{i ∈ H}}, transport systems read for Sonine polynomial order ξ ≥ 1
where P_{ξ} = {0, …, ξ − 1}. Likewise, for λ_{𝔢}, we have
where P_{ξ1} = {1, …, ξ − 1} for ξ ≥ 2, the matrices and are given up to third order as
and the righthand sides are given as
where are the socalled reduced collision integrals.
Parallel components of the anisotropic transport coefficients are related to the solutions of the real linear systems above,
Like wise, the perpendicular and transverse components are found using the solutions of the complex linear systems above, by replacing the left hand side of these equations with μ^{⊥}(ξ)+iμ^{⊙}(ξ) for μ = {D_{𝔢}, λ_{𝔢}, χ_{𝔢}, α_{𝔢i}, χ_{𝔢i}} and replacing (1) with (2).
C.2. Isotropic transport systems
The usual collision integral ratios are defined as
In addition, we define binary diffusion coefficients and pure species viscosities for heavy species as
First order heavyparticle multicomponent diffusion coefficients are found as the solution of 𝔫_{𝔥} linear systems of the form
where the matrix is given as
where the hat on mole and mass fractions denotes that they are only over heavy species, ie: for i ∈ H.
Shear viscosity and heavyparticle thermal conductivity are found to be the solution of a linear system of the form
where the matrices and are given as
and
Likewise, the heavy thermal diffusion ratios can be computed as
where
All Tables
Time scales hierarchy and macroscopic equations derived using the ChapmanEnskog method, Graille et al. (2009).
All Figures
Fig. 1. Mole fraction of the heliumhydrogen mixture S_{1}, for P = 10^{4} Pa (case A), as a function of temperature. 

Open with DEXTER  
In the text 
Fig. 2. Mole fraction of the heliumhydrogen mixture S_{1}, for P = 1 Pa (case A), as a function of temperature. 

Open with DEXTER  
In the text 
Fig. 3. Ionization degree the heliumhydrogen mixture of S_{1}, for case A (full line) and case B (dashed line) as function of the temperature. 

Open with DEXTER  
In the text 
Fig. 4. Parallel component of the electron thermal conductivity tensor for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER  
In the text 
Fig. 5. Perpendicular component of the electron thermal conductivity tensor for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER  
In the text 
Fig. 6. Heavy thermal conductivity λ_{𝔥} for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER  
In the text 
Fig. 7. Heavy particle viscosity η_{𝔥} for a fully ionized plasma S_{2}, as function of temperature. Dashed lines and full lines correspond to the transport coefficient from the model of Braginskii (1965), and from Graille et al. (2009) respectively. Bold lines correspond to the case A, the other lines correspond to the case B. 

Open with DEXTER  
In the text 
Fig. 8. Components of the electron thermal conductivity tensor for the isotropic case A, at the thirdorder Laguerre Sonine polynomials, for the heliumhydrogen mixture S_{1} as a function of temperature. 

Open with DEXTER  
In the text 
Fig. 9. Components of the electron thermal conductivity tensor for the anisotropic case B, at the thirdorder Laguerre Sonine polynomials, for the heliumhydrogen mixture S_{1} as a function of the temperature. 

Open with DEXTER  
In the text 
Fig. 10. Heavy thermal conductivity λ_{𝔥}, at the secondorder of Laguerre Sonine approximation, for case A and case B, for the heliumhydrogen mixture S_{1} as a function of temperature. 

Open with DEXTER  
In the text 
Fig. 11. Component of the total heat flux (33) as a function of the temperature for the isotropic case A for the heliumhydrogen mixture S_{1}. 

Open with DEXTER  
In the text 
Fig. 12. Parallel component of the electron/heavyparticle transport coefficient , for the isotropic case A, for the heliumhydrogen mixture S_{1}. 

Open with DEXTER  
In the text 
Fig. 13. Parallel component of the electron/heavyparticle transport coefficient , for the anisotropic case B, for the heliumhydrogen mixture S_{1}. 

Open with DEXTER  
In the text 
Fig. 14. Plasma beta coefficient β_{p} distribution from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 15. Temperature (K) distribution from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 16. Total mass density (kg m^{−3}) distribution from the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 17. Distribution of the heavy particle heat flux λ_{𝔥}∂_{x}T, for the heliumhydrogen mixture S_{1} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 18. Ratio distribution, computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 19. Ratio /λ_{𝔥} distribution, computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 20. Ratio distribution, computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{2} based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 21. Distribution of the resistive term (first term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 22. Distribution of the electron battery term (second term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 23. Distribution of the heavy battery term (third term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 24. Distribution of the Soret/Dufour term (fourth term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Fig. 25. Distribution of the Soret/Dufour term (last term of Eq. (25)), computed at the third order of the LaguerreSonine polynomials approximation, for the heliumhydrogen mixture S_{1}, based on the results of the radiative 3D MHD simulations of a pore by Kitiashvili et al. (2010). 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.