Pith. sign in

REVIEW 2 major objections 3 minor 66 references

Relativistic dissipative fluids in the trace-fixed particle frame: Strongly hyperbolic quasi-linear first-order evolution equations

T0 review · 2 major / 3 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read The paper proves that the trace-fixed particle frame gives a first-order relativistic dissipative fluid theory that is locally well-posed in the nonlinear regime, provided a single transport inequality and a non-crossing condition on…

desk verdict A real new first-order dissipative fluid theory with a sound local well-posedness theorem, though the abstract and conclusions quietly drop the non-crossing condition on characteristic speeds that the proof requires. read the letter →

arxiv 2412.03713 v2 pith:NL3LROSZ submitted 2024-12-04 gr-qc hep-thmath.AP

classification gr-qchep-thmath.AP MSC 35L4035L4576Y05 PACS 47.75.+f05.70.Ln
keywords relativisticdissipativefluidstrace-fixedparticleframestronghyperbolicityfirst-orderquasilinearsystemwell-posedCauchyproblemconstraintpropagationcharacteristicspeedstransportcoefficients
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

Relativistic dissipative fluid theories often fail to have a well-posed initial-value problem, which blocks both rigorous predictions and numerical simulation. This paper proves that a particular first-order theory, the trace-fixed particle frame, is locally well-posed in the full nonlinear regime: for sufficiently smooth initial data there is a unique solution that depends continuously on the data. The proof works by turning the mixed-order evolution equations into a constrained first-order quasilinear system, showing that system is strongly hyperbolic through an explicit symmetrizer, and proving that the auxiliary constraints propagate in time. The only conditions needed are the inequality $1+2\eta/\kappa \le e/(k_B T)$ and that the three characteristic speeds remain distinct. This makes the theory a concrete candidate for modeling strongly viscous relativistic fluids such as merging neutron stars or quark-gluon plasma.

What carries the argument

The central machinery is the first-order quasilinear system (61–68) and its principal symbol $\mathcal{A}(k,U)=\mathcal{A}^\mu(U)k_\mu$ for covectors $k$ orthogonal to the fluid four-velocity. The symbol is decomposed into scalar, vector, and tensor blocks along $k$; choosing the constraint-addition coefficients $\delta_1,\delta_2,\delta_3$ according to (107) makes the scalar and vector blocks diagonalizable with real eigenvalues. The proof of strong hyperbolicity then rests on a smooth symmetrizer $H(k,U)$ built from the eigenprojectors, whose smoothness follows because the eigenvalues $\mu_0,\mu_1,\mu_2$ are functions of temperature alone and are assumed distinct. Constraint propagation is shown by embedding the constraint fields in a larger system with a $k$-independent symmetrizer $H_c$, which is symmetric hyperbolic.

What would settle it

For a concrete gas model, such as the hard-sphere gas in $d=3$ or hard-disk gas in $d=2$ discussed in the companion letter, evaluate the functions $\mu_0,\mu_1,\mu_2$ from Eq. (38) over the full temperature range. If any two of them intersect at a finite temperature, the hypothesis of distinctness fails, and the smoothness of the symmetrizer would need to be re-examined; numerics could then test whether well-posedness still holds. Alternatively, discretize the constrained first-order system (61–68) with a standard stable method; if small constraint violations grow without bound as the grid is refined, the constraint-propagation proof contradicts the numerics.

Watch

Extended reading notes

Core claim

The central claim is Theorem 1: under smooth, strictly positive temperature-dependent transport coefficients $e,\eta,\kappa,\zeta,\Gamma_1$ with $0<c_v<d k_B$ and $\Gamma_1\neq 1$, with $\Gamma_2=h/(k_B T)$, inequality (1), and distinct characteristic speeds $\mu_0,\mu_1,\mu_2$, the nonlinear system (28–32) admits a unique local solution depending continuously on the initial data. The theorem is proved by embedding (28–32) in a first-order quasilinear system with auxiliary fields $N_\mu,T_\mu,B_{\mu\nu}$ and constraint fields $C^{(N)}_\mu,C^{(T)}_\mu,C^{(B)}_{\mu\nu}$. Constraint-violating terms are added off the constraint surface to make the principal symbol diagonalizable with real eigenvalues; a smooth symmetrizer is constructed from the eigenprojectors using the fact that the eigenvalues depend only on the temperature. Finally, the constraint fields themselves are shown to satisfy a symmetric hyperbolic system, so constraints imposed initially remain satisfied forever.

