🔢 mathematics

Steady compressible Navier-Stokes-Fourier system with general temperature dependent viscosities and hard sphere pressure law

This paper establishes the existence of solutions for the steady compressible Navier-Stokes-Fourier system in a three-dimensional bounded domain featuring temperature-dependent viscosities of the form $(1+\vartheta)^\alpha$ and a singular hard sphere pressure law, utilizing Bogovskii operator estimates to control the pressure under various boundary conditions.

Original authors: Zhengguang Guo, Milan Pokorný

Published 2026-06-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Zhengguang Guo, Milan Pokorný

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 a fluid, like air or water, flowing steadily inside a closed, three-dimensional room. This fluid isn't just moving; it's also getting hot and cold, changing its density, and pushing against the walls. The paper you are asking about is a mathematical "recipe" that proves such a flow can actually exist under very specific, realistic, and complicated conditions.

Here is a breakdown of what the authors, Zhengguang Guo and Milan Pokorný, achieved, using everyday analogies.

The Big Picture: The "Perfectly Realistic" Fluid

Most simple math models for fluids make life easy by assuming the fluid's "stickiness" (viscosity) and how well it conducts heat stay the same, no matter how hot it gets. But in the real world, honey flows differently when it's cold versus when it's hot.

This paper tackles the hard version of the problem:

Temperature-Dependent Stickiness: The fluid gets stickier or thinner depending on its temperature (specifically, following a rule like $(1 + \text{temperature})^\alpha$ ).
The "Hard Sphere" Pressure: Imagine the fluid is made of tiny, hard marbles. As you squeeze them together, they resist fiercely. If you try to pack them too tightly (reaching a specific density limit), the pressure shoots up to infinity, like trying to force a solid rock into a space that's already full. The paper uses a mathematical model for this "hard sphere" behavior.
The Room: The fluid is in a 3D box with specific rules for how it touches the walls (either sliding, sticking, or having heat flow in/out).

The Challenge: The "Pressure Puzzle"

The main difficulty in solving these equations is the pressure. In fluid dynamics, pressure is the "glue" that holds the solution together. If you can't prove the pressure behaves nicely, the whole mathematical solution falls apart.

The authors faced a tricky situation:

The fluid's density is capped (it can't get denser than the "hard spheres" allow).
The temperature can get very hot or very cold.
The "stickiness" changes with temperature.

To solve this, they had to prove that even with all these changing variables, the pressure doesn't go crazy. They did this using a mathematical tool called the Bogovskii operator.

The Analogy: Think of the Bogovskii operator as a specialized "pressure regulator" or a "traffic cop" for the math. It helps the authors take a messy, uneven distribution of fluid density and smooth it out just enough to prove that the pressure stays within a manageable range, preventing the math from breaking down.

The Four Scenarios (The "Problems")

The paper doesn't just solve one case; it solves four different "flavors" of the problem, depending on how the fluid interacts with the walls of the room:

Problem 1 (The Quiet Room): The fluid is stuck to the walls (no sliding), and heat flows out through the walls like a radiator cooling down a room.
Problem 2 (The Fixed Temperature Room): The fluid is stuck to the walls, but the walls themselves are kept at a specific, fixed temperature (like a thermostat).
Problem 3 (The Moving Wall): The walls are moving (pushing the fluid), but the fluid doesn't leak in or out. Heat flows out like a radiator.
Problem 4 (The Moving, Thermostated Wall): The walls are moving, and they are also kept at a fixed temperature. This is the most complex scenario.

The "Secret Sauce": Energy Balances

To prove the solution exists, the authors had to track the "energy" of the system. They used two main concepts:

Variational Entropy: Think of this as a "disorder meter." The authors proved that even though the fluid is chaotic, the total "disorder" (entropy) follows a predictable rule that prevents the system from exploding mathematically.
Ballistic Energy: This is a special accounting trick for the energy. It's like a "budget" that ensures the energy entering the system (from external forces or heat) balances out with the energy leaving or being used up by friction.

The Verdict

The authors successfully proved that solutions exist for all four scenarios, provided the "stickiness" and "heat conduction" follow certain mathematical rules (specifically, how fast they grow as the temperature rises).

Key Finding: They found the "Goldilocks zone" for the parameters. If the heat conduction grows too slowly or the stickiness changes too wildly, the math breaks. But if they stay within their calculated limits, a stable, steady flow is guaranteed to exist.

What They Did Not Do

It is important to note what this paper is not:

It does not simulate a specific engine or weather pattern.
It does not provide a computer code to run these simulations.
It does not suggest new medical treatments or industrial applications.

Instead, it is a theoretical foundation. It is the mathematical equivalent of an architect proving that a bridge can be built with these specific materials and forces before any construction crew ever picks up a hammer. It says, "Yes, the laws of physics allow this complex flow to exist; here is the proof."

Technical Summary: Steady Compressible Navier–Stokes–Fourier System with General Temperature-Dependent Viscosities and Hard Sphere Pressure Law

Problem Formulation
This paper investigates the existence theory for the steady compressible Navier–Stokes–Fourier system in a three-dimensional bounded domain $\Omega \subset \mathbb{R}^3$ . The system models the steady flow of a heat-conducting compressible fluid and is governed by the conservation of mass, momentum, and total energy. The study focuses on a specific equation of state known as the "hard sphere pressure law," defined as $p(\varrho, \vartheta) = \vartheta \varrho h(\varrho)$ , where $h(\varrho)$ is a non-decreasing function that becomes singular at a critical density $\varrho = a > 0$ . This pressure law formally corresponds to an adiabatic exponent $\gamma = \infty$ and vanishes as temperature $\vartheta \to 0$ .

