An asymptotic preserving scheme for the twofluid EulerPoisson model in the quasineutral limit
Abstract
This paper deals with the modeling of a plasma in the quasineutral limit using the twofluid EulerPoisson system. In this limit, explicit numerical schemes suffer from severe numerical constraints related to the small Debye length and large plasma frequency. Here, we propose an implicit scheme which reduces to a scheme for the quasineutral Euler model in the quasineutral limit. Such a property is referred to as 'asymptotic preservation'. One of the distinctive features of this scheme is that it has a comparable numerical cost to that of an explicit scheme: simply the Poisson equation is replaced by a different (but formally equivalent) elliptic problem. We present numerical simulations for two different onedimensional testcases. They confirm the expected stability of the scheme in the quasineutral limit. They also show that this scheme has some accuracy problems in the limit of small electron to ion mass ratio in reproducing the correct electron velocity. But this problem is already present in the results of the classical algorithm. Numerical simulations are also performed for a twodimensional problem of a plasma expansion in vacuum between two electrodes.
 Authors:
 CNES Centre de Toulouse, 18 av. Ed. Belin, 31401 Toulouse cedex 4 (France) and ONERA Centre de Toulouse, 2 av. Ed. Belin, 31055 Toulouse cedex 4 (France) and MIP, Univ. P. Sabatier, 118 rte de Narbonne, 31062 Toulouse cedex 4 (France). Email: crispel@mip.upstlse.fr
 MIP, Univ. P. Sabatier, 118 rte de Narbonne, 31062 Toulouse cedex 4 (France). Email: degond@mip.upstlse.fr
 MIP, Univ. P. Sabatier, 118 rte de Narbonne, 31062 Toulouse cedex 4 (France). Email: mhvignal@mip.upstlse.fr
 Publication Date:
 OSTI Identifier:
 20991572
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 223; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2006.09.004; PII: S00219991(06)004335; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; ALGORITHMS; COMPUTERIZED SIMULATION; DEBYE LENGTH; ELECTRONS; FLUIDS; LANGMUIR FREQUENCY; MATHEMATICAL MODELS; ONEDIMENSIONAL CALCULATIONS; PLASMA; PLASMA EXPANSION; POISSON EQUATION; TWODIMENSIONAL CALCULATIONS
Citation Formats
Crispel, Pierre, Degond, Pierre, and Vignal, MarieHelene. An asymptotic preserving scheme for the twofluid EulerPoisson model in the quasineutral limit. United States: N. p., 2007.
Web. doi:10.1016/j.jcp.2006.09.004.
Crispel, Pierre, Degond, Pierre, & Vignal, MarieHelene. An asymptotic preserving scheme for the twofluid EulerPoisson model in the quasineutral limit. United States. doi:10.1016/j.jcp.2006.09.004.
Crispel, Pierre, Degond, Pierre, and Vignal, MarieHelene. Tue .
"An asymptotic preserving scheme for the twofluid EulerPoisson model in the quasineutral limit". United States.
doi:10.1016/j.jcp.2006.09.004.
@article{osti_20991572,
title = {An asymptotic preserving scheme for the twofluid EulerPoisson model in the quasineutral limit},
author = {Crispel, Pierre and Degond, Pierre and Vignal, MarieHelene},
abstractNote = {This paper deals with the modeling of a plasma in the quasineutral limit using the twofluid EulerPoisson system. In this limit, explicit numerical schemes suffer from severe numerical constraints related to the small Debye length and large plasma frequency. Here, we propose an implicit scheme which reduces to a scheme for the quasineutral Euler model in the quasineutral limit. Such a property is referred to as 'asymptotic preservation'. One of the distinctive features of this scheme is that it has a comparable numerical cost to that of an explicit scheme: simply the Poisson equation is replaced by a different (but formally equivalent) elliptic problem. We present numerical simulations for two different onedimensional testcases. They confirm the expected stability of the scheme in the quasineutral limit. They also show that this scheme has some accuracy problems in the limit of small electron to ion mass ratio in reproducing the correct electron velocity. But this problem is already present in the results of the classical algorithm. Numerical simulations are also performed for a twodimensional problem of a plasma expansion in vacuum between two electrodes.},
doi = {10.1016/j.jcp.2006.09.004},
journal = {Journal of Computational Physics},
number = 1,
volume = 223,
place = {United States},
year = {Tue Apr 10 00:00:00 EDT 2007},
month = {Tue Apr 10 00:00:00 EDT 2007}
}

In this paper we consider a kineticfluid model for disperse twophase flows with uncertainty. We propose a stochastic asymptoticpreserving (sAP) scheme in the generalized polynomial chaos stochastic Galerkin (gPCsG) framework, which allows the efficient computation of the problem in both kinetic and hydrodynamic regimes. The sAP property is proved by deriving the equilibrium of the gPC version of the Fokker–Planck operator. The coefficient matrices that arise in a Helmholtz equation and a Poisson equation, essential ingredients of the algorithms, are proved to be positive definite under reasonable and mild assumptions. The computation of the gPC version of a translation operatormore »

Second order symmetrypreserving conservative Lagrangian scheme for compressible Euler equations in twodimensional cylindrical coordinates
In applications such as astrophysics and inertial confinement fusion, there are many threedimensional cylindricalsymmetric multimaterial problems which are usually simulated by Lagrangian schemes in the twodimensional cylindrical coordinates. For this type of simulation, a critical issue for the schemes is to keep spherical symmetry in the cylindrical coordinate system if the original physical problem has this symmetry. In the past decades, several Lagrangian schemes with such symmetry property have been developed, but all of them are only first order accurate. In this paper, we develop a second order cellcentered Lagrangian scheme for solving compressible Euler equations in cylindrical coordinates, basedmore » 
On the asymptotic preserving property of the unified gas kinetic scheme for the diffusion limit of linear kinetic models
The unified gas kinetic scheme (UGKS) of K. Xu et al. (2010) [37], originally developed for multiscale gas dynamics problems, is applied in this paper to a linear kinetic model of radiative transfer theory. While such problems exhibit purely diffusive behavior in the optically thick (or small Knudsen) regime, we prove that UGKS is still asymptotic preserving (AP) in this regime, but for the free transport regime as well. Moreover, this scheme is modified to include a time implicit discretization of the limit diffusion equation, and to correctly capture the solution in case of boundary layers. Contrary to many APmore » 
An asymptoticpreserving scheme for the semiconductor Boltzmann equation toward the energytransport limit
We design an asymptoticpreserving scheme for the semiconductor Boltzmann equation which leads to an energytransport system for electron mass and energy as mean free path goes to zero. As opposed to the classical driftdiffusion limit where the stiff collisions are all in one scale, new difficulties arise in the twoscale stiff collision terms because the simple BGK penalization [15] fails to drive the solution to the correct limit. We propose to set up a spatially dependent threshold on the penalization of the stiffer collision operator such that the evolution of the solution resembles a Hilbert expansion at the continuous level.more » 
Asymptotic Preserving schemes for highly oscillatory Vlasov–Poisson equations
This work is devoted to the numerical simulation of a Vlasov–Poisson model describing a charged particle beam under the action of a rapidly oscillating external field. We construct an Asymptotic Preserving numerical scheme for this kinetic equation in the highly oscillatory limit. This scheme enables to simulate the problem without using any time step refinement technique. Moreover, since our numerical method is not based on the derivation of the simulation of asymptotic models, it works in the regime where the solution does not oscillate rapidly, and in the highly oscillatory regime as well. Our method is based on a “twomore »