2D Incompressible Navier–Stokes (Projection · MAC) — SVG Heatmap v0.3.2

Simulation Controls

160
Cells along the nozzle axis (x). In a MAC (Marker‑and‑Cell) grid, higher values sharpen axial features but increase cost.
60
Cells across height (y). In a MAC grid, higher values resolve shear layers and near‑wall gradients.
0.020
Controls viscous diffusion (μ/ρ). Lower ν → higher Reynolds number at the throat.
0.250
Bulk inflow magnitude. With pulsatile mode: \(U_{in}(t)=U_0[1+A\sin(2\pi f t)]\).
Lc 0.30
Ld 0.30
0.35
With symmetry enabled, Ld mirrors Lc. “Cosine‑smooth” eases the converging/diverging ramps.
Term maps reveal dominant mechanisms: advection (jet), |∇p| (throat & diffuser), viscous diffusion (walls), vorticity (shear/recirculation).
Locking preserves min/max across frames so time evolution becomes visually obvious.
Bright white arrows with heads for visibility on the heatmap.
Using downstream‑referenced cross‑section mean clarifies pressure recovery along the diffuser.
Ready.

2D CFD Heatmap (velocity / pressure / shear / NS terms)

Velocity magnitude Kinematic pressure Shear & NS term maps
Low (blue) Moderate (red) High (yellow)
Throat is centered at mid‑domain (x = 0.5). MAC grid; explicit advection–diffusion and pressure projection (incompressible, constant ν).

Centerline |v|(x), pressure p(x), and upper‑wall τw(x)

p (arb.) & τw (arb.) |v| (arb.) x (0…1, throat @ 0.5)
Velocity |v|(x) — blue Pressure pc(x) — green Cross‑section mean p̄(x) — dashed green Violet dashed = upper‑wall τw(x)

Estimator (throat)

Estimator (throat)Value (arb./mmHg*)
Peak Doppler (4V2)
Modified (4(Vth2−Vin2))
Doppler + pressure‑recovery
CFD net Δp (up−down)
*CFD pressure is kinematic (per ρ); mmHg shown for relative comparison only.

Model, intuition, and clinical interpretation

Governing equations (incompressible, Newtonian)

\[ \nabla\!\cdot\!\mathbf{u}=0,\qquad \frac{\partial\mathbf{u}}{\partial t}+(\mathbf{u}\!\cdot\!\nabla)\mathbf{u} = -\nabla p + \nu \nabla^2 \mathbf{u}. \]

Navier–Stokes (two useful forms)

  1. Convective (non‑conservative) form \[ \nabla\!\cdot\!\mathbf{u}=0,\quad \frac{\partial\mathbf{u}}{\partial t}+(\mathbf{u}\!\cdot\!\nabla)\mathbf{u} = -\,\nabla p + \nu \nabla^2 \mathbf{u}. \]
    Observed via heatmap toggles: |(u·∇)u| (advection), |∇p| (pressure), |ν∇²u| (viscous). Advection dominates the throat jet; viscous terms and vorticity concentrate at walls and separation zones.
  2. Conservative (flux/stress) form \[ \frac{\partial(\rho\mathbf{u})}{\partial t} + \nabla\!\cdot\!\big(\rho\,\mathbf{u}\otimes\mathbf{u} + p\,\mathbf{I} - \mu(\nabla\mathbf{u}+\nabla\mathbf{u}^{\!\top})\big) = \mathbf{0},\quad(\rho{=}1,\ \mu{=}\nu). \]
    Momentum fluxes (convective + pressure + viscous stresses) balance over any control volume. The measured Δpnet reflects irreversible loss from viscous dissipation and separation; large recirculation and high τw predict reduced pressure recovery.

Reading pressure recovery

Enhancing pressure recovery

Interpreting the four throat estimates

Live interpretation helper (auto‑updates)