Compressible Navier-Stokes Equations: Density, Momentum, and Energy

The compressible Navier-Stokes equations describe flows where density changes. That makes the system larger than the incompressible equations: density becomes an unknown, pressure is thermodynamic, and energy usually has to be tracked.

This is the right framework for high-speed gas flow, acoustics, shock waves, combustion, and many atmospheric or astrophysical flows.

In conservation form, the compressible Navier-Stokes-Fourier system couples density $\rho$, velocity $u$, pressure $p$, and total specific energy $E=e+\tfrac12 |u|^2$.

Closure requires a constitutive law for viscous stress and heat flux, plus an equation of state such as $p=(\gamma-1)\rho e$ for an ideal gas.

The mass equation is

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

It says density changes when the flow compresses or expands fluid parcels. The momentum equation tracks how velocity changes under advection, pressure, viscosity, and external forces. The energy equation tracks heat, work, and kinetic energy.

A standard compressible viscous system is

$$\partial_t \rho + \nabla \cdot(\rho u)=0,$$

$$\partial_t(\rho u)+\nabla\cdot(\rho u\otimes u)+\nabla p=\nabla\cdot\tau+\rho f,$$

$$\partial_t(\rho E)+\nabla\cdot((\rho E+p)u)=\nabla\cdot(\tau u)+\nabla\cdot(\kappa\nabla\theta)+\rho f\cdot u.$$

For a Newtonian fluid, $\tau=\mu(\nabla u+\nabla u^T)+\lambda(\nabla\cdot u)I$. The second viscosity coefficient $\lambda$ and the common zero-bulk-viscosity assumption are explained under Stokes' hypothesis.

In compressible flow, pressure is tied to density and temperature through thermodynamics. In a simple ideal gas, $p=\rho R T$.

That is different from incompressible Navier-Stokes, where pressure mainly enforces the divergence-free constraint $\nabla\cdot u=0$.

The inviscid part of the compressible system is hyperbolic and has acoustic waves with finite sound speed $c=\sqrt{\partial p/\partial \rho|_s}$. The viscous and thermal terms add parabolic diffusion.

The inviscid Euler limit admits shock, rarefaction, and contact discontinuities. With nonzero viscosity and heat conduction, shocks and contacts are resolved into thin viscous or thermal transition layers under standard nondegenerate assumptions; this compressible wave structure still has no direct analogue in the incompressible model.

Compressibility usually matters when the Mach number is not small:

$$\mathrm{Ma}=\frac{|u|}{c}.$$

For many engineering flows, $\mathrm{Ma}<0.3$ is treated as effectively incompressible. Near and above the speed of sound, density changes and shock waves become central.

The incompressible equations arise formally as a low-Mach-number limit of the compressible equations under suitable preparation of the initial data. This limit is singular: acoustic waves become infinitely fast and the pressure law degenerates into the elliptic pressure constraint of incompressible flow.

The Clay Millennium Problem is not about the full compressible system. It is about the three-dimensional incompressible Navier-Stokes equations.

The compressible system has major open problems of its own, including global existence for large data and the interaction of shocks, viscosity, and vacuum states.

Compressible Navier-Stokes introduces difficulties absent from the Clay incompressible formulation: vacuum, acoustic modes, shock formation in inviscid limits, and thermodynamic closure. These are mathematically important but distinct from the 3D incompressible global regularity problem.

For the side-by-side comparison, read incompressible vs. compressible Navier-Stokes.