Next Article in Journal
Existence and Nonexistence of Positive Solutions for Perturbations of the Anisotropic Eigenvalue Problem
Next Article in Special Issue
Mapping GENERIC Hydrodynamics into Carter’s Multifluid Theory
Previous Article in Journal
Blow-Up Criterion and Persistence Property to a Generalized Camassa–Holm Equation
Previous Article in Special Issue
New Developments in Relativistic Magnetohydrodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Phenomenological Relativistic Second-Order Hydrodynamics for Multiflavor Fluids

by
Arus Harutyunyan
1,2,* and
Armen Sedrakian
3,4
1
Byurakan Astrophysical Observatory, National Academy of Sciences, Byurakan 0213, Armenia
2
Department of Physics, Yerevan State University, Yerevan 0025, Armenia
3
Frankfurt Institute for Advanced Studies, D-60438 Frankfurt am Main, Germany
4
Institute of Theoretical Physics, University of Wrocław, 50-204 Wrocław, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(2), 494; https://doi.org/10.3390/sym15020494
Submission received: 27 December 2022 / Revised: 3 February 2023 / Accepted: 9 February 2023 / Published: 13 February 2023

Abstract

:
In this work, we perform a phenomenological derivation of the first- and second-order relativistic hydrodynamics of dissipative fluids. To set the stage, we start with a review of the ideal relativistic hydrodynamics from energy–momentum and particle number conservation equations. We then go on to discuss the matching conditions to local thermodynamical equilibrium, symmetries of the energy–momentum tensor, decomposition of dissipative processes according to their Lorentz structure, and, finally, the definition of the fluid velocity in the Landau and Eckart frames. With this preparatory work, we first formulate the first-order (Navier–Stokes) relativistic hydrodynamics from the entropy flow equation, keeping only the first-order gradients of thermodynamical forces. A generalized form of diffusion terms is found with a matrix of diffusion coefficients describing the relative diffusion between various flavors. The procedure of finding the dissipative terms is then extended to the second order to obtain the most general form of dissipative function for multiflavor systems up to the second order in dissipative fluxes. The dissipative function now includes in addition to the usual second-order transport coefficients of Israel–Stewart theory also second-order diffusion between different flavors. The relaxation-type equations of second-order hydrodynamics are found from the requirement of positivity of the dissipation function, which features the finite relaxation times of various dissipative processes that guarantee the causality and stability of the fluid dynamics. These equations contain a complete set of nonlinear terms in the thermodynamic gradients and dissipative fluxes arising from the entropy current, which are not present in the conventional Israel–Stewart theory.

1. Introduction

Relativistic hydrodynamics has been widely applied in recent years in heavy-ion physics [1,2,3,4] and astrophysics [5,6,7]. For example, the experiments at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) motivated significantly the development of the relativistic hydrodynamics with the data being well-described in terms of a fluid with low shear viscosity; for reviews, see [8,9]. Relativistic hydrodynamics is also used in the general relativistic simulations of compact stars in isolation or binaries, particularly in the context of gravitational wave emission by neutron star mergers observed in 2017 [10]. A further astrophysical area of applications of relativistic hydrodynamics is the physics of compact star rotational dynamics, specifically glitches and their relaxations; for reviews, see [11,12].
The hydrodynamic description of fluids is valid close to the local thermal equilibrium. The hydrodynamic state of a relativistic fluid is described through its energy–momentum tensor and currents of conserved charges, which in the low-frequency and long-wavelength can be Taylor-expanded around their equilibrium values in thermodynamic gradients (so-called thermodynamic forces). The validity of such gradient expansion is guaranteed due to the clear separation between the typical microscopic and macroscopic scales of the system. The zeroth-order term in this expansion corresponds to the limit of the ideal hydrodynamics.
The truncation of the gradient expansion at the first order leads to the relativistic Navier–Stokes (NS) theory, which was worked out by Eckart [13] and Landau–Lifshitz [14]. It is known that the solutions of the relativistic NS equations are acausal and unstable [15,16,17,18]. The reason for the acausality is the parabolic structure of NS equations, which originated from the linear constitutive relations between the dissipative fluxes and the thermodynamic forces. Recent work demonstrated that the acausalities and instabilities in relativistic hydrodynamics are a consequence of the matching procedure to the local equilibrium reference state. More general matching conditions were used to render the theory causal and stable at first order [19,20,21].
The problem of acausality can be solved in the second-order theory, where additional terms appear that contain higher (second)-order derivatives in thermodynamic quantities. For non-relativistic fluids, the second-order theory was proposed by Müller [22], and then rediscovered and extended to relativistic systems by Israel and Stewart [23,24]. In these theories, the dissipative fluxes are treated as independent state variables that satisfy relaxation-type equations derivable from the entropy principle. The relaxation terms which appear in these equations recover the causality of the theory [18,25]. The relaxation equations for dissipative fluxes in the second-order theories and their numerical studies in simulations have been discussed extensively in the literature; see [26,27,28,29,30,31,32,33,34,35,36,37].
The relativistic second-order hydrodynamics can be obtained from moments of the Boltzmann equation for the distribution function [38,39,40]. This theory provides a systematic way of evaluating the new coefficients describing the relaxation of dissipative quantities at weak coupling and in the quasiparticle limit. An alternative approach valid in the strong coupling limit is Zubarev’s non-equilibrium statistical operator formalism, which was recently applied to obtain the second-order hydrodynamics and/or Kubo-type formulae for transport coefficients [41,42,43]. Related approaches based on field theory methods are given in Refs. [44,45].
Recent applications of relativistic hydrodynamics focused on the systems featuring spin, polarization, and vorticity. A systematic calculation of the corrections of the stress–energy tensor and currents up to the second order in thermal vorticity were given in Refs. [46,47] and the relevant Kubo formulae were derived. Similarly, the problem of relativistic hydrodynamics under strong spin was studied on the basis of entropy–current analysis in Ref. [48] where the seventeen transport coefficients of highly anisotropic relativistic hydrodynamics were identified. Ref. [49] formulated relativistic spin hydrodynamics of Dirac fermions in a torsionful curved background and derived the relevant Kubo formulae associated with the correlation functions. The spin hydrodynamics was also derived from the Boltzmann equation using the method of moments in Ref. [50] up to the second order. Plasmas with vorticity were studied in second-order dissipative hydrodynamics in the presence of a chiral imbalance in Ref. [51] and the dispersion relations of chiral vortical waves were obtained.
Specific applications of second-order hydrodynamics include, for example, the computation of the transport coefficients from a Chapman–Enskog-like expansion for a system of massless quarks and gluons [34]. The second-order hydrodynamics was recently applied to multicomponent systems with hard-sphere interactions in Ref. [52], where the multicomponent nature of the system was taken into account. The second-order hydrodynamics was applied to quark–gluon plasma in heavy ion collisions, for example, in Ref. [53], showing the coupling between diffusion and shear and bulk viscosities in a multi-component system. A discussion of the use of relativistic anisotropic hydrodynamics to study the physics of ultrarelativistic heavy-ion collisions was given in Ref. [54], using as the main ingredient the quasiparticle anisotropic hydrodynamics model for quark–gluon plasma. Anomalies in the field theory were included in the description of heavy-ion collisions in Ref. [55]. The strong gravity regime of second-order hydrodynamics, relevant for astrophysical applications, was studied in Ref. [56]. The first-order theory was applied in the context of the out-of-equilibrium dynamics of viscous fluids in a spatially flat cosmology [57].
Recent work also extended the current formulations of relativistic hydrodynamics to include the effect of magnetic fields. The first-order dissipative effects and relevant transport coefficients were derived, and a numerical method was developed using the method of moments in Ref. [58]. The effects of anomalies in the magnetohydrodynamics were included in Ref. [59]. Charge diffusion in the second-order theories was recently studied in Ref. [60]. Applications to neutron stars and relevant transport coefficients were derived, for example, in Refs. [61,62,63]. A comprehensive review of the new developments in this field is given in Ref. [64].
More formal recent developments of relativistic hydrodynamics include formulations that avoid specific frames. (The choice of frames in relativistic heavy ion collisions was discussed in Ref. [65].) First- and second-order viscous relativistic hydrodynamics were formulated without specific frame conditions in Refs. [21,66], respectively, with the causality and stability conditions explicitly formulated in the conformal regime in the second work. Earlier, the stability of first-order hydrodynamics was established in Ref. [20].
The purpose of this article is twofold. First, we review the phenomenological derivation of the second-order dissipative hydrodynamics from general principles of conservation laws and the second law of thermodynamics. The dissipative equations are obtained from expansions of the energy–momentum tensor and currents around their equilibrium values in thermodynamic gradients (so-called thermodynamic forces) order by order. Second, we consider the theory in the general case of a system with l independent flavors of conserved charges, keeping all possible second-order terms which arise from the entropy current and the second law of thermodynamics. In the case where the nonlinear terms in the thermodynamic gradients and dissipative fluxes are dropped, these equations are reduced to the original Israel and Stewart hydrodynamics [23,24].
This work is organized as follows. In Section 2, we present the ideal hydrodynamics of non-dissipative fluids. Dissipation in the fluids is introduced in Section 3. We discuss the first-order hydrodynamics and obtain the Navier–Stokes equations in Section 4. The second-order Israel–Stewart theory in the case of a system with l independent flavors of conserved charges is given in Section 5. A summary is provided in Section 6.

