pub fn compute_vorticity(grid: &MacGrid) -> Vec<[f64; 3]>
Compute curl (vorticity) of the velocity field at cell centers.
Returns a flat buffer of 3-vectors (ωx, ωy, ωz).