Pranab J. Deka,Ralf Kissmann,Lukas Einkemmer
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.
High Energy Astrophysical Phenomena,Instrumentation and Methods for Astrophysics,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the efficient numerical simulation method for the anisotropic diffusion of galactic cosmic rays (Cosmic Rays, CRs). Specifically, the author proposes a new numerical method - the exponential method - for dealing with linear anisotropic diffusion equations, especially in the presence of convection and time - dependent or non - time - dependent sources. These methods allow for larger time steps in simulations while maintaining high precision, thereby significantly increasing the speed of the simulation.
### Background and Problem Description
The propagation of cosmic rays in the Milky Way is mainly achieved through the diffusion mechanism, and anisotropic diffusion plays a crucial role in understanding the interaction of cosmic rays in the Milky Way, the heliosphere, and with the galactic magnetic field. This diffusion term makes the cosmic - ray transport equation highly stiff. Traditional numerical methods, such as implicit integrators (especially the Crank - Nicolson method), although they can use larger time steps, are still limited by the CFL condition, resulting in low computational efficiency.
### Proposed Method
The author proposes to use the exponential method to deal with linear anisotropic diffusion equations. These methods combine the technique of polynomial interpolation on Leja points and can use larger time steps while ensuring high precision, thus significantly accelerating the simulation process.
### Specific Contributions
1. **Proposal of the Exponential Method**: The author proposes the Exponential Quadrature Rules and details how to use Leja points for polynomial interpolation to calculate the matrix exponential function and the related \(\varphi_k\) functions.
2. **Performance Testing**: Through several different magnetic field configurations (such as concentric circles and spiral magnetic fields), the author shows the superiority of the proposed exponential method in dealing with anisotropic diffusion problems. In particular, in the presence of convection and time - dependent sources, the exponential Gauss quadrature method shows higher precision and computational efficiency.
3. **Comparison with Other Methods**: The author compares the performance of the exponential method with the traditional Crank - Nicolson method through the Work - Precision Diagram, and the results show that the exponential method has significant advantages at large time steps.
### Conclusions and Prospects
The author points out that the proposed exponential method is superior to the traditional Crank - Nicolson method in terms of both computational time and precision. Although this paper only considers two - dimensional problems, the author expects that these methods will also show similar performance improvements when dealing with the complete four - dimensional cosmic - ray transport equation. Future work will include integrating the LeXInt library into the P/i.pc/c.pc/a.pc/r.pc/d.pc code to study the cosmic - ray transport in the Milky Way and consider continuous energy loss and re - acceleration effects.
### Mathematical Formulas
The main equations involved in the paper are as follows:
1. **Cosmic - Ray Transport Equation**:
\[
\frac{\partial f}{\partial t}-\nabla\cdot(D\nabla f)+\nabla\cdot(\vec{v}f)-\frac{\partial}{\partial p}\left(\frac{p^2}{3}(\nabla\cdot\vec{v})f\right)+\left(\frac{1}{\tau_p}+\frac{1}{\tau_r}\right)f = Q(\vec{x},p,t)
\]
2. **Two - Dimensional Anisotropic Diffusion Equation**:
\[
\frac{\partial f}{\partial t}=\nabla\cdot(D\nabla f)+\nabla\cdot(\vec{v}f)+Q(x,y,t)
\]
3. **Diffusion Tensor**:
\[
D=\begin{bmatrix}
D_{xx}(x,y)&D_{xy}(x,y)\\
D_{yx}(x,y)&D_{yy}(x,y)
\end{bmatrix}
\]