Abstract

The accurate modeling of various features in high energy astrophysical scenarios requires the solution of the Einstein equations together with those of special relativistic hydrodynamics (SRHD). Such models are more complicated than the non-relativistic ones due to the nonlinear relations between the conserved and state variables. A high-resolution shock-capturing central upwind scheme is implemented to solve the given set of equations. The proposed technique uses the precise information of local propagation speeds to avoid the excessive numerical diffusion. The second order accuracy of the scheme is obtained with the use of MUSCL-type initial reconstruction and Runge-Kutta time stepping method. After a discussion of the equations solved and of the techniques employed, a series of one and two-dimensional test problems are carried out. To validate the method and assess its accuracy, the staggered central and the kinetic flux-vector splitting schemes are also applied to the same model. The scheme is robust and efficient. Its results are comparable to those obtained from the sophisticated algorithms, even in the case of highly relativistic two-dimensional test problems.

Highlights

  • The special relativistic hydrodynamical models can be used to simulate many high energy phenomena in astrophysics, including accretion flows, gamma-ray bursts, and jet flows [1, 2]

  • Free-electron laser technology, high energy particles beams and heavy-ion collisions can be modeled by using special relativistic hydrodynamics (SRHD) [3]

  • After integrating the two-dimensional SRHD Eq (31) over the control volume Cij, the two-dimensional extension of the scheme can be expressed as dWi;j dt

Read more

Summary

Introduction

The special relativistic hydrodynamical models can be used to simulate many high energy phenomena in astrophysics, including accretion flows, gamma-ray bursts, and jet flows [1, 2]. Several upwind high-resolution shock-capturing (HRSC) schemes were being applied to solve relativistic hydrodynamical models. The method can be implemented as a black box solver to any system of conservation laws This family of central schemes suffers from excessive numerical viscosity when a sufficiently small time step is imposed, e.g., due to the presence of degenerate diffusion terms. The central upwind scheme reduces large amount of numerical dissipation present in the NT central schemes This scheme has already been applied to different problems namely, two-layer shallow water equations [22] and Hamilton Jacobi equations [23]. The central upwind scheme is implemented for solving the SRHD equations in one and two space dimensions.

Description of Special Relativistic Hydrodynamics
Recovery of physical variables
One-dimensional central upwind scheme
One-dimensional central scheme
Two-dimensional central upwind scheme
Two-dimensional central scheme
One-dimensional test problems
Two-Dimensional Tests
Findings
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call