Navier–Stokes equations

Navier–Stokes equations are a set of partial differential equations that describe the motion of Newtonian fluids. They are fundamental to the field of continuum mechanics and are named for their principal developers, Claude-Louis Navier and George Gabriel Stokes. These equations govern a vast range of phenomena, from airflow over an aircraft wing to ocean currents and blood flow in the human circulatory system.

Mathematical formulation

The equations express the conservation of momentum and mass for a fluid continuum. In their most common form for an incompressible flow of a Newtonian fluid, they are written using vector calculus notation. The momentum equation is \(\rho \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \mathbf{f}\), where \(\mathbf{v}\) is the flow velocity, \(p\) is the pressure, \(\rho\) is the density, and \(\mu\) is the dynamic viscosity. The term \(\mathbf{f}\) represents external body forces, such as gravity. This is coupled with the continuity equation \(\nabla \cdot \mathbf{v} = 0\), which enforces mass conservation. The full system for compressible flow includes an additional equation of state, such as the ideal gas law, and an energy equation.

Physical interpretation

The equations balance the inertia of fluid particles against the forces acting upon them. The term \(\rho \frac{D \mathbf{v}}{Dt}\) represents the material derivative, accounting for the total acceleration of a fluid parcel. This is equated to forces from pressure gradients \(-\nabla p\), viscous stresses \(\mu \nabla^2 \mathbf{v}\), and any external influences like the Coriolis force in geophysical fluid dynamics. The viscous term models internal friction, originating from the no-slip condition at solid boundaries, a concept solidified by George Gabriel Stokes. This interplay determines complex behaviors like the transition from laminar flow to turbulence.

Derivation and assumptions

The derivation begins with Cauchy's momentum equation, a general statement of momentum conservation for a continuum. Applying constitutive equations specific to a Newtonian fluid, which linearly relates stress tensor to the strain rate tensor, yields the standard form. Key assumptions include the fluid being a continuum, obeying Newton's laws of motion, and having a viscosity that is constant and isotropic. The equations also assume the Stokes hypothesis regarding the bulk viscosity. For incompressible flow, the additional assumption of constant density is applied, greatly simplifying the continuity equation.

Solutions and existence

Analytical solutions are known only for simplified geometries and conditions, such as Poiseuille flow in a pipe or Couette flow between plates. For most practical situations, solutions are approximated using computational fluid dynamics methods like the finite element method or lattice Boltzmann method. The question of whether smooth, physically reasonable solutions always exist in three dimensions, known as the Navier–Stokes existence and smoothness problem, is one of the Millennium Prize Problems designated by the Clay Mathematics Institute. This profound mathematical challenge remains unsolved.

Applications

The equations are the cornerstone of modeling in numerous scientific and engineering disciplines. In aerospace engineering, they are solved to design airfoils and predict lift and drag forces. Meteorologists use them in global climate models to simulate atmospheric circulation and weather forecasting. They are essential in biomechanics for simulating hemodynamics in the aorta and in chemical engineering for reactor design. The Reynolds-averaged Navier–Stokes equations are a critical tool for modeling turbulent flow in industrial applications.

Special cases and simplifications

Under specific conditions, the full equations reduce to more tractable forms. For very slow, viscous flows, the Stokes equations neglect inertial terms entirely. The Euler equations describe inviscid flow by setting viscosity to zero, applicable to high-Reynolds number flows outside boundary layers. The boundary layer equations, developed by Ludwig Prandtl, simplify the analysis near surfaces. For potential flow, the velocity field is derived from a velocity potential, satisfying the Laplace equation. The Boussinesq approximation is used in natural convection problems to model density variations.