Skip to main content

Navier-Stokes equations for Newtonian fluid

Short answer

The Navier-Stokes equations follow from Newton's second law for a fluid element, with viscous stresses proportional to strain rates. The x-component contains inertia, pressure, gravity, and viscous terms. Setting viscosity to zero gives Euler's equation; setting velocity to zero yields the hydrostatic equation.

Explain how the Navier-Stokes equations are obtained for a Newtonian incompressible fluid. Write the x-component, identify each term, state the assumptions, and show how the equation reduces to (i) Euler's equation and (ii) the hydrostatic equation.

Basis

The Navier-Stokes (N-S) equations are Newton's second law written for a fluid element, with the viscous stresses related to the velocity gradients by Stokes' relations for a Newtonian fluid:

τxy=μ(∂u∂y+∂v∂x),σxx=−p+2μ∂u∂x\tau_{xy} = \mu\left(\frac{\partial u}{\partial y} + \frac{\partial v}{\partial x}\right), \qquad \sigma_{xx} = -p + 2\mu\frac{\partial u}{\partial x}

Steps

  1. Take a fluid element dx dy dzdx\,dy\,dz with normal and shear stresses on its faces.
  2. Net surface force in the xx-direction == (net normal stress gradient + shear stress gradients) × dx dy dz\times\ dx\,dy\,dz.
  3. Add the body force ρgx dx dy dz\rho g_x\,dx\,dy\,dz.
  4. Equate the sum to mass ×\times acceleration, ρ DuDt dx dy dz\rho\,\dfrac{Du}{Dt}\,dx\,dy\,dz.
  5. Substitute the Newtonian stress relations, and use incompressible continuity (∇⋅V⃗=0\nabla\cdot\vec V = 0) with constant μ\mu to simplify.

x-component (incompressible, constant viscosity)

ρ(∂u∂t+u∂u∂x+v∂u∂y+w∂u∂z)=ρgx−∂p∂x+μ(∂2u∂x2+∂2u∂y2+∂2u∂z2)\rho\left(\frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} + w\frac{\partial u}{\partial z}\right) = \rho g_x - \frac{\partial p}{\partial x} + \mu\left(\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2}\right)

In vector form: ρDV⃗Dt=ρg⃗−∇p+μ∇2V⃗\rho\dfrac{D\vec V}{Dt} = \rho\vec g - \nabla p + \mu\nabla^2\vec V.

TermMeaning
ρ ∂u/∂t\rho\,\partial u/\partial tLocal inertia force (unsteadiness)
ρ(V⃗⋅∇)u\rho(\vec V\cdot\nabla)uConvective inertia force
ρgx\rho g_xBody force (gravity)
−∂p/∂x-\partial p/\partial xPressure force
μ∇2u\mu\nabla^2uViscous force

The y and z components are written in the same way.

Assumptions

Newtonian fluid, incompressible flow, constant viscosity and a continuum. The equations are exact under these assumptions: they apply to laminar flow and to the instantaneous velocities in turbulent flow.

Reduction

  • Euler's equation: for an inviscid fluid, μ=0\mu = 0:
ρDV⃗Dt=ρg⃗−∇p\rho\frac{D\vec V}{Dt} = \rho\vec g - \nabla p
  • Hydrostatic equation: for a fluid at rest, V⃗=0\vec V = 0, and with zz upward, ∂p/∂z=−ρg\partial p/\partial z = -\rho g and ∂p/∂x=∂p/∂y=0\partial p/\partial x = \partial p/\partial y = 0.

The N-S equations have exact solutions only for simple cases, such as laminar flow between parallel plates and in circular pipes (Hagen-Poiseuille flow).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Fluid Mechanics chapter 4 questions →