Pith. sign in

REVIEW 1 major objections 5 minor 46 references

A simple family of solutions of relativistic viscous hydrodynamics for fireballs with Hubble flow and ellipsoidal symmetry

T0 review · 1 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read New exact analytic solutions of relativistic viscous hydrodynamics show that for Hubble-type ellipsoidal fireballs the shear viscosity cancels completely, leaving bulk viscosity as the only dissipative influence on the temperature…

desk verdict Solid exact viscous Hubble-flow solutions, but the advertised ellipsoidal symmetry is not realized in any physical solution — the thermodynamics ends up depending on τ alone. read the letter →

arxiv 1909.02498 v3 pith:4Q6ANZ5S submitted 2019-09-05 nucl-th hep-ph

classification nucl-thhep-ph PACS 20.2420.25
keywords relativisticviscoushydrodynamicsHubbleflowbulkviscosityexactanalyticsolutionsellipsoidalsymmetryheavy-ioncollisionsshearcancellationIsrael-Stewart
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper constructs new exact, analytic solutions of relativistic viscous hydrodynamics for expanding fireballs whose velocity is the Hubble-type profile $u^\mu = x^\mu/\tau$ and whose density and temperature profiles have ellipsoidal symmetry. In these solutions the shear-viscosity terms vanish identically, and when the temperature depends only on the proper time $\tau$, the heat-conduction terms cancel as well, so bulk viscosity is the only dissipative mechanism left. The paper reduces the full partial differential equations to ordinary differential equations for the pressure and solves them for six concrete assumptions about how the bulk viscosity coefficient depends on temperature, entropy density, particle density, or bulk pressure. The physical solutions all show that bulk viscosity slows the cooling relative to a perfect fluid, and different scenarios for $\zeta$ produce qualitatively different temperature histories, raising the possibility that measured energy-density and temperature evolution in heavy-ion collisions could constrain the bulk viscosity.

What carries the argument

The central object is the Hubble-type four-velocity profile $u^\mu = x^\mu/\tau$, with proper time $\tau=\sqrt{t^2-r_x^2-r_y^2-r_z^2}$. Its defining identities are zero acceleration, $u^\nu\partial_\nu u^\mu=0$, zero shear, $\sigma^{\mu\nu}=0$, and $u^\mu\partial_\mu S=0$ for the ellipsoidal scaling variable $S=r_x^2/X^2+r_y^2/Y^2+r_z^2/Z^2$. These identities collapse the divergences of the viscous stress tensor into equations that depend only on $\tau$; the machinery of the argument is the projection of $\partial_\mu T^{\mu\nu}=0$ into parts parallel and pseudo-orthogonal to $x^\mu$, reducing a multi-dimensional PDE problem to the ordinary differential equations (36)-(38).

What would settle it

Measure or simulate the velocity field of an expanding heavy-ion fireball and compute the shear tensor $\sigma^{\mu\nu}$: the paper's cancellation of shear viscosity holds only where $\sigma^{\mu\nu}=0$ exactly, so a single spacetime region with non-vanishing shear would break the claimed applicability there.

Watch

Extended reading notes

Core claim

For a relativistic fireball with Hubble-like velocity field $u^\mu=x^\mu/\tau$ and ellipsoidal symmetry, the energy-momentum conservation equations of first-order viscous hydrodynamics reduce to a single ordinary differential equation for the pressure, Eq. (36), in the Navier-Stokes case and to the pair (37)-(38) in the Israel-Stewart case. Because the shear tensor $\sigma^{\mu\nu}$ vanishes identically for this velocity field, all shear-viscosity effects cancel; if the temperature depends only on the proper time $\tau$, all heat-conduction terms cancel as well. The paper solves the reduced equations for six scenarios for the bulk viscosity coefficient, obtaining closed-form pressure and temperature histories. Four of these cases are thermodynamically physical, while two produce unphysical reheating or divergence and are rejected. In the physical cases bulk viscosity always slows the decrease of energy density and temperature relative to the perfect-fluid limit, and the shape of the cooling curve depends on the assumed scaling of $\zeta$.

Load-bearing premise

The entire construction rests on assuming that the expanding matter moves with a Hubble-type velocity profile whose shear and acceleration vanish; if a real fireball's flow deviates from this boost-invariant, acceleration-free form, the cancellation of shear viscosity and heat conduction fails and the solutions do not apply.

Editorial extensions

If this is right

  • Bulk viscosity changes the cooling law of the fireball: constant $\zeta$, $\zeta\propto s$, and $\zeta\propto n$ produce noticeably different temperature histories, so measured cooling curves can in principle distinguish these scenarios.
  • Shear viscosity drops out of these solutions for any $\eta(T)$ or $\eta/s$ function, making the solutions exact benchmarks for numerical viscous-hydro codes that must reproduce the same $\tau$-dependent evolution.
  • Two of the six scenarios, Case B (constant $\zeta$ with conserved particle number) and Case E ($\zeta\propto T^\kappa$ with $p=nT$), lead to unphysical reheating or divergence and are therefore excluded within this flow family.
  • Because thermal conductivity also cancels whenever $T=T(\tau)$, the solutions remain valid for arbitrary $\lambda$ in cases with no conserved charge (or with $V(S)=1$), so they isolate bulk viscosity as the only active dissipative transport coefficient.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If a real heavy-ion system ever approaches Hubble flow, shear viscosity would be hidden from the temperature profile by geometry rather than by smallness of $\eta$, so extracting bulk viscosity would require independent evidence that the flow is truly boost-invariant.
  • The technique of choosing a zero-shear velocity field to cancel shear viscosity could be carried over to other flow geometries, such as boosted or accelerating profiles, to construct analogues that isolate bulk viscosity in less symmetric settings.
  • The unphysical cases suggest a model-discrimination rule: a bulk-viscosity ansatz that grows too rapidly as the conserved density dilutes will produce runaway heating, so the observed absence of reheating in heavy-ion data can set a lower bound on how fast $\zeta$ may grow with $T$ or $n$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 5 minor

Summary. This paper constructs exact analytic solutions of relativistic viscous hydrodynamics for the Hubble-type velocity field u^mu = x^mu/tau. For this profile the shear tensor vanishes identically and, when thermal conductivity is neglected, energy-momentum conservation reduces to the ordinary differential equations (36)-(38). The authors solve these equations for six scenarios for the bulk viscosity (constant, proportional to s, n, T^kappa, or Pi) and for two equations of state (p = nT or epsilon = kappa p), in both the Navier-Stokes and one Israel-Stewart case, and they present explicit p(tau) and T(tau) formulas together with illustrative plots. They find that bulk viscosity slows the cooling of the fireball relative to the perfect-fluid Hubble flow, and they identify Case B and, conditionally, Case E as unphysical.

Significance. The central ODE reduction is sound, and the resulting explicit solutions are useful as benchmarks for numerical viscous hydrodynamics, particularly because the shear viscosity cancels identically for the Hubble profile and the thermal conductivity also cancels under the condition T = T(tau). Exact analytic viscous solutions with bulk viscosity are rare, so this is a genuine contribution. The paper is also honest in flagging the unphysical cases. However, the advertised ellipsoidal symmetry is not realized in the physical solutions, which reduces the significance relative to the title's promise.

major comments (1)
  1. [Title, abstract, and Section VI (Summary); see also Eqs. (21), (25), (29), (36) and Table I] The central claim that the solutions describe fireballs with ellipsoidal symmetry is not supported by the explicit solutions. In the physical Cases A, C, D, and F, all thermodynamic quantities are functions of tau alone: p(tau) and T(tau) are given by Eqs. (41), (43), (46), and (53), with no dependence on the scaling variable S of Eq. (21). The only solution that can retain an arbitrary V(S) in Eq. (25) is Case B, which the authors themselves reject as unphysical because T diverges as tau^(d-1) (Section IV, after Eq. (41), and Section V). In Cases D and E the paper sets V(S)=1 (Section III, after Eq. (25)), and Cases A, C, and F have no conserved density n, so no V(S) exists. The orthogonality condition Eq. (29) forces p - d zeta/tau to be a function of tau alone, and this is what eliminates S-dependence whenever zeta depends on T or n. Consequently, the title and abstract overstate the ellipsoidal-symmetry content of the new solutions; the physical viscous solutions presented are spherically symmetric Hubble-flow solutions with bulk viscosity. Please either construct a genuine ellipsoidal family (for example by allowing S-dependent p in the reduction leading to Eq. (36)) or revise the title, abstract, and Section VI to describe the solutions as Hubble-flow, spherically symmetric solutions.
minor comments (5)
  1. [Eq. (54)] Equation (54a) contains a likely typo: the right-hand side contains Gamma(B, tau/tau_Pi), which makes p_A tau-dependent and inconsistent with Eq. (53). From Eq. (53) and the boundary condition p(tau_0)=p_0, the argument should be Gamma(B, tau_0/tau_Pi).
  2. [Eq. (41)] The exponent of (tau_0/tau) in Eq. (41) is typeset ambiguously as "(d kappa+1)/kappa"; it should be d(kappa+1)/kappa to match the homogeneous solution of Eq. (36). Please clarify.
  3. [End of Section IV] The sentence stating that the solutions are valid for any type of thermal conductivity is too broad without the qualification given two paragraphs earlier: the cancellation requires T = T(tau), which in Cases B, D, and E forces V(S)=1. This condition should be stated in the summary sentence as well.
  4. [Figures 1-4 and Table I] The captions of Figures 1-4 do not state that the plotted solutions have no nontrivial S-dependence (V(S)=1 or no conserved density). Adding this information would make the ellipsoidal-symmetry limitation visible to the reader.
  5. [Section VI] There is a spelling error: "reseach directions" should be "research directions".

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the viscous solutions are derived explicitly from conservation laws and constitutive relations, with no fitted parameter renamed as a prediction.

full rationale

The paper carries out a direct analytic derivation: it postulates the Hubble-type velocity field (Eq. 23), writes the Eckart-frame energy-momentum tensor, imposes ε=κp, projects the conservation equations into Eqs. (28)-(31), and then solves the resulting ODEs (36)-(38) case by case. The integration constants p0, T0, τ0, ζ0, and Π0 are initial-condition scales, not parameters fitted to data, so the temperature and pressure curves in Figs. 1-4 are exact consequences of the stated assumptions rather than fits. The statement that shear viscosity cancels is also not a circular prediction: it follows immediately from σμν=0 for the chosen uμ=xμ/τ, and the paper explicitly presents this as a property of the ansatz. Citations to earlier Hubble-flow perfect-fluid solutions (Refs. 18, 19, 41, 42, 45) are used for comparison and for the standard continuity-equation solution n∝τ^{-d}V(S); they are not invoked to justify the new viscous results. The title's 'ellipsoidal symmetry' claim is arguably stronger than what the explicit solutions deliver, since in the physical cases p and T depend only on τ, but that is an overstatement of scope, not a circular reduction of the derivation to its inputs.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No new particles or forces are introduced. The constants appearing in the solutions are integration constants that set the initial conditions; they are not fitted to data. The model parameters kappa and zeta0 are chosen for illustration but the existence of the solutions does not depend on their specific values.

assumptions (4)
  • domain assumption Eckart frame with Navier-Stokes constitutive relations (Eqs. (4)-(7))
    The bulk viscosity is assumed to be given by Pi = -zeta partial_rho u^rho and the heat flux by a Fourier-like law; this is a first-order theory. The central equations (28)-(31) rely on this.
  • domain assumption Constant equation of state epsilon = kappa p (Eq. (10))
    Used throughout; a constant speed of sound is a simplification not valid near phase transitions.
  • domain assumption Hubble-type velocity profile u^mu = x^mu/tau (Eq. (23))
    This ansatz is the basis of all solutions; it makes the shear tensor vanish identically.
  • domain assumption Israel-Stewart relaxation equations for Case F (Eqs. (8)-(9))
    Second-order theory used only in Case F; requires a relaxation time tauPi and an independent bulk pressure scale Pi0.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A simple family of solutions of relativistic viscous hydrodynamics for fireballs with Hubble flow and ellipsoidal symmetry." pith.science (2026). https://pith.science/paper/4Q6ANZ5S

@misc{pith2026190902498,
  author       = {Pith},
  title        = {Pith review of: A simple family of solutions of relativistic viscous hydrodynamics for fireballs with Hubble flow and ellipsoidal symmetry},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4Q6ANZ5S}},
  note         = {Machine review of arXiv:1909.02498}
}
read the original abstract

New, analytic solutions of relativistic viscous hydrodynamics are presented, describing expanding fireballs with Hubble-like velocity profile and ellipsoidal symmetry, similar to fireballs created in heavy ion collisions. We find that with these specifications, one obtains solutions where the shear viscosity essentially does not influence the time evolution of the system, thus these solutions are particularly adept tools to study the effect of bulk viscosity alone, which always results in a slower decrease of energy density as well as temperature compared to the case of perfect fluid. We investigate different scenarios for the bulk viscosity and find qualitatively different effects on the time evolution which suggests that there is a possibility to infer the value of bulk viscosity from energy density and temperature measurements in high-energy heavy-ion collisions.

Figures

Figures reproduced from arXiv: 1909.02498 by the authors.

Figure 1
Figure 1. FIG. 1: (Color online) Temperature evolution in Case A: no [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: (Color online) Temperature evolution in Case C: no [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: (Color online) Temperature evolution in Case D: non [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4: (Color online) Temperature evolution in Case F: no [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

46 extracted references · 21 canonical work pages

  1. [1]

    S. A. Bass, M. Gyulassy, H. St¨ ocker, W. Greiner, J. Phys. G25, R1 (1999)

  2. [2]

    Gyulassy, and L

    M. Gyulassy, and L. McLerran, Nucl. Phys. A 750, 30-63 (2004)

  3. [3]

    Heinz, and R

    U. Heinz, and R. Snellings, Ann. Rev. Nucl. Part. Sci. 63, 123-151 (2013)

  4. [4]

    Song Huichao, S. A. Bass, U. Heinz, T. Hirano, and Shen Chun, Phys. Rev. Lett. 106, 192301 (2011)

  5. [5]

    Csan´ ad, T

    M. Csan´ ad, T. Cs¨ org˝ o, Ze Fang Jiang and C. B. Yang, Universe. 3 1, 9 (2017) [arXiv:1609.07176]

  6. [6]

    Ze Fang Jiang, C. B. Yang, M. Csan´ ad and T. Cs¨ org˝ o, Phys. Rev. C 97 (2018) no.6, 064906 [arXiv:1711.10740]

  7. [7]

    L. G. Pang, H. Petersen and X. N. Wang, Phys. Rev. C 97 (2018) no.6, 064918 [arXiv:1802.04449]

  8. [8]

    Xiang-Yu Wu and Long-Gang Pang and Guang-You Qin and Xin-Nian Wang, Phys. Rev. C 98, no. 2, 024913 (2018) [arXiv:1805.03762]

Show all 46 references
  1. [9]

    Csan´ ad, G

    Ze Fang Jiang, M. Csan´ ad, G. Kasza,C. B. Yang, and T. Cs¨ org˝ o, Acta Phys. Pol. B Proc. Suppl. vol. 12 (2), pp. 261–267 (2019) [arXiv:1806.05750]

  2. [10]

    W. Chen, S. Cao, T. Luo, L. G. Pang and X. N. Wang, Phys. Lett. B 777 (2018) 86 [arXiv:1704.03648]

  3. [11]

    L. D. Landau, Izv. Akad. Nauk Ser. Fiz. 17, 51 (1953)

  4. [12]

    R. C. Hwa, Phys. Rev. D 10, 2260 (1974)

  5. [13]

    J. D. Bjorken, Phys. Rev. D 27, 140 (1983)

  6. [14]

    S. S. Gubser, Phys. Rev. D82, 085027 (2010) [arXiv:1006.0006]

  7. [15]

    S. S. Gubser and A. Yarom, Nucl. Phys. B846, 469 (2011) [arXiv:1012.1314]

  8. [16]

    Hatta, J

    Y. Hatta, J. Noronha, and Bo-Wen Xiao, Phys. Rev. D89, 051702 (2014) [arXiv:1401.6248]

  9. [17]

    Hatta, J

    Y. Hatta, J. Noronha, and Bo-Wen Xiao, Phys. Rev. D89, 114011 (2014) [arXiv:1403.7693]

  10. [18]

    Cs¨ org˝ o, F

    T. Cs¨ org˝ o, F. Grassi, Y. Hama, T. Kodama, Phys. Lett. B565, 107 (2003) [arXiv:nucl-th/0305059]

  11. [19]

    Cs¨ org˝ o, L

    T. Cs¨ org˝ o, L. P. Csernai, Y. Hama and T. Kodama, Acta Phys. Hung. A 21, 73-84 (2004) [arXiv:nucl-th/0306004]

  12. [20]

    Cs¨ org˝ o, M

    T. Cs¨ org˝ o, M. I. Nagy, and M. Csan´ ad, Phys. Lett. B663, 306 (2008) [arXiv:nucl-th/0605070]

  13. [21]

    M. I. Nagy, T. Cs¨ org˝ o, and M. Csan´ ad, Phys. Rev. C77, 024908(2008) [arXiv:0709.3677]

  14. [22]

    Cs¨ org˝ o, G

    T. Cs¨ org˝ o, G. Kasza, M. Csan´ ad, Z. F. Jiang, Universe. 4 6, 69 (2018) [arXiv:1805.01427]

  15. [23]

    Cs¨ org˝ o, G

    T. Cs¨ org˝ o, G. Kasza, M. Csan´ ad, Z. F. Jiang, Acta Phys. Polon. B 50, 27 (2019) [arXiv:1806.06794]

  16. [24]

    Eckart, Phys

    C. Eckart, Phys. Rev. 58, 919 (1940)

  17. [25]

    L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- amon Press, London, 1959)

  18. [26]

    Israel, and J

    W. Israel, and J. M. Stewart, Annals Phys. 118, 341-372 (1979)

  19. [27]

    H. B. Meyer, Phys. Rev. Lett. 100, 162001 (2008) [arXiv:0710.3717]

  20. [28]

    Kharzeev, and K

    D. Kharzeev, and K. Tuchin, JHEP 09, 093 (2008) [arXiv:0705.4280]

  21. [29]

    Karsch, D

    F. Karsch, D. Kharzeev and K. Tuchin, Phys. Lett. B 663, 217 (2008) [arXiv:0711.0914]

  22. [30]

    Noronha-Hostler, J

    J. Noronha-Hostler, J. Noronha and C. Greiner, Phys. Rev. Lett. 103, 172302 (2009) [arXiv:0811.1571]

  23. [31]

    Harutyunyan and A

    A. Harutyunyan and A. Sedrakian, Particles 1, no. 1, 212 (2018) [arXiv:1808.06403]

  24. [32]

    S. Ryu, J. F. Paquet, C. Shen, G. Denicol, B. Schenke, S. Jeon and C. Gale, Phys. Rev. C 97, no. 3, 034910 (2018) [arXiv:1704.04216]

  25. [33]

    P. B. Arnold, C. Dogan and G. D. Moore, Phys. Rev. D 74, 085021 (2006) [arXiv:hep-ph/0608012]

  26. [34]

    J. E. Bernhard, J. S. Moreland, S. A. Bass, J. Liu and U. Heinz, Phys. Rev. C 94, no. 2, 024907 (2016) [arXiv:1605.03954]

  27. [35]

    J. E. Bernhard, J. S. Moreland and S. A. Bass, Nature Phys. 15, no. 11, 1113 (2019)

  28. [36]

    Rougemont, R

    R. Rougemont, R. Critelli, J. Noronha-Hostler, J. Noronha and C. Ratti, Phys. Rev. D 96, no. 1, 014032 (2017) [arXiv:1704.05558]

  29. [37]

    Alqahtani, M

    M. Alqahtani, M. Nopoush and M. Strickland, Prog. Part. Nucl. Phys. 101, 204 (2018) [arXiv:1712.03282]

  30. [38]

    Kovtun, Dan T

    P. Kovtun, Dan T. Son, Andrei O. Starinets, Phys. Rev. Lett. 94, 111601 (2005) [arXiv:hep-th/0405231]

  31. [39]

    Benincasa, A

    P. Benincasa, A. Buchel, A. Starinets, and O. An- drei, Nucl. Phys. B733, 160-187 (2006) [arXiv:hep- th/0507026]

  32. [40]

    I. M. Khalatnikov, Sov. Phys. JETP 2, 169 (1956)

  33. [41]

    Csan´ ad, T

    M. Csan´ ad, T. Cs¨ org˝ o, J. Phys. G30, S1079-S1082 (2004). [arXiv:nucl-th/0403074]

  34. [42]

    Csan´ ad, M

    M. Csan´ ad, M. I. Nagy, and S. L¨ ok¨ os, Eur. Phys. J. A48, 173 (2012) [arXiv:1205.5965]

  35. [43]

    Chojnacki, W

    M. Chojnacki, W. Florkowski and T. Cs¨ org˝ o, Phys. Rev. C 71, 044902 (2005) [arXiv:nucl-th/0410036]

  36. [44]

    Csan´ ad and M

    M. Csan´ ad and M. Vargyas, Eur. Phys. J. A44, 473-478 (2010) [arXiv:0909.4842]

  37. [45]

    Csan´ ad and A

    M. Csan´ ad and A. Szab´ o, Phys. Rev. C90, no.5, 054911 (2014) [arXiv:1405.3877]

  38. [46]

    New, multipole solutions of relativistic, viscous hydrodynamics

    T. Cs¨ org˝ o and G. Kasza, “New, multipole solutions of relativistic, viscous hydrodynamics” [arXiv:2003.08859]

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.