Numerical Model of Cosmic Ray Propagation in the Anisotropic Magnetic Field of the Galaxy

Abstract

The problem of diffusive propagation of cosmic rays (CRs) with a general tensor in a model homogeneous magnetic field of the Galaxy has been solved, and its stability has been investigated. An algorithm for generating stencils of second-order accuracy finite-difference schemes for a three-dimensional differential equation with tensor components is implemented. The components of the diffusion tensor are calculated with a specified spatial step and are determined based on the magnitudes of the large-scale magnetic field vector and the magnitude of the random field component, which is assumed to be 6 \(\mu\)G. The study demonstrates the fundamental possibility of determining the energy spectra of CRs in various regions of the Galaxy, where the free parameters are the distribution of CR sources and the structure of the magnetic field.