🔢 mathematics

Entropy stable finite difference schemes for One-Fluid Two-Temperature Euler Non-equilibrium Hydrodynamics

This paper presents the design and validation of higher-order entropy-stable finite difference schemes for the One-Fluid Two-Temperature Euler equations by reformulating the non-conservative terms to ensure zero entropy contribution, utilizing Tadmor's relation for conservative parts, and incorporating entropy-dissipation terms based on scaled eigenvectors.

Original authors: Chetan Singh, Harish Kumar

Published 2026-05-18

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Chetan Singh, Harish Kumar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Simulating a "Two-Temperature" Fluid

Imagine you are trying to simulate a super-hot gas, like the plasma inside a star or a fusion reactor. In normal air, the molecules bounce around and quickly share energy, so everything is at the same temperature. But in this specific type of gas (plasma), the heavy particles (ions) and the light particles (electrons) are like two different teams playing on the same field. They don't share their heat very well. The electrons might be scorching hot while the ions are just warm, or vice versa.

This paper is about building a computer simulation to track how this "two-temperature" gas moves and changes over time. The authors call this the "One-Fluid Two-Temperature Euler" model.

The Problem: The "Unruly" Math

The math describing this gas is tricky. It's like trying to drive a car where the steering wheel and the gas pedal are connected in a weird, non-standard way.

Non-Conservative Terms: In physics, "conservation" usually means things like mass or energy can't just disappear. But in this specific model, the math includes terms that don't fit the standard "conservation" rules perfectly. It's like a recipe that says "add some flour" but doesn't tell you how much unless you guess the path you took to get there. This makes it very hard to write a computer program that doesn't crash or give nonsense results.
Entropy (The Mess Factor): In physics, "entropy" is a measure of disorder or "messiness." The Second Law of Thermodynamics says that in a closed system, things naturally get messier (entropy goes up) or stay the same, but they never get perfectly organized on their own.
- The Risk: If a computer simulation is poorly designed, it might accidentally make the system less messy (decrease entropy) in a way that violates the laws of physics. This causes the simulation to blow up, creating wild, impossible numbers (like negative pressure) that ruin the whole calculation.

The Solution: Building a "Stable" Simulation

The authors wanted to build a computer code (a numerical scheme) that is Entropy Stable. Think of this as building a car with a "fail-safe" engine. No matter how rough the road gets (shocks, explosions, sudden changes), the car will never violate the laws of physics. It will always respect the rule that "things get messier, not cleaner."

Here is how they did it, step-by-step:

1. Rewriting the Rules (Reformulation)

First, they realized the original math was too messy to work with directly. So, they rewrote the equations.

The Analogy: Imagine you have a complex knot of ropes. It's hard to pull them apart. The authors found a way to untie the knot and rearrange the ropes so that the "messy" part (the non-conservative term) is moved to a place where it doesn't affect the "messiness" (entropy) of the system.
The Result: They created a new version of the equations where the "messy" part is neutralized. It's like putting a guard on a door so that the chaotic crowd outside can't mess up the order inside the room.

2. The Perfect Balance (Entropy Conservative)

Next, they built a basic version of the simulation that is perfectly balanced.