2. Relativistic Ideal Hydrodynamics

Relativistic hydrodynamics describes the state of the fluid employing its energy–momentum tensor T μ ν and currents of conserved charges N a μ , such as the baryonic, electric, etc. Here, we consider the general case of a system with l independent flavors of conserved charges, which are labeled by the index a. The equations of relativistic hydrodynamics are contained in the conservation laws for the energy–momentum tensor and the charge currents
μ T μ ν = 0 , μ N a μ = 0 , a = 1 , 2 , , l .
In dissipative hydrodynamics, also the entropy principle μ S μ 0 should be applied to close the system (1), where S μ is the entropy 4-current.
Ideal hydrodynamics corresponds to the zeroth-order expansion of the energy-momentum tensor and charge currents with respect to the thermodynamic forces. In this case, each fluid element maintains the local thermal equilibrium during its evolution [14,67]. The macroscopic state of the fluid is therefore fully described through the fields of the energy density ϵ ( x ) and the charge densities n a ( x ) . The fluid 4-velocity is defined as
u μ ( x ) = d x μ d τ = γ ( x ) , γ v i ( x ) ,
where x μ is the displacement 4-vector of a fluid element, d τ = g μ ν d x μ d x ν is the infinitesimal change in the proper time, v ( x ) is the fluid 3-velocity and γ ( x ) = ( 1 v 2 ) 1 / 2 is the relevant Lorentz factor. The 4-velocity given by Equation (2) is normalized by the condition u μ ( x ) u μ ( x ) = 1 and has only three independent components.
Because the thermal equilibrium is maintained locally, each fluid element can be assigned well-defined local values of the temperature T ( x ) , chemical potentials μ a ( x ) (conjugate to charge densities n a ( x ) ), entropy density s ( x ) and pressure p ( x ) . These quantities are related to the local energy and charge densities via an equation of state p = p ( ϵ , n a ) and the standard thermodynamic relations
T d s = d ϵ a μ a d n a , ( first law )
d p = s d T + a n a d μ a , ( Gibbs Duhem relation )
h = ϵ + p = T s + a μ a n a ,
where h is the enthalpy density.
The fluid 4-velocity is defined in such a way that in the fluid rest frame, the energy and charge flows vanish: N i = 0 and T 0 i = 0 , if u i = 0 . These conditions together with the spatial isotropy imply the following form for the energy–momentum tensor and charge currents [14,67]:
T 0 μ ν = ϵ u μ u ν p Δ μ ν , N a 0 μ = n a u μ ,
where the index 0 labels the quantities in ideal hydrodynamics; Δ μ ν = g μ ν u μ u ν is the projection operator onto the 3-space orthogonal to u μ and has the properties
u μ Δ μ ν = Δ μ ν u ν = 0 , Δ μ ν Δ ν λ = Δ λ μ , Δ μ μ = 3 .
From Equations (6) and (7), we obtain
ϵ = u μ u ν T 0 μ ν , p = 1 3 Δ μ ν T 0 μ ν , n a = u μ N a 0 μ .
In the fluid rest frame, u μ = ( 1 , 0 , 0 , 0 ) ; therefore, Δ μ ν = diag ( 0 , 1 , 1 , 1 ) and Δ i j = Δ i j = δ i j . In this case, Equation (8) simplifies to
ϵ = T 0 00 , p = 1 3 T 0 i i , n a = N a 0 0 .
We introduce also the energy 4-current or momentum density by the formula
P 0 μ = u ν T 0 μ ν + p u μ = h u μ ,
which is the relativistic generalization of the 3-momentum density (mass current). As seen from Equations (6) and (10), all charge currents are parallel to each other and to the energy flow, which is due to the possibility of a unique definition of the velocity field for ideal fluids.
The equations of ideal hydrodynamics are obtained by substituting the expressions (6) into the conservation laws (1)
μ N a 0 μ = μ ( n a u μ ) = 0 , μ T 0 μ ν = μ ( h u μ u ν p g μ ν ) = 0 .
Contracting the second equation in (11) once with u ν and then with the projector Δ α ν and taking into account Equation (7) and the condition u ν μ u ν = 0 , we obtain
D n a + n a θ = 0 , D ϵ + h θ = 0 , h D u α = α p ,
where we introduce the covariant time derivative, covariant spatial derivative, and velocity 4-divergence via D = u μ μ , α = Δ α β β and θ = μ u μ , respectively. The velocity 4-divergence θ quantifies how fast the fluid is expanding ( θ > 0 ) or contracting ( θ < 0 ); it is called the fluid expansion rate. In the case of an incompressible flow of the fluid, we have θ = 0 .
It is not difficult to recognize in the first two relations in (12) the covariant expressions for the charge conservation law and the energy conservation law, respectively. The third equation is nothing more than the relativistic generalization of the ordinary Euler equation familiar from nonrelativistic hydrodynamics. From these equations, one can deduce that in relativistic hydrodynamics, the role of rest mass density is taken over by the enthalpy density h, which thus provides the correct inertia measure for relativistic fluids. This fact illustrates the importance of the quantity (10) as the relativistic analog of the momentum flux.
The system (12) contains l + 4 equations for l + 5 variables ϵ , p, n a and u μ . To close the system, one still needs to specify an equation of state p = p ( ϵ , n a ) , which relates the pressure to the conserved thermodynamic variables.
It is easy to show that the equations of ideal hydrodynamics lead automatically to entropy conservation. The entropy flux can be written as
S 0 μ = s u μ .
Using the thermodynamic relations (3)–(5) and the equations of motion (12), we obtain
T μ S 0 μ = D s + s θ = D ϵ + h θ a μ a ( n a θ + D n a ) = 0 ,
which is the second law of thermodynamics for a non-dissipative system. Using Equations (5), (6), (10) and (13), we can rewrite the entropy current as
S 0 μ = β P 0 μ a α a N a 0 μ = p β μ + β ν T 0 μ ν a α a N a 0 μ ,
where we defined
β = T 1 , β ν = β u ν , α a = β μ a .
The expression (15) is the covariant form of the relation (5). To proceed further, it is convenient to modify Equations (3) and (4) using the definitions (16). We obtain for Equation (4)
d p = β 1 h d β + β 1 a n a d α a ,
where we used Equation (5) in the second step. Now the first law of thermodynamics and the Gibbs–Duhem relation can be written in an alternative form:
d s = β d ϵ a α a d n a , β d p = h d β + a n a d α a .
With the aid of Equations (6) and (15), these relations can be cast into a covariant form:
d S 0 μ = β ν d T 0 μ ν a α a d N a 0 μ ,
d ( p β μ ) = T 0 μ ν d β ν + a N a 0 μ d α a ,
where we used the second relation in Equation (18). One should note that, despite their vector form, these equations do not contain more information than the scalar thermodynamic relations (contraction of Equations (19) and (20) with the projector Δ μ ν leads to identities).

3. Relativistic Dissipative Hydrodynamics

3.1. Matching Conditions