Load-bearing premise

The proof requires that the three characteristic speeds $\mu_0,\mu_1,\mu_2$ never coincide for any temperature; if two of these curves cross, the eigenprojectors may lose smoothness and the explicit symmetrizer construction breaks down.

Editorial extensions

If this is right

  • The nonlinear evolution equations of the trace-fixed particle frame become locally well-posed; numerical simulations of dissipative relativistic fluids (e.g., neutron star mergers) can rely on a well-defined continuum problem.
  • Well-posedness, causality, and stability at equilibrium are achieved with a single transport inequality, simpler than the multi-parameter conditions of previous first-order theories.
  • The constrained first-order system (61–68) gives an explicit practical form for numerical implementation, with evolution equations for expansion, shear, and vorticity.
  • The proof provides a template for treating other first-order dissipative fluid theories: rewrite as a constrained system, make it strongly hyperbolic off the constraint surface by adding constraint terms, and prove constraint propagation.
  • Because the symmetrizer depends only on the temperature through the characteristic speeds, the result immediately extends to any fixed globally hyperbolic spacetime background.

Reading between the lines

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

  • If realistic equations of state produce crossings of the characteristic speeds, the theory might still be well-posed, but the present proof technique would need a more general symmetrizer that tolerates eigenvalue multiplicities or a direct argument showing crossings do not occur in the physically relevant regime.
  • The inequality $1+2\eta/\kappa \le e/(k_B T)$ and the bound $c_v < d k_B$ tie the theory's well-posedness to thermodynamics; testing these against tabulated transport coefficients for actual gases could map the theory's domain of validity.
  • The same constraint-addition strategy could be applied to the full Einstein-fluid system, which the article lists as an open problem; the fixed-background proof here provides the local machinery that such a coupled system would need.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 3 minor

Summary. This paper develops a first-order relativistic dissipative-fluid theory in the trace-fixed particle frame. The authors derive constitutive relations (22)-(24) with coefficients Gamma1 and Gamma2, rewrite the equations as a constrained first-order quasilinear system (61)-(68), and prove that, under inequality (1), the choice Gamma2 = h/(k_B T), and a distinctness condition on the characteristic speeds mu0, mu1, mu2, the principal symbol can be diagonalized by an explicit block decomposition with a smooth symmetrizer (Sec. V). They also show that the auxiliary constraints propagate by embedding them in a symmetric hyperbolic system (Sec. VI). The main result is Theorem 1, a local well-posedness statement for the Cauchy problem.

Significance. The conditional theorem is technically substantial. The block decomposition into scalar, vector, and tensor modes, the explicit choice of the constraint-addition coefficients delta_i, and the reduction of constraint propagation to symmetric hyperbolicity are concrete and checkable. The paper also makes a useful comparison with BDNK-type frame choices and identifies which inequalities are needed. However, the advertised result is broader than the theorem: the abstract and conclusions omit the non-crossing condition on mu0, mu1, mu2, and the theorem does not establish that condition for any physical equation of state, including the hard-sphere gases verified for inequality (1) in the companion letter. This gap, together with the closure issue for the electromagnetic term in Eq. (67), requires a major revision before the claims match the proof.

major comments (2)
  1. [Theorem 1; Secs. V.E, VII; Abstract] The distinctness assumption on mu0, mu1, mu2 is essential for the proof. In Sec. V.E the symmetrizer (112) is built from eigenprojectors Pj, and the cited smoothness result requires constant multiplicity; the same point is made in App. D. The theorem therefore only proves strong hyperbolicity and local well-posedness under this non-crossing condition. Yet the abstract and the conclusions state, without this qualification, that the full nonlinear system is strongly hyperbolic and yields a well-posed Cauchy problem. The companion-letter verification recalled in Sec. III(d) checks inequality (1) for hard-sphere and hard-disk gases, but it does not check that mu0, mu1 and mu2 remain distinct. Please either verify the non-crossing condition for the physical examples, extend the proof to eigenvalue crossings, or state the condition explicitly in the abstract and conclusions and restrict the advertised claims accordingly.
  2. [Eq. (67) and Sec. V.A] The first-order system (61)-(68) is claimed to be quasilinear in U, but Eq. (67) contains the term D_mu E_nu. Since E_nu = u^alpha F_{nu alpha}, this term involves D_mu u^alpha, i.e. a derivative of a state variable that is not among the first-order variables unless it is replaced by B_mu^alpha up to the constraint C(B)_(mu alpha). The principal-symbol computation in Sec. V.A, Eq. (77), contains no contribution of this type, so the system as written seems not to be closed in first-order form. Please clarify how D_mu E_nu is expressed in terms of U and its derivatives, or state explicitly that terms involving derivatives of the electromagnetic field are neglected because E_mu is O(gradient) as in footnote 1, and make the theorem's hypotheses on F explicit.
