Pith. sign in

REVIEW 2 major objections 5 minor 2 cited by

Conservative, pressure-equilibrium-preserving discontinuous Galerkin method for compressible, multicomponent flows

T0 review · 2 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A discontinuous Galerkin scheme that preserves velocity and pressure equilibria while conserving total energy for smooth multicomponent flows.

desk verdict New and useful DG scheme with clean algebra, but Theorem 17 doesn't cover the mixed uniform/non-uniform faces the actual scheme and tests use. read the letter →

arxiv 2501.12532 v1 pith:MRJTP7DH submitted 2025-01-21 math.NA cs.NAphysics.flu-dyn

Eric J. Ching , Ryan F. Johnson , Andrew D. Kercher This is my paper · ORCID
classification math.NAcs.NAphysics.flu-dyn MSC 65M6076N15
keywords discontinuousGalerkinmethodmulticomponentflowpressureequilibriumpreservationvelocitytotalenergyconservationcorrectiontermscompressibleEulerequationsthermallyperfectgas
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

This paper attacks a long-standing failure mode of conservative numerical methods for compressible multicomponent flows: at interfaces between gases with different specific heats, the computed pressure develops spurious oscillations and the solver can diverge. The authors prove that a standard discontinuous Galerkin (DG) discretization of the conservative Euler equations preserves a constant velocity under Lax-Friedrichs-type fluxes but cannot preserve a constant pressure. They therefore replace the total-energy conservation equation with a pressure-evolution equation, which preserves pressure equilibrium, and then add correction terms that restore semidiscrete total-energy conservation. A modified choice of auxiliary variables makes these corrections invisible to the velocity and pressure equations in constant-pressure, constant-velocity flows, and face-based corrections cover elementwise-constant states. The result is a DG scheme that, for smooth thermally perfect multicomponent flows initially at uniform pressure and velocity, keeps pressure and velocity flat while conserving total energy, and it stays stable for hundreds of advection periods in one, two, and three dimensions.

What carries the argument

The load-bearing object is the modified correction term $r_k = \alpha(\hat{z}_k - \bar{z})$, where $\alpha$ is chosen so that the corrected residual satisfies a discrete total-energy balance, and the auxiliary variable is $z = (v \sum_i W_i \partial(\rho e_t)/\partial C_i,\; P,\; \partial(\rho e_t)/\partial C_1, \ldots, \partial(\rho e_t)/\partial C_{n_s})^T$, with $v$ the velocity and $W_i$ the molar masses. Replacing the usual auxiliary variable $w$ by $z$ makes the correction drop out of the velocity and pressure equations whenever $P=P_0$ and $v=v_0$, because $\Delta z_{\rho v} = v_0 \sum_i W_i \Delta(\partial \rho e_t/\partial C_i)$ cancels identically, while $\Delta z_P = 0$. A face-based correction $\beta \llbracket z \rrbracket$ handles elementwise-constant states, and the nonnegativity of the denominator follows from $\Delta w^T \Delta z = \sum_i [\Delta(\partial \rho e_t/\partial C_i)]^2 \ge 0$ in equilibrium.

What would settle it

Simulate an elementwise-constant, two-species interface with equal pressure and velocity, or a smooth thermal bubble, using the modified scheme with third-order Runge-Kutta and systematically halved time steps; if pressure or velocity errors rise above round-off level, or if the global energy error does not shrink at third order with the time step, the central theorem is contradicted. A shock-tube version of the same setup tests whether the promised preservation extends past the stated smooth, constant-equilibrium assumptions.

Watch

Extended reading notes

Core claim

The central claim is Theorem 17: in a constant-pressure, constant-velocity flow, the DG scheme (3.7) with the modified correction term (3.17) preserves both velocity and pressure equilibria and, assuming continuity in time, conserves total energy at the semidiscrete level. The key is that the correction term is built from an auxiliary variable z whose momentum components are v Σ_i W_i ∂(ρe_t)/∂C_i and whose pressure component is P; in the equilibrium state the correction's contribution to the velocity equation cancels term-by-term and its pressure component is exactly zero. A face-based correction of the form β ⟦z⟧ extends the argument to elementwise-constant solutions with inter-element jumps, and the numerical total-energy flux is defined so that local conservation holds across mixed element types. The numerical experiments show pressure errors at finite-precision levels for the modified scheme over 100 advection periods in one, two, and three dimensions, with third-order temporal convergence of the total-energy error when paired with SSPRK3 time integration.

Load-bearing premise

