A New Time-Dependent Finite Difference Method for Relativistic Shock Acceleration
arXiv:1111.5795 · doi:10.1111/j.1365-2966.2011.20257.x
Abstract
We present a new approach to calculate the particle distribution function about relativistic shocks including synchrotron losses using the method of lines with an explicit finite difference scheme. A steady, continuous, one dimensional plasma flow is considered to model thick (modified) shocks, leading to a calculation in three dimensions plus time, the former three being momentum, pitch angle and position. The method accurately reproduces the expected power law behaviour in momentum at the shock for upstream flow speeds ranging from 0.1c to 0.995c (1 < Γ< 10). It also reproduces approximate analytical results for the synchrotron cutoff shape for a non-relativistic shock, demonstrating that the loss process is accurately represented. The algorithm has been implemented as a hybrid OpenMP--MPI parallel algorithm to make efficient use of SMP cluster architectures and scales well up to many hundreds of CPUs.
Accepted for publication in MNRAS
References in corpus (6)
- Particle acceleration by ultra-relativistic shocks: theory and simulations
- Particle acceleration in relativistic collisionless shocks: Fermi process at last?
- Magnetic field evolution in relativistic unmagnetized collisionless shocks
- A relativistic partially electromagnetic planar plasma shock
- Diffusive shock acceleration in radiation dominated environments
- Radiative losses and cut-offs of energetic particles at relativistic shocks