The lattice Boltzmann method in 25 lines of MATLAB code (LBM)

Here is two m-files showing the simplicity of LBM. A view of the density distributions is provided.

When "LBM_25_LINES.m" is started, the LBM solves the flow field in a "pseudo-porous" geometry with periodic boundary condition. A perturbation to the velocity field in the x-direction is applied so that the flow moves from the left to the right side of the computational domain.
When "LBM_VIEW.m" is started, the m-file script shows all LBM process at the mesoscopic level. At each time step, the collision, the bounce-back and the streaming step are showed with details and at the same time the macroscopic velocity field is plotted.

Note that these two m-files are not research code! Although, they might be use as a starting point.

Overall, I believe that these two m-files may be of a good addition for a general audience scientific presentation.

Enjoy!

If interested, my scientific publications may be found there:
http://scholar.google.ca/citations?user=_wUUbqwAAAAJ&hl=en