Algorithms for Forward and Backward Solution of the Fokker-Planck Equation in the Heliospheric Transport of Cosmic Rays

摘要

Motion of charged particles in an inhomogeneous turbulent medium as magnetic field is described by partial differential equations of the Fokker-Planck-Kolmogorov type. We present an algorithm of numerical solution of the four-dimensional Fokker-Planck equation in three-dimensional spherical coordinates system. The algorithm is based on Monte Carlo simulations of the stochastic motion of quasi-particles guided by the set of stochastic differential equations corresponding to the Fokker-Planck equation by the Ito formalism. We present the parallel algorithm in Julia programming language. We simulate the transport of cosmic rays in the heliosphere considering the full three-dimensional diffusion tensor. We compare forward- and backward-in-time solutions of the transport equation and discuss its computational advantages and disadvantages.