Numerical integration of stochastic differential equations: A parallel cosmic ray modulation implementation on Africa’s fastest computer

Three-dimensional studies of the transport and modulation of cosmic ray particles in turbulent astrospheres require large-scale simulations using specialized scientific codes. Essentially, a multi-dimensional Fokker-Planck type equation (a parabolic diffusion equation) must be integrated numerically. One such approach is to convert the relevant transport equation into a set of stochastic differential equations (SDEs), with the latter much easier to handle numerically. Due to the growing demand for high performance computing resources, research into the application of effective and suitable numerical algorithms to solve such equations is needed. We present a case study of the performance of a custom-written FORTRAN SDE numerical solver on the CHPC (Centre for High Performance Computing) Lengau cluster in South Africa for a realistic test problem with different set-ups. It is shown that SDE codes can scale very well on large parallel computing platforms. Finally, we consider an extremely computationally expensive application of the SDE approach to cosmic ray modulation, studying the behaviour of galactic cosmic ray proton latitude gradients and relative amplitudes in a physics-first manner. This is done using a modulation code that employs diffusion coefficients derived from first principles, which in turn are functions of turbulence quantities in reasonable agreement with spacecraft observations and modelled using a two-component turbulence transport model (TTM). We show that this approach leads to reduced latitude gradients qualitatively in line with spacecraft observations of the same, without making ad hoc assumptions as to anisotropic perpendicular diffusion coefficients as are often made in many cosmic ray modulation studies.