Incompressible vs. Compressible Navier-Stokes

"Incompressible vs compressible" boils down to density. Does it stay constant, or does it change?

Try it. Fill a syringe with water and push the plunger. The water moves, but under everyday conditions it doesn't compress noticeably. Water resists compression so strongly that treating it as incompressible is an excellent approximation. Incompressible. Now fill that syringe with air and seal the end. Push the plunger in and you'll feel the air give way, the same mass of air now packed into less volume as it compresses under your thumb. That's compressible flow.

In incompressible flow, the density $\rho$ is constant throughout the fluid, and every tiny parcel keeps its volume as it moves through space. Compressible flow is different. Density becomes a variable, free to change from place to place and moment to moment. Air around a jet engine, gas in an explosion, the atmosphere at large scales: all compressible, all driven by density variations.

Why care? Because this distinction reshapes what the Navier-Stokes equations look like, what they predict, and how difficult they are to analyze and solve.

The incompressible assumption sets the density field $\rho$ to a positive constant throughout the flow domain. Physically, this means that every material volume element preserves its volume under the flow map. The compressible setting promotes $\rho(x,t)$ to a full unknown, governed by its own evolution equation.

These aren't two notations for the same system. They differ in the number of unknowns, the structure of the constraint equations, and the character of the pressure. The incompressible system has $d+1$ scalar unknowns ($u$ and $p$ in $\mathbb{R}^d$); the compressible system typically has $d+2$ primary fields (velocity, density, and internal energy or temperature), with pressure determined through an equation of state.

The distinction also changes the PDE type of the pressure: elliptic in the incompressible model, where pressure is determined by an instantaneous spatial solve, versus mixed hyperbolic-parabolic in the compressible case, where pressure disturbances propagate at finite speed. This isn't a technicality. It controls how information moves through the fluid.

The incompressible Navier-Stokes equations describe fluids whose density is constant. They're the version that appears in the Clay Millennium Problem and the version this site focuses on.

The system has two parts. The momentum equation:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

and the incompressibility constraint:

$$\nabla \cdot u = 0$$

The constraint $\nabla \cdot u = 0$ says the velocity field is divergence-free: fluid neither piles up nor thins out anywhere. Whatever flows into a tiny region must flow out at the same rate. This single condition replaces the entire density equation. Density doesn't change, so you don't need an equation to track it.

Pressure plays a special role here. It isn't determined by a thermodynamic law (like the ideal gas law). Instead, it adjusts instantaneously everywhere to keep the flow divergence-free. Mathematically, $p$ solves a Poisson equation derived from the constraint. Pressure changes propagate infinitely fast. There's no "speed of sound" in incompressible flow.

The incompressible Navier-Stokes system has two unknown fields: velocity $u$ and pressure $p$. That simplicity is deceptive. The nonlinear term $(u \cdot \nabla)u$ still makes the system extremely difficult in three dimensions.

The incompressible Navier-Stokes system on $\mathbb{R}^3$ with kinematic viscosity $\nu > 0$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

The divergence-free condition $\nabla \cdot u = 0$ is a pointwise constraint, not an evolution equation. It encodes local volume preservation: the flow map $\Phi_t$ satisfies $\det(D\Phi_t) = 1$ for all $t$.

Applying the divergence operator to the momentum equation and using incompressibility yields the pressure Poisson equation:

$$-\Delta p = \partial_i \partial_j (u_i u_j) - \nabla \cdot f.$$

This is an elliptic equation for $p$ at each fixed time. The pressure isn't an independent thermodynamic variable; in the standard PDE interpretation, it acts as a Lagrange multiplier enforcing the divergence-free constraint, determined globally and instantaneously by the velocity field. Information propagates at infinite speed through the pressure, a structural difference from the compressible system that can't be papered over.

The unknowns are $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ and $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$. For the Clay formulation (Fefferman, 2000), $u_0 \in C^\infty(\mathbb{R}^3)$ is divergence-free and the question is whether $u$ remains in $C^\infty(\mathbb{R}^3 \times [0,\infty))$ with bounded energy.

The compressible Navier-Stokes equations govern flows where density varies. Bigger system. More unknowns. More equations.

You still have a momentum equation, but now density $\rho$ appears explicitly:

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

The constraint $\nabla \cdot u = 0$ is gone. In its place, you get a continuity equation that tracks how density evolves:

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

This says mass is conserved: density changes because the flow compresses or expands fluid parcels.

The system also needs an energy equation and an equation of state, a thermodynamic relation like $p = \rho R T$ (the ideal gas law) that ties pressure to density and temperature. Pressure is no longer a passive enforcer of a constraint. It has its own physics, its own dynamics, and it propagates at a finite speed: the speed of sound.

The compressible system is essential for aerodynamics at high speeds, astrophysical gas dynamics, combustion, and any flow where density changes matter. But it's a genuinely different mathematical object from the incompressible equations. More unknowns, more equations, different PDE structure entirely.

The compressible Navier-Stokes system couples the velocity $u(x,t)$, density $\rho(x,t)$, pressure $p(x,t)$, and specific internal energy $e(x,t)$ (or temperature $\theta$). In conservation form:

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

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

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

where $E = e + \tfrac{1}{2}|u|^2$ is the total specific energy, $\tau$ is the viscous stress tensor (for a Newtonian fluid, $\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$ with bulk viscosity $\lambda$), and $\kappa$ is thermal conductivity.

Closure requires an equation of state, e.g., $p = (\gamma - 1)\rho e$ for an ideal gas with adiabatic index $\gamma$.

