{"id":"3687c5c6-5491-4f0b-9c24-6804439f6679","arxiv_id":"2611.09582","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Two new flux-conservative formulations of BDNK viscous hydrodynamics are proved well-posed and validated on flat and curved spacetimes, including heavy-ion and black hole settings.","lead":"This paper constructs two new numerical ways to solve BDNK, a causal first-order theory of relativistic viscous fluids, on curved spacetimes, with proofs that the equations are well-posed. The simulations reproduce known flows and show when BDNK and the older Israel-Stewart approach agree or diverge for heavy-ion collisions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"Assumption IV.1 is the load-bearing restriction: the frame prescription does not guarantee it holds throughout nonconformal evolutions, so the simulated trajectories may leave the proven regime.","rationale":"The reader's weakest assumption and my primary concern coincide: Assumption IV.1 is genuinely load-bearing, and the available text does not establish that it is satisfied along the evolutions used to support the physical conclusions. The proofs use Assumption IV.1 in indispensable ways: it guarantees det M ≠ 0, it is used for real distinct roots of the principal symbol and strict hyperbolicity in Step 4 of Section IX and in Section X, and it underpins the equivalence arguments in Section VIII. The inequality conditions depend on ε and n through the EOS and transport coefficients, so satisfying them at t=0 does not imply preserving them in time. The frame-selection prescription is a real contribution, but saying that it works for 'a broad class of equations of state' is stronger than showing that the particular simulated trajectories remain admissible; the text even acknowledges restrictions. This makes the numerical results for QCD EOS and curved backgrounds conditional on an unmonitored assumption. The central mathematical theorems themselves are not undermined: if Assumption IV.1 holds, the proofs presented are detailed and internally consistent as far as the available text permits checking. I also considered the abstract's claim that both algorithms are causal despite superluminal cleaning speeds mentioned in Section III.D; that is a real overstatement, but it is a wording/scope issue that does not attack the LWP and equivalence core, whereas the IV.1 gap directly gates the applicability of those core results. The reader's CONDITIONAL verdict is therefore appropriate, and I see no basis to move it.","tokens_in":65349,"tokens_out":18452,"duration_ms":122201,"concrete_test":"Add a runtime diagnostic to the first-order and mixed-order codes that computes, at every cell and time step, c1², c+², c−², and Δ from Definition IV.2 using the actual EOS and frame parameters, and records the minima of (1−c_i²) and of |c_i²−c_j²| for i,j∈{1,+,−}. Run the TRENTo initialization at the largest shear viscosity reported and the Kerr-torus run to final time. If any c_i leaves (0,1), or any pairwise difference becomes zero, before final time, then that simulation lies outside the hypotheses of Theorems IX.1 and X.1 and the 'broad class of EOS' claim must be restricted accordingly. If the diagnostics remain uniformly positive throughout both runs, the concern is settled for those benchmarks.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theorems' conclusions — local well-posedness of the cleaning first-order system (74), local well-posedness of the mixed-order system (84), and strict equivalence to the original BDNK equations (55) — all require Assumption IV.1: positivity of τε, ρ, τQ, τP, Δ>0, c1,c±∈(0,1), and pairwise distinctness of c1,c+,c−. These are not identities satisfied by every EOS; they are state-dependent inequalities involving the equation of state and the transport coefficients. The paper offers a frame-selection prescription in Section XII.A, but the reviewed text does not prove that the inequalities are preserved by the evolution, nor does it report a diagnostic showing they remain satisfied in the TRENTo, QCD-EOS, Schwarzschild, or Kerr runs. Since ε and n evolve, c1,c± move; pairwise distinctness is an open condition that can fail at an isolated state, and any of the c_i can cross 1. If a trajectory leaves the admissible set, det M can degenerate or the principal symbol can lose the seven distinct real roots used in Step 4 of Section IX and in Theorem XI.2, so Theorems IX.1 and X.1, and the equivalence statements of Section VIII, no longer cover that solution. The paper itself qualifies the prescription as applying 'with certain restrictions' (Section I.A), so the unresolved question is whether the actual benchmarks stay inside those restrictions. This is not a defect in the conditional theorems, but it is the weakest point of the paper's broader claim of a practical admissible-frame prescription for nonconformal EOS relevant to heavy-ion collisions and astrophysics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This paper develops the mathematical foundations and numerical implementation of BDNK relativistic viscous hydrodynamics on fixed 3+1 background spacetimes, and presents two independent flux-conservative formulations: a first-order reduction with auxiliary gradient variables and div-curl cleaning fields (Section V, Equations (74)), and a cleaning-free mixed-order system that is first-order in time and second-order in space (Section VI, Equations (84)). Theorems IX.1 and X.1 establish local well-posedness in Sobolev spaces for the two systems under Assumptions IV.1-IV.2 together with real-analytic coefficients, suitable Sobolev data, and a compact-containment condition on the physical jet, and Section VIII proves strict equivalence of both formulations to the original second-order BDNK equations (55) for compatible initial data and vanishing cleaning fields. The proofs rest on a new general theorem (Section XI) giving sufficient conditions for reducing a mixed-order system with arbitrary-order equations to a strictly hyperbolic system, presented as a strengthening of Choquet-Bruhat's weak-hyperbolicity result. The numerical part validates the algorithms on smooth shocks, Kelvin-Helmholtz instabilities, Bjorken and Gubser flows, Michel accretion and Kerr tori, applies them to heavy-ion-like settings with a lattice-QCD equation of state and TRENTo initial conditions, and compares the results with a minimal Israel-Stewart theory.","tokens_in":65859,"tokens_out":18376,"duration_ms":112142,"significance":"The paper makes a substantial contribution if its claims hold. The well-posedness proofs are unusually detailed for a physics-oriented manuscript: explicit principal symbols (Equations (124), (129), (142)), characteristic polynomials, diagonalizability and root-distinctness checks, and reconstruction arguments are written out, and the hypotheses of the new reduction theorem (Section XI) are verified step by step in Sections IX and X. The equivalence proofs in Section VIII are constructive, and no fitted constants are promoted to predictions: the frame parameters are chosen to satisfy the causality inequalities, so circularity concerns do not arise. The paper is also notably honest about limitations: it states the one-derivative loss in the first-order existence proof (Sections III.B and IX), flags the nonphysical character and possibly superluminal speeds of the cleaning fields (Remark III.1), emphasizes that well-posedness is local, and qualifies the frame prescription as valid with certain restrictions.","major_comments":[{"comment":"The well-posedness theorems (Theorems IX.1 and X.1) and the equivalence statements of Section VIII are conditional on the physical jet remaining in the admissible set U_FO/U_MO, where Assumption IV.1 holds: positivity of tau_epsilon, rho, tau_Q, tau_P, Delta > 0, c_1, c_+, c_- in (0,1), and pairwise distinctness of c_1, c_+, c_-. These are state-dependent inequalities in (epsilon, n); as epsilon and n evolve, the characteristic speeds move, and the paper proves neither invariance of the admissible set under the flow nor that a simulated trajectory cannot cross a degeneracy Delta = 0 or a characteristic speed equal to 1 or to another speed. Section IV asserts that the causality conditions of [39, Theorem I] will be satisfied in the numerical simulations, and Section XII.A provides a frame-selection prescription explicitly qualified as applying with certain restrictions, but the reviewed text reports no diagnostic showing that Assumption IV.1 remained satisfied throughout the Bjorken, Gubser, Michel-accretion, Kerr-torus, or TRENTo runs, nor a check of the cleaning-speed bound of Section III.D on the curved-background runs. This is not a defect of the conditional theorems, but it is a gap in the paper's broader claim of practical admissibility: on any portion of a trajectory that leaves U_FO/U_MO, Theorems IX.1 and X.1 no longer cover the solution, and the equivalence to the original BDNK equations (55) established in Section VIII ceases to apply. I ask the authors to either (i) report, for each benchmark, the extrema and time histories of the quantities entering Assumption IV.1 (minima of tau_epsilon, rho, tau_Q, tau_P, Delta, and the ranges and pairwise separations of c_1, c_+, c_-) over the computational domain and integration interval, or (ii) explicitly delimit the maximal time interval on which the jet provably remains in the admissible set and restrict the claims about the simulations accordingly.","section":"IV (Assumption IV.1); IX-X; XII.A; XIII-XV","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The opening sentence of Section VIII reads \"In this Section we established that the new formulations...\", but the section establishes the statements; the tense should be corrected.","section":"VIII","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"In the displayed formula for Q^nu, the extracted text shows an isolated \"epsilon\" and a misplaced equation number at the end of the displayed expression; the typesetting of Equation (52) should be checked.","section":"IV, Eq. (52)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Remark IV.2 asserts without proof or a precise pointer that Assumption IV.1 is equivalent to the well-posedness hypotheses of [39, Theorem II]; since Assumption IV.1 carries all the hypotheses of the main theorems, an argument or an exact citation should be supplied.","section":"IV, Remark IV.2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Theorems IX.1 and X.1 require real-analytic coefficients and Sobolev data of order H^{ell+7}/H^{ell+8} with ell >= 4 or ell >= 3, whereas the numerical benchmarks include shocks, torus-atmosphere interfaces (Appendix E), and fluctuating TRENTo initial data of far lower regularity; a sentence stating explicitly that the regularity hypotheses are sufficient conditions for the continuum PDE results, and that numerical reliability is validated separately by the convergence tests of Section II.B, would prevent the theorems from being read as covering the discrete runs directly.","section":"IX, X vs. XIII-XV","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The sufficient bound on the cleaning speed c_h in Section III.D is derived in lapse-shift form for general metrics; reporting the actual values of kappa^A_1, kappa^A_2, lambda_1, lambda_2 used in the Milne, Schwarzschild, and Kerr runs, together with a check against this bound, would make the cleaning-parameter choices auditable.","section":"III.D, XII.C","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"The thing to know: this is a substantial paper, not another BDNK application. It builds the first general 3+1 curved-spacetime numerical framework for BDNK, gives two independent formulations, and backs both with real proofs. The mixed-order formulation is first-order in time and second-order in space, they prove it equivalent to the original second-order BDNK system, and the first-order formulation with div-curl cleaning is also proven locally well-posed, with a one-derivative loss that they state plainly. The new strictly hyperbolic reduction theorem in Section XI is a genuine strengthening of Choquet-Bruhat's weak-hyperbolicity result. I have not verified every algebra step, but the principal-symbol work in Section IX is explicit enough to be checked.\n\nWhat is good: the benchmark suite is broad—Bjorken, Gubser, shocks, Kelvin-Helmholtz, Michel accretion, Kerr tori, TRENTo comparisons with Israel-Stewart. The physical findings are mostly sanity checks rather than surprises, but the comparisons at larger viscosities are useful. The equivalence and constraint-propagation arguments are careful; they even point out that damping cleaning fields alone does not enforce the full gradient identities.\n\nThe soft spots are real but contained. The main one is Assumption IV.1. Every major theorem—well-posedness of both systems, equivalence to BDNK—assumes positivity of the frame parameters, positive discriminant, and causal, pairwise-distinct characteristic speeds. These are state-dependent inequalities. The paper gives a frame prescription for several equations of state but does not prove the inequalities are preserved by the evolution, and it does not report a diagnostic showing the TRENTo, QCD-EOS, Schwarzschild, or Kerr runs stay inside the admissible region. The authors do say the prescription works \"with certain restrictions,\" so this is not a hidden flaw in the conditional theorems. It is a gap between the theorems and the claim that the framework is practical for nonconformal equations of state. It should be fixed in revision, at minimum with a diagnostic, and ideally with a proof that the admissible set is forward invariant or a clear statement of what happens when it is left. The abstract also overstates causality for the cleaning-field system: cleaning speeds can be superluminal, which the text acknowledges is a permissible tradeoff for nonphysical fields.\n\nTwo smaller things: no code or data are provided, so the numerical claims are not independently reproducible from the paper alone; and the paper is long, with a large pedagogical Part I, which is fine but makes refereeing heavy.\n\nThe paper deserves a serious referee. It is aimed at numerical relativists, heavy-ion modelers, and PDE people working on relativistic viscous fluids. I would send it out, ask the referee to verify the symbol arithmetic in Sections IX and X, and require the admissibility diagnostic before accepting.","headline":"Serious, mostly sound framework paper for 3+1 BDNK numerics; the math is real and the main caveat is that the frame-admissibility condition underpins the theorems but is not proven to hold through the simulated runs.","tokens_in":66421,"tokens_out":2544,"duration_ms":18771},"pith_extraction":{"msc":["35L60","35L65","35Q75","76Y05"],"pacs":["47.75.+f","25.75.-q"],"model":"deepseek/deepseek-v4-flash","headline":"Two flux-conservative forms of BDNK viscous hydrodynamics are proven locally well-posed in general 3+1 spacetimes, one strictly equivalent to the original equations; both are validated from heavy-ion to black-hole settings.","keywords":["BDNK theory","relativistic viscous hydrodynamics","local well-posedness","strict hyperbolicity","constraint cleaning","mixed-order PDE systems","heavy-ion collisions","general-relativistic hydrodynamics"],"falsifier":"Take a nonconformal equation of state outside the paper's list, for instance a tabulated neutron-star-matter EOS, and a plausible merger or accretion trajectory; evaluate $\\Delta$ in definition (58c) and $c_1$, $c_+$, $c_-$ along it under the paper's frame prescription. If any speed leaves $(0,1)$, two speeds collide, or $\\Delta$ turns negative, then $\\det M$ vanishes or a characteristic root becomes complex, and Theorems IX.1 and X.1 plus the equivalence results cease to apply. A purely numerical check is to evolve the mixed-order code on heavy-ion initial conditions and compare directly with the second-order BDNK equations; disagreement while Assumption IV.1 holds signals a discretization bug, whereas disagreement coinciding with its violation locates the boundary of guaranteed validity.","tokens_in":69340,"feed_emoji":"🌊","tokens_out":5033,"duration_ms":106837,"temperature":0.7,"pith_summary":"BDNK is a first-order relativistic viscous hydrodynamic theory that remains causal and stable by redefining temperature, velocity, and chemical potential through the hydrodynamic frame. This paper's central claim is that BDNK can be written in two flux-conservative forms, both locally well-posed in Sobolev spaces on arbitrary 3+1 spacetimes, and one — the mixed-order form — is strictly equivalent to the original second-order equations, so simulations of it are unambiguously simulations of BDNK. The first-order form needs cleaning fields for its curl constraints, and the proof of well-posedness for that augmented system is included. The authors prove their results with a new theorem for reducing mixed-order systems to strictly hyperbolic form, and they validate both codes with convergence tests, Bjorken and Gubser flows, shock broadening, Kelvin–Helmholtz suppression, Michel accretion and Kerr tori, and heavy-ion comparisons with Israel–Stewart. A sympathetic reader would care because this provides the general-purpose numerical backbone needed to decide whether first-order viscous hydrodynamics describes quark-gluon plasma and neutron-star mergers.","feed_headline":"Two new schemes make BDNK hydrodynamics provably solvable","feed_subtitle":"Both are hyperbolic, causal, and proven equivalent to the original — ready for curved spacetimes and heavy-ion flows.","key_machinery":"The load-bearing object is the matrix $M^{\\mu\\beta}=H^{\\mu 00\\beta}$, the coefficient of $\\partial_t^2(\\varepsilon u^\\beta)$ in the conservation equations; Assumption IV.1 ensures $\\det M\\neq 0$, and the explicit inverse in Section V.A lets the time derivatives of the hydrodynamic fields be expressed algebraically through the conserved variables, converting a second-order system into flux-conservative form. The argument then rests on the principal symbol: for both formulations the characteristic determinant factors as $u_\\mu\\xi^\\mu$ times $(u_\\mu\\xi^\\mu)^2-c_1^2\\Delta_{\\mu\\nu}\\xi^\\mu\\xi^\\nu$, $(u_\\mu\\xi^\\mu)^2-c_+^2\\Delta_{\\mu\\nu}\\xi^\\mu\\xi^\\nu$, and $(u_\\mu\\xi^\\mu)^2-c_-^2\\Delta_{\\mu\\nu}\\xi^\\mu\\xi^\\nu$, so Assumption IV.1 yields seven distinct real characteristic speeds and strict hyperbolicity. The new mechanism is Theorem XI.1, a general result converting mixed-order systems with arbitrary equation orders into strictly hyperbolic systems, which supplies the Sobolev well-posedness for both formulations. For the first-order code, the div-curl cleaning fields — vector fields $\\varphi^{Ai}$ and $\\Psi^i$ with scalar partners $\\phi^A$ and $\\psi$ — control numerical violations of the curl constraints.","core_discovery":"On the paper's own terms, the result is that the BDNK system, which is second order in derivatives because the viscous terms enter through redefined $\\varepsilon$, $u^\\mu$, and $n$, admits two independent flux-conservative formulations whose Cauchy problem is locally well-posed. Theorem IX.1 gives this for the first-order system (74) together with div-curl cleaning fields, and the proof solves the cleaning fields first, yielding a one-derivative loss that the authors attribute to the method; Section VIII.A proves that compatible initial data force the constraints $X^{Ai}=\\partial_i U^A$ and $Z_i=\\partial_i n$ to propagate, so the system reproduces the original BDNK equations. Theorem X.1 gives local well-posedness for the mixed-order system (84) with unknowns $T^{0\\mu}$, $U^A$, $J^0$, and Section VIII.B proves strict equivalence: every mixed-order solution gives a solution of (55) and conversely, with first time derivatives fixed by algebraic identities involving the matrix $M^{\\mu\\beta}=H^{\\mu 00\\beta}$. The engine is Theorem XI.1, new sufficient conditions in terms of Leray indices and the characteristic determinant under which a mixed-order system reduces to a strictly hyperbolic system. Under Assumption IV.1, the principal symbol has seven distinct real roots, and the paper supplies frame-selection prescriptions for conformal, QCD, and ideal-gas equations of state.","pith_inferences":["The theorems guarantee well-posedness only while the solution stays within the region fixed by Assumption IV.1; nothing in the paper proves a trajectory cannot leave it, so monitoring $c_1$, $c_+$, $c_-$ during evolution is a natural safeguard.","The mixed-order formulation seems the safer default for production because it needs no cleaning-speed or damping parameters, whereas the first-order formulation requires choosing cleaning parameters within metric-dependent bounds.","If the two codes ever disagree, the equivalence theorem would place the difference in the discretization or cleaning rather than the underlying physics.","A testable extension is to apply the frame prescription to a tabulated neutron-star-matter equation of state and check Assumption IV.1 along merger trajectories; failure there would demarcate the theory's domain."],"forward_implications":["A general-purpose BDNK code can now run in any prescribed 3+1 spacetime without symmetry assumptions; the authors identify this as the main step toward fully dynamical BDNK–Einstein simulations in which viscosity may enter binary-neutron-star gravitational-wave signals.","Because both formulations are equal to the same original equations, agreement between the two codes becomes a built-in reliability check in regimes without analytic solutions.","The mixed-order system's strict equivalence means its nine-variable simulations carry no cleaning parameters: solving it is provably solving BDNK.","Theorem XI.1 applies beyond BDNK to any mixed-order system satisfying its Leray-index and symbol conditions, giving strict hyperbolicity and Sobolev well-posedness.","In heavy-ion simulations, BDNK and a minimal Israel–Stewart theory agree at small shear viscosity but differ substantially at larger viscosities, especially for early initialization with rapid longitudinal expansion."],"supporting_citations":[{"why":"Supplies the BDNK energy-momentum tensor, the constitutive relations, the causality and stability conditions behind Assumption IV.1, and the original first-order reformulation the new algorithms build on.","marker":"[39]"},{"why":"Cited jointly with [39] as the source of the first-order causal and stable formulation via hydrodynamic-frame freedom.","marker":"[70]"},{"why":"Earlier numerical implementations of BDNK that this work extends to arbitrary 3+1 metrics and to two independent schemes.","marker":"[72–74]"},{"why":"Supplies the GLM divergence-cleaning method that motivates the cleaning fields used in the first-order formulation.","marker":"[146]"},{"why":"Supplies the combined curl-plus-divergence cleaning system adopted in Equations (74f)–(74i) and (87).","marker":"[147]"},{"why":"The central scheme whose flux treatment naturally accommodates the mixed-order conservation form.","marker":"[156]"},{"why":"Cited for global existence of BDNK solutions for small data, the broader context for the local well-posedness theorems.","marker":"[126,127]"}],"fun_headline_variants":["Two new schemes make BDNK provably well-posed","BDNK hydrodynamics: two proof-backed solvers for curved spacetimes","First-order viscous hydro gets two proven flux-conservative schemes","BDNK numerical schemes proven hyperbolic, causal, and well-posed"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The whole edifice rests on Assumption IV.1: the characteristic speeds $c_1$, $c_+$, $c_-$ must stay real, distinct, and strictly between 0 and 1, with $\\tau_\\varepsilon$, $\\rho$, $\\tau_Q$, $\\tau_P$ positive, throughout the evolution; the paper's frame prescription realizes this for conformal, QCD, and ideal-gas equations of state, but does not prove the inequalities hold for every nonconformal equation of state or that a trajectory can never leave the admissible region.","fun_headline_variants_meta":{"raw":{"variants":["Two new schemes make BDNK provably well-posed","BDNK hydrodynamics: two proof-backed solvers for curved spacetimes","First-order viscous hydro gets two proven flux-conservative schemes","BDNK numerical schemes proven hyperbolic, causal, and well-posed"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.000131,"raw_usage":{"cost":0.0001312402,"is_byok":false,"cost_details":{"upstream_inference_cost":0.0001312402,"upstream_inference_prompt_cost":3.64182e-05,"upstream_inference_completions_cost":9.4822e-05},"total_tokens":1770,"prompt_tokens":1249,"completion_tokens":521,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":1061,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":457}},"tokens_in":1249,"tokens_out":521,"cost_source":"openrouter_usage","duration_ms":5395,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":457,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-09T05:33:44.482625+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"pipeline_version":"pith-pipeline@v1.4.1-daily-deepseek-or-medium","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"Take a nonconformal equation of state outside the paper's list, for instance a tabulated neutron-star-matter EOS, and a plausible merger or accretion trajectory; evaluate $\\Delta$ in definition (58c) and $c_1$, $c_+$, $c_-$ along it under the paper's frame prescription. If any speed leaves $(0,1)$, two speeds collide, or $\\Delta$ turns negative, then $\\det M$ vanishes or a characteristic root becomes complex, and Theorems IX.1 and X.1 plus the equivalence results cease to apply. A purely numerical check is to evolve the mixed-order code on heavy-ion initial conditions and compare directly with the second-order BDNK equations; disagreement while Assumption IV.1 holds signals a discretization bug, whereas disagreement coinciding with its violation locates the boundary of guaranteed validity.","supporting_citations":[{"cited_title":"First-Order General-Relativistic Viscous Fluid Dynamics","cited_arxiv_id":"2009.11388","evidence_quote":"Supplies the BDNK energy-momentum tensor, the constitutive relations, the causality and stability conditions behind Assumption IV.1, and the original first-order reformulation the new algorithms build on."},{"cited_title":"First-order relativistic hydrodynamics is stable","cited_arxiv_id":"1907.08191","evidence_quote":"Cited jointly with [39] as the source of the first-order causal and stable formulation via hydrodynamic-frame freedom."},{"cited_title":"D. Kr¨ oner, C.-D. Munz, T. Schnitzer, and M. Wesenberg, Journal of Computational Physics 175, 645","cited_arxiv_id":null,"evidence_quote":"Supplies the GLM divergence-cleaning method that motivates the cleaning fields used in the first-order formulation."},{"cited_title":"Journal of Computational Physics 404, 109088","cited_arxiv_id":null,"evidence_quote":"Supplies the combined curl-plus-divergence cleaning system adopted in Equations (74f)–(74i) and (87)."},{"cited_title":"Journal of Computational Physics 160, 241","cited_arxiv_id":null,"evidence_quote":"The central scheme whose flux treatment naturally accommodates the mixed-order conservation form."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"deepseek/deepseek-v4-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}