The lattice Boltzmann method in NumPy: vortex shedding behind a cylinder

Afterwards you can write a lattice Boltzmann solver in NumPy, set its relaxation time from the Reynolds number, and measure a vortex shedding frequency.