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 language||English (US)|
|Number of pages||17|
|State||Published - Sep 2005|
Bibliographical noteFunding 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).
- Particle acceleration