The guarantee is proven only for smooth flows sitting exactly at constant pressure and velocity; away from that state, at shocks, under nonuniform equilibria, or with real-fluid equations of state, the modified correction can lose its nonnegative denominator and the theorem's cancellations are not shown.

Editorial extensions

If this is right

  • Smooth interfacial flows of thermally perfect gas mixtures initialized at uniform pressure and velocity can be advected for many periods with only round-off-level pressure error and no artificial viscosity or limiting.
  • Total-energy conservation is restored to the pressure-based formulation at the semidiscrete level; with explicit third-order SSPRK time stepping the global energy error converges at third order in the time step.
  • The method works on unstructured and curved meshes in two and three dimensions, where the standard total-energy DG scheme with either colocated or overintegrated quadrature diverges in the tested cases.
  • The original correction term causes spurious production and destruction of zero-concentration species; the modified correction with zeroed components prevents this, which matters for reacting-flow extensions.
  • For polynomial degree p=1 the modified correction has a smaller denominator and can degrade accuracy or cause divergence on coarse meshes; p≥2 cases are largely unaffected.

Reading between the lines

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

  • Away from constant-pressure and constant-velocity equilibria, the modified correction's denominator can become negative, so a general-purpose implementation would likely need to switch back to the original correction term or add stabilization; the paper sketches such a hybrid but does not develop it.
  • The construction is not obviously restricted to the Euler equations: any system where a secondary conservation law is desired alongside a nonconservative primary formulation could use the same auxiliary-variable switching, though the paper does not explore this.
  • Because the method stores pressure instead of total energy, temperature is obtained directly from the state, which may simplify future real-fluid and chemically reacting extensions; the paper only demonstrates thermally perfect gas mixtures.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 5 minor

Summary. The paper addresses pressure- and velocity-equilibrium preservation in discontinuous Galerkin (DG) discretizations of the compressible, multicomponent Euler equations for smooth, inviscid flows of thermally perfect gas mixtures. The authors first analyze a standard conservative DG formulation and show that it can preserve velocity equilibrium under flux constraints but generally fails to preserve pressure equilibrium. They then introduce a pressure-based DG formulation, replacing the total-energy equation with a pressure-evolution equation, and restore semidiscrete total-energy conservation using the correction terms of Abgrall and co-workers. Because these corrections destroy the equilibrium-preservation property, the authors propose a modified correction term and a face-based correction, and they prove (Theorem 17) that the resulting scheme is velocity-equilibrium-preserving, pressure-equilibrium-preserving, and semidiscretely energy-conservative. Numerical experiments in one, two, and three dimensions, including curved-mesh cases, demonstrate stability, pressure-equilibrium preservation, and energy-conservation behavior for smooth thermal-bubble and Gaussian-wave problems, with an acknowledged robustness limitation for p=1 on coarse meshes.

Significance. If the results hold, the paper makes a valuable contribution to a long-standing problem: constructing fully conservative schemes that exactly preserve pressure equilibrium at fluid interfaces for thermally perfect multicomponent gases, a regime where prior conservative schemes either produced spurious oscillations or required calorically perfect gas assumptions. The paper's strengths include clean algebraic proofs of the equilibrium conditions and of the modified correction-term identities, transparent discussion of the tradeoffs involved (semidiscrete versus fully discrete conservation, exact-integration caveats, and p=1 fragility), and a substantial set of numerical experiments that demonstrate the practical benefit of the proposed P3 scheme. The independent verification on 1D/2D/3D thermal-bubble problems and on curved meshes provides convincing evidence that the proposed corrections are viable in realistic smooth-flow settings. The explicit statement of limitations and the clear distinction between proven and empirically observed behavior are commendable.

major comments (2)
  1. [Sec. 3.3.1 (Theorem 17) and Sec. 3.3 (Eqs. 3.19–3.20)] The proof of Theorem 17 covers two pure cases: elementwise-nonuniform solutions with only the elementwise correction (3.17), and elementwise-constant solutions with only the face-based correction (3.19). However, the P3 scheme as defined applies the face correction through the numerical flux at every face where it is activated, so at a face shared by a uniform and a non-uniform element the beta JbzK term enters the residual of the non-uniform element, whose elementwise correction is also active. The note after Eq. (3.20) explicitly relies on this mixed-face configuration to maintain local energy conservation, but neither case (a) nor case (b) of the proof establishes energy conservation or equilibrium preservation for it. Since the scheme and theorem are stated for the combined corrections, the central theoretical claim is incomplete for mixed faces. Please either provide a proof for the mixed-face case or restrict the theorem and algorithmic guarantee to configurations in which each face connects elements of the same type, and revise the text accordingly.
  2. [Sec. 3.3, Remark 19] The proposed fix for preserving zero species concentrations (zeroing components of r_k, and similarly for the face correction) is applied after alpha and beta have been determined from the energy-balance conditions (3.10) and (3.19). The paper does not analyze whether this post-hoc zeroing preserves the semidiscrete total-energy identity (3.15); in general, removing a nonzero component of the correction term will break the equality (3.10). Since the introduction advertises that the modifications 'enable exact preservation of pressure equilibrium, velocity equilibrium, and zero species concentrations (while maintaining semidiscrete total-energy conservation)', the compatibility of the zero-species fix with energy conservation should be either proved or explicitly qualified. The numerical test in Sec. 4.2 demonstrates zero-species preservation but does not report the energy error for that configuration.