minor comments (3)
  1. [Eq. (71)] In the vorticity evolution equation, the term beta2(alpha4 theta + alpha4 epsilon) should presumably read beta2(alpha4 theta + alpha5 epsilon), since dot T/T = alpha4 theta + alpha5 epsilon in Eq. (29).
  2. [Sec. III, property (d)] Property (d) says that the inequality and technical assumptions (i)-(iv) hold for hard-sphere gases, but the distinctness hypothesis of Theorem 1 is not listed there; this should be either verified and added, or explicitly identified as an open condition for those examples.
  3. [Sec. II, Eq. (25)] The comparison with BDNK theory around Eq. (25) would be easier to check if the derivation of Eq. (25) from Eqs. (11) of Ref. [17] were shown, since this equation is used to justify the difference between the two frames.

Circularity Check

0 steps flagged · score 0.0 of 10

The derivation is self-contained: the theorem follows from an explicit symmetrizer construction and constraint-propagation argument; reliance on the companion letter is ordinary citation, not circularity, and the distinctness caveat is a scope condition, not a circular reduction.

full rationale

The claimed result is a conditional mathematical theorem, not a fitted prediction. Starting from the Eckart constitutive relations, Sec. II obtains the TFP-frame relations (22)-(24) by a frame change and by adding combinations of the Euler equations, which is a legitimate representation freedom. The evolution system (28)-(32) is obtained by algebraically solving these relations for time derivatives, and the first-order reformulation (61)-(68) introduces gradients as new fields together with constraints. The strong-hyperbolicity proof in Sec. V constructs an explicit diagonalizer by choosing the constraint-addition coefficients δ1,δ2,δ3 as in Eq. (107); this is an engineering of free coefficients, not a fit of the theorem's conclusion. The constraint-propagation proof in Sec. VI is self-contained with the explicit symmetrizer (124). The eigenvalue bounds 0<µ1<µ2≤1 used in Sec. V.D are quoted from the authors' companion letter, but that is a separate, parameter-free result with stated assumptions (inequality (1) and Γ2=h/(kBT)); it does not assume the nonlinear well-posedness theorem being proved, so it is independent support rather than circular self-citation. The unverified hypothesis that µ0,µ1,µ2 remain distinct is a scope gap in the advertised theorem, but it is not a circular reduction of the conclusion to the assumptions.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The central claim rests on the chosen frame (matching J^µ and T^µ_µ), the ideal gas equation of state, the inequality (1), the non-crossing condition, and the standard theory of symmetric hyperbolic PDEs. The δi coefficients are construction parameters, not physical inputs, and are set to make the principal symbol diagonalizable.

free parameters (5)
  • Γ1 (representation coefficient) = unspecified; chosen large enough via Eq. (41) with Λ0
    Free dimensionless coefficient in constitutive relation (22). Does not affect the principal symbol; chosen to ensure stability in the companion letter.
  • Γ2 (representation coefficient) = h/(k_B T)
    Fixed by Eq. (37) to obtain strong hyperbolicity and simplify β2 = 0.
  • δ1 (constraint addition coefficient) = -2η/κ
    Chosen in Eq. (107) to make the vector block diagonalizable.
  • δ2 (constraint addition coefficient) = 2aη/κ, with a = α4 + 2(d-1)α5η
    Chosen to satisfy det Q_|| = 0, making the scalar block diagonalizable.
  • δ3 (constraint addition coefficient) = -2η
    Chosen in Eq. (97) to eliminate the off-diagonal term in the vector block.
