D2Q9 LBM simulation with bounce-back walls — draw obstacles!
D2Q9 lattice Boltzmann with BGK collision: f_i(t+1,x+e_i) = f_i + (f^eq_i - f_i)/τ. Bounce-back for no-slip walls. Click to draw obstacles. Parabolic Poiseuille profile emerges.