minor comments (5)
  1. [Sec. 4.1 and Fig. 4.1] The text states that the L2 error is computed at t = 0.2 (one advection period), while the caption of Figure 4.1 states t = 1; please correct the inconsistency.
  2. [Sec. 4.4] In the temporal-convergence paragraph, the results are said to be presented in Figure 4.2, but the correct reference appears to be Figure 4.7.
  3. [Sec. 2.1.1, Eq. (2.13)] The clause 'is the component corresponding to the ith species concentration,.' contains a stray comma; also, the 'if and only if' statement of Lemma 1 should carry the exact-integration qualification stated in Remark 4 directly in the lemma statement, rather than only in a later remark.
  4. [Sec. 3.3, Eq. (3.20)] The notation q wT F (y) y is nonstandard and not introduced; replacing it with the average notation used elsewhere, e.g., {{w^T F}}·n, would improve readability.
  5. [References] References [20] and [50] appear to be the same paper (Johnson & Kercher, JCP 423, 2020, 109826); please deduplicate the bibliography.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the scheme's energy conservation is enforced by construction of the correction coefficients, and the equilibrium-preservation claims are proven from the definitions; self-citations are present but not load-bearing.

full rationale

The paper does not derive its central claims from fitted data or from self-citation. The correction coefficients alpha and beta are deliberately constructed so that the corrected residual satisfies a discrete total-energy balance: alpha is defined by equation (3.10) so that the total-energy residual equals the boundary flux, and beta is defined by equation (3.19) so that the face correction satisfies the energy-jump relation (3.21). Stating Proposition 10 and the energy-conservation part of Theorem 17 after this construction is a report of the construction, not a circular prediction; the construction itself is an explicit design step. The equilibrium-preservation part of Theorem 17 is also obtained by direct algebra from the modified auxiliary variable z in equation (3.18): in a constant-pressure, constant-velocity state, z_P = P and the momentum component of z is v times the species-component sum, so the correction contributions to d(v)/dt and dP/dt cancel identically. This is a proof from the definitions, not an assumption of the result. The paper relies on prior work by Abgrall and others for the correction-term framework, but those are external references, and the authors' own self-citations (e.g., [3] for the smooth thermal-bubble test setup and [20,50] for the implementation framework) are contextual rather than load-bearing for the theorem. The proof of Theorem 17 does split into elementwise-nonuniform and elementwise-constant cases, and the treatment of a face connecting a uniform element to a nonuniform element is less explicit; however, that is a potential completeness or correctness gap in the proof, not a circularity. Overall, the derivation chain is self-contained: the claimed properties follow from the stated formulas and algebraic identities, and there is no fitted input disguised as a prediction.

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

No new physical entities are introduced. The central construction modifies the auxiliary vector w to z and adds face-based corrections, both numerical devices. The only hand-set numbers are two denominator thresholds (1e-7 and 1e-6). The main axioms are the thermally perfect-gas EOS, the DG/Lax-Friedrichs discretization, the integration-accuracy assumption in the equilibrium lemmas, and continuity in time for the semidiscrete energy-conservation statement.

free parameters (2)
  • Elementwise correction denominator cutoff for alpha = 1e-7
    Alpha is set to zero if the denominator sum_k (bw_k-w)^T(bz_k-z) is below this threshold. This avoids division by small denominators but can locally forfeit the energy-conservation guarantee.
  • Face-based correction denominator cutoff for beta = 1e-6
    Beta is set to zero if JbwK^T JbzK is below this threshold. A stabilization choice that can locally break the correction's effect.
