January 27, 2006
An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The HLLC solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper (Mignone & Bodo 2004) to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave. Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the contrained tranport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one and two dimensional numerical test problems.
Similar papers 1
November 17, 2017
We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the solution is approximated by two constant states flanked by two shocks separated by a contact wave. The accuracy of the new approximate solver is calibrated through one- and two-dimensional test problems.
November 28, 2005
This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the tradit...
November 17, 2021
We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers recover the solution of the Riemann problem by solving a simplified or reduced set of jump conditions, whose level of complexity depends on the intermediate modes that are included. Five different approaches - namely the HLL, HLLC, HLLD, H...
January 18, 2011
We describe a new Godunov algorithm for relativistic magnetohydrodynamics (RMHD) that combines a simple, unsplit second order accurate integrator with the constrained transport (CT) method for enforcing the solenoidal constraint on the magnetic field. A variety of approximate Riemann solvers are implemented to compute the fluxes of the conserved variables. The methods are tested with a comprehensive suite of multidimensional problems. These tests have helped us develop a hier...
November 10, 2008
We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for the equations of ideal MHD. The solution to the Riemann problem is approximated by a five wave pattern, comprised of two outermost fast shocks, two rotational discontinuities and a contact surface in the middle. The proposed scheme is conside...
June 17, 2005
We present an extension of the HLLC approximate Riemann solver by Toro, Spruce and Speares to the relativistic equations of fluid dynamics. The solver retains the simplicity of the original two-wave formulation proposed by Harten, Lax and van Leer (HLL) but it restores the missing contact wave in the solution of the Riemann problem. The resulting numerical scheme is computationally efficient, robust and positively conservative. The performance of the new solver is evaluated t...
August 11, 2021
We propose a new Harten-Lax-van Leer discontinuities (HLLD) approximate Riemann solver to improve the stability of shocks and the accuracy of low-speed flows in multidimensional magnetohydrodynamic (MHD) simulations. Stringent benchmark tests verify that the new solver is more robust against a numerical shock instability and is more accurate for low-speed, nearly incompressible flows than the original solver, whereas additional computational costs are quite low. The novel abi...
October 29, 2002
A third order shock-capturing numerical scheme for three-dimensional special relativistic magnetohydrodynamics (3-D RMHD) is presented and validated against several numerical tests. The simple and efficient central scheme described in Paper I (Del Zanna and Bucciantini, Astron. Astrophys., 390, 1177--1186, 2002) for relativistic hydrodynamics is here extended to the magnetic case by following the strategies prescribed for classical MHD by Londrillo and Del Zanna (Astrophys. J...
February 1, 2016
The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts of the mesh. In this paper we present a reconstruction strategy that overcomes this problem by making a single conservative to primitive transformation per cell followed by higher order WENO reconstruction on a carefully chosen set of prim...
June 22, 2005
We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic case and we show that the problem can be understood as a purely relativistic hydr...