Formal PDE Models for Electromagnetism, Plasma, and Fluids

4 Incompressible Navier-Stokes

The incompressible Navier-Stokes equations, Euler equations (inviscid limit), vorticity dynamics, and the pressure Poisson equation.

4.1 Equations of Motion

Definition 4.1 FluidConstants
#

Fluid constants: density \(\rho {\gt} 0\), dynamic viscosity \(\mu \ge 0\).

Lemma 4.2 FluidConstants.\(\rho \)_ne_zero
#

\(\rho \ne 0\).

Lemma 4.3 FluidConstants.\(\rho \)_pos
#

\(\rho {\gt} 0\).

Definition 4.4 FluidConstants.\(\nu \)
#

Kinematic viscosity \(\nu = \mu /\rho \).

Lemma 4.5 FluidConstants.\(\nu \)_nonneg
#

\(\nu \ge 0\).

Incompressible Navier-Stokes: \(\nabla \! \cdot \! \mathbf{v} = 0\) and \(\rho \frac{D \mathbf{v}}{D t} = -\nabla p + \mu \nabla ^2\mathbf{v} + \mathbf{f}\).

\(\frac{D \mathbf{v}}{D t} = -\frac{1}{\rho }\nabla p + \nu \nabla ^2\mathbf{v} + \frac{1}{\rho }\mathbf{f}\).

Theorem 4.8 Continuity from incompressible

\(\nabla \! \cdot \! (\rho \mathbf{v}) = 0\) since \(\rho \) is constant and \(\nabla \! \cdot \! \mathbf{v} = 0\).

4.2 Euler Equations

Theorem 4.9 Euler momentum

When \(\mu = 0\): \(\rho \frac{D \mathbf{v}}{D t} = -\nabla p + \mathbf{f}\).

Theorem 4.10 Euler momentum per unit mass

When \(\mu = 0\): \(\frac{D \mathbf{v}}{D t} = -\frac{1}{\rho }\nabla p + \frac{1}{\rho }\mathbf{f}\).

4.3 Vorticity

Definition 4.11 vorticity
#

Vorticity: \(\mathbf{\omega }(t,x) = \nabla \! \times \! \mathbf{v}(t,x)\).

Theorem 4.12 div vorticity = 0

\(\nabla \! \cdot \! \mathbf{\omega } = 0\): vorticity is solenoidal.

Theorem 4.13 curl pressure gradient = 0

\(\nabla \! \times \! (\nabla p) = 0\).

Definition 4.14 NSSmoothness

Smoothness bundle for vorticity/pressure derivations: \(C^3\) velocity, curl-time commutativity, and divergence-time commutativity.

Lemma 4.15 partialDerivComp_differentiable

Partial derivatives of \(C^2\) vector fields are differentiable.

Theorem 4.16 curl of advective term

\(\nabla \! \times \! ((\mathbf{v}\cdot \nabla )\mathbf{v}) = (\mathbf{v}\cdot \nabla )\mathbf{\omega } - (\mathbf{\omega }\cdot \nabla )\mathbf{v}\) for divergence-free \(\mathbf{v}\). The cornerstone identity for vorticity dynamics.

Theorem 4.17 curl commutes with Laplacian

\(\nabla \! \times \! (\nabla ^2\mathbf{v}) = \nabla ^2(\nabla \! \times \! \mathbf{v})\) for sufficiently smooth incompressible velocity fields.

\[ \frac{\partial \mathbf{\omega }}{\partial t} + (\mathbf{v}\cdot \nabla )\mathbf{\omega } = (\mathbf{\omega }\cdot \nabla )\mathbf{v} + \nu \nabla ^2\mathbf{\omega } + \frac{1}{\rho }\nabla \! \times \! \mathbf{f}. \]

The vortex stretching term \((\mathbf{\omega }\cdot \nabla )\mathbf{v}\) is unique to 3D.

4.4 Pressure Poisson Equation

Theorem 4.19 Pressure Poisson equation
\[ \nabla ^2_sp = -\rho \sum _{i,k}\frac{\partial v_k}{\partial x_i}\frac{\partial v_i}{\partial x_k} + \nabla \! \cdot \! \mathbf{f}. \]

The pressure is determined as a constraint enforcing incompressibility.