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

Abstract

We present a novel numerical scheme for the efficient and accurate solution of the isothermal two-fluid (electron and 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. 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. The scheme is assessed in two specifically designed test-cases: a quasi-neutral two-stream periodic perturbation with analytical solution and a low-temperature discharge that includes sheaths.

References in corpus (9)

• Theory of the Electron Sheath and Presheath
• Continuum Kinetic and Multi-Fluid Simulations of Classical Sheaths
• A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations
• Second order Implicit-Explicit Total Variation Diminishing schemes for the Euler system in the low Mach regime
• A new numerical strategy with space-time adaptivity and error control for multi-scale streamer discharge simulations
• Asymptotic-Preserving Schemes for Fluid Models of Plasmas
• Using the maximum entropy distribution to describe electrons in reconnecting current sheets
• A GPU-enabled implicit Finite Volume solver for the ideal two-fluid plasma model on unstructured grids
• Kinetic Solvers with Adaptive Mesh in Phase Space for Low-Temperature Plasmas

Cited by in corpus (4)

• Benchmarking between fluid and global models for low-pressure oxygen DC glow discharges
• 14-moment maximum-entropy modelling of collisionless ions for Hall thruster discharges
• Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling
• Comparison of high-order moment models for the ion dynamics in a bounded low-temperature plasma