Ideal hydrodynamics relies on the strong assumption of the local thermodynamic equilibrium of each fluid element. However, the local equilibrium cannot be maintained permanently since all real systems react to the non-uniformities of the fluid on finite time scales by generating irreversible (dissipative) fluxes. These fluxes tend to eliminate the local gradients to drive the system toward global thermal equilibrium, leading thereby to energy dissipation and entropy increase.
As in the case of ideal hydrodynamics, the equations of relativistic dissipative hydrodynamics can be obtained from the conservation laws given by Equation (1). The expressions (6) then need to be generalized by taking into account all possible effects of dissipation. The energy–momentum tensor obtains an anisotropic contribution because of irreversible momentum exchange between different fluid elements, as well as heat transfer due to the relative motion of energy and charge currents. The diffusion processes between different charge species in their turn introduce dissipative terms in the charge currents N a μ . As a result, the energy–momentum tensor and the charge currents for a dissipative fluid can be written in the following form:
T μ ν = T 0 μ ν + τ μ ν = ϵ u μ u ν p Δ μ ν + τ μ ν ,
N a μ = N a 0 μ + j a μ = n a u μ + j a μ ,
where j a μ and τ μ ν are dissipative terms, and the tensor τ μ ν is symmetric. Note that since the system is out of thermodynamic equilibrium, the thermodynamic parameters are not well-defined anymore, and one needs to impose additional conditions to specify what should be exactly understood under ϵ , n a and p in Equations (21) and (22). The thermodynamic variables in non-equilibrium states can be defined only by means of a fictitious equilibrium state, which should be constructed in such a way as to satisfy the thermodynamic relations (3)–(5) [24]. In order to construct such an equilibrium state for the given values of T μ ν and N a μ , we define first the energy and charge densities via the so-called matching (fitting) conditions:
ϵ ( u ) = u μ u ν T μ ν , n a ( u ) = u μ N a μ .
Equations (23) imply simply that ϵ and n a are the time-like eigenvalues of the energy–momentum tensor and the charge currents, respectively, measured by a local observer comoving with the fluid element. We remark that the quantities ϵ and n a depend in general on the choice of u μ , which is the consequence of the ambiguity of the definition of the velocity field for dissipative fluids (see Section 3.3 for details).
We define in the next step an equilibrium entropy density via an equation of state s = s ( ϵ , n a ) , which is chosen to be the same function of the parameters ϵ and n a , as it would be in full thermodynamic equilibrium [15,23,24,68]. The rest of the local thermodynamic quantities can be defined through standard thermodynamic relations (3)–(5), i.e.,
β = s ϵ n a , μ a = T s n a ϵ , n b n a , p = T s + a μ a n a ϵ ,
which implies that all thermodynamic parameters are defined in non-equilibrium states via the same functions of ϵ and n a as their equilibrium counterparts.
The thermodynamic parameters defined above form a fictitious equilibrium state for which the reversible thermodynamic relations are formally valid. However, it is worthwhile to stress that only the energy and charge densities can be ascribed a definite physical meaning, whereas the quantities s, p, T, and μ a do not retain their usual physical meaning when the system is out of thermal equilibrium. These quantities are mathematically convenient to use since they have the approximate physical meaning of their equilibrium counterparts for small departures from the local equilibrium [23,24].
For instance, the quantity p is not the actual thermodynamic pressure (i.e., the work done by a unit change of the fluid volume), but differs from the latter by an additional non-equilibrium term, which is of the first order in velocity gradients [14,24]. Similarly, the quantity s is not the actual non-equilibrium entropy density since it is defined through reversible thermodynamic relations. In other words, s is not the proper quantity which should increase in non-equilibrium processes according to the second law of thermodynamics. However, the first non-equilibrium correction, in this case, appears only at the second order in gradients, and, therefore, can be ignored in the first-order theory [14,23,24]. Indeed, the correction to s should be negatively defined since the entropy attains its maximum in equilibrium. The only scalar quantity which might contribute to the entropy density at the first order in gradients is the velocity 4-divergence θ ; therefore, the non-equilibrium entropy density should be s = s ( ϵ , n a ) a θ , with some thermodynamic coefficient a. Since θ can be of both signs, we conclude that a = 0 .

3.2. Decomposition in Different Dissipative Processes

To identify different dissipative processes, it is convenient to separate scalar, vector, and traceless parts of the tensor τ μ ν . We first note that the matching conditions (23) together with Equations (21) and (22) impose the following orthogonality conditions on the dissipative terms:
u μ u ν τ μ ν = 0 , u μ j a μ = 0 .
The tensor τ μ ν can be further decomposed into its irreducible components parallel and orthogonal to the fluid 4-velocity u μ . For that purpose, it is useful to introduce a fourth-rank traceless projector orthogonal to u μ via
Δ μ ν ρ σ = 1 2 Δ μ ρ Δ ν σ + Δ μ σ Δ ν ρ 1 3 Δ μ ν Δ ρ σ ,
which has the properties
Δ μ ν ρ σ = Δ ν μ ρ σ = Δ ρ σ μ ν , u μ Δ μ ν ρ σ = 0 , Δ α μ Δ μ ν ρ σ = Δ α ν ρ σ , Δ μ ρ σ μ = 0 , Δ ν μ σ μ = 5 3 Δ ν σ , Δ μ ν μ ν = 5 , Δ μ ν ρ σ Δ α β ρ σ = Δ μ ν α β .
The most general tensor decomposition of τ μ ν consistent with the orthogonality condition (25) can be now written as
τ μ ν = Π Δ μ ν + q μ u ν + q ν u μ + π μ ν ,
where we defined new dissipative quantities via
Π = 1 3 Δ μ ν τ μ ν = 1 3 τ μ μ , q μ = u α Δ β μ τ α β , π μ ν = Δ α β μ ν τ α β ,
and used the properties (7). From Equations (21), (22) and (28), we obtain the most general decompositions of the energy–momentum tensor and the charge currents into their irreducible components:
T μ ν = ϵ u μ u ν ( p + Π ) Δ μ ν + q μ u ν + q ν u μ + π μ ν ,
N a μ = n a u μ + j a μ .
The energy flow defined in Equation (10) can be generalized for dissipative fluids as
P μ = u ν T μ ν + p u μ = h u μ + q μ .
The dissipative terms j a μ , q μ , π μ ν and Π in Equations (30) and (31) are called charge diffusion fluxes, energy diffusion flux, shear stress tensor, and bulk viscous pressure, respectively. The shear stress tensor is the traceless spatial part of the energy–momentum tensor and describes the dissipation by anisotropic momentum flow, whereas the bulk viscous pressure is the non-equilibrium part of the pressure and is responsible for dissipation during isotropic expansion or compression. The dissipative currents satisfy the following conditions:
u ν q ν = 0 , u ν j a ν = 0 , u ν π μ ν = 0 ; π μ μ = 0 ,
the first three of which reflect the fact that the dissipation in the fluid should be spatial.
Finally, all quantities on the right-hand sides of Equations (21) and (22) can be obtained by the relevant projections of T μ ν and N a μ :
ϵ = u μ u ν T μ ν , n a = u μ N a μ , p + Π = 1 3 Δ μ ν T μ ν ,
π μ ν = Δ α β μ ν T α β , q μ = u α Δ β μ T α β , j a ν = Δ μ ν N a μ ,
as follows from Equations (21)–(23), (29) and the properties (7) and (27). In the fluid rest frame, we have
ϵ = T 00 , n a = N a 0 , p + Π = 1 3 T k k ,
π k l = δ k i δ l j 1 3 δ k l δ i j T i j , q i = T 0 i , j a i = N a i .
From these expressions, it can be seen that all dissipative quantities are, as expected, purely spatial in the rest state of the fluid. Each of the vectors q μ and j a μ thus has 3 independent components, and the shear stress tensor π μ ν has 5 independent components since its trace vanishes. Thus, the total number of independent quantities in Equations (36) and (37) is 4 l + 10 , which corresponds to the degrees of freedom of the energy–momentum tensor and the charge currents (remembering that p is determined by an equation of state). As for the velocity of the fluid, u μ should not be treated as an independent variable but should be associated with one of the physical currents. We devote the next subsection to the discussion of the possible definitions of the velocity field.

3.3. Definition of Flow Velocity

