Numerical simulation of the plasma of an electron cyclotron resonance ion source

Authors:
Alain Girard;Christian Lécot;Konstantin Serebrennikov
Affiliations:
Département de Recherche Fondamentale sur la Matière Condensée, SBT, CEA Grenoble, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France;Laboratoire de Mathématiques, Université de Savoia, Campus scientifique, 73376 Le Bourget-du-Lac Cedex, France;Département de Recherche Fondamentale sur la Matière Condensée, SBT, CEA Grenoble, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France
Venue:
Journal of Computational Physics
Year:
2003

Citing 7
Cited 0

Conservative finite-difference schemes for the Fokker-Planck equations not violating the law of an increasing entropy

Journal of Computational Physics
Random number generation and quasi-Monte Carlo methods

Random number generation and quasi-Monte Carlo methods
A control volume method to solve an elliptic equation on a two-dimensional irregular mesh

Computer Methods in Applied Mechanics and Engineering
Fast algorithms for numerical, conservative, and entropy approximations of the Fokker-Planck-Landau equation

Journal of Computational Physics
Conservative and entropy decaying numerical scheme for the isotropic Fokker-Planck-Landau equation

Journal of Computational Physics
An implicit energy-conservative 2D Fokker-Planck algorithm: I. difference scheme

Journal of Computational Physics
An implicit energy-conservative 2D Fokker-Planck algorithm: II. Jacobian-free Newton—Krylov solver

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

We propose a kinetic description of the particles in an electron cyclotron resonance ion source (ECRIS) for highly charged ion production. The electron distribution function (EDF) is solution of a Fokker-Planck equation in velocity space. A new finite-volume scheme for the computation of the EDF is presented. The system is azimuthally symmetric. A two-dimensional, structured grid of quadrangles is used to describe the geometry in spherical polar coordinates (v, µ). The temporal stiffness problems are addressed with a semi-implicit method of time integration. Evaluation of the coefficients of the collision operator requires derivatives ef the Rosenbluth potentials in every time step. The computation is accomplished through parallel implementation. Examples of the application of the algorithm are given and the role of the different design parameters is studied.