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

Title: ECsim-CYL: Energy Conserving Semi-Implicit 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 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 plasma 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.

Authors:
 [1];  [1]; ORCiD logo [1];  [2]; ORCiD logo [3]
  1. Univ. of Leuven (Belgium). Dept. of Mathematics
  2. Energy Matter Conversion Corporation (EMC2), San Diego, CA (United States)
  3. 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 Laboratory-National Energy Research Scientific Computing Center (NERSC)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1530608
DOE Contract Number:  
AC02-05CH11231
Resource Type:
Journal Article
Journal Name:
Computer Physics Communications
Additional Journal Information:
Journal Volume: 236; Journal Issue: C; Journal ID: ISSN 0010-4655
Country of Publication:
United States
Language:
English

Citation Formats

Gonzalez-Herrero, Diego, Micera, Alfredo, Boella, Elisabetta, Park, Jaeyoung, and Lapenta, Giovanni. ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates. United States: N. p., 2019. Web. doi:10.1016/j.cpc.2018.10.026.
Gonzalez-Herrero, Diego, Micera, Alfredo, Boella, Elisabetta, Park, Jaeyoung, & Lapenta, Giovanni. ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates. United States. doi:10.1016/j.cpc.2018.10.026.
Gonzalez-Herrero, Diego, Micera, Alfredo, Boella, Elisabetta, Park, Jaeyoung, and Lapenta, Giovanni. Fri . "ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates". United States. doi:10.1016/j.cpc.2018.10.026.
@article{osti_1530608,
title = {ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates},
author = {Gonzalez-Herrero, 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 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 plasma 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.},
doi = {10.1016/j.cpc.2018.10.026},
journal = {Computer Physics Communications},
issn = {0010-4655},
number = C,
volume = 236,
place = {United States},
year = {2019},
month = {3}
}