A dissipative extension to ideal hydrodynamics

Marcus J. Hatton,Ian Hawke

2024-05-16

Abstract:We present a formulation of special relativistic, dissipative hydrodynamics (SRDHD) derived from the well-established Müller- Israel-Stewart (MIS) formalism using an expansion in deviations from ideal behaviour. By re-summing the non-ideal terms, our approach extends the Euler equations of motion for an ideal fluid through a series of additional source terms that capture the effects of bulk viscosity, shear viscosity and heat flux. For efficiency these additional terms are built from purely spatial derivatives of the primitive fluid variables. The series expansion is parametrized by the dissipation strength and timescale coefficients, and is therefore rapidly convergent near the ideal limit. We show, using numerical simulations, that our model reproduces the dissipative fluid behaviour of other formulations. As our formulation is designed to avoid the numerical stiffness issues that arise in the traditional MIS formalism for fast relaxation timescales, it is roughly an order of magnitude faster than standard methods near the ideal limit.

General Relativity and Quantum Cosmology,High Energy Astrophysical Phenomena,Computational Physics,Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to effectively simulate hydrodynamic behaviors with dissipative effects within the framework of special relativity. Specifically, the author proposes a new method based on the Müller - Israel - Stewart (MIS) formalism to extend the ideal hydrodynamic equations so that they can capture the viscous and heat - conduction effects.

### Background and Motivation

1. **Binary Neutron Star Mergers**:
- Binary neutron star mergers are complex astrophysical phenomena, involving supranuclear - density matter, strong - gravitational spacetime, and the origin of heavy elements in the universe.
- These merger events will produce multi - messenger signals, including electromagnetic waves and gravitational waves, as confirmed by the GW170817 event.
- Observing these signals helps to constrain the mass, radius, and tidal deformability of neutron stars, thereby inferring their equations of state.

2. **The Need for Numerical Simulation**:
- Simulating these events requires highly nonlinear general - relativistic magnetohydrodynamics (GRMHD) models and combines numerical relativity for spacetime evolution.
- Usually, the fluid of neutron stars is regarded as "ideal", that is, the fluid stress is completely isotropic.
- However, recent studies have begun to focus on the influence of non - ideal fluid behaviors, such as the viscosity and heat - conduction effects generated by non - equilibrium processes, which are especially important shortly after neutron star mergers.

3. **Limitations of Existing Models**:
- The traditional MIS model will have numerical stiffness problems when dealing with fast relaxation time scales, which leads to inefficient simulations.
- When the dissipative relaxation time scale is small, the source terms become "stiff", requiring a significant reduction in the time step or the use of implicit time - integration methods, which greatly increases the computational cost.

### Main Contributions of the Paper

1. **Proposal of a New Model**:
- The author proposes a new dissipative - extension model (DEIFY), which extends the ideal hydrodynamic equations by resumming non - ideal terms.
- The new model captures the effects of bulk viscosity, shear viscosity, and heat flux through a series of additional source terms.

2. **Improvement of Numerical Efficiency**:
- By avoiding the numerical stiffness problem in the traditional MIS model, the new model is one order of magnitude faster in calculation speed than the standard method when approaching the ideal limit.
- Numerical simulation results show that the new model can reproduce the behaviors of other dissipative fluid models.

3. **Theoretical Basis**:
- This model is based on the Chapman - Enskog (CE) expansion method. By expanding the source terms of the MIS model, a new, non - rigid source - term form is obtained.
- Through this expansion method, the new model can significantly reduce the computational complexity while maintaining high accuracy.

### Conclusion

This paper solves the problem of efficiently simulating non - ideal fluid behaviors within the framework of special relativity by proposing a new dissipative - extension model. The new model not only has a solid theoretical foundation but also shows higher efficiency in numerical simulations, providing a powerful tool for future astrophysical research.