A 10-moment fluid numerical solver of plasma with sheaths in a Hall Effect Thruster

Joncquieres, Valentin; Péchereau, François; Laguna, Alejandro Alvarez; Bourdon, Anne; Vermorel, Olivier; Cuenot, Bénédicte

doi:10.2514/6.2018-4905

Cited by 15 publications

(19 citation statements)

References 39 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The work was an extension of previous attempts by Mikellides [227] to incorporate the ECDI in Hall2De simulations based on the hypothesis that the instability excites IAT which, in turn, enhances the effective collision frequency. Other larger scale hybrid approaches are being pursued, such as that by Joncquieres, et al [228] who are developing a 3-D unstructured massively parallel Particle-in-Cell/fluid solver. In this work, the fluid part for electrons and ions is based on a 10-moment model, while the kinetic simulations are used as a reference solver for the sheath boundary conditions.…”

Section: Fluid and Hybrid (Fluid-kinetic) Modeling Of Eb Dischargesmentioning

confidence: 99%

Physics of E × B discharges relevant to plasma propulsion and similar technologies

et al. 2020

Self Cite

View full text Add to dashboard Cite

This paper provides perspectives on recent progress in the understanding of the physics of devices where the external magnetic field is applied perpendicularly to the discharge current. This configuration generates a strong electric field, which acts to accelerates ions. The many applications of this set up include generation of thrust for spacecraft propulsion and the separation of species in plasma mass separation devices. These "E×B" plasmas are subject to plasma-wall interaction effects as well as various micro and macro instabilities, and in many devices, we observe the emergence of anomalous transport. This perspective presents the current understanding of the physics of these phenomena, state-of-the-art computational results, identifies critical questions, and suggests directions for future research.

show abstract

Section: Fluid and Hybrid (Fluid-kinetic) Modeling Of Eb Dischargesmentioning

confidence: 99%

Physics of E × B discharges relevant to plasma propulsion and similar technologies

et al. 2020

Self Cite

View full text Add to dashboard Cite

show abstract

“…Multifluid-Maxwell models have been rapidly gaining interest recently in laboratory, space, and astrophysical plasma phyiscs modeling [6,28,17,29,16,30,31,32,33,34,35,36,37]. In this paper, we have developed a locally implicit schemes to integrate the source terms for such systems.…”

Section: Discussionmentioning

confidence: 99%

“…Now we find two real, orthogonal solutions for λ = iω of different polarizations as in Eq. (34). In total, we recover all 2 × (N + 1) real, orthogonal solutions to the 2 (N + 1) order Eq.…”

Section: B33 Real Solution Basesmentioning

confidence: 99%

Exact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems

Wang,

Hakim,

et al. 2019

Preprint

View full text Add to dashboard Cite

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3S + 3) × (3S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. We present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

show abstract

“…If the Debye length is resolved, then the convective condition is sufficient to fulfil the condition ∆t < ω −1 pe . As shown in [38,70], when this scheme does not spatially resolve the Debye length, the simulation leads to large spurious charge separation errors that can excite plasma modes, leading to an erroneous solution. Similarly, the truncation error of the upwind discretization of the fluxes of the electrons leads to a large error in the flux of the electrons.…”

Section: Standard Upwind Finite Volume Discretizationmentioning

confidence: 99%

An asymptotic preserving well-balanced scheme for the isothermal fluid equations in low-temperature plasma applications

Laguna,

Pichard,

Magin

et al. 2019

Preprint

Self Cite

View full text Add to dashboard Cite

We present a novel numerical scheme for the efficient and accurate solution of the isothermal two-fluid (electron + ion) equations coupled to Poisson's equation for low-temperature plasmas. The model considers electrons and ions as separate fluids, comprising the electron inertia and charge separation. The discretization of this system with standard explicit schemes is constrained by very restrictive time steps and cell sizes related to the resolution of the Debye length, electron plasma frequency, and electron sound waves. Both sheath and electron inertia are fundamental to fully explain the physics in low-pressure and low-temperature plasmas. However, most of the phenomena of interest for fluid models occur at speeds much slower than the electron thermal speed and are quasi-neutral, except in small charged regions. A numerical method that is able to simulate efficiently and accurately all these regimes is a challenge due to the multiscale character of the problem. In this work, we present a scheme based on the Lagrange-projection operator splitting that preserves the asymptotic regime where the plasma is quasi-neutral with massless electrons. As a result, the quasi-neutral regime is treated without the need of an implicit solver nor the resolution of the Debye length and electron plasma frequency. Additionally, the scheme proves to accurately represent the dynamics of the electrons both at low speeds and when the electron speed is comparable to the thermal speed. In addition, a well-balanced treatment of the ion source terms is proposed in order to tackle problems where the ion temperature is very low compared to the electron temperature. The scheme significantly improves the accuracy both in the quasi-neutral limit and in the presence of plasma sheaths when the Debye length is resolved. In order to assess the performance of the scheme in low-temperature plasmas conditions, we propose two specifically designed testcases: a quasi-neutral two-stream periodic perturbation with analytical solution and a low-temperature discharge that includes sheaths. The numerical strategy, its accuracy, and computational efficiency are assessed on these two discriminating configurations.

show abstract

A 10-moment fluid numerical solver of plasma with sheaths in a Hall Effect Thruster

Cited by 15 publications

References 39 publications

Physics of E × B discharges relevant to plasma propulsion and similar technologies

Physics of E × B discharges relevant to plasma propulsion and similar technologies

Exact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems

An asymptotic preserving well-balanced scheme for the isothermal fluid equations in low-temperature plasma applications

Contact Info

Product

Resources

About