assumptions (5)
  • domain assumption Thermally perfect gas mixture equation of state: P = R0 T sum_i C_i, with species internal energies computed from NASA polynomial fits.
    All derivations and test cases assume this EOS; real-fluid equations of state are explicitly deferred to future work.
  • domain assumption The spatial discretization is a standard DG method with Lax-Friedrichs numerical flux and polynomial test space V_h^p.
    The equilibrium conditions are derived for this DG discretization and this numerical flux; other fluxes would require separate analysis.
  • domain assumption Exact or sufficiently accurate numerical integration of volume terms such as int_kappa dphi/dx P dx in Lemma 1 and Proposition 6.
    Remark 4 states this assumption explicitly. With standard overintegration, pressure-equilibrium preservation can fail, so the theoretical claim is conditional on integration accuracy.
  • standard math Continuity in time for the semidiscrete energy-conservation statement.
    Proposition 10 and Theorem 17 prove semidiscrete conservation; explicit Runge-Kutta integrators only give time-step-convergent energy error, as noted in Remark 12.
  • domain assumption Local constant pressure and velocity (P=P0, v=v0) when proving nonnegativity of the modified correction denominator.
    Remark 16 establishes denominator nonnegativity only in this regime; outside it, the modified alpha formula is not guaranteed to be well behaved.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Conservative, pressure-equilibrium-preserving discontinuous Galerkin method for compressible, multicomponent flows." pith.science (2026). https://pith.science/paper/MRJTP7DH

@misc{pith2026250112532,
  author       = {Pith},
  title        = {Pith review of: Conservative, pressure-equilibrium-preserving discontinuous Galerkin method for compressible, multicomponent flows},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MRJTP7DH}},
  note         = {Machine review of arXiv:2501.12532}
}
read the original abstract

This paper concerns preservation of velocity and pressure equilibria in smooth, compressible, multicomponent flows in the inviscid limit. First, we derive the velocity-equilibrium and pressure-equilibrium conditions of a standard discontinuous Galerkin method that discretizes the conservative form of the compressible, multicomponent Euler equations. We show that under certain constraints on the numerical flux, the scheme is velocity-equilibrium-preserving. However, standard discontinuous Galerkin schemes are not pressure-equilibrium-preserving. Therefore, we introduce a discontinuous Galerkin method that discretizes the pressure-evolution equation in place of the total-energy conservation equation. Semidiscrete conservation of total energy, which would otherwise be lost, is restored via the correction terms of [Abgrall, J. Comput. Phys., 372, 2018, pp. 640-666] and [Abgrall et al., J. Comput. Phys., 453, 2022, 110955]. Since the addition of the correction terms prevents exact preservation of pressure and velocity equilibria, we propose modifications that then lead to a velocity-equilibrium-preserving, pressure-equilibrium-preserving, and (semidiscretely) energy-conservative discontinuous Galerkin scheme, although there are certain tradeoffs. Additional extensions are also introduced. We apply the developed scheme to smooth, interfacial flows involving mixtures of thermally perfect gases initially in pressure and velocity equilibria to demonstrate its performance in one, two, and three spatial dimensions.

Figures

Figures reproduced from arXiv: 2501.12532 by the authors.

