A new HLLD Riemann solver with Boris correction for reducing Alfvén speed

Abstract

A new Riemann solver is presented for the ideal magnetohydrodynamics (MHD) equations with the so-called Boris correction. The Boris correction is applied to reduce wave speeds, avoiding an extremely small timestep in MHD simulations. The proposed Riemann solver, Boris-HLLD, is based on the HLLD solver. As done by the original HLLD solver, (1) the Boris-HLLD solver has four intermediate states in the Riemann fan when left and right states are given, (2) it resolves the contact discontinuity, Alfvén waves, and fast waves, and (3) it satisfies all the jump conditions across shock waves and discontinuities except for slow shock waves. The results of a shock tube problem indicate that the scheme with the Boris-HLLD solver captures contact discontinuities sharply and it exhibits shock waves without any overshoot when using the minmod limiter. The stability tests show that the scheme is stable when $|u|\lesssim 0.5c$ for a low Alfvén speed ($V_A\lesssim c$), where $u$, $c$, and $V_A$ denote the gas velocity, speed of light, and Alfvén speed, respectively. For a high Alfvén speed ($V_A\gtrsim c$), where the plasma beta is relatively low in many cases, the stable region is large, $|u|\lesssim (0.6-1)c$. We discuss the effect of the Boris correction on physical quantities using several test problems. The Boris-HLLD scheme can be useful for problems with supersonic flows in which regions with a very low plasma beta appear in the computational domain.

arXiv comment17 pages, 8 figures, Accepted by ApJ

References in corpus (11)

• The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics
• The Formation and Structure of a Strongly Magnetized Corona above Weakly Magnetized Accretion Disks
• Radiative MHD simulation of sunspot structure
• Local Study of Accretion Disks with a Strong Vertical Magnetic Field: Magnetorotational Instability and Disk Outflow
• Extension of the MURaM radiative MHD code for coronal simulations
• Directions of Outflows, Disks, Magnetic Fields, and Rotation of YSOs in Collapsing Molecular Cloud Cores
• Self-gravitational Magnetohydrodynamics with Adaptive Mesh Refinement for Protostellar Collapse
• Circumstellar disks and outflows in turbulent molecular cloud cores: possible formation mechanism for misaligned systems
• Collapse of Magnetized Singular Isothermal Toroids: I. Non-Rotating Case
• A Three-dimensional Simulation of a Magnetized Accretion Disk: Fast Funnel Accretion onto a Weakly Magnetized Star
• Global simulations of magnetorotational turbulence II: turbulent energetics

Cited by in corpus (3)

• Triple spiral arms of a triple protostar system imaged in molecular lines
• Non-linear evolutions of magnetised thick discs around black holes: dependence on the initial data
• Simulating a pulsed power-driven plasma with ideal MHD