View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
immiscible LB

5.0 | 3 ratings Rate this file 13 Downloads (last 30 days) File Size: 356 KB File ID: #24824 Version: 1.6

immiscible LB


G. Sken (view profile)


23 Jul 2009 (Updated )

Implements Immiscible Lattice Boltzmann (ILB, D2Q9) method for two phase flows

| Watch this File

File Information

The MatLab code uses lattice Boltzmann (LB) method to simulate two phase flows for immiscible fluids (blue & red fluids) in 2D according to the D2Q9 scheme.
The model is also referred to as color model or chromodynamic model or ILB (immiscible LB).
This MatLab implementation of ILB saves an AVI file in '' dir_avi='C:\' '' that records frames generated by the time evolution iterations.

The code implements in D2Q9 the model originally developed by Gunstensen [Gunstensen AK, Rothman D. Lattice Boltzmann model of immiscible fluids. Phys Rev A 1991;43(8):4320?4327], which is based on the original Rothman?Keller model [Rothman DH, Keller JM. Immiscible cellular automaton fluids. J Stat Phys 1988;52:1119?1127].
see also : Gunstensen, A., Rothman, D.H., Zaleski, S., and Zanetti, G., Physical Review A 43, 4320-4327, April 15, 1991.
page 190-191 of the reference book by S. Succi, The Lattice Boltzmann Equation, Oxford Univ. Press, 2001.
Only periodic boundary conditions are implemented in this version.
Also the method suggested by Dr Sergey Lishchuk for including surface tension is implemented, see :
S. V. Lishchuk, C. M. Care, I. Halliday. Lattice Boltzmann algorithm for surface tension with greatly reduced micro-currents. Phys. Rev. E, 2003,
In order to run this code you need MatLab 2009 or later with the image processing toolbox (IPT).
The files Run* are those to execute (eg. Run_Demo_ILBGK_D2Q9_Spinodal_Decomposition OR Run_ILB_Flow_inPM_Periodic_sigle_bubble).
If you stop the execution you can restart with >Restart_Time_Evolution

Required Products Image Processing Toolbox
MATLAB release MATLAB 7.9 (R2009b)
Other requirements IMAGE PROCESSING TOOLBOX tested on Windows XP 64
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (3)
09 May 2016 Zhenzhu MENG

03 Jul 2014 Isaac

Isaac (view profile)

14 Aug 2012 Victoria

28 Jul 2010 1.2

faster than previous submission

22 Sep 2010 1.5

new examples in Run* files

28 Sep 2010 1.6

corrects minor bug surface tension treatment

Contact us