Another important question in dissipative fluid dynamics is the proper definition of fluid velocity. The choice of the frame in the case of ideal hydrodynamics is simple, as it has the energy and charge current flowing parallel to each other. Then, the fluid rest frame is defined via the requirement that these currents vanish identically. The situation is different in the case of dissipative fluids because one is faced with simultaneous flows of energy and particle currents. There are two simple and also natural ways to define the fluid rest state in dissipative hydrodynamics, which we describe in this subsection in turn.
In the Landau frame (or L-frame) [14] the net energy flow vanishes, the definition of fluid velocity u μ in this frame is then given through the time-like eigenvector of T μ ν :
u L μ = u L ν T μ ν u L ν T μ ν u L λ T μ λ ,
which in combination with Equations (23), (30) and (33), leads to the following relations:
ϵ L = u L ν T μ ν u L λ T μ λ , u L ν T μ ν = ϵ L u L μ , q L μ = 0 ,
where the subscript L indicates that the quantities are evaluated according to the Landau definition of u μ . With such a choice of the velocity field, the energy diffusion flux is zero, whereas the heat transport is accommodated in the particle diffusion fluxes j L a μ . A different form of Equation (38) is obtained by noting that the energy flux (32) in the L-frame is parallel to the flow velocity:
u L μ = P μ P μ P μ .
Next, consider the velocity of fluid u μ in a generic frame. It can be related to u L μ upon noticing that in a generic fluid rest frame the current (32) is given by P μ = ( h , q i ) and, therefore, the boost velocity from an arbitrary rest state ( u i = 0 ) to the Landau rest frame ( u L i = 0 ) is given by v L i = q i / h = O 1 . We follow the convention of Refs. [23,24] to denote the n-th order quantities in deviation from equilibrium by the symbol O n . The transformation of the charge currents into the L-frame is then given by
N a i = N a i v L i N a 0 + O 2 .
Setting here N a i = j a i , N a i = j L a i and N a 0 = n a , we find the charge diffusion fluxes in the Landau rest state:
j L a i = j a i n a h q i + O 2 .
It is worthwhile to note that the current on the left-hand side j L a i is evaluated at a transformed coordinate denoted by x . However, the difference j L a i ( x ) j L a i ( x ) ( x x ) j L a i v L j L a i can be neglected, as it is of the third order in gradients. Because in the fluid rest frame j L a 0 = j a 0 = q 0 = 0 , we can express Equation (42) in a covariant form as
j L a μ = j a μ n a h q μ + O 2 ,
which is valid already in an arbitrary frame, i.e., not exclusively in the fluid rest frame. The 4-currents
J a μ = j a μ n a h q μ N a μ n a h P μ
are the charge diffusion fluxes with respect to the energy flow, i.e., these are the charge currents in the absence of energy flow. Note that the combinations (44) are invariant under first-order changes in u μ despite the fact that the energy-diffusion flux q μ and the charge-diffusion fluxes j a μ depend on the chosen velocity field. It will be demonstrated below (see in Section 4), that the energy dissipation in irreversible processes is associated precisely with these currents.
The remaining thermodynamic variables which appear in Equations (30) and (31) remain unchanged if we neglect the second-order changes in thermodynamic gradients induced by u μ [23,24]. Thus, we can summarize the relation between u and u L as
u L μ = u μ + q μ h + O 2 ,
which follows straightforwardly upon comparing Equations (32) and (40).
The Eckart frame (E-frame) is defined as the frame in which the velocity field is parallel to one of the conserved currents N a μ . For a fluid with a single conserved charge (for example, the net particle number), one has N μ = n u μ + j μ and, therefore, the 4-velocity is defined as [13]
u E μ = N μ N μ N μ ,
which in combination with Equations (23) and (31) gives
n E = N μ N μ , N μ = n E u E μ , j E μ = 0 ,
i.e., the particle diffusion flux vanishes. The subscript E in Equations (46) and (47) indicates that the quantities are evaluated in the E-frame.
The boost velocity from an arbitrary rest frame to the Eckart rest frame is v E i = j i / n . This implies that the velocities u μ and u E μ can be related by
u E μ = u μ + j μ n + O 2 .
We are now in a position to transform the energy flux into the Eckart rest frame; we find
P i = P i j i n P 0 + O 2 ,
It follows then that the energy diffusion flux in the E-frame reads
q E μ = q μ h n j μ + O 2 .
The quantity
h μ = q μ h n j μ = h n J μ
can be seen as the energy flow with respect to the particle flow; therefore, it can be identified with heat flux. The relation (51) demonstrates that heat conduction and particle diffusion correspond to the same phenomenon, which is observed from different reference frames. This is true only when the higher than the first-order deviations from equilibrium can be neglected. To order O 1 , the relation between the L-frame and the E-frame can be established from Equations (45) and (48) as
u L μ u E μ = h μ h = J μ n .
Finally, let us note that the extension of the discussion above to the case where multiple conserved charges are present is straightforward. This is achieved by attaching a reference frame to each of them, which formally then leads to the definition
u a μ = N a μ N a μ N a μ ,
which, in turn, enforces the vanishing of the corresponding diffusion flux, i.e., j a μ = 0 . In what follows below, we will keep the fluid velocity generic (if not stated otherwise) without specifying any particular reference frame.

3.4. Equations of Relativistic Dissipative Hydrodynamics

Equations of relativistic dissipative hydrodynamics are obtained by substituting the decompositions (30) and (31) into the conservation laws (1). Using the same technique as in Section 2 and recalling the properties (33), we obtain
D n a + n a θ + μ j a μ = 0 ,
D ϵ + ( h + Π ) θ + μ q μ q μ D u μ π μ ν σ μ ν = 0 ,
( h + Π ) D u α α ( p + Π ) + Δ α μ D q μ + q μ μ u α + q α θ + Δ α ν μ π μ ν = 0 ,
where we introduced the velocity stress tensor as σ μ ν = Δ μ ν α β α u β .
The system of Equations (54)–(56) contains l + 4 equations, as in the case of ideal hydrodynamics, while the number of independent variables is now 4 l + 10 . The additional unknown variables here are 3 l components of the diffusion fluxes, and 5 components of the shear stress tensor and the bulk viscous pressure (recall that the equilibrium pressure is given by the equation of state, and one of the diffusion fluxes can always be eliminated). Thus, in order to solve the system of Equations (54)–(56), we need additional 3 l + 6 equations for these dissipative quantities. These relations are described in the phenomenological theory via the second law of thermodynamics (the entropy principle).

4. Relativistic Navier-Stokes (First-Order) Theory

The entropy conservation law (14) derived in the framework of ideal fluid dynamics is no longer valid for dissipative fluids. In this case, it should be replaced by the second law of thermodynamics, which implies that the entropy production rate of an isolated system must be always non-negative:
μ S μ 0 ,
where equality holds only for reversible processes. By analogy with the decomposition of the charge currents (31), we can decompose S μ into contributions parallel and orthogonal to u μ
S 1 μ = s u μ + j s μ ,
where s is identified with the equilibrium entropy density (which is sufficient for the first-order accuracy, as explained in Section 3.1), and the vector j s μ satisfies the condition u μ j s μ = 0 . The index 1 denotes that (58) is only the first-order approximation to the entropy flux. For small departures from equilibrium, it is natural to assume that j s μ is a linear combination of the energy and charge diffusion fluxes
j s μ = β ˜ q μ a α ˜ a j a μ ,
where β ˜ and α ˜ a are functions of thermodynamic variables and should be determined from the condition (57). This formulation of relativistic dissipative fluid dynamics was proposed by Eckart [13] and Landau–Lifshitz [14], and leads to the relativistic version of the NS theory. Inserting Equation (59) into Equation (58) and introducing the dissipation function via R T μ S μ , we obtain
R = Π θ + π μ ν σ μ ν + q μ ( T μ β ˜ + D u μ ) T a j a μ μ α ˜ a + ( T β ˜ 1 ) μ q μ + a ( μ a T α ˜ a ) μ j a μ ,
where we employed the relations (3) and (5) and eliminated the terms D ϵ and D n a using Equations (54) and (55). Requiring R 0 , we identify from Equation (60) β ˜ = β , α ˜ a = α a . In the third and the fourth terms, we can replace μ μ due to the orthogonality conditions (33). Then, we have
R = Π θ + π μ ν σ μ ν + q μ ( T μ β + D u μ ) T a j a μ μ α a .
We can further simplify the third term by approximating D u μ = h 1 μ p from Equation (12) since q μ is already of the order O 1 . Using Equation (18), we obtain
T μ β + D u μ = T a n a h μ α a .
Substituting this result into Equation (61) and recalling the definition (44), we obtain finally
R = Π θ + π μ ν σ μ ν T a J a μ μ α a .
The second law of thermodynamics implies R 0 , which requires the right-hand side of Equation (63) to be a quadratic form of thermodynamic forces θ , π μ ν and μ α a . Assuming linear dependence of the dissipative fluxes on the thermodynamic forces, we obtain the constitutive relations
Π = ζ θ , π μ ν = 2 η σ μ ν , J a μ = b χ a b μ α b ,
where η and ζ are called the shear and the bulk viscosities, respectively, and χ a b is the matrix of the diffusion coefficients. These relations together with Equations (54)–(56) constitute a closed system of equations, which are known as relativistic NS equations.
From Equations (63) and (64), we obtain for the dissipation function
R = ζ θ 2 + 2 η σ μ ν σ μ ν T a b χ a b μ α a μ α b = 1 ζ Π 2 + 1 2 η π μ ν π μ ν T a b ( χ 1 ) a b J a μ J b μ .
The positivity of this expression is guaranteed by the positivity of the viscosity coefficients and the eigenvalues of the matrix χ a b (recall that the diffusion fluxes J a μ are spatial). We see from Equation (65) that the contribution of diffusion processes to the function R depends only on the relative diffusion currents J a μ , but not on the currents j a μ and q μ separately. This result, which was obtained in the framework of the relativistic NS (first-order) theory, is the direct consequence of the Lorentz invariance and indicates that the dissipation in the fluid is independent of the choice of the fluid velocity field, as expected.
The entropy current given by Equations (58) and (59) can now be written as
S 1 μ = s u μ + β q μ a α a j a μ = s h P μ a α a J a μ ,
where we used Equations (5), (16), (32) and (44) to obtain the second relation. One can give a simple interpretation of the second expression in Equation (66). Recalling that P μ / h = u L μ is the fluid velocity measured in the L-frame (see Equation (32)), we observe that the first term on the right-hand side of Equation (66) is the entropy current that is convected together with the energy. The second term arises as a result of the relative flow between the energy and the charges and, therefore, should be identified with the irreversible part of the entropy flow. Using Equations (5), (16), (32) and (44), we can rewrite Equation (66) also in the following form:
S 1 μ = β P μ a α a N a μ = p β μ + β ν T μ ν a α a N a μ ,
which formally coincides with its counterpart of ideal hydrodynamics (15).
If we have only one sort of conserved charge ( l = 1 ), then instead of the third relation in Equation (64), we have simply
J μ = χ μ α , h μ = κ n T 2 h μ α ,
where we used the definition of the heat current given by Equation (51) and introduced the coefficient of thermal conductivity via
κ = h n T 2 χ .
Equation (69) establishes the relation between the thermal conductivity and the diffusion coefficient χ . Thus, there are only three independent transport coefficients in the first-order theory which relate the irreversible fluxes to the corresponding thermodynamic forces. Employing Equation (62), we can write the heat flux (68) in the following way
h μ = κ μ T T h μ p ,
which in the fluid rest frame reads
h i = κ i T T h i p .
Equation (71) is the relativistic generalization of the well-known Fourier law h i = κ i T of the non-relativistic hydrodynamics [14].
The expressions (65)–(67) in the case of l = 1 reduce to
R = ζ θ 2 + 2 η σ μ ν σ μ ν T χ μ α 2 = 1 ζ Π 2 + 1 2 η π μ ν π μ ν 1 κ T h μ h μ ,
and
S 1 μ = s u μ + β q μ α j μ = s n N μ + β h μ .
The last relation in Equation (73) is the decomposition of the entropy flow into its reversible and irreversible components observed from the E-frame.
Note that in the case where there are no conserved charges, i.e., when l = 0 , the heat conduction and/or diffusion phenomena are absent [69].

