A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh
- We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh AREPO code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretized using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of AREPO. The interaction of CRs and gyroresonant Alfven waves is described by short time-scale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magnetohydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetized discontinuity, the hydrodynamic response to a CR overpressure,We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh AREPO code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretized using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of AREPO. The interaction of CRs and gyroresonant Alfven waves is described by short time-scale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magnetohydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetized discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magnetohydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions.…
Verfasserangaben: | Timon ThomasORCiDGND, Christoph PfrommerORCiDGND, Rüdiger PakmorORCiDGND |
---|---|
DOI: | https://doi.org/10.1093/mnras/stab397 |
ISSN: | 0035-8711 |
ISSN: | 1365-2966 |
Titel des übergeordneten Werks (Englisch): | Monthly notices of the Royal Astronomical Society |
Verlag: | Oxford University Press |
Verlagsort: | Oxford |
Publikationstyp: | Wissenschaftlicher Artikel |
Sprache: | Englisch |
Datum der Erstveröffentlichung: | 12.02.2021 |
Erscheinungsjahr: | 2021 |
Datum der Freischaltung: | 01.03.2024 |
Freies Schlagwort / Tag: | MHD; cosmic rays; hydrodynamics; methods: numerical |
Band: | 503 |
Ausgabe: | 2 |
Seitenanzahl: | 23 |
Erste Seite: | 2242 |
Letzte Seite: | 2264 |
Fördernde Institution: | European Research Council under ERCCoG grant CRAGSMAN [646955]; National Science FoundationNational Science Foundation (NSF) [NSF PHY1748958] |
Organisationseinheiten: | Mathematisch-Naturwissenschaftliche Fakultät / Institut für Physik und Astronomie |
DDC-Klassifikation: | 5 Naturwissenschaften und Mathematik / 52 Astronomie / 520 Astronomie und zugeordnete Wissenschaften |
Peer Review: | Referiert |