Navier-Stokes-Fourier Equations

• The Navier-Stokes-Fourier system is a mathematical model for compressible, viscous, and heat-conducting fluids governed by conservation laws.
• It couples mass, momentum, and energy equations with constitutive laws for pressure, stress, and heat flux to ensure thermodynamic consistency.
• Recent advances analyze the stability of strong solutions and extend measure-valued formulations to quantify uncertainty and address blow-up phenomena.

The Navier-Stokes-Fourier (NSF) system describes the evolution of a compressible, viscous, heat-conducting fluid under classical thermomechanical laws. The system couples conservation of mass, linear momentum, and internal energy via constitutive relations for pressure, viscous stress, and heat flux, typically modeled as Newtonian and Fourier laws. The framework is foundational in continuum thermomechanics, underpinning rigorous analysis of stability, statistical solutions, and uncertainty quantification. In particular, recent advances have established stability of strong solutions, robust pushforward representation of measure-valued (statistical) solutions, and analytic tools for extending solution concepts beyond blow-up time in the sense of Feireisl–Lukáčová-Medviďová (Feireisl et al., 2022). The NSF system serves as a canonical model for compressible flows subject to thermal effects, supporting deep connections to stochastic analysis, statistical mechanics, and numerical approximation theories.

1. Mathematical Formulation and Constitutive Laws

The NSF system is posed on a three-dimensional torus $\Omega = T^3$ with periodic boundary conditions. Primary fields are the density $\rho(t, x) > 0$, absolute temperature $\vartheta(t, x) > 0$, and velocity $u(t, x) \in \mathbb{R}^3$. The PDEs governing the system comprise:

• Continuity equation (mass conservation):

$\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$

