🔢 mathematics

Generalized solution and Weak-Strong uniqueness for a barotropic Euler-Riesz system

This paper establishes the existence of global-in-time dissipative solutions and proves their weak-strong uniqueness for the barotropic Euler-Riesz system on the torus by recasting the nonlocal Riesz force as a local stress tensor via the Caffarelli-Silvestre extension and adapting the relative energy method.

Original authors: Nilasis Chaudhuri

Published 2026-08-04

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Nilasis Chaudhuri

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

Imagine the universe as a giant, invisible dance floor where trillions of tiny particles are constantly bumping, swirling, and pushing against each other. Some of these particles, like gas in a star or electrons in a plasma, don't just bump; they also feel a mysterious, long-distance tug. It's like if everyone on the dance floor could feel a gentle pull or push from everyone else, even if they were on the opposite side of the room. This is the world of fluid dynamics, specifically the study of "compressible fluids"—stuff that can squish and stretch, like air or hot gas. Scientists use complex math equations to predict how these fluids move. Usually, they rely on two main forces: the pressure of the gas itself (like a spring pushing back) and the forces between particles.

For a long time, mathematicians have been great at predicting these dances when the forces are simple and local (only neighbors matter). But when the forces are "non-local"—meaning every particle feels the pull of every other particle across the entire room—the math gets incredibly messy. It's like trying to choreograph a dance where every dancer is connected to every other dancer by an invisible, stretchy rubber band. Sometimes, the math breaks down, and the equations stop giving a single, clear answer. This creates a puzzle: if the math allows for multiple different outcomes for the same starting point, how do we know which one is the "real" physical behavior? This paper dives into that exact puzzle, trying to prove that even when the math gets fuzzy, there is still only one true way the dance can unfold, provided the forces are pushing the particles apart (repulsive) rather than pulling them together.

The Paper's Big Idea: Taming the Ghostly Rubber Bands

This paper tackles a specific, tricky version of the fluid dance called the Euler–Riesz system. Imagine a crowd of people (the fluid) moving around a torus (which is just a fancy word for a donut-shaped space, like a video game map where if you walk off the right edge, you appear on the left). These people are pushing each other with a special kind of force called a Riesz kernel. Think of this force as a ghostly rubber band connecting every person to every other person. The strength of this band depends on how far apart they are.

The author is interested in the case where these rubber bands are repulsive—they push people apart. This is important because in the real world, this setup models things like charged plasmas (clouds of electrons) where particles naturally repel each other. When particles repel, the total energy of the system is "coercive," which is a fancy way of saying the system has a natural tendency to stay stable and not collapse into a singularity.

However, there's a catch. The math describing this "ghostly rubber band" force is non-local. In simple terms, to know how hard a particle is being pushed at one spot, you have to calculate the sum of pushes from every single other particle in the entire donut-shaped universe. This makes the equations incredibly hard to solve, especially when the fluid gets turbulent or develops sharp shocks (like a sonic boom).

The Magic Trick: The "Extra Dimension" Elevator

The paper's main breakthrough is a clever mathematical trick to turn this impossible, non-local problem into a manageable, local one. The author uses a tool called the Caffarelli–Silvestre extension.

Imagine the fluid is living on a 2D floor (the surface of the donut). The non-local force is a headache because it connects points that are far apart on this floor. The author's trick is to imagine an invisible elevator shaft (an extra dimension) going up from every point on the floor. They lift the problem into this 3D space (or 4D if you count time).

In this new, higher-dimensional world, the "ghostly rubber bands" disappear. Instead, the force becomes a simple, local stress tensor—a kind of pressure pushing on the walls of the elevator shaft. It's like turning a complicated global conversation where everyone shouts at everyone else into a simple, local conversation where you only talk to the person standing right next to you in the elevator. By solving the problem in this extra dimension and then looking at the "shadow" (the trace) it casts back down on the 2D floor, they can describe the non-local force using standard, local math.

The Main Result: One True Dance

The paper's central goal is to prove Weak-Strong Uniqueness. Here is the scenario:

Strong Solution: A perfect, smooth, classical solution to the equations. This is the "ideal" dance where everything is predictable and calm.
Weak (or Dissipative) Solution: A more general, "rough" solution. This allows for the math to get messy, like when the fluid develops sharp shocks or tiny, chaotic swirls that the standard equations can't quite pin down. These are called "defects" or "measure-valued solutions."

The big question is: If you start with the exact same initial conditions (the same starting positions and speeds for all particles), can the "rough" solution drift away from the "smooth" solution and become something completely different?

