Structure-preserving vs. standard particle-in-cell methods: The case of an electron hybrid model

•Structure-preserving discretization leads to noncanonical Hamiltonian system in time.•No spurious modes compared to standard numerical methods.•Better long-term energy conservation compared to standard numerical methods. We applied two numerical methods both belonging to the class of finite element...

Full description

Saved in:
Bibliographic Details
Published in:Journal of computational physics 2020-02, Vol.402, p.109108, Article 109108
Main Authors: Holderied, Florian, Possanner, Stefan, Ratnani, Ahmed, Wang, Xin
Format: Article
Language:English
Subjects:
Algorithms
B-splines
Computational physics
Differential equations
Electrons
Finite element exterior calculus
Finite elements
Mathematical models
Numerical methods
Ordinary differential equations
Particle in cell technique
Particle-in-cell
Phase velocity
Plasma simulation
Structure-preserving
Three dimensional models
Time integration
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:•Structure-preserving discretization leads to noncanonical Hamiltonian system in time.•No spurious modes compared to standard numerical methods.•Better long-term energy conservation compared to standard numerical methods. We applied two numerical methods both belonging to the class of finite element particle-in-cell methods to a four-dimensional (one dimension in real space and three dimensions in velocity space) hybrid plasma model for electrons in a stationary, neutralizing background of ions. Here, the term hybrid means that (energetic) electrons with velocities close to the phase velocities of the model's characteristic waves are treated kinetically, whereas electrons that are much slower than the phase velocity are treated with fluid equations. The two developed numerical schemes are based on standard finite elements on the one hand and on structure-preserving geometric finite elements on the other hand. We tested and compared the schemes in the linear and in the nonlinear stage. We show that the structure-preserving algorithm leads to better results in both stages. This can be related to the fact that the spatial discretization results in a large system of ordinary differential equations that exhibits a noncanonical Hamiltonian structure. To such systems special time integration schemes with good conservation properties can be applied.
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2019.109108