A guiding center implementation for relativistic particle dynamics in the PLUTO code
Abstract
We present a numerical implementation of the guiding center approximation to describe the relativistic motion of charged test particles in the PLUTO code for astrophysical plasma dynamics. The guiding center approximation (GCA) removes the time step constraint due to particle gyration around magnetic field lines by following the particle center of motion rather than its full trajectory. The gyration can be detached from the guiding center motion if electromagnetic fields vary sufficiently slow compared to the particle gyration radius and period. Our implementation employs a variable step-size linear multistep method, more efficient when compared to traditional one-step Runge Kutta schemes. A number of numerical benchmarks is presented in order to assess the validity of our implementation.
- Publication:
-
Computer Physics Communications
- Pub Date:
- April 2023
- DOI:
- arXiv:
- arXiv:2212.08064
- Bibcode:
- 2023CoPhC.28508625M
- Keywords:
-
- Numerical methods;
- Relativistic particles;
- Guiding center approximation;
- PLUTO code;
- Astrophysical plasma;
- Astrophysics - Instrumentation and Methods for Astrophysics
- E-Print:
- doi:10.1016/j.cpc.2022.108625