Using a viscosity matrix to construct a Riemann solve for the equations of special relativistic hydrodynamics
I.M. Kulikov
Institute of Computational Mathematics and Mathematical Geophysics of Siberian Branch of Russian Academy of Sciences, Novosibirsk, Russia
Keywords: numerical modeling, computational astrophysics, special relativistic magnetohydrodynamics
Abstract
Traditionally, to solve the hydrodynamic equations a Godunov method is used, whose main component is the solution of a Riemann problem to compute the fluxes of the conservative variables through the interfaces. Most numerical Riemann solvers are based on partial or full spectral decompositions of the Jacobian matrix with the spatial derivatives. However, when using complex hyperbolic models and various types of equations of state, even partial spectral decompositions are quite difficult to find analytically. Such hyperbolic systems include the equations of special relativistic magnetic hydrodynamics. In this paper, a numerical Riemann solver is constructed by means of a viscosity matrix on the basis of Chebyshev polynomials. This scheme does not require information about the spectral decomposition of the Jacobian matrix, while considering all types of waves in its design. To reduce the dissipation of the numerical solution, a piecewise parabolic reconstruction of the physical variables is used. The behavior of the numerical method is studied by using some classical test problems.
|
