Kelvin-Helmholtz instability at different Reynolds numbers with NGSolve

Simulation of a Kelvin-Helmholtz instability problem at different Reynolds numbers Re=1E2, Re=1E3 and Re=1E4 with NGSolve. Computed with a high-order exactly divergence-free H(div)-HDG method. Visualisation of the vorticity. More information available at