Figure 4.1
Figure 4.1. Convergence under grid refinement, with h = 0.02, for the multicomponent Gaussian-wave test. The L2 error of the normalized state with respect to the exact solution at t = 1 is computed. The dashed lines represent convergence rates of p + 1. 19 [PITH_FULL_IMAGE:figures/full_fig_p019_4_1.png] view at source ↗
Figure 4.2
Figure 4.2. Convergence of total-energy conservation with respect to time-step size, where [PITH_FULL_IMAGE:figures/full_fig_p020_4_2.png] view at source ↗
Figure 4.3
Figure 4.3. p = 3 solutions to the advection of a nitrogen/n-dodecane thermal bubble at v = 600 m/s. P1: uncorrected pressure–based DG scheme; P2: pressure-based DG scheme with original correction term (3.8); P3: proposed pressure-based DG scheme (Section 3.3); E1: total-energy-based DG scheme with overintegration; E2: total-energy-based DG scheme with colocated integration; Exact: exact solution. Next, to demonstrate the failu… view at source ↗
Figures from the paper (14 more)
Figure 4.4
Figure 4.4. Figure 4.4: YO2 profiles at t = 100τ obtained from p = 3 solutions to the advection of a nitrogen/n-dodecane thermal bubble at v = 600 m/s [PITH_FULL_IMAGE:figures/full_fig_p022_4_4.png]
Figure 4.5
Figure 4.5. Figure 4.5: p = 3 solution to the advection of a nitrogen/n-dodecane thermal bubble at v = 1 m/s. P1: uncorrected pressure– based DG scheme; P2: pressure-based DG scheme with original correction term (3.8); P3: proposed pressure-based DG scheme (Section 3.3); E1: total-energy-ba…
Figure 4.6
Figure 4.6. Figure 4.6: Initial p = 3 density and pressure fields, superimposed by the grid, for the two-dimensional advection of a nitrogen/n-dodecane thermal bubble. We first assess the temporal convergence of the pressure-based schemes in this two-dimensional setting. Four time-step size…
Figure 4.7
Figure 4.7. Figure 4.7: Convergence of total-energy conservation with respect to time-step size, where [PITH_FULL_IMAGE:figures/full_fig_p024_4_7.png]
Figure 4.8
Figure 4.8. Figure 4.8: Temporal variation of error in pressure for [PITH_FULL_IMAGE:figures/full_fig_p025_4_8.png]
Figure 4.9
Figure 4.9. Figure 4.9: p = 3 solution to two-dimensional advection of a nitrogen/n-dodecane thermal bubble at t = 100τ computed with the uncorrected DG scheme (P1). 25 [PITH_FULL_IMAGE:figures/full_fig_p025_4_9.png]
Figure 4.10
Figure 4.10. Figure 4.10: p = 3 solution to two-dimensional advection of a nitrogen/n-dodecane thermal bubble at t = 100τ computed with the proposed pressure-based DG scheme (P3). Curved elements. Finally, we recompute this case using curved elements of quadratic order with a CFL of 0.4. To …
Figure 4.11
Figure 4.11. Figure 4.11: Temporal variation of error in pressure for [PITH_FULL_IMAGE:figures/full_fig_p026_4_11.png]
Figure 4.12
Figure 4.12. Figure 4.12: p = 3 solution to two-dimensional advection of a nitrogen/n-dodecane thermal bubble at t = 100τ computed with the uncorrected DG scheme (P1) on a curved grid. 4.5. Three-dimensional, high-velocity thermal bubble Our final test case is a three-dimensional version of …
Figure 4
Figure 4. Figure 4: presents the temporal variation of percent error in pressure, sampled every [PITH_FULL_IMAGE:figures/full_fig_p027_4.png]
Figure 4.13
Figure 4.13. Figure 4.13: p = 3 solution to two-dimensional advection of a nitrogen/n-dodecane thermal bubble at t = 100τ computed with the proposed pressure-based DG scheme (P3) on a curved grid. shown that under certain constraints on the numerical flux, the scheme is velocity-equilibrium-…
Figure 4.14
Figure 4.14. Figure 4.14: T = 890 K isosurfaces colored by density at t = 0 for the three-dimensional advection of a nitrogen/n-dodecane thermal bubble. 0 20 40 60 80 100 t/τ 10−7 10−5 10−3 10−1 101 max x |P − P0| /P0 (%) P1 P2 P3 E1 E2 [PITH_FULL_IMAGE:figures/full_fig_p029_4_14.png]
Figure 4.15
Figure 4.15. Figure 4.15: Temporal variation of error in pressure for [PITH_FULL_IMAGE:figures/full_fig_p029_4_15.png]
Figure 4.16
Figure 4.16. Figure 4.16: T = 890 K isosurfaces colored by density at t = 100τ for the three-dimensional advection of a nitrogen/n-dodecane thermal bubble computed with the P1 and P3 schemes. nonlinear function of the state, although this requires further investigation. In certain cases, it …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Real Fluid Quasi-Conservative Method: Mechanical Equilibrium Mechanism and Liquid-Upwind Anomaly

    math-ph 2026-07 conditional novelty 6.0 of 10

    The RFQC real-fluid solver's liquid-upwind anomaly is a start-up failure caused by a jump in the affine slope xi and by thermodynamic re-projection canceling the pressure rise; smoothing the initial data fixes it.

  2. A Central Differential Flux with High-Order Dissipation for Robust Simulations of Transcritical Flows

    physics.flu-dyn 2025-08 conditional novelty 5.0 of 10

    CDHD is a central flux plus high-order DOTRS dissipation for primitive-variable Euler updates, giving second-order accuracy with O(Δx^4) dissipation for transcritical flows.

Reference graph

Works this paper leans on