The Analogy: Imagine a perfectly balanced scale. If you add weight to one side, you add the exact same amount to the other. In their simulation, they designed a "flux" (a way of moving data between grid points) that ensures the total "messiness" stays exactly the same if there are no shocks. It's a frictionless, perfect system.
The Trick: They used a specific mathematical recipe (Tadmor's relation) to ensure that when they calculated the flow of gas, they didn't accidentally create or destroy entropy.

3. Adding the Safety Net (Entropy Dissipation)

Real life isn't perfect. When gas hits a wall or creates a shockwave, it should get messier (entropy should increase). A perfectly balanced simulation would just keep the same level of mess, which isn't realistic for shocks.

The Analogy: Imagine driving a car on a perfect track. You need brakes for when you hit a sudden obstacle. The authors added a "braking system" (dissipation terms) to their simulation.
How it works: This braking system is smart. It only kicks in when the simulation detects a shock or a sudden jump in the data. It adds just enough "friction" to ensure the entropy goes up (or stays the same), but never down. This prevents the simulation from crashing.

The Results: Testing the Engine

The authors tested their new "car" (the simulation) in several ways:

Smooth Roads (Accuracy Tests): They ran the simulation on smooth, predictable flows. They checked if the results matched the known math.
- Result: They built versions that were 2nd, 3rd, and 4th order accurate. Think of this as the resolution of a camera. The 4th-order version is like a 4K camera, capturing tiny details that the 2nd-order version (standard definition) misses.
Rough Roads (Shock Tests): They simulated violent events, like:
- Double Rarefaction: Gas expanding rapidly in two directions.
- Sod and Lax Problems: Classic crash tests for gas simulations involving shockwaves.
- Shock-Bubble Interaction: A shockwave hitting a bubble of gas, squishing it and making it spin.
- Richtmyer-Meshkov Instability: A complex mixing of fluids caused by a shockwave, creating intricate, swirling patterns (like ink dropping in water).
The Outcome:
- The 2nd-order version worked but was "blurry." It smoothed out the details too much (like a low-resolution photo).
- The 3rd and 4th-order versions were sharp and clear. They captured the tiny swirls and sharp edges of the shockwaves perfectly.
- Crucially: In every single test, the simulation respected the entropy rules. It never crashed, and it never produced impossible physics. The "messiness" always behaved correctly.

Summary

The paper presents a new, robust way to simulate hot, two-temperature gases. The authors took a difficult, unstable set of equations, rewrote them to separate the "tricky" parts, and built a computer code that is mathematically guaranteed to respect the laws of thermodynamics. They proved that their high-resolution (4th-order) version can capture complex, violent fluid interactions with high accuracy without breaking the physics.

Technical Summary: Entropy Stable Finite Difference Schemes for One-Fluid Two-Temperature Euler Non-equilibrium Hydrodynamics

Problem Statement
The paper addresses the numerical simulation of the One-Fluid Two-Temperature Euler (OFTT-Euler) system, a model used for non-equilibrium hydrodynamics in fields such as astrophysics, inertial confinement fusion (ICF), and high-energy-density physics. This system describes radiative, high-energy-density plasmas containing ions and electrons, where each species maintains its own thermal equilibrium but possesses a distinct temperature.

Mathematically, the OFTT-Euler system is a set of nonlinear hyperbolic partial differential equations (PDEs) containing non-conservative products. The total pressure is decomposed into electron ( $p_e$ ) and ion ( $p_i$ ) components, leading to separate evolution equations for their internal energies. The presence of non-conservative terms (specifically involving $\nabla \cdot \mathbf{v}$ ) presents significant challenges for numerical stability. Standard methods often struggle to define weak admissible solutions without specifying a path, and numerical viscosity can adversely affect the solution. Furthermore, while entropy stability is a crucial theoretical nonlinear stability estimate for hyperbolic systems, designing entropy-stable schemes for non-conservative systems is particularly difficult.

Methodology
The authors propose a framework for designing higher-order entropy-stable finite difference schemes for the OFTT-Euler system. The methodology proceeds in three main stages:

Reformulation of the System:
Following the approach of [41, 42], the authors reformulate the OFTT-Euler equations. The goal is to modify the system such that the non-conservative terms do not contribute to the entropy evolution. By rewriting the fluxes and source terms, they derive a new system where the non-conservative matrices $C_d(U)$ satisfy the condition $\mathbf{w}^T C_d(U) = 0$ , where $\mathbf{w}$ represents the entropy variables. This reformulation ensures that the entropy evolution is governed solely by the conservative part of the system.
Entropy-Conservative Fluxes:
For the reformulated conservative part, the authors derive second-order entropy-conservative numerical fluxes using Tadmor's relation. These fluxes satisfy the condition $[\![\mathbf{w}]\!] \cdot \tilde{\mathbf{f}} = [\![\Psi]\!]$ , where $\Psi$ is the entropy potential. These second-order fluxes are then extended to higher orders (specifically 3rd and 4th order) using a linear combination of second-order fluxes, following the technique described in [48].
Entropy-Stable Dissipation:
To ensure stability in the presence of discontinuities (shocks), the authors add dissipation terms to the entropy-conservative fluxes. They construct a higher-order entropy diffusion operator using entropy-scaled right eigenvectors of the conservative part. The dissipation term is designed to be symmetric positive definite, based on Rusanov-type diffusion operators. To maintain high-order accuracy while preserving entropy stability, the authors employ a reconstruction procedure (using MinMod for 2nd-order and ENO-based reconstruction for higher orders) on scaled entropy variables. This ensures the sign-preserving property required for the stability of the numerical flux.

The time discretization is handled using explicit Strong Stability-Preserving (SSP) Runge-Kutta methods (SSP-RK) of orders 2, 3, and 4.

Key Contributions

Novel Reformulation: The paper presents a specific reformulation of the OFTT-Euler system that isolates the non-conservative terms such that they do not affect the entropy balance, a necessary step for applying entropy stability theory to this non-symmetrizable system.
Higher-Order Schemes: The authors successfully derive and implement 2nd, 3rd, and 4th-order entropy-stable finite difference schemes for this specific non-conservative hyperbolic system.
Rigorous Stability Proof: The work provides a theoretical proof that the proposed semi-discrete and fully discrete schemes satisfy the discrete entropy inequality, ensuring that the total entropy does not increase over time.

Numerical Results
The proposed schemes (denoted as O2es, O3es, and O4es) were tested on a variety of one-dimensional and two-dimensional problems:

Accuracy Tests: Smooth problems in 1D and 2D demonstrated that the schemes achieve their theoretically predicted orders of accuracy (2nd, 3rd, and 4th) for density and pressure variables.
Riemann Problems: Several 1D Riemann problems (Double rarefaction, Sod, and Lax) were simulated. The results showed that while all schemes captured the waves, the higher-order schemes (O3es and O4es) were less diffusive and produced sharper resolution of shocks and contact discontinuities compared to the 2nd-order scheme.
2D Riemann Problem: A complex 2D Riemann problem involving interacting waves showed that the higher-order schemes captured multi-scale structures significantly better than the 2nd-order scheme.
Shock-Bubble Interaction and Richtmyer-Meshkov Instability: These tests, relevant to ICF applications, highlighted the ability of the higher-order schemes to resolve small-scale features and vortex rings that were smeared out by the lower-order, more diffusive scheme.
Entropy Stability: In all test cases, the total entropy decay was monitored. The results confirmed that the schemes are entropy stable, with the expected entropy production at discontinuities and negligible decay in smooth regions.

Significance and Claims
The paper claims to provide the first entropy-stable finite difference schemes for the OFTT-Euler system, addressing the challenge posed by its non-conservative nature. The authors emphasize that their reformulation allows the system to be treated in a way that makes entropy stability achievable. The numerical results serve to demonstrate that the proposed schemes are not only theoretically sound but also practically accurate and stable for complex non-equilibrium hydrodynamic flows. The work establishes a foundation for robust, high-order simulations of two-temperature plasmas without the need for ad-hoc numerical viscosity that might violate physical entropy conditions.