← Circulation domains
PERIODIC NOZZLE · EULER FLOWBernoulli to Clebsch transport

One experiment · five layers

Why does a fluid
speed up here?

Start with a narrowing pipe. Keep the flux fixed, solve the slip-wall Euler problem, distinguish the velocity vector from its one-form, add vorticity with two scalar labels, and finally carry those labels with the flow.

CHAPTER 1 · CONSERVE VOLUMEthe same flux passes through every cross-section
area speed pressure
mean flux Q
flux spread
wide → throat speed
pressure drop
divergence RMS
vorticity RMS
harmonic circulation
tracked-label drift

The problem, then the representation

The constraint chooses φ.
Labels remember vorticity.

These are different jobs. Pressure enforces incompressibility in the Euler evolution; the spatial Clebsch reconstruction solves for a velocity potential φ. They are related by the constrained equations, but they are not the same scalar.

01 · POTENTIAL EULER BASE

The constriction chooses a harmonic representative.

closed, divergence-free, wall tangent

The solver begins with the closed one-form dx, Hodge-projects it against the curved walls, and keeps its period around the identified pipe. This is the two-dimensional potential-flow answer—not a resampled one-dimensional arrow field.

02 · BERNOULLI

Speed trades with pressure.

steady irrotational Euler

Pressure falls where the harmonic field accelerates and recovers when the channel widens. Because the model is inviscid, there is no permanent head loss.

03 · CLEBSCH RECONSTRUCTION

One global term, one exact term, two labels.

velocity one-form

At each time, the lab solves for φ so the reconstructed velocity is divergence-free and tangent to the wall. The coefficient cH is adjusted so the total harmonic period stays fixed; a single-valued φ cannot control that period.

04 · VORTICITY

Take one exterior derivative.

vorticity two-form

Both du and dα ∧ dβ are two-forms on this planar surface. Relative to the chosen area form they can be displayed as a signed scalar density; as an ambient picture, the associated vorticity vector points normal to the plane.

Vector or covector?

The arrow moves a particle.
The one-form measures a displacement.

VECTOR
particle velocity

Feed the tangent vector u to the advection ODE. In Chapter 3 it is drawn as a cyan arrow.

COVECTOR
circulation density

Feed the one-form u a test displacement. The orange bar is its kernel: moving along that bar contributes zero.

METRIC
lower the index

The Euclidean metric makes their components look identical here, but the objects still transform differently and do different jobs.

Chapter 5 · the algorithm

Carry α and β.
Re-solve φ.

The simulation repeats a material transport step followed by a spatial Hodge reconstruction. This makes the division of labor visible.

  1. 01
    Advect the scalar labels.
    material transport
  2. 02
    Build the vortical candidate.
    label one-form
  3. 03
    Solve the incompressibility constraint.
    Neumann Hodge solve
  4. 04
    Restore the prescribed global period.
    harmonic + exact + label pieces

Next: shallow water on a surface

Let the cross-section become state.

Here the wall shape A(x) is prescribed and divergence is projected away. In shallow water, the layer depth h evolves, horizontal convergence changes h, and pressure comes from the hydrostatic potential one half g h squared.

shallow-water continuity plus Clebsch transport
Continue to Clebsch shallow water →