The paper proves that the answer is NO.

As long as a smooth, strong solution exists, every rough, generalized solution with the same starting data must be identical to it. They are the same dance. The "defects" (the chaotic noise) vanish completely.

This is a huge deal because it means that even if the math gets messy and allows for many theoretical possibilities, the physical reality (represented by the smooth solution) is the only one that actually happens. The "rough" solutions don't get to invent their own reality; they are forced to follow the smooth one.

How They Did It: The Relative Energy Scorecard

To prove this, the author invented a new way to measure the "distance" between the rough solution and the smooth one. They call this the Relative Energy.

Think of it like a scorecard in a video game.

Kinetic Energy: How fast the particles are moving.
Internal Energy: How much the gas is squished or heated.
Interaction Energy: The energy stored in those ghostly rubber bands (now calculated via the extra dimension).

The author created a formula that adds up the differences between the rough solution and the smooth solution for all three of these energies. They then showed that this "distance" (the Relative Energy) cannot grow. In fact, they proved that if the distance starts at zero (because the solutions have the same start), it must stay at zero forever.

The proof relies on a clever inequality (a mathematical rule that says "A is less than or equal to B"). They showed that the rate at which the "distance" grows is controlled by the "distance" itself. Since the distance starts at zero, and the growth rate depends on the current distance, the distance can never escape zero. It's like a ball sitting at the very bottom of a bowl; no matter how you nudge it, if the bowl is shaped just right, it will roll right back to the center.

What This Means for the Future

The paper establishes that for this specific type of fluid (repulsive, barotropic pressure) in 2D or 3D, the math is robust. Even when the fluid gets wild and the equations get fuzzy, the "smooth" reality is the only valid outcome.

The author is careful to note that this works for a specific range of parameters (the "order" of the Riesz kernel is between 0 and 2, and the pressure exponent is greater than 1). They also clarify that this result is a proof, not just a simulation or a guess. They have rigorously demonstrated that the "defects" (the mathematical noise) are controlled by the energy of the system and cannot take over.

In short, this paper takes a chaotic, high-dimensional problem involving invisible long-range forces and uses a clever "elevator" trick to show that, despite the chaos, the universe has a strict rule: if you start with a smooth dance, you will always end up with a smooth dance. The rough, messy alternatives are mathematically possible but physically impossible to sustain.

Technical Summary: Generalized Solution and Weak-Strong Uniqueness for a Barotropic Euler-Riesz System

Problem Statement
The paper investigates the barotropic Euler-Riesz system on the torus $\mathbb{T}^d$ ( $d=2,3$ ), describing a compressible fluid with density $\varrho$ and velocity $u$ subject to a nonlocal interaction force. The system is governed by the continuity equation, the momentum equation, and a fractional Poisson equation linking the potential $\Phi$ to the density:
$\begin{aligned} \partial_t \varrho + \text{div}_x(\varrho u) &= 0, \\ \partial_t(\varrho u) + \text{div}_x(\varrho u \otimes u) + \nabla_x p(\varrho) + \varrho \nabla_x \Phi &= 0, \\ (-\Delta)^s \Phi &= \varrho - \bar{\varrho}, \end{aligned}$
where $p(\varrho) = a\varrho^\gamma$ ( $\gamma > 1$ ) and the interaction kernel corresponds to the Riesz potential of order $\beta = 2s \in (0, 2)$ . The study focuses on the repulsive case, where the interaction energy is coercive.

The central difficulty addressed is the rigorous definition and analysis of generalized (dissipative) solutions for this system. In the classical Euler-Poisson case ( $s=1$ ), the nonlocal force $\varrho \nabla \Phi$ can be rewritten as the divergence of a local stress tensor via the identity $(\varrho - \bar{\varrho})\nabla \Phi = -\text{div}_x T_c[\nabla \Phi]$ . However, for fractional orders $s \in (0, 1)$ , no such pointwise identity exists on the torus; the force is genuinely nonlocal, and the product $\varrho \nabla \Phi$ cannot be defined solely through standard energy bounds, as $\nabla \Phi$ lies in a distributional space.

Methodology
To overcome the lack of a local stress representation for fractional orders, the author employs the Caffarelli–Silvestre extension. This technique lifts the problem from the physical domain $\mathbb{T}^d$ to a half-cylinder $\Omega = \mathbb{T}^d \times (0, \infty)$ by introducing an auxiliary variable $z$ . The potential $\Phi$ is realized as the trace of an extended function $U(t, x, z)$ , which satisfies a local, degenerate elliptic equation:
$\text{div}_{x,z}(z^{1-2s} \nabla_{x,z} U) = 0 \quad \text{in } \Omega,$
with a Neumann-type boundary condition at $z=0$ recovering the fractional Laplacian.

