Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
157 tokens/sec
GPT-4o
8 tokens/sec
Gemini 2.5 Pro Pro
46 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

Efficient numerical methods for Anisotropic Diffusion of Galactic Cosmic Rays (2307.12276v2)

Published 23 Jul 2023 in astro-ph.HE, astro-ph.IM, cs.NA, and math.NA

Abstract: Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of the cosmic ray transport equation. To conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators (namely, Crank-Nicolson (CN)) have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods to treat the linear anisotropc diffusion equation in the presence of advection and time-independent and time-dependent sources. These methods allow us to take even larger step sizes that can substantially speed-up the simulations whilst generating highly accurate solutions. In or subsequent work, we will use these exponential solvers in the Picard code to study anisotropic cosmic ray diffusion and we will consider additional physical processes such as continuous momentum losses and reacceleration.

Citations (1)

Summary

We haven't generated a summary for this paper yet.