5. Israel–Stewart (Second-Order) Theory

The first-order theory described in the previous section turns out to be acausal, and, therefore, cannot be regarded as a consistent theory of relativistic dissipative fluids. The origin of acausality lies in the constitutive relations (64), which imply that the thermodynamic forces generate dissipative fluxes instantaneously [15,17,22,23]. In addition, this theory suffers also from instability, which is a consequence of acausality [15,16,17,18,25]. However, as pointed out in the introduction, acausalities, and instabilities are a consequence of the matching procedure to the local-equilibrium reference state which can be generalized to obtain causal and stable first-order dissipative hydrodynamics [19,20,21].
It turns out that to recover the causality, the entropy current S μ is required to be at least a quadratic function of the dissipative fluxes. This idea of an extension of the entropy current up to the second order was first proposed by Müller [22] for nonrelativistic fluids. For relativistic fluids, a similar second-order theory was developed by Israel and Stewart [23,24]. In this subsection, we review briefly the Israel–Stewart (IS) formulation of causal hydrodynamics, following mainly Ref. [23].
As a starting point, the entropy current given by Equations (66) and (67) is extended up to the second order in dissipative quantities j a μ , q μ , π μ ν and Π . It is worth stressing that, despite the frame dependence of these quantities, the entropy current S μ along with the energy–momentum tensor T μ ν and the charge currents N a μ should be regarded as a primary variable and should be therefore frame-independent up to the second order [23,24] (we note that the expressions (66) and (67) are frame-independent only at the first order in deviations from equilibrium).
We write now the entropy current in the following form:
S μ = S 1 μ R ¯ μ R μ ,
where S 1 μ is the first-order contribution given by Equations (66) and (67), and the terms R ¯ μ and R μ collect all possible second-order corrections.
The most general form of the vectors R μ and R ¯ μ in a generic frame is [23]
R μ = β u μ 2 β Π Π 2 + β π π α β π α β a b β J a b J a α J b α
+ β Π a α Π a J a μ + β π ν μ a α π a J a ν ,
R ¯ μ = β h 1 Π q μ π μ ν q ν 1 2 u μ q ν q ν ,
where the new coefficients β Π , β π , β J a b = β J b a , α Π a and α π a are unknown functions of ϵ and n a . The vector R μ is frame-independent up to the second order, and R ¯ μ collects the second-order contributions to S μ which are not frame-independent to the order O 2 (note that we use different metric signature from Ref. [23]). Because the contribution S 1 μ is frame-independent only to the first order, the term R ¯ μ is necessary to provide the frame-independence of the total entropy current (74) [23].
The terms in Equations (75) and (76) which are proportional to u μ are responsible for the second-order corrections to the equilibrium entropy density s. Indeed, the non-equilibrium entropy density is identified with S = S μ u μ , which can be found from Equations (74)–(76)
S = s β 2 β Π Π 2 + β π π α β π α β a b β J a b J a α J b α h 1 q α q α s ,
where we used Equations (66). The last inequality in Equation (77) requires the entropy density in non-equilibrium states to be smaller than its equilibrium value s. From here, we conclude that β Π 0 , β π 0 , and the matrix β J a b is positive-semidefinite. The term in Equation (77), which is proportional to q μ q μ represents the shift in the entropy density because of the change of the reference frame and is automatically negative. Those terms in Equations (75) and (76) which are orthogonal to u μ represent the irreversible entropy flux arising from couplings between the diffusion and viscous fluxes.
The phenomenological equations for the dissipative fluxes should be found again from the positivity condition of the dissipative function
R = T μ S μ = T μ S 1 μ T μ R ¯ μ T μ R μ 0 .
The first term in Equation (78) was already computed in Equation (61):
T μ S 1 μ = Π θ + π μ ν σ μ ν + q μ ( T μ β + D u μ ) T a j a μ μ α a .
Now we use Equation (56) to eliminate the acceleration term in Equation (79):
h D u μ = μ p + μ Π Π D u μ Δ μ ν D q ν q ν ν u μ q μ θ Δ μ ν α π α ν = h T μ β + T a n a μ α a + μ Π Π D u μ Δ μ ν D q ν q ν ν u μ q μ θ Δ μ ν α π α ν ,
where we used the second relation of Equation (18) to modify the pressure gradient in Equation (80). Equation (80) differs from Equation (62) by additional second-order terms, which now cannot be neglected. Substituting Equation (80) in the third term of Equation (79) and recalling the definitions (44), we obtain
T μ S 1 μ = Π θ + π μ ν σ μ ν T a J a μ μ α a + h 1 ( q μ μ Π Π q μ D u μ q μ D q μ q μ q ν ν u μ q μ q μ θ q ν μ π μ ν ) ,
where we used the notation μ = Δ μ ν ν . Using Equations (76) and (81), we obtain
T μ S 1 μ T μ R ¯ μ = Π θ + h 1 μ q μ + γ ¯ Π h 1 q μ D u μ + γ ˜ Π T q μ μ ( β h 1 ) + π μ ν σ μ ν + h 1 μ q ν + γ ¯ π T q ν μ ( β h 1 ) T a J a μ μ α a q μ [ h 1 q ν ν u μ + 1 2 h 1 q μ θ 1 2 T q μ D ( β h 1 ) + ( 1 γ ¯ Π ) h 1 Π D u μ + ( 1 γ ˜ Π ) T Π μ ( β h 1 ) ( 1 γ ¯ π ) T π μ ν ν ( β h 1 ) ] .
Here, we introduced three additional coefficients γ ¯ Π , γ ˜ Π and γ ¯ π because there is an ambiguity in “sharing” the terms involving Π q μ and π ν μ q ν between the diffusion and viscous fluxes [15]. Note that Ref. [23] assumed γ ¯ π = γ ¯ Π = γ ˜ Π = 0 .
Now we take the divergence of Equation (75):
μ R μ = β Π β Π D Π + 1 2 T Π μ ( β β Π u μ ) + a α Π a μ J a μ + γ Π T a J a μ μ ( β α Π a ) + β π μ ν β π D π μ ν + 1 2 T π μ ν λ ( β β π u λ ) + a α π a μ J a ν + γ π T a J a ν μ ( β α π a ) β a J a μ [ b β J a b D J b μ + 1 2 T b J b μ ν ( β β J a b u ν ) α Π a μ Π α π a ν π μ ν ( 1 γ Π ) T Π μ ( β α Π a ) ( 1 γ π ) T π μ ν ν ( β α π a ) ] ,
where we introduced again two additional coefficients γ Π and γ π to “share” the terms containing Π J a μ and π ν μ J a ν between the diffusion and viscous fluxes [15]. We kept also all terms that contain gradients of transport coefficients which were neglected in Ref. [23].
Combining Equations (82) and (83), we obtain for the dissipative function (78)
R = Π [ θ + β Π Π ˙ + a α Π a μ J a μ + γ Π T a J a μ μ ( β α Π a ) + 1 2 T Π μ ( β β Π u μ ) + h 1 μ q μ + γ ¯ Π h 1 q μ u ˙ μ + γ ˜ Π T q μ μ ( β h 1 ) ] + π μ ν [ σ μ ν β π π ˙ μ ν a α π a < μ J a ν > γ π T a J a < ν μ > ( β α π a ) 1 2 T π μ ν λ ( β β π u λ ) + h 1 < μ q ν > + γ ¯ π T q < ν μ > ( β h 1 ) ] a J a μ [ T μ α a b β J a b J ˙ b μ 1 2 T b J b μ ν ( β β J a b u ν ) + α Π a μ Π + α π a Δ α μ ν π ν α + ( 1 γ Π ) T Π μ ( β α Π a ) + ( 1 γ π ) T π μ ν ν ( β α π a ) ] q μ [ h 1 q ν ν u μ + 1 2 h 1 q μ θ 1 2 T q μ D ( β h 1 ) + ( 1 γ ¯ Π ) h 1 Π u ˙ μ + ( 1 γ ˜ Π ) T Π μ ( β h 1 ) ( 1 γ ¯ π ) T π μ ν ν ( β h 1 ) ] ,
where we introduced the short-hand notations
u ˙ μ = D u μ , Π ˙ = D Π , π ˙ μ ν = Δ μ ν ρ σ D π ρ σ , J ˙ a μ = Δ μ ν D J a ν ,
A < μ ν > = Δ μ ν α β A α β ,
and used the symmetries of the dissipative fluxes.
The expressions in the square brackets in Equation (84) are called generalized or extended thermodynamic forces. Requiring R 0 , fixing for simplicity the L-frame q μ = 0 and assuming linear relations between these forces and the dissipative fluxes, we obtain the following evolution equations:
Π = ζ [ θ + β Π Π ˙ + a α Π a μ J a μ + γ Π T a J a μ μ ( β α Π a ) + 1 2 T Π μ ( β β Π u μ ) ] ,
π μ ν = 2 η [ σ μ ν β π π ˙ μ ν a α π a < μ J a ν > γ π T a J a < ν μ > ( β α π a ) 1 2 T π μ ν λ ( β β π u λ ) ] ,
J a μ = b χ a b [ μ α b β c β J b c J ˙ c μ 1 2 c J c μ ν ( β β J b c u ν ) + β α Π b μ Π + β α π b Δ μ α ν π ν α + ( 1 γ Π ) Π μ ( β α Π b ) + ( 1 γ π ) π μ ν ν ( β α π b ) ] .
Equations (87)–(89) in this form for single type of conserved charges were derived in Ref. [15] and were written in the E-frame. In the original papers of Israel and Stewart [23,24,70], the terms which are nonlinear in thermodynamic forces and dissipative fluxes were dropped, although the general case of chemically reacting multicomponent mixtures was discussed in Ref. [23]. The current form generalizes the hydrodynamics equations derived in Refs. [15,23,24,70] to the case of multiple independent flavors of conserved charges, where all second-order terms arising from the entropy current are kept.
With the aid of Equations (87)–(89), the dissipative function (84) obtains the form
R = 1 ζ Π 2 + 1 2 η π μ ν π μ ν T a b J a μ ( χ 1 ) a b J b μ ,
which formally coincides with Equation (65). Defining relaxation times according to
τ Π = ζ β Π , τ π = 2 η β π , τ J a c = β b χ a b β J b c ,
we can write Equations (87)–(89) in the following form:
τ Π Π ˙ + Π = ζ θ ζ [ a α Π a μ J a μ + γ Π T a J a μ μ ( β α Π a ) + 1 2 T Π μ ( β β Π u μ ) ] ,
τ π π ˙ μ ν + π μ ν = 2 η σ μ ν 2 η [ a α π a < μ J a ν > + γ π T a J a < ν μ > ( β α π a ) + 1 2 T π μ ν λ ( β β π u λ ) ] ,
c τ J a c J ˙ c μ + J a μ = b χ a b μ α b + b χ a b [ 1 2 c J c μ ν ( β β J b c u ν ) + β α Π b μ Π + β α π b Δ μ α ν π ν α + ( 1 γ Π ) Π μ ( β α Π b ) + ( 1 γ π ) π μ ν ν ( β α π b ) ] .
The first terms on the right-hand sides of these equations represent the corresponding NS contributions to the dissipative fluxes; see Equation (64). The first terms on the left-hand sides incorporate the relaxation of the dissipative fluxes to their NS values on finite time scales given by Equation (91). Thus, these relaxation terms imply a delay in the response of the dissipative fluxes to thermodynamic forces and recover the causality of the theory [18,25]. The rest of the terms in Equations (92)–(94) are responsible for spatial inhomogeneities in the dissipative fluxes as well as nonlinear couplings between different dissipative processes.
We note that the derivation of the second-order hydrodynamics from the kinetic theory produces additional terms which are not obtained within the phenomenological theory [24]. The derivation of complete IS equations from kinetic theory is discussed in Refs. [30,38,71,72,73].
In the case of one conserved current, we have instead of Equations (92)–(94)
τ Π Π ˙ + Π = ζ θ ζ α Π μ J μ + γ Π T J μ μ ( β α Π ) + 1 2 T Π μ ( β β Π u μ ) ,
τ π π ˙ μ ν + π μ ν = 2 η σ μ ν 2 η [ α π < μ J ν > + γ π T J < ν μ > ( β α π ) + 1 2 T π μ ν λ ( β β π u λ ) ] ,
τ J J ˙ μ + J μ = χ μ α + χ [ β α Π μ Π + β α π Δ μ α ν π ν α 1 2 J μ ν ( β β J u ν ) + ( 1 γ Π ) Π μ ( β α Π ) + ( 1 γ π ) π μ ν ν ( β α π ) ] ,
where τ J = β χ β J .

