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:
Seb (2020). Lattice Boltzmann N-phase model with MRT and Enhanced Equilibrium (LBM) (https://www.mathworks.com/matlabcentral/fileexchange/49565-lattice-boltzmann-n-phase-model-with-mrt-and-enhanced-equilibrium-lbm), MATLAB Central File Exchange. Retrieved .