REVIEW 4 major objections 5 minor 29 references
Entropy-stable fluxes for high-order Discontinuous Galerkin simulations of high-enthalpy flows
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read This paper extends entropy-conserving high-order DG fluxes to multi-species reacting high-enthalpy flows by tabulating internal energies.
desk verdict A genuine multi-species extension of the authors' tabulated entropy-conservative flux with solid validation, but the presentation has gaps that need fixing before the method is fully reproducible. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the closed-form entropy integral obtained from piecewise-linear tabulated data. For each species, $\varepsilon^{\mathrm{int}}_c(T)$ and $c_{v,c}(T)$ are tabulated on a temperature grid of step $\Delta T$; linear reconstruction turns the entropy integral $\eta_c(T)=\int_0^T c_{v,c}(\tau)/\tau\,d\tau$ into Eq. (15), a sum of tabulated interval integrals $\eta_{c,i}$ plus one logarithm. This $\eta_c(T)$ feeds the entropy-conservative flux: logarithmic averages for species densities, geometric and arithmetic averages for temperature and internal-energy terms in the energy flux (Eq. (12)), so the discrete interface flux is consistent with the entropy pair and conserves kinetic energy. Precomputed prefix sums of $\eta_{c,i}$ keep the evaluation cost at one lookup and one logarithm per species.
What would settle it
Run the weak blast wave with the same physical O2/O mixture tabulated at a coarse step, say $\Delta T = 500$ K, and at a fine step, $\Delta T = 1$ K, on identical grids and time steps; if the coarse-table pressure and temperature fields differ by more than discretization error, or if entropy production leaves machine precision when temperatures approach the table boundary, the method's guarantee is tied to the interpolation table rather than to the physical model.
Extended reading notes
Core claim
The central claim is that a single numerical flux can conserve entropy and kinetic energy for a compressible gas mixture with arbitrary species internal-energy models. The construction replaces the true functions $c_{v,c}(T)$ and $\varepsilon^{\mathrm{int}}_c(T)$ by piecewise-linear tabulated reconstructions; the entropy integral $\eta_c(T)=\int_0^T c_{v,c}(\tau)/\tau\,d\tau$ then has the closed form of Eq. (15), so the entropy-conservative interface flux (Eqs. (9)-(12)) can be evaluated with table lookups and one logarithm per species. The paper verifies the flux on a weak blast wave, where discrete entropy production stays at machine precision, on homogeneous chemical relaxation, and on a Mach 10 reacting cylinder flow, matching an established finite-volume solver.
Load-bearing premise
The entropy-conservation guarantee holds for the piecewise-linear tabulated thermodynamic model, not for the true physical internal-energy functions, so a coarse temperature table would conserve the wrong gas's entropy.
Editorial extensions
If this is right
- Any species internal-energy model can enter the same flux by tabulation, with no analytic entropy integrals required.
- The cost per interface state is one table lookup and one logarithm per species, independent of thermodynamic model complexity.
- The discrete entropy production stays at machine precision in the non-reacting verification for polynomial orders 2, 4, and 6.
- Reacting high-enthalpy flows with strong shocks and adaptively refined non-conforming meshes run stably and match the benchmark solver.
- Chemical relaxation of an O2/O mixture reproduces the master-equation reference solution.
Reading between the lines
- A natural next step is to derive entropy-dissipative viscous and heat-conduction fluxes in the same tabulated framework, since the paper's conservation guarantee covers the inviscid Euler fluxes only.
- Table resolution $\Delta T$ becomes a user-facing accuracy parameter; for a given physical gas, practitioners will need a resolution criterion tied to the temperature range of the simulation.
- The same flux could ingest internal-energy tables from kinetic-theory or state-resolved chemistry data without a new derivation, coupling CFD directly to high-fidelity thermochemistry.
- In reacting simulations, the chemical source term is not entropy-conserving, so the overall scheme's stability depends on source-term treatment as well as on the flux.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper extends a previously proposed entropy-conservative flux construction for high-enthalpy single-species gas flows to multi-species reacting mixtures within a Discontinuous Galerkin Spectral Element Method (DGSEM). The flux uses piecewise-linear tabulated reconstructions of species specific heats and internal energies, with closed-form entropy integrals computed from the tabulated data, so that the cost of evaluating the flux is independent of the complexity of the underlying internal-energy expressions. The method is implemented in Trixi.jl and tested on three problems: a non-reacting weak blast wave (for which machine-precision entropy conservation is reported), a spatially homogeneous chemical relaxation of an O2/O mixture (compared against a master-equation reference), and a Mach 10 flow over a cylinder (compared against the DLR TAU solver). The paper claims excellent agreement in all cases and argues that the approach enables entropy-stable high-order simulations of high-enthalpy reacting flows with arbitrary internal energy models.
Significance. If the claims hold, this is a practically useful contribution: it provides a route to entropy-stable DGSEM for multi-species high-enthalpy flows with complex real-gas thermodynamics while keeping the per-element flux evaluation cheap and independent of the analytical form of the internal energy functions. The paper reports machine-precision entropy production in the blast-wave test, verifies the chemistry implementation against an independent master-equation solution, and validates the full reacting-flow setup against the established DLR TAU solver. The availability of a public reproducibility repository is a concrete strength. However, the manuscript is not fully self-contained: the flux derivation is delegated to a companion preprint, and two of the central printed formulas (Eqs. (14) and (15)) are corrupted or incorrect as written. These issues must be resolved before the central claim can be verified from the paper alone.
major comments (4)
- [Sec. 3, Eq. (14)] The printed definitions of c_{v,c}(T*) and c_{v,c}(T**) are garbled and cannot be read as mathematical formulas. In particular, the expression given for c_{v,c}(T**) contains extraneous characters and appears to mix an arithmetic mean with a jump in an undefined way. Since Eq. (12) is the central entropy-conservative energy flux, these definitions must be stated cleanly; otherwise the flux cannot be implemented from the manuscript and the claimed property cannot be checked.
- [Sec. 3, Eq. (15)] The closed-form expression for the reconstructed entropy integral η_c(T) is incorrect as printed. For a piecewise-linear c_{v,c}(T) on the interval [T_N, T_{N+1}], the coefficient multiplying ln(T/T_N) must involve the fixed interval slope (c_{v,c,N+1} - c_{v,c,N})/ΔT, but Eq. (15) instead contains (c_{v,c}(T) - c_{v,c,N}), which is not constant on the interval. As written, η_c(T) is not the exact integral of the linearly interpolated specific heat. This matters because Eq. (15) feeds directly into the flux through Eq. (14) and is the basis for the machine-precision entropy conservation reported in Sec. 4.1. The authors should correct the formula or clarify the intended notation.
- [Sec. 3, Eqs. (14)-(16) and Sec. 4.1] The paper does not state a consistency condition between the tabulated c_{v,c}(T) and ε_int,c(T) reconstructions. For the discrete entropy identity to hold, the temperature recovered from the total energy E must be compatible with the entropy variables computed from η_c; in particular, the reconstructed functions should satisfy dε_int,c/dT = c_{v,c}(T). When c_v and ε_int are linearly interpolated independently on the same ΔT grid, this equality fails away from the limit ΔT→0. The numerical entropy-production result in Sec. 4.1 demonstrates the property for the specific implementation, but it does not establish the general claim in the abstract for arbitrary tabulated data. The authors should either impose a consistency constraint on table generation (e.g., deriving ε_int from the same piecewise-linear c_v) or provide an analysis of how ΔT degrades entropy conservation and solution accuracy.
- [Sec. 3, flux derivation] The multi-species flux in Eqs. (9)-(12) is presented without derivation; the text only states that it follows the procedure of the authors' prior single-species work [17]. Since the central contribution of this manuscript is precisely this multi-species entropy-conservative flux, the paper should include the key steps of the Tadmor entropy-conservation condition for the reacting multi-species system, or at minimum state explicitly which equations and assumptions are being carried over from [17]. As it stands, a reader cannot verify the entropy and kinetic-energy preservation properties without consulting a separate preprint.
minor comments (5)
- [Sec. 1] There is a typo in the Introduction: "next-generation CDF solvers" should read "next-generation CFD solvers."
- [Sec. 2, Eq. (4)] The two cases for atoms and molecules are presented in a single run-on line; a proper case distinction would make the formula readable and avoid ambiguity about which expression applies to which species type.
- [Sec. 2, Eq. (8)] The physical entropy per unit mass and the mathematical entropy density are both denoted by the symbol s (with the mathematical entropy defined as −ρs). Using two different symbols (e.g., s and S̃) would prevent confusion in the subsequent derivation.
- [Sec. 4.1, Fig. 2] The y-axis label "|s|" is undefined; the plotted quantity appears to be the absolute value of the entropy production rate, so the label should read something like "|ds/dt|" or "entropy production rate."
- [Sec. 4.1] The paper claims the flux is kinetic-energy-preserving, but none of the numerical tests directly reports a kinetic-energy conservation error. Adding such a measurement would strengthen the empirical support for that part of the claim.
Circularity Check
No material circularity: the multi-species flux is constructed from the authors' prior self-cited derivation, but the entropy-conservation check is internal consistency and the validation rests on external master-equation and DLR TAU references.
full rationale
The paper's central contribution is a construction: given the entropy definition (8) and the piecewise-linear tabulated thermodynamic representation (14)-(16), the flux (9)-(12) is derived (via the authors' prior work [17]) so that the discrete entropy balance vanishes identically. The measured machine-precision entropy production in Fig. 2 therefore confirms implementation of this self-consistent construction rather than serving as an externally fitted prediction; this is not a circular reduction of a claim to a fitted input. No parameter of the method is tuned to reproduce the benchmarks: Arrhenius coefficients come from [24], initial conditions are stated independently, and the reference solutions are an independent master-equation solver and the DLR TAU code. The only self-citation is [17], which supplies the flux-derivation procedure; because the formulas are restated in the present paper and the resulting scheme is validated against external references, this self-citation is not load-bearing in the sense of an unverified premise that forces the conclusions. A genuine caveat is that entropy conservation is established for the tabulated interpolated cv,c and ε_int,c model, not for the original physical internal-energy functions, and Eq. (15) as printed contains an apparent interval-slope inconsistency; these are scope and correctness concerns, not circularity.
Assumptions & free parameters
free parameters (1)
- Temperature tabulation step ΔT =
not specified in paper
assumptions (6)
- domain assumption Ideal gas closure p = nkT with n = Σ_c ρ_c/m_c
- domain assumption Entropy definition (Eq. 8): s = Σ_c Y_c ∫_0^T cv,c/τ dτ - Σ_c Y_c (k/m_c) ln ρ_c
- domain assumption Internal-energy decomposition into translational, rotational, vibrational, and formation energies (Eq. 4)
- domain assumption Piecewise-linear table interpolation of cv,c and ε_int,c is an exact representation of the gas for the purposes of entropy computation (Eqs. 15-16)
- ad hoc to paper The multi-species flux expressions (9)-(12) are obtained by following the procedure of the authors' prior single-species derivation [17]
- domain assumption Arrhenius dissociation chemistry with coefficients from Park [24]
Cite this review
Pith. "Pith review of Entropy-stable fluxes for high-order Discontinuous Galerkin simulations of high-enthalpy flows." pith.science (2026). https://pith.science/paper/6ARJWLFB
@misc{pith2026241113168,
author = {Pith},
title = {Pith review of: Entropy-stable fluxes for high-order Discontinuous Galerkin simulations of high-enthalpy flows},
year = {2026},
howpublished = {\url{https://pith.science/paper/6ARJWLFB}},
note = {Machine review of arXiv:2411.13168}
}
read the original abstract
In the present work, we extend the Discontinuous Galerkin Spectral Element Method (DGSEM) to high-enthalpy reacting gas flows with internal degrees of freedom. An entropy- and kinetic energy-preserving flux function is proposed which allows for use of arbitrary expressions for the internal energies of the constituent gas species. The developed method is applied to simulation of several model problems and compared to the DLR TAU solver.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[17]
Entropy-conservative high-order methods for high-enthalpy gas flows
G. Oblapenko and M. Torrilhon, “Entropy-conservative high-order methods for high-enthalpy gas flows,” arXiv:2403.16882, 2024
work page Pith review arXiv 2024
-
[1]
S. Karl and T. Bykerk, “Sustainable space technologies—strategies toward a predictive aerothermal design of re-useable space transportation systems,” Rev. Sci. Instrum., vol. 95, no. 2, 2024
work page 2024
-
[2]
High-order CFD methods: current status and perspective,
Z. Wang, K. Fidkowski, R. Abgrall, F. Bassi, D. Caraeni, A. Cary, H. Deconinck, R. Hartmann, K. Hillewaert, H. Huynh, N. Kroll, G. May, P.-O. Persson, B. van Leer, and M. Visbal, “High-order CFD methods: current status and perspective,” Int. J. Numer. Meth. Fluids, vol. 72, pp. 811–845, July 2013
work page 2013
-
[3]
CFD vision 2030 study: a path to revolutionary computational aerosciences,
J. P. Slotnick, A. Khodadoust, J. Alonso, D. Darmofal, W. Gropp, E. Lurie, and D. J. Mavriplis, “CFD vision 2030 study: a path to revolutionary computational aerosciences,” Tech. Rep. CR–2014-218178, NASA, 2014
work page 2014
-
[4]
A brief survey on discontinuous Galerkin methods in computational fluid dynamics,
C.-W. Shu, “A brief survey on discontinuous Galerkin methods in computational fluid dynamics,”Advances in mechanics, vol. 43, no. 6, pp. 541–553, 2013. 6 Entropy-stable fluxes for high-order Discontinuous Galerkin simulations of high- enthalpy flows
work page 2013
-
[5]
G. J. Gassner and A. R. Winters, “A novel robust strategy for discontinuous Galerkin methods in computational fluid mechanics: Why? When? What? Where?,” Front. Phys., vol. 8, p. 500690, 2021
work page 2021
-
[6]
E. Tadmor, “Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems,” Acta Numer., vol. 12, pp. 451–512, 2003
work page 2003
-
[7]
On discretely entropy conservative and entropy stable discontinuous Galerkin methods,
J. Chan, “On discretely entropy conservative and entropy stable discontinuous Galerkin methods,” J. Comput. Phys., vol. 362, pp. 346–374, 2018
work page 2018
Show all 29 references
-
[8]
Review of entropy stable discontinuous Galerkin methods for systems of conservation laws on unstructured simplex meshes,
T. Chen and C.-W. Shu, “Review of entropy stable discontinuous Galerkin methods for systems of conservation laws on unstructured simplex meshes,” CSIAM Trans. Appl. Math., vol. 1, no. 1, pp. 1–52, 2020
2020
-
[9]
A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler equations,
S. Hennemann, A. M. Rueda-Ramírez, F. J. Hindenlang, and G. J. Gassner, “A provably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler equations,” J. Comput. Phys., vol. 426, p. 109935, 2021
2021
-
[10]
Subcell limiting strategies for discontinuous galerkin spectral element methods,
A. M. Rueda-Ramírez, W. Pazner, and G. J. Gassner, “Subcell limiting strategies for discontinuous galerkin spectral element methods,” Comput. Fluids, vol. 247, p. 105627, 2022
2022
-
[11]
Formulation of entropy-stable schemes for the multicomponent compressible Euler equations,
A. Gouasmi, K. Duraisamy, and S. M. Murman, “Formulation of entropy-stable schemes for the multicomponent compressible Euler equations,” Comput. Methods Appl. Mech. Eng., vol. 363, p. 112912, 2020
2020
-
[12]
A conservative discontinuous Galerkin discretization for the chemically reacting Navier-Stokes equations,
R. F. Johnson and A. D. Kercher, “A conservative discontinuous Galerkin discretization for the chemically reacting Navier-Stokes equations,” J. Comput. Phys., vol. 423, p. 109826, 2020
2020
-
[13]
Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case,
E. J. Ching, R. F. Johnson, and A. D. Kercher, “Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case,” J. Comput. Phys., vol. 505, p. 112878, 2024
2024
-
[14]
NASA Glenn coefficients for calculating thermodynamic properties of individual species,
B. J. McBride, “NASA Glenn coefficients for calculating thermodynamic properties of individual species,” tech. rep., National Aeronautics and Space Administration, John H. Glenn Research Center, 2002
2002
-
[15]
Entropy stable discretization of compressible flows in thermochemical nonequi- librium,
M. A. Hansen and T. C. Fisher, “Entropy stable discretization of compressible flows in thermochemical nonequi- librium,” tech. rep., Sandia National Laboratories, 2019
2019
-
[16]
High-order methods for hypersonic flows with strong shocks and real chemistry,
A. Peyvan, K. Shukla, J. Chan, and G. Karniadakis, “High-order methods for hypersonic flows with strong shocks and real chemistry,”J. Comput. Phys., p. 112310, 2023
2023
-
[18]
Simulation of real gas effects in supersonic methane jets using a tabulated equation of state with a discontinuous Galerkin spectral element method,
F. Hempert, S. Boblest, T. Ertl, F. Sadlo, P. Offenhäuser, C. Glass, M. Hoffmann, A. Beck, C.-D. Munz, and U. Iben, “Simulation of real gas effects in supersonic methane jets using a tabulated equation of state with a discontinuous Galerkin spectral element method,” Comput. Fl...
2017
-
[19]
Trixi.jl: Adaptive high-order numerical simulations of hyperbolic PDEs in Julia
M. Schlottke-Lakemper, G. J. Gassner, H. Ranocha, and A. R. Winters, “Trixi.jl: Adaptive high-order numerical simulations of hyperbolic PDEs in Julia.” https://github.com/trixi-framework/Trixi.jl, 08 2020
2020
-
[20]
Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing,
H. Ranocha, M. Schlottke-Lakemper, A. R. Winters, E. Faulhaber, J. Chan, and G. J. Gassner, “Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing,” arXiv:2108.06476, 2021
2021 arXiv
-
[21]
On high order strong stability preserving Runge-Kutta and multi step time discretizations,
S. Gottlieb, “On high order strong stability preserving Runge-Kutta and multi step time discretizations,”J. Sci. Comput., vol. 25, pp. 105–128, 2005
2005
-
[22]
Reproducibility repository for “Entropy-stable fluxes for high- order Discontinuous Galerkin simulations of high-enthalpy flows
G. Oblapenko, A. Tarnovskiy, and M. Torrilhon, “Reproducibility repository for “Entropy-stable fluxes for high- order Discontinuous Galerkin simulations of high-enthalpy flows”.” https://github.com/knstmrd/paper_ ec_trixi_chem, 2024
2024
-
[23]
Nagnibeda and E
E. Nagnibeda and E. Kustova, Non-equilibrium reacting gas flows: kinetic theory of transport and relaxation processes. Springer Science & Business Media, 2009
2009
-
[24]
Review of chemical-kinetic problems of future NASA missions. I-Earth entries,
C. Park, “Review of chemical-kinetic problems of future NASA missions. I-Earth entries,”J. Thermophys. Heat Transfer, vol. 7, no. 3, pp. 385–398, 1993
1993
-
[25]
p4est: Scalable algorithms for parallel adaptive mesh refinement on forests of octrees,
C. Burstedde, L. C. Wilcox, and O. Ghattas, “p4est: Scalable algorithms for parallel adaptive mesh refinement on forests of octrees,” SIAM J. Sci. Comp., vol. 33, no. 3, pp. 1103–1133, 2011
2011
-
[26]
Validation of the unstructured DLR-TAU-Code for hypersonic flows,
A. Mack and V . Hannemann, “Validation of the unstructured DLR-TAU-Code for hypersonic flows,” in32nd AIAA Fluid dynamics conference and exhibit, p. 3111, 2002
2002
-
[27]
A closely coupled experimental and numerical approach for hypersonic and high enthalpy flow investigations utilising the HEG shock tunnel and the DLR TAU code,
K. Hannemann, J. M. Schramm, A. Wagner, S. Karl, and V . Hannemann, “A closely coupled experimental and numerical approach for hypersonic and high enthalpy flow investigations utilising the HEG shock tunnel and the DLR TAU code,” Tech. Rep. RTO-EN-A VT-186, VKI, 2010. 7 Entrop...
2010
-
[28]
A sequel to AUSM: AUSM+,
M.-S. Liou, “A sequel to AUSM: AUSM+,” J. Comput. Phys., vol. 129, no. 2, pp. 364–382, 1996
1996
-
[29]
A flux splitting scheme with high-resolution and robustness for discontinuities,
Y . Wada and M.-S. Liou, “A flux splitting scheme with high-resolution and robustness for discontinuities,” in32nd Aerospace Sciences Meeting and Exhibit, 1994. 8
1994
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.