6. Summary

We provided a review of the phenomenological theory of second-order relativistic hydrodynamics for systems with multiple conserved charges with an extension to multiflavor fluids. The hydrodynamic state of the system is described through the energy–momentum tensor and the 4-currents of conserved charges. We reviewed the derivation and content of the equations at zeroth, first, and second order in gradient expansion of the energy–momentum tensor and currents, which led us to the ideal Navier–Stokes and Israel–Stewart hydrodynamics, respectively. From the positivity condition of the dissipative function, the most general set of dissipative processes was identified which contains also the relative diffusions between different conserved currents. We kept all the second-order gradient terms arising from the second law of thermodynamics, as well as the nonlinear terms in the thermodynamic gradients and dissipative fluxes, which were omitted in the canonical Israel–Stewart theory, but were kept in Ref. [15] in the case of one conserved flavor.
The hydrodynamics theory exposed in this article is phenomenological in nature. More formal but also complex derivations are available in the literature which utilize different concepts and approaches of statistical mechanics, for example, quasiparticle Boltzmann equation [38,39,40] or non-equilibrium statistical operator [41,43]. Nevertheless, the phenomenological theory can be reasonably expected to remain in the arsenal of theoretical tools that can be deployed in studies of fluid in new settings and/or under various external fields.
A few closing remarks are in order concerning the numerical implementations of the second-order dissipative hydrodynamics. First of all, the choice between the frames (i.e., Eckart vs. Landau) may be significantly influenced by the required computational cost, which in turn depends on such factors as system size, boundary conditions and the level of accuracy required. In general, the Landau frame is known to be computationally more expensive compared to the Eckart frame, but it also provides more insight into the behavior of the fluid. Eventually, the choice between the Landau and Eckart frames will depend on the specific problem at hand and the trade-off between accuracy and computational efficiency. A number of methods are available for the solution of the generalized (second-order) Navier–Stokes equations for multifluids, for example, finite difference method, finite volume method, spectral methods, and Lagrangian particle methods; see Ref. [9] Section 5 for a pedagogical discussion in the context of relativistic heavy ion collisions. Again, the choice of the method will be dictated by the specifics of the problem at hand, the computational cost, and required accuracy.

Author Contributions

A.H. and A.S. equally contributed to this research. All authors have read and agreed to the published version of the manuscript.

Funding

