An accurate treatment of scattering and diffusion in piecewise power-law models for cosmic ray and radiation/neutrino transport
arXiv:2202.05283 · doi:10.1093/mnras/stac3283
Abstract
A popular numerical method to model the dynamics of a 'full spectrum' of cosmic rays (CRs), also applicable to radiation/neutrino hydrodynamics (RHD), is to discretize the spectrum at each location/cell as a piecewise power law in 'bins' of momentum (or frequency) space. This gives rise to a pair of conserved quantities (e.g. CR number and energy) which are exchanged between cells or bins, that in turn give the update to the normalization and slope of the spectrum in each bin. While these methods can be evolved exactly in momentum-space (considering injection, absorption, continuous losses/gains), numerical challenges arise dealing with spatial fluxes, if the scattering rates depend on momentum. This has often been treated by either by neglecting variation of those rates 'within the bin,' or sacrificing conservation -- introducing significant errors. Here, we derive a rigorous treatment of these terms, and show that the variation within the bin can be accounted for accurately with a simple set of scalar correction coefficients that can be written entirely in terms of other, explicitly-evolved 'bin-integrated' quantities. This eliminates the relevant errors without added computational cost, has no effect on the numerical stability of the method, and retains manifest conservation. We derive correction terms both for methods which explicitly integrate flux variables (two-moment or M1-like) methods, as well as single-moment (advection-diffusion, FLD-like) methods, and approximate corrections valid in various limits.
10 pages, 1 figure. Updated to match published version
References in corpus (17)
- Cosmic Ray Propagation: Nonlinear Diffusion Parallel and Perpendicular to Mean Magnetic Field
- Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients
- A scheme for radiation pressure and photon diffusion with the M1 closure in RAMSES-RT
- Gamma-ray sky points to radial gradients in cosmic-ray transport
- Cosmic ray-driven galactic winds: streaming or diffusion?
- The dependence of cosmic ray driven galactic winds on halo mass
- Galactic cosmic-ray propagation in the light of AMS-02: I. Analysis of protons, helium, and antiprotons
- Monte-Carlo closure for moment-based transport schemes in general relativistic radiation hydrodynamics simulations
- First Predicted Cosmic Ray Spectra, Primary-to-Secondary Ratios, and Ionization Rates from MHD Galaxy Formation Simulations
- Standard Self-Confinement and Extrinsic Turbulence Models for Cosmic Ray Transport are Fundamentally Incompatible with Observations
- Implications of current nuclear cross sections on secondary cosmic rays with the upcoming DRAGON2 code
- Implications of Lithium to Oxygen AMS-02 spectra on our understanding of cosmic-ray diffusion
- Fermi-LAT Observations of High- and Intermediate-Velocity Clouds: Tracing Cosmic Rays in the Halo of the Milky Way
- The Spatially Uniform Spectrum of the Fermi Bubbles: the Leptonic AGN Jet Scenario
- Simulations of cosmic ray propagation
- A Consistent Reduced-Speed-of-Light Formulation of Cosmic Ray Transport Valid in Weak and Strong-Scattering Regimes
- Glimm-Godunov's Method for Cosmic-ray-hydrodynamics