MARCOS, a numerical tool for the simulation of multiple time-dependent non-linear diffusive shock acceleration

Gilles Ferrand (1) , Turlough Downes (2) , Alexandre Marcowith (3)

1 CESR - Centre d'étude spatiale des rayonnements
2 School of Mathematical Science
3 LPTA - Laboratoire de Physique Théorique et Astroparticules

We present a new code aimed at the simulation of diffusive shock acceleration (DSA), and discuss various test cases which demonstrate its ability to study DSA in its full time-dependent and non-linear developments. We present the numerical methods implemented, coupling the hydrodynamical evolution of a parallel shock (in one space dimension) and the kinetic transport of the cosmic-rays (CR) distribution function (in one momentum dimension), as first done by Falle. Following Kang and Jones and collaborators, we show how the adaptive mesh refinement technique (AMR) greatly helps accommodating the extremely demanding numerical resolution requirements of realistic (Bohm-like) CR diffusion coefficients. We also present the paral lelization of the code, which allows us to run many successive shocks at the cost of a single shock, and thus to present the first direct numerical simulations of linear and non-linear multiple DSA, a mechanism of interest in various astrophysical environments such as superbubbles, galaxy clusters and early cosmological flows.

Mots clés: cosmic-rays; acceleration of particles; numerical simulations

Domaines: Cosmologie et astrophysique extra-galactique [astro-ph.CO]

Gilles Ferrand, Turlough Downes, Alexandre Marcowith. MARCOS, a numerical tool for the simulation of multiple time-dependent non-linear diffusive shock acceleration. Monthly Notices of the Royal Astronomical Society, 2008, 383, pp.41. ⟨10.1111/j.1365-2966.2007.12511.x⟩. ⟨hal-00177684⟩