The viscosity coefficients are assumed to depend on the temperature $\vartheta$ in the form $\mu(\vartheta), \xi(\vartheta) \sim (1+\vartheta)^\alpha$ with $0 \le \alpha \le 1$ . The heat conductivity follows $\kappa(\vartheta) \sim (1+\vartheta)^m$ for $m > 0$ . The paper considers four distinct boundary value problems involving combinations of:

Velocity: Homogeneous Dirichlet ( $u=0$ ) or non-homogeneous Dirichlet ( $u=u_0$ with $u_0 \cdot n = 0$ ).
Temperature: Heat flux ( $q \cdot n = L(\vartheta - \vartheta_0)$ ) or Dirichlet ( $\vartheta = \vartheta_D$ ).

Additionally, the total mass of the fluid is prescribed, and the domain must satisfy a geometric constraint ensuring the density remains below the singularity $a$ .

Methodology
The authors employ a variational approach to establish the existence of weak and variational entropy solutions. The methodology proceeds through the following steps:

Approximation Scheme: An approximate system is constructed using regularization parameters ( $\varepsilon, \delta, \eta, R$ ). The density is truncated and regularized, and the pressure is modified to avoid the singularity at $\varrho = a$ during the approximation phase. The velocity is approximated using a Galerkin method in a finite-dimensional subspace.
A Priori Estimates: The core of the analysis involves deriving uniform estimates independent of the regularization parameters.
- Energy and Entropy: The authors utilize the total energy balance and the entropy inequality (or ballistic energy inequality for Dirichlet temperature conditions) to control the temperature and velocity fields. Special attention is paid to the behavior of the specific entropy near absolute zero, requiring a modified internal energy function $g(\vartheta)$ for Dirichlet boundary conditions.
- Pressure Control: A critical challenge is controlling the pressure in a space better than $L^1$ . The authors utilize the Bogovskii operator to test the momentum equation. By testing with $\Phi = B(p^\beta)$ , they derive estimates for the pressure in $L^{1+\beta}$ . This requires careful balancing of the exponents $\alpha$ and $m$ to ensure the integrability of nonlinear terms involving viscosity and temperature.
- Bogovskii Operator Estimates: The control of the pressure relies heavily on estimates involving the Bogovskii operator applied to the density and pressure terms, leveraging the boundedness of the density ( $\varrho < a$ ) and the specific structure of the hard sphere pressure.
Limit Passage: The authors pass to the limit as the regularization parameters tend to zero. This involves:
- Establishing the strong convergence of the density sequence, which is achieved by proving the effective viscous flux identity and utilizing the renormalized continuity equation.
- Handling the limit in the total energy balance, which imposes stricter conditions on the parameters $m$ and $\alpha$ compared to the momentum equation.

Key Contributions and Results
The paper establishes the existence of solutions for the four boundary value problems described above under specific constraints on the parameters $\alpha$ and $m$ . The main results are summarized as follows:

Theorem 1 (Problem 1: $u=0$ , Heat Flux): A variational entropy solution exists if $\alpha \in [0, 1]$ and $m > \max\{\frac{1+\alpha}{3}, \frac{1-\alpha}{2}\}$ . If $m > 1$ , the solution is a weak solution.
Theorem 2 (Problem 2: $u=0$ , Dirichlet Temperature): A variational entropy solution exists if $\alpha \in [0, 1]$ and $m > \max\{\frac{1+\alpha}{3}, 1-\alpha\}$ . If $m > 1$ , the solution is weak.
Theorem 3 (Problem 3: $u=u_0 \neq 0$ , Heat Flux):
- For $0 < \alpha < 1$ , existence holds if $m > \max\{\frac{1+\alpha}{3}, \frac{1-\alpha}{2}\}$ .
- For $\alpha = 0$ , existence of a weak solution requires $m > 1$ .
- For $\alpha = 1$ , existence requires $\|\nabla u_0\|_\infty$ to be sufficiently small and $m \ge \frac{2}{3}$ . If $m > 1$ , the solution is weak.
Theorem 4 (Problem 4: $u=u_0 \neq 0$ , Dirichlet Temperature): A weak entropy solution exists if $\alpha = 1$ , $m \ge 2$ , and $\|\nabla u_0\|_\infty$ is sufficiently small.

The paper distinguishes between "variational entropy solutions" (satisfying the entropy inequality and energy balance in a specific weak form) and "weak solutions" (satisfying the energy balance in the standard weak form). The transition from variational to weak solutions generally requires $m > 1$ to ensure sufficient integrability for the limit passage in the energy equation.

Significance
The paper extends the existing theory of steady compressible Navier–Stokes–Fourier systems in three dimensions. While previous works have addressed constant viscosities or specific power-law pressures (e.g., $p \sim \varrho^\gamma + \vartheta \varrho$ ), this study addresses the hard sphere pressure law, which is physically relevant for gases at high densities and formally corresponds to $\gamma = \infty$ .

The authors claim that their results represent an extension to the three-dimensional situation for this specific pressure law, which was previously only analyzed in two dimensions (where Trudinger-Moser inequalities played a key role). Furthermore, the paper generalizes the treatment of temperature-dependent viscosities and heat conductivity, providing existence results for a broad range of exponents $\alpha$ and $m$ under various boundary conditions, including the more challenging non-homogeneous velocity and Dirichlet temperature cases. The work relies on and refines techniques involving the Bogovskii operator and renormalized solutions to handle the singular nature of the pressure and the nonlinearity of the transport coefficients.