{"id":"c379b24c-d745-45f3-9c46-605959905a34","arxiv_id":"2508.21427","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive an entropy-conservative two-point flux for the ultra-relativistic Euler equations, prove its consistency, and validate an entropy-stable DG scheme against 1D radial reference solutions in 2D and 3D.","lead":"This paper derives the first entropy-conservative numerical flux for the ultra-relativistic Euler equations and builds a discontinuous Galerkin scheme around it, testing it in two and three dimensions. The simulations match an established radially symmetric reference solver on shock, rarefaction, and bubble tests, adding genuine 3D benchmark results.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The entropy-conservative flux derivation is the load-bearing step; it appears correct, but the printed algebra is compressed and no formal or CAS check is provided, and the method's p>0 restriction is unanalyzed for the low-pressure 3D pre-focus regime.","rationale":"The reader's verdict is ACCEPT with high confidence. I agree with the reader that the positivity/vacuum limitation is real, and I agree that the numerical comparisons are broadly consistent with the disclosed 3D under-resolution. However, I see the most load-bearing concern as the correctness of the algebraic derivation of Theorem 3.1 itself: the paper's central contribution is the new entropy-conservative flux, and the printed derivation is compressed in exactly the two places where the flux formulas (3.13) and (3.15) are selected. The entropy-rate test in Fig. 3 is supportive but would also pass for a slightly wrong flux that happened to conserve entropy on that particular smooth self-similar solution. The positivity gap is a genuine limitation but is explicitly acknowledged by the restriction p>0 and by the disclosure that the d=3 vacuum expansion is not considered; it does not undermine the stated theorem. A CAS-based residual check of (3.1) plus near-vacuum sampling would settle whether the derivation is fully correct and would also quantify the practical failure boundary near p=0. Therefore the verdict should remain UNCHANGED.","tokens_in":30930,"tokens_out":1617,"duration_ms":20465,"concrete_test":"Independently verify the Tadmor condition (3.1) for the flux (3.16) using a computer algebra system: sample random positive states (p_+, u_+, p_-, u_-) in d=2 and d=3 and check that the residual of (3.1) is zero to machine precision. Also run the same check for near-vacuum states, e.g. p_+=1e-3, p_-=1e-5, |u| up to 10, to confirm the algebraic identity still holds and to identify when roundoff degrades the residual.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is Theorem 3.1: the two-point flux (3.16) satisfies the Tadmor entropy condition (3.1). The printed derivation jumps from (3.12) to (3.13) and from (3.14) to (3.15) with the statement that each parenthesized term vanishes; I did not find a sign or indexing error, but the algebra is compressed and the equal-state consistency is asserted rather than demonstrated in the text. The entropy-rate tests in Fig. 3 are the only direct numerical evidence of the entropy-conservation property, and they only cover a smooth self-similar case; a small inconsistent flux could still yield near-machine-zero rates on that test. Additionally, the flux and the entropy variables (2.12) are only defined for p>0, and Examples 3--5 in 3D reach very low pressures before the focusing blow-up (e.g., Example 3 states the shock reaches the origin around t=4.165 and the pressure takes very low values). The paper applies the positivity-preserving limiter only for Example 2, not for Examples 3--5, and offers no analysis of whether the semidiscrete scheme preserves p>0. This is a real gap in the robustness claim, but it does not refute the algebraic EC-flux derivation itself.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an entropy-conservative two-point numerical flux for the multi-dimensional ultra-relativistic Euler equations, following the Tadmor framework. The authors compute the entropy variables and entropy potential, solve the Tadmor jump condition, and obtain explicit flux formulas (3.16). They embed this flux in a high-order DGSEM with flux differencing, Rusanov interface dissipation, finite-volume subcell shock capturing, AMR, and adaptive time stepping in Trixi.jl. They validate the method on five radially symmetric benchmark problems in 2D and 3D from Kunik et al., comparing against the dedicated 1D solver RadSymS and self-similar ODE reference solutions, and they report machine-precision entropy rates in one test.","tokens_in":31253,"tokens_out":23475,"duration_ms":225564,"significance":"The main claimed contribution is the first entropy-conservative two-point flux for the ultra-relativistic Euler equations, and the paper demonstrates its use in genuine 2D and 3D simulations of the recently proposed benchmarks. If the derivation is correct, this is a substantive, reusable building block for entropy-stable DG methods for relativistic flows. The paper is also strong on reproducibility: it provides a public Julia repository with code and data, and the numerical benchmarks are independently cross-checked against RadSymS and ODE-based self-similar solutions. The central algebraic derivation in Section 3 is checkable and not fitted, and the machine-precision entropy test, once precisely defined, is a meaningful falsifiable check.","major_comments":[{"comment":"Eq. (2.13) prints the entropy potential as ψ_k = -1/4 S u_k sqrt(1+|u|^2). With S=p^{3/4}sqrt(1+|u|^2), this equals -1/4 p^{3/4} u_k (1+|u|^2), whose jump is not the term [p^{3/4}u_k] used in (3.11). The derivation is consistent with ψ_k = -1/4 S u_k / sqrt(1+|u|^2) = -1/4 p^{3/4}u_k. Please correct (2.13) and adjust the sign/presentation in (3.1); as printed, the derivation cannot be followed.","section":"Eq. (2.13) and (3.11)"},{"comment":"The entropy variables (2.12) and all flux formulas (3.6)-(3.8) require p>0, and this domain of validity is never stated. The positivity-preserving limiter is applied only to Example 2 (Section 4), while Examples 3-5 reach very low pressures before focusing (e.g., Example 3 says p takes very low values around t=4.165 in 3D). No analysis or numerical check is provided that p remains positive for the computed solutions, so the scheme could, in principle, leave the domain of definition of the flux. State the p>0 assumption explicitly, verify it for the reported runs, or extend positivity protection/fallbacks to all tests.","section":"Section 3, eqs. (3.6)-(3.8); Section 5, Examples 3-5"},{"comment":"The proof of convexity of -S is incomplete. The text verifies that A1 is positive on vectors (v1,...,vd,0), A2 nonnegative on the same subspace, and A3 positive on (0,...,0,v_{d+1}), but this does not imply that A1+A2+A3 is positive definite on mixed vectors. The Hessian must be tested on arbitrary v=(y,z), including cross terms in A3. Either complete the argument or cite an existing proof of concavity of (2.9); the entropy-stable framework depends on this property.","section":"Appendix A, eqs. (A.6)-(A.10)"},{"comment":"The full scheme includes Rusanov interface dissipation and finite-volume subcell shock capturing, both entropy-dissipative. The statement that 'entropy is conserved up to machine precision' is therefore ambiguous for a test based on Example 1, which is a shock problem. Please define precisely what is plotted (e.g., the semidiscrete entropy residual of the volume EC flux only, or the total entropy balance), and reconcile any claim of exact conservation with the expected entropy production at shocks.","section":"Section 5, Fig. 3"}],"minor_comments":[{"comment":"In the displayed formula for the jump of the entropy variable, the expression for [p^{-1/2}] appears to contain a duplicated factor [p][p]; this typesetting issue should be cleaned up.","section":"Eq. (3.10)"},{"comment":"The statements 'This is clearly a consistent approximation' and 'This is a consistent approximation' are assertions; a one-line equal-state verification (p_L=p_R, u_L=u_R) would make the proof self-contained.","section":"Eqs. (3.13) and (3.15)"},{"comment":"Reference [69] has DOI 'TODO' and reference [13] has a truncated DOI. Please update before publication.","section":"References"},{"comment":"The caption should state how the 'entropy rate' is computed in space and time, and in which norm or discrete functional it is measured.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The central EC-flux derivation appears mathematically sound, but the manuscript needs a corrected Eq. (2.13), a completed positivity/domain discussion, a fixed convexity proof, and a precise definition of the reported entropy rate. These are local but load-bearing issues; they should be addressable in revision. The typo in (2.13) is likely an OCR/LaTeX artifact, but as printed it breaks the derivation chain."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the new EC flux in Theorem 3.1 is the real thing. I checked the equal-state limit for (3.16) and it gives the physical flux, and the Tadmor condition is satisfied by construction. The derivation is a bit compressed—the jumps from (3.12) to (3.13) and (3.14) to (3.15) skip some bookkeeping—but I did not find an error. The machine-precision entropy rates in Fig. 3 are a useful independent check, and the 2D agreement with the RadSymS reference and the ODE self-similar solutions is convincing. The 3D results are honestly presented as first attempts: the pressure peaks are under-resolved by factors of 10–30 and the authors say so explicitly. That is a limitation, not a defect.\n\nThe main soft spot is the p>0 restriction. The flux, the entropy variables, and the jump identities all require strictly positive pressure on both interface states. The paper only applies the Zhang–Shu positivity limiter in Example 2, and it offers no argument that the semidiscrete scheme preserves positivity for Examples 3–5, where the pressure gets very low before focusing. If an interface state ever hits p<=0, the flux formula blows up. This does not refute the EC flux derivation itself—that is an algebraic statement about positive states—but it does mean the method's robustness in near-vacuum regimes is unproven. A short discussion or a numerical experiment probing how close p gets to zero would have addressed it.\n\nThe citation pattern looks fine. The flux is genuinely new, not fitted; the entropy variables come from prior work on the same system; and the reproducibility repository is a plus, though the DOI in the references is still a TODO placeholder. The paper does what it claims: it derives the first EC flux for this system, demonstrates entropy conservation to machine precision, and delivers first 3D benchmark runs with clearly disclosed resolution limits. I would send this to a serious referee, and I would expect it to be accepted after minor revisions that acknowledge the positivity gap more prominently.","headline":"A genuinely new entropy-conservative flux for the ultra-relativistic Euler equations, derived cleanly and tested honestly; the positivity gap is real but does not sink the main claim.","tokens_in":31732,"tokens_out":2564,"would_cite":true,"duration_ms":27172,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M06","65M20","65M70"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives the first entropy-conservative numerical flux for the ultra-relativistic Euler equations and shows in two- and three-dimensional runs that a discontinuous Galerkin scheme built on it conserves the physical entropy to mach","keywords":["ultra-relativistic Euler equations","entropy-stable methods","entropy-conservative numerical flux","discontinuous Galerkin methods","flux differencing","summation-by-parts operators","radially symmetric solutions","adaptive mesh refinement"],"falsifier":"Two concrete checks decide the claims. (1) The theorem is algebraic: evaluate ⟦v⟧·F̃ − ⟦ψ⟧ for any two states with p_± > 0 (say p_+ = 2, p_− = 1, u_+ = (0.3, 0.1), u_− = (−0.2, 0.4)); if the result is not identically zero, the flux is not entropy-conservative. (2) The positivity boundary: run the self-similar expansion of Example 2 in 3D with initial radial velocity close enough to 1 that the flow expands into vacuum — the case the paper excludes; as p → 0 at the wave front, the flux's √p and p^{1/4} denominators should produce NaNs or unphysical states, confirming the p > 0 restriction is rea","tokens_in":30844,"feed_emoji":"⚛️","tokens_out":20677,"duration_ms":170160,"temperature":0.7,"pith_summary":"The paper sets out to give the ultra-relativistic Euler equations — the regime of relativistic gas dynamics where thermal energy dominates — a numerical flux that cannot create entropy by construction. It derives the entropy variables (the 'main field') and the entropy flux potential, then solves the entropy-conservation condition (3.1) algebraically to obtain an explicit two-point flux, stated as Theorem 3.1. This is the first entropy-conservative flux for this system; a discontinuous Galerkin scheme built on it via flux differencing, subcell shock capturing, and adaptive mesh refinement conserves the physical entropy to machine precision in both two and three dimensions. Applied to the five radially symmetric benchmark problems proposed in recent companion work, the scheme reproduces the reference solutions — straight shock formation, self-similar rarefactions, and pressure blow-up at the origin — and delivers the first genuine three-dimensional results for these benchmarks. The point of the construction is structural: the semidiscrete spatial discretization conserves entropy exactly, so entropy dissipation enters only where the user adds it, and other multi-dimensional solvers now have a structure-preserving benchmark to test against.","feed_headline":"First entropy-conserving flux for ultra-relativistic Euler","feed_subtitle":"New two-point flux keeps entropy at machine round-off in 2D and 3D shock and blow-up runs.","key_machinery":"The engine of the derivation is the entropy-variable transformation (the 'main field') v = ∇_w S (2.12), the entropy flux potential ψ = −¼ p^{3/4} u√(1+|u|²) (2.13), and the entropy-conservation condition ⟦v⟧·F̃ − ⟦ψ⟧ = 0 (3.1) that a two-point flux must satisfy exactly. The key algebraic step expresses every jump in that condition through two basic jumps, ⟦u_i p^{−1/4}⟧ and ⟦p⟧, using mean-value identities for products, squares, and roots ((3.4)–(3.8)); the condition becomes a linear system whose brackets can be forced to vanish independently, yielding the explicit flux (3.16) — 'affordable' because it uses differential means rather than integral averages. Flux differencing with summation-b","core_discovery":"Central claim (Theorem 3.1): there is an explicit two-point flux — momentum part 2(p−√p+ + p+√p−){{u_i p^{−1/4}}}{{u_j p^{−1/4}}} + {{p}}δ_ij, energy part 2(p−√p+ + p+√p−){{p^{−1/4}√(1+|u|²)}}{{u_j p^{−1/4}}}, with {{·}} the arithmetic mean across the interface — that is symmetric, consistent, and entropy-conservative for the ultra-relativistic Euler equations. It satisfies the entropy-conservation condition exactly, so a semidiscrete DG discretization conserves the physical entropy S = p^{3/4}√(1+|u|²) to machine precision in the runs. This is the first entropy-conservative flux for this system.","pith_inferences":["The same differential-averaging construction should transfer to neighbouring relativistic systems whose main field is known — the Synge-gas and general (non-ultra) relativistic Euler equations — since the only structural ingredients are the entropy variables, the flux potential, and two convenient jump variables.","The p > 0 restriction marks the method's real boundary: a 3D vacuum-expansion run (Example 2 with initial velocity near light speed) would immediately break the flux. A natural extension is a blended flux that switches to a positivity-preserving or kinetic solver when an interface pressure approaches zero.","The under-resolved 3D peaks are a resolution statement, not a sign of instability: re-running the focusing examples with increased refinement near the origin should push the computed peak toward the reference values (≈289, ≈3550, ≈5957), and an entropy-stable baseline should keep those runs stable throughout.","Machine-precision entropy conservation in smooth regions suggests the scheme can serve as a quantitative probe of the pressure blow-up: tracking peak growth on adaptively refined meshes could test whether the focusing singularity is self-similar, as the one-dimensional analysis predicts."],"forward_implications":["A semidiscrete DG discretization using the new flux conserves the physical entropy to machine precision; measured entropy rates in Fig. 3 sit at round-off in both 2D and 3D, so entropy production is controlled by the explicit interface dissipation and shock capturing rather than by accidental numerical error.","The scheme reproduces the radially symmetric reference solutions for the shock-plus-stationary example and the self-similar rarefaction example, in both two and three dimensions, matching ODE-based reference solutions.","The benchmark suite is run for the first time as genuine 3D computations: the bubble expansion, bubble collapse, and sine-velocity cases all produce the expected shock focusing and pressure blow-up near the origin, with the 3D peak pressures under-resolved relative to the one-dimensional reference (e.g. ≈9.7 versus ≈289 in the expansion example) because of computational cost.","Because the flux is an explicit, symmetric two-point flux, any conservative scheme that uses it as the volume flux inherits semidiscrete entropy conservation; extending it to higher order requires only flux differencing with summation-by-parts operators.","The solved benchmarks, with reference solutions and reproducing code, give other multi-dimensional relativistic solvers a structure-preserving test case."],"supporting_citations":[{"why":"Supplies the entropy-conservation condition (3.1) that defines an entropy-conservative numerical flux.","marker":"[68]"},{"why":"Supplies the differential-averaging algorithm used to turn the entropy-conservation condition into an explicit affordable flux.","marker":"[49]"},{"why":"Provides the model kinetic-energy-preserving and entropy-stable flux whose affordable two-point construction is adapted here.","marker":"[9]"},{"why":"Establishes the 'affordable' entropy-consistent flux paradigm based on differential means.","marker":"[29]"},{"why":"Defines the physical entropy (2.9) and the Rankine-Hugoniot shock conditions used for the Example 1 reference solution.","marker":"[35]"},{"why":"Introduces the five multi-dimensional benchmark problems and the radially symmetric quasi-1D formulation used for reference solutions.","marker":"[36]"},{"why":"Develops the radially symmetric solutions and the one-dimensional reference solver that the comparisons rely on.","marker":"[37]"},{"why":"Supplies the self-similar ODE system and jump conditions used to construct exact reference solutions in Examples 1 and 2.","marker":"[41]"},{"why":"Provides the finite-volume subcell shock-capturing technique that stabilizes the DG scheme at discontinuities.","marker":"[24]"},{"why":"Proves the positivity-preserving limiter used to keep pressure positive in Example 2.","marker":"[74]"}],"fun_headline_variants":["First entropy-conservative flux for ultra-relativistic Euler","Entropy preserved to machine round-off in Euler DG tests","Exact entropy conservation in ultra-relativistic Euler DG","New two-point flux conserves entropy in relativistic shocks","Entropy-conserving flux derived for ultra-relativistic Euler"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"Every flux formula divides by √p and p^{1/4}, and the entropy variables are defined only where pressure is strictly positive, so the method has no definition — and no safeguard — at vacuum or near-vacuum interfaces; the test cases chosen deliberately avoid the vacuum expansion that the 3D data could produce.","fun_headline_variants_meta":{"raw":{"variants":["First entropy-conservative flux for ultra-relativistic Euler","Entropy preserved to machine round-off in Euler DG tests","Exact entropy conservation in ultra-relativistic Euler DG","New two-point flux conserves entropy in relativistic shocks","Entropy-conserving flux derived for ultra-relativistic Euler"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001277,"raw_usage":{"total_tokens":5072,"prompt_tokens":771,"completion_tokens":4301,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":4218}},"tokens_in":515,"tokens_out":4301,"duration_ms":32877,"temperature":1.0,"reasoning_tokens":4218,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T14:18:54.355198+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Two concrete checks decide the claims. (1) The theorem is algebraic: evaluate ⟦v⟧·F̃ − ⟦ψ⟧ for any two states with p_± > 0 (say p_+ = 2, p_− = 1, u_+ = (0.3, 0.1), u_− = (−0.2, 0.4)); if the result is not identically zero, the flux is not entropy-conservative. (2) The positivity boundary: run the self-similar expansion of Example 2 in 3D with initial radial velocity close enough to 1 that the flow expands into vacuum — the case the paper excludes; as p → 0 at the wave front, the flux's √p and p^{1/4} denominators should produce NaNs or unphysical states, confirming the p > 0 restriction is rea","supporting_citations":[{"cited_title":"The numerical viscosity of entropy stable s chemes for systems of conservation laws. I","cited_arxiv_id":null,"evidence_quote":"Supplies the entropy-conservation condition (3.1) that defines an entropy-conservative numerical flux."},{"cited_title":"Comparison of some Entropy Conservative Numerical Fluxes for the Euler Equations","cited_arxiv_id":"1701.02264","evidence_quote":"Supplies the differential-averaging algorithm used to turn the entropy-conservation condition into an explicit affordable flux."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the physical entropy (2.9) and the Rankine-Hugoniot shock conditions used for the Example 1 reference solution."},{"cited_title":"Radially symme tric solutions of the ultra- relativistic Euler equations in several space dimensions","cited_arxiv_id":null,"evidence_quote":"Introduces the five multi-dimensional benchmark problems and the radially symmetric quasi-1D formulation used for reference solutions."},{"cited_title":"Radially symmetric so lutions of the ultra-relativistic Euler equations","cited_arxiv_id":null,"evidence_quote":"Develops the radially symmetric solutions and the one-dimensional reference solver that the comparisons rely on."},{"cited_title":"Self-similar solutions of the radially symmet ric relativistic Euler equations","cited_arxiv_id":null,"evidence_quote":"Supplies the self-similar ODE system and jump conditions used to construct exact reference solutions in Examples 1 and 2."},{"cited_title":"A provably entropy stable subcell shock capturing approach for high or der split form DG for the com- pressible Euler equations","cited_arxiv_id":null,"evidence_quote":"Provides the finite-volume subcell shock-capturing technique that stabilizes the DG scheme at discontinuities."}],"review_version":1}