A Reduced-Speed-of-Light Formulation of the Magnetohydrodynamic-Particle-in-Cell Method

Abstract

A reduced-speed-of-light (RSOL) approximation is a useful technique for magnetohydrodynamic (MHD)-particle-in-cell (PIC) simulations. With a RSOL, some "in-code" speed-of-light $\tilde{c}$ is set to much lower values than the true $c$, allowing simulations to take larger timesteps (which are restricted by the Courant condition given the large CR speeds). However, due to the absence of a well-formulated RSOL implementation from the literature, with naive substitution of the true $c$ with a RSOL, the CR properties in MHD-PIC simulations (e.g.\ CR energy or momentum density, gyro radius) vary artificially with respect to each other and with respect to the converged ($\tilde{c}\to c$) solutions with different choices of a RSOL. Here, we derive a new formulation of the MHD-PIC equations with a RSOL, and show that (1) it guarantees all steady-state properties of the CR distribution function and background plasma/MHD quantities are {\em independent} of the RSOL $\tilde{c}$ even for $\tilde{c}\ll c$, (2) ensures that the simulation can {\em simultaneously} represent the real physical values of CR number, mass, momentum, and energy density, (3) retains the correct physical meaning of various terms like the electric field, and (4) ensures the numerical timestep for CRs can always be safely increased by a factor $\sim c/\tilde{c}$. This new RSOL formulation should enable greater self-consistency and reduced CPU cost in simulations of CR-MHD interactions.

arXiv commentpublished by MNRAS

References in corpus (33)

• RAMSES-RT: Radiation hydrodynamics in the cosmological context
• Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients
• But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation
• Simulations of relativistic collisionless shocks: shock structure and particle acceleration
• Bayesian analysis of cosmic-ray propagation: evidence against homogeneous diffusion
• Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos
• Cosmic ray feedback in the FIRE simulations: constraining cosmic ray propagation with GeV gamma ray emission
• Models for Galactic cosmic-ray propagation
• Impact of Cosmic Ray Transport on Galactic Winds
• The Role of Cosmic Ray Transport in Shaping the Simulated Circumgalactic Medium
• Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks
• A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method
• Cosmic ray feedback heating of the intracluster medium
• Cosmic-ray hydrodynamics: Alfvén-wave regulated transport of cosmic rays
• The effects of cosmic rays on the formation of Milky Way-mass galaxies in a cosmological context
• Ion acceleration in non-relativistic astrophysical shocks
• Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies
• Role of Cosmic Rays in the Circumgalactic Medium
• The Impact of Cosmic Rays on Thermal Instability in the Circumgalactic Medium
• Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
• A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations
• Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast }$ Galaxies
• The failure of stellar feedback, magnetic fields, conduction, and morphological quenching in maintaining red galaxies
• On magnetic field amplification and particle acceleration near non-relativistic astrophysical shocks: Particles in MHD Cells simulations
• Magnetohydrodynamic-Particle-in-Cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution
• Virial shocks are suppressed in cosmic ray-dominated galaxy halos
• Heavy dust obscuration of z=7 galaxies in a cosmological hydrodynamic simulation
• A Consistent Reduced-Speed-of-Light Formulation of Cosmic Ray Transport Valid in Weak and Strong-Scattering Regimes
• Influence of Ion-Neutral Damping on the Cosmic-Ray Streaming Instability: Magnetohydrodynamic Particle-in-cell Simulations
• Towards First-principle Characterization of Cosmic-ray Transport Coefficients from Multi-scale Kinetic Simulations
• Chaotic cold accretion in giant elliptical galaxies heated by AGN cosmic rays
• Cosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments
• Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

Cited by in corpus (1)

• The Magnetohydrodynamic-Particle-In-Cell Module in Athena++: Implementation and Code Tests