An efficient numerical scheme for simulating particle acceleration in evolving cosmic-ray modified shocks

T. W. Jones, Hyesung Kang

Research output: Contribution to journalArticlepeer-review

17 Scopus citations

Abstract

We have developed a new, very efficient numerical scheme to solve the CR diffusion-convection equation that can be applied to the study of the nonlinear time evolution of CR modified shocks for arbitrary spatial diffusion properties. The efficiency of the scheme derives from its use of coarse-grained finite momentum volumes. This approach has enabled us, using ∼10-20 momentum bins spanning nine orders of magnitude in momentum, to carry out simulations that agree well with results from simulations of modified shocks carried out with our conventional finite difference scheme requiring more than an order of magnitude more momentum points. The coarse-grained, CGMV scheme reduces execution times by a factor approximately half the ratio of momentum bins used in the two methods. Depending on the momentum dependence of the diffusion, additional economies in required spatial and time resolution can be utilized in the CGMV scheme, as well. These allow a computational speed-up of at least an order of magnitude in some cases.

Original languageEnglish (US)
Pages (from-to)75-91
Number of pages17
JournalAstroparticle Physics
Volume24
Issue number1-2
DOIs
StatePublished - Sep 2005

Bibliographical note

Funding Information:
TWJ is supported by NSF grant AST03-07600, by NASA grants NAG5-10774, NNG05GF57G and by the University of Minnesota Supercomputing Institute. HK was supported by KOSEF through the Astrophysical Research Center for the Structure and Evolution of Cosmos (ARCSEC).

Keywords

  • Cosmic-rays
  • Particle acceleration

Fingerprint

Dive into the research topics of 'An efficient numerical scheme for simulating particle acceleration in evolving cosmic-ray modified shocks'. Together they form a unique fingerprint.

Cite this