56 extracted references · 51 canonical work pages · cited by 2 Pith papers

  1. [1]

    Abgrall, A general framework to construct schemes satisfying additional conservation relations

    R. Abgrall, A general framework to construct schemes satisfying additional conservation relations. application to entropy conservative and entropy dissipative schemes, Journal of Computational Physics 372 (2018) 640–666

    work page 2018

  2. [2]

    Abgrall, P

    R. Abgrall, P. Öffner, H. Ranocha, Reinterpretation and extension of entropy correction terms for residual distribution and discontinuous Galerkin schemes: application to structure preserving discretization, Journal of Computational Physics 453 (2022) 110955

    work page 2022

  3. [3]

    E. J. Ching, R. F. Johnson, A. D. Kercher, A note on reducing spurious pressure oscillations in fully conservative discon- tinuous Galerkin simulations of multicomponent flows, arXiv preprint arXiv:2310.17792 (2023)

    work page Pith review arXiv 2023

  4. [4]

    Abgrall, S

    R. Abgrall, S. Karni, Computations of compressible multifluids, Journal of Computational Physics 169 (2) (2001) 594 –

    work page 2001

  5. [5]

    Karni, Viscous shock profiles and primitive formulations, SIAM journal on numerical analysis 29 (6) (1992) 1592–1609

    S. Karni, Viscous shock profiles and primitive formulations, SIAM journal on numerical analysis 29 (6) (1992) 1592–1609

    work page 1992

  6. [6]

    Karni, Multicomponent flow calculations by a consistent primitive algorithm, Journal of Computational Physics 112 (1) (1994) 31 – 43.doi:https://doi.org/10.1006/jcph.1994.1080

    S. Karni, Multicomponent flow calculations by a consistent primitive algorithm, Journal of Computational Physics 112 (1) (1994) 31 – 43.doi:https://doi.org/10.1006/jcph.1994.1080

    arXiv 1994

  7. [7]

    Terashima, M

    H. Terashima, M. Koshi, Approach for simulating gas liquid-like flows under supercritical pressures using a high-order central differencing scheme, Journal of Computational Physics 231 (20) (2012) 6907 – 6923.doi:https://doi.org/10. 1016/j.jcp.2012.06.021

    work page 2012

  8. [8]

    Kawai, H

    S. Kawai, H. Terashima, H. Negishi, A robust and accurate numerical method for transcritical turbulent flows at super- critical pressure with an arbitrary equation of state, Journal of Computational Physics 300 (2015) 116–135

    work page 2015