A.H. and A.S. were funded by Volkswagen Foundation (Hannover, Germany) grant No. 96 839. A.S. was funded by Deutsche Forschungsgemeinschaft Grant No. SE 1836/5-2 and the Polish NCN Grant No. 2020/37/B/ST9/01937 at Wrocław University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaiswal, A.; Roy, V. Relativistic hydrodynamics in heavy-ion collisions: General aspects and recent developments. Adv. High Energy Phys. 2016, 2016, 9623034. [Google Scholar] [CrossRef]
  2. Florkowski, W.; Heller, M.P.; Spaliński, M. New theories of relativistic hydrodynamics in the LHC era. Rep. Prog. Phys. 2018, 81, 046001. [Google Scholar] [CrossRef] [PubMed]
  3. Denicol, G.S.; Rischke, D.H. Microscopic Foundations of Relativistic Fluid Dynamics; Springer: Berlin/Heidelberg, Germany, 2021. [Google Scholar]
  4. Elfner, H.; Müller, B. The exploration of hot and dense nuclear matter: Introduction to relativistic heavy-ion physics. arXiv 2022, arXiv:2210.12056. [Google Scholar]
  5. Font, J.A. Numerical hydrodynamics in general relativity. Living Rev. Rel. 2000, 3, 2. [Google Scholar] [CrossRef] [PubMed]
  6. Rezzolla, L.; Zanotti, O. Relativistic Hydrodynamics; Oxford University Press: Oxford, UK, 2013. [Google Scholar]
  7. Baiotti, L.; Rezzolla, L. Binary neutron star mergers: A review of Einstein’s richest laboratory. Rep. Prog. Phys. 2017, 80, 096901. [Google Scholar] [CrossRef]
  8. Busza, W.; Rajagopal, K.; van der Schee, W. Heavy Ion Collisions: The Big Picture and the Big Questions. Annu. Rev. Nucl. Part. Sci. 2018, 68, 339–376. [Google Scholar] [CrossRef]
  9. Jeon, S.; Heinz, U. Introduction to hydrodynamics. Int. J. Mod. Phys. 2015, 24, 1530010. [Google Scholar] [CrossRef]
  10. Baiotti, L. Gravitational waves from neutron star mergers and their relation to the nuclear equation of state. Prog. Part. Nucl. Phys. 2019, 109, 103714. [Google Scholar] [CrossRef]
  11. Haskell, B.; Sedrakian, A. Superfluidity and Superconductivity in Neutron Stars. Astrophys. Space Sci. Libr. 2018, 457, 401–454. [Google Scholar]
  12. Andersson, N.; Comer, G.L. Relativistic fluid dynamics: Physics for many different scales. Living Rev. Rel. 2021, 24, 3. [Google Scholar] [CrossRef]
  13. Eckart, C. The Thermodynamics of Irreversible Processes. III. Relativistic theory of the simple fluid. Phys. Rev. 1940, 58, 919–924. [Google Scholar] [CrossRef]
  14. Landau, L.; Lifshitz, E. Fluid Mechanics. In Fluid Mechanics; Landau, L., Lifshitz, E., Eds.; Butterworth-Heinemann: Oxford, UK, 1987. [Google Scholar]
  15. Hiscock, W.A.; Lindblom, L. Stability and causality in dissipative relativistic fluids. Ann. Phys. 1983, 151, 466–496. [Google Scholar] [CrossRef]
  16. Hiscock, W.A.; Lindblom, L. Generic instabilities in first-order dissipative relativistic fluid theories. Phys. Rev. D 1985, 31, 725–733. [Google Scholar] [CrossRef]
  17. Hiscock, W.A.; Lindblom, L. Linear plane waves in dissipative relativistic fluids. Phys. Rev. D 1987, 35, 3723. [Google Scholar] [CrossRef]
  18. Denicol, G.S.; Kodama, T.; Koide, T.; Mota, P. Stability and causality in relativistic dissipative hydrodynamics. J. Phys. G Nucl. Phys. 2008, 35, 115102. [Google Scholar] [CrossRef]
  19. Bemfica, F.S.; Disconzi, M.M.; Noronha, J. Nonlinear Causality of General First-Order Relativistic Viscous Hydrodynamics. Phys. Rev. D 2019, 100, 104020, Erratum in Phys. Rev. D 2022, 105, 69902. [Google Scholar] [CrossRef]
  20. Kovtun, P. First-order relativistic hydrodynamics is stable. J. High Energy Phys. 2019, 10, 34. [Google Scholar] [CrossRef]
  21. Noronha, J.; Spaliński, M.; Speranza, E. Transient Relativistic Fluid Dynamics in a General Hydrodynamic Frame. Phys. Rev. Lett. 2022, 128, 252302. [Google Scholar] [CrossRef] [PubMed]
  22. Müller, I. Zum Paradoxon der Wärmeleitungstheorie. Z. Phys. 1967, 198, 329–344. [Google Scholar] [CrossRef]
  23. Israel, W. Nonstationary irreversible thermodynamics: A causal relativistic theory. Ann. Phys. 1976, 100, 310–331. [Google Scholar] [CrossRef]
  24. Israel, W.; Stewart, J.M. Transient relativistic thermodynamics and kinetic theory. Ann. Phys. 1979, 118, 341–372. [Google Scholar] [CrossRef]
  25. Pu, S.; Koide, T.; Rischke, D.H. Does stability of relativistic dissipative fluid dynamics imply causality? Phys. Rev. D 2010, 81, 114039. [Google Scholar] [CrossRef] [Green Version]
  26. Baier, R.; Romatschke, P.; Thanh Son, D.; Starinets, A.O.; Stephanov, M.A. Relativistic viscous hydrodynamics, conformal invariance, and holography. J. High Energy Phys. 2008, 4, 100. [Google Scholar] [CrossRef]
  27. Betz, B.; Henkel, D.; Rischke, D.H. Complete second-order dissipative fluid dynamics. J. Phys. G Nucl. Phys. 2009, 36, 064029. [Google Scholar] [CrossRef]
  28. Romatschke, P. Relativistic viscous fluid dynamics and non-equilibrium entropy. Class. Quantum Gravity 2010, 27, 025006. [Google Scholar] [CrossRef]
  29. Tsumura, K.; Kunihiro, T. Second-order relativistic hydrodynamic equations for viscous systems; How does the dissipation affect the internal energy? Phys. Lett. B 2010, 690, 255–260. [Google Scholar] [CrossRef]
  30. Betz, B.; Denicol, G.S.; Koide, T.; Molnár, E.; Niemi, H.; Rischke, D.H. Second order dissipative fluid dynamics from kinetic theory. In Proceedings of the European Physical Journal Web of Conferences, Budapest, Hungary, 15–20 August 2011; Volume 13, p. 07005. [Google Scholar]
  31. Moore, G.D.; Sohrabi, K.A. Kubo Formulas for Second-Order Hydrodynamic Coefficients. Phys. Rev. Lett. 2011, 106, 122302. [Google Scholar] [CrossRef]
  32. Moore, G.D.; Sohrabi, K.A. Thermodynamical second-order hydrodynamic coefficients. J. High Energy Phys. 2012, 11, 148. [Google Scholar] [CrossRef]
  33. Jaiswal, A.; Bhalerao, R.S.; Pal, S. Complete relativistic second-order dissipative hydrodynamics from the entropy principle. Phys. Rev. C 2013, 87, 021901. [Google Scholar] [CrossRef]
  34. Jaiswal, A.; Friman, B.; Redlich, K. Relativistic second-order dissipative hydrodynamics at finite chemical potential. Phys. Lett. B 2015, 751, 548–552. [Google Scholar] [CrossRef]
  35. Florkowski, W.; Jaiswal, A.; Maksymiuk, E.; Ryblewski, R.; Strickland, M. Relativistic quantum transport coefficients for second-order viscous hydrodynamics. Phys. Rev. C 2015, 91, 054907. [Google Scholar] [CrossRef]
  36. Finazzo, S.I.; Rougemont, R.; Marrochio, H.; Noronha, J. Hydrodynamic transport coefficients for the non-conformal quark-gluon plasma from holography. J. High Energy Phys. 2015, 2, 51. [Google Scholar] [CrossRef]
  37. Tinti, L.; Jaiswal, A.; Ryblewski, R. Quasiparticle second-order viscous hydrodynamics from kinetic theory. Phys. Rev. D 2017, 95, 054007. [Google Scholar] [CrossRef]
  38. Denicol, G.S.; Niemi, H.; Molnár, E.; Rischke, D.H. Derivation of transient relativistic fluid dynamics from the Boltzmann equation. Phys. Rev. D 2012, 85, 114047. [Google Scholar] [CrossRef]
  39. Denicol, G.S.; Huang, X.G.; Koide, T.; Rischke, D.H. Consistency of field-theoretical and kinetic calculations of viscous transport coefficients for a relativistic fluid. Phys. Lett. B 2012, 708, 174–178. [Google Scholar] [CrossRef]
  40. Molnár, E.; Niemi, H.; Rischke, D.H. Derivation of anisotropic dissipative fluid dynamics from the Boltzmann equation. Phys. Rev. D 2016, 93, 114025. [Google Scholar] [CrossRef]
  41. Harutyunyan, A.; Sedrakian, A.; Rischke, D.H. Relativistic Dissipative Fluid Dynamics from the Non-Equilibrium Statistical Operator. Particles 2018, 1, 155–165. [Google Scholar] [CrossRef]
  42. Becattini, F.; Buzzegoli, M.; Grossi, E. Reworking Zubarev’s Approach to Nonequilibrium Quantum Statistical Mechanics. Particles 2019, 2, 197–207. [Google Scholar] [CrossRef]
  43. Harutyunyan, A.; Sedrakian, A.; Rischke, D.H. Relativistic second-order dissipative hydrodynamics from Zubarev’s non-equilibrium statistical operator. Ann. Phys. 2022, 438, 168755. [Google Scholar] [CrossRef]
  44. Tokarchuk, M.; Hlushak, P. Unification of Thermo Field Kinetic and Hydrodynamics Approaches in the Theory of Dense Quantum–Field Systems. Particles 2018, 2, 1–13. [Google Scholar] [CrossRef]
  45. Mirón Granese, N.; Kandus, A.; Calzetta, E. Field Theory Approaches to Relativistic Hydrodynamics. Entropy 2022, 24, 1790. [Google Scholar] [CrossRef] [PubMed]
  46. Buzzegoli, M.; Grossi, E.; Becattini, F. General equilibrium second-order hydrodynamic coefficients for free quantum fields. J. High Energy Phys. 2017, 2017, 91. [Google Scholar] [CrossRef] [Green Version]
  47. Becattini, F.; Buzzegoli, M.; Palermo, A. Exact equilibrium distributions in statistical quantum field theory with rotation and acceleration: Scalar field. J. High Energy Phys. 2021, 2021, 101. [Google Scholar] [CrossRef]
  48. Cao, Z.; Hattori, K.; Hongo, M.; Huang, X.G.; Taya, H. Gyrohydrodynamics: Relativistic spinful fluid with strong vorticity. Prog. Theor. Exp. Phys. 2022, 2022, 071D01. [Google Scholar] [CrossRef]
  49. Hongo, M.; Huang, X.G.; Kaminski, M.; Stephanov, M.; Yee, H.U. Relativistic spin hydrodynamics with torsion and linear response theory for spin relaxation. J. High Energy Phys. 2021, 2021, 150. [Google Scholar] [CrossRef]
  50. Weickgenannt, N.; Wagner, D.; Speranza, E.; Rischke, D.H. Relativistic second-order dissipative spin hydrodynamics from the method of moments. Phys. Rev. D 2022, 106, 096014. [Google Scholar] [CrossRef]
  51. Gorbar, E.V.; Rybalka, D.O.; Shovkovy, I.A. Second-order dissipative hydrodynamics for plasma with chiral asymmetry and vorticity. Phys. Rev. D 2017, 95, 096010. [Google Scholar] [CrossRef]
  52. Hu, J.; Shi, S. Multicomponent second-order dissipative relativistic hydrodynamics with binary reactive collisions. Phys. Rev. D 2022, 106, 014007. [Google Scholar] [CrossRef]
  53. Almaalol, D.; Dore, T.; Noronha-Hostler, J. Stability of multi-component relativistic viscous hydrodynamics from Israel-Stewart and reproducing DNMR from maximizing the entropy. arXiv 2022, arXiv:2209.11210. [Google Scholar]
  54. Alalawi, H.; Alqahtani, M.; Strickland, M. Resummed Relativistic Dissipative Hydrodynamics. Symmetry 2022, 14, 329. [Google Scholar] [CrossRef]
  55. Buzzegoli, M.; Kharzeev, D.E.; Liu, Y.C.; Shi, S.; Voloshin, S.A.; Yee, H.U. Shear-induced anomalous transport and charge asymmetry of triangular flow in heavy-ion collisions. Phys. Rev. C 2022, 106, L051902. [Google Scholar] [CrossRef]
  56. Lahiri, S. Second order causal hydrodynamics in Eckart frame: Using gradient expansion scheme. Class. Quantum Gravity 2020, 37, 075010. [Google Scholar] [CrossRef]
  57. Bemfica, F.S.; Disconzi, M.M.; Noronha, J.; Scherrer, R.J. Cosmological consequences of first-order general-relativistic viscous fluid dynamics. Phys. Rev. D 2023, 107, 023512. [Google Scholar] [CrossRef]
  58. Most, E.R.; Noronha, J. Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation. Phys. Rev. D 2021, 104, 103028. [Google Scholar] [CrossRef]
  59. Landry, M.J.; Liu, H. A systematic formulation of chiral anomalous magnetohydrodynamics. arXiv 2022, arXiv:2212.09757. [Google Scholar]
  60. Dash, A.; Shokri, M.; Rezzolla, L.; Rischke, D.H. Charge diffusion in relativistic resistive second-order dissipative magnetohydrodynamics. arXiv 2022, arXiv:2211.09459. [Google Scholar]
  61. Harutyunyan, A.; Sedrakian, A. Electrical conductivity of a warm neutron star crust in magnetic fields. Phys. Rev. C 2016, 94, 025805. [Google Scholar] [CrossRef]
  62. Harutyunyan, A.; Sedrakian, A. Electrical conductivity tensor of dense plasma in magnetic fields. In Proceedings of the Modern Physics of Compact Stars 2015 (MPCS2015), Yerevan, Armenia, 30 September–3 October 2015; p. 11. [Google Scholar]
  63. Shternin, P.; Ofengeim, D. Transport coefficients of magnetized neutron star cores. Eur. Phys. J. A 2022, 58, 42. [Google Scholar] [CrossRef]
  64. Hattori, K.; Hongo, M.; Huang, X.G. New Developments in Relativistic Magnetohydrodynamics. Symmetry 2022, 14, 1851. [Google Scholar] [CrossRef]
  65. Monnai, A. Landau and Eckart frames for relativistic fluids in nuclear collisions. Phys. Rev. C 2019, 100, 01490. [Google Scholar] [CrossRef]
  66. Bemfica, F.S.; Disconzi, M.M.; Noronha, J. First-Order General-Relativistic Viscous Fluid Dynamics. Phys. Rev. X 2022, 12, 021044. [Google Scholar] [CrossRef]
  67. Weinberg, S. Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity; Wiley: New York, NY, USA, 1972; p. 688. [Google Scholar]
  68. Weinberg, S. Entropy Generation and the Survival of Protogalaxies in an Expanding Universe. ApJ 1971, 168, 175. [Google Scholar] [CrossRef]
  69. Danielewicz, P.; Gyulassy, M. Dissipative phenomena in quark-gluon plasmas. Phys. Rev. D 1985, 31, 53–62. [Google Scholar] [CrossRef] [PubMed]
  70. Israel, W.; Stewart, J.M. Thermodynamics of nonstationary and transient effects in a relativistic gas. Phys. Lett. A 1976, 58, 213–215. [Google Scholar] [CrossRef]
  71. Betz, B.; Henkel, D.; Rischke, D.H. From kinetic theory to dissipative fluid dynamics. Prog. Part. Nucl. Phys. 2009, 62, 556–561. [Google Scholar] [CrossRef] [Green Version]
  72. Denicol, G.S.; Molnár, E.; Niemi, H.; Rischke, D.H. Derivation of fluid dynamics from kinetic theory with the 14-moment approximation. Eur. Phys. J. A 2012, 48, 170. [Google Scholar] [CrossRef]
  73. Denicol, G.S. Kinetic foundations of relativistic dissipative fluid dynamics. J. Phys. G Nucl. Phys. 2014, 41, 124004. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Harutyunyan, A.; Sedrakian, A. Phenomenological Relativistic Second-Order Hydrodynamics for Multiflavor Fluids. Symmetry 2023, 15, 494. https://doi.org/10.3390/sym15020494

AMA Style

Harutyunyan A, Sedrakian A. Phenomenological Relativistic Second-Order Hydrodynamics for Multiflavor Fluids. Symmetry. 2023; 15(2):494. https://doi.org/10.3390/sym15020494

Chicago/Turabian Style

Harutyunyan, Arus, and Armen Sedrakian. 2023. "Phenomenological Relativistic Second-Order Hydrodynamics for Multiflavor Fluids" Symmetry 15, no. 2: 494. https://doi.org/10.3390/sym15020494

APA Style

Harutyunyan, A., & Sedrakian, A. (2023). Phenomenological Relativistic Second-Order Hydrodynamics for Multiflavor Fluids. Symmetry, 15(2), 494. https://doi.org/10.3390/sym15020494

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop