ECsimCYL: Energy Conserving SemiImplicit particle in cell simulation in axially symmetric cylindrical coordinates
Abstract
Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsimCYL, 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 ECsimCYL discretizes the field equations in twodimensional 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 ECsimCYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasma inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. The ECsimCYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.
 Authors:

 Univ. of Leuven (Belgium). Dept. of Mathematics
 Energy Matter Conversion Corporation (EMC2), San Diego, CA (United States)
 Univ. of Leuven (Belgium). Dept. of Mathematics; Royal Observatory of Belgium, Brussels (Belgium). Solar Terrestrial Centre of Excellence (SIDC)
 Publication Date:
 Research Org.:
 Lawrence Berkeley National LaboratoryNational Energy Research Scientific Computing Center (NERSC)
 Sponsoring Org.:
 USDOE Office of Science (SC)
 OSTI Identifier:
 1530608
 DOE Contract Number:
 AC0205CH11231
 Resource Type:
 Journal Article
 Journal Name:
 Computer Physics Communications
 Additional Journal Information:
 Journal Volume: 236; Journal Issue: C; Journal ID: ISSN 00104655
 Country of Publication:
 United States
 Language:
 English
Citation Formats
GonzalezHerrero, Diego, Micera, Alfredo, Boella, Elisabetta, Park, Jaeyoung, and Lapenta, Giovanni. ECsimCYL: Energy Conserving SemiImplicit particle in cell simulation in axially symmetric cylindrical coordinates. United States: N. p., 2019.
Web. doi:10.1016/j.cpc.2018.10.026.
GonzalezHerrero, Diego, Micera, Alfredo, Boella, Elisabetta, Park, Jaeyoung, & Lapenta, Giovanni. ECsimCYL: Energy Conserving SemiImplicit particle in cell simulation in axially symmetric cylindrical coordinates. United States. doi:10.1016/j.cpc.2018.10.026.
GonzalezHerrero, Diego, Micera, Alfredo, Boella, Elisabetta, Park, Jaeyoung, and Lapenta, Giovanni. Fri .
"ECsimCYL: Energy Conserving SemiImplicit particle in cell simulation in axially symmetric cylindrical coordinates". United States. doi:10.1016/j.cpc.2018.10.026.
@article{osti_1530608,
title = {ECsimCYL: Energy Conserving SemiImplicit particle in cell simulation in axially symmetric cylindrical coordinates},
author = {GonzalezHerrero, Diego and Micera, Alfredo and Boella, Elisabetta and Park, Jaeyoung and Lapenta, Giovanni},
abstractNote = {Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsimCYL, 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 ECsimCYL discretizes the field equations in twodimensional 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 ECsimCYL and present a series of tests to validate the accuracy of the code including a wave spectrum in a homogeneous plasma inside a cylindrical waveguide and free expansion of a spherical plasma ball in vacuum. The ECsimCYL retains the stability properties of ECsim and conserves the energy within machine precision, while accurately describing the plasma behavior in the test cases.},
doi = {10.1016/j.cpc.2018.10.026},
journal = {Computer Physics Communications},
issn = {00104655},
number = C,
volume = 236,
place = {United States},
year = {2019},
month = {3}
}