assumptions (6)
  • domain assumption Ideal gas equation of state p = n k_B T
    Used in App. A.2, Eqs. (A30-A32), and throughout the derivation of constitutive relations; restricts the theory to ideal gases.
  • domain assumption Heat capacity condition c_v < d k_B
    Sec. II.B: required for the temperature to be uniquely determined from the trace-fixing condition.
  • domain assumption Inequality (1): 1 + 2η/κ ≤ e/(k_B T)
    Assumed in Theorem 1; used to prove causality and 0 < µ1 < µ2 ≤ 1. Shown in the companion letter for hard-sphere and hard-disk gases.
  • ad hoc to paper Characteristic speeds µ0, µ1, µ2 remain distinct for all T > 0
    Assumed before Theorem 1; needed for smooth eigenprojectors and a smooth symmetrizer in Sec. V.E via Kato's perturbation theory.
  • standard math Standard theory of symmetric hyperbolic PDE systems
    Invoked in Sec. V and App. D (Kreiss-Lorenz, Taylor) to conclude local well-posedness from strong hyperbolicity.
  • standard math Geroch-Reula covariant principal symbol formalism
    Used in Sec. V and App. D to analyze hyperbolicity without choosing a spacetime foliation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Relativistic dissipative fluids in the trace-fixed particle frame: Strongly hyperbolic quasi-linear first-order evolution equations." pith.science (2026). https://pith.science/paper/NL3LROSZ

@misc{pith2026241203713,
  author       = {Pith},
  title        = {Pith review of: Relativistic dissipative fluids in the trace-fixed particle frame: Strongly hyperbolic quasi-linear first-order evolution equations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NL3LROSZ}},
  note         = {Machine review of arXiv:2412.03713}
}
read the original abstract

In this paper we derive a new first-order theory of relativistic dissipative fluids by adopting the trace-fixed particle frame. Whereas in a companion letter we show that this theory is hyperbolic, causal and stable at global equilibrium states, here we prove that the full nonlinear system of equations can be cast into a first-order quasilinear system which is strongly hyperbolic. By rewriting the system in first-order form, auxiliary constraints are introduced. However, we show that these constraints propagate, and thus our theory leads to a well-posed Cauchy problem.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

66 extracted references · 61 canonical work pages

  1. [1]

    is satisfied, a technical difficulty arises when considering the nonlinear problem. This is due to the fact th at the velocity field is not hypersurface orthogonal in general , and thus when introducing local coordinates and writing the system as partial differential equations (PDEs), second or der time derivatives of the velocity field appear. To circumvent ...

  2. [2]

    together with ( 18), ( 19), coupled to the constitutive relations obtained in the previ ous subsection, i.e. Eqs. ( 22)-(24). The usual procedure con- sists in substituting the latter in Eqs. ( 2) to obtain an evolu- tion system for the state variables (n, T, u µ). This leads to a set of equations that involves second-order time derivati ves of T and uµ. ...

  3. [3]

    ( A9-A11) and the addition of multiples of the left-hand sides of Eqs

    Thermodynamic constraints Up to now, the number of coefficients involved in the gen- eral constitutive relations was reduced from 18 to 5 by means of considering two freedoms that are allowed within the first - order scheme: a change of frame given by Eqs. ( A9-A11) and the addition of multiples of the left-hand sides of Eqs. ( A30- A32). These transformati...

  4. [4]

    Furthermore, one can use the hatted transformations to achieve π1 = π2 = 0 and κ1 = 0

    Summary In order to shed light on the physical meaning of the in- variant coefficients defined so far, it is instructive to writ e down the constitutive relations in the Eckart frame, for whi ch νi = εi = 0 and γj = 0 (note that this frame can always be at- tained by means of the transformation ( A22– A27)), such that fi = πi and ℓj = −κj/h. Furthermore, on...

  5. [5]

    sufficiently regular

    implies β1 = −kBT /e, β2 = 0 , and β3 = −β1/κ, and it eliminates the gradient of the tem- perature on the right-hand side of Eq. ( 23), leading to Qµ = κ ïDµn n + e kBT aµ − q kBT Eµ ò . (42) The main objective of this article is to prove that, under the same hypothesis as in (a) and the assumption that the func- tions µ0, µ1, and µ2 defined in (b) do not ...

  6. [6]

    B. P . Abbott, R. Abbott, T. D. Abbott, and F. A. et al. (LIGO Scientific Collaboration, the Virgo Collaboration an d the KAGRA Collaboration), Observation of Gravitational Waves from Two Neutron Star–Black Hole Coalescences, The Astrophysical Journal Letters 915, L5 (2021)

  7. [7]

    Chabanov, L

    M. Chabanov, L. Rezzolla, and D. H. Rischke, General- relativistic hydrodynamics of non-perfect fluids: 3+1 cons erva- tive formulation and application to viscous black hole accretion, Monthly Notices of the Royal Astronomical Society 505, 5910 (2021)

  8. [8]

    M. J. Hatton and I. Hawke, A dissi- pative extension to ideal hydrodynamics, Monthly Notices of the Royal Astronomical Society 535, 47 (2024)

