A New Relativistic Hydrodynamic Code Dongsu Ryu, Indranil Chattopadhyay Department of Astronomy and Space Science, Chungnam National University, Daejeon 305-764, Korea: [email protected]nu.ac.kr, [email protected]u.ac.krand Eunwoo Choi Department of Physics and Astronomy, Georgia State University, P.O. Box 4106, Atlanta, GA 30302-4106, USA: [email protected]ABSTRACT Relativistic temperature of gas raises the issue of equation of state (EoS) in relat ivis tic hydrodyn amics. We present a code for relat ivistic hy drodyna mics with an EoS that is simple but approximates very closely the EoS of single- component perfect gas in the relativistic regime. Tests with a code based on the TVD scheme are presented to highlight differences induced by different EoS.
15
Embed
Dongsu Ryu and Indranil Chattopadhyay- A New Relativistic Hydrodynamic Code
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
8/3/2019 Dongsu Ryu and Indranil Chattopadhyay- A New Relativistic Hydrodynamic Code
Relativistic temperature of gas raises the issue of equation of state (EoS) inrelativistic hydrodynamics. We present a code for relativistic hydrodynamics
with an EoS that is simple but approximates very closely the EoS of single-
component perfect gas in the relativistic regime. Tests with a code based on the
TVD scheme are presented to highlight differences induced by different EoS.
8/3/2019 Dongsu Ryu and Indranil Chattopadhyay- A New Relativistic Hydrodynamic Code
where K 2 and K 3 are the modified Bessel functions of the second kind of order two and three,
respectively. In the extreme non-relativistic limit (Θ → 0), γ h → 5/2, and in the extreme
ultrarelativistic limit (Θ →∞), γ h → 4. However, using the above EoS comes with a price
of extra computational cost, since the thermodynamics of the fluid is expressed in terms of
the modified Bessel functions and no analytic expression can be written for LT.
In a recent paper, Mignone et al. (2005) proposed an EoS which fits RP well. The EoS,which is abbreviated as TM following Mignone et al. (2005), is given by
h =5
2Θ +
3
2
Θ2 +
4
9. (11)
With TM the expressions of n and cs become
n =3
2+
3
2
Θ Θ2 + 4/9
, c2s =5Θ
Θ2 + 4/9 + 3Θ2
12Θ
Θ2 + 4/9 + 12Θ2 + 2. (12)
TM was derived from the lower bound of the Taub’s inequality, (h − Θ)(h − 4Θ) = 1. Itproduces right asymptotic values for γ h.
In this poster we propose a new EoS, which is a simpler algebraic function of Θ and is
also a better fit of RP compared to TM. We abbreviate our proposed EoS as RC and give it
by
h = 26Θ2 + 4Θ + 1
3Θ + 2. (13)
With RC the expressions of n and cs become
n = 3
9Θ2 + 12Θ + 2
(3Θ + 2)2 , c2
s =
Θ(3Θ + 2)(18Θ2 + 24Θ + 5)
3(6Θ2 + 4Θ + 1)(9Θ2 + 12Θ + 2) . (14)
RC satisfies the Taub’s inequality, (h − Θ)(h − 4Θ) ≥ 1, for all Θ. It also produces right
asymptotic values for γ h. For both TM and RC, correctly c2s → 5Θ/3 in the extreme non-
relativistic limit, and c2s → 1/3 in the extreme ultrarelativistic limit, respectively.
In Figure 1, γ h, n, and cs are plotted with Θ to compare TM and RC to RP as well as
ID. One can see the RC is a much better fit of RP than TM with
|hTM − hRP|hRP
2%,|hRC− hRP|
hRP 0.8%. (15)
It is to be remembered that both γ h and n are independent of Θ, if ID is used.
8/3/2019 Dongsu Ryu and Indranil Chattopadhyay- A New Relativistic Hydrodynamic Code
The differences induced by different EoS are illustrated through a series of shock tube
tests, which were performed using the TVD code built with the EoS in §2 and the eigenvalues
and eigenvectors in §3. Two sets are considered.
For the first set with parallel velocity component only, two tests are presented:P1: ρL = 10, ρR = 1, pL = 13.3, pR = 10−6, and v p,L = v p,R = 0 initially, and tend = 0.45,
P2: ρL = ρR = 1, pL = 103, pR = 10−2, and v p,L = v p,R = 0 initially, and tend = 0.4.
For the second set with transverse velocity component, two tests, where different trans-
verse velocities were added to the test P2, are presented:
T1: initially vt,R = 0.99 to the right state, tend = 0.45,
T2: initially vt,L = 0.9 and vt,R = 0.99 to the left and right states, tend = 0.75.
The box covers the region of 0 ≤ x ≤ 1 in all the tests.
Figures 2, 3, 4, and 5 show the numerical solutions for RC and TM, but the analytic so-lutions for ID with γ = 5/3 and 4/3. For ID numerical solutions are almost indistinguishable
from analytic solutions, once they are calculated. The ID solutions are clearly different from
the RC and TM solutions. The ID solution with γ = 4/3 looks to match the RC and TM
solutions in P2, especially in the left region of contact discontinuity (hereafter CD) where
the flow is overall highly relativistic with Θ 1. But the difference is obvious in the region
between CD and shock, because Θ ∼ 1 there. On the other hand, the solutions of RC and
TM look very much alike. It reflects the similarity in the distributions of specific enthalpy
in equations (11) and (13). But yet there is a noticeable difference, especially in the density
in the region between CD and shock, and the difference reaches up to ∼ 5%.The most commonly used, ideal EoS, ID, can us used for entirely non-relativistic gas
(Θ 1) with γ = 5/3 or for entirely ultrarelativistic gas (Θ 1) with γ = 4/3. However,
if the transition from non-relativistic to relativistic with Θ ∼ 1 is involved, ID produces
incorrect results and using it should be avoid. The EoS proposed by Mignone et al. (2005),
TM, produces reasonably correct results with error of a few percent at most. The newly sug-
gested EoS, RC, which approximates the EoS of relativistic perfect gas, RP, most accurately,
produces thermodynamically the most accurate results. At the same time it is simple enough
to be implemented to numerical codes with minimum efforts and minimum computational
cost. The correctness and simplicity make RC suitable for astrophysical applications.
8/3/2019 Dongsu Ryu and Indranil Chattopadhyay- A New Relativistic Hydrodynamic Code