General-relativistic hydrodynamics of non-perfect fluids: 3+1 conservative formulation and application to viscous black-hole accretion

Michail Chabanov,Luciano Rezzolla,Dirk H. Rischke

DOI: https://doi.org/10.1093/mnras/stab1384

2023-02-25

Abstract:We consider the relativistic hydrodynamics of non-perfect fluids with the goal of determining a formulation that is suited for numerical integration in special-relativistic and general-relativistic scenarios. To this end, we review the various formulations of relativistic second-order dissipative hydrodynamics proposed so far and present in detail a particular formulation that is fully general, causal, and can be cast into a 3+1 flux-conservative form, as the one employed in modern numerical-relativity codes. As an example, we employ a variant of this formulation restricted to a relaxation-type equation for the bulk viscosity in the general-relativistic magnetohydrodynamics code $\texttt{BHAC}$. After adopting the formulation for a series of standard and non-standard tests in 1+1-dimensional special-relativistic hydrodynamics, we consider a novel general-relativistic scenario, namely, the stationary, spherically symmetric, viscous accretion onto a black hole. The newly developed solution $-$ which can exhibit even considerable deviations from the inviscid counterpart $-$ can be used as a testbed for numerical codes simulating non-perfect fluids on curved backgrounds.

General Relativity and Quantum Cosmology

What problem does this paper attempt to address?

This paper aims to solve some key problems in relativistic hydrodynamics of non - perfect fluids, especially to provide a suitable formula for numerical integration in special relativity and general relativity scenarios. Specifically, the goals of the paper include:

1. **Propose a comprehensive and causal formula**: The author reviews the existing formulas of relativistic second - order dissipative hydrodynamics and elaborates on one of the formulas that can be transformed into the 3 + 1 flux - conservation form. This form is compatible with the format used in modern numerical relativity codes.

2. **Apply to viscous black - hole accretion**: As an example, the author applies a variant of this formula to the bulk viscosity equation in the general - relativistic magnetohydrodynamics code BHAC. The validity of the formula is verified through a series of standard and non - standard tests.

3. **Study steady - state spherically - symmetric viscous accretion**: The author considers a new general - relativistic scenario, namely the accretion of steady - state, spherically - symmetric viscous matter onto a black hole. The newly developed solution shows significantly different characteristics from the non - viscous case and can serve as a test - bed for numerical codes simulating non - perfect fluids in a curvature background.

### Key problems

- **Numerical stability and causality**: Traditional first - order theories such as the Navier - Stokes equations have instabilities and non - causal behaviors in the relativistic framework. Therefore, this paper proposes a second - order theory to ensure its causality and stability.

- **Application background**: The study of extreme astrophysical phenomena such as binary neutron star mergers (BNS mergers) needs to consider non - ideal effects, such as bulk viscosity. These effects may significantly change the gravitational - wave signal and other observables after the merger.

### Summary of mathematical formulas

- **Energy - momentum tensor decomposition**:
\[
T^{\mu\nu}=e u^{\mu} u^{\nu}+(p + \Pi) h^{\mu\nu}+2 q^{(\mu} u^{\nu)}+\pi^{\mu\nu}
\]
where:
- \(e\) is the energy density.
- \(p\) is the isotropic pressure.
- \(\Pi\) is the bulk viscosity pressure.
- \(h^{\mu\nu}:=g^{\mu\nu}+u^{\mu} u^{\nu}\) is the projection operator orthogonal to \(u^{\mu}\).
- \(q^{\mu}\) is the heat flux, satisfying \(q^{\mu} u_{\mu} = 0\).
- \(\pi^{\mu\nu}\) is the shear - stress tensor, which is a symmetric, pure - spatial and traceless tensor.

- **Entropy flow**:
\[
S^{\mu}=s u^{\mu}+\frac{q^{\mu}}{T}-\left(\frac{\beta_{0} \Pi^{2}+\beta_{1} q^{\alpha} q_{\alpha}+\beta_{2} \pi^{\alpha\beta} \pi_{\alpha\beta}}{2T}\right) u^{\mu}+\frac{\alpha_{0} \Pi q^{\mu}+\alpha_{1} q^{\alpha} \pi_{\alpha}{}^{\mu}}{T}
\]

- **Dissipative current evolution equation**:
\[
\tau_{\Pi} \dot{\Pi}=\Pi_{NS}-\Pi-\frac{1}{2} \frac{\zeta}{T} \nabla_{\mu}(\tau_{\Pi} u^{\mu})+\alpha_{0} \zeta \nabla_{\mu} q^{\mu}+\gamma_{0} \zeta T q^{\mu} \nabla_{\mu} \left( \