Show all 66 references
  1. [9]

    V´ an, Nonequilibrium thermody- namics: emergent and fundamental, Philosophical Transactions of the Royal Society A: Mathema tical, Physical and

    P . V´ an, Nonequilibrium thermody- namics: emergent and fundamental, Philosophical Transactions of the Royal Society A: Mathema tical, Physical and

  2. [10]

    J. F. Salazar and T. Zannias, On extended thermo- dynamics: From classical to the relativistic regime, International Journal of Modern Physics D 29, 2030010 (2020)

  3. [11]

    Explicitly, they are given by µ1,2 = 1 2 Ä T r ± √ T r2 − 4D ä , (103) where D = 1 d kBT e , T r = 1 d + kBT e Å 1 + A0 κ ã

    and the inequality ( 1) hold. Explicitly, they are given by µ1,2 = 1 2 Ä T r ± √ T r2 − 4D ä , (103) where D = 1 d kBT e , T r = 1 d + kBT e Å 1 + A0 κ ã . (104) According to the proof of Lemma 1, two distinct nonzero eigenvalues of ˜M ‖ give rise to four linearly independent ...

  4. [12]

    G. S. Rocha, D. Wagner, G. S. Denicol, J. Noronha, and D. H. Rischke, Theories of relativistic dissipative fluid dynami cs, En- tropy 26, 10.3390/e26030189 (2024)

  5. [13]

    F. S. Bemfica, M. M. Disconzi, and J. Noronha, Causality and existence of solutions of relativistic viscous fluid dynami cs with gravity, Phys. Rev. D 98, 104064 (2018)

  6. [14]

    In particular, V = 0 is the unique solution with trivial initial data, and hence the constraints propag ate as desired

    is sym- metric hyperbolic. In particular, V = 0 is the unique solution with trivial initial data, and hence the constraints propag ate as desired. VII. CONCLUSIONS In this work and the companion letter [24] we presented and analyzed a novel theory for relativistic dissipative ...

  7. [15]

    Kovtun, First-order relativistic hydrodynamics is stab le, JHEP 10, 034

    P . Kovtun, First-order relativistic hydrodynamics is stab le, JHEP 10, 034. 19

  8. [16]

    A specific choice is referred to as a “frame” in the literature

    First-order changes of frame The determination of (n, T, uµ) from a given, off- equilibrium configuration is not unique. A specific choice is referred to as a “frame” in the literature. In the following, we consider a first-order change of the frame, that is, a transfo r- mation ...

  9. [17]

    Changes of representation As mentioned before, there is another freedom in the trans- port coefficients that should be taken into account. Indeed, the equations of motion ( 2) imply that ˙n n + θ = O(∂2), (A30) ˙T T + kB cv θ = O(∂2), (A31) aµ + 1 nh Dµp − q h F µν uν = O(∂2), ...

  10. [18]

    R. E. Hoult and P . Kovtun, Causal first-order hy- drodynamics from kinetic theory and holography, Phys. Rev. D 106, 066023 (2022)

  11. [19]

    G. S. Rocha, G. S. Denicol, and J. Noronha, Perturbative ap- proaches in relativistic kinetic theory and the emergence of first- order hydrodynamics, Phys. Rev. D 106, 036010 (2022)

  12. [20]

    time evolution vector field

    Strongly hyperbolic first-order systems We first recall some important definitions and results from Ref. [36] which are independent of a given spacetime folia- tion. Consider a first-order system of the form A µ ∇µU = F , (D1) where the components of the m × m matrix A µ and the m...

  13. [21]

    In contrast , we first perform a change of frame before adding these combina- tions

    directly to the constitutive relations in the Eckart frame. In contrast , we first perform a change of frame before adding these combina- tions. Moreover, when adding the combinations ( 21), we do it in such a fashion to preserve the TFP frame which leaves us with the two free ...

  14. [22]

    In order to do so we define A µ := uµI − Aµ(U ) with I denoting the identity ma- trix and the matrices Aµ(U ) satisfying Aµ(U )uµ = 0

    Strongly hyperbolic PDEs for the fluid system After the remarks made in the previous subsection we are ready to apply the theory to our system ( 72) for the fluid equa- tions, for which m = ( d + 1)(d + 3). In order to do so we define A µ := uµI − Aµ(U ) with I denoting the ident...

  15. [23]

    Causality Finally, we remark that the covectors ξµ for which A(ξ) has a nontrivial kernel are of the form ξµ = −λuµ + kµ, u µkµ = 0, k µkµ = 1, (D10) with λ the eigenvalues computed in Sec. V. Since the con- ditions for causality imply that |λ| ≤ 1, it follows that ξµξµ = 1 − ...

  16. [24]

    Derradi de Souza, T

    R. Derradi de Souza, T. Koide, and T. Kodama, Hy- drodynamic approaches in relativistic heavy ion reactions , Progress in Particle and Nuclear Physics 86, 35 (2016)

  17. [25]

    Du and U

    L. Du and U. Heinz, (3+1)-dimensional dissipative rel- ativistic fluid dynamics at non-zero net baryon density, Computer Physics Communications 251, 107090 (2020)

  18. [26]

    Shen and L

    C. Shen and L. Yan, Recent development of hy- drodynamic modeling in heavy-ion collisions, Nuclear Science and Techniques 31, 122 (2020)

  19. [27]

    B. P . Abbott, R. Abbott, T. D. Abbott, and F. A. et al. (LIGO Scientific Collaboration and Virgo Collaboration), GW1708 17: Observation of Gravitational Waves from a Binary Neutron Star Inspiral, Phys. Rev. Lett. 119, 161101 (2017)

  20. [28]

    B. P . Abbott, R. Abbott, T. D. Abbott, and F. A. et al, Gravitational Waves and Gamma-Rays from a Bi- nary Neutron Star Merger: GW170817 and GRB 170817a, The Astrophysical Journal Letters 848, L13 (2017)

  21. [29]

    Hawking and G

    S. Hawking and G. Ellis, The Large Scale Structure of Space Time (Cambridge University Press, Cambridge, 1973)

  22. [30]

    R. M. Wald, General Relativity (Chicago Univ. Pr., Chicago, USA, 1984)

  23. [31]

    G. F. R. Ellis, R. Maartens, and M. A. H. MacCallum, Relativis- tic Cosmology (Cambridge University Press, 2012)

  24. [32]

    Abalos and O

    F. Abalos and O. Reula, On necessary and sufficient con- ditions for strong hyperbolicity in systems with constrain ts, Class. Quant. Grav. 37, 185012 (2020)

  25. [33]

    Gavassino and M

    L. Gavassino and M. Antonelli, Unified extended irreversibl e thermodynamics and the stability of relativistic theories for dissipation, Frontiers in Astronomy and Space Sciences 8, 10.3389/fspas.2021.686344 (2021)

  26. [34]

    Abalos, O

    F. Abalos, O. Reula, and D. Hilditch, Hyperbolic extensions of constrained PDEs, (2024), arXiv:2410.18286 [math.AP]

  27. [35]

    Geroch, Partial differential equations of physics, Gene ral Relativity: Proceedings

    R. Geroch, Partial differential equations of physics, Gene ral Relativity: Proceedings. Edited by G.S. Hall and J.R. Pulha m. Edinburgh, IOP Publishing , 19 (1996), arXiv:gr-qc/9602055

  28. [36]

    F. S. Bemfica, M. M. Disconzi, and J. Noronha, Nonlinear causality of general first-order relativistic viscous hydr odynam- ics, Phys. Rev. D 100, 104020 (2019)

  29. [37]

    This means that global equilibrium solutions are mode sta- ble

    holds and Γ 1 = 1 + cv kB e2 kBT h Å kB cv − 1 d ã2 κ ζ Λ 0, (41) with large enough values of the constant Λ 0. This means that global equilibrium solutions are mode sta- ble. (d) The inequality (1) and the technical assumptions (i)-(iv) are automatically satisfied for a simple...

  30. [38]

    R. E. Hoult and P . Kovtun, Stable and causal relativistic Navier-Stokes equations, Journal of High Energy Physics 2020, 67 (2020)

  31. [39]

    F. S. Bemfica, M. M. Disconzi, and J. Noronha, First-order general-relativistic viscous fluid dynamics, Phys. Rev. X 12, 021044 (2022)

  32. [40]

    Taylor,Partial Differential Equations III: Nonlinear Equations , 2nd ed., Applied Mathematical Sciences, V ol

    M. Taylor,Partial Differential Equations III: Nonlinear Equations , 2nd ed., Applied Mathematical Sciences, V ol. 117 (Springer , New Y ork, 1996)

  33. [41]

    Sombras, lentes y ondas gravitatorias generadas por objetos compactos as- trof´ ısicos

    and (37), strong hyperbolicity and causality only require the fulfillment of the single inequality ( 1). Furthermore, this condition is indepen- dent of the only arbitrary parameter Λ 0 in our system, whose sole purpose is to achieve the stability property. Therefor e, we belie...

  34. [42]

    Disconzi, Recent developments in mathematical aspects o f relativistic fluids, Living Reviews in Relativity 27, 6 (2024)

    M. Disconzi, Recent developments in mathematical aspects o f relativistic fluids, Living Reviews in Relativity 27, 6 (2024)

  35. [43]

    Hiscock and L

    W. Hiscock and L. Lindblom, Generic instabilities in first-order dissipative relativistic fluid theories, Phys. Rev. D 31, 725 (1985)

  36. [44]

    Ciambelli and L

    L. Ciambelli and L. Lehner, Fluid-gravity correspondence and causal first-order relativistic viscous hydrodynamics , Phys. Rev. D 108, 126019 (2023)

  37. [45]

    (56) The derivatives of uµ that appear in k[µν] can be replaced by B[µν] using the constraint field C(B) µν

    one finds D[µC(N ) ν] = D[µNν] − k[µν] ˙n n , (54) D[µC(T ) ν] = D[µTν] − k[µν] ˙T T , (55) D[µC(B) ν]α = D[µBν]α − k[µν]aα − 1 2 ∆ µµ′ ∆ ν ν′ Rµ′ν′αβuβ. (56) The derivatives of uµ that appear in k[µν] can be replaced by B[µν] using the constraint field C(B) µν . This gives rise...

  38. [46]

    J. M. Stewart, Non-equilibrium relativistic kinetic theor y, in Non-Equilibrium Relativistic Kinetic Theory (Springer Berlin Heidelberg, Berlin, Heidelberg, 1971) pp. 1–113

  39. [47]

    J. F. Salazar, A. L. Garc´ ıa-Perciante, and O. Sar- bach, Relativistic dissipative fluids in the trace-fixed particle frame: Hyperbolicity, causality, and stability, Phys. Rev. D 111, L081501 (2025)

  40. [48]

    Gabarrete, A

    C. Gabarrete, A. L. Garc´ ıa-Perciante, and O. Sarbach, In prepa- ration

  41. [49]

    H. O. Kreiss and J. Lorenz, Initial-boundary value problems and the Navier-Stokes equations (Academic Press, San Diego, 1989)

  42. [50]

    G. Nagy, O. Ortiz, and O. Reula, Strongly hyperbolic second order Einstein’s evolution equations, Phys. Rev. D70, 044012 (2004)

  43. [51]

    Sarbach and M

    O. Sarbach and M. Tiglio, Continuum and Discrete Initial- Boundary-V alue Problems and Einstein’s Field Equations, Living Rev. Rel. 15, 9 (2012)

  44. [55]

    Abalos, On constraint preservation and strong hyperboli city, Class

    J. Abalos, On constraint preservation and strong hyperboli city, Class. Quant. Grav. 39, 215004 (2022)

  45. [57]

    Using the com- mutator identity ( C5) in App

    of Z (N ) µν and take a dot on both sides of this equation. Using the com- mutator identity ( C5) in App. C and Eqs. (113) and (115) one first obtains ˙Z (N ) µν = −N βD[µC(B) ν]β + k[µ αDν]C(N ) α − k[µ|α|DαC(N ) ν] + δ1D[µZν] + l.o., (116) 12 where here and from now on “ l.o....

  46. [58]

    Reula, Strongly hyperbolic systems in general relativit y, Diff

    O. Reula, Strongly hyperbolic systems in general relativit y, Diff. Eq. 01, 251 (2004)

  47. [59]

    Eckart, The thermodynamics of irreversible pro- cesses

    C. Eckart, The thermodynamics of irreversible pro- cesses. iii. relativistic theory of the simple fluid, Phys. Rev. 58, 919 (1940)

  48. [60]

    A. L. Garc´ ıa-Perciante, A. R. M´ endez, and O. Sar- bach, Existence of the Chapman-Enskog solution and its relation with first-order dissipative fluid theories, Journal of Non-Equilibrium Thermodynamics 50, 295 (2025)

  49. [61]

    Israel, Relativistic kinetic theory of a simple gas, J

    W. Israel, Relativistic kinetic theory of a simple gas, J. Ma th. Phys. 4, 1163 (1963)

  50. [63]

    Sarbach, E

    O. Sarbach, E. Barausse, and J. Preciado-L´ opez, Well- posed Cauchy formulation for Einstein-æther theory, Class. Quant. Grav. 36, 165007 (2019)

  51. [64]

    Kato, Perturbation Theory for Linear Operators (Springer- V erlag, New Y ork, 1980)

    T. Kato, Perturbation Theory for Linear Operators (Springer- V erlag, New Y ork, 1980)

  52. [65]

    Israel and J

    W. Israel and J. Stewart, Transient relativistic thermodyn amics and kinetic theory, Annals of Physics 118, 341 (1979)

  53. [66]

    Israel, Covariant fluid mechanics and thermodynamics: An introduction, in Relativistic Fluid Dynamics , edited by A

    W. Israel, Covariant fluid mechanics and thermodynamics: An introduction, in Relativistic Fluid Dynamics , edited by A. M. Anile and Y . Choquet-Bruhat (Springer Berlin Heidelberg, Berlin, Heidelberg, 1989) pp. 152–210

  54. [67]

    Conse- quently, only derivatives of the heat flux appear in this equa - tion, which simplifies the principal part of the equations

    become clear: they elim- inate the antisymmetric derivatives of Nµ and Tµ that would otherwise appear in the evolution equation for ωµν. Conse- quently, only derivatives of the heat flux appear in this equa - tion, which simplifies the principal part of the equations. V . STRONG...

  55. [72]

    For this, one consi d- ers the principal symbol, defined as A(k, U ) := Aµ(U )kµ, (73) where kµ is a given covector perpendicular to uµ

    without rewriting it explicitly as a PDE system, and this greatly simplifies the analysis. For this, one consi d- ers the principal symbol, defined as A(k, U ) := Aµ(U )kµ, (73) where kµ is a given covector perpendicular to uµ. The first- order system ( 72) is called (cf. Definiti...

  56. [74]

    implies that A(k, U ) is sym- metric with respect to the scalar product defined by H(k, U ) and hence it is diagonalizable and has only real eigenvalues . Conversely, if A(k, U ) is diagonalizable and has a real spec- trum, and if S(k, U ) denotes the matrix whose columns are t...

  57. [97]

    (98) Furthermore, the system is causal if 0 < β 3η ≤ 1

    holds, the vector block is diagonalizable with the purely real eigenvalues: 0, ± √ β3η. (98) Furthermore, the system is causal if 0 < β 3η ≤ 1. This is the same condition that was found in Eq. (20) of the compan- ion letter [24] when analyzing the vector block of the second - ...

  58. [122]

    is symmetric hyperbolic. Indeed, it is simple to verify that the symmetri c, positive definite matrix Hc(U0) := â 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 2 s ì (124) is a symmetrizer for Ac(U0, k), that is, Hc(U0)Ac(U0, k) is symmetric for all U0 a...

Pith tools

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