EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Klaus Steiniger; Rene Widera; Sergei Bastrakov; Michael Bussmann; Sunita Chandrasekaran; Benjamin Hernandez; Kristina Holsapple; Axel Huebl; Guido Juckeland; Jeffrey Kelling; Matt Leinhauser; Richard Pausch; David Rogers; Ulrich Schramm; Jeff Young; Alexander Debus

doi:10.1016/j.cpc.2023.108849

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Klaus Steiniger, Rene Widera, Sergei Bastrakov, Michael Bussmann, Sunita Chandrasekaran, Benjamin Hernandez, Kristina Holsapple, Axel Huebl, Guido Juckeland, Jeffrey Kelling, Matt Leinhauser, Richard Pausch, David Rogers, Ulrich Schramm, Jeff Young, Alexander Debus

Advanced Computing for Chemistry and Mat

Research output: Contribution to journal › Article › peer-review

Abstract

We present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

Original language	English
Article number	108849
Journal	Computer Physics Communications
Volume	291
DOIs	https://doi.org/10.1016/j.cpc.2023.108849
State	Published - Oct 2023

Funding

This research used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-00OR22725 . This work was partly funded by the Center for Advanced Systems Understanding (CASUS) which is financed by Germany's Federal Ministry of Education and Research (BMBF) and by the Saxon Ministry for Science, Culture and Tourism (SMWK) with tax funds on the basis of the budget approved by the Saxon State Parliament.

Funders	Funder number
Center for Advanced Systems Understanding
Saxon Ministry for Science, Culture and Tourism
Office of Science	DE-AC05-00OR22725
Bundesministerium für Bildung und Forschung
Sächsisches Staatsministerium für Wissenschaft und Kunst

Keywords

Charge conservation
Current deposition
EZ
PIConGPU
Particle-in-cell

Access to Document

10.1016/j.cpc.2023.108849

Cite this

Steiniger, K., Widera, R., Bastrakov, S., Bussmann, M., Chandrasekaran, S., Hernandez, B., Holsapple, K., Huebl, A., Juckeland, G., Kelling, J., Leinhauser, M., Pausch, R., Rogers, D., Schramm, U., Young, J., & Debus, A. (2023). EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations. Computer Physics Communications, 291, Article 108849. https://doi.org/10.1016/j.cpc.2023.108849

@article{ae1d4a1c618c4c189c533030fee2b896,

title = "EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations",

abstract = "We present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.",

keywords = "Charge conservation, Current deposition, EZ, PIConGPU, Particle-in-cell",

author = "Klaus Steiniger and Rene Widera and Sergei Bastrakov and Michael Bussmann and Sunita Chandrasekaran and Benjamin Hernandez and Kristina Holsapple and Axel Huebl and Guido Juckeland and Jeffrey Kelling and Matt Leinhauser and Richard Pausch and David Rogers and Ulrich Schramm and Jeff Young and Alexander Debus",

note = "Publisher Copyright: {\textcopyright} 2023 Elsevier B.V.",

year = "2023",

month = oct,

doi = "10.1016/j.cpc.2023.108849",

language = "English",

volume = "291",

journal = "Computer Physics Communications",

issn = "0010-4655",

publisher = "Elsevier",

}

Steiniger, K, Widera, R, Bastrakov, S, Bussmann, M, Chandrasekaran, S, Hernandez, B, Holsapple, K, Huebl, A, Juckeland, G, Kelling, J, Leinhauser, M, Pausch, R, Rogers, D, Schramm, U, Young, J & Debus, A 2023, 'EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations', Computer Physics Communications, vol. 291, 108849. https://doi.org/10.1016/j.cpc.2023.108849

TY - JOUR

T1 - EZ

T2 - An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

AU - Steiniger, Klaus

AU - Widera, Rene

AU - Bastrakov, Sergei

AU - Bussmann, Michael

AU - Chandrasekaran, Sunita

AU - Hernandez, Benjamin

AU - Holsapple, Kristina

AU - Huebl, Axel

AU - Juckeland, Guido

AU - Kelling, Jeffrey

AU - Leinhauser, Matt

AU - Pausch, Richard

AU - Rogers, David

AU - Schramm, Ulrich

AU - Young, Jeff

AU - Debus, Alexander

PY - 2023/10

Y1 - 2023/10

N2 - We present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

AB - We present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

KW - Charge conservation

KW - Current deposition

KW - EZ

KW - PIConGPU

KW - Particle-in-cell

UR - http://www.scopus.com/inward/record.url?scp=85165527540&partnerID=8YFLogxK

U2 - 10.1016/j.cpc.2023.108849

DO - 10.1016/j.cpc.2023.108849

M3 - Article

AN - SCOPUS:85165527540

SN - 0010-4655

VL - 291

JO - Computer Physics Communications

JF - Computer Physics Communications

M1 - 108849

ER -

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Abstract

Funding

Keywords

Access to Document

Other files and links

Fingerprint

Cite this