A fast numerical integrator for relativistic charged particle tracking

Qiang, Ji

doi:10.1016/j.nima.2017.04.015

Title: A fast numerical integrator for relativistic charged particle tracking

Abstract

In this paper, we report on a fast second-order numerical integrator to solve the Lorentz force equations of a relativistic charged particle in electromagnetic fields. This numerical integrator shows less numerical error than the popular Boris algorithm in tracking the relativistic particle subject to electric and magnetic space-charge fields and requires less number of operations than another recently proposed relativistic integrator.

Authors:

Qiang, Ji ^[1]

Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)

Publication Date:: 2017-09-01

Research Org.:: Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 1526495

Alternate Identifier(s):: OSTI ID: 1550495

Grant/Contract Number:: AC02-05CH11231

Resource Type:: Accepted Manuscript

Journal Name:: Nuclear Instruments and Methods in Physics Research. Section A, Accelerators, Spectrometers, Detectors and Associated Equipment

Additional Journal Information:: Journal Volume: 867; Journal Issue: C; Journal ID: ISSN 0168-9002

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY

Citation Formats


                    Qiang, Ji. A fast numerical integrator for relativistic charged particle tracking.  United States: N. p., 2017. 
        Web.  doi:10.1016/j.nima.2017.04.015.

Copy to clipboard


                    Qiang, Ji. A fast numerical integrator for relativistic charged particle tracking.  United States.  doi:https://doi.org/10.1016/j.nima.2017.04.015

Copy to clipboard


                    Qiang, Ji. Fri .  
        "A fast numerical integrator for relativistic charged particle tracking".  United States.  doi:https://doi.org/10.1016/j.nima.2017.04.015.  https://www.osti.gov/servlets/purl/1526495.

Copy to clipboard


                    
@article{osti_1526495,

  title        = {A fast numerical integrator for relativistic charged particle tracking},

  author       = {Qiang, Ji},

  abstractNote = {In this paper, we report on a fast second-order numerical integrator to solve the Lorentz force equations of a relativistic charged particle in electromagnetic fields. This numerical integrator shows less numerical error than the popular Boris algorithm in tracking the relativistic particle subject to electric and magnetic space-charge fields and requires less number of operations than another recently proposed relativistic integrator.},

  doi          = {10.1016/j.nima.2017.04.015},

  journal      = {Nuclear Instruments and Methods in Physics Research. Section A, Accelerators, Spectrometers, Detectors and Associated Equipment},

  number       = C,

  volume       = 867,

  place        = {United States},

  year         = {2017},

  month        = {9}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

DOI: https://doi.org/10.1016/j.nima.2017.04.015

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 2 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Works referencing / citing this record:

On the Boris solver in particle-in-cell simulation
journal, November 2018

Zenitani, Seiji; Umeda, Takayuki
Physics of Plasmas, Vol. 25, Issue 11
DOI: 10.1063/1.5051077

A Comprehensive Comparison of Relativistic Particle Integrators
journal, March 2018

Ripperda, B.; Bacchini, F.; Teunissen, J.
The Astrophysical Journal Supplement Series, Vol. 235, Issue 1
DOI: 10.3847/1538-4365/aab114

Similar Records in DOE PAGES and OSTI.GOV collections:

Generalized, Energy-conserving Numerical Simulations of Particles in General Relativity. II. Test Particles in Electromagnetic Fields and GRMHD

Journal Article Bacchini, F. ; Ripperda, B. ; Porth, O. ; ... - The Astrophysical Journal. Supplement Series (Online)

Observations of compact objects, in the form of radiation spectra, gravitational waves by LIGO/Virgo, and direct imaging with the Event Horizon Telescope, are currently the main information sources on plasma physics in extreme gravity. Modeling such physical phenomena Requires numerical methods that allow for the simulation of microscopic plasma dynamics in the presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018), we presented a detailed study of numerical techniques for the integration of free geodesic motion. In this paper, we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes.more »« less
Cited by 6
DOI: https://doi.org/10.3847/1538-4365/aafcb3

Full Text Available
On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

Journal Article Xu, Xinlu ; Li, Fei ; Tsung, Frank S. ; ... - Journal of Computational Physics

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwellmore »« less
Cited by 1
DOI: https://doi.org/10.1016/j.jcp.2020.109451

Full Text Available
Comment on “Symplectic integration of magnetic systems”: A proof that the Boris algorithm is not variational

Journal Article Ellison, C. L. ; Burby, J. W. ; Qin, H. - Journal of Computational Physics

One popular technique for the numerical time advance of charged particles interacting with electric and magnetic fields according to the Lorentz force law [1], [2], [3] and [4] is the Boris algorithm. Its popularity stems from simple implementation, rapid iteration, and excellent long-term numerical fidelity [1] and [5]. Excellent long-term behavior strongly suggests the numerical dynamics exhibit conservation laws analogous to those governing the continuous Lorentz force system [6]. Moreover, without conserved quantities to constrain the numerical dynamics, algorithms typically dissipate or accumulate important observables such as energy and momentum over long periods of simulated time [6]. Identification of themore »« less
Cited by 3
DOI: https://doi.org/10.1016/j.jcp.2015.09.007

Full Text Available
A semi-implicit, energy- and charge-conserving particle-in-cell algorithm for the relativistic Vlasov-Maxwell equations

Journal Article Chen, Guangye ; Chacon, Luis ; Yin, Lin ; ... - Journal of Computational Physics

Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in the light-wave modes. This can lead to intolerable simulation errors where accurate light propagation is needed (e.g. in laser-plasma interactions). In this study, we selectively combine the leap-frog and Crank-Nicolson methods to produce an exactly energy- and charge-conserving relativistic electromagnetic PIC algorithm. Specifically, we employ the leap-frog method for Maxwell's equations, and the Crank-Nicolson method for the particle equations. The semi-implicit formulation still features a timestep CFL, but facilitates exact globalmore »« less
DOI: https://doi.org/10.1016/j.jcp.2020.109228

Full Text Available
A Comprehensive Comparison of Relativistic Particle Integrators

Journal Article Ripperda, B. ; Bacchini, F. ; Teunissen, J. ; ... - The Astrophysical Journal. Supplement Series (Online)

We compare relativistic particle integrators commonly used in plasma physics, showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely, the Boris, Vay, and Higuera–Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Here, our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback onto the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields,more »« less
Cited by 11
DOI: https://doi.org/10.3847/1538-4365/aab114

Full Text Available

Similar Records

Title: A fast numerical integrator for relativistic charged particle tracking

Abstract

Citation Formats

On the Boris solver in particle-in-cell simulation journal, November 2018

A Comprehensive Comparison of Relativistic Particle Integrators journal, March 2018

On the Boris solver in particle-in-cell simulation
journal, November 2018

A Comprehensive Comparison of Relativistic Particle Integrators
journal, March 2018