arxiv: 2605.06131 · v1 · submitted 2026-05-07 · ⚛️ physics.plasm-ph · math-ph· math.MP

Recognition: unknown

Hugoniot Relation for Multi-Temperature Euler Equations of Compressible Plasma Flows

Zhifang Du , Aleksey Sikstel

Authors on Pith no claims yet

Pith reviewed 2026-05-08 04:27 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph math-phmath.MP

keywords multi-temperature Euler equationsHugoniot relationshock wavesplasma flowsnon-conservative termsshock admissibilitydiscontinuous solutionscompressible flows

0 comments

The pith

Multi-temperature Euler equations for compressible plasma flows admit two distinct admissible Hugoniot relations rather than one.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that shock solutions in multi-temperature Euler models become ambiguous once microscopic physics is averaged out and non-conservative terms appear in the equations. It derives two separate Hugoniot relations, one for the general multi-temperature case and one for two-temperature plasmas, then verifies that both satisfy the classical Courant-Friedrichs admissibility conditions. A reader would care because this means the jump conditions across a shock cannot be read off from the macroscopic PDEs alone; they must be supplied by experiments or first-principles calculations, limiting what continuum models can predict about discontinuous plasma flows without additional microscopic input.

Core claim

The authors derive two distinct Hugoniot relations for the multi-temperature Euler equations—one for the general case and one specialized to two-temperature plasma flows. Classical analysis following Courant and Friedrichs shows that both relations correspond to physically admissible shock solutions even though the system contains non-conservative terms. The result establishes a fundamental non-uniqueness: the Hugoniot relation is not fixed by the macroscopic equations and must be provided from external sources such as experiments or first-principles simulations.

What carries the argument

The two analytically derived Hugoniot relations that serve as explicit jump conditions, each satisfying admissibility criteria for the non-conservative multi-temperature system.

If this is right

Numerical schemes for discontinuous plasma flows can use the two relations as reference solutions to test structure-preserving and vanishing-viscosity methods.
Any macroscopic model of compressible plasma shocks must include an external rule to select or justify which Hugoniot relation applies.
Microscopic physics is required to resolve the ambiguity in shock structure for models reduced from multi-temperature descriptions.
Admissibility analysis can be performed directly on non-conservative multi-temperature systems to confirm multiple candidate solutions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Similar ambiguity may appear in other reduced fluid models that discard microscopic detail, suggesting the need for hybrid closures that carry partial microscopic information.
Different plasma regimes, distinguished by relaxation rates, may naturally select one relation over the other in practice.
The finding points toward numerical methods that dynamically incorporate microscopic data to choose the appropriate jump conditions during a simulation.

Load-bearing premise

Both derived Hugoniot relations correspond to physically admissible shocks when the classical Courant-Friedrichs criteria are applied to the multi-temperature Euler system that includes non-conservative terms.

What would settle it

A first-principles simulation or laboratory measurement of a shock propagating in a multi-temperature plasma whose measured jump conditions match neither of the two derived relations.

Figures

Figures reproduced from arXiv: 2605.06131 by Aleksey Sikstel, Zhifang Du.

**Figure 1.** Figure 1: The Hugoniot surface parameterized by the vanishing viscosity Hugoniot relation and the segment-path Hugoniot curve (red line). fig:surface view at source ↗

**Figure 2.** Figure 2: Test A. The ion pressure (left) and the electron pressure (right). fig:A_p -0.3 -0.2 -0.1 0 0.1 0.2 0.3 1 1.5 2 2.5 3 3.5 4 Ion Pressure structure preserving scheme, numerical solution structure preserving scheme, exact solution vanishing viscosity scheme, numerical solution vanishing viscosity scheme, exact solution -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.5 1 1.5 2 2.5 3 3.5 Electron Pressure structure preserving … view at source ↗

