This lattice Boltzmann code is an extension of one of my previous submission ( http://www.mathworks.com/matlabcentral/fileexchange/40607-lattice-boltzmann-immiscible-two-phase-model--lbm- ). It may give a better understanding of the multiphase RK-type lattice Boltzmann model. The numerical algorithm is constructed in such a way that the computer memory required is significantly reduced. With large N, the distribution functions of each color are not required [R1]. Variable density ratios are achieved with enhanced equilibrium distribution functions [R2]. The multiple-relaxation-time collision operator is also available [R3]. Two qualitative simulations are already pre-set:
1) An N-phase spinodal decomposition; and
2) a two-phase oscillating bubble.
This code is based on the following publications:

[R1] Leclaire, S., Reggio, M. and Trépanier, J.-Y. (2013) Progress and investigation on lattice Boltzmann modeling of multiple immiscible fluids or components with variable density and viscosity ratios. Journal of computational physics 246: 318-342.

[R2] Leclaire, S., Pellerin, N., Reggio, M. and Trépanier, J.-Y. (2013) Enhanced equilibrium distribution functions for simulating immiscible multiphase flows with variable density ratios in a class of lattice Boltzmann models. International Journal of Multiphase Flow 57: 159-168.

[R3] Leclaire, S., Pellerin, N., Reggio, M. and Trépanier, J.-Y. (2014) Unsteady immiscible multiphase flow validation of a multiple-relaxation-time lattice Boltzmann method. Journal of Physics A: Mathematical and Theoretical 47(10): 105501.

If these m-files were useful in your work, please cite:
1) The above three references R1, R2 and R3 accordingly.

Others related publications may be found there:
http://scholar.google.ca/citations?user=_wUUbqwAAAAJ&hl=en