Analysis of the performance of numerical integration methods for the tracking of ultra-high energy cosmic rays

被引:3
|
作者
Costa Jr, R. P. [1 ]
de Oliveira, M. A. Leigui [1 ]
机构
[1] Univ Fed ABC, Ctr Ciencias Nat & Humanas, Av Estados 5001, BR-09210580 Santo Andre, SP, Brazil
基金
巴西圣保罗研究基金会;
关键词
Particle tracking; Larmor radius; Magnetic fields; Ultra high energy cosmic rays;
D O I
10.1016/j.jcp.2019.04.058
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We analyzed the performance of several numerical integration methods in the ultra-relativistic regime. The integration methods include the fourth order Runge-Kutta, Boris, Vay, and a new method (called New Euler). For a proton traveling in a circular trajectory through a uniform perpendicular magnetic field, we have calculated the relative percentage error of Larmor radius of the particle's trajectory, relative efficiency between the methods, and relative percentage error of the particle's energy as a function of CPU time. In these calculations, we have investigated the influence of the following parameters: particle's energy, magnetic field strength, and the integration step. On calculating the Larmor radius, the four methods examined have all presented great precision. Among all studied cases, the most efficient method was the Boris method, which completes within a given CPU time more than twice cycles (complete circles, or periods of revolution) given by the Runge-Kutta one. The least efficient was the New Euler method (0.08 cycles), followed by Vay one (0.21 cycles). On the calculation of the particle's energy, New Euler method has exhibited the best properties, with the maximum relative percentage error of the particle's energy of roughly 10(-13)% remaining constant during the entire simulation. Except in the cases with the biggest simulation step, when for the Runge-Kutta method the particle's energy diverges asymptotically from its initial value, and with the smallest simulation step and lowest energy, when the Boris method has presented less precision than the Runge-Kutta one. The Runge-Kutta and Boris methods have presented approximately the same precision, with the maximum relative percentage error of similar to 10(-11)%. The Vay method has presented the worst precision, with the particle's energy diverging asymptotically for all studied cases, with the relative percentage error of the particle's energy reaching approximately 10(-8)%, in ten seconds of CPU time. (C) 2019 Elsevier Inc. All rights reserved.
引用
收藏
页码:432 / 443
页数:12
相关论文
共 50 条
  • [1] ULTRA-HIGH ENERGY COSMIC RAYS
    Watson, A. A.
    ACTA PHYSICA POLONICA B, 2019, 50 (12): : 2035 - 2056
  • [2] Ultra-high energy cosmic rays
    Department of Physics and Astronomy, University of Leeds Leeds, LS2 9JT, United Kingdom
    Acta Phys Pol B, 1600, 12 (2035-2056):
  • [3] Ultra-high energy cosmic rays
    Harari, Diego
    PHYSICS OF THE DARK UNIVERSE, 2014, 4 : 23 - 30
  • [4] On ultra-high energy cosmic rays
    Lemoine, Martin
    NEUTRINO OSCILLATION WORKSHOP, NOW2018, 2019,
  • [5] Ultra-high energy cosmic rays
    Kim, Hang Bae
    JOURNAL OF THE KOREAN PHYSICAL SOCIETY, 2021, 78 (10) : 912 - 917
  • [6] Ultra-high energy cosmic rays
    Sokolsky, P
    MODERN PHYSICS LETTERS A, 2004, 19 (13-16) : 959 - 966
  • [7] Ultra-high energy cosmic rays
    Hang Bae Kim
    Journal of the Korean Physical Society, 2021, 78 : 912 - 917
  • [8] Ultra-high energy cosmic rays
    Werner, K
    Chernatkin, V
    Engel, R
    Kalmykov, NN
    Ostapchenko, S
    Pierog, T
    IX HADRON PHYSICS AND VII RELATIVISTIC ASPECTS OF NUCLEAR PHYSICS, 2004, 739 : 385 - 396
  • [9] Ultra-High Energy Cosmic Rays and Neutrinos
    Nagataki, Shigehiro
    XIV MEXICAN SCHOOL ON PARTICLES AND FIELDS, 2011, 287
  • [10] The anisotropy of the ultra-high energy cosmic rays
    Evans, NW
    Ferrer, F
    Sarkar, S
    ASTROPARTICLE PHYSICS, 2002, 17 (03) : 319 - 340