This localization allows the nonlocal force to be expressed as the trace of a local stress tensor $T_{ER}[U]$ in the extended space. Specifically, the force term $\varrho \nabla \Phi$ is replaced by an integral involving $T_{ER}[U] = z^{1-2s}(\nabla_x U \otimes \nabla_x U - \frac{1}{2}|\nabla_{x,z} U|^2 \text{Id})$ .

Building on this structure, the paper introduces a notion of dissipative solutions. Unlike classical weak solutions, these are defined as a septuplet $(\varrho, m, U, R_{\text{fluid}}, R_{\text{ext}}, E_{\text{fluid}}, E_{\text{ext}})$ , where:

$(\varrho, m, U)$ are the fluid density, momentum, and extended potential.
$R_{\text{fluid}}$ and $R_{\text{ext}}$ are defect measures (turbulent stresses) arising from the weak limits of nonlinear terms in the fluid and extension equations, respectively.
$E_{\text{fluid}}$ and $E_{\text{ext}}$ are energy defect measures accounting for the loss of energy in the limit.

Crucially, the definition includes compatibility conditions linking the momentum defects to the energy defects. For the extension part, the defect measure $R_{\text{ext}}$ is not sign-definite but is dominated by the energy defect $E_{\text{ext}}$ via a Radon–Nikodým density with operator norm at most 1.

Key Contributions and Results

Definition of Dissipative Solutions: The paper establishes a consistent framework for generalized solutions for the Euler-Riesz system with $s \in (0, 1)$ . It proves that this class is sequentially stable, meaning that limits of consistent approximating sequences (e.g., vanishing viscosity limits) yield dissipative solutions. The defect measures are shown to arise naturally as the weak limits of the nonlinear terms.
Compatibility with Classical Solutions: It is proven that any dissipative solution that possesses sufficient regularity (specifically, if the density is strictly positive and fields are $C^1$ ) coincides with a classical solution, and all defect measures vanish. This ensures the definition does not introduce spurious solutions in the smooth regime.
Weak-Strong Uniqueness: The main result (Theorem 2) establishes that if a strong solution $(r, v)$ exists on a time interval $[0, T]$ , then any dissipative solution with the same initial data coincides with the strong solution on that interval. Consequently, all defect measures vanish, and the energy is conserved.
- The proof utilizes the relative energy method. A relative energy functional $E_{\text{rel}}$ is constructed, combining the kinetic/internal energy difference and the extension energy difference.
- A key technical achievement is the derivation of a relative energy inequality where the nonlocal force term is absorbed into the relative stress tensor of the extension. This allows the force to be controlled algebraically by the relative energy, without requiring fractional interpolation inequalities that would couple the pressure exponent $\gamma$ and the singularity order $\beta$ .

Significance and Claims
The paper claims that its primary contribution is filling a structural gap in the theory of nonlocal fluid equations: the lack of a local divergence structure for fractional Riesz interactions. By utilizing the Caffarelli–Silvestre extension, the author provides a mechanism to control the singular nonlocal force through a local stress tensor in an extended dimension.

A significant aspect of the result is the decoupling of the pressure exponent $\gamma$ and the interaction order $\beta$ . Previous works on weak-strong uniqueness for Euler-Riesz systems (e.g., Alves et al.) required a coupling condition such as $\gamma \geq 2 - \beta/d$ . In contrast, this paper establishes weak-strong uniqueness for every $\gamma > 1$ and every $\beta \in (0, 2)$ independently. The author notes that the restriction $\beta \in (0, 2)$ is inherent to the first-order nature of the Caffarelli–Silvestre extension used, while the independence from $\gamma$ is a direct consequence of controlling the nonlocal force via the extension energy rather than the pressure.

The work extends the theory of generalized solutions (initiated for compressible Euler by DiPerna, Majda, and later Breit, Feireisl, and Hofmanová) to systems with singular nonlocal interactions, providing a robust framework for analyzing large-data solutions and their stability.

The Paper's Big Idea: Taming the Ghostly Rubber Bands

The Magic Trick: The "Extra Dimension" Elevator

The Main Result: One True Dance

How They Did It: The Relative Energy Scorecard

What This Means for the Future

More like this