**Figure 3.** Figure 3: Test B. The ion pressure (left) and the electron pressure (right). fig:B_p 5.3.2. Double rarefaction wave Riemann problem The last case involves two rarefaction waves. The initial data is (ρ, u, pi , pe) = ( (1, −2, 1 3 , 0.2), x < 0, (1, 2, 1 3 , 0.2), x > 0. (84) The thermodynamical parameters are γi = 1.4, γe = 5 3 , µi = 1, µe = 100. The final time is 0.1. Both structure preserving scheme and vanishing… view at source ↗

**Figure 4.** Figure 4: Test C. The ion pressure (left) and the electron pressure (right). fig:C_p -0.6 -0.4 -0.2 0 0.2 0.4 0.6 -0.1 -0.05 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 Ion Pressure structure preserving scheme, numerical solution structure preserving scheme, exact solution vanishing viscosity scheme, numerical solution vanishing viscosity scheme, exact solution -0.6 -0.4 -0.2 0 0.2 0.4 0.6 -0.05 0 0.05 0.1 0.15 0.2 Electron P… view at source ↗

**Figure 5.** Figure 5: Double rarefaction. The ion pressure (left) and the electron pressure (right). fig:D_p 6. Discussion sec:discussion The analysis in this paper shows that the ambiguity of shock solutions for multi-temperature Euler equations is not a numerical artifact, but a structural feature of the non-conservative model. The Rankine-Hugoniot relations for mass, momentum, and total energy determine only the Hugoniot sub… view at source ↗

read the original abstract

Shock solutions for multi-temperature Euler equations are inherently ambiguous due to the loss of microscopic physical detail during model reduction and occurrence of non-conservative terms. This paper presents a detailed analytical study of shock structures in such models. We derive two distinct Hugoniot relations, each corresponding to a physically admissible shock solution: one for the general multi-temperature case and one for two-temperature plasma flows. Through classical analysis \`a la Courant--Friedrichs, we demonstrate that both satisfy admissibility conditions, revealing a fundamental non-uniqueness in shock structures. By relating these solutions to existing numerical schemes, the structure preserving and vanishing viscosity approaches, we provide physically justified references for constructing and evaluating discontinuous numerical approximations. In particular, we emphasize that the Hugoniot relation is not uniquely determined by the macroscopic PDEs alone, but must be supplied from external sources such as experiments or first-principles simulations. This insight demonstrates the essential role of microscopic physics in resolving shock ambiguity and contributes to the theoretical foundation for modeling discontinuous plasma flows.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper derives two explicit Hugoniot relations for multi-temperature plasma Euler shocks that both pass classical admissibility, but the non-conservative terms raise questions about whether both can be admissible under one consistent definition.

read the letter

The core point is that the Hugoniot relation for these multi-temperature flows cannot be pinned down from the macroscopic equations alone and must come from outside, like experiments or micro-scale simulations. The authors construct two distinct relations—one for the general case and one for two-temperature plasmas—and show each satisfies the Courant-Friedrichs admissibility conditions via standard characteristic counting. They also connect the relations to structure-preserving and vanishing-viscosity numerical schemes, which gives modelers concrete references to check against when building discontinuous approximations.

Referee Report

1 major / 2 minor

Summary. The manuscript derives two distinct Hugoniot relations for shock solutions in the multi-temperature Euler equations of compressible plasma flows—one for the general multi-temperature case and one specialized to two-temperature plasmas. It applies classical Courant-Friedrichs analysis to show both satisfy admissibility conditions, concludes that shock structures are non-unique, and argues that the Hugoniot relation cannot be fixed by the macroscopic PDEs alone but must be supplied from external sources such as experiments or first-principles simulations. The work relates these relations to structure-preserving and vanishing-viscosity numerical schemes to provide references for discontinuous approximations.

Significance. If the derivations hold and the admissibility analysis is robust, the result would clarify an important source of ambiguity in modeling discontinuous multi-temperature plasma flows. It explicitly credits the necessity of microscopic physics input and supplies concrete references for evaluating numerical schemes, which could strengthen theoretical foundations in plasma hydrodynamics. The parameter-free character of the non-uniqueness conclusion (no free parameters introduced inside the model) is a positive feature.

major comments (1)

[Admissibility analysis (post-derivation of Hugoniot relations)] The admissibility section (following the derivation of the two Hugoniot relations): the claim that both relations satisfy classical Courant-Friedrichs characteristic counting does not address the path-dependence of Rankine-Hugoniot conditions in the presence of non-conservative products. Different regularizations (e.g., different viscosity matrices or paths) generally produce different jump relations; the counting itself can change with the chosen path. The manuscript must specify the implicit regularization underlying each derived relation and demonstrate that both remain admissible under a single fixed definition of the weak solution, rather than under separate regularizations.

minor comments (2)

[Governing equations] Notation for the multi-temperature variables and the non-conservative terms should be introduced with explicit definitions in the governing-equations section to avoid ambiguity when the two Hugoniot relations are stated.
[Abstract] The abstract states that the relations 'pass admissibility tests' but does not preview the concrete form of either relation; adding one-line expressions for the two Hugoniot relations would improve readability.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the insightful comment on the admissibility analysis in the presence of non-conservative products. We address this point below and will revise the manuscript to provide the requested clarifications.

read point-by-point responses

Referee: The admissibility section (following the derivation of the two Hugoniot relations): the claim that both relations satisfy classical Courant-Friedrichs characteristic counting does not address the path-dependence of Rankine-Hugoniot conditions in the presence of non-conservative products. Different regularizations (e.g., different viscosity matrices or paths) generally produce different jump relations; the counting itself can change with the chosen path. The manuscript must specify the implicit regularization underlying each derived relation and demonstrate that both remain admissible under a single fixed definition of the weak solution, rather than under separate regularizations.

Authors: We acknowledge the referee's concern regarding the path-dependence of Rankine-Hugoniot conditions for non-conservative systems. The two Hugoniot relations in the manuscript are derived under different physical assumptions that implicitly correspond to distinct regularizations of the underlying multi-fluid model. The general multi-temperature relation assumes decoupled temperature equations, consistent with a vanishing-viscosity regularization using a diagonal viscosity matrix. The two-temperature relation includes inter-species energy exchange, corresponding to a regularization that incorporates relaxation terms in the limit. In the revised version, we will explicitly specify these implicit regularizations in the admissibility section and relate them to the structure-preserving and vanishing-viscosity schemes mentioned in the paper. We will further show that both relations can be admissible under a single fixed weak solution definition by adopting a common regularization path, for instance the one from the vanishing viscosity limit with uniform viscosity scaling. Under this path, the Courant-Friedrichs characteristic counting confirms admissibility for shocks satisfying the respective Lax entropy conditions derived from each Hugoniot relation. This supports the manuscript's conclusion that the macroscopic equations alone do not uniquely determine the Hugoniot relation, which must instead be informed by external data. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation explicitly requires external input for uniqueness

full rationale

The paper derives two distinct Hugoniot relations from the multi-temperature Euler system (one general, one for two-temperature plasma) and applies classical Courant-Friedrichs characteristic analysis to confirm admissibility for both. It then states outright that the Hugoniot relation is not uniquely fixed by the macroscopic PDEs alone and must be supplied from experiments or first-principles simulations. No derivation step reduces by construction to a fitted parameter, self-definition, or self-citation chain; the non-uniqueness is presented as an inherent feature resolved externally rather than internally manufactured. The analysis remains self-contained against the stated PDEs and standard admissibility criteria without load-bearing reliance on prior author work.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard hyperbolic PDE theory and the existence of non-conservative terms in the multi-temperature reduction; no new free parameters, invented entities, or ad-hoc axioms are introduced in the abstract.

axioms (1)

domain assumption Multi-temperature Euler equations contain non-conservative terms arising from model reduction that render shock solutions ambiguous.
Stated in the abstract as the root cause of non-uniqueness.

pith-pipeline@v0.9.0 · 5479 in / 1400 out tokens · 63583 ms · 2026-05-08T04:27:03.079901+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

27 extracted references · 3 canonical work pages

[1]

1948 , doi=

Supersonic Flow and Shock Waves , author =. 1948 , doi=

1948
[2]

, title =

Hugoniot, H. , title =. Journal de l'. 1887 , pages =
[3]

, title =

Hugoniot, H. , title =. Journal de l'. 1889 , pages =
[4]

Splitting shock heating between ions and electrons in an ionized gas , author =. Comput. & Fluids , year =
[5]

Christophe Chalons and Frederic Coquel , journal =. THE
[6]

2006 , pages =

Christophe Chalons and Frederic Coquel , journal =. 2006 , pages =

2006
[7]

Structure-preserving operators for thermal-nonequilibrium hydrodynamics , author =. J. Comput. Phys. , volume =. 2018 , pages =

2018
[8]

Breil and S

J. Breil and S. Galera and P.-H. Maire , journal =. Multi-material. 2011 , pages =

2011
[9]

Hugoniot equation of state of rock materials under shock compression , author =. Phil. Trans. R. Soc. A , volume =. 2016 , pages =

2016
[10]

Hugoniot and sound velocity measurements of bismuth in the range of 11–70

Ye Tan and Yuying Yu and Chengda Dai and Ke Jin and Qingsong Wang and Jianbo Hu and Hua Tan , journal =. Hugoniot and sound velocity measurements of bismuth in the range of 11–70. 2013 , pages =

2013
[11]

Optical technique for determining rarefaction wave velocities at very high pressures , author =. Rev. Sci. Instrum. , volume =. 1982 , pages =

1982
[12]

Fowles, G. R. and Isbell, W. M. , title = ". J. of Appl. Phys. , volume =. 1965 , doi =

1965
[13]

Equations of State of Metals , author =

Shock-Wave Compressions of Twenty-Seven Metals. Equations of State of Metals , author =. Phys. Rev. , volume =. 1957 , doi =

1957
[14]

and Yarger, Frederick L

Christian, Russell H. and Yarger, Frederick L. , title = ". J. Chem. Phys. , volume =. 1955 , doi =

1955
[15]

McQueen, R. G. and Marsh, S. P. , title = ". J. of Appl. Phys. , volume =. 1960 , doi=

1960
[16]

Modelling and numerical approximation for the nonconservative bitemperature

Aregba. Modelling and numerical approximation for the nonconservative bitemperature. ESAIM: Mathematical Modelling and Numerical Analysis , pages =. 2018 , publisher =. doi:10.1051/m2an/2017007 , mrnumber =

work page doi:10.1051/m2an/2017007 2018
[17]

SIAM journal on numerical analysis , volume=

Time-implicit approximation of the multipressure gas dynamics equations in several space dimensions , author=. SIAM journal on numerical analysis , volume=. 2010 , publisher=

2010
[18]

Hantke, Maren and M. Baer-. ESAIM: Proceedings and Surveys , volume=. 2025 , publisher=

2025
[19]

Nonlinear projection methods for multi-entropies

Berthon, Christophe and Coquel, Fr. Nonlinear projection methods for multi-entropies. Mathematics of computation , volume=. 2007 , doi=

2007
[20]

International Journal on Finite Volumes , volume=

Numerical model of a compressible multi-fluid fluctuating flow , author=. International Journal on Finite Volumes , volume=
[21]

Communications in Mathematical Sciences , year=

About viscous approximations of the bitemperature Euler system , author=. Communications in Mathematical Sciences , year=
[22]

Journal de math

Definition and weak stability of nonconservative products , author =. Journal de math
[23]

Numerische Mathematik , volume=

Navier-Stokes equations with several independent pressure laws and explicit predictor-corrector schemes , author=. Numerische Mathematik , volume=. 2005 , publisher=

2005
[24]

Chalons, Christophe and Coquel, Fr. The. Journal of Hyperbolic Differential Equations , volume=. 2005 , publisher=

2005
[25]

2013 , address=

Riemann solvers and numerical methods for fluid dynamics: a practical introduction , author=. 2013 , address=

2013
[26]

First-principles and classical molecular dynamics simulation of shocked polymers , author =. Phys. Rev. B , volume =. 2010 , month =. doi:10.1103/PhysRevB.81.054103 , url =

work page doi:10.1103/physrevb.81.054103 2010
[27]

First-principles equations of state for simulations of shock waves in silicon , author =. Phys. Rev. B , volume =. 2001 , month =. doi:10.1103/PhysRevB.64.214107 , url =

work page doi:10.1103/physrevb.64.214107 2001