Show all 56 references
  1. [9]

    Karni, Hybrid multifluid algorithms, SIAM Journal on Scientific Computing 17 (5) (1996) 1019–1039

    S. Karni, Hybrid multifluid algorithms, SIAM Journal on Scientific Computing 17 (5) (1996) 1019–1039

    1996

  2. [10]

    Fedkiw, X.-D

    R. Fedkiw, X.-D. Liu, S. Osher, A general technique for eliminating spurious oscillations in conservative schemes for multiphase and multispecies Euler equations, International Journal of Nonlinear Sciences and Numerical Simulation 3 (2) (2002) 99–106

    2002

  3. [11]

    Y. Lv, M. Ihme, Discontinuous Galerkin method for multicomponent chemically reacting flows and combustion, Journal of Computational Physics 270 (2014) 105 – 137.doi:https://doi.org/10.1016/j.jcp.2014.03.029

    2014 doi

  4. [12]

    B. Boyd, D. Jarrahbashi, A diffuse-interface method for reducing spurious pressure oscillations in multicomponent trans- critical flow simulations, Computers & Fluids 222 (2021) 104924

    2021

  5. [13]

    Gaburro, W

    E. Gaburro, W. Boscheri, S. Chiocchetti, M. Ricchiuto, Discontinuous Galerkin schemes for hyperbolic systems in non- conservativevariables: quasi-conservativeformulationwithsubcellfinitevolumecorrections, ComputerMethodsinApplied Mechanics and Engineering 431 (2024) 117311

    2024

  6. [14]

    W. H. Reed, T. Hill, Triangular mesh methods for the neutron transport equation, Tech. rep., Los Alamos Scientific Lab., N. Mex.(USA) (1973)

    1973

  7. [15]

    Bassi, S

    F. Bassi, S. Rebay, High-order accurate discontinuous finite element solution of the 2D Euler equations, Journal of Computational Physics 138 (2) (1997) 251–285

    1997

  8. [16]

    F.Bassi, S.Rebay, Ahigh-orderaccuratediscontinuousfiniteelementmethodforthenumericalsolutionofthecompressible Navier–Stokes equations, Journal of Computational Physics 131 (2) (1997) 267–279

    1997

  9. [17]

    Cockburn, C.-W

    B. Cockburn, C.-W. Shu, The Runge–Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems, Journal of Computational Physics 141 (2) (1998) 199–224

    1998

  10. [18]

    Cockburn, G

    B. Cockburn, G. Karniadakis, C.-W. Shu, The development of discontinuous Galerkin methods, in: Discontinuous Galerkin Methods, Springer, 2000, pp. 3–50

    2000

  11. [19]

    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, M. Visbal, High-order CFD methods: Current status and perspective, International Journal for Numerical Methods in...

    2013 doi

  12. [21]

    Bando, Towards high-performance discontinuous Galerkin simulations of reacting flows using Legion, Ph.D

    K. Bando, Towards high-performance discontinuous Galerkin simulations of reacting flows using Legion, Ph.D. thesis, Stanford University (2023)

    2023

  13. [22]

    Johnsen, T

    T. Johnsen, T. Colonius, Implementation of WENO schemes in compressible multicomponent flow problems, Journal of Computational Physics 219 (2) (2006) 715 – 732.doi:https://doi.org/10.1016/j.jcp.2006.04.018

    2006 doi

  14. [23]

    Franchina, M

    N. Franchina, M. Savini, F. Bassi, Multicomponent gas flow computations by a discontinuous Galerkin scheme using L2-projection of perfect gas EOS, Journal of Computational Physics 315 (2016) 302–322

    2016

  15. [24]

    Fujiwara, Y

    Y. Fujiwara, Y. Tamaki, S. Kawai, Fully conservative and pressure-equilibrium preserving scheme for compressible multi- component flows, Journal of Computational Physics 478 (2023) 111973

    2023

  16. [25]

    Terashima, N

    H. Terashima, N. Ly, M. Ihme, Approximately pressure-equilibrium-preserving scheme for fully conservative simulations of compressible multi-species and real-fluid interfacial flows, Journal of Computational Physics 524 (2025) 113701

    2025

  17. [26]

    Chen, C.-W

    T. Chen, C.-W. Shu, Review of entropy stable discontinuous Galerkin methods for systems of conservation laws on unstructured simplex meshes, CSIAM Transactions on Applied Mathematics 1 (1) (2020) 1–52

    2020

  18. [27]

    Gaburro, P

    E. Gaburro, P. Öffner, M. Ricchiuto, D. Torlo, High order entropy preserving ADER-DG schemes, Applied Mathematics and Computation 440 (2023) 127644

    2023

  19. [28]

    Mantri, P

    Y. Mantri, P. Öffner, M. Ricchiuto, Fully well-balanced entropy controlled discontinuous Galerkin spectral element method for shallow water flows: Global flux quadrature and cell entropy correction, Journal of Computational Physics 498 (2024) 112673

    2024

  20. [29]

    Alberti, E

    L. Alberti, E. Carnevali, A. Colombo, A. Crivellini, An entropy conserving/stable discontinuous Galerkin solver in entropy variables based on the direct enforcement of entropy balance, Journal of Computational Physics 508 (2024) 113007. 31

    2024

  21. [30]

    Abgrall, S

    R. Abgrall, S. Busto, M. Dumbser, A simple and general framework for the construction of thermodynamically compatible schemes for computational fluid and solid mechanics, Applied Mathematics and Computation 440 (2023) 127629

    2023

  22. [31]

    M. J. Castro, P. G. LeFloch, M. L. Muñoz-Ruiz, C. Parés, Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes, Journal of Computational Physics 227 (17) (2008) 8107–8129

    2008

  23. [32]

    Abgrall, S

    R. Abgrall, S. Karni, A comment on the computation of non-conservative products, Journal of Computational Physics 229 (8) (2010) 2759–2763

    2010

  24. [33]

    Giovangigli, Multicomponent flow modeling, Birkhauser, Boston, 1999

    V. Giovangigli, Multicomponent flow modeling, Birkhauser, Boston, 1999

    1999

  25. [34]

    B. J. McBride, S. Gordon, M. A. Reno, Coefficients for calculating thermodynamic and transport properties of individual species (1993)

    1993

  26. [35]

    B. J. McBride, M. J. Zehe, S. Gordon, NASA Glenn coefficients for calculating thermodynamic properties of individual species (2002)

    2002

  27. [36]

    E. J. Ching, R. F. Johnson, S. Burrows, J. Higgs, A. D. Kercher, Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Navier-Stokes equations, arXiv preprint arXiv:2310.17637 (2023)

    2023 arXiv

  28. [37]

    Toro, Riemann solvers and numerical methods for fluid dynamics: A practical introduction, Springer Science & Business Media, 2013

    E. Toro, Riemann solvers and numerical methods for fluid dynamics: A practical introduction, Springer Science & Business Media, 2013

    2013

  29. [38]

    Gottlieb, C

    S. Gottlieb, C. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM review 43 (1) (2001) 89–112

    2001

  30. [39]

    Y. Lv, Development of a nonconservative discontinuous Galerkin formulation for simulations of unsteady and turbulent flows, International Journal for Numerical Methods in Fluids 92 (5) (2020) 325–346

    2020

  31. [40]

    Chan, On discretely entropy conservative and entropy stable discontinuous Galerkin methods, Journal of Computational Physics 362 (2018) 346–374

    J. Chan, On discretely entropy conservative and entropy stable discontinuous Galerkin methods, Journal of Computational Physics 362 (2018) 346–374

    2018

  32. [41]

    Atkins, C

    H. Atkins, C. Shu, Quadrature-free implementation of discontinuous Galerkin methods for hyperbolic equations, ICASE Report 96-51, 1996, Tech. rep., NASA Langley Research Center, nASA-CR-201594 (August 1996)

    1996

  33. [42]

    H. L. Atkins, C.-W. Shu, Quadrature-free implementation of discontinuous Galerkin method for hyperbolic equations, AIAA Journal 36 (5) (1998) 775–782

    1998

  34. [43]

    Nordström, T

    J. Nordström, T. Lundquist, Summation-by-parts in time, Journal of Computational Physics 251 (2013) 487–499

    2013

  35. [44]

    Ranocha, M

    H. Ranocha, M. Sayyari, L. Dalcin, M. Parsani, D. I. Ketcheson, Relaxation Runge–Kutta methods: Fully discrete explicit entropy-stable schemes for the compressible Euler and Navier–Stokes equations, SIAM Journal on Scientific Computing 42 (2) (2020) A612–A638

    2020

  36. [45]

    Ranocha, L

    H. Ranocha, L. Lóczi, D. I. Ketcheson, General relaxation methods for initial-value problems with application to multistep schemes, Numerische Mathematik 146 (4) (2020) 875–906

    2020

  37. [46]

    C. Wang, X. Zhang, C.-W. Shu, J. Ning, Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations, Journal of Computational Physics 231 (2) (2012) 653–665

    2012

  38. [47]

    Zhang, C.-W

    X. Zhang, C.-W. Shu, On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes, Journal of Computational Physics 229 (23) (2010) 8918–8934

    2010

  39. [48]

    E. J. Ching, R. F. Johnson, A. D. Kercher, Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case, Journal of Computational Physics (2024) 112881

    2024

  40. [49]

    R. F. Johnson, A. D. Kercher, The nondimensionalization of the chemically reacting Navier-Stokes equations, Tech. Rep. NRL/6041/MR–2022/4, U.S. Naval Research Laboratory (2022)

    2022

  41. [50]

    R. F. Johnson, A. D. Kercher, A conservative discontinuous Galerkin discretization for the chemically reacting Navier– Stokes equations, Journal of Computational Physics 423 (2020) 109826.doi:10.1016/j.jcp.2020.109826

    2020

  42. [51]

    Trojak, T

    W. Trojak, T. Dzanic, Positivity-preserving discontinuous spectral element methods for compressible multi-species flows, Computers & Fluids 280 (2024) 106343

    2024

  43. [52]

    P. C. Ma, Y. Lv, M. Ihme, An entropy-stable hybrid scheme for simulations of transcritical real-fluid flows, Journal of Computational Physics 340 (2017) 330–357

    2017

  44. [53]

    A. E. Magzumov, I. A. Kirillov, V. D. Rusanov, Effect of small additives of ozone and hydrogen peroxide on the induction- zone length of hydrogen-air mixtures in a one-dimensional model of a detonation wave, Combustion, Explosion and Shock Waves 34 (1998) 338–341

    1998

  45. [54]

    Crane, X

    J. Crane, X. Shi, A. V. Singh, Y. Tao, H. Wang, Isolating the effect of induction length on detonation structure: Hydrogen– oxygen detonation promoted by ozone, Combustion and Flame 200 (2019) 44–52

    2019

  46. [55]

    E. J. Ching, R. F. Johnson, Effect of ozone sensitization on the reflection patterns and stabilization of standing detonation waves induced by curved ramps, arXiv preprint arXiv:2410.13845https://arxiv.org/abs/2410.13845 (2024)

    2024 arXiv

  47. [56]

    Geuzaine, J.-F

    C. Geuzaine, J.-F. Remacle, Gmsh: A three-dimensional finite element mesh generator with built-in pre- and post- processing facilities, International Journal for Numerical Methods in Engineering (79(11)) (2009) 1310–1331. Appendix A. Compressible vortex transport: Grid converg...

    2009

  48. [623]

    doi:https://doi.org/10.1006/jcph.2000.6685

    2000

Pith tools

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