An Efficient Numerical Scheme for Simulating Particle Acceleration in Evolving Cosmic-Ray Modified Shocks
arXiv:astro-ph/0506212 · doi:10.1016/j.astropartphys.2005.05.006
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 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.
Accepted for publication in Astroparticle Physics; 19 pages and 5 figures
References in corpus (12)
- A semi-analytical approach to non-linear shock acceleration
- Analytic solution for nonlinear shock acceleration in the Bohm limit
- Numerical Studies of Cosmic Ray Injection and Acceleration
- Simulating Electron Transport and Synchrotron Emission in Radio Galaxies: Shock Acceleration and Synchrotron Aging in Axis-Symmetric Flows
- COSMOCR: A Numerical Code for Cosmic Ray Studies in Computational Cosmology
- Radio Emission from a Young Supernova Remnant Interacting with an Interstellar Cloud: MHD Simulation with Relativistic Electrons
- Spectral universality of strong shocks accelerating charged particles
- Diffusive Shock Acceleration Simulations: Comparison with Particle Methods and Bow Shock Measurements
- Diffusive Shock Acceleration in Oblique MHD Shocks: Comparison with Monte Carlo Methods and Observations
- Time Evolution of Cosmic-Ray Modified Plane Shocks
- Time dependent cosmic-ray shock acceleration with self-consistent injection
- Computation of diffusive shock acceleration using stochastic differential equations
Cited by in corpus (24)
- Detecting shock waves in cosmological smoothed particle hydrodynamics simulations
- Spectral evolution of superluminal components in parsec-scale jets
- Cosmological MHD Simulations of Galaxy Cluster Radio Relics: Insights and Warnings for Observations
- Numerical Studies of Diffusive Shock Acceleration at Spherical Shocks
- MHD Simulations of AGN Jets in a Dynamic Galaxy Cluster Medium
- First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations
- Comparison of Different Methods for Nonlinear Diffusive Shock Acceleration
- Spectrally resolved cosmic ray hydrodynamics -- I. Spectral scheme
- Interactions Between Radio Galaxies and Cluster Shocks - 1: Jet Axes Aligned with Shock Normals
- Simulations of cosmic ray propagation
- Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. I. Spectrally resolved propagation of CR electrons on Eulerian grids
- CRESCENDO: An on-the-fly Fokker-Planck Solver for Spectral Cosmic Rays in Cosmological Simulations
- A Fresh Look at Narrow-Angle Tail Radio Galaxy Dynamics, Evolution and Emissions
- Glimm-Godunov's Method for Cosmic-ray-hydrodynamics
- Simulations of Precessing Jets and the Formation of X-shaped Radio Galaxies
- Plasma physics of the intracluster medium
- Simulated Interactions Between Radio Galaxies and Cluster Shocks -- 2: Jet Axes Orthogonal to Shock Normals
- Shocked Narrow-Angle Tail Radio Galaxies: Simulations and Emissions
- MARCOS, a numerical tool for the simulation of multiple time-dependent non-linear diffusive shock acceleration
- An accurate treatment of scattering and diffusion in piecewise power-law models for cosmic ray and radiation/neutrino transport
- Non-thermal emissions from a head-tail radio galaxy in 3D magnetohydrodynamic simulations
- Observations and Simulations of Radio Emission and Magnetic Fields in Minkowski's Object
- Time Evolution of Cosmic Ray MHD Shocks and Their Emissions
- Simulating cosmic ray electron spectra and radio emission from an AGN jet outburst in a cool-core cluster