English
 
Help Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh

MPS-Authors
/persons/resource/persons4732

Pakmor,  R.
Stellar Astrophysics, MPI for Astrophysics, Max Planck Society;

External Resource
No external resources are shared
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available
Citation

Thomas, T., Pfrommer, C., & Pakmor, R. (2021). A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh. Monthly Notices of the Royal Astronomical Society, 503(2), 2242-2264. doi:10.1093/mnras/stab397.


Cite as: http://hdl.handle.net/21.11116/0000-0008-F371-6
Abstract
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 Alfvén 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.