No Arabic abstract
Relativistic temperature of gas raises the issue of the equation of state (EoS) in relativistic hydrodynamics. We study the EoS for numerical relativistic hydrodynamics, and propose a new EoS that is simple and yet approximates very closely the EoS of the single-component perfect gas in relativistic regime. We also discuss the calculation of primitive variables from conservative ones for the EoSs considered in the paper, and present the eigenstructure of relativistic hydrodynamics for a general EoS, in a way that they can be used to build numerical codes. Tests with a code based on the Total Variation Diminishing (TVD) scheme are presented to highlight the differences induced by different EoSs.
The characteristic formulation of the relativistic hydrodynamic equations (Donat et al 1998 J. Comput. Phys. 146 58), which has been implemented in many relativistic hydro-codes that make use of Godunov-type methods, has to be slightly modified in the case of evolving barotropic flows. For a barotropic equation of state, a removable singularity appears in one of the eigenvectors. The singularity can be avoided by means of a simple renormalization which makes the system of eigenvectors well defined and complete. An alternative strategy for the particular case of barotropic flows is discussed.
The dynamics of Quark-gluon plasma (QGP) as a lump of deconfined free quarks and gluons is elaborated. Based on the first principal we construct the Lagrangian that represents the dynamics of QGP. To induce a hydrodynamics approach, we substitute the gluon fields with flow fields. As a result, the derived equation of Motion (E.O.M) for gluon dominated QGP shows the form that similar to Euler equation, and the energy momentum tensor also represents explicitly the system of ideal fluid. Combining the E.O.M and energy momentum tensor, the pressure and energy density distribution as the equation of states are analytically derived.
We present the first numerical solutions of the causal, stable relativistic Navier-Stokes equations as formulated by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). For this initial investigation we restrict to plane-symmetric configurations of a conformal fluid in Minkowski spacetime. We consider evolution of three classes of initial data: a smooth (initially) stationary concentration of energy, a standard shock tube setup, and a smooth shockwave setup. We compare these solutions to those obtained with the Muller-Israel-Stewart (MIS) formalism, variants of which are the common tools used to model relativistic, viscous fluids. We find that for the two smooth initial data cases, simple finite difference methods are adequate to obtain stable, convergent solutions to the BDNK equations. For low viscosity, the MIS and BDNK evolutions show good agreement. At high viscosity the solutions begin to differ in regions with large gradients, and there the BDNK solutions can (as expected) exhibit violation of the weak energy condition. This behavior is transient, and the solutions evolve toward a hydrodynamic regime in a way reminiscent of an approach to a universal attractor. For the shockwave problem, we give evidence that if a hydrodynamic frame is chosen so that the maximum characteristic speed of the BDNK system is the speed of light (or larger), arbitrarily strong shockwaves are smoothly resolved. Regarding the shock tube problem, it is unclear whether discontinuous initial data is mathematically well-posed for the BDNK system, even in a weak sense. Nevertheless we attempt numerical solution, and then need to treat the perfect fluid terms using high-resolution shock-capturing methods. When such methods can successfully evolve the solution beyond the initial time, subsequent evolution agrees with corresponding MIS solutions, as well as the perfect fluid solution in the limit of zero viscosity.
This paper extends the second-order accurate BGK finite volume schemes for the ultra-relativistic flow simulations [5] to the 1D and 2D special relativistic hydrodynamics with the Synge equation of state. It is shown that such 2D schemes are very time-consuming due to the moment integrals (triple integrals) so that they are no longer practical. In view of this, the simplified BGK (sBGK) schemes are presented by removing some terms in the approximate nonequilibrium distribution at the cell interface for the BGK scheme without loss of accuracy. They are practical because the moment integrals of the approximate distribution can be reduced to the single integrals by some coordinate transformations. The relations between the left and right states of the shock wave, rarefaction wave, and contact discontinuity are also discussed, so that the exact solution of the 1D Riemann problem could be derived and used for the numerical comparisons. Several numerical experiments are conducted to demonstrate that the proposed gas-kinetic schemes are accurate and stable. A comparison of the sBGK schemes with the BGK scheme in one dimension shows that the former performs almost the same as the latter in terms of the accuracy and resolution, but is much more efficiency.
We study the conditions for the formation and propagation of Korteweg-de Vries (KdV) solitons in nuclear matter. In a previous work we have derived a KdV equation from Euler and continuity equations in non-relativistic hydrodynamics. In the present contribution we extend our formalism to relativistic fluids. We present results for a given equation of state, which is based on quantum hadrodynamics (QHD).