Compressible Navier–Stokes flow around a slip-wall bluff body — CPU translation
of the p5.js Bluff Body Slip sketch (N = 100, ~340 lines).
The body is a 20×20 square excluded from the domain with slip
walls on all four faces. Drag/lift are computed by integrating the
pressure difference across the body.
Equations: ∂ρ/∂t + ∇·m = 0,
∂m/∂t + ∇·(m u) + ∇p = ∇·(ν∇u),
∂e/∂t + ∇·(e u) + p∇·u = ν|∇u|² with
p = γe, u = m/ρ.
Cross-cuts at j = N/2 — Cyan: ρ. Yellow-orange: ux. Red: p. Green: T = e/ρ − 1.
Blue dots: particle traces.
Blue square: bluff body.
Dcum = ∫∫ν|∇u|² dx dt — cumulative dissipation.