Mass
U₁ = ρ
Density is conserved mass per unit volume. Its flux is ρu.
Build a discontinuity, let mass, momentum and energy evolve, and inspect the density, velocity and pressure solution. Adapted directly from the supplied Python solver.
Density is conserved mass per unit volume. Its flux is ρu.
Momentum flux combines advection and pressure: ρu² + p.
Energy flux is u(E+p): internal and kinetic energy plus pressure work.
At every interface, average the physical flux from both sides, then subtract a dissipative term proportional to the state jump and maximum local wave speed.
The CFL condition chooses Δt = CFL Δx / max(|u|+a). Each cell is updated from its outgoing and incoming interface fluxes.