Abstract
We present a numerical algorithm for the incorporation of the active cosmic ray transport into the ZEUS-3D magnetohydrodynamical code. The cosmic ray transport is described by the diffusion-advection equation. The applied form of the diffusion tensor allows for anisotropic diffusion of cosmic rays along and across the magnetic field direction, which is controlled by two parameters: the parallel and perpendicular diffusion coefficients. The implemented numerical algorithm is tested by comparison of the diffusive transport of cosmic rays to analytical solutions of the diffusion equation. Our method is numerically stable for a wide range of diffusion coefficients, including the realistic values inferred from the observational data for the Milky Way of about 6x10^28 cm^2 s^(-1) . The presented algorithm is applied for exemplary simulations of the Parker instability triggered by cosmic rays injected by a single SN remnant.
Dokumententyp: | Zeitschriftenartikel |
---|---|
Keywords: | ISM: magnetic fields / ISM: cosmic rays |
Fakultät: | Physik |
Themengebiete: | 500 Naturwissenschaften und Mathematik > 520 Astronomie
500 Naturwissenschaften und Mathematik > 530 Physik |
ISSN: | 0004-6361 |
Sprache: | Englisch |
Dokumenten ID: | 71525 |
Datum der Veröffentlichung auf Open Access LMU: | 02. Apr. 2020, 10:14 |
Letzte Änderungen: | 04. Nov. 2020, 13:52 |