skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An asymptotic preserving scheme for the two-fluid Euler-Poisson model in the quasineutral limit

Abstract

This paper deals with the modeling of a plasma in the quasineutral limit using the two-fluid Euler-Poisson 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 one-dimensional test-cases. 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 two-dimensional problem of a plasma expansion in vacuum between two electrodes.

Authors:
 [1];  [2];  [3]
  1. 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). E-mail: crispel@mip.ups-tlse.fr
  2. MIP, Univ. P. Sabatier, 118 rte de Narbonne, 31062 Toulouse cedex 4 (France). E-mail: degond@mip.ups-tlse.fr
  3. MIP, Univ. P. Sabatier, 118 rte de Narbonne, 31062 Toulouse cedex 4 (France). E-mail: mhvignal@mip.ups-tlse.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: S0021-9991(06)00433-5; 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; ONE-DIMENSIONAL CALCULATIONS; PLASMA; PLASMA EXPANSION; POISSON EQUATION; TWO-DIMENSIONAL CALCULATIONS

Citation Formats

Crispel, Pierre, Degond, Pierre, and Vignal, Marie-Helene. An asymptotic preserving scheme for the two-fluid Euler-Poisson model in the quasineutral limit. United States: N. p., 2007. Web. doi:10.1016/j.jcp.2006.09.004.
Crispel, Pierre, Degond, Pierre, & Vignal, Marie-Helene. An asymptotic preserving scheme for the two-fluid Euler-Poisson model in the quasineutral limit. United States. doi:10.1016/j.jcp.2006.09.004.
Crispel, Pierre, Degond, Pierre, and Vignal, Marie-Helene. Tue . "An asymptotic preserving scheme for the two-fluid Euler-Poisson 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 two-fluid Euler-Poisson model in the quasineutral limit},
author = {Crispel, Pierre and Degond, Pierre and Vignal, Marie-Helene},
abstractNote = {This paper deals with the modeling of a plasma in the quasineutral limit using the two-fluid Euler-Poisson 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 one-dimensional test-cases. 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 two-dimensional 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 kinetic-fluid model for disperse two-phase flows with uncertainty. We propose a stochastic asymptotic-preserving (s-AP) scheme in the generalized polynomial chaos stochastic Galerkin (gPC-sG) framework, which allows the efficient computation of the problem in both kinetic and hydrodynamic regimes. The s-AP 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 » that arises in the inversion of the Fokker–Planck operator is accelerated by a spectrally accurate splitting method. Numerical examples illustrate the s-AP property and the efficiency of the gPC-sG method in various asymptotic regimes.« less
  • In applications such as astrophysics and inertial confinement fusion, there are many three-dimensional cylindrical-symmetric multi-material problems which are usually simulated by Lagrangian schemes in the two-dimensional 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 cell-centered Lagrangian scheme for solving compressible Euler equations in cylindrical coordinates, basedmore » on the control volume discretizations, which is designed to have uniformly second order accuracy and capability to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. The scheme maintains several good properties such as conservation for mass, momentum and total energy, and the geometric conservation law. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of accuracy, symmetry, non-oscillation and robustness. The advantage of higher order accuracy is demonstrated in these examples.« less
  • 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 » schemes, this method is based on a standard finite volume approach, it does neither use any decomposition of the solution, nor staggered grids. Several numerical tests demonstrate the properties of the scheme.« less
  • We design an asymptotic-preserving scheme for the semiconductor Boltzmann equation which leads to an energy-transport system for electron mass and energy as mean free path goes to zero. As opposed to the classical drift-diffusion limit where the stiff collisions are all in one scale, new difficulties arise in the two-scale 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 » Formal asymptotic analysis and numerical results confirm the efficiency and accuracy of our scheme.« less
  • 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 » scale” reformulation of the initial equation, with the introduction of an additional periodic variable.« less