A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks

Abstract

Cosmic rays (CRs) are frequently modeled as an additional fluid in hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations of astrophysical flows. The standard CR two-fluid model is described in terms of three conservation laws (expressing conservation of mass, momentum and total energy) and one additional equation (for the CR pressure) that cannot be cast in a satisfactory conservative form. The presence of non-conservative terms with spatial derivatives in the model equations prevents a unique weak solution behind a shock. We investigate a number of methods for the numerical solution of the two-fluid equations and find that, in the presence of shock waves, the results generally depend on the numerical details (spatial reconstruction, time stepping, the CFL number, and the adopted discretization). All methods converge to a unique result if the energy partition between the thermal and non-thermal fluids at the shock is prescribed using a subgrid prescription. This highlights the non-uniqueness problem of the two-fluid equations at shocks. From our numerical investigations, we report a robust method for which the solutions are insensitive to the numerical details even in absence of a subgrid prescription, although we recommend a subgrid closure at shocks using results from kinetic theory. The subgrid closure is crucial for a reliable post-shock solution and also its impact on large scale flows because the shock microphysics that determines CR acceleration is not accurately captured in a fluid approximation. Critical test problems, limitations of fluid modeling, and future directions are discussed.

18 pages, 12 figures, 2 tables; MNRAS accepted version

References in corpus (11)

• PLUTO: a Numerical Code for Computational Astrophysics
• Feedback Heating by Cosmic Rays in Clusters of Galaxies
• Detecting shock waves in cosmological smoothed particle hydrodynamics simulations
• The Role of Cosmic Ray Transport in Shaping the Simulated Circumgalactic Medium
• A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations
• Interaction of Cosmic Rays with Cold Clouds in Galactic Halos
• Constraining cosmic ray acceleration in young star clusters using multi-wavelength observations
• Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping
• Lack of thermal energy in superbubbles: hint of cosmic rays?
• Approximate Riemann Solvers for the Cosmic Ray Magnetohydrodynamical Equations
• Fluid Simulations of Cosmic Ray Modified Shocks

Cited by in corpus (13)

• Simulations of cosmic ray propagation
• Dynamical Effects of Cosmic Rays on the Medium Surrounding Their Sources
• Cosmic-ray generated bubbles around their sources
• Constraining Cosmic-ray Transport with Observations of the Circumgalactic Medium
• The Impact of Cosmic Rays on the Kinematics of the Circumgalactic Medium
• Role of cosmic rays in the early stages of galactic outflows
• Cosmic rays from massive star clusters : A close look at Westerlund 1
• Steady states of the Parker instability
• Entropy-Conserving Scheme for Modeling Nonthermal Energies in Fluid Dynamics Simulations
• Capturing Turbulence with Numerical Dissipation: a Simple Dynamical Model for Unresolved Turbulence in Hydrodynamic Simulations
• IllustrisTNG + Cosmic Rays with a Simple Transport Model: From Dwarfs to L${}^{\star }$ Galaxies
• Comparing energy and entropy formulations for cosmic ray hydrodynamics
• Radio Halo of NGC 4631: Comparing Observations and Simulations