Approximate Riemann Solvers for the Cosmic Ray Magnetohydrodynamical Equations

Abstract

We analyze the cosmic-ray magnetohydrodynamic (CR MHD) equations to improve the numerical simulations. We propose to solve them in the fully conservation form, which is equivalent to the conventional CR MHD equations. In the fully conservation form, the CR energy equation is replaced with the CR "number" conservation, where the CR number density is defined as the three fourths power of the CR energy density. The former contains an extra source term, while latter does not. An approximate Riemann solver is derived from the CR MHD equations in the fully conservation form. Based on the analysis, we propose a numerical scheme of which solutions satisfy the Rankine-Hugoniot relation at any shock. We demonstrate that it reproduces the Riemann solution derived by Pfrommer et al. (2006) for a 1D CR hydrodynamic shock tube problem. We compare the solution with those obtained by solving the CR energy equation. The latter solutions deviate from the Riemann solution seriously, when the CR pressure dominates over the gas pressure in the post-shocked gas. The former solutions converge to the Riemann solution and are of the second order accuracy in space and time. Our numerical examples include an expansion of high pressure sphere in an magnetized medium. Fast and slow shocks are sharply resolved in the example. We also discuss possible extension of the CR MHD equations to evaluate the average CR energy.

arXiv commentAccepted to MNRAS, 16 pages, 14 figures, 7 tables

References in corpus (14)

• Solar Modulation of Cosmic Rays
• Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for One-Dimensional Flow
• Detecting shock waves in cosmological smoothed particle hydrodynamics simulations
• Cosmic-ray Driven Outflows in Global Galaxy Disk Models
• Cosmic ray feedback in hydrodynamical simulations of galaxy formation
• The Fermi Bubbles: Supersonic AGN Jets with Anisotropic Cosmic Ray Diffusion
• Cosmic rays can drive strong outflows from gas-rich high-redshift disk galaxies
• Nonthermal particles and photons in starburst regions and superbubbles
• Incorporation of cosmic ray transport into the ZEUS MHD code. Application for studies of Parker instability in the ISM
• Forecasts for the detection of the magnetised cosmic web from cosmological simulations
• Modelling injection and feedback of Cosmic Rays in grid-based cosmological simulations: effects on cluster outskirts
• An implicit scheme for solving the anisotropic diffusion of heat and cosmic rays in the RAMSES code
• Nonlinear Parker Instability with the Effect of Cosmic-Ray Diffusion
• Numerical simulations of buoyancy instabilities in galaxy cluster plasmas with cosmic rays and anisotropic thermal conduction

Cited by in corpus (5)

• Dynamics of the Parker-Jeans instability of gaseous disks including the effect of cosmic rays
• A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks
• Entropy-Conserving Scheme for Modeling Nonthermal Energies in Fluid Dynamics Simulations
• Fluid Simulations of Cosmic Ray Modified Shocks
• Comparing energy and entropy formulations for cosmic ray hydrodynamics