A fast numerical integrator for relativistic charged particle tracking

doi:10.1016/j.nima.2017.04.015

Title: A fast numerical integrator for relativistic charged particle tracking

Full Record
Other Related Research

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 Laboratory (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center

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: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: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: 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

Similar Records in DOE PAGES and OSTI.GOV collections:

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: 10.1016/j.jcp.2015.09.007

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 6
DOI: 10.3847/1538-4365/aab114

Full Text Available
A high-order relativistic two-fluid electrodynamic scheme with consistent reconstruction of electromagnetic fields and a multidimensional Riemann solver for electromagnetism

Journal Article Balsara, Dinshaw S., E-mail: dbalsara@nd.edu ; Amano, Takanobu ; Garain, Sudip ; ... - Journal of Computational Physics

In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equationsmore »« less
DOI: 10.1016/J.JCP.2016.05.006
Energy Loss by Radiation in Many-Particle Numerical Simulation With Lorentz-Dirac Equation

Journal Article Zacek, Martin - AIP Conference Proceedings

We studied the possibilities for numerical integration of Lorentz-Dirac equation that is the equation describing the motion of a charged point particle when radiation reaction is taken into account. In numerical modelling based on particle models usually the equations of motion without radiation force are used and the corrections for radiation are used consequently, expressed by laws given by averaged particle parameters as the temperature or particle density. If the complete equation of motion concluding the radiation would be used, the corrections for radiation reaction force could be used for every charged particle individually from more fundamental laws. Thus themore »« less
DOI: 10.1063/1.2168812
Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

Journal Article Godfrey, Brendan B., E-mail: brendan.godfrey@ieee.org ; Vay, Jean-Luc - Journal of Computational Physics

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the commonmore »« less
DOI: 10.1016/J.JCP.2013.04.006

Similar Records