• Momentum balance (Newton's law):

$$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$

• Internal energy balance (first law, with heat conduction):

$$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$

where $g(t,x)$ is a given volume force, $Q(t, x)$ an external heating, and the closure relations for a Newtonian caloric ideal gas are:

• Equation of state: $p(\rho, \vartheta) = \rho \vartheta$, $\rho(t, x) > 0$0 with $\rho(t, x) > 0$1.
• Viscous stress: $\rho(t, x) > 0$2 with shear viscosity $\rho(t, x) > 0$3, bulk viscosity $\rho(t, x) > 0$4.
• Heat conduction: $\rho(t, x) > 0$5 with $\rho(t, x) > 0$6.

Admissibility of initial data $\rho(t, x) > 0$7 requires $\rho(t, x) > 0$8, $\rho(t, x) > 0$9, $\vartheta(t, x) > 0$0, $\vartheta(t, x) > 0$1, $\vartheta(t, x) > 0$2 (Feireisl et al., 2022).

2. Well-Posedness and Strong Solution Theory

Let $\vartheta(t, x) > 0$3 and trajectory space $\vartheta(t, x) > 0$4 built from regularity classes of initial data and time-evolved fields. A triple $\vartheta(t, x) > 0$5 is a strong solution on $\vartheta(t, x) > 0$6 if:

• $\vartheta(t, x) > 0$7, $\vartheta(t, x) > 0$8, $\vartheta(t, x) > 0$9.
• $u(t, x) \in \mathbb{R}^3$0, $u(t, x) \in \mathbb{R}^3$1.
• $u(t, x) \in \mathbb{R}^3$2, $u(t, x) \in \mathbb{R}^3$3, $u(t, x) \in \mathbb{R}^3$4.

The classical result (cf. Cho–Kim) assures existence of a unique strong solution up to a maximal time $u(t, x) \in \mathbb{R}^3$5, with blow-up detected by divergence of relevant Sobolev norms:

$u(t, x) \in \mathbb{R}^3$6

(Feireisl et al., 2022)

3. Stability of Solutions and Data Dependence

A central result demonstrates stability of strong solutions with respect to data perturbations. For sequences $u(t, x) \in \mathbb{R}^3$7 in the product space $u(t, x) \in \mathbb{R}^3$8:

• $u(t, x) \in \mathbb{R}^3$9.
• For every $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$0, solution trajectories $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$1 weak-* in $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$2.

This yields convergence in $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$3 on any subinterval $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$4 prior to possible singularity formation (Feireisl et al., 2022). The proof exploits uniform energy-type bounds and Sobolev compactness structures.

4. Statistical Solutions and Measure-Valued Framework

Statistical solutions arise by randomization of data $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$5, distributed according to a Borel probability measure $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$6 on $\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$7. Denote

$\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$8

and define the time-$\partial_t \rho + \nabla_x \cdot (\rho u) = 0.$9 marginal as push-forward $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$0. The statistical solution is then the family $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$1 of probability measures on phase space $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$2, satisfying a Markov semigroup property:

$$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$3

Crucial properties:

• $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$4 puts zero mass on blow-up events for $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$5, by lower semicontinuity of $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$6.
• The semigroup structure persists in the autonomous case (time-independent $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$7).
• Disintegration of $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$8 with respect to physical parameters yields decomposable family of statistical solutions.

Push-forward construction avoids further selection arguments—Polish metric topology and continuity of solution operator suffice for compactness and tightness (Feireisl et al., 2022).

5. Extension Beyond Blow-Up and the Statistical Measure Approach

The analytic machinery supports extension of individual strong solutions beyond blow-up time by setting $$\partial_t (\rho u) + \nabla_x \cdot (\rho u \otimes u) + \nabla_x p(\rho, \vartheta) = \nabla_x \cdot \mathbb{S}(\nabla_x u) + \rho g.$$9 for $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$0, and augmenting the phase space $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$1 with a zero element. The endowment of a Polish metric $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$2 enforces $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$3 as $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$4. Solution flow $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$5 remains continuous in this topology. This enables rigorous definition of statistical solutions "beyond blow-up", applicable to uncertainty quantification and numerical schemes (Feireisl et al., 2022).

Monte Carlo methods: Empirical distributions from Monte Carlo samples of random $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$6 converge almost surely in the sense of the Banach-valued strong law (Ledoux–Talagrand) to the statistical solution $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$7, linking theory with probabilistic numerical methods.

6. Analytical and Probabilistic Implications

This framework establishes a mathematically sound basis for statistical solutions of compressible NSF systems—measure-valued, consistent for random or uncertain data, supporting robust propagation up to (and beyond) singularity formation. The tightness and continuity properties facilitate rigorous uncertainty quantification and reliability analysis in physical or engineering simulations. The methods also provide a template for extending similar constructions to more complex multi-physics PDE systems, stochastic and measure-valued formulations, and for connecting deterministic and probabilistic solution concepts in nonlinear thermomechanics. The statistical approach is compatible with empirical algorithms, and avoids delicate selection issues present in other measure-valued PDE theories.

7. Conceptual Table: Solution Structures and Extensions

Solution Type	Definition	Extension Beyond Blow-Up
Strong solution	Pointwise, regular in Sobolev spaces	Not defined beyond $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$8
Statistical solution	Pushforward of data measure via solution operator	Set $$\partial_t (\rho e(\rho, \vartheta)) + \nabla_x \cdot (\rho e(\rho, \vartheta) u) + \nabla_x \cdot q(\nabla_x \vartheta) = \mathbb{S}(\nabla_x u) : \nabla_x u - p(\rho, \vartheta) \nabla_x \cdot u + Q,$$9 after blow-up
Empirical measure	Monte Carlo sample average of deterministic solutions	Converges to statistical soln

The statistical solution formalism enables interpretation and analysis of NSF dynamics under uncertainty, accommodates loss of regularity, and ensures stability with respect to input data.

---

The developments in (Feireisl et al., 2022) rigorously unify classical stability, modern measure-valued/statistical solution theory, and uncertainty quantification for compressible Navier-Stokes-Fourier systems, establishing a robust analytic foundation for theoretical and applied investigations.