ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates
arXiv:1810.08433 · doi:10.1016/j.cpc.2018.10.026
Abstract
Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsim-CYL, to be used for axially symmetric problems. The main motivation for the development of the cylindrical version is to greatly improve the computational speed by utilizing cylindrical symmetry. The ECsim-CYL discretizes the field equations in two-dimensional cylindrical coordinates using the finite volume method . For the particle mover, it uses a modification of ECsim's mover for cylindrical coordinates by keeping track of all three components of velocity vectors, while only keeping radial and axial coordinates of particle positions. In this paper, we describe the details of the algorithm used in the ECsim-CYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasmas inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. The ECsim-CYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.
arXiv commentSubmitted to Computer Phys. Comm
References in corpus (7)
- An Energy- and Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm
- The Energy Conserving Particle-in-Cell Method
- Exactly Energy Conserving Semi-Implicit Particle in Cell Formulation
- CPIC: a Curvilinear Particle-In-Cell code for plasma-material interaction studies
- Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations
- Multiple scale kinetic simulations with the energy conserving semi implicit particle in cell (PIC) method
- Performance analysis and implementation details of the Energy Conserving Semi-Implicit Method code (ECsim)
Cited by in corpus (3)
- Advances in the Implementation of the Exactly Energy Conserving Semi-Implicit (ECsim) Particle in Cell Method
- Fully Kinetic Simulations of Proton-Beam-Driven Instabilities from Parker Solar Probe Observations
- Particle-In-Cell simulations of the parallel proton firehose instability influenced by the electron temperature anisotropy in solar wind conditions