Statement

For a smooth solution of the , set e=∣u∣2/2e=|u|^2/2. The local kinetic-energy balance is

∂te+∇⋅((e+p)u)=νΔe−ν∣∇u∣2+f⋅u.\partial_t e+\nabla\cdot((e+p)u)=\nu\Delta e-\nu|\nabla u|^2+f\cdot u.

Here pressure and force are normalized by the constant density.

Derivation

Take the scalar product of the momentum equation with uu. The identities u⋅∂tu=∂teu\cdot\partial_tu=\partial_te, u⋅(u⋅∇)u=∇⋅(eu)u\cdot(u\cdot\nabla)u=\nabla\cdot(eu), u⋅∇p=∇⋅(pu)u\cdot\nabla p=\nabla\cdot(pu), and u⋅Δu=Δe−∣∇u∣2u\cdot\Delta u=\Delta e-|\nabla u|^2 give the equation. The divergence identities use ∇⋅u=0\nabla\cdot u=0.

For a rough distributional solution, these multiplications require justification; this smooth equality does not automatically persist.