In the compressible system, acoustic disturbances propagate at finite speed, with characteristic sound speed $c = \sqrt{\partial p / \partial \rho |_s}$. This contrasts with the elliptic, instantaneous pressure determination in the incompressible model. The compressible equations support shock waves, rarefactions, and contact discontinuities that have no analogue in incompressible flow.

When does compressibility matter? One number decides: the Mach number.

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$ is flow speed. $c$ is the speed of sound. Their ratio tells you how fast the flow moves compared to the speed at which pressure disturbances can propagate through the medium, and that comparison determines whether you can safely ignore density changes or whether they'll dominate the physics.

When $\text{Ma} < 0.3$, density changes by less than about 5%. Incompressible equations work. Air in a room, water in a pipe, wind around a building: all low-Mach flows where pressure disturbances travel so much faster than the flow itself that density barely budges.

Above $\text{Ma} \approx 0.3$, compressibility starts to bite, and around $\text{Ma} \approx 1$ you hit the transonic regime where local supersonic pockets appear and shock waves form. Fighter jets. Rocket nozzles. Re-entering spacecraft.

Not a binary switch. Most everyday fluid flows, and the Clay Millennium Problem, sit firmly in the low-Mach regime where the incompressible equations apply.

The Mach number $\text{Ma} = |u|/c$ parametrizes the importance of compressibility, where $c = \sqrt{\partial p / \partial \rho |_s}$ is the isentropic sound speed. Formally, the incompressible equations arise as the low-Mach limit of the compressible system.

The asymptotic expansion in powers of $\text{Ma}^2$ (see Klainerman & Majda, 1981, 1982; Schochet, 1986) shows that as $\text{Ma} \to 0$ with suitable initial data, the compressible solutions converge to the incompressible solution. In standard low-Mach nondimensionalizations with well-prepared data, one often writes pressure as a nearly spatially uniform thermodynamic background plus a smaller dynamic correction that enforces the incompressibility constraint in the limit.

This is a singular limit: the sound speed $c \to \infty$ and the compressible system's hyperbolic character degenerates to the elliptic pressure equation of incompressible flow. The acoustic modes become infinitely fast and decouple from the vortical dynamics.

The regime $\text{Ma} < 0.3$ is an engineering heuristic. It reflects the empirical observation that relative density variation $\delta\rho / \rho \sim \text{Ma}^2 / 2$ stays below $\sim$5% in this range. The mathematical justification is the convergence theorem for the low-Mach limit, which requires well-prepared initial data and compatible boundary conditions.

The Clay Millennium Problem asks a precise question: given a smooth, divergence-free initial velocity on $\mathbb{R}^3$, does the incompressible Navier-Stokes system always produce a smooth solution that exists for all time?

Why incompressible specifically? Three reasons.

First, it's already hard enough. The incompressible 3D equations have resisted proof of global regularity since Leray's foundational work in 1934. Adding variable density, thermodynamics, and shock waves would make the problem vastly harder, not more tractable.

Second, the difficulty is pure fluid mechanics. The incompressible system isolates the core mathematical challenge, the competition between nonlinear advection $(u \cdot \nabla)u$ and viscous dissipation $\nu \Delta u$, without thermodynamic or acoustic complications. It's the cleanest arena to ask the regularity question.

Third, the physics is clean. The incompressible equations model the most common everyday flows. Whether they can produce singularities from smooth data is a fundamental question about the mathematical consistency of classical fluid mechanics.

The compressible system has its own deep open problems (existence of global solutions with large data, formation and interaction of shocks), but those are different problems with different structures. The Clay prize targets the incompressible case because that's the specific regularity question Fefferman formulated for 3D Navier-Stokes.

The official Clay formulation (Fefferman, 2000) specifies the incompressible system on $\mathbb{R}^3$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

with $u_0 \in C^\infty(\mathbb{R}^3)$ divergence-free, and the question is whether $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ with $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ bounded for all $t \geq 0$.

The choice of the incompressible system is mathematically motivated. The key open problem, the gap between Leray-Hopf weak solutions (which exist globally but may not be unique or smooth) and classical smooth solutions (which exist locally but may blow up), is specific to the incompressible 3D equations. In 2D, global regularity for smooth incompressible Navier-Stokes solutions is known; the 3D case remains open.

The compressible system introduces qualitatively different difficulties: shock formation (which occurs even for Euler equations with smooth data), vacuum states ($\rho \to 0$), and the coupling between vorticity and acoustic modes. These are important open problems, but they're structurally distinct from the incompressible regularity question.

The incompressible problem isolates the competition between the energy-supercritical nonlinearity and viscous dissipation. The natural energy estimate gives $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$, which falls short of the scaling-critical space $L^\infty_t \dot{H}^{1/2}_x$ by half a derivative in 3D. Closing this gap, or proving it can't be closed, is the heart of the Millennium Problem.

Start here. Want every term in the incompressible system pulled apart, with the physical meaning and mathematical role of each piece explained from scratch? What Are the Navier-Stokes Equations?

Where does this system come from? Derivation of the Navier-Stokes Equations.

Drop viscosity and you get the Euler equations, which are a century older, look simpler on the page, and in some ways are even harder to understand mathematically because you lose the smoothing effect of the diffusion term. Euler vs. Navier-Stokes.

The prize. The Navier-Stokes Existence and Smoothness Problem.

Suggested paths from here:

• What Are the Navier-Stokes Equations?, the full incompressible system with term-by-term analysis
• Derivation of the Navier-Stokes Equations, from continuum mechanics axioms to the PDE system
• Euler vs. Navier-Stokes, the inviscid limit $\nu \to 0$ and the role of viscosity in regularity
• The Navier-Stokes Existence and Smoothness Problem, the precise Clay formulation, Leray-Hopf theory, and the gap between weak and strong solutions